Response and Modeling of Hybrid Maize Seed Vigor to Water Deficit at Different Growth Stages

Research is imperative to predict seed vigor of hybrid maize production under water deficit in arid areas. Field experiments were conducted in 2018 and 2019 in arid areas of northwestern China to investigate the effects of different irrigation strategies at various growth stages with drip irrigation under film mulching on grain yield, kernel weight, seed protein content, and seed vigor of hybrid maize (Zea mays L.). Water deficit at vegetative, flowering, and grain-filling stages was considered and a total of 16 irrigation treatments was applied. A total of 12 indices of germination percentage, germination index (GI), shoot length (SL), and root length (RL) under different germination conditions (standard germination and accelerated aging); electrical conductivity (EC) of the leachate; and activities of peroxidase, catalase, and superoxide dismutase in seeds were measured and analyzed using the combinational evaluation method (CEM). Furthermore, five water production functions (Blank, Stewart, Rao, Jensen, and Minhas) were used to predict seed vigor evaluated by CEM under water deficit. The results showed that leachate EC was higher under water deficit than that under sufficient irrigation. The SL, RL, and GI of different germination conditions increased under water deficit at the flowering stage. The Rao model was considered the best fitted model to predict the vigor of hybrid maize seeds under water deficit, and an appropriate water deficit at the flowering stage is recommended to ensure high seed vigor of hybrid maize production with drip irrigation under film mulching. Our findings would be useful for reducing crop water use while ensuring seed vigor for hybrid maize production in arid areas.


Introduction
Hybrid maize was commercially grown for the first time in the early 1930s in the United States [1]. Hybrid maize lines show increased yield and disease resistance compared with parental lines, and over 95% of the corn-growing regions worldwide cultivate hybrid maize varieties [2]. With the expansion of the area under cultivation of this crop, the production of hybrid maize seeds has rapidly increased globally. Notably, hybrid maize seed production requires suitable environmental conditions, such as temperature and moisture, as well as technical and professional expertise for industrial production [3]. Over 60% of the total hybrid maize seed production in the world occurs in China (23%), India (15%), the United States (9%), Brazil (5%), Mexico (5%), and Argentina (4%) [4].
Seed vigor involves the combination of all traits that determine the potential activity and performance of seeds during germination and establishment [5]. High seed vigor significantly increases grain yield [6]. Seed vigor tests are divided into physiological methods, which measure germination characteristics, and biochemical methods, which measure specific biochemical reactions

Experimental Design
The growth of hybrid maize over the cropping cycle and seed production was divided into five stages: seedling, vegetative, flowering, grain-filling, and maturity stage. The vegetative stage lasted from the six-leaf stage to before the tasseling of inbred female lines, and the flowering stage lasted from the tasseling of inbred female lines to the end of pollination. A total of 16 treatments were delivered in each year (Table 2). A randomized block design with three replicates per irrigation treatment was used. The experimental plot was 4.8 m wide and 46 m long, and a buffer area (6 m wide and 46 m long) Water 2020, 12, 3289 4 of 20 separated different treatments. The inbred female lines were planted in the buffer zones. In this study, we assumed that the irrigation efficiency is 1, and the irrigation depth per event of sufficient irrigation treatment (I) was calculated as follows: where ET 0i is the reference evapotranspiration at the ith day in a given irrigation period (mm) (calculated according to the United Nations Food and Agriculture Organization (FAO) Penman-Monteith equation [40]); n is the length of the irrigation period in days; P e is the effective precipitation (difference between precipitation per rainfall event and 5 mm) in a given irrigation period (mm); and K c is the crop coefficient in different growth stages.  (4) In Equation (1), K c values were 0.33, 0.90, 1.10, 1.02, and 0.80 at the seedling, vegetative, flowering, grain-filling, and maturity stages, respectively [41]. And ET 0i was calculated by meteorological data in 2018 and 2019. Under the control treatment (V2F2G2), the maize inbred lines were irrigated every 10 days starting from 5 June 2018 and 4 June 2019, which followed local practices. The two irrigation levels applied at the vegetative stage were 100% I and 50% I, and the three irrigation levels applied at the flowering and grain-filling stages were 100% I, 50% I, and 0 I. Because the local hybrid maize cultivars are no irrigated at the seedling and maturity stages, no irrigation treatments were applied at these two stages. The irrigation treatments, depth, and times during different growth stages are presented in Table 2.
Meteorological data were gathered by an automatic weather station (Hobo, Onset Computer Corp., Bourne, MA, USA) about 50 m away from the experimental field and logged every 15 min. The meteorological data of all growth stages in both years are presented in Table 3. Table 3. Daily average solar radiation (R s ), air temperature (T a ), relative humidity (RH), wind speed (S w ), vapor pressure deficit (VPD), reference evapotranspiration (ET 0 ), and total precipitation (P tot ) at all growth stages in 2018 and 2019.
Year In each plot, 50 plants of female inbred lines were randomly harvested for yield and kernel weight (KW) measurements. Grains were dried at 105 • C for 30 min and then at 75 • C to a constant weight. A total of 300 grains were randomly selected from each plot and weighted to a constant KW. Seed yield was calculated based on yield per plant and female plant density. Seed yield and KW were converted to weight with kernel water content of 13%. Seed protein content (SPC, Nitrogen × 6.25) was determined using the micro-Kjeldahl method [42].

Evapotranspiration
Soil water content was measured using Diviner 2000 (Sentek Pty Ltd., Adelaide, Australia) within a depth of 1.0 m at 0.1 m intervals every 3-5 d. In addition, soil water content was measured before and after each irrigation and rainfall events. In each plot, one PVC tube was installed under bare soil and another under mulching. A calibration of the Diviner 2000 readings against gravimetric measurements of soil water content was performed for all the tubes separately. Crop evapotranspiration was estimated assuming no water uptake below 1 m, and using the soil water balance method as follows: where ET is evapotranspiration (mm); I is the total irrigation depth (mm); R is surface runoff (mm); D is deep drainage (mm); CR is capillary rise from a shallow water table towards the root zone (mm); ∆SF is difference between subsurface flow in and out of the root zone (mm); and ∆S is the variation in soil water storage (mm). ∆S can be calculated as follows: where W t1 and W t2 are the mean water contents in the root zone at times t 1 and t 2 (cm −3 cm −3 ); and H is soil depth of the planned moisture layer, and here considered as 1 m.
Because some fluxes such as SF, D and CR are difficult to assess and short time periods cannot be considered [40]. The experimental site is flat, and surface runoff is not intensive; therefore, R is considered as 0. Therefore, Equation (2) can be simplified as follows: Water 2020, 12, 3289 6 of 20

Germination Test
For SG test, seeds (100 seeds × 4 replicates) were placed in the germination chamber at 25 • C under 100% relative humidity and 12 h of light per day. Deterioration using AA test is associated with emergence and storage potential. Seeds (100 seeds × 4 replicates) were aged for 72 h at 45 • C under 100% relative humidity. After aging, the seeds were air-dried. GP of SG and AA were evaluated every day until the 7th day according to the recommendations of ISTA [8].

Seedling Growth Test
SLs and RLs of 10 representative seedlings per replicate were measured with a ruler on the 7th day of different germination tests. These indices reflect the robustness of the seedlings.

Germination Index (GI)
GI reflects the germination rate of a seed lot. A higher GI indicates a faster germination rate and can be calculated as follows: where GI is the germination index; Gt is the number of germinated seeds per day; and Dt is the number of days.

Leachate EC
Approximately 40 g of seeds per replication were soaked in 250 mL of deionized water for 24 h at 25 • C, and the beaker was sealed with plastic wrap. Leachate EC of sample was measured using a conductivity meter (S230, Mettler-Toledo International Inc., Greifensee, Switzerland), and the measurement was repeated four times. The leachate EC per gram of seed weight for each sample can be calculated as follows: where EC is the electrical conductivity of leachate per unit mass [µS (cm g) −1 ]; n is the number of repetitions; S i is the electrical conductivity of the samples (µS cm −1 ); S 0 is the electrical conductivity of deionized water (µS cm −1 ); and m i is the mass of seeds (g).

Enzyme Activity Assay
Two grams of seeds per replication preserved in liquid nitrogen were ground in a mortar and homogenized with 0.4 g polyvinylpyrrolidone in 20 mL of 0.1 M phosphate buffer (pH 7.8) containing 2 mM dithiothreitol, 0.1 mM EDTA and 1.25 mM PEG-4000. The homogenate was centrifuged at 16,000× g for 20 min. The supernatant obtained was used as an enzyme extract. All steps of the extraction procedure were carried out at 1-4 • C. The specific activities of the three enzymes were expressed as units of enzyme activity per milligram of protein, and replicated 15 times per treatment for each enzyme.
SOD activity (U mg −1 protein) was assayed by measuring the inhibition of photochemical reduction of nitroblue tetrazolium (NBT) [43]. The reaction mixture contained 75 µM NBT, 13 mM methionie, 2 µM riboflavine, 10 µM EDTA in 0.1 M phosphate buffer (pH 7.8), and 0.1 mL of enzyme extract to a final volume of 3.3 mL. Glass test tubes containing the mixture were placed in a light incubator for 20 min at 25 • C, and the tubes that were not illuminated served as blanks. Absorbance at 560 nm was recorded using a spectrophotometer (UV-2600, Shimadzu Corporation, Kyoto, Germany). A single unit of SOD was defined as the level of enzyme activity required to inhibit the photoreduction of NBT to blue formazan by 50%.
CAT activity (U mg −1 protein) was determined according to the method reported by Bailly et al. [13]. Hydrogen peroxide (H 2 O 2 ) can be decomposed by CAT, and the absorbance of the reaction solution Water 2020, 12, 3289 7 of 20 decreases with increasing reaction time. The enzyme assay contained 0.1 mM H 2 O 2 in 50 mM phosphate buffer (pH 7.8) and 0.2 mL of enzyme extract to a final volume of 3 mL. The change in absorbance was recorded using the spectrophotometric method at 240 nm at 30 s intervals. A single unit of CAT activity was defined as a 0.1 change in optical density (OD) per minute.
POD activity (U mg −1 protein) was assayed according to the method reported by Shannon et al. [44]. Under POD catalysis, hydrogen peroxide oxidizes guaiacol to form a brown product. The enzyme assay contained 19 µL guaiacol, 30% H 2 O 2 of 28 µL in 0.2 mM phosphate buffer (pH 6.0) and 1mL enzyme extract to a final volume of 4 mL. The change in absorbance was recorded using spectrophotometry at 470 nm at 30 s intervals. A single unit of POD activity was defined as a 0.01 change in OD per minute.
The protein content of the enzyme extraction solution was measured using Coomassie brilliant blue staining [45], and bovine serum albumin was used as the standard for setting up the calibration curves.

Assessment Methods
In this study, a total of 12 seed vigor indices (GP, GI, SL, and RL under SG and AA conditions, EC, POD, SOD and CAT) were evaluated by GRA, PCA and CEM. PCA and CEM were performed described elsewhere [31,36], and GRA is calculated by the following steps.
(1) Normalize the seed vigor index. The normalized value of the-smaller-the-better index (EC) is expressed as follows: The normalized values of the-greater-the-better indices (remaining 11 indices) are expressed as follows: where x ij is normalized value of seed vigor index; y ij is the original seed vigor index; i = 1, 2, . . . , m; j = 1, 2, . . . , n; y jmax is the maximum value of the jth index; and y jmin is the minimum value of the jth index. In this study, m = 16 and n = 12.
(2) The grey relational coefficient (ξ ij ) is calculated as follows: where x j 0 is the ideal normalized result for the jth index; and ρ = 0.5.
(3) The weighted grey relational grade for the ith treatment (R i * ) is calculated as follows: where W X is the weight of each seed vigor index, derived from the entropy method. And the entropy method has been detailed in Ding et al. [34].

Water-Seed Vigor Models
The crop water production functions using the Blank [46], Stewart [47], Rao [48], Jensen [49], and Minhas [50] models were applied to simulate the seed vigor of hybrid maize under water deficit in the present study.
Blank model: Water 2020, 12, 3289 8 of 20 Stewart model: Rao model: Jensen model: Minhas model: where V i is the seed vigor under the ith treatment evaluated by CEM; V CK is the seed vigor under sufficient irrigation (V2F2G2) evaluated by CEM; a = 1, 2, and 3 indicate the vegetative, flowering, and grain-filling stages, respectively; ET ia is evapotranspiration under the ith treatment at the ath growth stage (mm); ET CKa is evapotranspiration under sufficient irrigation at the ath growth stage (mm); and λ a , δ a , A a , B a , and C a indicate water sensitivity indices at the ath growth stage.

Data Analysis
Because of high rainfall at the grain-filling stage in 2018, the water-seed vigor models were calibrated using 12 irrigation treatments (V2F2G2, V2F1G0, V2F0G2, V1F2G0, V1F1G1, and V1F0G2 in each year) and validated using the remaining 20 irrigation treatments in 2018 and 2019. SPSS 21.0 (SPSS Inc., USA) was used for the one-way analysis of variance, verification normality and homogeneity of variance, and comparison of mean values with the Duncan multiple range tests at a 5% probability level. SAS 9.3 (SAS Institute, USA) was used for PCA and non-linear fitting. Model performance was evaluated based on relative root-mean-square error (RRMSE), coefficient of determination (R 2 ), average relative error (ARE), modeling efficiency (EF), and agreement index (d IA ). The equations for computing these indicators are the following: Water 2020, 12, 3289 9 of 20 where n is the sample size; P i is the predicted value; M i is the measured value; P is the mean of the predicted values; and M is the mean of the measured values.

Results and Discussion
3.1. Soil Water Content, Evapotranspiration, Yield, Kernel Weight, and Seed Protein Content Figure 1 shows the trends of mean soil volumetric water content in the 0-1 m soil layer under different irrigation treatments in 2018 and 2019. In this study, the lower limit of readily available water was 0.20 cm −3 cm −3 [40]. Soil water contents in all treatments were significantly affected by irrigation and rainfall, and fluctuated differently between the wilting point and field capacity, and decreased with the reduction of irrigation depth at different growth stages.
There were no significant differences on all the indices among the different irrigation treatments at the grain-filling stage in 2018 due to heavy rainfall (Table 3). ET, yield, KW, and SPC of hybrid maize seeds under different treatments are presented in Table 4.   Irrigation of V2F2G2; rainfall; field water capacity; wilting point; lower limit of readily available water, and here considered as 0.20 cm −3 cm −3 [40]. Irrigation of V2F2G2; rainfall; field water capacity; wilting point; lower limit of readily available water, and here considered as 0.20 cm −3 cm −3 [40]. Irrigation of V2F2G2; rainfall; field water capacity; wilting point; lower limit of readily available water, and here considered as 0.20 cm −3 cm −3 [40].   [40].  lower limit of readily available water, and here considered as 0.20 cm −3 cm −3 [40].
KW is an important component of maize yield. In this study, water deficit at the flowering stage significantly increased KW in both 2018 and 2019. This result is consistent with that of the study conducted by NeSmith and Ritchie [51]. In contrast, some studies showed that water deficit at the flowering stage decreased KW [52][53][54], while others found no effect [24,55,56]. Water deficit at the flowering stage significantly affected flowering characteristics and decreased kernel number [24]. The post-flowering assimilate availability per kernel under water deficit at the flowering stage may be the reason for increased KW. On the other hand, water deficit at vegetative and grain-filling stage decreased KW as a result of reduced the post-flowering assimilate availability per kernel.
SPC usually decreases with increasing total irrigation volume [57][58][59][60]. However, severe water stress delays the growth of leaf and stem cells, affects nutrient transport, and decreases SPC [61][62][63]. In this study, water deficit at the flowering stage did not significantly affect SPC in 2019, although water deficit at the vegetative and grain-filling stages significantly reduced SPC in both years. Furthermore, compared with the value under V2F2G2, SPC under V2F0G2 increased by 4.5% in 2018. Similar to KW, SPC is typically affected by changes in post-flowering assimilate availability per kernel [64,65]. Therefore, increased post-flowering assimilate availability per kernel under water deficit at the flowering stage may be another reason for increased SPC.

Seed Vigor
Cell membrane damage caused by water stress may increase the EC of seed leachate. Values of leachate EC under different irrigation treatments are shown in Figure 2. The results showed that EC significantly increased under water deficit at different growth stages. This is consistent with findings reported in maize [29] and soybean [18,19,66] significantly increased under water deficit at different growth stages. This is consistent with findings reported in maize [29] and soybean [18,19,66]. Compared with the value under V2F2G2, EC increased by 47.4% and 47.3% under V1F0G1 in 2018 and 2019, respectively.  Table 5. Water deficit at the flowering stage did not affect GPSG, except in V1F0G2 and V1F0G1; however, water deficit at the vegetative and grain-filling stages significantly decreased GPSG. Conversely, SLSG, RLSG and GISG decreased with the intensity of water deficit at the vegetative and grain-filling stages but increased with the intensity of water deficit at the flowering stage, particularly under V2F0G2.   Table 5. Water deficit at the flowering stage did not affect GP SG , except in V1F0G2 and V1F0G1; however, water deficit at the vegetative and grain-filling stages significantly decreased GP SG . Conversely, SL SG , RL SG and GI SG decreased with the intensity of water deficit at the vegetative and grain-filling stages but increased with the intensity of water deficit at the flowering stage, particularly under V2F0G2. There was no significant change in GP SG under water deficit at the flowering stage, except under V1F0G2 and V1F0G1. However, Takele and Farrant [28] found that water deficit at the flowering stage decreased GP SG by 80%. Lieffering et al. [67] reported that the rate and degree of seed imbibition were significantly related to the colloidal properties of seeds. Since proteins are an important colloidal constituent, the germination rate increased with increase in SPC [68]. Moreover, seeds with higher SPC showed a faster rate of transfer of dry matter and N reserves from the endosperm to the embryo to enhance seedling vigor [69][70][71]. On the other hand, larger seeds tend to produce more vigorous seedlings than smaller seeds [72][73][74]. Therefore, increases in SPC and KW caused by water deficit at the flowering stage may have increased GI SG , SL SG and RL SG in both years.
Deterioration using AA test is associated with emergence and storage potential. The GP, GI, SL, and RL of AA test under different irrigation treatments in 2018 and 2019 are provided in Table 6. The results showed that water deficit at vegetative and grain-filling stages significantly decreased GP AA , GI AA , SL AA , and RL AA . This is consistent with findings reported in previous researches [28,75,76]. However, water deficit at the flowering stage significantly increased GP AA in both years. Pedrinho et al. [77] suggested that GP AA was positively correlated with SPC, and in this study, water deficit at the flowering stage increased GP AA by increasing the SPC.
The activities (U mg −1 protein) of SOD, CAT, and POD of hybrid maize seeds under different irrigation treatments are listed in Table 7. The results showed that water deficit at the vegetative and grain-filling stages significantly decreased SOD, CAT, and POD. Compared with its effects on SOD activity, water deficit at the flowering stage showed limited effects on CAT and POD activities. Anwar et al. [63] found that the maximum increase in CAT activity in maize grains was observed under drought stress at the grain-filling and silking stages. Eisvand et al. [78] suggested that late drought decreased CAT activity but did not affect POD and SOD activities in wheat grains. Water deficit at the vegetative stage significantly reduced antioxidant enzyme activities, possibly due to the inhibition of enzyme synthesis or change in the assembly of enzyme subunits under water stress [78]. Table 6. Germination percentage (GP), germination index (GI), seedling length (SL), and root length (RL) on the 7th day of accelerated aging (AA) test under different irrigation treatments and weight determined using the entropy weight method (W X ) in 2018 and 2019. The weight of seed vigor indices determined using the entropy weight method ranged from 0.056 to 0.161 (Figure 2 and Tables 5-7) in both years. The weighted grey relational grade determined by GRA (R i * ), seed vigor calculated by PCA (P i * ) and CEM (V i ) are presented in   1 V, vegetative; F, flowering; G, grain-filling; 2, sufficient irrigation; 1, 50% sufficient irrigation; 0, no irrigation. 2 Numbers in parentheses indicate weights calculated by CEM, and described previously [36]. 3 R i * is weighted grey relational grade, and calculated by Equation (10). 4 P i * is the seed vigor calculated by PCA, and described previously [31]. 5 V i is the seed vigor calculated by a combinational value of each method, as was described by Jiang et al. [36].
Data analysis (Table 8) revealed that the highest seed vigor calculated by CEM was observed under treatments with V2F1G2 in 2018 and V2F0G2 in 2019, but decrease in yield by 6.3% in 2018 and 38.0% in 2019. The results showed that seed vigor and yield of hybrid maize cannot be improved at the same time under water deficit. Therefore, to ensure high yield and seed vigor of hybrid maize under drip irrigation, the optimal irrigation depth at different growth stages must be investigated in detail.

Water-Seed Vigor Model
The crop water production function is often used to simulate crop yield based on the linear relationship between relative crop yield and relative evapotranspiration [36,48]. Figure 3b shows a linear relationship between relative seed vigor of hybrid maize evaluated by CEM (V i /V CK ) and relative seasonal ET (ET i /ET CK ), which is similar to relative yield in response to ET i /ET CK (Figure 3a). Therefore, in this study, we used the crop water production functions to simulate the vigor of hybrid maize seeds under water deficit at different growth stages.
The crop water production function is often used to simulate crop yield based on the linear relationship between relative crop yield and relative evapotranspiration [36,48]. Figure 3b shows a linear relationship between relative seed vigor of hybrid maize evaluated by CEM (Vi/VCK) and relative seasonal ET (ETi/ETCK), which is similar to relative yield in response to ETi/ETCK (Figure 3a). Therefore, in this study, we used the crop water production functions to simulate the vigor of hybrid maize seeds under water deficit at different growth stages. The calibration results of the water-seed vigor models are summarized in Table 9. The calibration results of the water-seed vigor models are summarized in Table 9. For different models, R 2 and RRMSE were 0.65−0.88 and 0.132−0.175, respectively. The Rao model showed the best calibration results of R 2 = 0.88 and RRMSE = 0.132, and the water sensitivity indices at the vegetative, flowering, and grain-filling stages were 0.772, −0.343, and 0.916, respectively. Table 9. The calibration results of seed vigor water sensitivity indices at the vegetative, flowering, and grain-filling stages for different models, calibrated using 12 irrigation treatments in 2018 and 2019. R 2 , determination coefficients; RRMSE, relative root-mean-square error. The results of the water-seed vigor models validated using data of 20 irrigation treatments are presented in Figure 4. The results showed that althought the Rao model underestimated the measured value by 3%, it had the highest R 2 , EF, and d IA , and the lowest RRMSE and ARE. Therefore, the Rao model is recommended for predicting seed vigor of hybrid maize under water deficit with drip irrigation under film mulching in arid regions.

Model
Because seed vigor is a complex trait, it is imperative to establish a model to predict seed vigor under water deficit. Compared with additive models, the multiplicative models are more physiologically meaningful. In this study, the Rao model is recommended for predicting seed vigor of hybrid maize under water deficit. The Rao model showed that water deficit at the grain-filling and vegetative stages reduced seed vigor, and seed vigor was more sensitive to water deficit at the grain-filling (C 3 = 0.916) than at the vegetative stage (C 1 = 0.772), whereas water deficit at the flowering stage increased seed vigor (C 2 = −0.343). Based on the results of previous studies, seed vigor is positively correlated with KW and SPC [70,72,76,77]. In this study, water deficit at the flowering stage increased KW and SPC (Table 4), thereby improving seed vigor. On the contrary, water deficit at the vegetative and grain-filling stages decreased of KW and SPC, thereby reducing seed vigor.
The water-seed vigor model can be used to optimize irrigation water management to ensure seed vigor of hybrid maize production under limited water availability. The results of the water-seed vigor models validated using data of 20 irrigation treatments are presented in Figure 4. The results showed that althought the Rao model underestimated the measured value by 3%, it had the highest R 2 , EF, and dIA, and the lowest RRMSE and ARE. Therefore, the Rao model is recommended for predicting seed vigor of hybrid maize under water deficit with drip irrigation under film mulching in arid regions. Because seed vigor is a complex trait, it is imperative to establish a model to predict seed vigor under water deficit. Compared with additive models, the multiplicative models are more physiologically meaningful. In this study, the Rao model is recommended for predicting seed vigor of hybrid maize under water deficit. The Rao model showed that water deficit at the grain-filling and vegetative stages reduced seed vigor, and seed vigor was more sensitive to water deficit at the grainfilling (C3 = 0.916) than at the vegetative stage (C1 = 0.772), whereas water deficit at the flowering stage

Conclusions
The present study assessed the vigor of hybrid maize seeds that were produced under water deficit at different growth stages. Water deficit at different growth stages significantly decreased yield. Water deficit at the flowering stage increased SL, RL, and GI under different germination conditions. The seed vigor ranking of combinational evaluation method showed that V2F1G2 (sufficient irrigation at the vegetative and grain-filling stages and 50% of sufficient irrigation at the flowering stage) and V2F0G2 (sufficient irrigation at the vegetative and grain-filling stages and no irrigation at the flowering stage) produced higher seed vigor and reduced total irrigation volume by 17.7% and 19.3% in 2018 and 2019, respectively. The Rao model was superior to other models in predicting the vigor of hybrid maize seeds under water deficit. The results showed that an appropriate water deficit at the flowering stage could ensure seed vigor of hybrid maize production under limited water availability. Our research created a foundation for a water-saving irrigation strategy for hybrid maize production with drip irrigation under film mulching in an arid area.