
Conventional farming practices not only constrained food security due to low yield but also threatened the ecosystem by causing groundwater decline and groundwater nitrate contamination. A twoear field experiment was conducted at the research station of North China University of Water Resources and Electric Power, Zhengzhou. The WHCNS model was used to simulate grain yield, water and nitrogen fertilizer use efficiencies (WUE and FNUEs) of spring maize under border irrigation method, drip irrigation, and rainfed conditions. In addition, a scenario analysis was also performed on different dry and rainy seasons to assess the long-term impact of rainfall variability on spring maize from 2000–2017. The result showed that the model precisely simulated soil water content, N concentration, crop biomass accumulation, and grain yield. The maximum and minimum range of relative root mean squire error (RRMSE) values were 0.5–36.0% for soil water content, 14.0–38.0% for soil nitrate concentrations, 19.0–24.0% for crop biomass and 1.0–2.0% for grain yield, respectively under three irrigation methods. Both the index of agreement (IA) and Pearson correlation coefficient (r) values were close 1. We found the lowest grain yield from the rainfed maize, whereas the drip irrigation method increased grain yield by 14% at 40% water saving than border irrigation method for the two years with the 11% lower evaporation and maintained transpiration rate. Moreover, the drip irrigated maize had a negligible amount of drainage and runoff, which subsequently improved WUE by 27% in the first growing season and 16% in the second rotation than border irrigation. The drip irrigated maize also showed 24% higher FNUE. The reason of lower WUE and FNUEs under the border irrigation method was increased drainage amounts and N leaching rates. Furthermore, scenario analysis indicated that the dry season could result in a 30.8% yield decline as compared to rainy season.
Citation: Shu Xu, Yichang Wei, Abdul Hafeez Laghari, Xianming Yang, Tongchao Wang. Modelling effect of different irrigation methods on spring maize yield, water and nitrogen use efficiencies in the North China Plain[J]. Mathematical Biosciences and Engineering, 2021, 18(6): 9651-9668. doi: 10.3934/mbe.2021472
[1] | Favour Adenugba, Sanjay Misra, Rytis Maskeliūnas, Robertas Damaševičius, Egidijus Kazanavičius . Smart irrigation system for environmental sustainability in Africa: An Internet of Everything (IoE) approach. Mathematical Biosciences and Engineering, 2019, 16(5): 5490-5503. doi: 10.3934/mbe.2019273 |
[2] | Chenyang Zhan, Min Tian, Yang Liu, Jie Zhou, Xiang Yi . A novel greedy adaptive ant colony algorithm for shortest path of irrigation groups. Mathematical Biosciences and Engineering, 2022, 19(9): 9018-9038. doi: 10.3934/mbe.2022419 |
[3] | Edwin Duque-Marín, Alejandro Rojas-Palma, Marcos Carrasco-Benavides . A soil water indicator for a dynamic model of crop and soil water interaction. Mathematical Biosciences and Engineering, 2023, 20(8): 13881-13899. doi: 10.3934/mbe.2023618 |
[4] | Todorka Glushkova, Stanimir Stoyanov, Lyubka Doukovska, Jordan Todorov, Ivan Stoyanov . Modeling of an irrigation system in a virtual physical space. Mathematical Biosciences and Engineering, 2021, 18(5): 6841-6856. doi: 10.3934/mbe.2021340 |
[5] | Gianluca D'Antonio, Paul Macklin, Luigi Preziosi . An agent-based model for elasto-plastic mechanical interactions between cells, basement membrane and extracellular matrix. Mathematical Biosciences and Engineering, 2013, 10(1): 75-101. doi: 10.3934/mbe.2013.10.75 |
[6] | Zixiao Xiong, Xining Li, Ming Ye, Qimin Zhang . Finite-time stability and optimal control of an impulsive stochastic reaction-diffusion vegetation-water system driven by L$ {\rm \acute{e}} $vy process with time-varying delay. Mathematical Biosciences and Engineering, 2021, 18(6): 8462-8498. doi: 10.3934/mbe.2021419 |
[7] | Alejandro Rincón, Fredy E. Hoyos, John E. Candelo-Becerra . Comparison, validation and improvement of empirical soil moisture models for conditions in Colombia. Mathematical Biosciences and Engineering, 2023, 20(10): 17747-17782. doi: 10.3934/mbe.2023789 |
[8] | Salamida Daudi, Livingstone Luboobi, Moatlhodi Kgosimore, Dmitry Kuznetsov . Dynamics for a non-autonomous fall armyworm-maize interaction model with a saturation functional response. Mathematical Biosciences and Engineering, 2022, 19(1): 146-168. doi: 10.3934/mbe.2022008 |
[9] | Chenxi Dai, ZhiWang, Weiming Wang, Yongqin Li, Kaifa Wang . Epidemics and underlying factors of multiple-peak pattern on hand, foot and mouth disease inWenzhou, China. Mathematical Biosciences and Engineering, 2019, 16(4): 2168-2188. doi: 10.3934/mbe.2019106 |
[10] | Sunmi Lee, Romarie Morales, Carlos Castillo-Chavez . A note on the use of influenza vaccination strategies when supply is limited. Mathematical Biosciences and Engineering, 2011, 8(1): 171-182. doi: 10.3934/mbe.2011.8.171 |
Conventional farming practices not only constrained food security due to low yield but also threatened the ecosystem by causing groundwater decline and groundwater nitrate contamination. A twoear field experiment was conducted at the research station of North China University of Water Resources and Electric Power, Zhengzhou. The WHCNS model was used to simulate grain yield, water and nitrogen fertilizer use efficiencies (WUE and FNUEs) of spring maize under border irrigation method, drip irrigation, and rainfed conditions. In addition, a scenario analysis was also performed on different dry and rainy seasons to assess the long-term impact of rainfall variability on spring maize from 2000–2017. The result showed that the model precisely simulated soil water content, N concentration, crop biomass accumulation, and grain yield. The maximum and minimum range of relative root mean squire error (RRMSE) values were 0.5–36.0% for soil water content, 14.0–38.0% for soil nitrate concentrations, 19.0–24.0% for crop biomass and 1.0–2.0% for grain yield, respectively under three irrigation methods. Both the index of agreement (IA) and Pearson correlation coefficient (r) values were close 1. We found the lowest grain yield from the rainfed maize, whereas the drip irrigation method increased grain yield by 14% at 40% water saving than border irrigation method for the two years with the 11% lower evaporation and maintained transpiration rate. Moreover, the drip irrigated maize had a negligible amount of drainage and runoff, which subsequently improved WUE by 27% in the first growing season and 16% in the second rotation than border irrigation. The drip irrigated maize also showed 24% higher FNUE. The reason of lower WUE and FNUEs under the border irrigation method was increased drainage amounts and N leaching rates. Furthermore, scenario analysis indicated that the dry season could result in a 30.8% yield decline as compared to rainy season.
Maize is a major crop of the North China Plain (NCP), which accounts for approximately more than 31% of the national grain production. To increase yield per unit area, farmers in this region typically apply 5–7 times of irrigation to crop via the conventional border irrigation method, with each time pumping of more than 60 mm water. As a result, annual groundwater withdrawal ranges from 600 to 800 mm yr–1 [1]. According to a report, over 7.6 million tubewells are used to pump groundwater across the region [2]. Consequently, excessive water pumping, continuously causing groundwater to drop at a rate of more than 1 m yr–1 [3]. The first aquifer is completely depleted, and the second aquifer is about to be depleted, as farmers pump water from the ground faster than it can refill. Furthermore, because of little precipitation, annual natural water replenishment is uncertain [4]. In the near future, the situation may become more difficult. For example, the water table in the counties surrounding Beijing is declining year after year, and rapid urbanisation has increased demand for water while reducing its availability [5–7]. Presently, the water shortage drastically limiting crop yield [8].
On the other hand, excessive use of N fertilizer in crop production posed a hazard to both terrestrial and aquatic ecosystems. A typical plant can only utilize half of the applied N fertilizer and remaining enters in water bodies via drainage [9]. In an estimate, 249 kg N ha–1 year–1 exceeds the crop N uptake capacity [10]. The N fertilizer is highly water–soluble, it transports with water on the wide range of soil types [11]. Many researchers showed that about 20% of N fertilizer loss takes place due to leaching [12]. Therefore, fertilizer N use efficiency in the region is often less than 20 to 25%, because farmers apply N fertilizer through broadcasting, which is an inefficient method as compared to modern techniques, such as application of water and nutrients through drip irrigation system [13,14].
The drip irrigation system consists of small diameter plastic pipes with outlets that drip a low rate of water into the soil. Unlike traditional border irrigation system, which saturate the entire soil, water is supplied near to the root-zone of plants, wetting only the region of the soil in which the roots grow. Water is applied more frequent with the drip irrigation system, ensuring full utilization of water stored in the soil profile and a very favourable high moisture content in the soil for plants to flourish. Despite the fact that drip irrigation system is only practiced on a small scale, it has a great potential for increasing production [15] improving photosynthesis rate [16], reducing evaporation, deep drainage, and save irrigation when compared to border irrigation method [14,17]. Drip irrigation significantly reduces evaporation without affecting transpiration rate of wheat crop [18]. Drip irrigation also enhances WUE of different crops and reduces water loss via drainage and runoff [19,20]. There are also many other advantages of drip irrigation system. Drip irrigation, in particular, maintains soil aeration, which promotes plant root growth, increases leaf area index, and improves grain quality. Despite these advantages, the drip is not well adopted in the agricultural system of the region. Due to water shortages in most areas, maize is frequently cultivated under rainfed conditions after border irrigation.
Many researchers have focused on spring maize and examined its yield potentiality under the impact of different agronomic methods to offset the effect of water shortage on crop production in the study area [1]. Still, modelling based information on water and nitrogen balance of spring maize under different irrigation methods is limited. In present study, we used WHCNS model as a tool to determine the N fate and water consumption, and to answer how much the WUE and FNUE changes under border irrigation method, drip irrigation and rainfed conditions. The model has been successfully applied in the NCP [21]. Wang et al. [22] reported that simulated soil water content, soil N content, crop biomass and N uptake matched well with the measured data.
The objectives were to; (i) compare grain yield, water consumption, N fates, WUE and FNUE of spring maize under different irrigation methods, (ii) evaluate simulation performance of WHCNS model, (iii) study effect of long-term dry and rainy seasons on spring maize.
The research station is located in Yuanyang county (34.47° N, 113.38° E) of Zhengzhou, Henan province, China. The study area experiences warm temperate, sub–humid, and continental monsoon weather conditions (Figure 1). The mean annual temperature is 14.4 ℃ and rainfall is about 641 mm. Maximum rainfall (60–70%) occurs during spring maize growing season. Soil of experimental site has sandy loam texture (USDA classification), pH 8.2, organic matter 11.2 g kg–1, nitrate (N) 76.5 g kg–1, phosphorus (P) 8.2 g kg–1 and exchangeable (K) 67.3 g kg–1 [23,24]. The layer-wise basic soil properties are given in Table 1.
Depth (cm) |
BD (g cm–3) |
Particle fraction (%) | Soil texture (USDA) |
pH | θs (cm3 cm–3) |
θfc (cm3 cm–3) |
θwp (cm3 cm–3) |
Ks (cm d–1) |
||
Sand | Silt | Clay | ||||||||
0–20 | 1.31 | 52 | 23 | 25 | Sandy clay loam | 8.2 | 0.42 | 0.31 | 0.14 | 42.2 |
20–40 | 1.34 | 54 | 25 | 21 | Sandy clay loam | 8.0 | 0.40 | 0.32 | 0.15 | 33.7 |
40–60 | 1.30 | 57 | 23 | 20 | Sandy clay loam | 7.9 | 0.41 | 0.31 | 0.13 | 31.4 |
60–80 | 1.35 | 58 | 22 | 20 | Sandy clay loam | 7.5 | 0.40 | 0.28 | 0.14 | 50.5 |
80–100 | 1.37 | 57 | 22 | 21 | Sandy clay loam | 7.8 | 0.39 | 0.27 | 0.14 | 51.2 |
Note: BD; Bulk density, θs; Soil saturated water content, θfc; Field capacity, θwp; Permanent wilting point and Ks; Saturated hydraulic conductivity. |
A two-year field experiment was conducted from 2018–2019. The experimental design was a randomized complete block design (RCBD), with four replications. There were three treatments, including conventional border irrigation method, drip irrigation and rainfed conditions. The plot size was 108 m2 (12 m long × 9 m wide). Each drip irrigation plot consisted of lateral pipes that delivered water directly to the root zones of the crops, while the sub–main pipeline system supplied water to the laterals. With a flow rate of 1.2 L h-1, the inner diameter of lateral drip pipes was 1.60 cm. For border irrigation and drip irrigation methods, there were separate water outlets and water meters to record the amount of irrigation. The irrigation system was controlled by a computer and a smartphone application [25].
Spring maize variety Luobai–3 was sown on 13 April in 2018 and 14 April in 2019 at a row spacing of 60 cm with a seed rate of 30 kg ha–1. Total amount of irrigation under drip irrigation method was 90 mm (20–30 mm each time). Irrigation was applied when the field moisture level fell below 50% of the field capacity and irrigated to 80% of the field capacity. In the border irrigation method, the total irrigation amount was 150 mm (50 mm each time). For the strong seedling establishment, the full dose of phosphorus (25 kg P ha−1) and potassium (21 kg K ha−1) fertilizers was applied as a basal fertilizer. Superphosphate was used as a source of P fertilizer, and potassium chloride was used as a source of K fertilizer. For each treatment, the total amount of nitrogen (N) was 200 kg ha–1, from which 65% of the N fertilizer was incorporated in the soil before sowing. The plant density was 6–8 plants m2 at harvest time. Weed control, insect pest control, and disease control measures were all according to the local farmer practices. Manual hand hoeing was performed during the 2–3 leaf stage and mechanical plough operation was performed at 4–6 leaf stage. The crop was harvested on 24–25 September.
The soil bulk density and soil particle fraction were determined by the cutting ring method [26] and using hydrometer [27], respectively. Soil saturated and residual water contents were determined with the water retention curve [28]. Soil saturated hydraulic conductivity was measured using KSAT instrument. While soil water content was measured from different soil depths (0–100 cm at 20 cm intervals) in each experimental plot using on–field installed WiTu technology tubular moisture sensors. This sensor is new, modern, and most accurate to determine soil water content [29]. While, for determination of soil nitrate concentrations, soil samples were sieved with 4 mm sieves, extracted with 1 mol L–1 KCL (1:10, soil/water), and filtered. Ready samples were analyzed by a continuous flow analyzer (AA3, Bran and Luebbe, Germany). The grain yields and crop biomass were measured at harvest and different growth stages [30].
Integrated WHCNS (Soil-Water-Heat-Carbon-Nitrogen-Simulator) model simulates crop growth and yield as a function of the soil–plant–atmosphere dynamics. The model consists of five primary modules, including (i) soil water, (ii) soil temperature (iii), carbon, (iv) nitrogen, and (v) crop growth. In the functional system of model, the simulation methods for soil water and temperature are imported from HYDRUS1D [31] and RZWQM models [32]. The FAO Penman–Monteith equation is used to compute the reference ET0 [33]. While water infiltration and redistribution are modelled via modified Green–Ampt [34] and Richard's [35] equations, respectively. Runoff is estimated using the SCS curve method [36]. Soil carbon and nitrogen cycles are adopted from the Daisy model [37] and crop growth and yields were projected using an improved PS123 model [38].
The model runs on a daily time step and requires latitude and longitude of site, soil saturated water content (θs), water content at field capacity (θfc), permanent wilting point (θwp), saturated hydraulic conductivity (Ks), soil C/N ratio (kg h–1), root, shoot and leaf weight (g), plant type (C3 or C4), sowing date, seed rate, and fertilizer dose, harvesting date, as well as daily rainfall, relative humidity, maximum, minimum and mean air temperature. The soil data, which were calibrated, is presented in Table 1, and basic crop parameters are given in Table 2. The meteorological data used in the model were acquired from an on–field established weather station.
Parameters and definitions of parameters | Values |
Tbase defined base temperature for crop growth | 8 ℃ |
Tsum defined as effective accumulated temperature | 2249℃ |
Ke define as extinction coefficient (0–1) | 0.7 |
Kc_ini defined crop coefficients at initial stage (0–1) | 0.64 |
Kc_mid defined as crop coefficients at middle stage (0–1) | 1.26 |
Kc_end defined as crop coefficients at final stage (0–1) | 0.7 |
SLA_max defined as maximum specific leaf area (m2 kg–1) | 29 |
SLA_min defined as minimum specific leaf area (m2 kg–1) | 20 |
AMAX defined as maximum assimilation rate (kg ha–1) | 50 |
R_max defined as maximum soil depth contributing to N transformation (cm) | 142 |
Vn defined as maximum nitrification rate (mg L–1 d–1) | 100 |
Kn defined as half saturation constant for nitrification (g m−3) | 200 |
Kd defined as empirical proportionality factor for denitrification (−) | 0.2 |
Ad defined as first–order kinetic coefficient for denitrification (d−1) | 0.05 |
Kv defined as first–order kinetic coefficient for ammonia volatilization denitrification (d−1) | 0.03 |
To evaluate the WHNCS model simulations goodness of fit, the Relative Root Mean Square Error (RRMSE), Index of Agreement (IA) and Pearson Correlation Coefficient (r) were used.
RRMSE=RMSEˉO×100 | (1) |
IA=1−∑ni=1(Oi−Pi)2∑ni=1(|Oi−ˉO|+|Pi−ˉO|)2 | (2) |
r=∑ni=1(Oi−ˉO)×(Pi−ˉPi)√∑ni=1(Oi−ˉO)2√∑ni=1(Pi−ˉPi)2 | (3) |
where Oi is observed (or measured) value, Piis predicted value, and
Statistically, RRMSE is a good indicator of model accuracy, lower the RRMSE (%), better the model performance [39], and IA valueranges from 0–1, close to 1, IA value represents a perfect match of simulated data with measured, and 0 specify no agreement [40]. Whereas Pearson correlation coefficient r value equal to either 1 or –1 illustrates perfect correlation between simulated and measured data [41].
In addition, the analysis of variance (ANOVA) was applied using Statistix 8.1 software (Analytical Software, Tallahassee, FL, USA) to calculate the overall effect of border, drip and rainfed conditions on spring maize. The significant differences among irrigation methods were determined using Fisher's LSD at P < 0.05 through the All–pairwise comparison method [42]. Furthermore, water balance and water use efficiency (WUE), nitrogen balance and nitrogen fertilizer use efficiency (FNUE), savings, and yield variations (%) were calculated using the following equations.
Wbal=(I+P)−(ET+D+R) | (4) |
WUE=GY(kgha−1)ET(mm) | (5) |
WS(%)=WU(BI−DIorRF)WU(BI)×100 | (6) |
FNUE=GY(kgha−1)AppliedN(kgha−1) | (7) |
YV(%)=GYBI−GYDIorRFGYBI | (8) |
where, Wbal is water balance, I is Irrigation (mm); P is precipitation (mm); ET is evapotranspiration (mm) and D is drainage (mm), R is runoff (mm); WUE is water use efficiency; GY is grain yield (kg ha–1); WS is water-saving; WU is water use (mm); BI is border irrigation; DI is drip irrigation; RF is rainfed [43]; FNUE is fertilizer nitrogen use efficiency; YV is yield variation.
To extensively explore the effect of long-term dry and rainy seasons on spring maize yield, water movement, and N fate, we performed 18 years scenario analysis from 2000–2017. The border irrigation method was used as a basic treatment. The dry and rainy years were differentiated according to the method used by Murray-Tortarolo et al. [44] and Gao et al. [45]. Years with annual seasonal rainfall greater than the mean annual seasonal rainfall over 2000–2017 were identified as rainy years, while those that were less than the mean annual seasonal rainfall were separated as dry years. There was total 10 dry years and 8 rainy years (Figure 2). There is no broadly recognized criteria to justify rainy and dry years; therefore, separation is relative [45]. To avoid complexity, we used the long term mean annual season rainfall as a threshold level in the study. Figure 2 shows the differentiated dry and rainy seasons.
The model was calibrated using measured data (season 2018) of soil water content (cm3 cm–3), N concentrations (mg kg–1), crop biomass (kg ha–1), and grain yield under border irrigation method, drip irrigation and rainfed conditions, and then validated using data from the next year (season 2019). Highly sensitive crop and soil indicators were adjusted using a trial–and–error method during calibration to reduce the chance of error between simulated and measured data [46]. When compared to prior investigations, statistical evaluation of simulated data using absolute values revealed good model performance [47,48]. The model precisely predicted soil water content, N concentrations, crop biomass, and grain yield. The measured data points fluctuated up–down, close to the projected curve with the acceptable ranges of RRMSE values. The RRMSE varied under effect of different irrigation method and variables. Thus, top to deep soil layer (0–100 cm) RRMSE ranged from 5.0–36.0% for soil water content, 14.0–38.0% for soil nitrate concentrations, 19.0–24.0% for crop biomass and 1.0–2.0% for grain yield. In general, there is no threshold limit of RRMSE value. The model's performance considered good as the RRMSE or RMSE value decreases [49]. The maximum range of RRMSE for soil variables was relatively larger than crop growth and yield variables due to relatively high uncertainty and complex soil processes. Whereas IA values were between 0.53–0.91 for soil water content, 0.75–0.96 for soil nitrate concentrations, 0.97–0.98 for crop biomass, and 0.55–0.64 for grain yield. In the case of correlation coefficients, the r equal to 1 or close 1 existed between simulated and measured crop biomass and grain yields. Similarly, simulated soil water content and N concentrations had moderate to strong correlation with a minimum and maximum range of r corresponding to 0.43–0.98 (Table 3 and Figures 3 and 4). A strong correlation (r > 0.85) between measured and simulated for grain yield and biomass also reported by Shi et al. [46] with simulation error less than 10%.
SWC (cm3 cm–3) | Border | Drip | Rainfed | ||||||
RRMSE | IA | r | RRMSE | IA | r | RRMSE | IA | r | |
0–20 cm | 24.0 | 0.76 | 0.70 | 30.0 | 0.62 | 0.58 | 24.0 | 0.80 | 0.73 |
20–40 cm | 21.0 | 0.62 | 0.47 | 9.0 | 0.77 | 0.70 | 5.0 | 0.91 | 0.84 |
40–60 cm | 26.0 | 0.58 | 0.43 | 35.0 | 0.53 | 0.47 | 14.0 | 0.88 | 0.80 |
60–80 cm | 27.0 | 0.64 | 0.62 | 36.0 | 0.56 | 0.46 | 21.0 | 0.77 | 0.65 |
80–100 cm | 21.0 | 0.61 | 0.69 | 14.0 | 0.73 | 0.58 | 12.0 | 0.82 | 0.78 |
SNC (mg kg–1) | |||||||||
0–20 cm | 27.0 | 0.84 | 0.76 | 22.0 | 0.93 | 0.89 | 36.0 | 0.76 | 0.66 |
20–40 cm | 14.0 | 0.92 | 0.93 | 21.0 | 0.75 | 0.61 | 24.0 | 0.80 | 0.67 |
40–60 cm | 22.0 | 0.86 | 0.77 | 26.0 | 0.81 | 0.77 | 24.0 | 0.88 | 0.81 |
60–80 cm | 21.0 | 0.90 | 0.96 | 27.0 | 0.85 | 0.87 | 33.0 | 0.87 | 0.84 |
80–100 cm | 18.0 | 0.96 | 0.94 | 38.0 | 0.83 | 0.79 | 28.0 | 0.94 | 0.89 |
Biomass (kg ha–1) | 19.0 | 0.98 | 0.98 | 20.0 | 0.98 | 0.97 | 24.0 | 0.98 | 0.97 |
Grain yield (kg ha–1) | 1.0 | 0.64 | 1.00 | 2.0 | 0.61 | 1.00 | 2.0 | 0.55 | 1.00 |
Note: SWC: Soil water content; SNC: Soil nitrate concentrations. |
The farmer-based border irrigation method consumed a high amount of irrigation (300 mm) as compared to drip irrigation and showed a total grain yield of 20184 kg ha–1 (10021 kg ha–1 in 2018 and 10,163 kg ha–1 in 2019). Significantly (P < 0.05) high total grain yield 22949 kg ha–1 (11567 kg ha–1 in 2018 and 11382 kg ha–1 in 2019) obtained from drip irrigated maize at 40% decreased water (180 mm) for the two years. Thus, drip irrigation increased grain yield by 14 and 73% than border irrigation and rainfed maize, respectively. The rainfed maize exhibited total grain yield of 13249 kg ha–1 (6666 kg ha–1 in 2018 and 6584 kg ha–1 in 2019), which was 42% lower than drip irrigated maize because the rainfed plot was highly under water–stress. Throughout the experimental period, the distribution of precipitation was uneven. While vigorous plant growth and development, which promoted root growth, increased biomass accumulation, grain number, and cob length, was the key reason for better production from the drip irrigated maize. Figures 5 and 6 indicate the changes in yield from season to season as well as the physiological response of spring maize under different irrigation methods.
The simulated soil evaporation, crop transpiration, field drainage, and runoff as well as water balance and water use efficiencies under different irrigation methods, are present in Table 4. Both drainage and runoff were negligible under the drip irrigation method, and rainfed conditions and border irrigation had higher average drainage by 22.2 mm, and evaporation rates were decreased in the rainfed and drip fields. The rainfed and drip plots had an average evaporation rate of 139.5 and 153.5 mm, respectively. While border irrigation method exhibited 173.4 mm evaporation, which was 13% increased evaporation as compared to drip plot. In the case of crop transpiration, the drip and rainfed plots had 6–24% decreased water loss through transpiration than border irrigation method. The border irrigated plot had 209.4 mm transpiration. Minimum evapotranspiration from the drip plot indicates increased crop water use [50]. Furthermore, the drip field increased water balance by 46% as compared to border irrigation method. The border irrigation method showed an average 58.25 mm water balance. The lower water balance was under rainfed condition with 48.38 mm. Due to minimum water loss, the WUE of drip irrigated maize was 3.23 kg m–3, significantly higher than border irrigation method and rainfed maize. In the first season, the WUE under drip irrigation method was 3.41 kg m–3, and the subsequent year was 3.15 kg m–3. Thus, the order of WUE was drip irrigation > border irrigation > rainfed condition.
Water use | Ire (mm) |
Pr (mm) |
Eva (mm) |
Tr (mm) |
D (mm) |
R (mm) |
Wbal (mm) | WUE (kg m–3) |
|
2018 | Border | 150 | 370 | 170.7 | 216.3 | 24.5 | 13.6 | 75.2 | 2.61 |
Drip | 90 | 370 | 152.0 | 195.0 | 2.0 | 0.5 | 110.5 | 3.31 | |
Rainfed | 0.0 | 370 | 137.0 | 160.0 | 1.0 | 0.6 | 71.3 | 2.31 | |
2019 | Border | 150 | 326 | 176.1 | 202.5 | 19.9 | 11.5 | 41.3 | 2.72 |
Drip | 90 | 326 | 155.0 | 197.0 | 3.0 | 1.0 | 60.0 | 3.15 | |
Rainfed | 0.0 | 326 | 142.0 | 157.0 | 1.0 | 0.6 | 25.3 | 2.19 | |
Note: Ire is Irrigation amount; Pr is precipitation; Eva is evapotration; Tr is transpiration rate; D is drainage; R is runoff; Wbal is water balance and WUE is water use efficiency. |
The summary of nitrogen fate under different irrigation methods is given in Table 5. Results showed that the N leaching was negligible under drip irrigation method and rainfed condition and higher in border method. The border irrigation treatment had an average N leaching of 21.8 kg ha–1. Whereas the ammonia volatilization was higher in both drip and rainfed plots and lowered in border irrigation method. The average ammonia volatilization was 19.5 kg ha–1 in drip irrigation method and 30 kg ha–1 in rainfed condition, significantly higher than border irrigation method. The N denitrification rates were minor in all treatments. The maximum N leaching under border irrigation method could be attributed to excessive application of water, and the reason of increased ammonia volatilization from rainfed and drip fields could be due to low moisture at the soil surface layer. In the case of N mineralization and crop N uptake. The border irrigated maize had an average 54.1 kg N ha–1 mineralization, and drip plot had 41 kg N mineralization ha–1. Despite this difference in N mineralization, the drip irrigated maize had almost equal N utilization. The average crop N uptake of the drip irrigated crop was 191 kg ha–1, which was only 1.2% lower than border irrigated maize and 4% higher than rainfed maize. Moreover, the drip irrigated maize field had significantly good N balance as compared to border irrigated and rainfed maize. The drip field had 27.2 kg ha–1 N balance, 51.3–63.7% higher than border and rainfed fields. Result showed that the drip irrigated maize also had higher FNUE. The FNUE of drip irrigated maize increased by 24.1%. Thus, the order of FNUE was drip irrigation > border irrigation > rainfed condition.
Nitrogen fate | Nfer (kg ha–1) |
Nmin (kg ha–1) |
Nvol (kg ha–1) |
Nden (kg ha–1) |
Nlea (kg ha–1) |
Nup (kg ha–1) |
Nbal (kg ha–1) |
FNUE (kg kg–1) |
|
2018 | Border | 200 | 56.9 | 13.5 | 8.4 | 21.6 | 196.5 | 16.7 | 41.71 |
Drip | 200 | 43.0 | 18.0 | 2.0 | 0.7 | 190.0 | 32.3 | 54.51 | |
Rainfed | 200 | 36.0 | 31.0 | 3.0 | 0.8 | 187.0 | 14.2 | 30.98 | |
2019 | Border | 200 | 53.3 | 13.0 | 9.0 | 22.0 | 190.0 | 19.3 | 43.43 |
Drip | 200 | 39.0 | 21.0 | 3.0 | 0.8 | 192.0 | 22.2 | 51.14 | |
Rainfed | 200 | 33.0 | 29.0 | 5.0 | 0.9 | 179.0 | 19.1 | 30.61 | |
Note: Nfer is fertilization; Nmin is net mineralization; Nvol is ammonia volatilization; Nden is denitrification; Nlea is nitrogen leaching; Nup is crop Nitrogen uptake, Nbal is nitrogen balance and FNUE is Nitrogen use efficiency [51]. |
To evaluate the impact of climate change on spring maize yield, water, and nitrogen fate, we ran the model for the long-term period from 2000–2017. We selected conventional border irrigation method as a treatment. The model well-captured effect of dry and rainy seasons. The results showed that the rainfall anomaly significantly influences spring maize yield, water, and N balances. The yield bars were fluctuated high-low according to the rainfall variability. The lower the rate of rainfall, the lower the spring maize yield [46]. The dry season produced 30.8% low yield as compared to rainy season. The average yield for rainy years simulated as 10,386 kg ha–1 and for dry years as 7188 kg ha–1. Among the rainy years, maximum yields were 12,147 and 12,360 kg ha–1 in 2004 and 2006, respectively, whereas simulated yields in 2013 and 2014 rainy years were lowest. The reason of low yield during these seasons was rainfall uncertainty. Over 70% of rainfall occurred after 60 days of sowing. Consequently, both WUE and FNUE also affected. In dry and rainy seasons, the average WUEs were 2.06 and 2.49 kg m–3, and FNUEs were 30.71 and 32.97 kg kg–1, respectively. Furthermore, we found that N leaching was highly correlated with drainage amounts. Increased and decreased amount of drainage resulted in increased and decreased N leaching. The maximum average drainage was 50.6 mm, which resulted in 27.2 kg ha–1 N leaching. The year–wise simulated results are shown in Figure 6.
Crop water productivity is strongly related to irrigation methods [52]. Farmers' based border irrigation methods tend to induce excessive soil water content in the topsoil layer, which possibly caused water-stress and subsequently affected grain yield [53]. The soil of the study area is predominantly sandy clay loam. It has hard surface layer due to low organic content [54]. Therefore, in current study, the border irrigation method showed 13% higher evaporation rate as compared to drip irrigation method. The drip irrigated maize not only reduced evaporation, but also maintained plant transpiration, thus increased yield by 14% with 40% water conservation. Compared to rainfed maize, the drip irrigated maize highly increased yield. The increase of yield under drip irrigation method, could be attributed to improved WUE. The results of yield increase, water saving and enhanced WUE of drip irrigated maize, agree with the findings of Zhou et al. [55] study, which showed that keeping soil moisture at an optimal level via reducing irrigation amounts and increasing number of irrigations, maximizes crop yield as compared to increasing soil water content by apply full irrigation [56]. In our study, drip irrigation primarily promoted growth of roots, biomass accumulation and then grain–filling. Figure 5 shows the effects of border irrigation method, drip and rainfed conditions on spring maize. Previously, also demonstrated by Sandhu et al. [57] and Yang et al. [18], they found that drip irrigation method improved maize yield by > 13% while reducing ET by 20 mm with significantly reduction in runoff rates and deep–drainage water losses. Similarly, Leghari et al. [21] compared yield and water-saving of spring maize under different irrigation methods and stated that drip irrigation increased yield by 9.8% at 62.5% reduced water consumption, improved WUE by 24.5% than border irrigation method. The increase of WUE was mainly due to a significant decrease in water loss. Overall, results suggesting that the drip irrigation method would ensure food security by increasing yield and maintain ecosystem.
The loss of N due to volatilization, denitrification, and leaching is not only yield affecting factor, but it is also the main cause of maximal N deposition in terrestrial and aquatic systems. Because N in the form of NO3− is highly water–soluble, it moves freely with excessive water [58]. Plants can uptake half of the applied N and remaining resulting in environmental contamination. Wang et al. [22] reported that border irrigation method resulted in 97 to 186 kg N ha−1 leaching, while drip irrigation method reduced N leaching by 68%. The N emission losses including ammonia volatilization and denitrification were also higher under border irrigation method. Similarly, our study showed that border irrigation method is inefficient as compared to drip irrigation method, which is efficient due to high N productivity. The drip irrigated maize produced maximum yield while reducing N leaching loss by 97%, improving N utilization, and enhancing FNUE by 24%. These findings are in line with those of Li et al. [59] who found less than 40 kg and more than 50 kg ha–1 N leaching when irrigation application was less and greater than 300 mm ha–1. Reduced N losses in drip-irrigated maize could be due to direct N supplementation to the plant root zone in the required amount. Drip irrigation enhanced soil moisture first, then N uptake and thus improved nutrient management [60].
Maize is a drought-tolerance crop (ABA accumulates in leaf tissue, which induces stomatal closure). However, prolonged drought stress could adversely affect vegetative growth, biomass accumulation, formation of reproductive organs and yield attributes [48–61]. Deficit water-stress mainly reduces plant height due to decrease in photosynthesis and osmotic potential [62]. Scenario analysis on the long-term impact of climate change exhibited that rainfall distribution during the growing season is critical. For example, simulated yields in 2012 and 2017 seasons were 3943 and 5018 kg ha–1, respectively, lower than the rest of growing seasons due to rainfall irregularity from sowing to harvesting. On average, the dry years reduced yield by 30.8% as compared to rainy season. Similar yield response of spring maize according to the season to season variation in local weather conditions, also observed by Shi et al. [46] who found that the crop yields were decreased as rainfall decreased in 2011 and 2012. During the dry season, the supplemental application of water is recommended at the critical growth stage or when the crop is profoundly suffering from water stress and finally wilting up [63].
In this study, we simulated effect of border irrigation method, drip irrigation and rainfed conditions under farmer-based fertilizer application rate. The model was evaluated for crop growth and soil hydraulic parameters, except for N leaching and drainage, as well as evaporation and transpiration due to the unavailability of measured data. In future works, the model should be assessed for these factors under different irrigation methods. Furthermore, the exhibited results of model to the dry and rainy years, should be examined under long term field experiments for better understanding of the impact of dry and rainy seasons on crop yield, water, and N balances.
The statistical indicators exhibited that the WHCNS model could be potentially used in the region to analyze the water movement and N cycling quantitatively. Based on the measured and modelled results, we found that the drip irrigation method would significantly improve crop yield as compared to conventional border irrigation method and rainfed. The drip irrigation not only increased grain yield but also conserved water by decreasing evaporation, drainage, and runoff. Moreover, the drip irrigation had a negligible amount of N leaching. Both WUE and FNUEs were higher in drip irrigated maize. Whereas border irrigation method resulted in increased evaporation as compared to drip irrigation, and this was a major pathway of water loss under the conventional border irrigation method. Despite drip irrigation system is costly, the water saving, and yield increase under drip irrigation suggesting suitability of this system in the NCP. Furthermore, considering rainfall variability based on long–term dry and rainy seasons, the scenario analysis showed that although maize is a drought tolerance crop, the rainfall uncertainty during the crop growth period could result in low yield.
The study was funded by China Agriculture Research System of MOF and MARA (CARS-02), National Key Research and Development Program of China (2017YFD0301102–5) and Leading Talent Project of Science and Technology Innovation of Central Plains (194200510008).
The authors have declared that no competing interests exist.
[1] | Z. Liu, A. Qin, B. Zhao, S. T. Ata-Ul-Karim, J. Xiao, J. Sun, et al., Yield response of spring maize to inter-row subsoiling and soil water deficit in Northern China, PLoS One, 11 (2016). |
[2] |
W. Feng, M. Zhong, J. M. Lemoine, R. Biancale, H. T. Hsu, J. Xia, Evaluation of groundwater depletion in North China using the Gravity Recovery and Climate Experiment (GRACE) data and ground-based measurements, Water Resour. Res., 49 (2013), 2110–2118. doi: 10.1002/wrcr.20192
![]() |
[3] |
X. Yang, Y. Chen, S. Pacenka, W. Gao, L. Ma, G. Wang, et al., Effect of diversified crop rotations on groundwater levels and crop water productivity in the North China Plain, J. Hydrol., 522 (2015), 428–438. doi: 10.1016/j.jhydrol.2015.01.010
![]() |
[4] |
L. Min, Y. Shen, H. Pei, Estimating groundwater recharge using deep vadose zone data under typical irrigated cropland in the piedmont region of the North China Plain, J. Hydrol., 527 (2015), 305–315. doi: 10.1016/j.jhydrol.2015.04.064
![]() |
[5] | J. Liu, G. Cao, C. Zheng, Sustainability of groundwater resources in the North China Plain, Sustaining Groundwater Res., (2011), 69–87. |
[6] | X. Yang, Y. Chen, S. Pacenka, W. Gao, M. Zhang, P. Sui, T. S. Steenhuis, Recharge and groundwater use in the north china plain for six irrigated crops for an eleven year period, PLoS One, 10 (2015), 1–17. |
[7] |
H. Gong, Y. Pan, L. Zheng, X. Li, L. Zhu, C. Zhang, et al., Long-term groundwater storage changes and land subsidence development in the North China Plain (1971–2015), Hydrogeol. J., 26 (2018), 1417–1427. doi: 10.1007/s10040-018-1768-4
![]() |
[8] | L. Wang, J. Xie, Z. Luo, Y. Niu, J. A. Coulter, R. Zhang, et al., Forage yield, water use efficiency, and soil fertility response to alfalfa growing age in the semiarid Loess Plateau of China, Agric. Water Manage., 243 (2021). |
[9] | M. N. Khan, M. Mobin, Z. K. Abbas, S. A. Alamri, Fertilizers and their contaminants in soils, surface and groundwater, Encycl. Anthropocene, 5 (2018), 225–240. |
[10] |
Z. Cui, X. Chen, F. Zhang, Current Nitrogen management status and measures to improve the intensive wheat-maize system in China, Ambio, 39 (2010), 376–384. doi: 10.1007/s13280-010-0076-6
![]() |
[11] |
J. Letey, P. Vaughan, Soil type, crop and irrigation technique affect nitrogen leaching to groundwater, Calif. Agric., 67 (2013), 231–241. doi: 10.3733/ca.E.v067n04p231
![]() |
[12] |
X. Ju, C. Zhang, Nitrogen cycling and environmental impacts in upland agricultural soils in North China: A review, J. Integr. Agric., 16 (2017), 2848–2862. doi: 10.1016/S2095-3119(17)61743-X
![]() |
[13] |
Z. Wang, J. Li, Y. Li, Simulation of nitrate leaching under varying drip system uniformities and precipitation patterns during the growing season of maize in the North China Plain, Agric. Water Manage., 142 (2014), 19–28. doi: 10.1016/j.agwat.2014.04.013
![]() |
[14] |
Ö. Çetin, E. Akalp, Efficient use of water and fertilizers in irrigated agriculture: drip irrigation and fertigation, Acta Hortic. Regiotect., 22 (2019), 97–102. doi: 10.2478/ahr-2019-0019
![]() |
[15] |
D. Zaccaria, M. T. Carrillo-Cobo, A. Montazar, D. H. Putnam, K. Bali, Assessing the viability of sub-surface drip irrigation for resource-efficient alfalfa production in central and Southern California, Water, 9 (2017), 837. doi: 10.3390/w9110837
![]() |
[16] |
M. Umair, T. Hussain, H. Jiang, A. Ahmad, J. Yao, Y. Qi, et al., Water-saving potential of subsurface drip irrigation for winter wheat, Sustainability, 11 (2019), 2978. doi: 10.3390/su11102978
![]() |
[17] |
H. M. Darouich, C. M. G. Pedras, J. M. Gonçalves, L. S. Pereira, Drip vs. surface irrigation: A comparison focussing on water saving and economic returns using multicriteria analysis applied to cotton, Biosyst. Eng., 122 (2014), 74–90. doi: 10.1016/j.biosystemseng.2014.03.010
![]() |
[18] | D. Yang, S. Li, S. Kang, T. Du, P. Guo, X. Mao, et al., Effect of drip irrigation on wheat evapotranspiration, soil evaporation and transpiration in Northwest China, Agric. Water Manage., 232 (2020). |
[19] | Y. Wei, L. Sun, S. Wang, Y. Wang, Z. Zhang, Q. Chen, et al., Effects of different irrigation methods on water distribution and nitrate nitrogen transport of cucumber in greenhouse, Nongye Gongcheng Xuebao/Trans. Chinese Soc. Agric. Eng., 26 (2010), 67–72. |
[20] |
J. E. Ayars, A. Fulton, B. Taylor, Subsurface drip irrigation in California–Here to stay?, Agric. Water Manage., 157 (2015), 39–47. doi: 10.1016/j.agwat.2015.01.001
![]() |
[21] |
S. J. Leghari, K. Hu, Y. Wei, T. C. Wang, T. A. Bhutto, M. Buriro, Modelling water consumption, N fates and maize yield under different water-saving management practices in China and Pakistan. Agric. Water Manage., 255 (2021), 107033. doi: 10.1016/j.agwat.2021.107033
![]() |
[22] | Y. Wang, S. Li, H. Liang, K. Hu, S. Qin, H. Guo, Comparison of water-and nitrogen-use efficiency over drip irrigation with border irrigation based on a model approach, Agronomy, 10 (2020), 1890. |
[23] |
Z. Zhao, X. Qin, Z. Wang, E. Wang, Performance of different cropping systems across precipitation gradient in North China Plain, Agric. For. Meteorol., 259 (2018), 162–172. doi: 10.1016/j.agrformet.2018.04.019
![]() |
[24] |
X. L. Yang, Y. Q. Chen, T. S. Steenhuis, S. Pacenka, W. S. Gao, L. Ma, et al., Mitigating groundwater depletion in North China plain with cropping system that alternate deep and shallow rooted crops, Front. Plant Sci., 8 (2017), 980. doi: 10.3389/fpls.2017.00980
![]() |
[25] | S. Lorvanleuang, Y. Zhao, Automatic irrigation system using android, OALib, 5 (2018). |
[26] |
J. Wang, P. Wang, Q. Qin, H. Wang, The effects of land subsidence and rehabilitation on soil hydraulic properties in a mining area in the Loess Plateau of China, Catena, 159 (2017), 51–59. doi: 10.1016/j.catena.2017.08.001
![]() |
[27] | A. N. Beretta, A. V. Silbermann, L. Paladino, D. Torres, D. Bassahun, R. Musselli, et al., Soil texture analyses using a hydrometer: modification of the Bouyoucos method, Cienc. Invest. Agrar., 41 (2014), 263–271. |
[28] |
M. T. van Genuchten, A closed-form equation for predicting the hydraulic conductivity of unsaturated soils, Soil Sci. Soc. Am. J., 44 (1980), 892–898. doi: 10.2136/sssaj1980.03615995004400050002x
![]() |
[29] |
S. Zhang, Y. Wei, M. N. Kandhro, F. Wu, Evaluation of FDR MI2X and new WiTu technology sensors to determine soil water content in the corn and weed field, AIMS Agric. Food, 5 (2020), 169–180. doi: 10.3934/agrfood.2020.1.169
![]() |
[30] |
S. J. Leghari, K. Hu, H. Liang, Y. Wei, Modeling water and nitrogen balance of different cropping systems in the North China Plain, Agronomy, 9 (2019), 696. doi: 10.3390/agronomy9110696
![]() |
[31] | J. Šimůnek, A. M. Šejna, H. Saito, M. Sakai, M. T. V. Genuchten, The HYDRUS-1D software package for simulating the movement of water, heat, and multiple solutes in variably saturated media, Software Ser., 2013. |
[32] |
J. D. Hanson, L. R. Ahuja, M. D. Shaffer, K. W. Rojas, D. G. DeCoursey, H. Farahani, et al., RZWQM: Simulating the effects of management on water quality and crop production, Agric. Syst., 57 (1998), 161–195. doi: 10.1016/S0308-521X(98)00002-X
![]() |
[33] | R. G. Allen, L. S. Pereira, D. Raes, M. Smith, Crop evapotranspiration: Guidelines for computing crop requirements, Irrig. Drain. Pap., 56 (1980) |
[34] | A. Y. Hachum, J. F. Alfaro, Rain infiltration into layered soils: prediction, J. Irrig. Drain. Div. ASCE., 106 (1980). |
[35] |
Y. Pachepsky, D. Timlin, W. Rawls, Generalized Richards' equation to simulate water transport in unsaturated soils, J. Hydrol., 272 (2003), 3–13. doi: 10.1016/S0022-1694(02)00251-2
![]() |
[36] |
J. W. Jung, K. S. Yoon, D. H. Choi, S. S. Lim, W. J. Choi, S. M. Choi, et al., Water management practices and SCS curve numbers of paddy fields equipped with surface drainage pipes, Agric. Water Manage., 110 (2012), 78–83. doi: 10.1016/j.agwat.2012.03.014
![]() |
[37] | S. Hansen, P. Abrahamsen, C. T. Petersen, M. Styczen, Daisy: model use, calibration, and validation, 55 (2012), 1315–1333. |
[38] | H. Liang, K. Hu, W. D. Batchelor, Z. Qi, B. Li, An integrated soil-crop system model for water and nitrogen management in North China, Sci. Rep., 6 (2016). |
[39] |
D. N. Moriasi, J. G. Arnold, M. W. V. Liew, R. L. Bingner, R. D. Harmel, T. L. Veith, Model evaluation guidelines for systematic quantification of accuracy in watershed simulations. Trans., Asabe., 50 (2007), 885–900. doi: 10.13031/2013.23153
![]() |
[40] |
C. J. Willmott, On the validation of models, Phys. Geogr., 2 (1981), 184–194. doi: 10.1080/02723646.1981.10642213
![]() |
[41] | M. M. Mukaka, Statistics corner: A guide to appropriate use of correlation coefficient in medical research, Malawi Med. J., 24 (2012), 69–71. |
[42] | M. Razaq, P. Zhang, H. L. Shen, Salahuddin, Influence of nitrogen and phosphorous on the growth and root morphology of Acer mono, PLoS One, 12 (2017). |
[43] |
A. A. Siyal, A. S. Mashori, K. L. Bristow, M. T. V. Genuchten, Alternate furrow irrigation can radically improve water productivity of okra, Agric. Water Manage., 173 (2016), 55–60. doi: 10.1016/j.agwat.2016.04.026
![]() |
[44] | G. Murray-Tortarolo, V. J. Jaramillo, M. Maass, P. Friedlingstein, S. Sitch, The decreasing range between dry- and wet-season precipitation over land and its effect on vegetation primary productivity, PLoS One, 12 (2017). |
[45] |
X. Gao, X. Chen, T. W. Biggs, H. Yao, Separating wet and dry years to improve calibration of SWAT in Barrett watershed, Southern California, Water, 10 (2018), 274. doi: 10.3390/w10030274
![]() |
[46] | X. Shi, W. D. Batchelor, H. Liang, S. Li, B. Li, K. Hu, Determining optimal water and nitrogen management under different initial soil mineral nitrogen levels in northwest China based on a model approach, Agric. Water Manage., 234 (2020). |
[47] |
Y. He, H. Liang, K. Hu, H. Wang, L. Hou, Modeling nitrogen leaching in a spring maize system under changing climate and genotype scenarios in arid Inner Mongolia, China, Agric. Water Manage., 210 (2018), 316–323. doi: 10.1016/j.agwat.2018.08.017
![]() |
[48] |
H. Liang, K. Hu, W. Qin, Q. Zuo, L. Guo, Y. Tao, et al., Ground cover rice production system reduces water consumption and nitrogen loss and increases water and nitrogen use efficiencies, F. Crop. Res., 233 (2019), 70–79. doi: 10.1016/j.fcr.2019.01.003
![]() |
[49] | Q. Xu, K. Hu, H. Liang, S. J. Leghari, M. T. Knudsen, Incorporating the WHCNS model to assess water and nitrogen footprint of alternative cropping systems for grain production in the North China Plain, J. Clean. Prod., 263 (2020). |
[50] |
Y. Liang, Y. Jiang, M. Du, B. Li, L. Chen, M. Chen, et al., ZmASR3 from the maize ASR gene family positively regulates drought tolerance in transgenic arabidopsis, Int. J. Mol. Sci., 20 (2019), 2278. doi: 10.3390/ijms20092278
![]() |
[51] |
H. Liang, K. Hu, W. Qin, Q. Zuo, Y. Zhang, Modelling the effect of mulching on soil heat transfer, water movement and crop growth for ground cover rice production system, F. Crop. Res., 201 (2017), 97–107. doi: 10.1016/j.fcr.2016.11.003
![]() |
[52] |
S. K. Jha, T. S. Ramatshaba, G. Wang, Y. Liang, H. Liu, Y. Gao, et al., Response of growth, yield and water use efficiency of winter wheat to different irrigation methods and scheduling in North China Plain, Agric. Water Manage., 217 (2019), 292–302. doi: 10.1016/j.agwat.2019.03.011
![]() |
[53] | S. Kang, L. Zhang, Y. Liang, X. Hu, H. Cai, B. Gu, Effects of limited irrigation on yield and water use efficiency of winter wheat in the Loess Plateau of China, Agric. Water Manage., 55 (2002), 203–216. |
[54] |
G. Wang, Q. Fang, B. Wu, H. Yang, Z. Xu, Relationship between soil erodibility and modeled infiltration rate in different soils, J. Hydrol., 528 (2015), 408–418. doi: 10.1016/j.jhydrol.2015.06.044
![]() |
[55] |
B. Zhou, X. Sun, Z. Ding, W. Ma, M. Zhao, Multisplit nitrogen application via drip irrigation improves maize grain yield and nitrogen use efficiency, Crop Sci., 57 (2017), 1687–1703. doi: 10.2135/cropsci2016.07.0623
![]() |
[56] |
S. E. El-Hendawy, E. M. Hokam, U. Schmidhalter, Drip irrigation frequency: The effects and their interaction with nitrogen fertilization on sandy soil water distribution, maize yield and water use efficiency under Egyptian conditions, J. Agron. Crop Sci., 194 (2008), 180–192. doi: 10.1111/j.1439-037X.2008.00304.x
![]() |
[57] |
O. S. Sandhu, R. K. Gupta, H. S. Thind, M. L. Jat, H. S. Sidhu, Yadvinder-Singh, Drip irrigation and nitrogen management for improving crop yields, nitrogen use efficiency and water productivity of maize-wheat system on permanent beds in north-west India, Agric. Water Manage., 219 (2019), 19–26. doi: 10.1016/j.agwat.2019.03.040
![]() |
[58] | S. J. Leghari, N. A. Wahocho, G. M. Laghari, A. HafeezLaghari, G. MustafaBhabhan, K. H. Talpur, et al., Role of nitrogen for plant growth and development : A review, Adv. Environ. Biol., 10 (2016), 209–218. |
[59] |
Y. Li, H. Liu, G. Huang, R. Zhang, H. Yang, Nitrate nitrogen accumulation and leaching pattern at a winter wheat: summer maize cropping field in the North China Plain, Environ. Earth Sci., 75 (2016), 118. doi: 10.1007/s12665-015-4867-8
![]() |
[60] |
D. Wu, X. Xu, Y. Chen, H. Shao, E. Sokolowski, G. Mi, Effect of different drip fertigation methods on maize yield, nutrient and water productivity in two-soils in Northeast China, Agric. Water Manage., 213 (2019), 200–211. doi: 10.1016/j.agwat.2018.10.018
![]() |
[61] |
A. Gadédjisso-Tossou, T. Avellán, N. Schütze, Potential of deficit and supplemental irrigation under climate variability in northern Togo, West Africa, Water, 10 (2018), 1803.. doi: 10.3390/w10121803
![]() |
[62] |
W. Zhao, L. Liu, Q. Shen, J. Yang, X. Han, F. Tian, et al., Effects of Water Stress on Photosynthesis, Yield, and Water Use Efficiency in Winter Wheat, Water, 12 (2020), 2127. doi: 10.3390/w12082127
![]() |
[63] |
Y. Liu, H. Yang, J. Li, Y. Li, H. Yan, Estimation of irrigation requirements for drip-irrigated maize in a sub-humid climate, J. Integr. Agric., 17 (2018), 677–692. doi: 10.1016/S2095-3119(17)61833-1
![]() |
1. | Biao Jia, Jiangpeng Fu, Huifang Liu, Zhengzhou Li, Yu Lan, Xue Wei, Yongquan Zhai, Bingyuan Yun, Jianzhen Ma, Hao Zhang, Estimation of Critical Nitrogen Concentration Based on Leaf Dry Matter in Drip Irrigation Spring Maize Production in Northern China, 2022, 14, 2071-1050, 9838, 10.3390/su14169838 | |
2. | Meenakshi Rawat, Rintu Sen, Ikenna Onyekwelu, Travis Wiederstein, Vaishali Sharda, Modeling of Groundwater Nitrate Contamination Due to Agricultural Activities—A Systematic Review, 2022, 14, 2073-4441, 4008, 10.3390/w14244008 | |
3. | Mengjie Liu, Fei Liang, Quansheng Li, Guodong Wang, Yuxin Tian, Hongtao Jia, Enhancement growth, water use efficiency and economic benefit for maize by drip irrigation in Northwest China, 2023, 13, 2045-2322, 10.1038/s41598-023-35611-9 | |
4. | Xiaomeng Huang, Xinpeng Xu, Qichao Zhu, Yitao Zhang, Optimizing water and nitrogen inputs for sustainable wheat yields and minimal environmental impacts, 2024, 220, 0308521X, 104061, 10.1016/j.agsy.2024.104061 | |
5. | Kai Liu, Huan Liao, Haibo Hao, Zhenan Hou, Water and nitrogen supply at spatially distinct locations improves cotton water productivity and nitrogen use efficiency and yield under drip irrigation, 2024, 296, 03783774, 108808, 10.1016/j.agwat.2024.108808 | |
6. | Hui Guo, Sien Li, A Review of Drip Irrigation’s Effect on Water, Carbon Fluxes, and Crop Growth in Farmland, 2024, 16, 2073-4441, 2206, 10.3390/w16152206 | |
7. | Qing Sun, Hongxiang Zhang, Xuejie Li, Zixuan Zhao, Zengxu Li, Peiyu Zhang, Shutang Liu, Wen Jiang, Xuefang Sun, Effects of Drip Tape Layout and Flow Rate on Water and Nitrogen Distributions within the Root Zone and Summer Maize Yield in Sandy Tidal Soil, 2023, 13, 2073-4395, 2689, 10.3390/agronomy13112689 | |
8. | Yu Cheng, Tibin Zhang, Xinglu Hu, Zhenyuan Liu, Qing Liang, Sihui Yan, Hao Feng, Kadambot H.M. Siddique, Drip fertigation triggered by soil matric potential reduces residual soil nitrate content and improves maize nitrogen uptake and yield stability in an arid area, 2023, 150, 11610301, 126932, 10.1016/j.eja.2023.126932 | |
9. | Shah Jahan Leghari, Kelin Hu, Yichang Wei, Tongchao Wang, Yaseen Laghari, Modelling the effects of cropping systems and irrigation methods on water consumption, N fates and crop yields in the North China Plain, 2024, 218, 01681699, 108677, 10.1016/j.compag.2024.108677 | |
10. | Danyue Zhao, Wei Wang, Xiangxiang Ji, Pute Wu, La Zhuo, Nitrogen cycling and associated grey water footprint in croplands under different irrigation practices, 2024, 480, 09596526, 144081, 10.1016/j.jclepro.2024.144081 | |
11. | Jin Tong, Yulei Xiong, Yu Lu, Wen Li, Wen Lin, Jianfu Xue, Min Sun, Yuechao Wang, Zhiqiang Gao, Drip Fertigation with Moderate Nitrogen Topdressing Rate Achieves High Nitrogen and Water Use Efficiencies for Irrigated Wheat, 2025, 15, 2073-4395, 259, 10.3390/agronomy15020259 | |
12. | Yuyan Fan, Wen Li, Limin Zhang, Jinxiao Song, Depeng Wang, Jianfu Xue, Yuechao Wang, Zhiqiang Gao, Wheat Nitrogen Use and Grain Protein Characteristics Under No-Tillage: A Greater Response to Drip Fertigation Compared to Intensive Tillage, 2025, 15, 2073-4395, 588, 10.3390/agronomy15030588 | |
13. | Zhendong Liu, Mingliang Gao, Fuying Liu, Shiyu Wei, Haoze Zhang, Huabin Yuan, Quanqi Li, Rui Zong, Evapotranspiration partitioning and crop water productivity under varying planting densities of maize in north China, 2025, 22145141, 10.1016/j.cj.2025.03.005 | |
14. | Khushkarandeep Singh, Vinay Kumar Sindhu, Ajmer Singh Brar, J.S. Lamba, J.S. Deol, R.S. Jakhar, Ranjeet Singh, Optimizing water use in spring maize through integrated agro-techniques in groundwater-stressed regions of northwestern India, 2025, 317, 03783774, 109635, 10.1016/j.agwat.2025.109635 |
Depth (cm) |
BD (g cm–3) |
Particle fraction (%) | Soil texture (USDA) |
pH | θs (cm3 cm–3) |
θfc (cm3 cm–3) |
θwp (cm3 cm–3) |
Ks (cm d–1) |
||
Sand | Silt | Clay | ||||||||
0–20 | 1.31 | 52 | 23 | 25 | Sandy clay loam | 8.2 | 0.42 | 0.31 | 0.14 | 42.2 |
20–40 | 1.34 | 54 | 25 | 21 | Sandy clay loam | 8.0 | 0.40 | 0.32 | 0.15 | 33.7 |
40–60 | 1.30 | 57 | 23 | 20 | Sandy clay loam | 7.9 | 0.41 | 0.31 | 0.13 | 31.4 |
60–80 | 1.35 | 58 | 22 | 20 | Sandy clay loam | 7.5 | 0.40 | 0.28 | 0.14 | 50.5 |
80–100 | 1.37 | 57 | 22 | 21 | Sandy clay loam | 7.8 | 0.39 | 0.27 | 0.14 | 51.2 |
Note: BD; Bulk density, θs; Soil saturated water content, θfc; Field capacity, θwp; Permanent wilting point and Ks; Saturated hydraulic conductivity. |
Parameters and definitions of parameters | Values |
Tbase defined base temperature for crop growth | 8 ℃ |
Tsum defined as effective accumulated temperature | 2249℃ |
Ke define as extinction coefficient (0–1) | 0.7 |
Kc_ini defined crop coefficients at initial stage (0–1) | 0.64 |
Kc_mid defined as crop coefficients at middle stage (0–1) | 1.26 |
Kc_end defined as crop coefficients at final stage (0–1) | 0.7 |
SLA_max defined as maximum specific leaf area (m2 kg–1) | 29 |
SLA_min defined as minimum specific leaf area (m2 kg–1) | 20 |
AMAX defined as maximum assimilation rate (kg ha–1) | 50 |
R_max defined as maximum soil depth contributing to N transformation (cm) | 142 |
Vn defined as maximum nitrification rate (mg L–1 d–1) | 100 |
Kn defined as half saturation constant for nitrification (g m−3) | 200 |
Kd defined as empirical proportionality factor for denitrification (−) | 0.2 |
Ad defined as first–order kinetic coefficient for denitrification (d−1) | 0.05 |
Kv defined as first–order kinetic coefficient for ammonia volatilization denitrification (d−1) | 0.03 |
SWC (cm3 cm–3) | Border | Drip | Rainfed | ||||||
RRMSE | IA | r | RRMSE | IA | r | RRMSE | IA | r | |
0–20 cm | 24.0 | 0.76 | 0.70 | 30.0 | 0.62 | 0.58 | 24.0 | 0.80 | 0.73 |
20–40 cm | 21.0 | 0.62 | 0.47 | 9.0 | 0.77 | 0.70 | 5.0 | 0.91 | 0.84 |
40–60 cm | 26.0 | 0.58 | 0.43 | 35.0 | 0.53 | 0.47 | 14.0 | 0.88 | 0.80 |
60–80 cm | 27.0 | 0.64 | 0.62 | 36.0 | 0.56 | 0.46 | 21.0 | 0.77 | 0.65 |
80–100 cm | 21.0 | 0.61 | 0.69 | 14.0 | 0.73 | 0.58 | 12.0 | 0.82 | 0.78 |
SNC (mg kg–1) | |||||||||
0–20 cm | 27.0 | 0.84 | 0.76 | 22.0 | 0.93 | 0.89 | 36.0 | 0.76 | 0.66 |
20–40 cm | 14.0 | 0.92 | 0.93 | 21.0 | 0.75 | 0.61 | 24.0 | 0.80 | 0.67 |
40–60 cm | 22.0 | 0.86 | 0.77 | 26.0 | 0.81 | 0.77 | 24.0 | 0.88 | 0.81 |
60–80 cm | 21.0 | 0.90 | 0.96 | 27.0 | 0.85 | 0.87 | 33.0 | 0.87 | 0.84 |
80–100 cm | 18.0 | 0.96 | 0.94 | 38.0 | 0.83 | 0.79 | 28.0 | 0.94 | 0.89 |
Biomass (kg ha–1) | 19.0 | 0.98 | 0.98 | 20.0 | 0.98 | 0.97 | 24.0 | 0.98 | 0.97 |
Grain yield (kg ha–1) | 1.0 | 0.64 | 1.00 | 2.0 | 0.61 | 1.00 | 2.0 | 0.55 | 1.00 |
Note: SWC: Soil water content; SNC: Soil nitrate concentrations. |
Water use | Ire (mm) |
Pr (mm) |
Eva (mm) |
Tr (mm) |
D (mm) |
R (mm) |
Wbal (mm) | WUE (kg m–3) |
|
2018 | Border | 150 | 370 | 170.7 | 216.3 | 24.5 | 13.6 | 75.2 | 2.61 |
Drip | 90 | 370 | 152.0 | 195.0 | 2.0 | 0.5 | 110.5 | 3.31 | |
Rainfed | 0.0 | 370 | 137.0 | 160.0 | 1.0 | 0.6 | 71.3 | 2.31 | |
2019 | Border | 150 | 326 | 176.1 | 202.5 | 19.9 | 11.5 | 41.3 | 2.72 |
Drip | 90 | 326 | 155.0 | 197.0 | 3.0 | 1.0 | 60.0 | 3.15 | |
Rainfed | 0.0 | 326 | 142.0 | 157.0 | 1.0 | 0.6 | 25.3 | 2.19 | |
Note: Ire is Irrigation amount; Pr is precipitation; Eva is evapotration; Tr is transpiration rate; D is drainage; R is runoff; Wbal is water balance and WUE is water use efficiency. |
Nitrogen fate | Nfer (kg ha–1) |
Nmin (kg ha–1) |
Nvol (kg ha–1) |
Nden (kg ha–1) |
Nlea (kg ha–1) |
Nup (kg ha–1) |
Nbal (kg ha–1) |
FNUE (kg kg–1) |
|
2018 | Border | 200 | 56.9 | 13.5 | 8.4 | 21.6 | 196.5 | 16.7 | 41.71 |
Drip | 200 | 43.0 | 18.0 | 2.0 | 0.7 | 190.0 | 32.3 | 54.51 | |
Rainfed | 200 | 36.0 | 31.0 | 3.0 | 0.8 | 187.0 | 14.2 | 30.98 | |
2019 | Border | 200 | 53.3 | 13.0 | 9.0 | 22.0 | 190.0 | 19.3 | 43.43 |
Drip | 200 | 39.0 | 21.0 | 3.0 | 0.8 | 192.0 | 22.2 | 51.14 | |
Rainfed | 200 | 33.0 | 29.0 | 5.0 | 0.9 | 179.0 | 19.1 | 30.61 | |
Note: Nfer is fertilization; Nmin is net mineralization; Nvol is ammonia volatilization; Nden is denitrification; Nlea is nitrogen leaching; Nup is crop Nitrogen uptake, Nbal is nitrogen balance and FNUE is Nitrogen use efficiency [51]. |
Depth (cm) |
BD (g cm–3) |
Particle fraction (%) | Soil texture (USDA) |
pH | θs (cm3 cm–3) |
θfc (cm3 cm–3) |
θwp (cm3 cm–3) |
Ks (cm d–1) |
||
Sand | Silt | Clay | ||||||||
0–20 | 1.31 | 52 | 23 | 25 | Sandy clay loam | 8.2 | 0.42 | 0.31 | 0.14 | 42.2 |
20–40 | 1.34 | 54 | 25 | 21 | Sandy clay loam | 8.0 | 0.40 | 0.32 | 0.15 | 33.7 |
40–60 | 1.30 | 57 | 23 | 20 | Sandy clay loam | 7.9 | 0.41 | 0.31 | 0.13 | 31.4 |
60–80 | 1.35 | 58 | 22 | 20 | Sandy clay loam | 7.5 | 0.40 | 0.28 | 0.14 | 50.5 |
80–100 | 1.37 | 57 | 22 | 21 | Sandy clay loam | 7.8 | 0.39 | 0.27 | 0.14 | 51.2 |
Note: BD; Bulk density, θs; Soil saturated water content, θfc; Field capacity, θwp; Permanent wilting point and Ks; Saturated hydraulic conductivity. |
Parameters and definitions of parameters | Values |
Tbase defined base temperature for crop growth | 8 ℃ |
Tsum defined as effective accumulated temperature | 2249℃ |
Ke define as extinction coefficient (0–1) | 0.7 |
Kc_ini defined crop coefficients at initial stage (0–1) | 0.64 |
Kc_mid defined as crop coefficients at middle stage (0–1) | 1.26 |
Kc_end defined as crop coefficients at final stage (0–1) | 0.7 |
SLA_max defined as maximum specific leaf area (m2 kg–1) | 29 |
SLA_min defined as minimum specific leaf area (m2 kg–1) | 20 |
AMAX defined as maximum assimilation rate (kg ha–1) | 50 |
R_max defined as maximum soil depth contributing to N transformation (cm) | 142 |
Vn defined as maximum nitrification rate (mg L–1 d–1) | 100 |
Kn defined as half saturation constant for nitrification (g m−3) | 200 |
Kd defined as empirical proportionality factor for denitrification (−) | 0.2 |
Ad defined as first–order kinetic coefficient for denitrification (d−1) | 0.05 |
Kv defined as first–order kinetic coefficient for ammonia volatilization denitrification (d−1) | 0.03 |
SWC (cm3 cm–3) | Border | Drip | Rainfed | ||||||
RRMSE | IA | r | RRMSE | IA | r | RRMSE | IA | r | |
0–20 cm | 24.0 | 0.76 | 0.70 | 30.0 | 0.62 | 0.58 | 24.0 | 0.80 | 0.73 |
20–40 cm | 21.0 | 0.62 | 0.47 | 9.0 | 0.77 | 0.70 | 5.0 | 0.91 | 0.84 |
40–60 cm | 26.0 | 0.58 | 0.43 | 35.0 | 0.53 | 0.47 | 14.0 | 0.88 | 0.80 |
60–80 cm | 27.0 | 0.64 | 0.62 | 36.0 | 0.56 | 0.46 | 21.0 | 0.77 | 0.65 |
80–100 cm | 21.0 | 0.61 | 0.69 | 14.0 | 0.73 | 0.58 | 12.0 | 0.82 | 0.78 |
SNC (mg kg–1) | |||||||||
0–20 cm | 27.0 | 0.84 | 0.76 | 22.0 | 0.93 | 0.89 | 36.0 | 0.76 | 0.66 |
20–40 cm | 14.0 | 0.92 | 0.93 | 21.0 | 0.75 | 0.61 | 24.0 | 0.80 | 0.67 |
40–60 cm | 22.0 | 0.86 | 0.77 | 26.0 | 0.81 | 0.77 | 24.0 | 0.88 | 0.81 |
60–80 cm | 21.0 | 0.90 | 0.96 | 27.0 | 0.85 | 0.87 | 33.0 | 0.87 | 0.84 |
80–100 cm | 18.0 | 0.96 | 0.94 | 38.0 | 0.83 | 0.79 | 28.0 | 0.94 | 0.89 |
Biomass (kg ha–1) | 19.0 | 0.98 | 0.98 | 20.0 | 0.98 | 0.97 | 24.0 | 0.98 | 0.97 |
Grain yield (kg ha–1) | 1.0 | 0.64 | 1.00 | 2.0 | 0.61 | 1.00 | 2.0 | 0.55 | 1.00 |
Note: SWC: Soil water content; SNC: Soil nitrate concentrations. |
Water use | Ire (mm) |
Pr (mm) |
Eva (mm) |
Tr (mm) |
D (mm) |
R (mm) |
Wbal (mm) | WUE (kg m–3) |
|
2018 | Border | 150 | 370 | 170.7 | 216.3 | 24.5 | 13.6 | 75.2 | 2.61 |
Drip | 90 | 370 | 152.0 | 195.0 | 2.0 | 0.5 | 110.5 | 3.31 | |
Rainfed | 0.0 | 370 | 137.0 | 160.0 | 1.0 | 0.6 | 71.3 | 2.31 | |
2019 | Border | 150 | 326 | 176.1 | 202.5 | 19.9 | 11.5 | 41.3 | 2.72 |
Drip | 90 | 326 | 155.0 | 197.0 | 3.0 | 1.0 | 60.0 | 3.15 | |
Rainfed | 0.0 | 326 | 142.0 | 157.0 | 1.0 | 0.6 | 25.3 | 2.19 | |
Note: Ire is Irrigation amount; Pr is precipitation; Eva is evapotration; Tr is transpiration rate; D is drainage; R is runoff; Wbal is water balance and WUE is water use efficiency. |
Nitrogen fate | Nfer (kg ha–1) |
Nmin (kg ha–1) |
Nvol (kg ha–1) |
Nden (kg ha–1) |
Nlea (kg ha–1) |
Nup (kg ha–1) |
Nbal (kg ha–1) |
FNUE (kg kg–1) |
|
2018 | Border | 200 | 56.9 | 13.5 | 8.4 | 21.6 | 196.5 | 16.7 | 41.71 |
Drip | 200 | 43.0 | 18.0 | 2.0 | 0.7 | 190.0 | 32.3 | 54.51 | |
Rainfed | 200 | 36.0 | 31.0 | 3.0 | 0.8 | 187.0 | 14.2 | 30.98 | |
2019 | Border | 200 | 53.3 | 13.0 | 9.0 | 22.0 | 190.0 | 19.3 | 43.43 |
Drip | 200 | 39.0 | 21.0 | 3.0 | 0.8 | 192.0 | 22.2 | 51.14 | |
Rainfed | 200 | 33.0 | 29.0 | 5.0 | 0.9 | 179.0 | 19.1 | 30.61 | |
Note: Nfer is fertilization; Nmin is net mineralization; Nvol is ammonia volatilization; Nden is denitrification; Nlea is nitrogen leaching; Nup is crop Nitrogen uptake, Nbal is nitrogen balance and FNUE is Nitrogen use efficiency [51]. |