Next Article in Journal
Inactivation Effect of Antibiotic-Resistant Gene Using Chlorine Disinfection
Next Article in Special Issue
Incorporating the Vadose Zone into the Budyko Framework
Previous Article in Journal
Long Term Trends and Dynamics of Dissolved Organic Carbon (DOC) in a Subtropical Reservoir Basin
Previous Article in Special Issue
HYDRUS Simulation of Sustainable Brackish Water Irrigation in a Winter Wheat-Summer Maize Rotation System in the North China Plain
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influencing Factors and Simplified Model of Film Hole Irrigation

1
College of Water Resources and Architectural Engineering, Northwest A&F University, Yangling 712100, China
2
Key Laboratory of Agricultural Soil and Water Engineering in Arid and Semiarid Areas, Ministry of Education, Northwest A&F University, Yangling 712100, China
3
College of Energy and Power Engineering, Lanzhou University of Technology, Lanzhou 730050, China
*
Author to whom correspondence should be addressed.
Water 2017, 9(7), 543; https://doi.org/10.3390/w9070543
Submission received: 24 May 2017 / Revised: 28 June 2017 / Accepted: 17 July 2017 / Published: 20 July 2017
(This article belongs to the Special Issue Water and Solute Transport in Vadose Zone)

Abstract

:
Film hole irrigation is an advanced low-cost and high-efficiency irrigation method, which can improve water conservation and water use efficiency. Given its various advantages and potential applications, we conducted a laboratory study to investigate the effects of soil texture, bulk density, initial soil moisture, irrigation depth, opening ratio (ρ), film hole diameter (D), and spacing on cumulative infiltration using SWMS-2D. We then proposed a simplified model based on the Kostiakov model for infiltration estimation. Error analyses indicated SWMS-2D to be suitable for infiltration simulation of film hole irrigation. Additional SWMS-2D-based investigations indicated that, for a certain soil, initial soil moisture and irrigation depth had the weakest effects on cumulative infiltration, whereas ρ and D had the strongest effects on cumulative infiltration. A simplified model with ρ and D was further established, and its use was then expanded to different soils. Verification based on seven soil types indicated that the established simplified double-factor model effectively estimates cumulative infiltration for film hole irrigation, with a small mean average error of 0.141–2.299 mm, a root mean square error of 0.177–2.722 mm, a percent bias of −2.131–1.479%, and a large Nash–Sutcliffe coefficient that is close to 1.0.

1. Introduction

Severe global water shortages have attracted ample attention because water is crucial for human beings, particularly in agriculture [1,2,3]. Notably, the arid and semiarid regions of the western part of China have encountered serious water scarcity problems because of limited rainfall and great soil-moisture evaporation [4,5]. The use of advanced water-saving irrigation methods, such as sprinkler irrigation, microsprinkler irrigation, surface drip irrigation, and subsurface drip irrigation, is encouraged, particularly in these arid and semiarid areas of China with high rates of evaporation. However, the cost of advanced irrigation systems is relatively high for farmers in underdeveloped areas [6]. Therefore, a low-cost and efficient irrigation method is urgently required.
Currently, plastic mulching has become a globally applied agricultural practice for high yields and increased water-use efficiency. This method is primarily used to protect seedlings in arid climates, prevent evaporation, maintain or slightly increase soil temperature and humidity, prevent weed growth, and reduces herbicide and fertilizer use [7,8,9,10,11,12]. These prospects have increased the immediate relevance of plastic film technology; currently, plastic films constitute the largest proportion of agricultural surface coverings [13,14]. According to the National Bureau of Statistics of China [15], China is the world’s largest film mulch consumer and boasts the world’s largest film mulch coverage area. The area covered by film mulch considerably increased to 20 million hm2 in 2014, and its use increased fourfold (from 642 to 2580 megatons) from 1991 to 2014. Film hole irrigation has also developed along with film-mulched techniques for crops with wide line spacing [16]. Film hole irrigation is a relatively new irrigation method that involves completely covering a bordered field with plastic film with holes of uniform size [17]. Water penetrates into the soil through the holes during irrigation, and seedlings sprout through these holes on germination. Compared with the traditional surface irrigation method, this technique considerably reduces water losses and improves the uniformity of irrigation along the long direction of the border [6,16,18,19].
Mathematical models and software, such as winSRFR [20], have been widely used in the design of irrigation systems to improve the efficiency of water application and the uniformity of water distribution. Film hole irrigation is somewhat similar to point-source irrigation, in that water infiltration occurs in the region directly around the film hole. Unlike other point source irrigation systems, in which the water is transported by tubes, in film mulch irrigation, water is applied to the top of the border or furrow and it flows above the applied film mulch to the end of the border or furrow under the influence of gravity, similar to surface irrigation. In the design of a surface irrigation system, a zero-inertia model and winSRFR software are the most useful tools. In the model, the infiltration and roughness are the key parameters to be determined. As with surface irrigation, the infiltration characteristics of film hole irrigation are fundamental for determining a field film hole irrigation scheme with high application efficiency and distribution uniformity. Therefore, studying a simple and easily-estimated infiltration model of film hole irrigation is essential. Numerous laboratory studies and some field studies on film hole infiltration have been conducted in the last two decades in China; such studies have mainly focused on wetting patterns, empirical modeling, and infiltration characteristics, including single and multipoint source infiltration [18,21,22]. However, the models of all of these studies have been empirical descriptions of some specific soils; therefore, a more universally applicable model must be developed.
Numerical simulation is often used in soil research. From a theoretical point of view, film hole irrigation involves three-dimensional point-source infiltration under a low-pressure water head (i.e., irrigation depth). Similar to subsurface drip irrigation, film hole infiltration is affected by many factors, such as soil texture, bulk density (γd), irrigation depth, film hole diameter (D), and hole spacing (i.e., distance between the centers of two neighboring holes). With the development of computer simulation techniques, numerical simulations based on the theory of unsaturated soil water movement are being increasingly used to study soil water infiltration. Several programs, such as HYDRUS and SWMS-2D, are often used to simulate soil water movement, and these have been effectively and accurately applied to subsurface drip irrigation [23,24,25,26,27,28] for predicting wetting patterns and infiltration characteristics, yielding more generally applicable results [29].
Therefore, the objectives of this study were: (1) to assess the feasibility of SWMS-2D for simulation of the cumulative infiltration of film hole irrigation through a laboratory experiment; (2) to investigate the effects of various influencing factors on cumulative infiltration in film hole irrigation and then to select the dominant factors; and (3) to propose and verify a simplified double-factor model that estimates the infiltration of film hole irrigation.

2. Materials and Methods

2.1. Laboratory Experiments

Given the realities that constrain film hole irrigation in the field, the film hole diameters and spacing are generally determined by the actual situation of field crops. The diameters of film holes and the spaces between them are typically 3–8 and 12–30 cm, respectively. The opening ratio, ρ, is defined as the ratio of the area of the open holes to the total area under the plastic mulching; in general, ρ is 2–5%; the irrigation amount is 225–450 m3·hm2; and the irrigation depth relative to the film mulch is kept constant within a range of 4–6 cm [30,31]. The soil is usually irrigated when the soil water content (SWC) is at 40–60% field capacity. Considering all of these irrigation variables, we designed eight treatments for our experiments (Table 1).
Figure 1 illustrates the laboratory setup for the experiments. Considering the symmetry characteristics of film hole infiltration, two soil bins, composed of 10-mm-thick transparent acrylic material, were designed; one of these bins measured 10 cm in length, 10 cm in width, and 60 cm in depth; the other was 15 cm in length, 15 cm in width, and 60 cm in depth. The bottom of each soil bin had numerous 2-mm parallel air vents for ventilation, and on both sides of the bins were side holes for soil sampling for SWC measurement; the diameter of each hole and the spaces between the holes were 1.5 and 5 cm, respectively. Before the soil was filled into the soil bin, transparent adhesive tape was placed on both sides of the soil bin to prevent the soil from falling out of the box. A Marriotte vessel was used to maintain a constant hydraulic head.
Silt loam and sandy loam (based on the USDA Soil Taxonomy System) from Yangling District, China, were used for the experiments. The particle size distributions of the silt loam were 13.52% 0–0.002 mm, 79.23% 0.002–0.02 mm, and 7.25% 0.02–2 mm, while those of the sandy loam were 3.10% 0–0.002 mm, 42.63% 0.002–0.02 mm, and 54.27% 0.02–2 mm. The soil samples were collected from the 20–60-cm depths of a field. Soil samples were air-dried, sieved through a 2-mm mesh, and then compacted into the soil bin at bulk densities of 1.30, 1.35, 1.40, and 1.45 g·cm−3 to simulate the in situ bulk density. Some water was added to the dry soil and mixed with the soil before filling the soil bin to reach the experimental desired initial SWC values of the eight treatments. The soil was loaded and compacted in the bin in 5-cm layers to obtain a homogeneous soil profile. To maintain zero evaporation, the soil surface was covered with a polyethylene sheet. Based on the design of the experiments, the irrigation volume of a single hole could be determined; accordingly, an infiltration depth of 4–6 cm was maintained by the Marriotte bottle. Cumulative infiltration was recorded during the infiltration. Finally, soil samples were collected from side holes and SWC was determined by recording the weight loss of the samples after oven drying at 105 °C for 24 h.

2.2. Numerical Model Establishment and Calibration

Because film hole irrigation can be simplified as a process of point-source water infiltration, water movement during infiltration can be considered as an axisymmetric three-dimensional infiltration process. The following partial differential equation based on the Richards equation governs water flow through variably unsaturated media:
θ t = r [ K ( φ ) φ r ] + K ( φ ) r φ r + z [ K ( φ ) ( φ z 1 ) ] ,
where θ is the soil water content (cm3·cm−3); φ is the potential head (cm); K ( φ ) is the hydraulic conductivity (cm·min−1); z is the depth from the soil surface (cm), measured with positive values in the downward direction; r is the radial coordinate (cm); and t is the time (min).
The infiltration space can be described as a three-dimensional axisymmetric domain, as shown in Figure 2a. With the infiltration process, water moves into the inner soil; the wetting pattern is shown in Figure 2b.
SWMS-2D under the axisymmetric three-dimensional flow model was used to solve the point source infiltration process of film hole irrigation [23]. The initial conditions of the system are as follows:
φ = φ 0 ,   0 r 1 / 2   S ,   0 z H ,   t = 0 ,
where φ 0 is the initial water potential before irrigation (cm−1); R and S are the radius and space of the film hole (cm), respectively; and H is the depth of simulation domain (cm).
The boundary AB is the surface through which water enters with a water head of h0 (Equation (3)). The boundary BC, covered with plastic mulch, with no infiltration or evaporation, is the zero-flux boundary, as shown in Equation (4). The boundaries AE and CD are partially symmetric boundaries without exchange of water flow in the r direction; they are also zero-flux boundaries, as shown in Equations (5) and (6). The boundary ED constitutes the bottom and is not affected by infiltration; therefore, it has a fixed initial soil water potential, as shown in Equation (7).
φ = h 0 ,   z = 0 ,   0 r R ,   t 0 ,
K ( φ ) [ φ z 1 ] = 0 ,   z = 0 ,   R r 1 / 2   S ,   t 0 ,
K ( φ ) φ r = 0 ,   0 z H ,   r = 0 ,   t 0 ,
K ( φ ) φ r = 0 ,   0 z H ,   r = 1 / 2   S ,   t 0 ,
φ = φ 0 ,   z = H ,   0 r 1 / 2   S ,   t 0 ,
The water retention curve and hydraulic conductivity curve are important soil characteristics [32]. A centrifugal machine (SCR-20) was used to obtain a water retention curve under rotational speeds of 900, 1700, 2200, 2800, 3100, 5300, 6900, and 8100 rpm; each speed was maintained for 1 h. Then, SWC was measured for each rotational speed, whereupon the capillary head was calculated from Equation (8) [33].
h = 1.118 × 10 5 × R 0 × ( rpm ) 2
where R0 is the radial distance to the midpoint of the soil sample (cm).
The van Genuchten–Mualem model was used to describe the curves of soil water retention, θ(φ), and of hydraulic conductivity, K(θ):
θ ( φ ) = θ r + θ s θ r ( 1 + | α φ | n ) m ,
K ( θ ) = K s S e l [ 1 ( 1 S e 1 m ) m ] 2 ,
S e = θ θ r θ s θ r ,
D ( θ ) = K ( θ ) d φ d θ
where θr and θs are the residual and saturated water contents (cm3·cm−3), respectively; K(θ) is the saturated hydraulic conductivity (cm·min−1); D(θ) is the soil diffusivity, which is obtained by horizontal column infiltration method (cm2·min−1); α is an empirical parameter that is inversely related to the air-entry pressure value (cm−1); n is an empirical parameter related to the pore-size distribution; l is an empirical shape parameter; m = 1 − 1/n; and Se is the effective saturation.
The experimental data of soil water retention and diffusivity were fitted to the van Genuchten–Mualem model by using the RETC code [34,35]. The parameters of the van Genuchten–Mualem model are listed in Table 2. A comparison of fitted soil water retention and diffusivity curves by model and measured data is shown in Figure 3.
The van Genuchten–Mualem model parameters for 11 selected types of soil, including silt loam, sandy loam, sand, sandy loam, loam, loamy sand, silt, silt loam, sandy clay loam, clay loam, and silty clay loam, were obtained using RETC software [35] and are listed in Table 2.

2.3. Simplified Kostiakov Model for Film Hole Infiltration

The cumulative infiltration and duration were described using the Kostiakov model given as Equation (13). In this study, six influencing factors for film hole irrigation were analyzed. Moreover, dominant factors were selected, and a simplified equation based on the Kostiakov model for film hole infiltration would be further proposed.
I = k t α
where I is the cumulative infiltration (mm), t is the infiltration duration (min), k is an infiltration coefficient (mm min−1), and α is the infiltration index.

2.4. Error Analysis

Four indicators, namely mean absolute error (MAE), root mean squared error (RMSE), percent bias (PBIAS), and Nash–Sutcliffe coefficient (NS), were selected to make error analyses between the measured and simulated values of SWC and cumulative infiltration.
MAE = 1 n | Y i s i m Y i o b s | n
RMSE = 1 n 1 i = 1 n ( Y i s i m Y i o b s ) 2 , for   i = 1 , 2 n
PBIAS = i = 1 n ( Y i s i m Y i o b s ) i = 1 n Y i o b s
NS = 1 i = 1 n ( Y i s i m Y i o b s ) 2 i = 1 n ( Y i o b s Y ¯ i o b s ) 2
where n is the total number of data points in each case,. Y i sim is the ith simulated datum, and Y i o b s is the ith measured datum. The MAE can potentially identify the presence of bias. The RMSE provides an overall measure of the degree to which the data differ from the model predictions, whereas the PBIAS is the deviation of data being evaluated, expressed as a percentage. If PBIAS within ±10%, the PBIAS values are considered to be within a very accurate range. A value of 1 for the NS means the simulated value is as accurate as the measured value [36].

3. Results and Discussion

3.1. Application of SWMS-2D to Film Hole Irrigation and Comparison with the Experimental Results

Figure 4 shows the measured and simulated SWC distribution data for the eight treatments in Table 1. The measured SWC is represented by black dots, whereas simulation results obtained through SWMS-2D are shown as contour lines. The MAE, RMSE, PBIAS, and NS values for measured and simulated SWC are presented in Table 3; the results indicated the SWMS-2D model to be suitable for describing soil water distributions resulting from film hole irrigation, because it had extremely small MAE, RMSE, and PBIAS values (1.091–2.600 cm3·cm−3, 0.013–0.034 cm3·cm−3, and 0.461–3.648%, respectively), and a large NS (very close to 1.0).
Figure 5 compares simulated and measured cumulative infiltration values for different treatments. The simulated values agreed reasonably well with the measured values, with small PBIAS values from −0.197% to 0.914% and large NS values from 0.989 to 0.999, which indicated that the SWMS-2D model can effectively simulate cumulative infiltration resulting from film hole irrigation.

3.2. Different Factors Affecting Cumulative Infiltration of Film Hole Irrigation

3.2.1. Effect of Initial SWC on Cumulative Infiltration

Simulations of three soil types, namely silt loam (γd = 1.30 g·cm−3), sandy loam (γd = 1.35 g·cm−3), and loamy sand (simulated on the basis of information obtained from RETC software), were executed at different initial SWC levels with a film hole diameter of 5 cm, spacing of 20 cm, irrigation depth of 6 cm, and irrigation amount of 450 m3·hm−2. The field capacities of the silt loam, sandy loam, and loamy sand were 0.243, 0.126, and 0.160 cm3·cm−3, respectively. The cumulative infiltration curves of the film holes at different initial SWC levels are shown in Figure 6; the initial SWC had little effect on the cumulative infiltration dynamic of a single film hole. As the SWC increased, the water potential gradient slightly decreased, leading to a slight decrease in cumulative infiltration. Therefore, the impacts of initial SWC could be ignored in film hole irrigation research.

3.2.2. Effect of Irrigation Depth on Cumulative Infiltration

Figure 7 illustrated the effects of different irrigation depths in the silt loam, sandy loam, and loamy sand. All simulations were conducted at a film hole diameter of 5 cm, film hole spacing of 20 cm, initial SWC of 50% field water capacity, and irrigation amount of 450 m3·hm−2. The irrigation depth had little effect on the film hole infiltration characteristics (Figure 7). As the irrigation depth increased, the film hole infiltration rate slightly increased. However, because the value of gravity potential caused by irrigation depth was so small, relative to matrix potential that changes of gravity potential at different irrigation depths was smaller, that it can be ignored. Accordingly, because of the slight changes in film hole infiltration at different irrigation depths, the infiltration was disregarded in later investigations.

3.2.3. Effect of Layered-Soil Depth on Cumulative Infiltration

To consider the actual soil conditions in the field, and to analyze the effect of upper soil depth on film hole infiltration, this study conducted simulations of six cases of layered soil. Different soil types and bulk densities were obtained at various upper soil depths. Homogeneous soils and upper soils with 15-, 20-, and 30-cm-deep soil layers were simulated for all cases. Case 1 involved sandy loam (1.35 g·cm−3) on top with silt loam (1.30 g·cm−3) underneath. Case 2 involved silt loam (1.40 g·cm−3) on top with sandy loam (1.45 g·cm−3) underneath. Case 3 involved loamy sand on top with loam underneath. Case 4 involved loam on top with loamy sand underneath. Cases 5 and 6 involved combinations of the same soil types but with different bulk densities: Case 5 considered silt loam soil with bulk densities of 1.40 and 1.30 g·cm−3, whereas Case 6 considered sandy loam with bulk densities of 1.35 and 1.45 g·cm−3. As illustrated in Figure 8, the hydraulic characteristics of the soil at the lower layers (at depths beyond 15 cm) had no significant effect on cumulative infiltration, which was predominantly affected by the characteristics of the soil at the upper layers. In other words, the upper soil characteristics predominantly affected film hole infiltration characteristics under film hole irrigation with the normal amount of water. Thus, the soil at depths of 0–15 cm has the dominant influence on the cumulative infiltration of layered soil.

3.2.4. Effect of Film Hole Diameter and Opening Ratio on Cumulative Infiltration

Silt loam (γd = 1.30 g·cm−3), sandy loam (γd = 1.45 g·cm−3), loam, and sand were selected to simulate cumulative infiltration of film hole irrigation under different D and ρ values. The simulations were conducted at a film hole spacing of 20 cm, an irrigation water depth of 6 cm, an initial SWC of 50% field capacity, and an irrigation amount of 450 m3·hm−2. As shown in Figure 9, for all treatments, the cumulative infiltration increased as ρ increased, but decreased as D increased under the same ρ; in addition, sand had the highest infiltration rate, followed by loam and sandy loam. From the above analysis, the effects of D and ρ should be taken into account in film hole irrigation research.

3.3. Establishment and Universality of a Simplified Model

3.3.1. Establishment of a Simplified Model

This study selected 11 soil types, including silt loam (γd = 1.30 and 1.40 g·cm−3) and sandy loam (γd = 1.35 and 1.45 g·cm−3) as well as other sandy, loam, and silt soils, to analyze the characteristics of film hole infiltration. Cumulative infiltration was simulated at a film hole diameter of 5 cm, spacing of 20 cm, irrigation water depth of 6 cm, and irrigation amount of 450 m3·hm−2.
Equation (13) was used to fit the relationship between k and α, listed in Table 4, for different soils. The correlation coefficients (R2) for all soils were all larger than 0.99, indicating that the Kostiakov model can adequately describe the relationship between cumulative infiltration and duration.
Based on the preceding analysis, with the same soil texture and the same bulk density as actual field film hole irrigation, the initial SWC and irrigation depth weakly affected cumulative infiltration. Therefore, D and ρ were viewed as the two main influencing factors in model establishment for a given soil. The values of k increased with ρ, but decreased with increased D under the same ρ. The values of α decreased with ρ, and had no significant relationship with D. After extensive verification, Equations (18) and (19) were proposed based on the relationships among k, α, ρ, and D.
k = A ρ m D n ,
α = B η ρ ,
where k is the infiltration coefficient (mm min−1); ρ is the opening ratio (%); D is the film hole diameter (cm); and A, m, n, B, and η are the fitting parameters.
Based on the k and α data for all soil samples listed in Table 4, the model parameters A, B, m, n, and η were fitted using MATLAB, and the values of m, n, and η were 1.131, 0.698, and 0.997 (where 0.997 can be considered as effectively equal to 1.0), respectively; and A and B varied with soil texture and bulk density. Therefore, Equation (13) was reconstructed using a combination of D and ρ as follows:
I = A ( ρ 1.131 D 0.698 ) t B ρ

3.3.2. Evaluation of the Simplified Model Universality

Table 4 compares parameters k and α for different soil texture and bulk densities; their estimated values based on Equations (18) and (19) are shown in Figure 10. The simulated values of k and α were in good agreement with almost all estimated values, indicating that the simplified model, Equation (20), is suitable for infiltration estimation for various soils with different D and ρ values. In addition, α was observed to be numerically stable at approximately 0.9; however, this requires further analysis.
The experiment results, including four soil types located in different regions on the Loess Plateau, China, which were loam from Wugong (γd = 1.30 g·cm−3), silt loam from Wei River Bench (γd = 1.45 g·cm−3), loam from Luochuan (γd = 1.40 g·cm−3) and sandy loam from Ansai (γd = 1.35 g·cm−3) [37], were used to verify the universality of the simplified model. The infiltration models were obtained based on Equation (20). As presented in Figure 11, the calculated values obtained from Equations (21)–(24) were well matched with the experimental values, signifying that the simplified model of Equation (20) can effectively describe the characteristics of film hole infiltration.
I = 42.94 ( ρ 1.131 / D 0.698 ) t 0.8635 ρ
I = 50.92 ( ρ 1.131 / D 0.698 ) t 0.8809 ρ
I = 48.24 ( ρ 1.131 / D 0.698 ) t 0.8677 ρ
I = 65.75 ( ρ 1.131 / D 0.698 ) t 0.9362 ρ
Three published data sets were then next selected for a more thorough evaluation of the simplified model established in this study. The first data set was obtained from Li [21], and the experimental soil was sandy loam with dry bulk density of 1.30 g·cm−3, which was taken from the Yulin region, located in the North Loess Plateau, China. The cumulative infiltration was analyzed with four different combinations of the ρ and D values. The D of the experiments was 2–5 cm. The second data set came from Hu [38], whose experimental soil was silt loam with bulk density of 1.35 g·cm−3, selected from Handan City, located in the North China Plain. The D of the experiments was 3–5 cm and the spacing between holes was 17.5 cm. The third data set obtained by Wu [18] with sand as the experimental soil and bulk density of 1.52 g·cm−3, which taken from the Yulin region. The film hole diameters in the experiment were 3–6 cm, and the spacing between holes was 12 cm. The details from these experiments have been provided in their respective publications.
A comparison between the calculated and measured values of the cumulative infiltration of a single film hole within the controlling area under different D and ρ values is illustrated in Figure 12. Equations (25)–(27) are the simplified model equations for each soil type motioned above; the MAE, RMSE, PBIAS, and NS values for measured and calculated cases are presented in Table 5. The MAE, RMSE, and PBIAS values ranged from 0.141 to 2.299 mm, 0.177 to 2.722 mm, and −2.131% to 1.479%, respectively; meanwhile, the NS values were very close to 1.0. Notably, all results were in good agreement, indicating that the model can effectively describe the characteristics of film hole infiltration.
I = 62.64 ( ρ 1.131 / D 0.698 ) t 0.7415 ρ
I = 51.05 ( ρ 1.131 / D 0.698 ) t 0.9217 ρ
I = 480.53 ( ρ 1.131 / D 0.698 ) t 0.9245 ρ

4. Conclusions

Film hole infiltration is a low-cost and high-efficiency irrigation method that is widely used in China. In this study, the SWMS-2D-simulated wetting patterns and cumulative infiltration observed for film hole irrigation were in good agreement with experimental observations for silt loam and sandy loam at different soil bulk densities. Therefore, the SWMS-2D model is suitable for simulating film hole irrigation. In addition, the initial SWC and irrigation depth had little effect on cumulative infiltration during film hole irrigation, whereas the D and ρ significantly affected the cumulative infiltration. Furthermore, the observed cumulative infiltration increased with ρ, but it decreased with increased D with the same ρ. Finally, we proposed a simplified model for film hole irrigation, and this model has extensive applicability. In conclusion, cumulative infiltration can be adequately described with D and ρ for a soil with certain texture and bulk density.

Acknowledgments

This research was jointly supported by grant of National Natural Science Foundation of China (No. 51279167 and No. 51579205), Special Fund for Agro-scientific Research in the Public Interest (No. 201503124) and Specialized Research Fund for the Doctoral Program of Higher Education of China (No. 20120204110023). A special and sincere thanks to the Key Laboratory of Agricultural Soil and Water Engineering in Arid and Semiarid Areas, Ministry of Education for the support for the experiments.

Author Contributions

Xiao-Yi Ma and Yan-Wei Fan conceived and designed the experiments; Yi-Bo Li and Yan-Wei Fan performed the experiments; Yi-Bo Li and Ye Liu analyzed the data; Xiao-Yi Ma contributed reagents/materials/analysis tools; Yi-Bo Li wrote the paper; and Xiao-Yi Ma revised the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Matti, K.; Ward, P.J.; Moel, H.D.; Varis, O. Is physical water scarcity a new phenomenon? Global assessment of water shortage over the last two millennia. Environ. Res. Lett. 2010, 5. [Google Scholar] [CrossRef] [Green Version]
  2. Wallace, J.S. Increasing agricultural water use efficiency to meet future food production. Agric. Ecosyst. Environ. 2001, 82, 105–119. [Google Scholar] [CrossRef]
  3. Swatuk, L.; Mcmorris, M.; Leung, C.; Zu, Y. Seeing “invisible water”: Challenging conceptions of water for agriculture, food and human security. Can. J. Dev. Stud. 2015, 36, 24–37. [Google Scholar] [CrossRef]
  4. Deng, X.P.; Shan, L.; Zhang, H.; Turner, N.C. Improving agricultural water use efficiency in arid and semiarid areas of China. Agric. Water Manag. 2006, 80, 23–40. [Google Scholar] [CrossRef]
  5. Li, X.Y.; Gong, J.D.; Wei, X.H. In-situ rainwater harvesting and gravel mulch combination for corn production in the dry semi-arid region of China. J. Arid Environ. 2000, 46, 371–382. [Google Scholar] [CrossRef]
  6. Camp, C.R. Subsurface drip irrigation: A review. Trans. ASAE 1998, 41, 1353–1367. [Google Scholar] [CrossRef]
  7. Tarara, J.M. Microclimate modification with plastic mulch. Hortscience 2000, 35, 169–180. [Google Scholar]
  8. McKenzie, C.; Duncan, L.W. Landscape fabric as a physical barrier to neonate diaprepes abbreviatus (Coleoptera: Curculionidae). FL Entomol. 2001, 84, 721–722. [Google Scholar] [CrossRef]
  9. Chalker-Scott, L. Impact of mulches on landscape plants and the Environment—A review. J. Environ. Hort. 2007, 25, 239–249. [Google Scholar]
  10. Espí, E.; Salmerón, A.; Fontecha, A.; García, Y.; Real, A.I. Plastic films for agricultural applications. J. Plast. Film Sheet 2006, 22, 85–102. [Google Scholar] [CrossRef]
  11. Lamont, W.J. Plastic mulches for the production of vegetable crops. Horttechnology 1993, 3, 35–39. [Google Scholar]
  12. Scarasciamugnozza, G.; Sica, C.; Russo, G. Plastic materials in European agriculture: Actual use and perspectives. J. Agric. Eng. 2011, 42, 15–28. [Google Scholar] [CrossRef]
  13. Steinmetz, Z.; Wollmann, C.; Schaefer, M.; Buchmann, C.; David, J.; Tröger, J.; Muñoz, K.; Frör, O.; Schaumann, G.E. Plastic mulching in agriculture. Trading short-term agronomic benefits for long-term soil degradation? Sci. Total Environ. 2016, 550, 690–705. [Google Scholar] [CrossRef] [PubMed]
  14. Kader, M.A.; Senge, M.; Mojid, M.A.; Ito, K. Recent advances in mulching materials and methods for modifying soil environment. Soil Till. Res. 2017, 168, 155–166. [Google Scholar] [CrossRef]
  15. National Bureau of Statistics of China. China Statistical Year Book 2014; National Bureau of Statistics of China: Beijing, China, 2014.
  16. Xu, S. Fim hole irrigation technology. Xinjiang Water Res. 1994, 83, 22–26. (In Chinese) [Google Scholar]
  17. Saeed, M.; Mahmood, S. Application of film hole irrigation on borders for water saving and sunflower production. Arab J. Sci. Eng. 2013, 38, 1347–1358. [Google Scholar] [CrossRef]
  18. Wu, J.; Fei, L.; Wang, W. Study on the infiltration characteristics of single filmed hole and its mathematical model under filmed hole irrigation. Adv. Water Sci. 2001, 9, 307–311. (In Chinese) [Google Scholar]
  19. Li, Y.; Tian, J. Research advances irrigating technique of film hole irrigation. J. Ningxia Agric. Coll. 2003, 24, 96–100. (In Chinese) [Google Scholar]
  20. Bautista, E.; Clemmens, A.J.; Strelkoff, T.S.; Schlegel, J. Modern analysis of surface irrigation systems with WinSRFR. Agric. Water Manag. 2009, 96, 1146–1154. [Google Scholar] [CrossRef]
  21. Li, F.; Zhang, X.; Fei, L. Experimental study on influential factors in film hole bilateral interference infiltration. J. Irrig. Drain. 2002, 22, 44–48. (In Chinese) [Google Scholar]
  22. Dong, Y.; Wang, B.; Jia, L.; Fei, L. Effects of different irrigation treatments on maize water consumption, growth and yield under film hole irrigation. Appl. Mech. Mater. 2014, 501–504, 1986–1992. [Google Scholar] [CrossRef]
  23. Šimůnek, J.; Huang, K.; Genuchten, M.T.V. The SWMS-2D Code for Simulating Water Flow and Solute Transport in Two-Dimensional Variably-Saturated Media; Version 1.0; US Salinity Laboratory, US Department of Agriculture, Agricultural Research Service: Riverside, CA, USA, 1995.
  24. Skaggs, T.H.; Trout, T.J.; Simunek, J.; Shouse, P.J. Comparison of HYDRUS-2D simulations of drip irrigation with experimental observations. J. Irrig. Drain. Eng. 2014, 130, 304–310. [Google Scholar] [CrossRef]
  25. Zhou, Q.Y.; Kang, S.Z.; Zhang, L. Comparison of APRI and Hydrus-2D models to simulate soil water dynamics in a vineyard under alternate partial root zone drip irrigation. Plant Soil 2007, 291, 211–223. [Google Scholar] [CrossRef]
  26. Lazarovitch, N.; Poulton, M.; Furman, A.; Warrick, A.W. Water distribution under trickle irrigation predicted using artificial neural networks. J. Eng. Math. 2009, 64, 207–218. [Google Scholar] [CrossRef]
  27. Mguidiche, A.; Provenzano, G.; Douh, B.; Khila, S.; Rallo, G.; Boujelben, A. Assessing hydrus-2D to simulate soil water content (SWC) and salt accumulation under an SDI system: Application to a potato crop in a semi-arid area of Central Tunisia. Irrig. Drain. 2015, 64, 263–274. [Google Scholar] [CrossRef]
  28. Hou, L.; Zhou, Y.; Bao, H.; Wenninger, J. Simulation of maize (Zea mays L.) water use with the HYDRUS-1D model in the semi-arid Hailiutu River catchment, Northwest China. Hydrol. Sci. J. 2016, 62, 1–11. [Google Scholar] [CrossRef]
  29. Kandelous, M.M.; Šimůnek, J. Numerical simulations of water movement in a subsurface drip irrigation system under field and laboratory conditions using HYDRUS-2D. Agric. Water Manag. 2010, 97, 1070–1076. [Google Scholar] [CrossRef]
  30. Jiang, L. Technique element of film hole infiltration on cotton. China Rural Water Hydropower 1995, 5, 3–6. (In Chinese) [Google Scholar]
  31. Zhang, G.; Zhu, L. Study on irrigating technique of film hole irrigation in Weibei dry plateau. Res. Soil Water Conserv. 1999, 6, 60–63. (In Chinese) [Google Scholar]
  32. Xing, X.; Kang, D.; Ma, X. Differences in loam water retention and shrinkage behavior: Effects of various types and concentrations of salt ions. Soil Till. Res. 2017, 167, 61–72. [Google Scholar] [CrossRef]
  33. Yao, W.W.; Ma, X.Y.; Li, J.; Parkes, M. Simulation of point source wetting pattern of subsurface drip irrigation. Irrig. Sci. 2011, 29, 331–339. [Google Scholar] [CrossRef]
  34. Maulem, Y. A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Resour. Res. 1976, 12, 513–522. [Google Scholar] [CrossRef]
  35. Van Genuchten, M.T.; Leij, F.J.; Yates, S.R. The RETC Code for Quantifying the Hydraulic Functions of Unsaturated Soils; Research Report No. EPA/600/2-91/065; US Salinity Laboratory, US Department of Agriculture, Agricultural Research Service: Riverside, CA, USA, 1991. Available online: http://www.pc-progress.com/documents/programs/retc.pdf (accessed on 1 December 1991).
  36. Moriasi, D.; Arnold, J.; Van Liew, M.W.; Veith, T.L. Model evaluation guidelines for systematic quantification of accuracy in watershed simulations. Trans. ASABE 2007, 50, 885–900. [Google Scholar] [CrossRef]
  37. Ma, X.; Fan, Y.; Wang, S.; Sun, X.; Wang, B. Simplified Infiltration Model of Film Hole Irrigation and Its Validation. Trans. Chin. Soc. Agric. Mach. 2009, 40, 67–73. (In Chinese) [Google Scholar]
  38. Hu, H. Soil Water and Nitrogen Transport Characteristic Experiment and Numerical Simulation of Film Hole Irrigation under Facilities. Ph.D. Dissertation, Wuhan University, Wuhan, China, 2010. (In Chinese). [Google Scholar]
Figure 1. Simplified experimental setup (h0 represents irrigation depth).
Figure 1. Simplified experimental setup (h0 represents irrigation depth).
Water 09 00543 g001
Figure 2. Simplified graph of film hole irrigation model ((a) represents the three-dimensional axisymmetric domain; (b) represents the wetting pattern during the infiltration process).
Figure 2. Simplified graph of film hole irrigation model ((a) represents the three-dimensional axisymmetric domain; (b) represents the wetting pattern during the infiltration process).
Water 09 00543 g002
Figure 3. Measured (dots) and fitted (line) water retention and diffusivity curves of silt loam and sandy loam using the van Genuchten–Mualem Model.
Figure 3. Measured (dots) and fitted (line) water retention and diffusivity curves of silt loam and sandy loam using the van Genuchten–Mualem Model.
Water 09 00543 g003
Figure 4. Comparison of simulated and measured values of soil wetting patterns (Measured data are represented by black dots, whereas simulated values are shown using water content as contour lines).
Figure 4. Comparison of simulated and measured values of soil wetting patterns (Measured data are represented by black dots, whereas simulated values are shown using water content as contour lines).
Water 09 00543 g004aWater 09 00543 g004b
Figure 5. Comparison of simulated and measured values of cumulative infiltration under different treatments (cumulative infiltration represents the irrigation volume of the controlling area of a single hole in the experimental setup (see Figure 1); measured values are represented by dots, whereas simulated values are represented by lines).
Figure 5. Comparison of simulated and measured values of cumulative infiltration under different treatments (cumulative infiltration represents the irrigation volume of the controlling area of a single hole in the experimental setup (see Figure 1); measured values are represented by dots, whereas simulated values are represented by lines).
Water 09 00543 g005
Figure 6. Effect of initial SWC on film hole infiltration: (a) silt loam γd = 1.30 g·cm−3; (b) sandy Loam γd = 1.35 g·cm−3; and (c) loamy sand (the percentage is the percentage of field capacity).
Figure 6. Effect of initial SWC on film hole infiltration: (a) silt loam γd = 1.30 g·cm−3; (b) sandy Loam γd = 1.35 g·cm−3; and (c) loamy sand (the percentage is the percentage of field capacity).
Water 09 00543 g006
Figure 7. Infiltration pattern curves under different irrigated water depths: (a) silt loam, γd = 1.30 g·cm−3; (b) sandy loam, γd = 1.35 g·cm−3; and (c) loamy sand.
Figure 7. Infiltration pattern curves under different irrigated water depths: (a) silt loam, γd = 1.30 g·cm−3; (b) sandy loam, γd = 1.35 g·cm−3; and (c) loamy sand.
Water 09 00543 g007
Figure 8. Effect of layered soil depth on film hole infiltration. (a) Experimental soils: 1–4 sandy loam (γd = 1.35 g·cm−3) + silt loam (γd = 1.30 g·cm−3); 5–8 silt loam (γd = 1.40 g·cm−3) + sandy loam (γd = 1.45 g·cm−3); (b) The soils selected from RETC software: 1–4 loamy sand + loam, 5–8 loam + loamy sand; (c) Double bulk densities: 1–4 silt loam, 1.40 g·cm−3 + 1.30 g·cm−3; 5–8 sandy loam, 1.35 g·cm−3 + 1.45 g·cm−3. For Nos. 1–8: 1 and 5—upper depth 15 cm; 2 and 6—upper depth 20 cm; 3 and 7—upper depth 30 cm; 4 and 8—homogeneity.
Figure 8. Effect of layered soil depth on film hole infiltration. (a) Experimental soils: 1–4 sandy loam (γd = 1.35 g·cm−3) + silt loam (γd = 1.30 g·cm−3); 5–8 silt loam (γd = 1.40 g·cm−3) + sandy loam (γd = 1.45 g·cm−3); (b) The soils selected from RETC software: 1–4 loamy sand + loam, 5–8 loam + loamy sand; (c) Double bulk densities: 1–4 silt loam, 1.40 g·cm−3 + 1.30 g·cm−3; 5–8 sandy loam, 1.35 g·cm−3 + 1.45 g·cm−3. For Nos. 1–8: 1 and 5—upper depth 15 cm; 2 and 6—upper depth 20 cm; 3 and 7—upper depth 30 cm; 4 and 8—homogeneity.
Water 09 00543 g008
Figure 9. Infiltration pattern curve under different film hole spacing: (a1,a2) silt loam, γd = 1.30 g·cm−3; (b1,b2) sandy loam, γd = 1.45 g·cm−3; (c1,c2) loam; and (d1,d2) sand. In the legend, the percentage values are equal to ρ, and the following values are equal to D.
Figure 9. Infiltration pattern curve under different film hole spacing: (a1,a2) silt loam, γd = 1.30 g·cm−3; (b1,b2) sandy loam, γd = 1.45 g·cm−3; (c1,c2) loam; and (d1,d2) sand. In the legend, the percentage values are equal to ρ, and the following values are equal to D.
Water 09 00543 g009
Figure 10. Comparison of simulated and estimated k and α. Nos. 1–13 indicate 13 types of tested soil form experiments performed using RETC software: 1 = silt loam (1.30 g·cm−3); 2 = silt loam (1.30 g·cm−3); 3 = sandy loam (1.35 g·cm−3); 4 = sandy loam (1.45 g·cm−3); 5 = loam; 6 = sand; 7 = sandy loam; 8 = loamy sand; 9 = silt; 10 = silt loam; 11 = sandy clay loam; 12 = clay loam; and 13 = silt clay loam.
Figure 10. Comparison of simulated and estimated k and α. Nos. 1–13 indicate 13 types of tested soil form experiments performed using RETC software: 1 = silt loam (1.30 g·cm−3); 2 = silt loam (1.30 g·cm−3); 3 = sandy loam (1.35 g·cm−3); 4 = sandy loam (1.45 g·cm−3); 5 = loam; 6 = sand; 7 = sandy loam; 8 = loamy sand; 9 = silt; 10 = silt loam; 11 = sandy clay loam; 12 = clay loam; and 13 = silt clay loam.
Water 09 00543 g010
Figure 11. Comparison of four soil types between calculated values and measured values of cumulative infiltration.
Figure 11. Comparison of four soil types between calculated values and measured values of cumulative infiltration.
Water 09 00543 g011
Figure 12. Comparison of three infiltration models (cited from published studies) between calculated and measured values of cumulative infiltration.
Figure 12. Comparison of three infiltration models (cited from published studies) between calculated and measured values of cumulative infiltration.
Water 09 00543 g012
Table 1. Test scheme for film hole irrigation.
Table 1. Test scheme for film hole irrigation.
TreatmentSoil Type *Bulk Density (g/cm3)Distance between Film Holes (cm)Hole Diameter (cm)Initial Water Content ** (%)Irrigation Depth (cm)Irrigation Amount (m3/hm2)Irrigation Volume *** (mL)
1Silt loam1.30204404225900
21.402046064501800
31.303066062252025
41.403064044504050
5Sandy loam1.353046044504050
61.453044062252025
71.352064064501800
81.45206604225900
Notes: * Based on the USDA Soil Taxonomy System; ** Initial water content is the percentage of field capacity; *** Irrigation volume represents the irrigation volume of water supplied within the controlling area of a single hole in the experimental setup.
Table 2. Parameters for the van Genuchten–Mualem model of experimental soils.
Table 2. Parameters for the van Genuchten–Mualem model of experimental soils.
Soil TypeBulk Density (g/cm3)θr (cm3/cm3)θs (cm3/cm3)α (1/cm)nlKs (cm/min)
Silt loam1.300.1130.4920.0141.720.50.0208
1.400.1230.4560.0171.550.50.0154
Sandy loam1.350.0570.380.0251.760.50.0409
1.450.0490.3620.0241.690.50.0257
Sand-0.0450.4300.1452.680.50.4950
Sandy loam-0.0650.4100.0751.890.50.0737
Loam-0.0780.4300.0361.560.50.0173
Loamy sand-0.0570.0400.1242.280.50.2432
Silt-0.0340.4600.0161.370.50.0042
Silt loam-0.0670.4500.0201.410.50.0075
Sandy clay loam-0.1000.3900.0591.480.50.0218
Clay loam-0.0950.4100.0191.310.50.0043
Silty clay loam-0.0890.4300.0011.230.50.0012
Table 3. Correlation between simulated and measured SWC levels.
Table 3. Correlation between simulated and measured SWC levels.
TreatmentMAE (cm3/cm3)RMSE (cm3/cm3)PBIAS (%)NS
10.0130.0162.9470.979
20.0220.0282.7110.998
30.0260.0343.6480.997
40.0110.0131.9620.995
50.0200.0233.5080.968
60.0110.0131.4380.998
70.0130.0160.4610.998
80.0260.0282.7620.999
Table 4. Fitted infiltration parameter values.
Table 4. Fitted infiltration parameter values.
Opening Ratio (%)Film Hole Diameter (cm)Silt LoamSilt LoamSandy LoamSandy LoamLoamSandSandy Loam
(γd = 1.30 g/cm3)(γd = 1.40 g/cm3)(γd = 1.35 g/cm3)(γd = 1.45 g/cm3)
kαkαkαkαkαkαkα
2.1830.25530.8770.16360.8580.29190.9010.1990.890.12470.8661.59570.960.31660.923
40.21040.8710.1350.8530.23250.9020.1610.8880.1050.8611.2810.9550.26580.912
50.18110.8660.11550.8520.19870.9010.1380.8860.09280.8541.08950.950.23270.904
Mean0.21560.8710.1380.8550.2410.9010.1660.8880.10750.861.32210.9550.27170.913
3.1420.52420.8580.33540.8390.61070.8750.4140.8650.25030.8493.27320.9540.61460.923
30.38250.8610.24660.8410.43450.8850.2980.8730.1860.8522.32250.9530.46270.915
40.31640.8560.20520.8360.34690.8880.2430.8730.15820.8461.86650.9480.39120.903
Mean0.40770.8580.26240.8390.4640.8830.3180.870.19820.8492.48740.9510.48950.913
4.9130.62650.8350.40620.8150.70980.8570.4480.8460.3030.833.69610.9370.74050.899
40.51520.8340.3350.8150.57320.8630.3970.8510.25580.8282.95350.9370.62320.891
50.44950.830.2940.8130.49590.8620.3460.8490.23210.822.52650.9320.55430.881
Mean0.53040.8330.34510.8140.5930.8610.3970.8490.26360.8263.05870.9350.63930.89
7.0760.61290.8060.40410.7880.66520.840.4670.8270.32080.7963.27530.9160.74750.859
13.651.4080.8550.93450.7371.54220.7761.0870.7680.71970.7597.37730.8631.65530.828
Opening Ratio (%)Film Hole Diameter (cm)Loamy SandSiltSilt LoamSandy Clay LoamClay LoamSilt Clay Loam
kαkαkαkαkαkα
2.1830.77570.9670.03110.9050.04450.9230.07890.9460.02360.9220.09890.805
40.60490.9690.02580.8980.04260.8960.06190.9460.02470.8860.10710.752
50.49960.9730.01940.9150.03140.9170.04910.9560.01520.930.04410.861
Mean0.62680.9690.02550.9060.03950.9120.06330.950.02120.9130.08330.806
3.1421.67260.9440.10240.8270.13840.8470.17310.9240.06830.8610.35560.686
31.22720.9520.06050.8640.08770.8780.12770.9310.04740.8790.20410.737
40.90040.960.03580.9030.05550.910.09410.9370.03290.8960.11710.791
Mean1.26680.9520.06620.8650.09390.8790.13160.9310.04950.8790.22560.738
4.9131.89150.9320.10620.8330.150.8480.21230.9040.08010.8510.38780.678
41.47470.9410.07090.8640.11040.8670.16620.9110.06340.860.24640.726
51.22570.9470.05810.8720.09140.8740.13550.9190.05180.8680.17810.754
Mean1.53060.940.07840.8560.11730.8630.17130.9120.06510.860.27080.719
7.0761.61890.9290.08560.8470.13360.8480.18580.9020.07620.8430.25810.720
13.654.06320.8460.32170.7380.45360.7460.50990.8380.24120.7630.93970.580
Table 5. Correlation between measured and calculated values of cumulative infiltration.
Table 5. Correlation between measured and calculated values of cumulative infiltration.
SoilMAE (mm)RMSE (mm)PBIAS (%)NS
Loam from Wugong0.1900.190−2.1310.995
Silt Loam from Weihe Bench0.3910.552−1.0820.997
Loam from Luochuan2.2992.722−1.7560.989
Sandy Loam from Ansai0.5080.830−1.3750.985
Sandy Loam from Yulin0.4310.5661.4790.999
Silt Loam from Handan0.1410.177−1.5690.982
Sand from Yulin0.3900.466−1.1700.998

Share and Cite

MDPI and ACS Style

Li, Y.-B.; Fan, Y.-W.; Liu, Y.; Ma, X.-Y. Influencing Factors and Simplified Model of Film Hole Irrigation. Water 2017, 9, 543. https://doi.org/10.3390/w9070543

AMA Style

Li Y-B, Fan Y-W, Liu Y, Ma X-Y. Influencing Factors and Simplified Model of Film Hole Irrigation. Water. 2017; 9(7):543. https://doi.org/10.3390/w9070543

Chicago/Turabian Style

Li, Yi-Bo, Yan-Wei Fan, Ye Liu, and Xiao-Yi Ma. 2017. "Influencing Factors and Simplified Model of Film Hole Irrigation" Water 9, no. 7: 543. https://doi.org/10.3390/w9070543

APA Style

Li, Y. -B., Fan, Y. -W., Liu, Y., & Ma, X. -Y. (2017). Influencing Factors and Simplified Model of Film Hole Irrigation. Water, 9(7), 543. https://doi.org/10.3390/w9070543

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop