Next Article in Journal
Compensation for Agricultural Economic Losses Caused by Restoration of Healthy Eco-Hydrological Sequences of Rivers
Previous Article in Journal
Occurrence of Micropollutants in Wastewater and Evaluation of Their Removal Efficiency in Treatment Trains: The Influence of the Adopted Sampling Mode
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Water Use Effectiveness Is Enhanced Using Film Mulch Through Increasing Transpiration and Decreasing Evapotranspiration

1
Center for Agricultural Water Research in China, China Agricultural University, Beijing 100083, China
2
Wuwei Experimental Station for Efficient Water Use in Agriculture, Ministry of Agriculture and Rural Affairs, Wuwei 733000, China
*
Author to whom correspondence should be addressed.
Water 2019, 11(6), 1153; https://doi.org/10.3390/w11061153
Submission received: 19 April 2019 / Revised: 5 May 2019 / Accepted: 10 May 2019 / Published: 1 June 2019
(This article belongs to the Section Water Use and Scarcity)

Abstract

:
Water shortage is a main limitation of crop growth and yield in drought northwest China, which is an important area of seed maize growth. Plastic film mulch is widely adopted to reduce soil evaporation (E) and conserve water resources, which changes evapotranspiration (ET) and its components, E and transpiration (Tr) and crop growth. The AquaCrop model, one of widely used crop models powered by water, can well simulate crop ET components and growth. However, there are few studies that examine ET partitioning and growth with and without plastic film mulch. The calibrated AquaCrop model was used to partition ET and simulate growth of seed maize with and without plastic film mulch in a drought region of northwest China in 2014 and 2015. The AquaCrop model can well simulate canopy cover curve (CC), and the dynamic and accumulated courses of ET and ET components. Plastic film mulch could advance the growth stage of seed maize and reduce seasoned ET. The initial stage with plastic film mulch was 37–42 days, while it was 46–48 days for no-mulch. Plastic film mulch increased Tr by 14.16% and 14.48% and significantly decreased E by 57.25% and 34.28% in 2014 and 2015, respectively, resulting in the reduction of seasonal total ET. Plastic film mulch increased averaged mid-season crop coefficient for transpiration (Kc Tr) by 0.88% and decreased soil evaporation coefficient (Ke) by 62.50%. Collectively, the results suggest that, in comparison with no-mulch, plastic film mulch advanced crop growth, and decreased total ET and increased Tr related with crop production, i.e., improve water use effectiveness.

1. Introduction

Water resources scarcity is the main cause of crop failure [1]. The trend is further worsened by climate change, since climate change projections indicate that vapor pressure deficit (VPD), and thus, evaporation demand, will become an expected increase [2]. Crop evapotranspiration (ET) is the largest proportion of ET, precipitation and irrigation in cropland water balance component [3,4] because more than 90% of water used in agriculture is lost through soil evaporation and crop transpiration, referred to as evapotranspiration [5]. Soil evaporation cannot promote the increase of crop yield and should be controlled by management measures (e.g., plastic film mulch) [6,7]. On the contrary, transpiration occurs simultaneously through stomata and photosynthesis in leaves and is closely related to crop growth and production [8], whereas the two independent processes Tr and E occur simultaneously and are difficult to distinguish [9]. Therefore, in order to better understand the water consumption law of seed maize and formulate accurate irrigation scheduling, so as to improve the water use effectiveness (EUW, means effective use of water, usually represents a process, and not water use efficiency (WUE)) and save water resources [10], ET needs to be accurately partitioned and growth simulated through models.
There are many crop models which have been widely used to estimate ET and its components and crop growth, such as DSSAT [11], WOFOST [12], APSIM [13], EPIC [14,15], WAVES [16], and Crop-Syst [17]. However, these models are complex and often require a large number of hard-to-obtain input variables and calibration parameters. The AquaCrop model is the newest international water-driven crop model developed by FAO to estimate ET and crop growth and yield [18,19]. The AquaCrop model can simulate canopy cover and estimate Tr and E depending upon canopy cover curve. The AquaCrop model was widely used in annual field crops, e.g., Stricevic et al., Nyakudya and Stroosnijder, Abedinpour et al., Hsiao et al., and Paredes et al. [20,21,22,23,24] for maize; Mkhabela and Bullock, Wang et al., Andarzian et al., Iqbal et al., and Toumi et al. [25,26,27,28,29] for wheat; Linker et al., and Garcia-Vila et al. [30,31] for cotton; Maniruzzaman et al. [32] for rice; Todorovic et al., and Stricevic et al. [20,33] for sunflowers; Araya et al. [34] for barley, and others.
In the drought region of northwest China, a major yield limiting factor is water shortage, where seed maize plays an important role in major cereal crops because of its high price and income. To save water resources by reducing soil evaporation in the field, plastic film mulch practice has been widely used in this area since this century [35,36,37]. Many studies indicate that plastic film mulch can not only reduce soil evaporation, but also accelerate crop growth due to increased soil temperature and weeds control [37,38]. Several studies used the AquaCrop model to partition ET into E and Tr [24,39,40,41]. However, there are few studies that used the AquaCrop model to partition ET and investigate crop growth with and without plastic film mulch.
The objectives of this study were to quantify the differences of ET partitioning and maize growth, and to analyze the causes of the differences between mulched and unmulched plastic film field using the AquaCrop model.

2. Material and Methods

2.1. Experimental Site and Design

The experimental observation occurred, respectively, in 2014 and 2015 from April to September in northwest China’s Gansu province Wuwei city the Shiyanghe Experimental Station for Water-saving in Agriculture and Ecology of China Agricultural University (N 37°51′, E 102°52′, and 1581 m elevation). This region is typical of arid continental temperate climate, where it is rich in light and heat resources. The average annual sunshine time is more than 3000 h, and the hottest monthly daily temperature is 29 °C in July, with more than 150 days of frost-free period. However, the region is short of water, with an average annual precipitation of 164 mm and an average annual evaporation from a free water surface of 2000 mm. The groundwater table is 30–40 m below the ground surface. For sandy loam soil and corresponding mean soil dry bulk density of 1.40 g⸱cm−3, field capacity of 30%, and permanent wilting point of 11% for the 0–100 cm layers, the average saturated water content of 0.41 cm3⸱cm−3. Saturated hydraulic conductivity is 500 mm⸱d−1.
The experimental design was a completely random design, including two treatments with and without plastic film mulch (M, N). Each treatment was repeated three times with a length of 9 m from north to south and 12.6 m from east to west. The bare soil surface was mulched with a transparent plastic thin film strip every 20 cm. The plastic film thickness was 0.04 mm and the width was 140 cm. Each film strip was mulched with four rows of seed maize. For seed maize, holes with a diameter of 5.0 cm, a plant spacing of 25 cm, and a row spacing of 40 cm were opened on the plastic film to sowing of single row male parents and five rows female parents, and there were no maize plants in the bare ground. Therefore, due to the existence of bare soil and holes, the actual mulching rate was 0.70. In 2014 and 2015, seed maize was sown on April 15 and 17, and harvested on September 20.
The experimental irrigation method was border-irrigation with uniform water infiltration according to the local actual agricultural production. The irrigation schedule in the experiment was estimated by subtracting precipitation from preceding ET calculated by crop coefficient method, and the water deficit was calculated by the difference between root area, field capacity and observed soil water content (SWC). The actual irrigation water in the experiment was measured and recorded by turbine flow meters. The precipitation and irrigation schedule in 2014 and 2015 are shown in Table 1. According to the local practice, the ratio of fertilizers application for 420 kg⸱ha−1 N in 2014 and 390 kg⸱ha−1 N in 2015. At the same time, 240 kg⸱ha−1 P2O5 and 50 kg⸱ha−1 K2O were applied in each year. Full doses of potassium and phosphate were applied prior to seeding, with half of the N used as a base treatment at seeding and the rest applied during the first irrigation.

2.2. Data Measurements

2.2.1. Meteorological Data

Meteorological data, including solar radiation (Rs, W⸱m−2), precipitation (P, mm), air temperature (T, °C) relative humidity (RH, %), and wind speed (u2, m⸱s−1), were continuously measured at a standard automatic weather station (Hobo, Onset Computer Corp, Bourne, MA, USA), 2.0 m above the ground. The data were collected every 5 s, and the data recorder recorded an average of 15 min. The FAO Penman-Monteith method [38] was used to calculate daily reference crop evapotranspiration (ET0, mm). The daily maximum, minimum temperatures, precipitation, and ET0 for each growing season are shown in Figure 1.

2.2.2. Canopy Cover

In 2015, canopy cover (CC) was directly measured by an overhead digital camera, and image analysis software was used to calculate the ratio of weekly green leaf pixels in each plot to the whole picture. In 2014, it was derived from the leaf area index (LAI) by the following empirical equation, given by Hsiao et al. [23]:
C C = 1.005 [ 1 e x p ( 1 e x p ( 0.6 L A I ) ] 1.2
Five maize plants were randomly selected and the plant height of each treatment plot was measured every 10 days. The leaf area index (LAI) was taken to be leaf area per unit of area.

2.2.3. Soil Water Content

Volumetric soil water content (SWC) was measured by a tubular soil water probes (Trime-TDR, Imko, Gmbh, Ettlingen, Germany). The maximum soil depth was 1.0 m, with an interval of 0.2 m per week. Each plot was fitted with two detection tubes, one in the mulched area and the other in the bare soil area. The frequency of the SWC measurements would increase before and after each irrigation and precipitation.

2.2.4. Evapotranspiration

In this study, the water balance equation was used to estimate ET.
E T = P + I + W D R Δ S
where P is the precipitation (mm), I is the irrigation quantity (mm), W is the groundwater flow into the root zone (mm), D is the drainage (mm), R is the surface runoff (mm), ΔS is the change of soil moisture storage. The terrain of the experimental site was flat and the precipitation was not strong, so surface runoff and deep drainage could be ignored. W could also be ignored, as the groundwater table was 30–40 m below the ground surface, which is not enough to affected the growth of seed maize. Therefore, the Equation (2) can be simplified as:
E T = P + I Δ S
Δ S = 1000 × Z × ( θ 2 θ 1 )
where Z is the depth of root zone (m) with the change of plant growth, and θ2 and θ1 are average soil water contents in the root zone at two measurement times (m3⸱m−3).

2.3. Model Description

2.3.1. Model Principles and Algorithms

The AquaCrop model is a crop model powered by water. It first calculates transpiration and then converts it into biomass through water productivity, which combines four sub-models: (1) The atmosphere sub-model, treatment of rainfall, evaporative demand (reference evapotranspiration, ET0), and CO2 concentration; (2) the development, growth and yield of crops; (3) the management sub-model, including irrigation and fertilization; and (4) the balance with soil water. The AquaCrop model use a “Kc-ET0” approach to estimate Tr and E, depending on the daily time step used for canopy cover. Tr is calculated as:
T r = C C * K c T r , x E T 0
where ET0 is reference evapotranspiration (mm) calculated using the FAO Penman-Monteith method (mm⸱d−1), KcTr,x is the maximum standard crop transpiration coefficient (non-dimensional) (well-watered soil and complete canopy, canopy cover = 100%), and CC* is the actual crop canopy cover (%) adjusted for micro-advective effects. The actual transpiration is obtained by adjusting Tr to soil water stress conditions using the water stress coefficient Ks (0–1), i.e., T r a d j = K s Δ T r . The coefficient Ks describes the effects of soil water stress on the following crop growth processes including water logging (Ksaer), stomatal closure (Kssto) and soil salinity stress (Kssalt). There was no water stress in this study, so the effect of Ks was not considered.
CC* is calculated as:
C C * = 1.72 C C C C 2 + 0.30 C C 3
Under optimal conditions, the development and senescence of the green canopy is described as:
C C = C C 0 e t   C G C   for   C C C C x / 2
C C = C C x 0.25 ( C C x ) 2 C C 0 e t   C G C   for   C C > C C x / 2
C C = C C x [ 1 0.05 ( e C D C C C x t 1 ) ]   for   C C   decline
where CC0 is the initial canopy cover at 90% emergence, CCx is the maximum canopy cover, CGC is the canopy growth coefficient (d−1), and t is time (d).
CGC and CDC are the two most important parameters for building up the CC, which are modified by soil water stress, as follows:
C G C a d j = K s e x p C G C
C D C a d j = ( 1 K s 8 s e n ) C D C
where CGCadj is CGC adjusted for water stress (d−1), Ksexp is the water stress coefficient for canopy expansion (0–1), CDCadj is CDC adjusted for water stress (d−1), and Kssen is the water stress coefficient for early canopy senescence (0–1).
The soil evaporation is also obtained by CC* and ET0:
E = K r ( 1 C C * ) K e x E T 0
The calculation equation of soil evaporation under plastic film mulch is as follows:
E = K r   ( 1 f m f r   m u l c h ) ( 1 C C * ) K e x E T 0
where Kex is the maximum soil evaporation coefficient (non-dimensional) of completely wet and unshaded soil surface, and AquaCrop model defaults to 1.10. Kr is the evaporation reduction coefficient (0–1), and when the surface soil moisture is insufficient, Kr < 1, in response to atmospheric evaporation requirements. fm is a regulator of the influence of mulch on soil evaporation, fr mulch is a fraction of the mulch (fr mulch = 0.70).

2.3.2. Parameters Sensitivity Analysis

Accurate estimation of the AquaCrop model parameters is imperative for the accuracy of model outputs. A sensitivity analysis (SA) quantifies the impact of each uncertain factor (parameter or driving variable) on the variability of the model output and helps to understand the model’s response to the changing behavior of these factors [42,43]. The SA methods are typically classified as local (i.e., derivative-based) or global [44]. A local SA only studies the influence of the change of one parameter on the output of the model and fixes all other parameters to an arbitrary value [42]. For these reasons, local methods are criticized as not suitable for nonlinear models [45]. A global SA (GSA) studies the average response of model output to changes in all parameters within a defined range. GSA technique fully considers parameter interactions or nonlinear responses, but requires repeated analysis of the changing parameter values in the parameter space, so it has high requirements for calculation and workload [46,47,48]. There are many existing GSA techniques, and this study used a combination of screening and variance-based methods. First, the Morris screening method [49]( was used to identify the parameters with the greatest influence. Then, using the Sobol’ [50] sensitivity analysis method, the first-and higher-order effects of these most influential model parameters on the selected model output were quantified. In this study, the Morris SA method was used to screen the 29 input parameters affecting the AquaCrop model to select the first 15 input parameters that are most sensitive to E and Tr, and then the Sobol’s SA method was used to analyze the global sensitivity of the selected 15 parameters. The 15 parameters were base temperature (BaseTemp), upper temperature (UpperTemp), soil water depletion factor for canopy expansion—upper threshold (PexpUpper), soil water depletion factor for canopy expansion— lower threshold (PexpLower), shape factor for water stress coefficient for canopy expansion (SFexp), soil water depletion factor for canopy senescence—upper threshold (PsenUpper), shape factor for water stress coefficient for canopy senescence (SFsen), crop coefficient when canopy is complete but prior to senescence (KcTr,x), minimum effective rooting depth (MinRootDepth), maximum effective rooting depth (MaxRootDepth), soil surface covered by an individual seedling at 90% emergence (CC0), canopy growth coefficient (CGC), canopy decline coefficient (CDC), the time from sowing to emergence (Emer), and the time from sowing to maturity (Mat), respectively. This combinatorial approach has been applied successfully in other studies (e.g., Kendall et al. [51]).

2.3.3. Model calibration and Performance Evaluation

The model was calibrated with the data during the field experiment in 2014 and validated in 2015, including ET and canopy cover. In order to evaluate the “goodness-of-fit” of the model, a variety of statistical methods were used, such as the previous study [24,52]. These indicators are computed from the pairs of measured and simulated values Mi and Si (i = 1, 2,…, n) whose means are, respectively, M ¯ and S ¯ , thus:
R 2 = { i = 1 n ( M i M ¯ ) ( S i S ¯ ) [ i = 1 n ( M i M ¯ ) 2 ] 0.5 [ i = 1 n ( S i S ¯ ) 2 ] 0.5 } 2
R M S E = [ i = 1 n ( M i S i ) 2 n ] 0.5
N R M S E = 100 M ¯ [ i = 1 n ( M i S i ) 2 n ] 0.5
E F = 1.0 i = 1 n ( M i S i ) 2 i = 1 n ( M i M ¯ ) 2
d = 1 i = 1 n ( S i M i ) 2 i = 1 n ( | S i M ¯ | + | M i M ¯ | ) 2
It indicates that when d and R2 values are close to 1, the simulated values and the measured values are statistically close. It was reported that when the value of R2 is >0.5, the simulation results are acceptable in the basin simulation. For the value of RMSE (root mean squared error), when the value of RMSE is close to zero, the consistency between simulation results and measured results is good. The value of NRMSE (normalized root mean squared error) may indicate the extent of agreement as a whole and has been applied in several studies. To test the quality of the modelling approach the Nash and Sutcliffe [53] modelling efficiency (EF, non-dimensional) was used. Willmott’s [54] index of agreement (d). However, none of the parameters can be generally accepted by most people as the evaluation criteria. We also recognize the complexity of partitioning evapotranspiration into evaporation and transpiration. Considering only one indicator is not comprehensive enough, so we considered all indicators including d, R2, RMSE, NRMSE, and EF to indicate the efficiency and accuracy of the models used in this study.

3. Results and Discussion

3.1. Differences in Model Parameters Sensitivity and Calibration with and without Mulch

The parameters sensitivity analysis of the AquaCrop model was shown in Figure 2. Since soil evaporation accounts for a small proportion in total ET and is less affected by crops, the values of parameters sensitivity were small and the sensitivity difference was small. Compared with the results of the unmulched treatment, plastic film mulch significantly reduced the sensitivity of 15 parameters, especially the global sensitivity, namely the interaction between parameters. As shown in Figure 2a, Mat is the most sensitive parameter under the condition of plastic film mulch. The GDDays (growing degree days) from sowing to maturity significantly influenced on E through affecting maize growth stage. Plastic film mulch can reduce external environment influence on soil temperature, resulting in conservative GDDays for the whole growth period. As a result, the effect obviously reduced the global sensitivity of the parameters for the impact of E due to the external environment change.
As shown in Figure 2b, for the Tr of seed maize, the sensitivity values of parameters were larger than E due to the large proportion of Tr in total ET, and the sensitivity difference between parameters was large. The most sensitive parameter to Tr was CC0. For seed maize, CC0 was mainly related to planting density and canopy architecture, which was determined by field management and genotype, and was relatively conservative. So, although CC0 was sensitive, it did not need to be adjusted in the later parameter calibration. The CGC was sensitive parameter next to CC0. The CGC determines the shape of the CC curve of seed maize, and the simulation performance of ET by AquaCrop model depends on the simulation of CC curve. Therefore, we tuned the CGC during the parameter calibration process. The BaseTemp was also a relatively sensitive parameter, which affected Tr through the GDDays, thus affecting the whole growth process of seed maize. The sensitivity of most parameters (except CGC) for Tr was reduced by plastic film mulch, which was consistent with the effect of E. These results suggest that plastic film mulch could provide a relatively stable growth environment for seed maize and reduce the interference of the external environment on the growth of seed maize. The improved CGC’s sensitivity by plastic film mulch might be because of ET directly dependence on CC determined by CGC.
The default and calibrated parameters of the AquaCrop model for both M and N treatments were shown in Table 2. The calibration of parameters was based on the above sensitivity analysis. Our study aims to use the AquaCrop model to partition ET on the basis of accurate simulation of ET, so we need to simulate ET as accurately as possible. The simulation accuracy of ET depends on the simulation of CC. Because CGC, CDC, and CCx affect the shape of CC curve and the change process of CC, we gave different values of CDC, CGC, and CCx for each year and each treatment, respectively (Table 2). Plastic film mulch increased CGC from 1.15% GDD−1 to 1.34% GDD−1 in 2014, and from 0.94% GDD−1 to 1.18% GDD−1 in 2015, respectively. The larger the CGC, the faster growth of canopy, i.e., the less time it takes to reach a specific maximum CC. This corresponds to the fact that plastic film mulch can increase GDDays, advance the growth stage and shorten the time required for maturity. For CDC, plastic film mulch slightly improved the CDC, indicating it could accelerate the senescence and maturity of seed maize.
Another important influencing factor is Kc Tr,x for accurately simulate Tr according to Equation (5). Our calibrated Kc Tr,x was 1.25, which was slightly higher than default (1.05) and previous values. The Kc Tr,x was 1.08 for maize in Alpiarca, central Portugal by Paredes et al. [24]; the Kc Tr,x was 1.20 for seed maize in the same experimental site by Ran et al. [41]. In fact, Kc Tr,x is the transpiration coefficient when CC = 100%. For different years and different climatic conditions, the values of Kc Tr,x are different. For our study, the experimental site was the same as Ran et al. [41], but the experimental years were different, so the meteorological factors and other environmental conditions were different, resulting in slightly higher Kc Tr,x due to higher evaporative demand. The Kc Tr,x is similar to the Kcb full in FAO 56, the basic crop coefficient in the middle of growth when the plant completely covers the ground or when LAI > 3. The reference value of Kcb full is 1.20, which should be adjusted to local value according to relative humidity, wind speed, and plant height. The adjusted Kcb full was higher than the default value for dry and hot climatic region. Our Kc Tr,x was close to the adjusted Kcb full, indicating that the calibrated Kc Tr,x was reasonable.
The PstoUpper and PsenUpper were slightly less than the default value. Raes et al. [19] indicated that these two parameters will be reduced correspondingly if there were some stresses. For our study, the research area is located in the arid area of northwest China, where the climate is hot and dry, and there may be some stress in the growth period of seed maize.

3.2. Differences in Canopy Cover and Growth With and without Mulch

As shown in Figure 3, a good agreement was found between the simulated and measured CC for both M and N treatments for the two years. The goodness-of-fit indicators presented higher R2 (>0.9) and lower RMSE of 3.6% and 4.1% in 2014, and 4.1% and 4.0% in 2015, respectively, for the two treatments (Table 3), which was slightly lower than previous studies, e.g., Hsiao et al. [23] [4.8–13.6%] and Zeleke et al. [55] [8.4–12.4%]. Farahani et al. [56] stated that the most logical pathway for a systematic calibration of AquaCrop is first and foremost to ensure a sound prediction of the canopy cover. Because field management could influence CGC and CDC [24,57], it was required to perform distinct parameterizations of the CC curves for different field management. As mentioned above, in order to accurately simulate CC, we used different CGC and CDC values for different years and different treatment methods (Table 2), which was the same as the previous studies on maize and soybean [24,39]. Plastic film mulch advanced the growth stage of seed maize. This was clearly reflected in the length of the growth period, the initial stage for the M treatment was 37–42 days, while it was 46–48 days for the N treatment. These results show that plastic film mulch can accelerate the growth of maize canopy, improve the maximum canopy cover, and increase CCx, thus increasing the effective utilization of water in seed maize.

3.3. Differences in Evapotranspiration with and without Mulch

Although the simulated seasonal ET totals were slightly lower than the measurements for the two treatments of the two years, there was a good agreement between measured and simulated accumulated courses (Figure 4). The bias might be due to errors of measurement calculated using the water balance method [38]. The goodness-of-fit indicators were better, with the linear regression R2 and RMSE of 0.98 and 17.827%, and 0.98 and 19.471% in 2014, 1.00 and 10.740%, and 0.99 and 15.198% in 2015, respectively, for the two treatments (Table 3), which was higher previous R2 of 0.812 in Ran et al. [41]. In 2014, the simulated ET total was 512.6 mm for the M treatment, and 529.2 mm for the N treatment, respectively (Figure 4a,b). In 2015, the simulated ET total was 526.4 mm for the M treatment, and 534.5 mm for the N treatment, respectively (Figure 4c,d). Compared with the total ET of the two treatments in two years, the ET value under the plastic film mulch condition was significantly lower than that under the non-mulched condition. These results indicate that plastic film mulch could decrease ET, which was consistent with previous study [58].

3.4. Differences in Transpiration and Evaporation with and without Mulch

Seasonal variation of daily Tr was similar for the M and N treatments both years, but the M value was slightly higher the N treatment (Figure 5). Daily Tr increased rapidly from less than 1.0 mm⸱d−1 at the initial growth stage to 7.0–8.0 mm⸱d−1 at the mid-season growth stage. Subsequently, daily Tr decreased to less than 1.0 mm⸱d−1 at the end of season for the two treatments in 2014 and 2015. The maximum daily Tr of M treatment is 3.2 mm larger than that of N treatment. E was larger due to lower CC at the initial and late-season stages before the development and after the senescence of maize canopy. After each precipitation or irrigation, there was a significant increase in E. In the mid-season with high CC, E was lower or near to zero. Seasonal variation of daily E was significantly different for the M and N treatments both years. The N E was distinctly higher than the M treatment, especially during the initial growth stage (Figure 5), which was attributed to the reduction of plastic film mulch. According to Equations (12) and (13), soil evaporation was not only dependent on the CC, but also related to the available water in the upper soil (Kr) and the effect of film mulch (fm and fr mulch).
Table 4 presented total and daily averaged Tr, E, and ET during the different growth stages and whole period. Total Tr accounted for most of total ET, with a highest value of 89.31% for the M treatment in 2014. The maximum of total E was 32.80% for the N treatment in 2015. Total Tr during the mid-season stage was the largest proportion, with 98.91% of total ET for the M treatment in 2014. During the initial stage, E accounted for the largest proportion, with 91.65% of total ET for the N treatment in 2014 when precipitation was much. Compared to the N treatment, the ratio of Tr to ET increased from 75.77% to 89.31% in 2014, and 67.20% to 78.12% in 2015 for the M treatment. The lower ratio of Tr to ET in 2015 was due to pest and diseases. In conclusion, in 2014 and 2015, the ET was increased by 33 mm and 25 mm, respectively, by 6.30% and 4.55%.
Higher M Tr indicates that plastic film mulch improved the effective water use of maize or the ratio of Tr to ET, and then biomass production and yield [58]. In another rain-fed maize field, Zhou et al. [59] reported that partly (fm = 0.7) and whole mulching with plastic film reduced total ET by 11.1% and 16.6%, respectively, comparing with no mulching. Averaged ratio of E to ET for the two years was 28.93% under no plastic film mulch, which was close to the ratios of no-mulch maize field reported by other researchers with 26.0% and 30.3% [16,60]. One of the main reasons for ET reduced by mulching was that the proportion of soil evaporation (E) was reduced. According to the Table 4, film mulch shortened the initial and development stages of maize. That is to say, film mulch could shorten the crop season duration and accelerate the growth of maize.
Equation (5) shows that the product of CC* and Kc Tr,x was basal crop coefficient (Kc Tr) as presented in the FAO 56 [38]. We gave the similar table with the FAO 56 basal crop coefficients (Table 17 in the FAO 56). As shown in the red solid line in Figure 6, the M Kc Tr was significantly higher than that for the N treatment, which was consistent with the variation of CC curve (Figure 3). Both Kc Tr of two conditions were lower in the initial stage, and gradually increased with the advance of the growth stage until reaching the maximum in the mid-season. The adjusted Ke was shown in the dotted line in Figure 6. The change trend of Ke was relatively consistent with that of E.
Table 5 presented stage-averaged Kc Tr, Ke, and Kc. Kc Tr was the largest in the mid-season (1.14) and the smallest in the initial stage (0.04). Ke is the largest at the initial stage (0.48) and the smallest at the mid-season stage (0.03). The total Kc reached the maximum at the mid-season stage (1.21) and the minimum at the initial stage (0.38). As shown in Table 5, in the initial stage, there was no difference in Kc Tr between the two treatments, and Ke was reduced by 29.17% with plastic film mulch. In the development stage, the advantages of laminating were gradually shown. Laminating increased Kc Tr by 3.33% and decreased Ke by 63.89%. In the middle-season, the Kc Tr was increased by 0.88% and the Ke was reduced by 62.50% by the plastic film mulch. In the late-season Ke was reduced by 32.26% by plastic film mulch, while for Kc Tr, there were pests and diseases in 2015, and plastic film mulch did not improve Kc Tr.

4. Conclusions

The AquaCrop model performed well for simulating CC and ET, and can be used to partition ET and investigate crop growth of seed maize for treatments with and without plastic film mulch. Plastic film mulch increased the GDDays, changed the phenology, advanced the growth stage, increased beneficial Tr, and reduced E, which was not conducive to seed maize and ET during the growth period of seed maize, resulting in improved water use effectiveness (EUW).

Author Contributions

Designed the experiments and analyzed the data: Q.S. and R.D. Contributed reagents/materials/analysis tools: R.D., T.D., L.T. and S.L. Wrote the manuscript: Q.S., R.D., T.D., L.T. and S.L.

Funding

This research was partially supported by the National Natural Science Foundation of China (51679237, 51790534 and 51621061), the Government Public Research Funds for Projects (201503125), and China Agriculture Research System (CARS-03-35). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Acknowledgments

We thank the two anonymous reviewers for their constructive comments that have helped to improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sinclair, T.R.; Devi, J.; Shekoofa, A.; Choudhary, S.; Sadok, W.; Vadez, V.; Riar, M.; Rufty, T. Limited-transpiration response to high vapor pressure deficit in crop species. J. Cit. Rep. 2017, 260, 109–118. [Google Scholar] [CrossRef] [PubMed]
  2. Lobell, D.B.; Roberts, M.J.; Schlenker, W.; Braun, N.; Little, B.B.; Rejesus, R.M.; Hammer, G.L. Greater sensitivity to drought accompanies maize yield increase in the US Midwest. Science 2014, 344, 516–519. [Google Scholar] [CrossRef]
  3. Leuning, R.; Zhang, Y.Q.; Rajaud, A.; Cleugh, H.; Tu, K. A simple surface conductance model to estimate regional evaporation using MODIS leaf area index and the Penman-Monteith equation. Water Resour. Res. 2008, 44, W10419. [Google Scholar] [CrossRef]
  4. Liu, H.J.; Yang, H.Y.; Zheng, J.H.; Jia, D.D.; Wang, J.; Li, Y.; Huang, G.H. Irrigation scheduling strategies based on soil matric potential on yield and fruit quality of mulched-drip irrigated chili pepper in Northwest China. Agric. Water Manag. 2012, 115, 232–241. [Google Scholar] [CrossRef]
  5. Rana, G.; Katerji, N. Measurement and estimation of actual evapotranspiration in the field under Mediterranean climate: A review. Eur. J. Agron. 2000, 13, 125–153. [Google Scholar] [CrossRef]
  6. 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]
  7. Ding, R.S.; Kang, S.Z.; Li, F.S.; Zhang, Y.Q.; Tong, L. Evapotranspiration measurement and estimation using modified Priestley-Taylor model in an irrigated maize field with mulching. Agric. For. Meteorol. 2013, 168, 140–148. [Google Scholar] [CrossRef]
  8. Pieruschka, R.; Huber, G.; Berry, J.A. Control of transpiration by radiation. In Proceedings of the National Academy of Sciences of The United States of America, Stanford, CA, USA, 16 June 2010; Volume 107, pp. 13372–13377. [Google Scholar]
  9. Er-Raki, S.; Chehbouni, A.; Boulet, G.; Williams, D.G. Using the dual approach of FAO-56 for partitioning ET into soil and plant components for olive orchards in a semi-arid region. Agric. Water Manag. 2010, 97, 1769–1778. [Google Scholar] [CrossRef] [Green Version]
  10. Katerji, N.; Campi, P.; Mastrorilli, M. Productivity, evapotranspiration, and water use efficiency of corn and tomato crops simulated by AquaCrop under contrasting water stress conditions in the Mediterranean region. Agric. Water Manag. 2013, 130, 14–26. [Google Scholar] [CrossRef]
  11. Jones, J.W.; Hoogenboom, G.; Porter, C.H.; Boote, K.J.; Batchelor, W.D.; Hunt, L.A.; Wilkens, P.W.; Singh, U.; Gijsman, A.J.; Ritchie, J.T. The DSSAT cropping system model. Eur. J. Agron. 2003, 18, 235–265. [Google Scholar] [CrossRef]
  12. Diepen, C.A.; Wolf, J.; Keulen, H.; Rappoldt, C. WOFOST: A simulation model of crop production. Soil Use Manag. 1989, 5, 16–24. [Google Scholar] [CrossRef]
  13. Keating, B.A.; Carberry, P.S.; Hammer, G.L.; Probert, M.E.; Robertson, M.J.; Holzworth, D.; Huth, N.I.; Hargreaves, J.N.G.; Meinke, H.; Hochman, Z.; et al. An overview of APSIM, a model designed for farming systems simulation. Eur. J. Agron. 2003, 18, 267–288. [Google Scholar] [CrossRef] [Green Version]
  14. Sharpley, A.N.; Williams, J.R. EPIC-Erosion/Productivity Impact Calculator: 1. Model Documentation; Technical Bulletin-United States Department of Agriculture; Agricultural Research Service: Beltsville, MA, USA, 1990. [Google Scholar]
  15. Ko, J.; Piccinni, G.; Steglich, E. Using EPIC model to manage irrigated cotton and maize. Agric. Water Manag. 2009, 96, 1323–1331. [Google Scholar] [CrossRef]
  16. Kang, S.; Zhang, L.; Liang, Y.; Dawes, W. Simulation of winter wheat yield and water use efficiency in the Loess Plateau of China using WAVES. Agric. Syst. 2003, 78, 355–367. [Google Scholar] [CrossRef]
  17. Stöckle, C.O.; Donatelli, M.; Nelson, R. CropSyst, a cropping systems simulation model. Eur. J. Agron. 2003, 18, 289–307. [Google Scholar] [CrossRef]
  18. Steduto, P.; Hsiao, T.C.; Raes, D.; Fereres, E. AquaCrop—The FAO crop model to simulate yield response to water: I. Concepts and underlying principles. Agron. J. 2009, 101, 426–437. [Google Scholar] [CrossRef]
  19. Raes, D.; Steduto, P.; Hsiao, T.C.; Fereres, E. AquaCrop-the FAO crop model to simulate yield response to water: II. Main algorithms and software description. Agron. J. 2009, 101, 438–447. [Google Scholar] [CrossRef]
  20. Stricevic, R.; Cosic, M.; Djurovic, N.; Pejic, B.; Maksimovic, L. Assessment of the FAO AquaCrop model in the simulation of rainfed and supplementally irrigated maize, sugar beet and sunflower. Agric. Water Manag. 2011, 98, 1615–1621. [Google Scholar] [CrossRef]
  21. Nyakudya, I.W.; Stroosnijder, L. Effect of rooting depth: Plant density and planting date on maize (zea mays L.) yield and water use efficiency in semi-arid Zimbabwe: Modelling with AquaCrop. Agric. Water Manag. 2014, 146, 280–296. [Google Scholar] [CrossRef]
  22. Abedinpour, M.; Sarangi, A.; Rajput, T.B.S.; Singh, M.; Pathak, H.; Ahmad, T. Performance evaluation of AquaCrop model for maize crop in a semi-arid environment. Agric. Water Manag. 2012, 110, 55–66. [Google Scholar] [CrossRef]
  23. Hsiao, T.C.; Heng, L.; Steduto, P.; Rojas-Lara, B.; Raes, D.; Fereres, E. AquaCrop-the FAO crop model to simulate yield response to water: III. Parameterization and testing for maize. Agron. J. 2009, 101, 448–459. [Google Scholar] [CrossRef]
  24. 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]
  25. Mkhabela, M.S.; Bullock, P.R. Performance of the FAO AquaCrop model for wheat grain yield and soil moisture simulation in Western Canada. Agric. Water Manag. 2012, 110, 16–24. [Google Scholar] [CrossRef]
  26. Wang, X.; Wang, Q.; Fan, J.; Fu, Q. Evaluation of the AquaCrop model for simulating the impact of water deficits and different irrigation regimes on the biomass and yield of winter wheat grown on China’s Loess Plateau. Agric. Water Manag. 2013, 129, 95–104. [Google Scholar]
  27. Andarzian, B.; Bannayan, M.; Steduto, P.; Mazraeh, H.; Barati, M.E.; Barati, M.A.; Rahnama, A. Validation and testing of the AquaCrop model under full and deficit irrigated wheat production in Iran. Agric. Water Manag. 2011, 100, 1–8. [Google Scholar] [CrossRef]
  28. Iqbal, M.A.; Shen, Y.; Stricevic, R.; Pei, H.; Sun, H.; Amiri, E.; Penas, A.; Rio, S. Evaluation of the FAO AquaCrop model for winter wheat on the North China Plain under deficit irrigation from field experiment to regional yield simulation. Agric. Water Manag. 2014, 135, 61–72. [Google Scholar] [CrossRef]
  29. Toumi, J.; Er-Raki, S.; Ezzahar, J.; Khabba, S.; Jarlan, L.; Chehbouni, A. Performance assessment of AquaCrop model for estimating evapotranspiration, soil water content and grain yield of winter wheat in Tensift Al Haouz (Morocco): Application to irrigation management. Agric. Water Manag. 2016, 163, 219–235. [Google Scholar] [CrossRef]
  30. Linker, R.; Ioslovich, I.; Sylaios, G.; Plauborg, F.; Battilani, A. Optimal model-based deficit irrigation scheduling using AquaCrop: A simulation study with cotton, potato and tomato. Agric. Water Manag. 2016, 163, 236–243. [Google Scholar] [CrossRef]
  31. Garcia-Vila, M.; Fereres, E.; Mateos, L.; Orgaz, F.; Steduto, P. Deficit irrigation optimization of cotton with AquaCrop. Agron. J. 2009, 101, 477–487. [Google Scholar] [CrossRef]
  32. Maniruzzaman, M.; Talukder, M.S.U.; Khan, M.H.; Biswas, J.C.; Nemes, A. Validation of the AquaCrop model for irrigated rice production under varied water regimes in Bangladesh. Agric. Water Manag. 2015, 159, 331–340. [Google Scholar] [CrossRef]
  33. Todorovic, M.; Albrizio, R.; Zivotic, L.; Saab, M.T.A.; Stöckle, C.; Steduto, P. Assessment of AquaCrop, CropSyst, and WOFOST models in the simulation of sunflower growth under different water regimes. Agron. J. 2009, 101, 509–521. [Google Scholar] [CrossRef]
  34. Araya, A.; Keesstra, S.D.; Stroosnijder, L. Simulating yield response to water of Teff (Eragrostis tef) with FAO’s AquaCrop model. Field Crops Res. 2010, 116, 196–204. [Google Scholar] [CrossRef]
  35. Dong, H.Z.; Li, W.J.; Tang, W.; Zhang, D.M. Early plastic mulching increases stand establishment and lint yield of cotton in saline fields. Field Crops Res. 2009, 111, 269–275. [Google Scholar] [CrossRef]
  36. Zhao, W.; Liu, B.; Zhang, Z. Water requirements of maize in the middle Heihe River basin, China. Agric. Water Manag. 2010, 97, 215–223. [Google Scholar] [CrossRef]
  37. Hou, X.Y.; Wang, F.X.; Han, J.J.; Kang, S.Z.; Feng, S.Y. Duration of plastic mulch for potato growth under drip irrigation in an arid region of Northwest China. Agric. For. Meteorol. 2010, 150, 115–121. [Google Scholar] [CrossRef]
  38. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop evapotranspiration. guide-lines for computing crop water requirements. In FAO Irrigation and Drainage Paper 56; Food and Agriculture Organization of United Nations: Rome, Italy, 1998. [Google Scholar]
  39. Paredes, P.; Wei, Z.; Liu, Y.; Xu, D.; Xin, Y.; Zhang, B.; Pereira, L.S. Performance assessment of the FAO AquaCrop model for soil water, soil evaporation, biomass and yield of soybeans in North China plain. Agric. Water Manag. 2015, 152, 57–71. [Google Scholar] [CrossRef]
  40. Pereira, L.S.; Paredes, P.; Rodrigues, G.C.; Neves, M. Modeling malt barley water use and evapotranspiration partitioning in two contrasting rainfall years. Assessing AquaCrop and SIMDualKc models. Agric. Water Manag. 2015, 159, 239–254. [Google Scholar] [CrossRef]
  41. Ran, H.; Kang, S.Z.; Li, F.S.; Tong, L.; Ding, R.S.; Du, T.S.; Li, S.E.; Zhang, X.T. Performance of AquaCrop and SIMDualKc models in evapotranspiration partitioning on full and deficit irrigated maize for seed production under plastic film-mulch in an arid region of China. Agric. Syst. 2017, 151, 20–32. [Google Scholar] [CrossRef]
  42. Campolongo, F.; Cariboni, J.; Saltelli, A. An effective screening design for sensitivity analysis of large models. Environ. Model. Softw. 2007, 22, 1509–1518. [Google Scholar] [CrossRef]
  43. Confalonieri, R.; Bellocchi, G.; Tarantola, S.; Acutis, M.; Donatelli, M.; Genovese, G. Sensitivity analysis of the rice model WARM in Europe: Exploring the effects of different locations, climates and methods of analysis on model sensitivity to crop parameters. Environ. Model. Softw. 2010, 25, 479–488. [Google Scholar] [CrossRef]
  44. Saltelli, A.; Ratto, M.; Andres, T.; Campolongo, F.; Cariboni, J.; Gatelli, D.; Saisana, M.; Tarantola, S. Global Sensitivity Analysis: The Primer; John Wiley: Chichester, UK, 2008; p. 292. [Google Scholar]
  45. Saltelli, A.; Annoni, P. How to avoid a perfunctory sensitivity analysis. Environ. Model. Softw. 2010, 25, 1508–1517. [Google Scholar] [CrossRef]
  46. Cariboni, J.; Gatelli, D.; Liska, R.; Saltelli, A. The role of sensitivity analysis in ecological modelling. Ecol. Model. 2007, 203, 167–182. [Google Scholar] [CrossRef]
  47. Elsawwaf, M.; Willems, P.; Feyen, J. Assessment of the sensitivity and prediction uncertainty of evaporation models applied to Nasser Lake, Egypt. J. Hydrol. 2010, 395, 10–22. [Google Scholar] [CrossRef]
  48. Saltelli, A.; Tarantola, S.; Campolongo, F. Sensitivity analysis as an ingredient of modeling. Stat. Sci. 2000, 15, 377–395. [Google Scholar]
  49. Morris, M.D. Factorial sampling plans for preliminary computational experiments. Technometrics 1991, 33, 161–174. [Google Scholar] [CrossRef]
  50. Sobol, I.M. Sensitivity estimates for nonlinear mathematical models. Math. Model. Comput. Exp. 1993, 1, 407–414. [Google Scholar]
  51. Kendall, C.D.; James, C.A. II.; Mehdi, A.; Allan, A.A.; Mazdak, A. Global sensitivity and uncertainty analysis of a dynamic agroecosystem model under different irrigation treatments. Ecol. Model. 2012, 231, 113–125. [Google Scholar]
  52. 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]
  53. Nash, J.E.; Sutcliffe, J.V. River flow forecasting through conceptual models part I—A discussion of principles. J. Hydrol. 1970, 10, 282–290. [Google Scholar] [CrossRef]
  54. Willmott, C.J. On the validation of models. Phys. Geogr. 1981, 2, 184–194. [Google Scholar] [CrossRef]
  55. Zeleke, K.T.; Luckett, D.; Cowley, R. Calibration and testing of the FAO AquaCrop model for canola. Agron. J. 2011, 103, 1610–1618. [Google Scholar] [CrossRef]
  56. Farahani, H.J.; Izzi, G.; Oweis, T.Y. Parameterization and evaluation of the AquaCrop model for full and deficit irrigated cotton. Agron. J. 2009, 101, 469–476. [Google Scholar] [CrossRef]
  57. Abi Saab, N.T.; Todorovic, M.; Albrizio, R. Comparing AquaCrop and CropSyst models in simulating barley growth and yield under different water and nitrogen regimes. Does calibration year influence the performance of crop growth models? Agric. Water Manag. 2015, 147, 21–33. [Google Scholar] [CrossRef]
  58. Fan, Y.Q.; Ding, R.S.; Kang, S.Z.; Hao, X.M.; Du, T.S.; Tong, L. 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]
  59. Zhou, L.M.; Li, F.M.; Jin, S.L.; Song, Y.J. 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]
  60. Liu, C.; Zhang, X.; Zhang, Y. Determination of daily evaporation and evapotranspiration of winter wheat and maize by large-scale weighing lysimeter and micro-lysimeter. Agric. For. Meteorol. 2002, 111, 109–120. [Google Scholar] [CrossRef]
Figure 1. Variations in daily meteorological data against days after planting (DAP) of seed maize in 2014 and 2015, including maximum (Tmax, —) and minimum (Tmin, – – – –) air temperature, precipitation (▐), and reference evapotranspiration (ET0, ––◦––).
Figure 1. Variations in daily meteorological data against days after planting (DAP) of seed maize in 2014 and 2015, including maximum (Tmax, —) and minimum (Tmin, – – – –) air temperature, precipitation (▐), and reference evapotranspiration (ET0, ––◦––).
Water 11 01153 g001
Figure 2. The Sobol’ global sensitivity analysis (GSA, M, Water 11 01153 i001; N, Water 11 01153 i002), and first order sensitivity analysis (SA, M, Water 11 01153 i003; N, Water 11 01153 i004) indices for the AquaCrop model output of evaporation (E) and transpiration (Tr). Model input parameters were those sensitive by the Morris screening method.
Figure 2. The Sobol’ global sensitivity analysis (GSA, M, Water 11 01153 i001; N, Water 11 01153 i002), and first order sensitivity analysis (SA, M, Water 11 01153 i003; N, Water 11 01153 i004) indices for the AquaCrop model output of evaporation (E) and transpiration (Tr). Model input parameters were those sensitive by the Morris screening method.
Water 11 01153 g002
Figure 3. Simulated canopy cover (CC) using default (- - - - -) and calibrated (—) parameters and measurements (●) of seed maize in 2014 and 2015.
Figure 3. Simulated canopy cover (CC) using default (- - - - -) and calibrated (—) parameters and measurements (●) of seed maize in 2014 and 2015.
Water 11 01153 g003
Figure 4. Simulated accumulated ET using default (- - - - -) and calibrated (—) parameters and measurements (●) of seed maize in 2014 and 2015.
Figure 4. Simulated accumulated ET using default (- - - - -) and calibrated (—) parameters and measurements (●) of seed maize in 2014 and 2015.
Water 11 01153 g004
Figure 5. Variations in daily evaporation (E, - - - - - - -) and transpiration (Tr, —) against DAP of seed maize for treatments with and without plastic film mulch (M and N) in 2014 and 2015.
Figure 5. Variations in daily evaporation (E, - - - - - - -) and transpiration (Tr, —) against DAP of seed maize for treatments with and without plastic film mulch (M and N) in 2014 and 2015.
Water 11 01153 g005
Figure 6. Variations in daily crop coefficient components against DAP of seed maize partitioned for treatments with and without plastic film mulch (M and N) in 2014 to 2015, including crop coefficient Kc (- - - - - - -), soil evaporation coefficient Ke (–⸱⸱–⸱⸱–⸱⸱–) and transpiration coefficient Kc Tr (—).
Figure 6. Variations in daily crop coefficient components against DAP of seed maize partitioned for treatments with and without plastic film mulch (M and N) in 2014 to 2015, including crop coefficient Kc (- - - - - - -), soil evaporation coefficient Ke (–⸱⸱–⸱⸱–⸱⸱–) and transpiration coefficient Kc Tr (—).
Water 11 01153 g006
Table 1. Irrigation schedules for treatments with and without plastic film mulch (M and N) and precipitation during the growth season of seed maize in 2014 and 2015.
Table 1. Irrigation schedules for treatments with and without plastic film mulch (M and N) and precipitation during the growth season of seed maize in 2014 and 2015.
YearsTreatmentsIrrigation Schedules and Irrigation Depth (mm)Irrigation (I)Precipitation (P)I + P
6/86/136/277/77/127/197/257/318/68/158/208/219/2(mm)(mm)(mm)
2014M 80.6 80.6 80.6 30.7 272.5241.0513.5
N 86.4 98.6 98.6 53.8 337.4241.0578.4
2015M69.1 53.6 72.9 87.2 51.9 90.6425.3150.6575.9
N69.1 75.7 77.0 78.5 72.2 102.9475.4150.6626.0
Table 2. Default and calibrated values of parameters of seed maize in the AquaCrop model.
Table 2. Default and calibrated values of parameters of seed maize in the AquaCrop model.
DescriptionDefault ValuesCalibrated Values
Base temperature, °C88
Upper temperature, °C3030
Leaf growth threshold (PexpUpper)0.140.14
Leaf growth threshold (PexpLower)0.720.72
Leaf growth stress coefficient curve shape (SFexp)2.92.9
Minimum effective rooting depth (MinRootDepth), m0.30.3
Maximum effective rooting depth (MaxRootDepth), mUp to 2.800.8
Senescence stress coefficient curve shape (SFsen)2.72.7
Canopy cover per seedling at 90% emergence (CC0), cm2 per plant6.56.5
Stomatal conductance threshold (PstoUpper)0.690.45
Senescence stress coefficient (PsenUpper)0.690.5
Crop coefficient for transpiration at CC = 100% (KcTr,x)1.051.25
Time from sowing to emergence (Emer), GDD, °C d50–10030–77
Time from sowing to maturity (Mat), GDD, °C dTSE a + 1450 − 17501551–1691
14M14N15M15N
Canopy growth coefficient (CGC), % (°C d)−11.2–1.31.341.151.180.94
Canopy decline coefficient (CDC), % (°C d)−11.00.290.290.680.66
Maximum canopy cover (CCx), %65–9996899087
a TSE is the time from sowing to emergence. Bold values are the specific calibrated parameters in this study.
Table 3. Goodness-of-fit indicators between measurements and simulations for CC and accumulated evapotranspiration (ET) of seed maize for treatments with and without plastic film mulch (M and N) in 2014 and 2015.
Table 3. Goodness-of-fit indicators between measurements and simulations for CC and accumulated evapotranspiration (ET) of seed maize for treatments with and without plastic film mulch (M and N) in 2014 and 2015.
TreatmentsNR2RMSE (%/mm)NRMSE (%)EFd
CC (%)
2014 M101.003.65.50.991.00
2014 N100.994.17.40.991.00
2015 M110.994.16.10.981.00
2015 N111.004.06.70.981.00
ET (mm)
2014 M120.9817.8278.10.991.00
2014 N120.9819.4718.30.991.00
2015 M211.0010.7404.41.001.00
2015 N210.9915.1986.30.991.00
Table 4. ET, Tr, E, and the percentage of Tr and E to ET in different growth stages for treatments with and without plastic film mulch (M and N) in 2014 and 2015.
Table 4. ET, Tr, E, and the percentage of Tr and E to ET in different growth stages for treatments with and without plastic film mulch (M and N) in 2014 and 2015.
TreatmentsStageDAPDaysTr (mm)Daily Tr (mm⸱d−1)E (mm)Daily E (mm⸱d−1)ET (mm)Daily ET (mm⸱d−1)Tr/ETE/ET
14 MInitial1–42 d426.30.15 44.71.06 511.21 12.35%87.65%
Development43–57 d15453.00 4.90.33 49.93.33 90.18%9.82%
Mid-season58–150 d93391.14.21 4.30.05 395.44.25 98.91%1.09%
Late-season151–158 d815.41.93 0.90.11 16.32.04 94.48%5.52%
Whole1–158 d158457.82.90 54.80.35 512.63.24 89.31%10.69%
14 NInitial1–48 d487.50.16 82.31.71 89.81.87 8.35%91.65%
Development49–68d2048.22.41 24.41.22 72.63.63 66.39%33.61%
Mid-season69–147 d79323.64.10 18.70.24 342.34.33 94.54%5.46%
Late-season148–158 d1121.71.97 2.80.25 24.52.23 88.57%11.43%
Whole1–158 d1584012.54 128.20.81 529.23.35 75.77%24.23%
15 MInitial1–37d375.90.16 54.41.47 60.31.63 9.78%90.22%
Development38–59 d2258.12.64 18.90.86 773.50 75.45%24.55%
Mid-season60–126 d673014.49 12.30.18 313.34.68 96.07%3.93%
Late-season127–156d3046.21.54 29.60.99 75.82.53 60.95%39.05%
Whole1–156 d156411.22.64 115.20.74 526.43.37 78.12%21.88%
15 NInitial1–46 d468.40.18 72.51.58 80.91.76 10.38%89.62%
Development47–71 d2561.52.46 40.71.63 102.24.09 60.18%39.82%
Mid-season72–124 d532364.45 20.20.38 256.24.83 92.12%7.88%
Late-season125–156 d3253.31.67 41.91.31 95.22.98 55.99%44.01%
Whole1–156 d156359.22.30 175.31.12 534.53.43 67.20%32.80%
Table 5. Mean crop coefficient Kc, transpiration coefficient Kc Tr, soil evaporation coefficient Ke of seed maize in different growth stages for treatments with and without plastic film mulch (M and N).
Table 5. Mean crop coefficient Kc, transpiration coefficient Kc Tr, soil evaporation coefficient Ke of seed maize in different growth stages for treatments with and without plastic film mulch (M and N).
ComponentsTreatmentsInitialDevelopmentMid-SeasonLate-Season
KeM0.34 0.13 0.03 0.21
N0.48 0.36 0.08 0.31
Kc TrM0.04 0.62 1.14 0.58
N0.04 0.60 1.13 0.63
KcM0.38 0.76 1.17 0.79
N0.52 0.96 1.21 0.94

Share and Cite

MDPI and ACS Style

Shen, Q.; Ding, R.; Du, T.; Tong, L.; Li, S. Water Use Effectiveness Is Enhanced Using Film Mulch Through Increasing Transpiration and Decreasing Evapotranspiration. Water 2019, 11, 1153. https://doi.org/10.3390/w11061153

AMA Style

Shen Q, Ding R, Du T, Tong L, Li S. Water Use Effectiveness Is Enhanced Using Film Mulch Through Increasing Transpiration and Decreasing Evapotranspiration. Water. 2019; 11(6):1153. https://doi.org/10.3390/w11061153

Chicago/Turabian Style

Shen, Qianxi, Risheng Ding, Taisheng Du, Ling Tong, and Sien Li. 2019. "Water Use Effectiveness Is Enhanced Using Film Mulch Through Increasing Transpiration and Decreasing Evapotranspiration" Water 11, no. 6: 1153. https://doi.org/10.3390/w11061153

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