Abstract
Agriculture is vulnerable to drought indicating that the increasing climate crisis requires the necessity of sustainable crop production. In this study, we developed the Irrigation Schedule and Management (ISM) model based on a simulation–optimization (Soil Water Atmosphere Plant-SWAP model with Genetic Algorithm-GA) framework. The ISM model finds an optimal combination of Irrigation Water Amount (IWA) and Irrigation Interval (II) by adjusting Water Stress (WS) responding to environmental conditions (weather, soils, crops and bottom boundary conditions) throughout growing periods. By conditioning the crop (WS) and water management (IWA and II) variables, ISM improves the sustainability of optimal crop productions under different climatic-land surface conditions. The Regional Agromet Center (RAC) site in Faisalabad (at Punjab, Pakistan) was selected to test the proposed ISM model for the field validation/synthetic numerical experiments with various crops (Wheat, Corn and Potato) and soils. We demonstrated that the ISM model that reflects the relationship between crop and water management variables improved the sustainability of crop productions and Water Productivity (WP) compared to those of the conventional irrigation method at the RAC site under various environment conditions. Additionally, the ISM-based long-term crop productions showed the variations along the yearly precipitation changes indicating that optimal combinations of the crop and water management variables are considerably influenced by environmental conditions. Although uncertainties exist, our proposed ISM model can contribute to the establishment of efficient irrigation schedule/management plans under agricultural drought.
Similar content being viewed by others
Introduction
The frequency and intensity of natural disasters (e.g., floods, droughts, etc.) have been gradually increasing across the world as a result of climate changes (Lettenmaier et al. 1996; Aswathanarayana 2001; Gobiet et al. 2014; Dai et al. 2018). Mishra and Singh (2009) reported that global/local regions might be exposed to more severe droughts in the future (2001–2050). Drought causes widespread and severe economic losses in hydrology, agriculture, and ecosystem across various regions (FEMA 1995; Svoboda et al. 2002). Agriculture is highly vulnerable to drought due to the meteorological conditions. A large amount of available fresh water across the world (about 70% of the existing freshwater supplies) is utilized for agriculture (Evans and Sadler 2008). Therefore, an increase in water scarcity might be inevitable in the future indicating that the agricultural sector needs to be prepared to deal with drought.
Soil moisture is a key factor in crop production because plant roots directly absorb soil moisture for crop growth (Morita and Abe 1996). Deficient soil moisture contents disturb normal root growth (e.g., reduction of the number of roots, total root length, limiting of crop growth and yield) while excessive soil moisture conditions cause notable stimulation (Pardales et al. 2002; Sugimoto et al. 2000; Lv et al. 2011). A few studies reported that frequent irrigations with low water amount can improve crop productions (Jin et al. 2014). Appropriate irrigation supplies to wheat crop during the critical growth stage contribute to the improvement of photosynthetic rate, dry matter production, and transportation indicating that crop productions can be significantly increased compared to those under the rain-fed condition (Boughdiri et al. 2014; M’hamed et al. 2015; Saeedipour and Moradi 2011). Dilsher and Khalaf (2014) reported that crop productions were increased by 37.5% and 84.5% with irrigation applied done once/twice were increased compared with that of non-irrigation. Especially, timely irrigation significantly contributes to the improvement of the water use efficiency (WUE) of crops (Sun et al. 2006). Several studies on water-saving concepts (irrigation methods for crops) have been conducted for the efficient use of water in agriculture (Seckler 1996; Molden 1997; Ines et al. 2002). These water-saving concepts contributed to the improvement of water productivity in irrigation (Bos and Nugteren 1990; Bos 1997), but the efficiency indicators only adjusted the water flow ratio for the reuse of water in the irrigation system.
Crop growth highly depends on both climatic-land surface conditions (solar radiation, temperature, rainfall, soil physical characteristics, etc.) and artificial crop/water (Water Stress-WS, Irrigation Water Amount-IWA and Irrigation Interval-II) management (Olesen and Bindi 2002). Considering that artificial crop/water management are only adjustable for outdoor cultivation, irrigation is the key factor that can control crop growth, development and yield (Gornall et al. 2002; Prasad et al. 2008). Thus, the requirement of timely efficient irrigation schedule and management is inevitable for the sustainability of crop production against rapidly deteriorating drought in future. Water and crop management during crop growing seasons can improve the performance of the irrigation system (Ines et al. 2006), while Shin and Jung (2014) suggested the drought index-based irrigation water management approach that can adjust irrigation water based on drought condition. Basically, crop status (WS) under limited water use is highly related to irrigation water (IWA and II) management (Ines et al. 2006; Lee and Shin 2016; Domínguez et al. 2017). In other word, the degree of WS affecting by regional climatic-land surface characteristics significantly influences on determining IWA and II in crop cultivation and water efficiency under irrigation water deficit. Various research works have been done for optimal irrigation decisions, but few studies with respect to timely irrigation schedule and management that can consider the optimal combination (considering both crop and water management) of IWA and II by adjusting WS responding to environmental factors (weather, soil and crop) throughout crop growing seasons have been conducted for the sustainability of crop production. Considering that food security is becoming a global issue due to climate crisis, rapidly increasing world population, continental-/local-disputes, etc., an optimal combination of crop and water management is highly required for the sustainability of optimal (maximum) crop production and water efficiency.
The main goal of this study is to improve the sustainability of optimal crop production by considering timely efficient irrigation management under agricultural drought. The objectives of this research were as follows: (1) to develop an Irrigation Schedule and Management (ISM) model that can determine an optimal combination of IWA and II by adjusting WS responding to environmental factors (weather, soil and crop) based on a simulation–optimization framework and (2) to provide an efficient irrigation management plan for the sustainability of optimal crop production under different climatic-land surface conditions.
Materials and methods
Concept of the ISM model
In this study, we developed the ISM model for improving the sustainability of optimal crop production by considering timely efficient irrigation schedule/management plans under irrigation water deficit. The ISM model adapting the simulation (Soil Water Atmosphere Plant—SWAP, Kroes et al. 1999; van Dam et al. 1997)–optimization (Genetic Algorithm—GA, Goldenberg 1989; Holland 1975) framework is capable of sustaining optimal crop productions at field regions by determining the optimal combination of IWA and II by adjusting WS responding to environment conditions (weather, soil and crop) in Fig. 1. The ISM model incorporated with GA randomly generates the WS values and determines the IWA and II values based on the generated WS. The physically-based hydrological model of SWAP incorporated within ISM estimates the daily root zone (0–30 cm) soil moisture dynamics. Here, SWAP requires the input data of weather data (solar radiation, max. and min. temperature, wind speed, precipitation, humidity), soil hydraulic properties (α, n, θres, θsat, Ksat), and crop information for simulation. Also, WOFOST (WOrld FOod STudies) (van Dam et al. 1997; van Dam 2000) adapted within the SWAP model can estimate crop production (kg ha−1) for various crops at field-scales under different hydro-climatic conditions. To obtain the soil hydraulic properties, the ISM model was combined with the near-surface soil moisture data assimilation scheme (Ines and Mohanty 2008) based on Inverse Modeling (IM). The data assimilation scheme extracts the effective soil hydraulic properties that can represent Remotely Sensed (RS)/measured (in situ) soil moisture values by minimizing the differences between the measured and simulated soil moisture contents. The SWAP model also calculates the actual/potential evaporation (Ea/Ep) and transpiration (Ta/Tp) based on the Penman–Monteith equation. Basically, the ISM model determines the IWA and II variables based on the WS ([1−Ta/Tp] × 100) levels. Considering that the WS levels are allowable between 5 and 40% in fields before irrigation (Berkoff and Huppert 1987; Bandaragoda 1998; Polinova et al. 2019; Bonfante et al. 2019; Tolomio and Casa 2020), the ISM model randomly generates the WS values in the range of 10 (wet) to 40% (dry) based on GA. Then, GA derives the IWA and II values based on the generated WS until GA finds the best combination (solution, k = {IWAi=1,…,I, II, WS}) of IWA, II and WS using the objective function (Z) adapting the maximization function (N∈M) in Eq. (1) along with the given GA generations (1 to G). Mathematically, the best combination is determined by obtaining the better Crop Production (CP) in the current generation (G = g) compared to those of the previous generations (G = 1,…,g−1) based on the objective function. Here, the Irrigation Schedule Date (ISD) value was initially determined based on the WS (e.g., [1−Ta/Tp] × 100 > 40%) value across the daily time-scale. When the dry days are continued since the ISD is initially determined with the water scarcity (e.g., WS > 40%) for the individual irrigation events, ISM determines II within the constrained ranges (14–60 days used in this study) to trigger irrigation (flood irrigation method used in this study) for removing drought damages to crops. GA usually finds a single solution from the search space indicating that GAs have limitations that can be trapped at the local regions (minima). For this reason, we integrated an ensemble scheme with GA, called the Ensemble-based Genetic Algorithm (EGA), to overcome GA limitations for better searching solutions and determining uncertainty ranges of searched solutions.
Here, Z is the objective function, e is the running index of ensembles (E = 10), E is the total number of ensembles, CP is the crop production, g is the running index of EGA, G is the total number of generations in the GA process, N is the individual population, and M is the total number of populations in EGA.
Soil–water–atmosphere–plant (SWAP)
The ISM model adapts the one-dimensional (1-D) SWAP model to simulate soil water flow through the soil profile by using the Richards’ equation (Eq. 2) with the input variables of weather, soil, and crop information. The soil hydraulic functions can be described using the analytical expression (van Genuchten 1980) based on the relationship of the soil water content (θ), pressure head (h), and unsaturated hydraulic conductivity (K):
where θ is the water content in soil (cm3 cm−3), K is the hydraulic conductivity (cm d−1), h is the pressure head (-cm), z is the soil depth (cm), t is the time domain (d), C is the differential water capacity (cm−1), and S(h) is the actual soil moisture rate extracted by the crop root (cm3 cm−3 d−1), as shown in Eq. (3):
In the equation, Tpot is the potential transpiration (cm d−1), Zr is the rooting depth (cm), and αw is the reduction factor by the pressure head (h), and it accounts for water deficit/oxygen stress (Feddes et al. 1978). The finite difference approach incorporated in Richards’ Eq. (2) (Belmans et al. 1983) can simulate the soil water flow using the soil hydraulic properties and various management scenarios. The soil hydraulic functions of θ(h) and K(h) are derived based on the analytical expressions of van Genuchten (1980) and Mualem (1976) as below,
where Se is the relative saturation (−); θres is the water content (cm3 cm−3) of soil in residual, and θsat is the water content (cm3 cm−3) of soil in saturation, α (cm−1), n(−), m(−), λ(−) are the shape parameters for the soil hydraulic functions, Ksat is the hydraulic conductivity (cm d−1) of soil in saturation, and m = 1 − \(\frac{{1}}{{\text{n}}}\).
Various top (weather variables) and bottom boundary (free-drainage and shallow ground water table depths) conditions need to be set manually to simulate the soil moisture dynamics (van Dam et al. 1997). SWAP is comprised of three crop modules: (1) a simple module to consider weather conditions, soil texture, and crop type, (2) a detailed module (WOFOST), and (3) the WOFOST module for simulating grass growth. Moreover, SWAP combines the water management modules to consider irrigation and drainage processes for irrigation events and schedules (van Dam et al. 1997; van Dam 2000). The SWAP model calculates potential evapotranspiration (ETp) based on the Penman–Monteith equation. Soil evaporation (Ep) and potential transpiration (Tp) are calculated based on the leaf area and the soil cover fraction. ETp is converted into actual ET (ETa) by adjusting soil evaporation and potential transpiration using empirical relationships as the soil becomes dry. The SWAP model performs well with various meteorological and environmental criteria; detailed information is included in the references (Wesseling and Kroes 1998; Sarwar et al. 2000; Droogers et al. 2000; Singh et al. 2006).
Ensemble-based genetic algorithm (EGA)
Genetic algorithms (GAs) are one of the various optimization algorithms that solve complex problems (Holland 1975; Goldberg 1989). The performance of GAs is highly affected by initial random generator seeds (e.g., idum: − 3000, − 2000, − 1000, etc.). Also, GAs usually have limitations that can be trapped at local regions (minima) at the unknown search space. To minimize the above-mentioned drawbacks of GAs, the Ensemble-based Genetic Algorithm (EGA) was developed for better-searching solutions and assessing uncertainties of searched solutions (Shin and Mohanty 2013). EGA derives the fittest chromosomes (k = {IWAi=1,…,I, II, WS}) in the individual population via the total generation (G). The derived chromosomes have new genetic values via the GA operators, which comprise the selection, crossover and mutation processes and search for more solution spaces. Moreover, EGA restarts from other search spaces to find solutions when the reproduced chromosomes converge to one region before all generations are completed. Using the restarting function, EGA can derive new chromosomes via the creep/jump mutation processes (Ines and Honda 2005). Furthermore, EGA remembers the best (elite) chromosomes in the previous generation (g − 1) and delivers them to the next generation (g) (Ines and Mohanty 2008). To search more unknown search spaces, a random resampling (ensemble) algorithm (IBM Programmers’ Guide) (Efron 1982) was integrated with EGA. The constrained ranges of solutions by GA are shown in Table 1.
Numerical Experiments
The field validation and synthetic numerical experiments were conducted to test the performance of ISM model for improving the sustainability of crop productions. The numerical experiments were as follows: (1) Case Study 1: field validation experiments and (2) Case Study 2: synthetic numerical experiments, respectively. The Regional Agromet Center (RAC) site located in the district of Faisalabad at Punjab (Pakistan) under the semi-arid/arid climate condition was selected for conducting the Case Studies (Fig. 2). The wheat crop was cultivated during the winter periods (Nov. to Apr.) of 2007–2008 and 2008–2009 with sandy loam soil (predominant). The soil hydraulic properties (α: 0.018, n: 1.260, θres: 0.044, θsat: 0.340, Ksat: 145.0) at 0–40 cm depth measured by Ahmad (2002) and Time Domain Reflectometry (TDR)-based in situ root zone (0–30 cm) soil moisture measurements (Hanif 2008) were used for the field validation of near-surface soil moisture data assimilation scheme that simulates the soil moisture dynamics. The in situ soil moisture data were measured at the multiple soil depths of 5, 10, 20, and 30 cm at the RAC site during the winter seasons (2007–2008 and 2008–2009). The averaged root zone soil moisture measurements in the multiple soil depths were used for the field validation experiments. The weather data (e.g., precipitation, maximum and minimum temperature, wind speed, humidity and solar radiation) were collected from the Pakistan Meteorological Department for the RAC site in Faisalabad (Pakistan). Also, the Field-scale Crop and Irrigation Management Information-FCIMI (Hanif 2008) such as crop growth period, sowing/harvesting dates, irrigation amounts (mm), irrigation type (flood), number of irrigations and measured crop production (kg ha−1) for wheat are shown in Table 2.
In Case Study 1, the near-surface soil moisture data assimilation scheme was calibrated for extracting the effective soil hydraulic properties (α, n, θres, θsat, Ksat) from the TDR-based in situ soil moisture measurements (Hanif 2008) during the winter season (2007–2008). The soil hydraulic functions of θ(h) and K(h) based on the extracted (calibrated) soil hydraulic parameters were compared to those of the measured soil hydraulic parameters (Ahmad 2002). Based on the calibrated soil hydraulic parameters (2007–2008), the root zone (0–30 cm) soil moisture dynamics were predicted and compared with the TDR-based measurements during 2008–2009 for validation of the near-surface soil moisture data assimilation. Then, the crop productions estimated by the SWAP (only)/ISM (integrated with SWAP) models were compared with the measured crop productions (Hanif 2008) during the winter seasons for calibration (2007–2008) and validation (2008–2009).
In Case Study 2, the synthetic numerical experiments were conducted to analyze the impacts of various soils/crops (corn, potato and wheat) on crop productions. Various reference-based soil hydraulic properties (Clay Loam-CL, Loam-L, Sandy Clay Loam-SCL, Silty Clay Loam-SiCL, Silt Loam-SiL, and Sandy Loam-SL) were used from the UNSODA databases (Leij et al. 1999). Figure 3 shows the soil hydraulic functions of θ(h) and K(h) derived from the UNSODA database. The Bottom Boundary Condition (BBC) of free-drainage and Shallow Ground Water (SGW) table depths have a high influence on crop growth and production at field-scales. Thus, free-drainage and different SGW table depths (− 200, − 150, and− 100 cm from the soil surface) were tested for the ISM model with various soils and crops at the RAC site. In addition, the long-term crop productions were predicted by the ISM model from 1989–1990 to 2008–2009 based on the calibrated soil hydraulic properties (2007–2008). We assumed that the soil profile is homogeneous at the soil depth of 0–200 cm from the soil surface (Fig. 4), and the soil profile was discretized into the 33 soil layers (1–10th layer: the intervals of 1 cm, 11–20th layer: the intervals of 5 cm, 21–32nd layer: the intervals of 10 cm, and 33rd layer: the interval of 20 cm). The RAC site has the ground water table depth of 800 ~ 1000 cm from the soil surface (Irrigation Department, http://irrigation.punjab.gov.pk/). For this reason, the BBC was set as the free-drainage condition under the assumption that the initial condition is in equilibrium with the presence of BBC under the modeling condition. Additionally, the Water Productivity (WP = CP/IWA, kg ha−1 mm−1, French and Schultz 1984; Bessembinder et al. 2005; Passioura 2006) values were estimated for evaluating irrigation water use efficiency. To determine uncertainties of the model outputs, we set the number of ensemble (E = 10) for individual populations in EGA (Table 2). Also, the Pearson’s correlation (R) and root-mean-square error (RMSE) were used to evaluate the model performance as shown in Eqs. (6) and (7),
where \(\theta_{{\text{t}}}^{{{\text{obs}}}}\): the observed soil moisture with time (t), \(\overline{\theta }_{{\text{t}}}^{{{\text{obs}}}}\): the average of observed soil moisture with time (t), \(\theta_{{\text{t}}}^{{{\text{sim}}}}\): the simulated soil moisture with time (t), and \(\overline{\theta }_{{\text{t}}}^{{{\text{sim}}}}\): the average of simulated soil moisture with time (t), T is the time domain and t is the running index of time, respectively.
Results and discussion
Case study 1: Field validation experiment
In this study, we developed the ISM model that can improve the sustainability of crop productions with the optimal combination of IWA and II determined by adjusting WS responding to the environmental conditions (weather, soils and crops). The effective soil hydraulic properties as the input variables to SWAP were derived using the near-surface soil moisture data assimilation scheme. Note that the actual irrigation events (individually 75 mm for 1st: Dec. 11th in 2007 and 2nd: Jan. 8th, and 3rd: Mar. 2nd, 4th: Mar. 25th in 2008) at the field site were included in the data assimilation process. Overall, the RAC site has the less precipitation amounts during the crop growing seasons (Nov. to Apr.) in the years of calibration (2007–2008) and validation (2008–2009) while the relatively higher precipitation amounts and frequencies were generated in May to June. Figure 5a shows the comparisons of measured and near-surface data assimilation (ensemble-averaged)-based soil moisture dynamics for calibration at the RAC site during the winter period (2007–2008). The calibrated soil moisture dynamics (R: 0.908 and RMSE: 0.014) matched well with the measurements. The good agreement between the soil hydraulic functions of θ(h) and K(h) based on the measured (Ahmad 2002)/calibrated soil hydraulic properties support the robustness of near-surface soil moisture data assimilation scheme in Figs. 5b, c. For validation (2008–2009), the soil moisture dynamics based on the calibrated soil hydraulic properties (2007–2008) were estimated using the SWAP model in a forward mode and compared with the TDR measurements in Fig. 6. Although uncertainties were slightly increased compared to the calibration results (Fig. 5a), the estimated root zone soil moisture dynamics (R: 0.721 and RMSE: 0.035) were highly comparable with the in situ measurements during the winter period of 2008–2009.
To test the model performance of SWAP (WOFOST module) for simulating crop productions, the wheat crop productions for the periods of 2007–2008 and 2008–2009 based on the effective soil hydraulic properties calibrated during the period of 2007–2008 were also estimated using SWAP in a forward mode and compared with the measurements at the RAC site in Figs. 7a, b. The SWAP-based crop production (3139 kg ha−1) using the actual IWA and II (as shown in Table 2) matched well with the measured data (3212 kg ha−1) during the period of 2007–2008 while the estimated crop production (3260 kg ha−1) for the period of 2008–2009 was slightly underestimated compared to the measured wheat production (3707 kg ha−1) with increased uncertainties. Although uncertainties slightly increased in the period of 2008–2009, these results support the good performance of SWAP for simulating crop productions at field-scales.
Before testing the performance of ISM model, the ISM-based II with the fixed WS ([1−Ta/Tp] × 100%) of 10 to 40% were estimated and compared to the field-scale ISD (1st: Dec. 11th 2007, 2nd: Jan. 08th 2008, 3rd: Mar. 2nd 2008 in Table 2) determined by farm managers manually (Fig. 8). With the multiple WS values from 10 to 40%, we confirmed that the ISM-derived irrigation required (dry) days were increased, respectively. In general, the conventional irrigation method empirically determines the IWA and II by considering the increasing WS due to agricultural drought. Compared to the conventional method, we confirmed that the ISM model determines the ISD (indicating irrigation events) right after the generations of WS throughout the crop growing period under the fixed WS (10, 20, 30 and 40%) conditions. These results supported that the ISM model can respond immediately to changes of the environmental conditions (weather, soils and crops) throughout the growing period.
In Table 3, the ISM-based wheat crop productions were estimated and compared with the measurements for the years of 2007–2008 and 2008–2009. The measured wheat crop productions of 3212/3707 kg ha−1 were shown in Table 2 with the precipitation amounts (108 mm for 2007–2008 and 96 mm for 2008–2009) during the crop growing periods, respectively. Overall, the ISM-based crop productions (4289 kg ha−1 for 2007–2008 and 6297 kg ha−1 for 2008–2009) were significantly improved, while the estimated IWA values slightly increase compared to the measured data (FCIMI in Table 2) by the conventional irrigation method (manually operated by farm mangers). It is inferred that the increased IWA might be due to the prompt responses of ISM model to the generated WS at the field site. The decreased II values (6.1 and 4.4 days for 2007–2008 and 2008–2009) compared to those (about monthly irrigation intervals) of the conventional method (Table 2) support this inference. Also, the WP values of ISM for the periods of 2007–2008 (12.2 kg ha−1 mm−1) and 2008–2009 (22.7 kg ha−1 mm−1) were higher than those (11.8 and 13.6 kg ha−1 mm−1) of the conventional method (Table 2), respectively. Especially, the ISM-based WP value for 2008–2009 was significantly increased indicating that the ISM model contributed to both the optimal crop production and irrigation water use efficiency. These results supported that the ability that can consider the relationship between crop (WS) and water management (IWA and II) variables allows the ISM model to respond immediately to environmental changes.
Case study 2: Synthetic numerical experiment
Figure 9 presents the rain-fed/ISM-based crop productions and IWA with the combinations of different soil types (CL, L, SCL, SiCL, SiL, and SL) and crops (wheat, corn and potato) at the RAC site during the crop growing period (2007–2008) under the synthetic conditions. The crop productions under the rain-fed condition were significantly affected by different soil types. Physical characteristics of the L (38.3 cm/day), SiL (30.5 cm/day) and SL (41.6 cm/day) soils in Fig. 3 weakens the crop productions compared to those of the other soils (1.8 ~ 9.78 cm/day for CL, SiCL, and SiL) indicating that crops with the predominant L, SiL and SL soils having the relatively high infiltration at field regions could be highly vulnerable to water deficit under agricultural drought. Overall, the ISM-based crop productions for the whole soils and crops were highly improved, especially for the L, SCL and SL soils. We also confirmed that the ISM-based IWA values were proportional to the increased crop productions, but the estimated crop productions using the ISM model were still fluctuated by the physical soil characteristics. These trends for wheat, corn and potato were similarly shown in Fig. 9a–c, but the differences ([rain-fed CP/ISM-based CP] × 100%) between the rain-fed and ISM-based CP values showed variations for the combination of soils and crops in Fig. 9. It is inferred that alternative crop selections by considering physical soil characteristics at local regions can improve to the sustainability of crop productions. Additionally, the SL soil having the highest infiltration rate in Fig. 3 was extremely vulnerable to water deficit condition, but the crop productions with the optimal IWA and II were dramatically increased with comparing to those of the other soils. The soil pore space for the SL soil is restively larger than the other soils indicating that air circulation and diffusion through large pore space actively occur in the soil matrix and are favorable to crop growth (Goorahoo et al. 2001; Du et al. 2018). The similar trends were shown for the L and SCL soils, respectively. Although the crop productions have variations based on the physical characteristics of soils, the ISM model sustained the relatively maximum levels of crop productions for the whole soils and crops.
In Fig. 10, we tested the impacts of BBCs such as free-drainage, SGW − 200 cm, SGW − 150 cm and SGW − 100 cm on estimating the wheat crop production at the RAC site during the crop growing period (2007 − 2008) under the synthetic condition, respectively. The crop production was only 4 kg ha−1 for free-drainage under the rain-fed condition, while SGW − 200 cm showed the crop production of 923 kg ha−1. When the SGW table depths of − 150 and − 100 cm become closer to the soil surface, the crop productions (3286/3467 kg ha−1) was considerably increased, respectively. These findings indicated that the presence of SGW that becomes closer to the soil surface in the soil profile can be more easily accessible to crop root based on the root moisture extraction force and capillary phenomena. When we applied the ISM model to the wheat crop with the different BBC, the whole crop productions reached to the maximum levels of 4289, 4439, 4279 and 3465 kg ha−1 for free-drainage, SGW − 200 cm, SGW − 150 cm and SGW − 100 cm with the IWA of 435, 272, 147 and 59 mm, respectively. We confirmed that the estimated irrigation water amounts responded to the levels of SGW table depths (free-drainage to SGW − 100 cm), respectively. When the SGW table depth of − 100 cm was set in the modeling condition, the crop production was slightly underestimated with the relatively low IWA compared to those of free-drainage, SGW − 200 cm and SGW − 150 cm. This underestimation might be inferred that the soil matrix from 0 (soil surface) to − 100 cm can retain the wet conditions (indicating out of the water stress ranges) because of the upward (to the soil surface) fluxes of groundwater induced by the root moisture extraction force and capillary phenomena. Capillary phenomena characteristics that govern soil water flow through the soil profile are highly varied based on various soil textures indicating that investigations of soil characteristics at field regions are required for establishing efficient irrigation schedule and management plans. Considering that field-scale investigations are time-consuming/expensive and requires intensive labors for large regions, our proposed ISM model can be expanded to large regions by combining RS soil moisture products such as advanced microwave scanning radiometer for Advanced Microwave Scanning Radiometer for EOS-AMSR-E (Njoku 2008), Soil Moisture Ocean Salinity-SMOS (Kerr et al. 2001), Soil Moisture Active and Passive-SMAP (Entekhabi et al. 2010), etc. Thus, these findings showed the potential application of ISM model in large spatial domains by combining the RS soil moisture data and near-surface soil moisture data assimilation schemes.
We estimated the long-term wheat crop productions from 1989–1990 to 2008–2009 (18 years) at the RAC site with the long-term weather conditions in Fig. 11. Note that the years of 2001–2002 and 2005–2006 were excluded because of the weather data missing. Also, the same sowing/harvesting dates in the winter season of 2007–2008 were used for the whole years. The sowing/harvesting dates need to be adjusted based on the weather conditions, but these are not considered for the long-term crop production estimations owing to the lack of data. Overall, the precipitation amounts highly influence on the yearly crop productions compared to the relatively monotonous changes of yearly (growing period) averaged solar radiation (15,561.4/1412.5 for average/SD, KJ/m2) and temperature (18.5/3.1 for average/SD, °C) values. Especially, the ISM-derived IWA, II and WS values showed the yearly variations and were inversely correlated with the occurrence days of effective precipitation (more than 5 mm). These findings indicated that optimal combinations of the crop and water management variables are considerably affected by yearly environmental conditions. Although uncertainties exist in these results, the ISM model improved the sustainability of crop productions with the optimal IWA, II and WS variables responding to the environmental conditions. Thus, these results showed the possibility that the ISM model can be used to establish long-term crop and irrigation management plans combining with climate change scenarios in future.
Conclusions
In this study, the ISM model was developed to improve the sustainability of optimal crop productions by providing timely efficient irrigation schedules and management under agricultural drought. The ISM model finds the optimal combination of IWA and II by adjusting WS responding to the environmental conditions based on the simulation–optimization scheme. The field validation and synthetic numerical experiments were conducted to test the proposed ISM model at the field region. The numerical experiments were conducted as follows: (1) Case Study 1: field validation experiments and (2) Case Study 2: synthetic numerical experiments, respectively. The RAC site in Faisalabad (at Punjab, Pakistan) was selected to conduct the Case Studies.
In Case Study 1, the ISM model improved both the sustainability of crop productions and water productivity during the periods of 2007–2008 and 2008–2009 compared to those of the conventional irrigation method. Especially, the estimated WP value was significantly increased for 2008–2009 indicating that the ISM model contributed to both the sustainability of optimal crop production and irrigation water use efficiency. These findings demonstrated that the ability that can consider the relationship between crop (WS) and water management (IWA and II) variables allows the ISM model to respond immediately to environmental changes.
In Case Study 2, the ISM model was tested and improved the crop productions under various environment conditions (weather, soils, crops and BBCs). The L, SCL and SL soils having the high infiltration rate were extremely vulnerable to water deficit, but air circulation and diffusion through large pore space for these soils (L, SCL and SL) are favorable to crop growth with the appropriate crop and water management plans suggested by ISM. Also, the CP values were improved with the presence of SGW becoming closer to the soil surface in the soil profile, because crop root can access easily to SGW due to the root moisture extraction force and capillary phenomena under the rain-fed condition. However, the ISM-based crop productions usually reached to the maximum levels with the presences of BBC (free-drainage and SGW − 200 ~ − 100 cm). Additionally, the long-term wheat crop productions from 1989–1990 to 2008–2009 were estimated and showed the variations along the yearly precipitation changes. Especially, the yearly crop productions were inversely correlated with the occurrence days of effective precipitation (more than 5 mm). These results supported that optimal combinations of the crop and water management variables are highly influenced by environmental conditions.
In this study, the applicability of newly developed ISM model was tested to improve the crop productions for various environment conditions under different hydro-climatic conditions. Although the field validation experiments were conducted at one region (the RAC site) in Faisalabad (Pakistan), our findings demonstrated the robustness of the ISM model for improving the sustainability of optimal (maximum) crop productions at field regions. Thus, our proposed ISM model can be used for establishing efficient irrigation schedule and management plans under agricultural drought.
References
Ahmad MUD, Bastiaanssen WG, Feddes RA (2002) Sustainable use of groundwater for irrigation: a numerical analysis of the subsoil water fluxes irrigation and drainage. J Int Comm Irrig Drain 51(3):227–241. https://doi.org/10.1002/ird.59
Aswathanarayana U (2001) Water resources management and the environment. Balkema, Rotterdam, The Netherlands
Bandaragoda DJ (1998) Design and practice of water allocation rules: lessons from warabandi in Pakistan’s Punjab. Research Report No. 17. International Water Management Institute, Colombo, Sri lanka
Belmans C, Wesseling JG, Feddes RA (1983) Simulation of the water balance of a cropped soil. SWATRE J Hydrol 63:271–286
Berkoff J, Huppert W (1987) Matching crop water requirements in large systems with a variable water supply: Experiments in India, OID-IIMI Irrigation Management Network Paper 87/3d. Overseas Development Institute, London
Bessembinder JJE, Leffelaar PA, Dhindwal AS, Ponsioen TC (2005) Which crop and which drop, and the scope for improvement of water productivity. Agric Water Manag 73:113–130. https://doi.org/10.1016/j.agwat.2004.10.004
Bonfante A, Monaco E, Manna P, De Mascellis R, Basile A, Buonanno M et al (2019) LCIS DSS—An irrigation supporting system for water use efficiency improvement in precision agriculture: a maize case study. Agric Syst 176:102646. https://doi.org/10.1016/j.agsy.2019.102646
Bos MG (1997) Performance indicators for irrigation and drainage. Irrig Drain Syst 11(2):119–137
Bos MG, Nugteren J (1990) On irrigation efficiencies. (4th ed. p. 141), ILRI Publication 19, Int. Institute for Land Reclamation and Improvement Wageningen, The Netherlands
Boughdiri A, Daghari H, Saidi A (2014) The water use efficiency for different varieties of wheat and the effect of supplemental irrigation in the semi-arid regions of Tunisia. IJAIR 3:326–330
Dai A, Zhao T, Chen J (2018) Climate change and drought: a precipitation and evaporation perspective. Curr Clim Change Rep 4:301–312. https://doi.org/10.1007/s40641-018-0101-6
van Dam JC, Huygen J, Wesseling JG, Feddes RA, Kabat P, van Waslum PEV, Groenendjik P, van Diepen CA (1997) Theory of SWAP version 2.0: Simulation of water flow and plant growth in the soil water-atmosphere-plant environment. Tech. Doc. 45, Wageningen Agric.Univ., and DLO Winand Staring Cent., Wageningen, Netherlands
van Dam JC (2000) Field-scale water flow and solute transport. SWAP model concepts, parameter estimation and case studies. Ph.D Dissertation, Wageningen Univ., Wageningen, The Netherlands
Dilsher SA, Khalaf AS (2014) Effect of preceding crops and supplementary irrigation on yield and yield components of two varieties of common wheat (Triticum aestivum L.). Am J Exp Agric 4:1944–1957. https://doi.org/10.9734/AJEA/2014/10889
Domínguez A, Martínez-Navarro A, López-Mata E, Tarjuelo JM, Martínez-Romero A (2017) Real farm management depending on the available volume of irrigation water (part I): financial analysis. Agric Water Manag 192:71–84. https://doi.org/10.1016/j.agwat.2017.06.022
Droogers P, Bastiaanssen WGM, Beyazgül M, Kayam Y, Kite GW, Murray-Rust H (2000) Distributed agro-hydrological modeling of an irrigation system in Western Turkey. Agric Water Manag 43(2):183–202. https://doi.org/10.1016/S0378-3774(99)00055-4
Du Y, Niu W, Gu X, Zhang Q, Cui B, Zaho Y (2018) Crop yield and water use efficiency under aerated irrigation: A meta-analysis. Agric Water Manag 210:158–164. https://doi.org/10.1016/j.agwat.2018.07.038
Efron B (1982) The jackknife, the bootstrap and other resampling plans. Society for industrial and applied mathematics, Philadelphia, PA
Entekhabi D, Njoku EG, O’Neill PE, Kellogg KH, Crow WT, Edelstein WN, Entin JK, Goodman SD, Jackson TJ, Johnson J, Kimball J, Piepmeier JR, Koster RD, Martin N, McDonald KC, Moghaddam M, Moran S, Reichle R, Shi JC, Spencer MW, Thurman SW, Tsang L, Zyl JV (2010) The soil moisture active passive (SMAP) mission. Proc IEEE 98(5):704–716. https://doi.org/10.1109/jproc.2010.2043918
Evans RG, Sadler EJ (2008) Methods and technologies to improve efficiency of water use. Water Resour Res 44:1–15. https://doi.org/10.1029/2007WR006200
Feddes RA, Kowalik PJ, Zarandy H (1978) Simulation of field water use and crop yield. Simulation Monographs. Pudoc., Wageningen
Federal Emergency Management Agency (FEMA) (1995) National mitigation strategy: partnerships for building safer communities. Washington, DC, USA
French RJ, Schultz JE (1984) Water use efficiency of wheat in a Mediterranean type environment. I. The relation between yield, water use and climate. Aust J Agric Res 35:743–764
Gobiet A, Kotlarski S, Beniston M, Heinrich G, Rajczak J, Stoffel M (2014) 21st century climate change in the European Alps-A review. Sci Total Environ 493:1138–1151. https://doi.org/10.1016/j.scitotenv.2013.07.050
Goldberg DE (1989) Genetic algorithms in search and optimization and machine learning. Addison-Wesley Publishing Co Inc, USA
Goorahoo D, Carstensen G, Mazzei A (2001) A pilot study on the impact of air injected into water delivered through subsurface drip irrigation tape on the growth and yield of bell peppers. California Agricultural Technology Institute, USA
Gornall J, Betts R, Burke E, Clark R, Camp J, Willett K, Wiltshire A (2002) Implications of climate change for agricultural productivity in the early twenty-first century. Philos Trans R Soc B Biol Sci 365:2973–2989. https://doi.org/10.1098/rstb.2010.0158
Hanif M (2008) Weather and wheat development during 2007–08 at R.A.M.C. Faisalabad. National Agromet Center Islamabad, Faisalabad, pp 1–12
Holland JH (1975) On quantifying agricultural and water management practices from low spatial resolution RS data using genetic algorithms: a numerical study for mixed pixel environment. Adv Water Resour 28:856–870
Ines AVM, Honda K (2005) On quantifying agricultural and water management practices from low spatial resolution RS data using genetic algorithms: a numerical study for mixed pixel environment. Adv Water Resour 28:856–870. https://doi.org/10.1016/j.advwatres.2004.11.015
Ines AVM, Mohanty BP (2008) Near-surface soil moisture assimilation to quantify effective soil hydraulic properties using genetic algorithm I conceptual modeling. Water Resour Res. https://doi.org/10.1029/2007WR005990
Ines AVM, Gupta AD, Loof R (2002) Application of GIS and crop growth models in estimating water productivity. Agric Water Manag 54(3):205–225. https://doi.org/10.1016/S0378-3774(01)00173-1
Ines AVM, Honda K, Gupta AD, Droogers P, Clemente RS (2006) Combining remote sensing-simulation modeling and genetic algorithm optimization to explore water management options in irrigated agriculture. Agric Water Manag 83(3):221–232. https://doi.org/10.1016/j.agwat.2005.12.006
Jin X, Feng H, Zhu X, Li Z, Song S, Song X, Yang G, Xu X, Guo WS (2014) Assessment of the aqua crop model for use in simulation of irrigated winter wheat canopy cover, biomass, and grain yield in the North China Plain. PLoS ONE 9:1–11. https://doi.org/10.1371/journal.pone.0086938
Kerr YH, Waldteufel P, Wigneron JP, Martinuzzi JM, Font J, Berger M (2001) Soil moisture retrieval from space: the soil moisture and ocean salinity (SMOS) mission. IEEE Trans Geosci Remote Sens 39(8):1729–1735. https://doi.org/10.1109/36.942551
Kroes JG, van Dam JC, Huygen J, Vervoort RW (1999) User’s guide of SWAP version 2.0; Simulation of water, solute transport, and plant growth in the soil-atmosphere-plant environment. Wageningen Agricultural University, Report 81, DLO Winand Staring Centre. Technical Document, pp 53–128
Lee T, Shin Y (2016) Estimation of irrigation water amounts for farm products based on various soil physical properties and crops. J Korean Soc Agric Eng 58(6):1–8. https://doi.org/10.5389/KSAE.2016.58.6.001
Leij FJ, Alves WJ, van Genuchten MT, Williams JR (1999) The UNSODA unsaturated soil hydraulic database. In: Van Genuchten MT et al (eds) Characterization and measurement of the hydraulic properties of unsaturated porous media. University of California, Riverside, USA, pp 1269–1281
Lettenmaier DP, McCabe G, Stakhiv EZ (1996) Global climate change: effects on hydrologic cycle. In: Part V, Mays LW (eds) Chapter 29 water resources handbook. McGraw-Hill, New York, p 29.3-29.33
Lv LH, Wang HJ, Jia XL, Wang ZM (2011) Analysis on water requirement and water-saving amount of wheat and corn in typical regions of the North China Plain. Front Mech Eng China 5:556–562. https://doi.org/10.1007/s11703-011-1149-4
M’hamed HC, Rezig M, Naceur MB (2015) Deficit irrigation of durum wheat (triticum durum desf): Effects on total dry matter production, light interception and radiation use efficiency under different nitrogen rates. Sustain Agric Res 4:26–40. https://doi.org/10.5539/sar.v4n1p26
Mishra AK, Singh VP (2009) Analysis of drought severity-area-frequency curves using a general circulation model and scenario uncertainty. J Geophys Res 114(D6):D06120. https://doi.org/10.1029/2008JD010986
Molden D (1997) Accounting for water use and productivity, SWIM Paper. Int. Irrigation Management Institute Colombo, Sri Lanka, 1, 16
Morita S, Abe J (1996) Development of root systems in wheat and rice. In: Ito O, Katayama K, Johansen C, Kumar Rao JVDK, Adu-Gyamfi JJ, Rego TJ (eds) Roots and nitrogen in cropping systems of the semi-arid tropics. JIRCAS, Tsukuba, pp 185–198
Mualem Y (1976) A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Resour Res 12:513–522
Njoku E (2008) AMSR-E/Aqua Daily L3 Surface soil moisture, interpretive parameters, and QC EASE-Grids V002, [list the dates of the data used]. Boulder, Colorado USA: National Snow and Ice Data Center. Digital media
Olesen JE, Bindi M (2002) Consequences of climate change for European agricultural productivity, land use and policy. Eur J Agron 16(4):239–262. https://doi.org/10.1016/S1161-0301(02)00004-7
Pardales JR, Dionisio MB, Akira Y, Morio I, Esquibel CB (2002) The effect of fluctuations of soil moisture on root development during the establishment phase of sweetpotato. Plant Prod Sci 3(2):134–139. https://doi.org/10.1626/pps.3.134
Passioura J (2006) Increasing crop productivity when water is scarce—from breeding to field management. Agric Water Manag 80:176–196. https://doi.org/10.1016/j.agwat.2005.07.012
Polinova M, Salinas K, Bonfante A, Brook A (2019) Irrigation optimization under a limited water supply by the integration of modern approaches into traditional water management on the cotton fields. Remote Sens 11(18):2127. https://doi.org/10.3390/rs11182127
Prasad PVV, Staggenborg SA, Ristic Z (2008) Impacts of drought and/or heat stress on physiological, developmental, growth, and yield processes of crop plants. In: Response of crops to limited water: Understanding and modeling water stress effects on plant growth processes. American Society of Agronomy, Crop Science Society of America, Soil Science Society of America, Madison, WI, Vol 1, pp. 301–356. https://doi.org/10.2134/advagricsystmodel1.c11
Saeedipour S, Moradi F (2011) Effect of drought at the post-anthesis stage on remobilization of carbon reserves and some physiological changes in the flag leaf of two wheat cultivars differing in drought resistance. J Agric Sci 3:81–92. https://doi.org/10.5539/jas.v3n3p81
Sarwar A, Bastiaanssen WGM, Boers MT, van Dam JC (2000) Evaluating drainage design parameters for the fourth drainage project, Pakistan by using SWAP model: I Calibration. Irrig Drain Syst 14(4):257–280. https://doi.org/10.1023/A:1006468905194
Seckler D (1996) The new era of water resources management: from ‘dry’ to ‘wet’ water savings. Int. Water Management Institute, Colombo, Sri Lanka, 17
Shin Y, Jung Y (2014) Development of irrigation water management model for reducing drought severity using remotely sensed soil moisture footprints. J Irrig Drain Eng 140(7):04014021. https://doi.org/10.1061/(ASCE)IR.1943-4774.0000736
Shin Y, Mohanty BP (2013) Development of a deterministic downscaling algorithm for remote sensing soil moisture footprint using soil and vegetation classifications. Water Resour Res. https://doi.org/10.1002/wrcr.20495
Singh R, Jhorar RK, van Dam JC, Feddes RA (2006) Distributed ecohydrological modeling to evaluate irrigation system performance in Sirsa district, India: II. Impact of viable water management scenarios. J Hydrol 329:714–723. https://doi.org/10.1016/j.jhydrol.2006.03.016
Sugimoto H, Koesmaryono Y, Nakano R (2000) Effects of excess moisture in the soil at different stages of development on the growth and seed yield of soybean. Pak J Biol Sci 3(9):1465–1467. https://doi.org/10.3923/pjbs.2000.1465.1467
Sun HY, Liu CM, Zhang XY, Shen YJ, Zhang YQ (2006) Effects of irrigation on water balance, yield and WUE of winter wheat in the North China Plain. Agric Water Manag 85:211–218. https://doi.org/10.1016/j.agwat.2006.04.008
Svoboda M et al (2002) The drought monitor. Bull Am Meteorol Soc 83:1181–1190
Tolomio M, Casa R (2020) Dynamic crop models and remote sensing irrigation decision support systems: a review of water stress concepts for improved estimation of water requirements. Remote Sens 12(23):3945. https://doi.org/10.3390/rs12233945
van Genuchten MT (1980) A closed-form equation foe predicting the hydraulic conductivity of unsaturated soils. Soil Sci Soc Am J 44:892–898
Wesseling JG, Kroes JG (1998) A global sensitivity analysis of the model SWAP. Rep. 160, DLO Winand Staring Centre, Wageningen, The Netherlands
Funding
Korea Ministry of Environment, The SS (Surface Soil conservation, Yongchul Shin, management) projects, Yongchul Shin, 2019002820002, Yongchul Shin
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Lee, T., Jang, W., Chun, B. et al. Development of irrigation schedule and management model for sustaining optimal crop production under agricultural drought. Paddy Water Environ 21, 31–45 (2023). https://doi.org/10.1007/s10333-022-00911-9
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10333-022-00911-9