Wheat is the main crop in the world ranks after rice and the largest grain source of Pakistan. Among several reasons for diminishing wheat yield in Pakistan, water stress throughout the growing season decreases crop production because of the short life span. Two years (2015–16 and 2016–17) of field experiments were conducted to assess the impact of various water regimes (full irrigation, irrigation at 45, 60, and 75 mm potential soil moisture deficit (PSMD)) on the growth and yield of wheat. Maximum crop growth rate was recorded by application of irrigation at 45 mm PSMD. Application of irrigation at 45 mm PSMD ensured maximum radiation use efficiency regarding total dry matter production and grain yield. The maximum number of productive tillers, spike length, and grain yield were recorded under 45 mm PSDM treatment. The present results show that the effect of water is more pronounced regarding the growth and productivity of wheat. Application of irrigation at 45 mm PSMD ensures higher economical yield.

  • Physiological availability of water to plants is a key factor in crop production in the world.

  • Potential soil moisture deficit (PSMD) is a useful approach to save water.

  • Application of irrigation at 45 mm PSMD ensures higher economical yield.

Wheat is the largest grain source of Pakistan and ranks after rice as the main crop. It accounts for 80% of the cultivated area of about nine million hectares during Rabi season, therefore, it has an important role in formulation of new agricultural policies. Wheat flour is a key product with a big share of about 72% of Pakistan's everyday caloric consumption with per capita intake of about 124 kg/year, which is the highest in the world (Blackshaw et al. 2006; USDA 2017). The average yield of wheat crop in Pakistan is 2,827 kg ha−1 which is very low as compared to other countries (GOP 2020). In Pakistan, wheat yield could not achieve its maximum in comparison with different wheat producing countries in the world. Several ecological factors are responsible for limiting yield of crops (Bali & Sidhu 2020), and a small change in these from their normal array results in a decrease in plant growth and yield. Environmental variations result in numerous morphological, physiological, physiochemical, and molecular changes in plants. The world agriculture cereal production is also decreasing with climate change (Raza et al. 2019; Mubarik et al. 2021).

Pakistan has a good irrigation infrastructure, but climate change would result in inconsistency of the amount of water taken out at canal headworks and unpredictability in the frequency and intensity of rainfall. Furthermore, the average precipitation during July to September (monsoon season) is 138 mm, while in 2011–12 it received 237 mm precipitation which is more than normal precipitation. However, in the winter season of 2012, the amount of precipitation was 34 mm which was 51% less than the normal 71 mm, and consequently a decline in water amount taken out at canal headworks (GOP 2012). Hence, a decrease in canal water stores and unpredictability in precipitation are major threats to agriculture. The average need of water in rabi season (October to April) is 44.9 × 109 m3. However, severe scarcities of water will mean these requirements are not met (GOP 2012). Variability in precipitation is a burning issue for Pakistan, and it is a huge threat to the masses' food security in the long run. Water stress is a globally concerning issue and considered a complex hazard that directly influences the water balance of especially arid and semi-arid regions (Sur & Lunagaria 2020; Wang et al. 2020). Water shortage has a critical role in the final yield of crops (Holzman et al. 2018). Basal & Szabó (2020) reported that water stress is one of the most hazardous abiotic stresses that adversely influences crop yield. One of the national and universal concerns is to articulate different plans to manage existing water resources for agricultural use (Smith 2000).

For optimum growth and better yield of wheat recurrent irrigations should be applied because it is a rather sensitive crop to water stress (Alderfasi & Neilsen 2001). A better option or alternative is crop cultivation under a deficit irrigation scheme with less evapotranspiration to conserve water and other resources (Jalota et al. 2006). The ultimate objective of irrigation to crops is to meet water requirements and to boost the yield, in spite of the fact that water shortage is a severe issue for crop production in tropical areas owing to less or erratic rainfall patterns and climate change having a huge influence on tropical wheat production (Messmer et al. 2009). Irrigation water management is vital in all crops, especially in wheat because under- or overirrigation results in growth retardation of wheat and a decrease in yield. Most areas of Punjab often experience a dry weather environment, particularly during Rabi cropping season (Amin et al. 2020). Keeping in view the above discussion, the present study was planned (i) to evaluate the performance of wheat crop under different potential soil moisture deficit (PSMD) strategies and (ii) to explore PSMD as a helpful approach for irrigation scheduling in wheat.

Crop husbandry

The present study was carried out at Agronomic Farm, Department of Agronomy, University of Agriculture, Faisalabad-Pakistan. Soaking irrigation (Rouni) was applied before the sowing of the crop to maintain soil moisture at field capacity (25% v/v) during the wheat crop growing season of 2015–16 and 2016–17. Wheat cultivar Sahar-2006 was sown using a seed rate of 100 kg ha−1 in rows 20 cm apart using a manual single row hand drill. Physio-chemical properties of the soil at the experimental site (Table S1 in Supplementary Material) were analyzed before the crop was sown.

Recommended nutrient dose of N:P:K at 120:90:60 kg ha−1 was applied using urea, diammonium phosphate (DAP), and sulphate of potash (SOP) fertilizers as the source of nutrient. Nitrogen was applied in three equal parts, one-third of the total nitrogen was applied during seed bed preparation while the two remaining equal parts were applied at first and second irrigation. All phosphate and potash fertilizers were applied during seed bed preparation. Uniform weed management and plant protection measures were adopted throughout the experimentation cycle. Treatments included:

  • Full irrigation/irrigation at all growth stages (control)

  • Irrigation at 45 mm potential soil moisture deficit

  • Irrigation at 60 mm potential soil moisture deficit

  • Irrigation at 75 mm potential soil moisture deficit.

Before each irrigation, soil moisture percentage was determined by gravimetric method and depth of irrigation for the next irrigation was estimated as follows:
formula
where SMC = soil moisture content, Ww = weight of moist soil (as collected from the field at 45 cm depth), and Wd = weight of oven dried soil.
SMC (% on weight basis) was converted to SMC (% on volume basis) by the formula:
formula
Depth of irrigation (Di) was determined by the formula used by Fahad et al. (2019):
formula
where Di = depth of irrigation (cm) or crop water requirement in depth (cm), FC = field capacity (% on volume basis), SMC = soil moisture content (% on volume basis), BD = bulk density (g cm−3), and Dr = depth of root zone (cm).
A measured quantity of irrigation was applied according to the requirement (Di) to replenish the depleted soil moisture using the equation:
formula
where T = time in seconds for pre-determined amount of irrigation, A = area to be irrigated (m2), Di = depth of irrigation (m), Q = discharge of water channel at field head (flow of water through a point at m3 sec−1).

For measurement of discharge, a cutthroat flume was installed in the water channel and readings were observed at every irrigation. These readings were used to derive the water flow rate. The amount (mm) and date of irrigation are presented in Table 1. Data regarding prevalent weather condition throughout the crop season were obtained from a field meteorological laboratory situated in the vicinity of the experimental site (Figure S1).

Table 1

Date of irrigation water applied to crop during course of experimentation

2015–16
2016–17
DateI1-CI2-45I3-60I4-75I1-CI2-45I3-60I4-75Date
12/09/2015 70    70    12/12/2016 
12/16/2015  60    60   12/19/2016 
12/16/2015   60    60  12/19/2016 
01/06/2016 85   60 85   60 01/09/2017 
01/27/2016  60    60   01/30/2017 
01/30/2016 60    60    02/02/2017 
02/08/2016   60    60  02/11/2017 
02/22/2016 60 60  60 60 60  60 02/25/2017 
03/09/2016 60 60 60  60 60 60  03/12/2017 
03/16/2016    60    60 03/19/2017 
Total (mm) 335 240 180 180 335 240 180 180 Total (mm) 
2015–16
2016–17
DateI1-CI2-45I3-60I4-75I1-CI2-45I3-60I4-75Date
12/09/2015 70    70    12/12/2016 
12/16/2015  60    60   12/19/2016 
12/16/2015   60    60  12/19/2016 
01/06/2016 85   60 85   60 01/09/2017 
01/27/2016  60    60   01/30/2017 
01/30/2016 60    60    02/02/2017 
02/08/2016   60    60  02/11/2017 
02/22/2016 60 60  60 60 60  60 02/25/2017 
03/09/2016 60 60 60  60 60 60  03/12/2017 
03/16/2016    60    60 03/19/2017 
Total (mm) 335 240 180 180 335 240 180 180 Total (mm) 

I1 = Full Irrigation/Irrigation at all stages (control), I2 = irrigation at 45 mm potential soil moisture deficit, I3 = irrigation at 60 mm potential soil moisture deficit, I4=irrigation at 75 mm potential soil moisture deficit.

Growth, yield and their attributes

To determine the growth parameters, each plot was divided into two sub-plots. One of them was used for destructive biomass sampling and the other half was kept intact for final grain yield measurements. From each plot, a 900 cm2 area was harvested at ground level at intervals of every 15 days leaving appropriate borders. Fresh and dry weights of constituent fractions of the plant (leaf and stem) were determined. A sub-sample from each fraction was taken to dry in an oven to a constant weight. Crop growth rate was calculated as proposed by Hunt (1978):
formula
where W1 and W2 are the dry weights harvested at time intervals of t1 and t2, respectively.
From the measurements of leaf area and dry weights the following parameters were calculated. Leaf area was measured by using a leaf area meter (Model CI-202, CID, Inc.). Five gram of leaf laminae from each experimental unit was used for the estimation of leaf area and then converted to total leaf area of harvested samples. Leaf area index (LAI) was then calculated as the ratio of leaf area to land area (Watson 1952):
formula
Leaf area duration (LAD) was estimated as suggested by Hunt (1978):
formula
Net assimilation rate (NAR) was determined according to Hunt (1978):
formula
where ATDM = final total dry matter and LAD = leaf area duration.
The fraction of intercepted radiation (Fi) was calculated from measurements of LAI using the exponential equation suggested by Monteith & Elston (1983):
formula
where K = extinction coefficient for total solar radiation (Monteith & Elston 1983). A ‘k’ value of 0.45 was used for wheat.
The amount of intercepted light (Sa) was determined by multiplying Fi with incident PAR (Si) during the season as:
formula
The photosynthetically active radiation (PAR) was assumed as 50% of the total incident radiation.
Radiation use efficiency for TDM (RUE-TDM) and grain yield (RUE-GY) was calculated as the ratio of total biomass and grain yield to cumulative intercepted PAR (ΣSa):
formula
formula
Water use efficiency (WUE) on TDM basis (WUETDM) and grain yield basis (WUEGY) were calculated as the ratio of total biomass and grain yield to the total amount of irrigation applied and rainfall:
formula
formula
Fertile tillers per unit area, randomly selected from each experimental unit, were recorded manually at maturity. Ten plants were randomly selected from each treatment plot for measurement of plant height and spike length using a meter rod. The number of grains per spike and spikelets per spike were also counted manually. An electronic balance was used to record the 1,000 grain weight. Total wheat biomass was recorded with the help of a weighing balance. Grain yield from each plot was recorded by using an electronic weighing balance. Harvest index (HI) on a percentage basis was recorded for each treatment as ratio of grain yield to total biomass produced:
formula

Statistical analysis

Data regarding all parameters were analyzed using Fisher's analysis of variance technique and the LSD test (least significance difference) was carried out at a probability level of 0.05 for comparison of variances among treatment means (Steel et al. 1997).

Growth, yield and their attributes

Crop growth rate (CGR) was significantly (P < 0.05) influenced by different irrigation regimes (Figure 1). The results showed that maximum CGR (9.98 g m−2 d−1) was recorded by application of irrigation at 45 mm PSMD (Figure 1). Figure 2 depicts the effect of treatments on the fraction of intercepted radiation (Fi). The Fi reached a maximum value of 0.92 at 90 days after sowing, then it (0.61) gradually decreased up to maturity (Figure 2). Maximum net assimilation rate (NAR) was recorded under irrigation at 60 mm PSMD with PAR under full irrigation regimes (Table 2). A significant difference was observed in water use efficiency (WUE) on total dry matter (TDM) and grain yield (GY) basis under different irrigation regimes (Table 2). Maximum value of RUE on TDM (3.94 g MJ−1) and GY (1.18 g MJ−1) basis was noted from wheat crop irrigated at 45 mm PSMD (Tables 2 and 3). Maximum tillers per unit area were recorded by irrigation at 45 mm potential soil moisture deficit (330 m−2). The wheat crop, irrigated at 45 mm PSMD, resulted in maximum plant height (99.2 cm) (Table 3).

Table 2

Influence of irrigation regimes on net assimilation rate (NAR), photosynthetic active radiation (PAR), water use efficiency (WUE) by total dry matter (TDM) and grain yield (GY), and radiation use efficiency (RUE) by TDM

TreatmentsNAR (g m−2 d−1)
PAR (MJ m−2)
WUE by TDM (g m−2 mm−1)
WUE by GY (g m−2 mm−1)
RUE-TDM (g MJ−1)
2015–162016–17Mean2015–162016–17Mean2015–162016–17Mean2015–162016–17Mean2015–162016–17Mean
I1 5.29 5.28 5.29B 600 598 599A 3.68 3.66 3.67D 1.19 1.17 1.18D 3.75 3.73 3.74B 
I2 5.20 5.18 5.19C 589 588 588B 5.19 5.21 5.20C 1.83 1.80 1.82B 3.95 3.93 3.94A 
I3 5.55 5.57 5.56A 582 581 582BC 6.37 6.42 6.39A 1.97 1.96 1.97A 3.67 3.67 3.67C 
I4 5.49 5.50 5.50A 577 576 577C 5.53 5.53 5.53B 1.63 1.64 1.63C 3.45 3.44 3.45D 
Mean 5.38 5.37  587 586  5.19 5.21  1.65 1.64  3.71 3.69  
LSD Y = ns, I = 0.0887, Y × I = ns Y = ns, I = 6.0814, Y × I = ns Y = ns, I = 0.0675, Y × I = ns Y = ns, I = 0.0331, Y × I = ns Y = ns, I = 0.0212, Y × I = ns 
TreatmentsNAR (g m−2 d−1)
PAR (MJ m−2)
WUE by TDM (g m−2 mm−1)
WUE by GY (g m−2 mm−1)
RUE-TDM (g MJ−1)
2015–162016–17Mean2015–162016–17Mean2015–162016–17Mean2015–162016–17Mean2015–162016–17Mean
I1 5.29 5.28 5.29B 600 598 599A 3.68 3.66 3.67D 1.19 1.17 1.18D 3.75 3.73 3.74B 
I2 5.20 5.18 5.19C 589 588 588B 5.19 5.21 5.20C 1.83 1.80 1.82B 3.95 3.93 3.94A 
I3 5.55 5.57 5.56A 582 581 582BC 6.37 6.42 6.39A 1.97 1.96 1.97A 3.67 3.67 3.67C 
I4 5.49 5.50 5.50A 577 576 577C 5.53 5.53 5.53B 1.63 1.64 1.63C 3.45 3.44 3.45D 
Mean 5.38 5.37  587 586  5.19 5.21  1.65 1.64  3.71 3.69  
LSD Y = ns, I = 0.0887, Y × I = ns Y = ns, I = 6.0814, Y × I = ns Y = ns, I = 0.0675, Y × I = ns Y = ns, I = 0.0331, Y × I = ns Y = ns, I = 0.0212, Y × I = ns 

Means sharing the same letter did not differ significantly at P = 0.05.

I1 = Full irrigation/Irrigation at all stages (control), I2=irrigation at 45 mm potential soil moisture deficit, I3 = irrigation at 60 mm potential soil moisture deficit, I4 = irrigation at 75 mm potential soil moisture deficit, Y = sowing year, I = irrigation treatments, Y × I = interaction, ns = non-significant.

Table 3

Influence of irrigation regimes on RUE by GY, number of tillers, plant height, spike length and number of spikelets per spike

TreatmentsRUE-GY (g MJ−1)
Number of tillers (m−2)
Plant height (cm)
Spike length (cm)
Spikelets per spike
2015–162016–17Mean2015–162016–17Mean2015–162016–17Mean2015–162016–17Mean2015–162016–17Mean
I1 1.07 1.08 1.08B 311 302 306B 94.0 93.3 93.7B 11.9 11.8 11.8B 16.3 16.2 16.3B 
I2 1.19 1.17 1.18A 334 326 330A 99.7 98.7 99.2A 14.1 14.5 14.3A 17.3 17.5 17.4A 
I3 0.97 0.96 0.97C 297 289 293B 89.4 89.9 89.7C 11.0 11.5 11.3C 15.1 14.9 15.0C 
I4 0.84 0.83 0.84D 259 251 255C 84.5 85.3 84.9D 10.4 10.2 10.3D 13.6 13.4 13.5D 
Mean 1.02 1.01  300 292  91.9 91.8  11.9 12.0  15.6 15.5  
LSD Y = ns, I = 0.0138, Y × I = ns Y = ns, I = 22.120, Y × I = ns Y = ns, I = 1.0651, Y × I = ns Y = ns, I = 0.3433, Y × I = ns Y = ns, I = 0.1870, Y × I = ns 
TreatmentsRUE-GY (g MJ−1)
Number of tillers (m−2)
Plant height (cm)
Spike length (cm)
Spikelets per spike
2015–162016–17Mean2015–162016–17Mean2015–162016–17Mean2015–162016–17Mean2015–162016–17Mean
I1 1.07 1.08 1.08B 311 302 306B 94.0 93.3 93.7B 11.9 11.8 11.8B 16.3 16.2 16.3B 
I2 1.19 1.17 1.18A 334 326 330A 99.7 98.7 99.2A 14.1 14.5 14.3A 17.3 17.5 17.4A 
I3 0.97 0.96 0.97C 297 289 293B 89.4 89.9 89.7C 11.0 11.5 11.3C 15.1 14.9 15.0C 
I4 0.84 0.83 0.84D 259 251 255C 84.5 85.3 84.9D 10.4 10.2 10.3D 13.6 13.4 13.5D 
Mean 1.02 1.01  300 292  91.9 91.8  11.9 12.0  15.6 15.5  
LSD Y = ns, I = 0.0138, Y × I = ns Y = ns, I = 22.120, Y × I = ns Y = ns, I = 1.0651, Y × I = ns Y = ns, I = 0.3433, Y × I = ns Y = ns, I = 0.1870, Y × I = ns 

Means sharing the same letter did not differ significantly at P = 0.05.

I1 = Full irrigation/Irrigation at all stages (control), I2 = irrigation at 45 mm potential soil moisture deficit, I3 = irrigation at 60 mm potential soil moisture deficit, I4 = irrigation at 75 mm potential soil moisture deficit, Y = sowing year, I = irrigation treatments, Y × I = interaction, ns = non-significant.

Figure 1

Effect of deficit irrigation on crop growth rate of wheat. I1 = Full irrigation/Irrigation at all stages, I2 = irrigation at 45 mm potential soil moisture deficit, I3 = irrigation at 60 mm potential soil moisture deficit, I4 = irrigation at 75 mm potential soil moisture deficit.

Figure 1

Effect of deficit irrigation on crop growth rate of wheat. I1 = Full irrigation/Irrigation at all stages, I2 = irrigation at 45 mm potential soil moisture deficit, I3 = irrigation at 60 mm potential soil moisture deficit, I4 = irrigation at 75 mm potential soil moisture deficit.

Close modal
Figure 2

Fraction of intercepted radiation (Fi) as affected by deficit irrigations. I1 = Full irrigation/Irrigation at all stages, I2 = irrigation at 45 = mm potential soil moisture deficit, I3 = irrigation at 60 mm potential soil moisture deficit, I4 = irrigation at 75 mm potential soil moisture deficit, DAS = Days after sowing.

Figure 2

Fraction of intercepted radiation (Fi) as affected by deficit irrigations. I1 = Full irrigation/Irrigation at all stages, I2 = irrigation at 45 = mm potential soil moisture deficit, I3 = irrigation at 60 mm potential soil moisture deficit, I4 = irrigation at 75 mm potential soil moisture deficit, DAS = Days after sowing.

Close modal

Application of irrigation at 45 mm PSMD resulted in maximum spike length of 14.3 cm. Spikelets per spike were observed to be maximum in the treatment where wheat crop was irrigated at 45 mm PSMD (17.4), which was followed by the treatment which was irrigated at all growth stages of wheat (16.3) (Table 3). Maximum number of grains per spike (43) was observed in wheat crop irrigated at all growth stages, which was statistically similar to the treatment irrigated at 45 mm PSMD (42.1) (Table 4). Similarly, irrigation application at 45 mm PSMD produced maximum 1,000 grain weight (44.3 g) followed by the treatment irrigated at all growth stages (41.6 g). Maximum grain yield was observed in the wheat crop by application of irrigation at 45 mm PSMD. The biological yield was observed as maximum under the treatment which received irrigation at 45 mm PSMD (12,908 kg ha−1). Wheat crop, irrigated at a moisture deficit of 45 mm PSMD, resulted in maximum harvest index (30.08%) (Table 4).

Table 4

Influence of irrigation regimes on number of grains per spike, 1,000-grain weight, grain yield, biological yield and harvest index

TreatmentsGrains per spike
1,000 Grain weight (g)
Grain yield (kg ha−1)
Biological yield (kg ha−1)
Harvest index (%)
2015–162016–17Mean2015–162016–17Mean2015–162016–17Mean2015–162016–17Mean2015–162016–17Mean
I1 43.3 42.7 43.0A 41.8 41.5 41.6B 3,511 3,493 3502B 12,284 12,272 12,278B 28.58 28.46 28.52B 
I2 42.6 41.7 42.1A 44.5 44.2 44.3A 3,881 3,885 3883A 12,900 12,916 12,908A 30.09 30.08 30.08A 
I3 40.2 41.7 40.9B 37.5 37.2 37.4C 3,177 3,164 3171C 12,103 12,099 12,101C 26.25 26.15 26.20C 
I4 37.7 36.3 37.0C 35.6 35.3 35.5D 2,806 2,792 2799D 11,349 11,329 11,339D 24.72 24.65 24.68D 
Mean 41.0 40.6  39.9A 39.5B  3344A 3334B  12,159 12,154  27.41A 27.34B  
LSD Y = ns, IM = 0.9471, Y × I = ns Y = 0.226; I = 0.319, Y × I = ns Y = 7.301; I = 10.32, Y × I = ns I = 8.2864, Y × I = ns Y = 0.066; I = 0.093, Y × I = ns 
TreatmentsGrains per spike
1,000 Grain weight (g)
Grain yield (kg ha−1)
Biological yield (kg ha−1)
Harvest index (%)
2015–162016–17Mean2015–162016–17Mean2015–162016–17Mean2015–162016–17Mean2015–162016–17Mean
I1 43.3 42.7 43.0A 41.8 41.5 41.6B 3,511 3,493 3502B 12,284 12,272 12,278B 28.58 28.46 28.52B 
I2 42.6 41.7 42.1A 44.5 44.2 44.3A 3,881 3,885 3883A 12,900 12,916 12,908A 30.09 30.08 30.08A 
I3 40.2 41.7 40.9B 37.5 37.2 37.4C 3,177 3,164 3171C 12,103 12,099 12,101C 26.25 26.15 26.20C 
I4 37.7 36.3 37.0C 35.6 35.3 35.5D 2,806 2,792 2799D 11,349 11,329 11,339D 24.72 24.65 24.68D 
Mean 41.0 40.6  39.9A 39.5B  3344A 3334B  12,159 12,154  27.41A 27.34B  
LSD Y = ns, IM = 0.9471, Y × I = ns Y = 0.226; I = 0.319, Y × I = ns Y = 7.301; I = 10.32, Y × I = ns I = 8.2864, Y × I = ns Y = 0.066; I = 0.093, Y × I = ns 

Means sharing the same letter did not differ significantly at P = 0.05.

I1 = Full irrigation/Irrigation at all stages (control), I2 = irrigation at 45 mm potential soil moisture deficit, I3 = irrigation at 60 mm potential soil moisture deficit, I4 = irrigation at 75 mm potential soil moisture deficit, Y = sowing year, I = irrigation treatments, Y × I = interaction, ns = non-significant.

The relationship between different growth and yield parameters is calculated and presented in Figure 3. The relationship between crop growth rate and total dry matter was recorded as linear and was strongly positive giving an R2 value of 0.95. Regression relation of the number of productive tillers with grain yield was positive and linear giving R2 values of 0.80. The spike length had a good positive and linear relationship with grain yield, and the R2 value was 0.96. Number of grains per spike gave a regression value of 0.84 in relation to grain yield. The relationship between grain yield and harvest index was linear and positive, and gave a regression value of 0.84. The regression relation of radiation use efficiency with grain yield was strongly positive and linear, and the value of regression was 0.99 (Figure 3).

Figure 3

Co-relation among various attributes of wheat cultivated under different irrigation regimes.

Figure 3

Co-relation among various attributes of wheat cultivated under different irrigation regimes.

Close modal

Appropriate irrigation frequency and intensity is a critical factor for the optimum production of field crops. Various irrigation regimes significantly affected the growth and yield contributing parameters of wheat crop. Irrigation water management is vital in all crops, especially in wheat because under- or overirrigation resulted in growth retardation of wheat and a decrease in yield. A better option or alternative is crop cultivation under a deficit irrigation scheme with less evapotranspiration to conserve water and other resources (Jalota et al. 2006). In the present study, application of irrigation at 45 PSMD produced the maximum yield and was found more responsive regarding the growth and yield of wheat crop. Zhang et al. (2004) stated that soil water deficit played a significant role in crop production. Their findings are in line with outcomes of the present experimentation. They stated that severe soil water deficit significantly reduced the economical yield of wheat crop as compared to slight soil water deficit. They also stated that evapotranspiration of crop also depended on irrigation amount. Qiu et al. (2008) also observed WUE on the basis of photosynthesis and biomass production and found a positive correlation between WUE and winter wheat grain yield. Application of irrigation at 45 mm PSMD produced the highest grain yield and its components in the present study. These outcomes are also supported by Bashir et al. (2016). PSMD at 45 mm was found to be the most appropriate throughout the experimentation because it appears responsible for maximum plant height, number of grains per spike, thousand grain weight, productive tillers, and economic yield.

The current study results revealed that deficit irrigation at 45 mm PSMD is more water-saving than conventional practices of irrigation of traditional farmers. The consequences of the current experiment are comparable with the findings of Li et al. (2005), who found that the highest number of tillers were observed when wheat crop was irrigated at critical growth stages and the least number of tillers when irrigation was applied at longer intervals. Application of water deficit irrigation at 45 PSMD produced taller plants comparatively, as recorded in the present study. The reason for taller plants might be the application of irrigation at all crucial stages, i.e., tillering, elongation, booting, anthesis and grain development stages, while similar plant height attained when irrigation at grain filling stage was skipped was due to the fact that the vegetative growth stage of the plant was completed at that stage. Shorter plants developed when irrigation was skipped at tillering and anthesis stage which resulted in lack of adequate moisture for stem elongation. Plant height was favored as irrigation was provided to wheat crop at all critical stages of crop growth (Qiu et al. 2008). Kazemi et al. (2021) also reported that water stress was responsible in the significant reduction of yield attributing traits.

Water deficit conditions resulted in the reduced number of grains per spike as well as in productive tillers per unit area in wheat crop (Sayyah et al. 2015). Our results are in agreement with Moayedi et al. (2010) and Plaut et al. (2004), who reported reduced 1,000 grain weight and grain yield in wheat due to water deficit during grain filling stage. According to Shehzad et al. (2012), more grain weight was recorded under 50 mm PSMD treatment, and increasing the water deficit level up to 75 mm PSMD significantly reduced the grain weight. A similar trend of 1,000 grain weight was found in relation to deficit irrigations in the current study. Salsinha et al. (2021) stated that water stress affected various growth parameters negatively which ultimately reduced the economic yield.

The correlation between productive tillers and grain yield is of prime importance (Figure 3). A positive correlation between the number of productive tillers and grain yield was observed in wheat under drought stress conditions (Mohammadi et al. 2011). Our findings of positive correlation between productive tillers and grain yield are also in line with the results of Peymaninia et al. (2012) and Chalabi & Rashidi (2012). Taiz & Zeiger (2006) reported that measurement of biomass is the key factor to determine drought stress in plants. Therefore, biomass has a pivotal role in increasing grain yield. A positive correlation between grain yield and harvest index was observed (Figure 3). The study of Kirigwi et al. (2004) supports our findings as it reported a positive correlation of grain yield with biological yield and harvest index in varying water deficit levels. The current study results revealed that deficit irrigation at 45 mm PSMD is more water-saving than conventional practices of irrigation of traditional farmers. PSMD at 45 mm is also found most appropriate throughout the experimentation because it appears responsible for maximum plant height, number of grains per spike, 1,000 grain weight, productive tillers, and economic yield.

The physiological availability of water to plants is a key factor in crop production throughout the world. Data suggest that there is still huge scope to increase wheat grain yield by using deficit irrigation sensibly. Under a climate change-induced water shortage scenario, lack of adequate moisture at early growth stages has a more detrimental effect than at later growth stages. Maximum growth and yield of wheat crop in the field can be achieved by applying irrigation at 45 mm PSMD. Hence, irrigation at 45 mm PSMD is suggested for achieving high productivity.

The authors extend their appreciation to the researchers supporting project (No. RSP. 2021/173), King Saud University, Riyadh, Saudi Arabia. This research did not receive any specific funding. The authors declare no conflict of interest. The data that support the outcomes of current experimentation are available from the corresponding author upon reasonable request.

All relevant data are included in the paper or its Supplementary Information.

Amin
M.
Khan
M. R.
Hassan
S. S.
Khan
A. A.
Imran
M.
Goheer
M. A.
Hina
S. M.
Perveen
A.
2020
Monitoring agricultural drought using geospatial techniques: a case study of Thal region of Punjab, Pakistan
.
Journal of Water and Climate Change
11
(
S1
),
203
216
.
Bali
A. S.
Sidhu
G. P.
2020
Growth and morphological changes of agronomic crops under abiotic stress
. In:
Agronomic Crops
(
Hassanuzzaman
M.
, ed.).
Springer
,
Singapore
, pp.
1
11
.
Basal
O.
Szabó
A.
2020
Physiology, yield and quality of soybean as affected by drought stress
.
Asian Journal of Agriculture and Biology
8
(
3
),
247
252
.
https://doi.org/10.35495/ajab.2019.11.505
.
Bashir
M. U.
Wajid
S. A.
Ahmad
A.
Iqbal
M.
2016
Potential soil moisture deficit: an alternative approach for irrigation scheduling in wheat
.
International Journal of Agriculture and Biology
18
,
16
22
.
Blackshaw
R. E.
Odonovan
J. T.
Harker
K. N.
Clayton
G. W.
Stougaard
R. N.
2006
Reduced herbicide doses in field crops: a review
.
Weed Biology and Management
6
,
10
17
.
Chalabi
Y. S.
Rashidi
V.
2012
Selection indices in the improvement of wheat grain yield on drought stress conditions
.
African Journal of Agricultural Research
7
,
1177
1183
.
Fahad
M.
Wajid
S. A.
Ahmad
A.
Cheema
M. J. M.
2019
Response of wheat cultivars to deficit irrigation under semiarid conditions of Faisalabad
.
International Journal of Agriculture and Biology
21
,
1004
1012
.
GOP
.
2012
Economic Survey of Pakistan
.
Government of Pakistan
,
Islamabad, Pakistan
, pp.
21
24
.
GOP
.
2020
Economic Survey of Pakistan
.
Government of Pakistan
,
Islamabad, Pakistan
, pp.
21
22
.
Holzman
M. E.
Carmona
F.
Rivas
R.
Niclòs
R.
2018
Early assessment of crop yield from remotely sensed water stress and solar radiation data
.
ISPRS Journal of Photogrammetry and Remote Sensing
145
,
297
308
.
Hunt
R.
1978
Plant Growth Analysis
.
Edward Arnold
,
London, UK
, pp.
26
38
.
Jalota
S. K.
Sood
A.
Chachal
G. B. S.
Choudhury
B. U.
2006
Crop water productivity of cotton wheat system as influenced by deficit irrigation, soil texture and precipitation
.
Agricultural Water Management
84
,
137
146
.
Kazemi
S.
Zakerin
A.
Abdossi
V.
Moradi
P.
2021
Fruit yield and quality of the grafted tomatoes under different drought stress conditions
.
Asian Journal of Agriculture and Biology
2021
(
1
).
https://doi.org/10.35495/ajab.2020.03.164
Kirigwi
F. M.
Ginkel
M. V.
Trethowan
R.
Sears
R. G.
Rajaram
S.
Paulsen
G. M.
2004
Evaluation of selection strategies for wheat adaptation across water regimes
.
Euphytica
135
,
361
371
.
Li
J.
Inanga
S.
Li
Z.
Eneji
A.
2005
Optimizing irrigation scheduling for winter wheat in the North China Plain
.
Agricultural Water Management
76
,
8
23
.
Messmer
R.
Fracheboud
Y.
Banziger
M.
Vargas
M.
Stamp
P.
Ribaut
M.
2009
Drought stress and wheat QTL by environment interactions and stability of QTLs across environments for yield components and secondary traits
.
Theoretical and Applied Genetics
199
,
913
930
.
Moayedi
A. A.
Boyce
A. N.
Barakbah
S. S.
2010
The performance of durum and bread wheat genotypes associated with yield and yield components under different water deficit conditions
.
Australian Journal of Basic and Applied Sciences
4
,
106
113
.
Mohammadi
M.
Karimizadeh
R.
Shefazadeh
M. K.
Sadeghzadeh
B.
2011
Statistical analysis of durum wheat yield under semi-warm dry land condition
.
Australian Journal of Basic and Applied Sciences
5
,
1292
1297
.
Monteith
J. L.
Elston
J. F.
1983
Performance and productivity of foliage in the field
. In:
The Growth and Functioning of Leaves
(
Dale
J. E.
Milthorpe
F. L.
, eds).
Cambridge University Press
,
Butterworths, London, UK
, pp.
499
518
.
Mubarik
M. S.
Khan
S. H.
Sajjad
M.
Raza
A.
Hafeez
M. B.
Yasmeen
T.
Rizwan
M.
Ali
S.
Arif
M. S.
2021
A manipulative interplay between positive and negative regulators of phytohormones: a way forward for improving drought tolerance in plants
.
Physiologia Plantarum
.
https://doi.org/10.1111/ppl.13325
Peymaninia
Y.
Valizadeh
M.
Shahryari
R.
Ahmadizadeh
M.
Habibpour
M.
2012
Relationship among morpho-physiological traits in bread wheat against drought stress at presence of a leonardite derived humic fertilizer under greenhouse condition
.
International Research Journal of Applied and Basic Sciences
3
,
822
830
.
Qiu
G. Y.
Wang
L.
He
X.
Zhang
X.
Chen
S.
Chen
J.
Yang
Y.
2008
Water use efficiency and evapotranspiration of winter wheat and its response to irrigation regime in the north China plain
.
Agricultural and Forest Meteorology
148
,
1848
1859
.
Raza
A.
Razzaq
A.
Mehmood
S. S.
Zou
X.
Zhang
X.
Lv
Y.
Xu
J.
2019
Impact of climate change on crops adaptation and strategies to tackle its outcome: a review
.
Plants
8
(
2
),
34
.
Salsinha
Y. C. F.
Maryani
I. D.
Purwestri
Y. A.
Rachmawati
D.
2021
Morphological and anatomical characteristics of Indonesian rice roots from East Nusa Tenggara contribute to drought tolerance
.
Asian Journal of Agriculture and Biology
2021
(
1
).
https://doi.org/10.35495/ajab.2020.05.304.
Sayyah
S. S.
Ghobadi
M.
Mansoorifar
S.
Zebarjadi
A. R.
2015
The yield of wheat genotypes associated with yield components under irrigated and drought stress after anthesis
.
Archives of Agronomy and Soil Science
61
,
1743
1755
.
Shehzad
M. A.
Maqsood
M.
Iqbal
S.
Saleem
M.
Hassan
M. U.
Ahmad
W.
2012
Impact of nitrogen nutrition and moisture deficits on growth, yield and radiation use efficiency of wheat (Triticum aestivum L.)
.
African Journal of Biotechnology
11
,
13980
13987
.
Steel
R. G. D.
Torrie
J. H.
Dickey
D. A.
1997
Principles and Procedures of Statistics. A Biometrical Approach
, 3rd edn.
McGraw Hill
,
New York, USA
, pp.
400
428
.
Taiz
L.
Zeiger
E.
2006
Plant Physiology
, 4th edn.
Sinauer Associates
,
Sunderland, MA, USA
.
USDA
.
2017
GAIN Report Number PK1704
.
USDA Foreign Agricultural Service
,
Washington DC, USA
.
Watson
D. J.
1952
The physiological basis of variation in yield
.
Advances in Agronomy
4
,
101
145
.
Zhang
Y.
Kendy
E.
Qiang
Y.
Changming
L.
Yanjun
S.
Hongyong
S.
2004
Effect of soil water deficit on evapotranspiration, crop yield and water use efficiency in the North China Plain
.
Agricultural Water Management
64
,
107
122
.
This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY 4.0), which permits copying, adaptation and redistribution, provided the original work is properly cited (http://creativecommons.org/licenses/by/4.0/).

Supplementary data