Next Article in Journal
Deriving an Optimal Operation Plan for Hydraulic Facilities with Complex Channels through Unsteady Flow Simulations
Next Article in Special Issue
Analysis and Design of Monopile Foundations for Offshore Wind and Tidal Turbine Structures
Previous Article in Journal
Backward Trajectory Analysis Based on Specific Humidity Correction for the Influence of Moisture Sources on Precipitation Isotopes in the Western Loess Plateau, China
Previous Article in Special Issue
Investigation on the Impact of Different Absorber Materials in Solar Still Using CFD Simulation—Economic and Environmental Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Sensitivity Analysis and Determination of the Optimal Level of Water Use Efficiency for Winter Wheat and Barley under Different Irrigation Scenarios Using the AquaCrop Model in Arid and Semiarid Climatic Conditions (Case Study: Dehloran Plain, Iran)

Department of Water Science Engineering, Islamic Azad University of Arak Branch, Arak 38135-567, Iran
*
Author to whom correspondence should be addressed.
Water 2022, 14(21), 3455; https://doi.org/10.3390/w14213455
Submission received: 26 September 2022 / Revised: 23 October 2022 / Accepted: 25 October 2022 / Published: 29 October 2022
(This article belongs to the Special Issue Renewable Energy Systems Flexibility for Water Desalination)

Abstract

:
The AquaCrop model is one of the most recent models that can simulate the growth rates and yields of various crops based on water consumption levels. To determine the optimal irrigation level, data measured in two crop years (2018–2019) in different irrigation scenarios (full irrigation or 100% water requirement and 90, 80, 70, 60, and 50% irrigation) were calibrated and validated for arid and semiarid climatic conditions using the AquaCrop model. The model was calibrated to simulate the grain yields of winter wheat and barley using R2, RMSE, d, EF, and NRMSE statistical indicators. The obtained values of these indices were, respectively, 0.97, 3, 0.98, 0.94, and 4 for winter wheat and 0.98, 4, 0.92, 0.89, and 7 for barley. The model efficiency was also validated using crop harvest data in the crop year 2019. For grain yield simulation, the calculated values of R2, RMSE, d, EF, and NRMSE statistical indicators were, respectively, 0.99, 4, 0.97, 0.93, and 4.4 for winter wheat and 0.97, 7, 0.94, 0.91, and 9 for barley. The data of field and modeled samples were analyzed by analysis of variance (ANOVA) using the F-test, and significant results were obtained for both crops in all applied scenarios at the 95% level.

1. Introduction

Nowadays, pioneer scientists around the world are seeking to quantitatively examine plant growth and performance using mathematical models and represent the results of existing research practically in broader dimensions. The growth process of crops can be studied at different levels of accuracy and detail. Since various field survey studies and experiments require costly and time-consuming procedures, computer sciences have made it possible to discover the relationships between growth processes and modeling. As such, the growth rate and yield of various crops can be estimated in variable climatic and managerial situations. These models reduce geographical and environmental limitations and can be generalized to various crop cultivars.
Modeling makes it possible to investigate quantitative and reasonable viewpoints about concurrent processes that interact with each other. The growth process and performance of plants can be simulated using theories proposed in various studies, thereby examining their accuracy and richness. It is noteworthy that models are necessarily a simple demonstration of the reality and behavior of the real system and, therefore, cannot completely express the behavior of the real system in real situations. Accordingly, the AquaCrop model was developed by the Food and Agriculture Organization (FAO) of the United Nations, which is a general conceptual model that strikes a balance between simplicity, accuracy, and strength. Compared to other models, the AquaCrop model is simple to apply and allows for crop yield simulations in multiple scenarios of field management. It is a water-based model that has primarily been developed to simulate crop responses to water and irrigation management strategies. The model requires a limited set of input parameters, most of which are easily accessible [1].
The AquaCrop model has been in use in different regions of the world with success. A brief review of its application for various crops in different scenarios of management and operations showed that one of the main applications of this model is to simulate different scenarios of water irrigation [2]. The AquaCrop model was used to simulate the irrigation management of winter wheat in an arid region [3]. Researchers conducted an experiment with five irrigation strategies to evaluate the water stress tolerance of wheat [4]. Model predictions of the root-zone soil water content (SWC), coverage (CC), grain yield (GY), and underground biomass (BM) produced good results [5,6]. In a study involving a two-year experiment, the authors calibrated the AquaCrop model and examined its accuracy in simulating barley crop indices in the Darab region. The model was calibrated by comparing the results of field studies and simulations. The amount of dry matter in the aerial part of the simulated plant and the amount of grain yield predicted by the model were investigated with d and NRSME statistical indices in the first and second years of cultivation. The results of statistical tests showed that the AquaCrop model was highly accurate [7,8]. Another study was carried out to assess water productivity under conditions of water stress (40% water-holding capacity) compared with normal conditions (75% water-holding capacity) as predicted by the AquaCrop model [9,10]. In a study entitled Simulating Soil Water Content, Evapotranspiration, and Yield of Variably Irrigated Grain Sorghum Using AquaCrop in the Aqualala Falls Aquifer Area, the authors concluded that the overall performance of the AquaCrop model shows which technique can be used as an effective tool to assess the effects of variable irrigation levels on grain yield within the study area. The effect of the planting density was found to be negligible [11]. Energy deficit improvements should be made using appropriate economic analysis. Increasing the use of renewable energy just to eliminate energy shortages may not be economically justified [12,13]. The crop growth trends were observed to be related to the cover crop canopy instead of the leaf area index. This was a major discovery for the development and application of the model. However, the model may also lead to some inaccuracies in predicting water balance. In effect, the AquaCrop model is suitable for supporting crop production but not for predicting changes in soil moisture and evaporation under drip irrigation [14]. An evaluation of the model was performed using the coefficient of the root-mean-square error (RMSE), normalized RMSE, and R2. The AquaCrop model performed well in simulating grain yield and final biomass production, with R2 > 0.90 and RMSE and normalized RMSE values less than 10 [15,16]. With calibration and validation, the AquaCrop and APSIM models can be used to derive the best management strategies in terms of N fertilizer and the water regime for wheat conditions [17]. Model analysis was carried out for the purpose of investigating opportunities to maximize wheat yield and resource efficiency, including irrigation water use efficiency (WUE), nitrogen use efficiency (NUE), and solar radiation use efficiency (RUE) [18]. The authors reported an even greater wheat water productivity of 12.49 kg ha−1 mm−1 with a yield threshold of 132 mm in northeast Colorado. The difficulty in measuring the components of the soil water balance encourages the use of simulation models to investigate the processes involved [19]. Crop-water production functions, which express the relationship between the crop yield (biomass yield, grain yield, and lint yield) and crop water use during the growing season, are important tools for quantifying the effects of water scarcity on agricultural production [20].
Therefore, the main objectives of this study can be summarized as follows:
(i)
To evaluate the AquaCrop model for simulating the yield responses of winter wheat and barley at different irrigation levels by comparing the model results with those of field surveys.
(ii)
To calibrate the AquaCrop models for winter wheat and barley in a large agricultural area in Ilam Province in central Iran.
The evaluation results can help in adopting optimal irrigation management strategies for winter wheat and barley in arid and semiarid climatic conditions.

2. Materials and Methods

2.1. AquaCrop Model Structure

This model is the most recent model for plant growth simulation based on water consumption [1,21]. It requires fewer parameters compared to other models to simulate plant growth. The basic equation used is that described by Doorenbos and Kassam in 1979, published in FAO Publication No. 33. This method has been the main reference for the analysis of the response to the amount of water used in agriculture for more than 20 years. Although the structure of the model is simple and requires few parameters, the results of research have shown that this model produces highly accurate results and is powerful in simulating the growth and yield of plants [1].
1 Y Y x = K y   ( 1 E T E T x )
where Yx and Y are the maximum and actual yields, (1 − Y/Yx) is the relative yield decline, ETx and ET are the maximum and actual evapotranspiration, (1 − ET/ETx) is the relative water stress, and Ky is the proportionality factor between the relative yield decline and the relative reduction in evapotranspiration.
The fundamental changes introduced in Equation (1) for the simulation of plant growth and yield by AquaCrop models are:
  • The consideration of the effect of the harvest index (HI) for the calculation of the final biomass yield estimate.
  • The consideration of the amount of evapotranspiration separately from soil evaporation (Es), transpiration (Ta), and ground crevices with the final performance calculation. Since there is little plant cover for plant growth in the first stage, the amount of evaporation from the soil surface is significant, and it is not necessary to consider when calculating the amount of water consumed by plants. In this model, the rate of evaporation from the soil surface is calculated using Ritchie’s equation (1972).
The model inputs include four categories: climate, plants, soil, and management. The input data in the region include the maximum and minimum air temperature (°C), the amount of rainfall (mm·day−1), reference evapotranspiration (ETo) (mm·day−1), and the average annual concentration of carbon dioxide (CO2) (ppm).
The input data can be those obtained on a daily basis, ten-day basis, or monthly basis. To obtain the reference evapotranspiration calculated from these parameters, the ETo Calculator was used. The data required by this calculator include the minimum and maximum temperatures, relative humidity, wind speed, and the number of hours of sunshine.
The most important input data in the plant category include: 1. the type of plant (fruit or seed crops, leafy vegetables, or root or tuberous plants and C3 or C4); 2. the planting method (seed or seedling); 3. the planting date; 4. the growth period (days); 5. plant development information (early-stage vegetation (CCo), plant celestial development, flowering and seed formation, plant root depth development, temperature, and plant density estimation); 6. plant information in conditions without water stress, fertilizer, and salinity, which includes water use efficiency data (WP*) and the crop harvest index (HI), where the reference harvest index (HIo) is presented as a representative of the harvest index; 7. water (coefficients and the pattern of water absorption); 8. the pattern of water extraction by roots; and 9. the parameters of water irrigation under stress conditions. Data entry for each category above was in accordance with the defined projects.

2.2. The study Area

The study area of the Dehloran plain is 54,755 hectares in Dehloran city, located in Ilam Province (with a longitude of 47′17 to 47′22 and a latitude of 32′32 to 32′37). In this region, the long-term average rainfall is 262 mm, with a maximum monthly rainfall of 50 mm of precipitation in January and a minimum rainfall of 0 mm during the summer. The average maximum monthly temperature in this region is 42.6 °C in August, and the average minimum monthly temperature is 8.3 °C in January.
To validate the model, highly transparent data were used from twenty observed fields (Figure 1).

2.3. Plant Data

The plant data and crop characteristics include constant parameters and user-specific data. The values of constant plant parameters do not change over time or among geographic situations and are present in the model by default for important crops, including wheat. These constant parameters are calibrated with no limitations using plant growth data in favorable conditions and are used with the stress response factor under water stress conditions. In addition to constant parameters, some other data required for the simulation are dependent on the crop species and/or have different values according to various managerial and environmental conditions. Some of these parameters include the cultivation date and density, the timing of plant phenological stages, and maximum root depth, which are variable depending on the specific conditions of each region and crop and are determined by the user (Table 1).

2.4. Irrigation and Field Management

These data include two parts, irrigation and field management, and the latter consists of selecting the soil fertility level, the method of soil water balance, such as mulch to reduce soil evaporation or mounds to save water in the field, and the plowing method. For irrigation management, rain-fed or irrigated crops are selected by the user. An irrigation system was selected in this research, and part of the wet and defined level was applied to detect the water quality, time, and level of irrigation. Basin irrigation was selected based on the commonly used irrigation method in the region. No mounds or mulch was used in field management. Six irrigation treatments, namely, full irrigation (100% (T100)) and 90% (T90), 80% (T80), 70% (T70), 60% (T60), and 50% (T50) water requirements (Table 2). The experimental design used in this research is considered a completely randomized block design; the output of Cropwat results is annual, but the calculations pertain to 2 years, which exactly corresponds to the 2-year period of crop cultivation.

2.5. Meteorological Data

In the AquaCrop model, the required meteorological data include precipitation, minimum/maximum temperatures, reference evapotranspiration, and atmospheric CO2 concentration. Data obtained from the Ilam meteorological station were used in this study. The default atmospheric CO2 concentration from 1902 to 2099 was used in the model. The reference plant evapotranspiration was calculated using the ETo Calculator tool in the AquaCrop model [1] (Table 3).

2.6. Soil Properties

In this section, the required data include the texture and layers of soil, electrical conductivity (EC) of saturated soil, volumetric water content, field capacity, and permanent wilting point. The model can estimate the EC of unsaturated soil using soil water content data (Table 4).

2.7. Model Calibration

Calibration involves changing some of the parameters and coefficients of the model so that the data simulated by the model are in acceptable agreement with observed or field data. In this study, the model was calibrated based on data from the first year of cultivation (2018), and output data simulated by the model were compared with observed data or actual changes and expressed as a model coefficient. This was carried out until the observed and simulated data had acceptable error estimates. Climate data measured for plants, soil, and farm management were entered into the model. Then, some of the best-fit model coefficients between observed and simulated yield values were reduced or increased.
One of the critical factors in plant model calibration is the calibration of the canopy cover (CC). Furthermore, the growth rate of the crop shell or cover is calculated by the model after entering the phenological dates of the plant, such as the time of emergence, the time when reaching the maximum canopy coverage (CCx), the time of the onset of aging, and the maturity and ripening of the crop. The canopy growth coefficient (CGC), loss coverage ratio (CGC), and water stress profiles are the main factors in estimating vegetation. The parameters of the water stress curve in relation to default data within the model can be changed to best fit observed and simulated data (Table 5).

2.8. Model Validation and Performance Evaluation

In the validation phase, the calibration model performance was evaluated. For this study, field data from the second year of cultivation (2019) were used. The output data from the model were compared with the observed data. To evaluate the performance of the model, statistical methods, such as the best fit of the regression line between observed and simulated data and comparison with the 1:1 line, the coefficient of determination (R2), root-mean-square error (RMSE), normalized root-mean-square error (NRMES), and Wilmot index, were used.
The statistical indicators of the quality of the model (EF) are dimensionless, and an indicator of the relative magnitude of the measured residual variance with a value close to zero or negative indicates that the model provides a good prediction from the observed data.
Given that the RMSE value depends on the unit of measurement, the evaluation of different models for two variables with different units will not be correct. Hence, the RMSE values are divided into dependent variable ranges. This is named the normalized RMSE (NRSME). This criterion is appropriate for the comparison of different models. To compare field data with those of the model, the testing of the model parameters and the output results is required. For this purpose, and in order to compare overlapping models, we used the F-test with a degree of freedom of n − 1 and a significance level of 95%. In the F-test, when the absolute computational value is less than the critical absolute value, we can accept the null hypothesis (H0) and reject hypothesis 1 (H1).
RMSE = ( 1 n i = 1 n ( S i O i ) 2 )
NRMSE = 1 O a v g ( 1 n i = 1 n ( S i O i ) 2 ) × 100
EF = 1 i = 1 n ( O i S i ) 2 i = 1 n ( O i O a v g ) 2
d = 1 i = 1 n ( S i O i ) 2 i = 1 n ( | S i O a v g | + | O i O a v g | ) 2
where “n” is the number of data, “Si” is the simulated data, “Oi” is the observed data, and “Oavg” is the observed average data.

2.9. Sensitivity Analysis of Input Data

The simulated results were subjected to sensitivity analysis to determine the effect of the input data used in the model calibration test. Grain yield, biomass, vegetation, and plant water use efficiency (WUE) were considered the input data (Table 6). Accordingly, the model input data were categorized according to their effects on the model inputs based on the degree of sensitivity (high, moderate, and low). In this study, three classes were used: the first is a 15% change in the output values of the model, defined as very sensitive data; the second is a change from 5% to 15% in the output values, defined as semi-sensitive; and the third is a change of less than 5% in the output data, considered data with a low sensitivity coefficient [1,18,21].

3. Results

3.1. Model Performance Evaluation

The simulation results for grain yield, biomass, and water use efficiency produced by the model compared with the data observed in the calibration test using different irrigation levels for winter wheat plants are listed in Table 6. Table 7 lists the same values for winter barley plants. The results show that the minimum and maximum errors for winter wheat are 1.90% when estimating the grain yield with 90% irrigation and 4.92% when estimating water use efficiency with 50% deficit irrigation. For barley, the minimum and maximum errors are 1.40% for the estimation of water use efficiency with 80% irrigation and 4.90% for the grain yield estimation with 80% deficit irrigation.
Using the functions represented in Figure 1 and Figure 2, the test calibration of the model performance to simulate grain yield resulted in values of R2, RMSE, d, EF, and NRMSE equal to 0.97, 3, 0.98, 0.94, and 4 for winter wheat, respectively, and 0.98, 4, 0.92, 0.89, and 7 for barley; thus, it can be concluded that the model is ideal. According to the results, it can be concluded that the accuracy of the model decreases with increasing water stress. The calibration results of the model for simulating grain yield, biomass, and water use efficiency under different irrigation treatments are presented in Figure 2 for the winter wheat crop and in Figure 3 for the barley crop.
Production performance curves were generated to determine the irrigation water productivity index, and the performance of the functions were used. The functions of the relationship between productivity and the water used during the irrigation season are shown. When interpreting the points on the curve and production functions, a limitation to bear in mind is that the equation is true only in a certain range. In addition, each crop of any plant in any climate has its own function, and it is necessary to apply the results only to the same climate or the same conditions, with each plant having its own curve.
Figure 4 illustrates the production yield of wheat, and Figure 5 shows the production yield of barley with the amount of water applied to the land under different irrigation treatments. For each scenario, the required water depth was measured through Microflumes installed in the field. The water requirements of winter wheat and barley were calculated to be 415 and 391.3 mm, respectively. The levels of deficit irrigation applied to the plant during different growth periods were calculated using the product method [22]. As accurately represented by the curves, the slopes are initially high and then gradually decrease. This implies that water usage efficiency is higher for deficit irrigation strategies. When increasing the amount of irrigation, the corresponding curve finally reaches its peak. After that, further irrigation does not significantly change productivity. Due to the use of water during the germination and early growth of the plant, the curves do not start from the origin of the plot axis.
By examining and comparing the coefficients of efficiency for the performance of the product for both winter wheat and barley, including a common fuel utilization factor, among the scenarios of irrigation, scenarios with more irrigation were found to have a higher yield (Figure 6 and Figure 7).

3.2. Model Validation Results

Tests to validate the model were performed using data from the second year of cultivation (2019), and calibration was performed using planting data from the first year (2018). The results show that the minimum and maximum error estimates for the winter wheat grain yield amount to 1.05% when irrigation is 80% and 39.98% when irrigation is 50%. Furthermore, the minimum and maximum error estimates for the amount of biomass in the atmosphere are 0.95% when irrigation is 80% and 38.46% when irrigation is 50% (Table 8 and Table 9).
In addition, the efficiency of the model in the validation experiment was obtained for simulating grain yield, where R2, RMSE, d, EF, and NRMSE are equal to 0.99, 4, 0.97, 0.93, and 4.4 for winter wheat, respectively, and 0.97, 7, 0.94, 0.91, and 9 for barley.
Using the data in Figure 8 and Figure 9 for all scenarios of irrigation, the F-test was performed in order to compare the variance between the field data and plant models for both winter wheat and barley products at a 95% significance level (Table 10).
Figure 10 and Figure 11 display the comparisons of the observed and simulated performance (ha) for winter wheat and barley crops during the years 2018 and 2019 at different levels of water stress.

4. Discussion

Field measurements over a large area with arid and semiarid climatic conditions were used in this study. To validate the model, highly transparent data were used on twenty observed fields. On the one hand, a large area of Iran has an arid or semiarid climate. On the other hand, the water that agriculture needs is a limiting factor. Therefore, it seems that the AquaCrop model used in this study is suitable for studying the production function of two strategic products that are cultivated at a very high level in Iran, along with an approach to increasing water use efficiency per yield. In this study, the AquaCrop models accurately simulated crop canopy cover (CC), biological yield (BY), and grain yield (GY) under both regular and deficit irrigation. In water-limited environments, the AquaCrop model is a potentially valuable tool to use for maximizing winter wheat and barley yields in the region.
This study’s results illustrate that the current scenario-based policies for agriculture with 100% irrigation increase crop yield in existing conditions. Furthermore, with less irrigation, not only is productivity not reduced, but the percentage deficit in the amount of irrigation to save water instead increases productivity and the water efficiency factor by a notable degree.
AquaCrop models enable the simulation of transpiration during the growing season, yield, crop water requirements, and water use efficiency based on required water, climate, soil, and plants. Due to the high cost of operations and field measurements, as well as time and place limitations, in this study, the AquaCrop model was calibrated by using measurement data taken on the farm. Finally, after calibration, the results could be used to estimate performance in the future; therefore, although the AquaCrop model may have some problems, it seems that it can be used after being calibrated with a series of initial data in the field, which can be easily measured. The AquaCrop model also makes it easy to compute scheduling irrigation requirements, water consumption, and irrigation intervals, and compared to other plant models, the AquaCrop model requires less input data. Based on existing statistics, it is capable of accurately evaluating and simulating plant growth. The calibration of model parameters related to water stress, the crop canopy, and the reference harvest index at default values is required, as changes in these variables can be caused by the use or climate of the study area.

5. Conclusions

In reviewing the results of the calibration and validation of field data using the model, it is concluded that the percentage error in the estimations of grain yield, biomass, and water use efficiency in deficit irrigation scenarios increases with the decreasing amount of water required by the plant. The results of calibration and validation errors during the crop years 2018 and 2019 for yield, biomass, and water use efficiency were close in magnitude.
To determine an optimal irrigation level, data measured in different irrigation scenarios (full irrigation or 100% water requirement and 90, 80, 70, 60, and 50% irrigation) were calibrated and validated for arid and semiarid climatic conditions using the AquaCrop model. The model was calibrated to simulate the grain yields of winter wheat and barley using R2, RMSE, d, EF, and NRMSE statistical indicators. The obtained values of these indices were, respectively, 0.97, 3, 0.98, 0.94, and 4 for winter wheat and 0.98, 4, 0.92, 0.89, and 7 for barley. The model efficiency was also validated using crop harvest data in the crop year 2019. For the grain yield simulation, the calculated values of R2, RMSE, d, EF, and NRMSE statistical indicators were, respectively, 0.99, 4, 0.97, 0.93, and 4.4 for winter wheat and 0.97, 7, 0.94, 0.91, and 9 for barley. The data of field and modeled samples were analyzed by analysis of variance (ANOVA) using the F-test, and significant results were obtained for both crops in all applied scenarios at the 95% level. The AquaCrop model can be used to predict and optimize winter wheat and barley yields and water productivity under different irrigation treatments in Illam.

Author Contributions

Conceptualization, A.M.K.; methodology, M.N.; software, A.M.K.; validation, A.M.K.; formal analysis, A.M.K. and M.N.; investigation, A.M.K., M.N. and R.J.; resources A.M.K.; Project administration S.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

In this paper all data were gathered in the field by authors and they were presented in the main body of manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Abedinpour, M. Agricultural Water Management with AquaCrop Model, 1st ed.; Isfahan University Jahad Publications: Isfahan, Iran, 2019; pp. 11–18, 81–86, 210–219. [Google Scholar]
  2. Saab, M.A.; Albrizio, R.; Nangia, V.; Karam, F.; Rouphael, Y. Developing scenarios to assess sunflower and soybean yield under different sowing dates and water regimes in the Bekaa valley (Lebanon): Simulations with Aquacrop. Int. J. Plant Prod. 2014, 8, 1735–6814. [Google Scholar] [CrossRef]
  3. Hamidreza, S.; Mohd, A.M.S.; Teang, S.L.; Sayed, F.M.; Arman, G.; Mohd, K.Y. Application of AquaCrop model in deficit irrigation management of Winter wheat in arid region. Afr. J. Agric. Res. 2011, 6, 2204–2215. [Google Scholar]
  4. Shirshahi, F.; Babazadeh, H.; Ebrahimipak, N.; Zeraatkish, Y. Calibration and evaluation of the performance of Aquacrop model in managing the amount and time of deficit Irrigation in wheat. Iran. J. Irrig. Sci. Eng. 2017, 41, 31–44. (In Persian) [Google Scholar]
  5. Sema, K.A.L.E.; Madenoğlu, S. Evaluating AquaCrop Model for Winter Wheat under Various Irrigation Conditions in Turkey. J. Agric. Sci. 2017, 24, 205–217. [Google Scholar]
  6. Ahmadi, S.H.; Mosallaeepour, E.; Kamgar-Haghighi, A.A.; Sepaskhah, A.R. Modeling Maize Yield and Soil Water Content with AquaCrop Under Full and Deficit Irrigation Managements. Water Resour. Manag. 2015, 29, 2837–2853. [Google Scholar] [CrossRef]
  7. Ramezani, M.; Babazadeh, H.; Sarai Tabrizi, M. Simulating Barley Yield under Different Irrigation Levels by using AquaCrop Model. Irrig. Sci. Eng. (JISE) 2019, 41, 161–172. [Google Scholar] [CrossRef]
  8. Goosheh, M.; Pazira, E.; Gholami, A.; Andarzian, B.; Panahpour, E. Improving Irrigation Scheduling of Wheat to Increase Water Productivity in Shallow Groundwater Conditions Using Aquacrop. Irrig. Drain. 2018, 67, 738–754. [Google Scholar] [CrossRef]
  9. Hellal, F.; Mansour, H.; Abdel-Hady, M.; El-Sayed, S.; Abdelly, C. Assessment water productivity of barley varieties under water stress by AquaCrop model. AIMS Agric. Food 2019, 4, 501–517. [Google Scholar] [CrossRef]
  10. Khoshravesh, M.; Mostafazadeh-Fard, B.; Heidarpour, M.; Kiani, A.-R. AquaCrop model simulation under different irrigation water and nitrogen strategies. Water Sci. Technol. 2013, 67, 232–238. [Google Scholar] [CrossRef] [PubMed]
  11. Masasi, B.; Taghvaeian, S.; Gowda, P.H.; Warren, J.; Marek, G. Simulating Soil Water Content, Evapotranspiration, and Yield of Variably Irrigated Grain Sorghum Using AquaCrop. JAWRA J. Am. Water Resour. Assoc. 2018, 55, 976–993. [Google Scholar] [CrossRef]
  12. Hoseinzadeh, S.; Astiaso Garcia, D. Techno-economic assessment of hybrid energy flexibility systems for islands’ decarbonization: A case study in Italy. Sustain. Energy Technol. Assess. 2022, 51. [Google Scholar] [CrossRef]
  13. Hoseinzadeh, S.; Ghasemi, M.H.; Heyns, S. Application of hybrid systems in solution of low power generation at hot seasons for micro hydro systems. Renew. Energy. 2020, 160, 323–332. [Google Scholar] [CrossRef]
  14. Mousavizadeh, S.; Honar, T.; Ahmadi, S. Assessment of the AquaCrop Model for simulating Canola under different irrigation managements in a semiarid area. Int. J. Plant Prod. 2016, 10, 425–445. [Google Scholar] [CrossRef]
  15. He, Q.; Li, S.; Hu, D.; Wang, Y.; Cong, X. Performance assessment of the AquaCrop model for film-mulched maize with full drip irrigation in Northwest China. Irrig. Sci. 2020, 39, 277–292. [Google Scholar] [CrossRef]
  16. Amiri, A.; Bahrani, S.; Irmak, N.; Mohammadiyan, R. Evaluation of irrigation scheduling and yield response for wheat cultivars using the AquaCrop model in an arid climate. Water Supply. 2021, 22, 602–614. [Google Scholar] [CrossRef]
  17. Kheir, A.M.S.; Alkharabsheh, H.M.; Seleiman, M.F.; Al-Saif, A.M.; Ammar, K.A.; Attia, A.; Zoghdan, M.G.; Shabana, M.M.A.; Aboelsoud, H.; Schillaci, C. Calibration and Validation of AQUACROP and APSIM Models to Optimize Wheat Yield and Water Saving in Arid Regions. Land 2021, 10, 1375. [Google Scholar] [CrossRef]
  18. Kheir, A.M.; Hoogenboom, G.; Ammar, K.A.; Ahmed, M.; Feike, T.; Elnashar, A.; Liu, B.; Ding, Z.; Asseng, S. Minimizing trade-offs between wheat yield and resource-use efficiency in the Nile Delta–A multi-model analysis. Field Crop. Res. 2022, 287. [Google Scholar] [CrossRef]
  19. Shen, X.; Wang, G.; Zeleke, K.T.; Si, Z.; Chen, J.; Gao, Y. Crop Water Production Functions for Winter Wheat with Drip Fertigation in the North China Plain. Agronomy 2020, 10, 876. [Google Scholar] [CrossRef]
  20. Soomro, K.B.; Alaghmand, S.; Shahid, M.R.; Andriyas, S.; Talei, A. Evaluation of Aquacrop Model in Simulating Bitter Ground Water Productivity under Saline Irrigation. Irrig. Drain. J. 2019, 69, 63–73. [Google Scholar] [CrossRef]
  21. Abedinpour, M.; Sarangi, A.; Rajput, T.; 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]
  22. Najarchi, M.; Kaveh, F.; Babazadeh, H.; Manshouri, M. Determination of the yield response factor for field crop deficit irrigation. Afr. J. Agric. Res. 2011, 6, 3700–3705. [Google Scholar] [CrossRef]
Figure 1. Location map of the study area in Iran—Ilam Province.
Figure 1. Location map of the study area in Iran—Ilam Province.
Water 14 03455 g001
Figure 2. The results of model calibration for simulating grain yield, biomass, and water use efficiency for the winter wheat crop under different irrigation conditions.
Figure 2. The results of model calibration for simulating grain yield, biomass, and water use efficiency for the winter wheat crop under different irrigation conditions.
Water 14 03455 g002aWater 14 03455 g002b
Figure 3. The results of model calibration for simulating grain yield, biomass, and water use efficiency for the barley crop under different irrigation conditions.
Figure 3. The results of model calibration for simulating grain yield, biomass, and water use efficiency for the barley crop under different irrigation conditions.
Water 14 03455 g003aWater 14 03455 g003b
Figure 4. The crop production yield versus the irrigation level for the winter wheat crop.
Figure 4. The crop production yield versus the irrigation level for the winter wheat crop.
Water 14 03455 g004
Figure 5. The crop production yield versus the irrigation level for the barley crop.
Figure 5. The crop production yield versus the irrigation level for the barley crop.
Water 14 03455 g005
Figure 6. The productivity coefficient versus crop yield for the winter wheat crop.
Figure 6. The productivity coefficient versus crop yield for the winter wheat crop.
Water 14 03455 g006
Figure 7. The productivity coefficient versus crop yield for the barley crop.
Figure 7. The productivity coefficient versus crop yield for the barley crop.
Water 14 03455 g007
Figure 8. The results of model calibration for simulating the crop yield of the winter wheat crop under different irrigation conditions.
Figure 8. The results of model calibration for simulating the crop yield of the winter wheat crop under different irrigation conditions.
Water 14 03455 g008
Figure 9. The results of model calibration for simulating the crop yield of the barley crop under different irrigation conditions.
Figure 9. The results of model calibration for simulating the crop yield of the barley crop under different irrigation conditions.
Water 14 03455 g009
Figure 10. Comparison of the yield (t/ha) simulated by the AquaCrop model with values observed for the winter wheat crop at different water stress levels.
Figure 10. Comparison of the yield (t/ha) simulated by the AquaCrop model with values observed for the winter wheat crop at different water stress levels.
Water 14 03455 g010
Figure 11. Comparison of the yield (t/ha) simulated by the AquaCrop model with values observed for the barley crop at different water stress levels.
Figure 11. Comparison of the yield (t/ha) simulated by the AquaCrop model with values observed for the barley crop at different water stress levels.
Water 14 03455 g011
Table 1. Growth period duration and the planting/growing dates of the crop.
Table 1. Growth period duration and the planting/growing dates of the crop.
CropGrowth Period (Days) in Different StagesPlanting DateDate of Harvest
Vegetative GrowthFloweringSeedingGranulationSenescence
Winter wheat30607020180The second half of NovemberThe second half of May
Barley20607020170The second half of NovemberThe third half of April
Table 2. Estimation of the water requirements of the planting patterns.
Table 2. Estimation of the water requirements of the planting patterns.
CropMaximum Depth of Root Development (cm)Net Volume of Irrigation Water (m3)Gross Volume of Irrigation Water (m3)Net Need for Irrigation Cycle (mm)Gross Need for Irrigation Cycle (mm)Maximum Hydromodule (lit·s−1) Number of Irrigations
Winter wheat12037896888631151.156
Barley10030925622621131.045
Table 3. Climatic characteristics over a 20-year statistical period (1999–2019) in Dehloran.
Table 3. Climatic characteristics over a 20-year statistical period (1999–2019) in Dehloran.
MonthET
(mm·day−1)
Sunny HoursWind Speed (Km·day−1)Humidity (%)Maximum Temperature (°C)Minimum Temperature (°C)Effective Rainfall (mm)Rainfall (mm)
January1.936.31645918.18.833.541.9
February3.017.21814822.712.124.230.3
March4.667.72254228.416.628.335.4
April5.988.21813235.422.918.322.9
May10.0110.73022242.729.10.10.1
June10.8311.32942045.631.700
July10.2810.92592046.63200
August10.2810.72512143.428.71.11.4
September7.059.22162637.4242.53.1
October4.217.11814328.116.922.728.4
November2.36.31475820.410.839.549.4
December1.696.21476317.48.339.549.4
Table 4. Soil input data for the model.
Table 4. Soil input data for the model.
SoilValue
Soil typeLumi Sandy
Saturated hydraulic conductivity (mm·day−1)1200
Saturated moisture (V %)41
Crop capacity point (V %)22
Permanent wilting point (V %)10
Thickness of soil layer (m)2.5
Soil penetration coefficient46
Bulk density (gr·cm−3)1.4
Table 5. Plant parameters used in the model for wheat and barley simulation.
Table 5. Plant parameters used in the model for wheat and barley simulation.
Crop CharacteristicsCropUnit
Winter WheatBarley
Crop typeRootRoot
Planting methodSowingSowing
Category of the plant in terms of carbonC3C3
Cropping period22 October22 October
Length of growing cycle180170Days
Canopy developmentCanopy growth coefficient (CGC)18.4011% inc. in CC relative to existi. CC per GDD
Canopy decline coeff. (CDC) at senescence8.709.40%; decrease in CC relative to CCx per GDD
Canopy cover (CCo)66% at 90% emergence
Maximum canopy cover9393CCx (%)
Shading surface during germination1.501.50cm2
Growing cycleGermination3020day
Flowering6060day
Granulation period9070day
Senescence11090day
Root deepeningMin.0.200.20m
Max.1.201.00m
Time to reach maximum root depth7070Day
TemperatureBase temperature00°C
Cut-off temperature3015°C
Minimum degree of pollination55°C
Maximum degree of pollination3535°C
Harvest Index50%
Soil water drainage deduction for vegetation developmentP(upper)0.250.25At this amount, vegetative growth stops
P(lower)0.600.60At this amount, vegetative growth stops
Upper threshold of stomatal conductanceP(upper)0.650.65Above this, stomata begin to close
Upper threshold of senescence stressP(upper)0.650.65
Canopy growth factor53
Stomatal control method factor2.503
TranspirationTranspiration coefficient at maximum coverage1.151
Effect of canopy on reducing evaporation at the end of growth5050%
Percentage decrease in Kc with age0.150.13%
Irrigation methodFurrow irrigation
Available water7070
Number of irrigations65
Salinity stressUpper limit threshold of salinity stress1515ds·m−1
Salinity threshold decreases yield67ds·m−1
Table 6. The results of model calibration for grain yield, biomass, and water use efficiency for the winter wheat crop under different irrigation conditions (2018).
Table 6. The results of model calibration for grain yield, biomass, and water use efficiency for the winter wheat crop under different irrigation conditions (2018).
TreatmentYield (ton/ha)PeBiomass (ton/ha)PeWP (kg·m−3)Pe
Obs.Sim.(±%)Obs.Sim.(±%)Obs.Sim.(±%)
T1004.504.694.228.8378.6681.910.660.6823.33
T904.0374.214.298.3578.1981.900.710.743.04
T803.6773.793.077.4827.7273.270.790.8148.86
T703.2003.292.816.9517.1813.300.850.872.35
T602.9183.013.156.3156.5032.970.920.9513.37
T502.5832.652.595.5845.8594.921.131.084.42
Note: Pe: Simulation error.
Table 7. The results of model calibration for grain yield, biomass, and water use efficiency for the barley crop under different irrigation conditions (2018).
Table 7. The results of model calibration for grain yield, biomass, and water use efficiency for the barley crop under different irrigation conditions (2018).
TreatmentYield (ton/ha)PeBiomass (ton/ha)PeWP (kg·m−3)Pe
Obs.Sim.(±%)Obs.Sim.(±%)Obs.Sim.(±%)
T1002.802.924.2911.7111.343.160.710.744.23
T902.5482.643.6110.8410.5992.220.760.793.95
T802.3352.413.219.7219.8571.400.840.822.38
T702.0022.14.908.8639.1122.810.870.903.45
T601.791.874.477.9478.1832.971.021.063.92
T501.551.63.236.8827.1954.551.061.103.77
Note: Pe: Simulation error.
Table 8. The results of model validation for grain yield, biomass, and water use efficiency for the winter wheat crop under different irrigation conditions (2019).
Table 8. The results of model validation for grain yield, biomass, and water use efficiency for the winter wheat crop under different irrigation conditions (2019).
TreatmentYield (ton/ha)PeBiomass (ton/ha)PeWP (kg·m−3)Pe
Obs.Sim.(±%)Obs.Sim.(±%)Obs.Sim.(±%)
T1004.84.476.879.148.971.860.730.752.74
T904.13.982.938.878.701.910.760.793.95
T803.83.761.057.937.673.280.840.873.57
T703.343.097.487.487.243.210.910.954.40
T603.132.878.37.076.872.830.981.035.10
T502.732.518.066.626.324.531.181.223.39
Table 9. The results of model validation for grain yield, biomass, and water use efficiency for barley crop under different irrigation conditions (2019).
Table 9. The results of model validation for grain yield, biomass, and water use efficiency for barley crop under different irrigation conditions (2019).
TreatmentYield (ton/ha)PeBiomass (ton/ha)PeWP (kg·m−3)Pe
Obs.Sim.(±%)Obs.Sim.(±%)Obs.Sim.(±%)
T1003.13.254.8412.312.002.440.740.7785.14
T902.812.923.9111.4611.281.570.790.823.8
T802.482.563.2310.4710.370.950.870.914.6
T702.422.503.319.869.652.130.910.9625.71
T602.112.193.798.778.572.281.11.1524.73
T501.831.882.737.627.363.411.171.224.27
Table 10. F-test results of field data and model for winter wheat and barley crops.
Table 10. F-test results of field data and model for winter wheat and barley crops.
Absolute Value of the Critical Point of Winter Wheat: 3.43
Absolute Value of the Critical Point of Barley: 3.17
Hypothesis
H0: Yij = Yji
H1: Yij ≠ Yji
Significance Level: 0.05
Irrigation scenariosAbsolute value of F statisticAbsolute value of F statistic
Winter WheatBarley
50%1.151.06
60%1.181.31
70%1.101.28
80%1.131.39
90%1.151.30
100%1.061.15
Note: Source: Research Findings.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Khoshsirat, A.M.; Najarchi, M.; Jafarinia, R.; Mokhtari, S. Sensitivity Analysis and Determination of the Optimal Level of Water Use Efficiency for Winter Wheat and Barley under Different Irrigation Scenarios Using the AquaCrop Model in Arid and Semiarid Climatic Conditions (Case Study: Dehloran Plain, Iran). Water 2022, 14, 3455. https://doi.org/10.3390/w14213455

AMA Style

Khoshsirat AM, Najarchi M, Jafarinia R, Mokhtari S. Sensitivity Analysis and Determination of the Optimal Level of Water Use Efficiency for Winter Wheat and Barley under Different Irrigation Scenarios Using the AquaCrop Model in Arid and Semiarid Climatic Conditions (Case Study: Dehloran Plain, Iran). Water. 2022; 14(21):3455. https://doi.org/10.3390/w14213455

Chicago/Turabian Style

Khoshsirat, Amir Mahyar, Mohsen Najarchi, Reza Jafarinia, and Shahroo Mokhtari. 2022. "Sensitivity Analysis and Determination of the Optimal Level of Water Use Efficiency for Winter Wheat and Barley under Different Irrigation Scenarios Using the AquaCrop Model in Arid and Semiarid Climatic Conditions (Case Study: Dehloran Plain, Iran)" Water 14, no. 21: 3455. https://doi.org/10.3390/w14213455

APA Style

Khoshsirat, A. M., Najarchi, M., Jafarinia, R., & Mokhtari, S. (2022). Sensitivity Analysis and Determination of the Optimal Level of Water Use Efficiency for Winter Wheat and Barley under Different Irrigation Scenarios Using the AquaCrop Model in Arid and Semiarid Climatic Conditions (Case Study: Dehloran Plain, Iran). Water, 14(21), 3455. https://doi.org/10.3390/w14213455

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