Next Article in Journal
Detailed Trend Analysis of Extreme Climate Indices in the Upper Geum River Basin
Next Article in Special Issue
Searching for Sustainable-Irrigation Issues of Clementine Orchards in the Syrian Akkar Plain: Effects of Irrigation Method and Canopy Size on Crop Coefficients, Transpiration, and Water Use with SIMDualKc Model
Previous Article in Journal
Development of a New 8-Parameter Muskingum Flood Routing Model with Modified Inflows
Previous Article in Special Issue
Estimation of Evapotranspiration and Crop Coefficient of Chinese Cabbage Using Eddy Covariance in Northwest China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evapotranspiration and Quantitative Partitioning of Spring Maize with Drip Irrigation under Mulch in an Arid Region of Northwestern China

1
Center for Agricultural Water Research in China, China Agricultural University, Beijing 100083, China
2
National Field Scientific Observation and Research Station on Efficient Water Use of Oasis Agriculture, Wuwei 733009, China
*
Author to whom correspondence should be addressed.
Water 2021, 13(22), 3169; https://doi.org/10.3390/w13223169
Submission received: 12 October 2021 / Revised: 29 October 2021 / Accepted: 29 October 2021 / Published: 10 November 2021

Abstract

:
To examine evapotranspiration (ETc), soil evaporation (Es), and transpiration (Tr), and partitioning of ETc, a two-year field experiment was carried out in a maize field with drip irrigation under mulch in an arid region of northwestern China in 2017 and 2018. In the experiment we designed two treatments with full irrigation (T1) and growth stage-based strategic regulated deficit irrigation (T2). The applied irrigation of T2 was 40% of the T1 during both late vegetative and reproductive growth stages. Based on the measurements of soil water content (SWC) and Tr, a dual crop coefficient model (SIMDualKc) was calibrated and validated, and daily ETc, Es, and Tr were estimated. The model can simulate well the dynamic variations of SWC and Tr. The calibrated basic crop coefficient at the initial, mid-season, and end growth stages was 0.2, 1.15, and 0.75, respectively. The ETc was 507.9 and 519.1 mm for the T1 treatment, and 428.9 and 430.9 mm for the T2 treatment. The ratios of Tr to ETc were higher for the two treatments, ~90%, for two years. Collectively, both drip irrigation under mulch and strategic deficit irrigation after canopy covering of the ground can significantly reduce the ineffective proportion of ETc and Es.

1. Introduction

Crop evapotranspiration (ETc) is one of the key indicators of field water management, crop irrigation scheduling, and planning and design of farmland water conservancy projects [1]. ETc is divided into two parts, soil evaporation (Es) and plant transpiration (Tr). Among them, Es, known as ineffective water consumption for crop growth and yield, can be decreased by ground coverage or proper irrigation management [2,3]. Tr, associated with photosynthetic carbon fixation through leaf pores, directly decides crop growth and the final yield [4]. However, as two water consumption processes in the farmland, Tr and Es occur simultaneously, so it is difficult to carry out quantitative partitioning. Therefore, accurate determination of crop evapotranspiration and its components is of great significance for guiding field irrigation and improving the water use efficiency.
The FAO-56 dual crop coefficient approach is widely used because it can be used to accurately estimate crop evapotranspiration and realize quantitative partitioning of daily Es and Tr [5]. Fan and Cai [6] and Lu et al. [7] demonstrated that ETc can be accurately estimated by the dual crop coefficient approach. A micro-lysimeter can be used to measure Es, but owing to the limited measuring accuracy of the instrument, the accuracy can merely be controlled within 15–20% [8]. Rosa et al. [9] developed a dual crop coefficient model (SIMDualKc) based on the dual crop coefficient approach, making it easier to partition ETc. Many studies showed that the model has a highly accurate estimation of ETc and its components for wheat, maize, forage, tomato, chili, pea, cucumber, etc., in Brazil, Uruguay, Portugal, Spain, and North China [10,11,12,13,14,15,16,17].
Agricultural irrigation is a large water user in arid regions of Northwest China, which is short of water resources, so the use of a new effective water-saving irrigation technology is of great strategic significance for ensuring the water resources security and ecological safety of Northwest China [18]. Drip irrigation under mulch, which is a new type of water-saving technology integrating the advantages of the mulch film, such as soil temperature conservation, soil moisture conservation, yield increase, and the water-saving advantage of drip irrigation, can be used to decrease Es and increase water use efficiency utilization during the initial stage of crop growth [19]. Thus, it has been widely used in arid regions of Northwest China. Previous studies have indicated that Es was reduced by ~50% with plastic film mulch over the whole growing season [20,21,22]. Fan et al. [23] indicated that plastic mulch decreases the available energy and ETc of maize in an arid region of northwest China, and thus the crop coefficient (Kc). Ding et al. [20] introduced a ground-mulching factor to modify the original soil evaporation coefficient in order to account for the reduction of the evaporation area by plastic film mulch. Zhang et al. [24] found that maize ETc with drip irrigation under mulch was reduced by 2.8–5.2%, with reduced soil evaporation by 45.2% and increased transpiration by 8.9% in Northeastern China. However, there remain very few studies on ETc and its components related to the use of drip irrigation under mulch in arid regions of Northwest China.
Crop regulated deficit irrigation (RDI) is a water-saving and high-yield irrigation technology based on the relationship between crops and water. Moderate water deficit in the growth stage of crops can reduce crop water consumption but has a small impact on the final grain yield, thereby improving water use efficiency [25]. RDI reduces crop water consumption mainly by reducing crop growth and leaf area or canopy coverage, but a reduction in canopy coverage will increase the area of bare soil and increase soil evaporation. For example, water deficit in the seedling or early growth period would delay crop growth and canopy cover time, increasing the proportion of ineffective soil evaporation [26]. After the canopy covers the ground (or the leaf area index is greater than 3.0 m2 m−2), the implementation of the strategic stage of deficit adjustment can ensure the reduction of crop water consumption without increasing the proportion of soil evaporation [27]. Therefore, the timing of RDI is very important to reduce crop water consumption without increasing ineffective soil evaporation.
In this study, a two-year field experiment of maize with drip irrigation under mulch was carried out, and two water treatments were set up, namely full irrigation (T1) and strategic stage regulated deficit irrigation in the late growth and reproductive periods (T2). The SIMDualKc model was used to estimate the ETc and Es and Tr of maize during the whole growth period. The objectives were: (1) to quantify the proportion of ETc and Es and Tr of maize with drip irrigation under mulch, and (2) to compare the differences in water use between the two treatments. These results provide a novel approach for efficient water management by strategic growth stage-based RDI in field maize.

2. Materials and Methods

2.1. Experimental Area

The experiment was conducted at the Shiyanghe Experiment Station, China Agricultural University in 2017 and 2018. The station is located in Liangzhou District, Wuwei, Gansu Province, northwest China (37°51′ N, 102°52′ E, at an altitude of 1581 m). The area has a typical continental temperate climate (arid inland desert climate) and abounds in photothermal resources. The annual sunshine duration exceeds 3000 h; the frost-free season lasts for more than 150 d; the annual average temperature is 8 °C and accumulated temperature above 0 °C is higher than 3550 °C; the multi-year average wind speed is 1.3 m s−1; the multi-year average precipitation is 164 mm; the groundwater depth is greater than 30 m. The soil in the experimental area is light sandy loam. The average dry bulk density in the 100 cm soil layer of the root zone is 1.38 g cm−3, with an average field capacity (θFC) of 0.32 cm3 cm−3 and permanent wilting point (θWP) of 0.13 cm3 cm−3.

2.2. Experiment Design

A randomized block experiment was used, with two irrigation treatments, i.e., full irrigation (T1) and regulated deficit irrigation (T2). Each treatment had three replicates, and there were six plots in total. Each plot had a size of 7 × 4.5 m, and the plots stayed unchanged in terms of size and treatment location in the two years. The planting crop was spring maize (Xianyu 335), which was sown on 29 April 2017, and the harvest date was 24 September, and the growth period was 148 days; in 2018, the planting was carried out on 26 April, the harvest date was 24 September, and the length of the growth period was 151 days. We used drip irrigation under the film, and each plot was laid with three white transparent films (Figure 1). The film width was 1.4 m, and each film had three drip irrigation tapes. The seeds were sown on one side of the drip irrigation tapes under the mulch, with a pore diameter of 5 cm, row spacing of 50 cm, and plant spacing of 25 cm. The film coverage rate was one minus the sum of the bare soil area per unit area and the film hole area, which was about 80%. The dripper flow rate was 2.5 L h−1, the dripper spacing was 30 cm, and the working pressure was 0.1 MPa. Nitrogen fertilizer of 250 kg ha−1, phosphate fertilizer of 60 kg ha−1, and potassium fertilizer of 139 kg ha−1 were applied during the whole growth period. Nitrogen fertilizer of 60 kg ha−1 was applied before sowing, and the remaining nitrogen fertilizer was applied four times. Other agronomic measures were consistent with local field management.

2.3. Irrigation Management

For the T1 treatment, irrigation scheduling was designed based on both the water requirements of the crop estimated by the FAO-56 approach and on the measured soil water content. The irrigation amount was set to 100% ETc or θFC. The water amount for the T2 treatment was 40% of that for T1 during both the late vegetative and reproductive growth stages and irrigated to θFC both at the seedling stage and the filling stage. For T1 treatment, unified irrigation was performed before the soil water content decreased below the level of readily available water (RAW). ETc was determined according to the reference evapotranspiration (ET0) and crop coefficient (Kc), while Kc was determined by canopy cover (fc) calculation. Table 1 shows the irrigation time and amount for T1 and T2 in 2017 and 2018. The total irrigation amount for T1 and T2 was 433 and 337 mm in 2017, and 382 and 347 mm in 2018, respectively.

2.4. Data Measurements

The meteorological data were measured by a 2 m-high automatic weather station (Hobo, Onset Computer Corporation, Cape Cod, MA, USA) at the Experimental Station. The data included solar radiation (Rs), air temperature (Ta), relative humidity (RH), 2 m wind speed (u2), and precipitation (P) recorded every 15 min. ET0 was calculated using the FAO-56 Penman–Monteith equation [5]. The average wind speed during the growth period was 0.7 m s−1 in 2017 and 0.66 m s−1 in 2018. The average Rs during the growth period was 223.54 W m−2 in 2017 and 213.4 W m−2 in 2018. Figure 2 shows the ET0, P, and maximum and minimum Ta (Tmax, and Tmin) in 2017 and 2018.
The volumetric soil water content (SWC, cm3 cm−3) was measured in 10 cm increments in depths of 0–200 cm using a neutron probe (CPN-503 Hydroprobe, InstroTek, San Francisco, CA, USA). One neutron tube was installed at the center of each plot. SWC was measured every 7–10 d, and an additional measurement was made before and after irrigation and after rain. The soil drying method was used for calibration.
Maize transpiration was measured by the wrapped sap flowmeter Flow32-1k (Dynamax Inc., Houston, TX, USA). Three uniformly growing maize plants in each plot were selected for wrapping. Before wrapping, the stem diameter of maize at the wrapping site was measured with a vernier caliper with an accuracy of 0.01 mm. An average value was used to calculate the cross-sectional area of the maize stalk and then the cross-sectional area was converted into the transpiration of the plot based on the leaf area index as follows:
T r = 1 N i = 1 n Q di LA i LAI
where Tr is the transpiration rate of the plot (mm d−1); Qdi is the sap flow per plant of the i-th plant (L d−1); LAi is the leaf area of the i-th plant (m2); LAI is the leaf area index (m2 m−2).
The crop height (hc) was measured with a ruler every 10–15 d. The canopy coverage (fc) was measured by photographing above the crop perpendicular to the ground. The ratio of the green area to the total area in the photo was equal to fc. The root zone depth (Zr) was measured at each growth stage by root drilling.

2.5. Quantitative Partitioning of ETc Using the SIMDualKc Model

The SIMDualKc model calculates daily crop ETc by considering both Es and Tr based on the soil water balance and dual Kc method [9,28]. In the model, actual crop ETc is computed as follows:
ETc = (Ks·Kcb + Ke)ET0
where Kcb is the basal crop coefficient, Ke is the soil evaporation coefficient, Ks is the water stress coefficient [0, 1], and ET0 is the reference evapotranspiration. The SIMDualKc model was used to calculate ETc and its components by simulating the dynamic variations of SWC in the root zone. The input data of the model included soil data (field water holding capacity, withering coefficient, saturated moisture content), meteorological data, crop growth data (start and end dates of each growth stage, root depth, plant height, canopy coverage), and irrigation data (irrigation amount and date). The model also considers the effects of mulching film coverage, groundwater recharge, surface runoff, and deep percolation on Tr. Before running the model, the total evaporable water (TEW), readily evaporable water (REW), depth of evaporation layer (Ze), basic crop coefficient (Kcb), and soil water depletion fraction (p) were calibrated.
To calibrate the model parameters, according to the FAO-56 method [5,28], the whole growth period of maize was divided into the initial stage (from seed sowing to fc = 10%), development stage (10% < fc < 80%), mid-season stage (from fc = 80% to maturing) and late-season stage (from maturing to harvest). The average growth indicators of maize in 2017 and 2018 are shown in Table 2 for each treatment. The parameters were calibrated by the trial-and-error method. The simulated soil water content was compared with the measured value. When the error between the simulated and the measured values reached a minimum, the parameter calibration process ended [28,29]. In this study, the measured SWC of 2017 was used for parameter calibration while the data of 2018 were used for verification. The initial values of TEW, REW, Ze, Kcb, and p were set to be equal to the values recommended by Allen et al. [5] and corrected according to the local meteorological conditions and crop factors. Because drip irrigation under mulch was used, the irrigation water–soil wetting ratio (fw) was 0.4 and the film mulching rate was 0.6. The irrigation amount did not exceed the water capacity of the root layer, so deep-water seepage or deep percolation was not taken into consideration. Surface runoffs were not detected in the two years. A simulation was performed using the given Kcb and p. Since the T2 treatment caused some limitations on the growth of maize, the fc of the T2 treatment decreased somewhat at the mid-season and late-season stages compared with T1. Therefore, Kcb was adjusted according to the mid-and late-season stages’ measured values of fc.
Model performance was assessed using the regression coefficient (b), determination coefficient (R2), root mean square error (RMSE), maximum absolute error (Emax), average absolute error (AAE), Willmott index of agreement (dIA), and Nash and Sutcliffe modeling efficiency (EF) between the simulated value and the measured value [13,30,31,32,33]. Among them, b, R2, dIA, and EF were closer to 1.0, and RMSE, Emax, and AAE were closer to 0, indicating that the fitting effect was better.

3. Results and Discussion

Table 3 shows the initial and calibration values of the main model parameters. After calibration, the Kcb of maize with drip irrigation under mulch at the initial stage (Kcb-ini), mid-season stage (Kcb-mid), and late season stage (Kcb-end) were equal to 0.2, 1.15, and 0.55, respectively. The values of Kcb obtained in this study were similar to those in the existing studies and sit within the reviewed and updated range of Kcb for field maize based on accurate crop ETc measurement and FAO56 method by Pereira et al. [34]. Chauhdary et al. [35] presented Kcb-mid = 0.93, Kcb-end = 0.47 for dripped maize with high grain moisture; they used the SALTMED model and gravimetric SWC measurements in Faisalabad, Pakistan. The experimental results achieved by Gimenez et al. [11] in western Uruguay showed that Kcb-ini = 0.15, Kcb-mid = 1.05, and Kcb-end = 0.3. Martins et al. [36] studied maize with sprinkling irrigation and drip irrigation under organic film in southern Brazil and showed that Kcb-ini = 0.2, Kcb-mid = 1.12, and Kcb-end = 0.2. Rodrigues et al. [37] conducted a study on maize under full irrigation and deficit drip irrigation in Portugal and found that Kcb-ini = 0.15, Kcb-mid = 1.15, and Kcb-end = 0.4. Paredes et al. [38], in Portugal, showed by using the AquaCrop model that KcTr,x = 1.18. Paredes et al. [12] in 2014 showed that Kcb-ini = 0.15, Kcb-mid = 1.15, and Kcb-end = 0.3. Yan et al. [39] studied summer maize under different drip irrigation conditions using the SIMDualKc model in Yangling, Shaanxi, concluding that Kcb-ini = 0.15, Kcb-mid = 1.13, and Kcb-end = 0.2. Zhao et al. [40] studied summer maize in Beijing, concluding that Kcb-ini = 0.2, Kcb-mid = 1.1, and Kcb-end = 0.45. Li et al. [25] studied maize by drip irrigation under mulch in northeastern Inner Mongolia, concluding that Kcb-ini = 0.15, Kcb-mid = 1.05, and Kcb-end = 0.4. The slightly higher Kcb-end might be due to the incomplete senescence of maize.
The measured and simulated SWC in the root zone of the two treatments in 2017 and 2018 are shown in Figure 3. The goodness-of-fit statistic of calibration and verification are shown in Table 4. The simulated value and measured SWC fit well. The simulated SWC can capture a dynamic process in which the SWC increased in a short period with irrigation or rainfall, and then gradually decreased due to ETc. The regression coefficient b was 0.96–1.07, R2 was 0.84–0.95, RMSE was 0.005–0.008 cm3 cm−3, AAE 0.01 was cm3 cm−3, Emax 0.025 was cm3 cm−3, and dIA reached up to 0.96, which was better than the results of the study of rain-fed maize in Inner Mongolia by Wu et al. [41]. These results were slightly lower than those found by Zhao et al. [39] on summer maize in Beijing (b = 0.91–1.01, R2 = 0.87–0.93), but the relative error of SWC in this study was lower than 10%, suggesting that the SIMDualKc model was accurately able to calculate SWC and can be used to calculate ETc of maize and its partitioning [9].
The Es, Tr, and ETc of maize were estimated using the calibrated and verified SIMDualKc model. Daily Ke, Kcb, and Kcbadj, as well as Es, Tr, and ETc, and measured Tr for T1 and T2 in 2017 and 2018 are shown in Figure 4 and Figure 5, respectively. The goodness-of-fit statistics of the measured and simulated Tr are presented in Table 5. The simulated and measured Tr had the same changing trend during the growth period. The b was 0.91–1.04, R2 was 0.91–0.97, RMSE was 0.366–0.389 mm d−1, AAE < 0.5 mm d−1, Emax was 1.163 mm d−1, dIA > 0.95, and EF 0.80–0.91. Although Tr was only verified during the mid-to-late growth period, we concluded that the model can estimate Tr throughout the growth period since it accurately simulated SWC throughout the growth period. Qiu et al. [42] compared tomato ETc measured by a lysimeter with SIMDualKc simulations and found that b was 0.91–1.13 and R2 was 0.55–0.82. Yan et al. [17] compared measured Tr values of greenhouse cucumber with simulations and demonstrated that the R2 was 0.89–0.92 and RMSE was 0.36–0.51 mm d−1. Our results were similar to those of previous studies. Overall, after being calibrated, the SIMDualKc model can better simulate the changes in ETc of maize with drip irrigation under mulch during the growth period.
Es and Tr values and their ratios to ETc in different growth stages of maize are shown in Table 6. In 2017, the ETc for T1 and T2 was 507.9 and 428.9 mm, Es was 32.0 and 43.6 mm, and Tr was 476.0 and 385.3 mm, respectively during the whole growth period of maize. In 2018, the ETc for T1 and T2 was 519.1 and 430.9 mm, Es was 35.2 and 43.4 mm, and Tr was 484.0 and 387.5 mm, respectively during the whole growth period of maize. There were large differences in ETc, Es, and Tr between T1 and T2. In particular, there was a difference of 90.7–96.5 mm in Tr, which occurred in the middle growth period. The pattern was similar for two years, which suggests that drip irrigation with film mulching can significantly reduce soil evaporation regardless of whether full or regulated deficit irrigation are used.
Tr was the major component of ETc, with the Tr/ETc ratio of 93.7% and 89.8% for T1 and T2 in 2017, and 93.2% and 89.9% in 2018, respectively. Although Tr and ETc decreased for T2, the Tr/ETc ratio did not decrease significantly, suggesting that the growth-based RDI strategy maintains a higher percentage of crop effective transpiration. The Es/ETc ratio obtained for T1 in the two years was 6.3% and 6.8%, while it was 10.1% and 10.2% for T2, respectively. T2 caused higher evaporation than T1 for the reason that T2 restricted the growth of maize and the fc for T2 was lower than that for T1 at the mid-season stage and late-season stage, causing an increase in the exposed soil area, thus increasing the Es. In the early stage of growth, the fc of maize was very low, with Es as the major active component, and the Es/ETc ratio was highest, in the range of 39–49.4%. At the development stage, the evaporation ratio was 9.9–12.2% in 2017 and 1–1.6% in 2018. Such a large difference was due to a decrease of 18.4 mm in rainfall and a decrease in irrigation volume of 7.2 mm in the same period in 2018. The Es in the mid-growth period in 2017 was smaller than that in the late-growth period, and the opposite was true in 2018. This was because the rainfall in the mid-growth period in 2018 was 140.2 mm, which was 57.6 mm more than in 2017, which led to an increase in soil evaporation. These results indicated that soil evaporation is greatly affected by the degree and coefficient of soil surface moisture and canopy coverage. For efficient crop water management practices, inefficient water consumption can be minimized by covering the ground with the canopy as soon as possible before performing deficit irrigation.
Previous studies have shown that drip irrigation under mulch can effectively reduce soil evaporation, thus improving the effective water use efficiency of crops or increasing Tr/ETc, thereby promoting the growth of biomass and yield [23]. Ding et al. [20] found that for maize for seed under film conditions (fm = 0.7) in arid regions of Northwest China, Es decreased by 55.7% compared to film-free conditions, while Tr was higher. Martins et al. [36] found that the Es/ETc ratio in a maize field was 8–9% under drip irrigation with straw mulch. Li et al. [19] found that that the maize Es/ETc ratio was 19.85–20.29% with film-mulched treatment but 26.15–27.23% without mulch in northeastern Inner Mongolia. Kang et al. [43] studied irrigated maize without mulch in the Guanzhong area, concluding that the Es/ETc ratio was 26%. In this study, the Es/ETc ratio of the two treatments under the condition of mulching drip irrigation were 6.3–10.2%, which is lower than the results of previous studies, indicating that drip irrigation under mulching mainly increases the effective transpiration rate of crops by reducing soil evaporation to save water and increase yield.
The Es/ETc ratios were 10.1% and 10.2% for T2 for the two years, respectively, which is slightly higher than those of T1, at 6.3–6.8%. We started to implement water deficits in the late growth period after the canopy covered the ground, which might cause leaf curling, reduce the canopy coverage, and increase the area of bare soil and evaporated surface. Comas et al. also found that in addition to reducing crop growth and leaf area, water deficit also increased the proportion of rolled leaves, thereby reducing canopy coverage [27]. In this study, due to the use of drip irrigation under the mulch, the area of irrigated wetness and bare soil was small. Even though RDI reduced the canopy coverage and increased the bare soil area, the actual wet soil evaporation area did not increase, so there was no significant increase in Es. These results indicate that in the practices of efficient water management for crops, sufficient irrigation in the early stage of growth can be used to quickly cover the ground in the canopy and then implement the strategic stage of RDI. At the same time, combined with high-efficiency water-saving irrigation methods such as drip irrigation under mulch, it can reduce water use but does not increase the proportion of effectless water.
Although our study area is arid and cold with an annual average temperature of 8 °C, our methods and result patterns can be extended to other areas. The purpose of our study was to estimate ETc and its components to support irrigation scheduling using the SIMDualKc model based on daily soil water balance. The estimation accuracy can be improved if ones take into account soil water infiltration together with the root water uptake [44,45,46]. Further work will be needed to incorporate the two processes into dynamic soil water equations, e.g., using the Richards equation, for accurate partitioning of ETc and soil water flow.

4. Conclusions

A two-year experiment of full irrigation and regulated deficit irrigation of maize with drip irrigation under mulch was conducted in an arid region of Northwest China. The daily evapotranspiration (ETc), soil evaporation (Es), and transpiration (Tr) of maize during its whole growth period and their ratios were calculated using the calibrated dual crop coefficient model SIMDualKc. Then, the differences in ETc and its components between the two treatments were analyzed, drawing the following conclusions: (1) The SIMDualKc model can well simulate the dynamic variations of soil water content and plant transpiration in the maize field with drip irrigation under mulch, and can be used to calculate the evapotranspiration, soil evaporation, and transpiration of maize during its whole growth period; (2) a local basic crop coefficient was obtained for maize with drip irrigation under mulch in an arid region of Northwest China; (3) drip irrigation under mulch can significantly reduce the proportion of soil evaporation, and increase the proportion of plant transpiration that is effective for crop production. Growth-based strategic RDI can reduce crop water use without significantly increasing the proportion of ineffective soil evaporation.

Author Contributions

Experiment design: C.X. and R.D. Data analysis: C.X., J.S. and Y.L. Contributed reagents/materials/analysis tools: C.X., R.D., J.S. and Y.L. Manuscript writing: C.X., R.D., J.S. and Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (51790534 and 52179051) and the China Agriculture Research System of MOF and MARA (CARS-03).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study is available on request from the corresponding author.

Acknowledgments

The authors would like to express their gratitude to the funding agencies, the editor, and reviewers for leveraging the quality of this work, and the students who participated in the fieldwork and laboratory work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhao, N.; Liu, Y.; Cai, J.; LS, P. Validation and application of dual crop coefficient model SIMDual_Kc. Trans. Chin. Soc. Agric. Eng. 2011, 27, 89–95. [Google Scholar]
  2. Allen, R.G. Using the FAO-56 dual crop coefficient method over an irrigated region as part of an evapotranspiration intercomparison study. J. Hydrol. 2000, 229, 27–41. [Google Scholar] [CrossRef]
  3. Zhao, L.; Ji, X. Quantification of transpiration and evaporation over agricultural field using the FAO-56 dual crop coefficient approach-a case study of the maize field in an oasis in the middlestream of the Heihe River basin in Northwest China. Sci. Agric. Sin. 2010, 43, 4016–4026. [Google Scholar]
  4. Pieruschka, R.; Huber, G.; Berry, J.A. Control of transpiration by radiation. Proc. Natl. Acad. Sci. USA 2010, 107, 13372–13377. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration-Guidelines for Computing Crop Water Requirements-FAO Irrigation and Drainage Paper 56; FAO: Rome, Italy, 1998; Volume 300, p. D05109. [Google Scholar]
  6. Fan, Y.; Cai, H. Comparison of crop water requirements computed by single crop coefficient approach and dual crop coefficient approach. J. Hydraul. Eng. 2002, 3, 50–54. [Google Scholar]
  7. Lu, X.; Duan, S.; Ma, X.; Bai, S. A comparative study between single crop coefficient method and double crop coefficient method in calculation of water requirement of maize. Water Sav. Irrig. 2012, 11, 18–21. [Google Scholar]
  8. Shi, J.; Gong, D.; Mei, X.; Hao, W.; Ma, X.; Hu, X. Research progress in surface evapotranspiration measuring and partitioning methods. China Rural Water Hydropower 2012, 2, 49–53,59. [Google Scholar]
  9. Rosa, R.D.; Paredes, P.; Rodrigues, G.C.; Alves, I.; Fernando, R.M.; Pereira, L.S.; Allen, R.G. Implementing the dual crop coefficient approach in interactive software. 1. Background and computational strategy. Agric. Water Manag. 2012, 103, 8–24. [Google Scholar] [CrossRef]
  10. Fandiño, M.; Olmedo, J.L.; Martínez, E.M.; Valladares, J.; Paredes, P.; Rey, B.J.; Mota, M.; Cancela, J.J.; Pereira, L.S. Assessing and modelling water use and the partition of evapotranspiration of irrigated hop (Humulus Lupulus), and relations of transpiration with hops yield and alpha-acids. Ind. Crop Prod. 2015, 77, 204–217. [Google Scholar] [CrossRef]
  11. Giménez, L.; Petillo, M.G.; Paredes, P.; Pereira, L.S. Predicting maize transpiration, water use and productivity for developing improved supplemental irrigation schedules in western Uruguay to cope with climate variability. Water 2016, 8, 309. [Google Scholar] [CrossRef] [Green Version]
  12. Paredes, P.; Rodrigues, G.C.; Alves, I.; Pereira, L.S. Partitioning evapotranspiration, yield prediction and economic returns of maize under various irrigation management strategies. Agric. Water Manag. 2014, 135, 27–39. [Google Scholar] [CrossRef]
  13. Qiu, R.; Du, T.; Kang, S.; Chen, R.; Wu, L. Assessing the SIMDualKc model for estimating evapotranspiration of hot pepper grown in a solar greenhouse in Northwest China. Agric. Syst. 2015, 138, 1–9. [Google Scholar] [CrossRef]
  14. Wei, Z.; Paredes, P.; Liu, Y.; Chi, W.W.; Pereira, L.S. Modelling transpiration, soil evaporation and yield prediction of soybean in North China Plain. Agric. Water Manag. 2015, 147, 43–53. [Google Scholar] [CrossRef]
  15. Wu, Y.; Liu, T.; Paredes, P.; Duan, L.; Wang, H.; Wang, T.; Pereira, L.S. Ecohydrology of groundwater-dependent grasslands of the semi-arid Horqin sandy land of inner Mongolia focusing on evapotranspiration partition. Ecohydrology 2016, 9, 1052–1067. [Google Scholar] [CrossRef]
  16. Wang, Z.; Cai, H.; Yu, L.; Wang, X.; Shi, X. Estimation of evapotranspiration and soil evaporation of winter wheat in arid region of Northwest China based on SIMDualKC model. Trans. Chin. Soc. Agric. Eng. 2016, 32, 126–136. [Google Scholar]
  17. Yan, H.; Wu, H.; Zhang, C.; Acquah, S.; Zhao, B.; Huang, S. Estimation of greenhouse cucumber evapotranspiration in different seasons based on modified dual crop coefficient model. Trans. Chin. Soc. Agric. Eng. 2018, 34, 117–125. [Google Scholar]
  18. Zhang, Y.; Wang, J.; Gong, S.; Xu, D.; Sui, J.; Wu, Z. Analysis of water saving and yield increasing mechanism in maize field with drip irrigation under film mulching based on transpiration estimated by sap flow meter. Trans. Chin. Soc. Agric. Eng. 2018, 34, 89–97. [Google Scholar]
  19. Li, R.; Zhao, J.; Shi, H.; Wang, N.; Qi, Y.; Feng, Y.; Jia, Q. Simulation of evapotranspiration of corn under mulched drip irrigation based on SIMDual_Kc model in Tongliao area of Inner Mongolia. Trans. Chin. Soc. Agric. Eng. 2018, 34, 127–134. [Google Scholar]
  20. Ding, R.; Kang, S.; Zhang, Y.; Hao, X.; Tong, L.; Du, T. Partitioning evapotranspiration into soil evaporation and transpiration using a modified dual crop coefficient model in irrigated maize field with ground-mulching. Agric. Water Manag. 2013, 127, 85–96. [Google Scholar] [CrossRef]
  21. Hou, X.; Wang, F.; Han, J.; Kang, S.; Feng, S. Duration of plastic mulch for potato growth under drip irrigation in an arid region of Northwest China. Agric. Forest Meteor. 2010, 150, 115–121. [Google Scholar] [CrossRef]
  22. Zhou, L.; Li, F.; Jin, S.; Song, Y. How two ridges and the furrow mulched with plastic film affect soil water, soil temperature and yield of maize on the semiarid Loess Plateau of China. Field Crops Res. 2009, 113, 41–47. [Google Scholar] [CrossRef]
  23. Fan, Y.; Ding, R.; Kang, S.; Hao, X.; Du, T.; Tong, L.; Li, S. Plastic mulch decreases available energy and evapotranspiration and improves yield and water use efficiency in an irrigated maize cropland. Agric. Water Manag. 2017, 179, 122–131. [Google Scholar] [CrossRef]
  24. Zhang, Y.; Wang, J.; Gong, S.; Xu, D.; Sui, J.; Wu, Z.; Mo, Y. Effects of film mulching on evapotranspiration, yield and water use efficiency of a maize field with drip irrigation in Northeastern China. Agric. Water Manag. 2018, 205, 90–99. [Google Scholar] [CrossRef]
  25. Cai, H.; Kang, S.; Zhang, Z.; Chai, H.; Hu, X.; Wang, J. Proper growth stages and deficit degree of crop regulated deficit irrigation. Trans. Chin. Soc. Agric. Eng. 2000, 16, 24–27. [Google Scholar]
  26. Trout, T.J.; DeJonge, K.C. Water productivity of maize in the US high plains. Irrig. Sci. 2017, 35, 251–266. [Google Scholar] [CrossRef] [Green Version]
  27. Comas, L.H.; Trout, T.J.; DeJonge, K.C.; Zhang, H.; Gleason, S.M. Water productivity under strategic growth stage-based deficit irrigation in maize. Agric. Water Manag. 2019, 212, 433–440. [Google Scholar] [CrossRef]
  28. Rosa, R.D.; Paredes, P.; Rodrigues, G.C.; Fernando, R.M.; Alves, I.; Pereira, L.S.; Allen, R.G. Implementing the dual crop coefficient approach in interactive software: 2. Model testing. Agric. Water Manag. 2012, 103, 62–77. [Google Scholar] [CrossRef]
  29. Zhao, N.; Liu, Y.; Cai, J.; Paredes, P.; Rosa, R.D.; Pereira, L.S. Dual crop coefficient modelling applied to the winter wheat–summer maize crop sequence in North China Plain: Basal crop coefficients and soil evaporation component. Agric. Water Manag. 2013, 117, 93–105. [Google Scholar] [CrossRef]
  30. Legates, D.R.; McCabe, G.J., Jr. Evaluating the use of “goodness-of-fit” Measures in hydrologic and hydroclimatic model validation. Water Resour. Res. 1999, 35, 233–241. [Google Scholar] [CrossRef]
  31. Moriasi, D.N. Model evaluation guidelines for systematic quantification of accuracy in watershed simulations. Trans. ASABE 2007, 50, 885–900. [Google Scholar]
  32. Xu, D.; Mermoud, A. Modeling the soil water balance based on time-dependent hydraulic conductivity under different tillage practices. Agric. Water Manag. 2003, 63, 139–151. [Google Scholar] [CrossRef]
  33. Zhang, B.; Liu, Y.; Xu, D.; Zhao, N.; Lei, B.; Rosa, R.D.; Paredes, P.; Paço, T.A.; Pereira, L.S. The dual crop coefficient approach to estimate and partitioning evapotranspiration of the winter wheat–summer maize crop sequence in North China Plain. Irrig. Sci. 2013, 31, 1303–1316. [Google Scholar] [CrossRef]
  34. Pereira, L.S.; Paredes, P.; Hunsaker, D.J.; López-Urrea, R.; Mohammadi Shad, Z. Standard single and basal crop coefficients for field crops. Updates and advances to the FAO56 crop water requirements method. Agric. Water Manag. 2021, 243, 106466. [Google Scholar] [CrossRef]
  35. Chauhdary, J.N.; Bakhsh, A.; Ragab, R.; Khaliq, A.; Engel, B.A.; Rizwan, M.; Shahid, M.A.; Nawaz, Q. Modeling corn growth and root zone salinity dynamics to improve irrigation and fertigation management under semi-arid conditions. Agric. Water Manag. 2020, 230, 105952. [Google Scholar] [CrossRef]
  36. Martins, J.D.; Rodrigues, G.C.; Paredes, P.; Carlesso, R.; Oliveira, Z.B.; Knies, A.E.; Petry, M.T.; Pereira, L.S. Dual crop coefficients for maize in southern Brazil: Model testing for sprinkler and drip irrigation and mulched soil. Biosyst. Eng. 2013, 115, 291–310. [Google Scholar] [CrossRef]
  37. Rodrigues, G.C.; Paredes, P.; Gonçalves, J.M.; Alves, I.; Pereira, L.S. Comparing sprinkler and drip irrigation systems for full and deficit irrigated maize using multicriteria analysis and simulation modelling: Ranking for water saving vs. farm economic returns. Agric. Water Manag. 2013, 126, 85–96. [Google Scholar] [CrossRef]
  38. Paredes, P.; de Melo-Abreu, J.P.; Alves, I.; Pereira, L.S. Assessing the performance of the FAO AquaCrop model to estimate maize yields and water use under full and deficit irrigation with focus on model parameterization. Agric. Water Manag. 2014, 144, 81–97. [Google Scholar] [CrossRef] [Green Version]
  39. Yan, S.; Zhang, F.; Wu, Y.; Qiang, S.; Zou, H.; Xiang, Y.; Fan, J.; Tian, J. Estimation of drip irrigated summer maize soil water content and evapotranspiration based on SIMDualKc model. Trans. Chin. Soc. Agric. Eng. 2017, 33, 152–160. [Google Scholar]
  40. Zhao, N.; Liu, Y.; Cai, J.; Yu, F.; Li, C. Research on soil evaporation of summer maize by field measurement and model simulation. Trans. Chin. Soc. Agric. Eng. 2012, 28, 66–73. [Google Scholar]
  41. Wu, Y.; Liu, T.; Paredes, P.; Duan, L.; Pereira, L.S. Water use by a groundwater dependent maize in a semi-arid region of Inner Mongolia: Evapotranspiration partitioning and capillary rise. Agric. Water Manag. 2015, 152, 222–232. [Google Scholar] [CrossRef] [Green Version]
  42. Qiu, R.; Du, T.; Chen, R. Application of the dual crop coefficient model for estimating tomato evapotranspiration in greenhouse. J. Hydaul. Eng. 2015, 46, 678–686. [Google Scholar]
  43. Kang, S.; Gu, B.; Du, T.; Zhang, J. Crop coefficient and ratio of transpiration to evapotranspiration of winter wheat and maize in a semi-humid region. Agric. Water Manag. 2003, 59, 239–254. [Google Scholar] [CrossRef]
  44. Difonzo, F.V.; Masciopinto, C.; Vurro, M.; Berardi, M. Shooting the numerical solution of moisture flow equation with root water uptake models: A python tool. Water Resour. Manag. 2021, 35, 2553–2567. [Google Scholar] [CrossRef]
  45. De Melo, M.L.A.; de Jong van Lier, Q. Revisiting the Feddes reduction function for modeling root water uptake and crop transpiration. J. Hydrol. 2021, 603, 126952. [Google Scholar] [CrossRef]
  46. Wu, X.; Zuo, Q.; Shi, J.; Wang, L.; Xue, X.; Ben-Gal, A. Introducing water stress hysteresis to the Feddes empirical macroscopic root water uptake model. Agric. Water Manag. 2020, 240, 106293. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of maize planting with drip irrigation under mulch (a) and photo (b).
Figure 1. Schematic diagram of maize planting with drip irrigation under mulch (a) and photo (b).
Water 13 03169 g001
Figure 2. Daily variations of reference evapotranspiration (ET0), precipitation (P), daily minimum relative humidity (RHmin), and maximum and minimum air temperature (Tmax and Tmin) with days after planting (DAP) during the whole growth period of maize in 2017 (a,b) and 2018 (c,d).
Figure 2. Daily variations of reference evapotranspiration (ET0), precipitation (P), daily minimum relative humidity (RHmin), and maximum and minimum air temperature (Tmax and Tmin) with days after planting (DAP) during the whole growth period of maize in 2017 (a,b) and 2018 (c,d).
Water 13 03169 g002
Figure 3. Measured and simulated seasonal soil water content (SWC) for two water treatments (T1 and T2) with days after planting (DAP) during the whole growth period of maize in 2017 (a,b) and 2018 (c,d).
Figure 3. Measured and simulated seasonal soil water content (SWC) for two water treatments (T1 and T2) with days after planting (DAP) during the whole growth period of maize in 2017 (a,b) and 2018 (c,d).
Water 13 03169 g003
Figure 4. Dynamic variations of basic crop coefficient (Kcb), actual adjustment Kcb (Kcbadj) and soil evaporation coefficient (Ke) for two water treatments (T1 and T2) with days after planting (DAP) during the whole growth period of maize in (a,b) and 2018 (c,d).
Figure 4. Dynamic variations of basic crop coefficient (Kcb), actual adjustment Kcb (Kcbadj) and soil evaporation coefficient (Ke) for two water treatments (T1 and T2) with days after planting (DAP) during the whole growth period of maize in (a,b) and 2018 (c,d).
Water 13 03169 g004
Figure 5. Seasonal variations of simulated evapotranspiration (ETc), transpiration (Tr), and soil evaporation (Es), and measured Tr for two water treatments (T1 and T2) with days after planting (DAP) during the whole growth period of maize in (a,b) and 2018 (c,d).
Figure 5. Seasonal variations of simulated evapotranspiration (ETc), transpiration (Tr), and soil evaporation (Es), and measured Tr for two water treatments (T1 and T2) with days after planting (DAP) during the whole growth period of maize in (a,b) and 2018 (c,d).
Water 13 03169 g005
Table 1. Irrigation scheduling for maize with drip irrigation under mulch for two water treatments (T1 and T2) during the whole growth period of maize in 2017 and 2018.
Table 1. Irrigation scheduling for maize with drip irrigation under mulch for two water treatments (T1 and T2) during the whole growth period of maize in 2017 and 2018.
YearsDatesIrrigation Depth (mm)
T1T2
20175/33030
5/311716
6/18–6/194041
6/30–7/15522
7/10–7/115321
7/22–7/24120120
8/84357
8/317530
20184/273030
5/92121
6/105051
6/236025
6/304117
7/93716
7/18–7/1968100
8/15–8/167587
Table 2. Growth traits for two water treatments (T1 and T2) during the whole growth period of maize in 2017 and 2018.
Table 2. Growth traits for two water treatments (T1 and T2) during the whole growth period of maize in 2017 and 2018.
TraitsYearsTreatmentsGrowth Stages
InitialDevelopmentMid-SeasonLate-SeasonWhole Season
Growth length (d)2017T127276430148
T228295833148
2018T132256430151
T232285932151
Plant height (m)2017T10.31.52.93.1
T20.291.52.42.4
2018T10.41.42.93.1
T20.41.22.72.7
Root depth (m)2017T10.10.40.740.74
T20.20.50.650.65
2018T10.20.440.70.7
T20.250.50.70.7
Canopy cover2017T10.10.970.930.6
T20.10.90.850.56
2018T10.10.950.90.6
T20.10.880.850.55
Table 3. Initial and calibrated values of key parameters for the SIMDualKc model.
Table 3. Initial and calibrated values of key parameters for the SIMDualKc model.
ParametersInitial ValuesCalibrated
Crop parameters
Kcb-ini0.150.2
Kcb-mid1.151.15
Kcb-end0.500.55
pini0.550.55
pmid0.550.55
pend0.550.55
Soil parameters
REW (mm)1012
TEW (mm)3030
Ze (m)0.120.15
Note: Kcb and p are the maize basal crop coefficient and the soil–water depletion fraction, respectively, for no stress at the initial (ini), mid-season (mid) and late-season (end) stages; REW and TEW are readily and total evaporable water, respectively; and Ze is the depth of the soil evaporation layer. The emboldened values are calibrated parameters that are different from the initial ones.
Table 4. Statistical indicators of goodness-of-fit between measured and simulated seasonal soil water content (SWC) for the two treatments (T1 and T2) in 2017 and 2018.
Table 4. Statistical indicators of goodness-of-fit between measured and simulated seasonal soil water content (SWC) for the two treatments (T1 and T2) in 2017 and 2018.
YearsTreatmentsbR2RMSE (cm3·cm−3)AAE (cm3·cm−3)Emax (cm3·cm−3)dIAEF
2017T10.980.840.0080.0060.0130.960.82
T20.960.900.0080.0060.0220.970.89
2018T11.070.950.0050.0040.0100.990.94
T21.000.950.0080.0060.0250.990.94
Note: b, linear regression coefficient; R2, coefficient of determination; RMSE, root mean square error; AAE, average absolute error; Emax, maximum absolute error; dIA, Willmott index of agreement; and EF, the Nash and Sutcliffe modeling efficiency.
Table 5. Statistical indicators of goodness-of-fit between measured and simulated seasonal plant transpiration (Tr) for the two treatments (T1 and T2) in 2017 and 2018.
Table 5. Statistical indicators of goodness-of-fit between measured and simulated seasonal plant transpiration (Tr) for the two treatments (T1 and T2) in 2017 and 2018.
YearsTreatmentsbR2RMSE (mm·d−1)AAE (mm·d−1)Emax (mm·d−1)dIAEF
2017T10.99 0.95 0.366 0.294 1.0600.970.88
T21.00 0.91 0.379 0.293 1.1630.950.80
2018T10.91 0.97 0.367 0.310 0.7090.980.91
T21.04 0.95 0.389 0.346 0.6490.960.82
Note: b, linear regression coefficient; R2, coefficient of determination; RMSE, root mean square error; AAE, average absolute error; Emax, maximum absolute error; dIA, Willmott index of agreement; and EF, the Nash and Sutcliffe modeling efficiency.
Table 6. Soil evaporation (Es), transpiration (Tr), evapotranspiration (ETc) and ratios of evaporation and transpiration to evapotranspiration for the two treatments (T1 and T2) at different growth stages of maize in 2017 and 2018.
Table 6. Soil evaporation (Es), transpiration (Tr), evapotranspiration (ETc) and ratios of evaporation and transpiration to evapotranspiration for the two treatments (T1 and T2) at different growth stages of maize in 2017 and 2018.
Growth StagesYearsEs (mm)Tr (mm)ETc (mm)Es/ETc (%)Tr/ETc (%)
T1T2T1T2T1T2T1T2T1T2
Initial201717.821.518.322.036.243.649.349.450.750.6
201821.021.031.932.952.954.039.739.060.361.0
Development20177.810.070.571.878.381.89.912.290.187.8
20180.81.277.373.378.174.51.01.699.098.4
Mid-season20173.97.9302.9214.2306.9222.11.33.598.796.5
20185.58.6306.4218.3311.9226.91.83.898.296.2
Late-season20172.44.284.277.286.681.42.85.297.294.8
20187.912.668.362.976.375.510.416.789.683.3
Whole season201732.043.6476.0385.3507.9428.96.310.293.789.8
201835.243.4484.0387.5519.1430.96.810.193.289.9
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xuan, C.; Ding, R.; Shao, J.; Liu, Y. Evapotranspiration and Quantitative Partitioning of Spring Maize with Drip Irrigation under Mulch in an Arid Region of Northwestern China. Water 2021, 13, 3169. https://doi.org/10.3390/w13223169

AMA Style

Xuan C, Ding R, Shao J, Liu Y. Evapotranspiration and Quantitative Partitioning of Spring Maize with Drip Irrigation under Mulch in an Arid Region of Northwestern China. Water. 2021; 13(22):3169. https://doi.org/10.3390/w13223169

Chicago/Turabian Style

Xuan, Chenggong, Risheng Ding, Jie Shao, and Yanshuo Liu. 2021. "Evapotranspiration and Quantitative Partitioning of Spring Maize with Drip Irrigation under Mulch in an Arid Region of Northwestern China" Water 13, no. 22: 3169. https://doi.org/10.3390/w13223169

APA Style

Xuan, C., Ding, R., Shao, J., & Liu, Y. (2021). Evapotranspiration and Quantitative Partitioning of Spring Maize with Drip Irrigation under Mulch in an Arid Region of Northwestern China. Water, 13(22), 3169. https://doi.org/10.3390/w13223169

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