Phenotypic Traits, Grain Yield and Yield Components of Maize Cultivars Under Combinations of Management Practices in Semi-arid Conditions of Iran

In semi-arid regions, selecting cultivars and planning management practices are critical issues for improving yields and reducing risks of maize cultivation in the short summer cropping season. This study aimed to assess the responses of grain yield (GY) and important phenotypic characteristics of maize cultivars from different maturity groups under various irrigation regimes, planting dates, and nitrogen rates, in a 2-year experiment. According to the clusters identified in the loading plot, stover yield, radiation use efficiency (RUE), height, leaf greenness index, stem weight, and ear size during flowering (VT/R1) were strongly correlated with GY, yield components and harvest index. Based on analysis of variance, more irrigation or N, or their interaction often increased GY, rows ear−1, grains row−1 and hundred-grain weight. Late planting decreased GY of KSC704 (late maturity cultivar). KSC260 (early maturity cultivar) had greater flexibility in planting time, required less days for maturity, and had less water consumption. The findings highlight the physiological basis of the relationship between the different phenotypic characteristics and how they affect GY and its components. They were in line with the established theories that higher RUE, biomass and sink activity (e.g. grain weight and numbers, and larger ear size) are associated with better genetic gains to produce high GY. Although the results were not consistent between years, findings suggested the good performance of recently released early maturity cultivars for use during the summer growing season.


Introduction
Semi-arid areas are characterized by low precipitation and high spatial and temporal variability leading to negative impacts of severe temperature and drought on crop production (Murungweni et al., 2016). Iran is a country in the midlatitude belt of the arid and semi-arid zones (Modarres et al., 2007). Maize in Iran is one of the most important crops that is widely cultivated as a second and summer crop after the winter cereals like wheat and barley (Rafiee & Kalhor, 2016). During the summer season, stress impacts on crop performance can be more common and persistent. Strong effects of genotype-by-environment interaction, extreme temperatures, and the short cropping time are among the most important factors that impair the development and functioning of morpho-physiological traits and grain yield of maize crops during summer (Adnan et al., 2020;Moradi et al., 2013). In addition, challenges of inadequate water and low fertilizer use efficiency because of volatilization or leaching are other important reasons for maize yield loss in these areas (Moradi et al., 2013(Moradi et al., , 2014Nasielski et al., 2019).
There are several strategies recommended and adopted by breeders and farmers to reduce the yield gap in the arid environmental zone. Optimal combinations of sowing window and cultivar are among the most important strategies for minimizing the impact of drought and heat during summer on anthesis and autumn cold on grain filling in that environment (Chen et al., 2020;Padovan et al., 2020).
For summer cropping in semi-arid parts of Iran, while early planting is better for yield maximization by taking advantage of spring rainfalls for better germination and crop establishment, however, farmers are advised to postpone cultivation to prevent synchronizing of the maize reproductive stage with extreme temperatures (Moradi et al., 2014). In addition, farmers are recommended to use early maturity maize cultivars (Moradi et al., 2013(Moradi et al., , 2014, so that they can reach maturity before the onset of autumn cold. Additionally, early maturity cultivars consume less water, due to a shorter maturation cycle (Jahangirlou, 2015). However, early maturing cultivars are not well received by farmers due to their high yield loss under water and N stresses, compared with late maturing ones (Jahangirlou, 2015).
In recent years, grain yield and agronomic aspects of new maize cultivars have been improved due to breeding efforts. Screening the phenotypic characteristics, grain yield and yield components of the best performing cultivars under multi-environment trials and stress conditions is an important strategy for investigating their behavior to genetic by environment interactions.
Grain yield, as a complex variable, can reflect the interaction of the environment and management with the growth and development processes that occur throughout the crop's maturation cycle (Studnicki et al., 2016). In addition, yield component traits adjust their expressions to determine grain yield under different environmental and agronomic conditions (Kumar et al., 2016(Kumar et al., , 2017. Since the reduction in grain yield and yield components of cultivars under stress conditions is often correlated with changes of some phenotypic expressions (Kumar et al., 2017), the evaluation of phenotypic traits are more important. As an example, variation in intercepted photosynthetically active radiation (IPAR) due to changes in location, year and planting date is the main cause of variation in the phenology, kernel set and harvest index (HI) of maize crop (Tsimba et al., 2013). Grain yield reductions due to drought stress at different growth stages are well documented, and the causes are associated with reduced tissue expansion and biomass, cell turgor, plant height elongation, ear elongation, number of fertilized ovaries, and reduced kernel number (Ruiz et al., 2019). The reduction of grain yield can also be associated with the reduction in IPAR due to changed phenology, LAI, and stress intensity (Sarlikioti et al., 2011;Shi et al., 2015;Villegas et al., 2016). In contrast, sufficient water and nitrogen (N) rates increase dry matter accumulation in maize grain due to enhancing leaf area and chlorophyll content, absorption of solar radiation, and biomass partitioning to individual organs, especially stems as reservoirs to remobilize assimilates to grain (Boomsma et al., 2010;Zhou et al., 2018).
Breeders have developed high yielding early maturity maize cultivars in the recent years, but understanding their responses to diverse systems of cultivation is limited. Identifying correlated phenotypic characteristics with grain yield and their response to challenging agronomic factors or interactions is key for developing strategies required to enhance and stabilize grain yields. Therefore, this study aimed to (i) evaluate the responses of locally popular maize cultivars from different maturity classes to combinations of planting times, irrigation, and N regimes to gain a quantitative understanding of yield formation, and (ii) identify important phenotypic characteristics affecting grain yield and its components. A motivation of the study was to provide useful knowledge, more efficiently use money and time defining and conducting future research, and develop guidelines to enhance maize grain yield and reduce production risks in Iran.

Experiment Site
The field experiments were performed in 2018 and 2019 in Pakdasht city (35.4669° N, 51.6861° E), Tehran province, Iran. The geographical location of Pakdasht city in Iran is shown in Fig. 1. According to the de Martonne index climate classification, this region has a semi-arid climate with relatively cold winters and hot dry summers (Emadodin et al., 2019). Irrigated summer maize in this region is often cultivated after winter cereals like wheat and barley. Summaries of minimum and maximum average temperatures and precipitation during the summer maize growing season for the two experimental years are presented in Fig. 2. The highest average temperature occurred in July and the total growing season precipitation was less than 1.0 mm in each year.
Soil sampling was done before the first cultivation year. The samples were air dried, passed through a 2-mm sieve, and tested for organic carbon using the Walkley-Black method (Walkley & Black, 1934), total N using the Kjeldahl method (Bremner, 1960), available phosphorus (P) by the Olsen procedure (Olsen et al., 1954) and available Potassium (K) using a flame photometer (Mehlich, 1953). Based on particle size analysis, the texture was characterized as clay-loam. The initial physico-chemical characteristics in the 0−30 cm soil layer were as follows: organic carbon, 1.09%; total N, 0.12%; P, 90.17 mg kg −1 ; K, 453 mg kg −1 .

Experimental Design
In mid-June 2018 experimental strips and plots were designated at a farm where no crops had been cultivated for the previous 3 years. The field was divided into two strips (24-m length) to apply different irrigation timings (12-day and 6-day intervals). Two planting dates (20 June and 21 July) were randomized within the irrigation treatments. Two ditches crossed the strips, partitioning each of them into three plots. Thus, the experiment comprised twelve main plots. The four combinations of cultivars (KSC704 and KSC260) and N rates (0 and 184 kg ha −1 ) were randomized to subplots within each plot. Thus, a total of 12 main plots containing a total of 48 sub-plots were obtained from three replications in each cultivation year. Individual sub-plots were 6 m in length and consisted of 6 rows sown at a density of six and eight plants per square meter for KSC704 and KSC260, respectively. Treatments were applied to the same plots in each year.
KSC704 is a locally popular high-yield late maturity hybrid which has been cultivated in Iran since 1980 and accounts for 80% of the maize growing areas of Iran (Estakhr et al., 2015). KSC260, an early maturity maize cultivar, was introduced in 2008, and has shown good performance in various experiments comparing new early-maturing cultivars (Jahangirlou, 2015). Harvesting date, maturity period and the total growing degree days (GDDs) of cultivars in present study are shown in Table 2. GDD was calculated using GDD calculator program (Soufizadeh, 2012). For this study, 10 °C as the base temperature and 34 and 40 °C as the optimum and maximum temperature thresholds, respectively, were adapted from Cutforth and Shaykewich, (1990). The two planting dates were chosen to create different temperature conditions during grain filling. The variability of daily average temperature during the reproductive period (listed in Table 1) for both planting dates in 2018 and 2019 is shown in Fig. 3. There were clear differences between the recorded temperatures for the planting dates in each year, and also differences between the 2 years. The second year was warmer year during the reproductive period. An irrigation regime with 6-day intervals and N rate of 184 kg ha −1 were specified as typical non-stressed growing conditions. Furrow irrigation was used for irrigation. The closed-end furrows were constructed with a ditcher. The ridge and furrow widths were 40 cm and 20 cm respectively, and the depth of the furrow was 15 cm. Based on the flow rate per furrow (∼1.3 L/s.), water was delivered until the furrows were completely full and the duration of each irrigation event was recorded. The known flow rate of the irrigation pump was then used to estimate irrigation volumes. In total, the KSC704 cultivar received about 9200 m 3 ha −1 of water with 6-day irrigation timing and about 5800 m 3 ha −1 with 12-day irrigation timing. KSC260 received less water due to a shorter growing period and less irrigation events: 7800 m 3 ha −1 with 6-day irrigation timing and 4800 m 3 ha −1 with 12-day irrigation timing. The Food and Agriculture Organization (FAO) estimates the water requirement of latematuring maize cultivars in arid and semi-arid regions to be around 7000 to 8000 m 3 ha −1 (Faostat, 2017). Urea was used as the N source, as the common source of N fertilizer used by farmers and top-dressed in equal proportions at two stages: pre-planting and V4-V6 (4-6 leaves with visible collars) stage of maize growth. Fresh irrigation water was used from the main ditch for each replication, and excess water was not re-used.

Sampling and Measurements
In both years of cultivation, a set of phenological and morpho-physiological variables were measured. Approximately every week, ten marked plants of each subplot were checked to determine important phenological stages including days to germination, tasseling, silking, physiological maturity, and grain filling period. Destructive sampling was conducted for six and eight stages, including VT (tasseling) and R1 (silking), in the first and second year, respectively. Each sampling time, 4 maize plants were randomly selected and harvested from each subplot. After measuring plant height and leaf area index (LAI), the organs were separated and dried in a laboratory oven at 60 °C until a constant weight (∼48 h) and subsequently leaf and stem weight were measured. The trend of changes in the desired characteristics was measured using destructive sampling and their maximum values, which often occurred in the VT and R1 stage, were recorded. A leaf area   meter (LI-COR LI-3100C model, USA) was used to measure LAI. At the R6 (physiological maturity) stage 15 plants in an area of ∼11 m 2 (4-m row length × 0.7-m row spacing × 4 rows) were cut from each subplot to measure stover yield, grain yield and harvest index (HI). The organs were separated and the following were measured: weights of organs, ear length, barren length, ear and cob diameter, rows ear −1 , grains row −1 , grains ear −1 , and hundred grain weight. After appearance of buds on the stem that will form ears, a Minolta SPAD-502chlorophyll meter (Minolta Camera Co., Ltd, Japan) was used for estimation of greenness index of the main ear leaf at R1 and R4 (dough) stages. Greenness index measurements were made on undamaged ear leaves of each of the 10 marked plants in each sub plot at the same time close to solar noon. To assess and characterize the ability of the varieties to photosynthetically utilize radiation for biomass production under different growing conditions, radiation use efficiency (RUE) was calculated (Dingkuhn et al., 2015;Zhang et al., 2009). Similar to Dingkuhn et al., (2015), the photosynthetically active radiation (PARi) was calculated from the dynamics of LAI by fitting two or three regression models that had the best coefficient of determination (R 2 ) to interpolate daily LAI values from planting to maturity (supplementary 1, Fig. S1 and Table S1). The fraction of light intercepted was estimated using Lambert-Beer's law (Monsi & Saeki, 1953), with PARi = PAR (1 -e −kLAI ); where k is the extinction coefficient. The extinction coefficient was set at 0.49 for KSC704 and 0.52 for KSC260, adopted from the study of Jahangirlou (2015). Daily global solar radiation was converted into PAR by multiplying by 0.48 (Dingkuhn et al., 2015); then PARi was cumulated from sowing to maturity and used to estimate RUE from the final above ground dry weight, separately for each plot.

Statistical Analysis
Principal component scores were derived using Minitab, v.19 (Minitab Inc., 2006) after standardizing the variables by using the correlation matrix. Results from the two cultivation years were analyzed separately, using mixed model procedures with PROC MIXED (SAS Institute, 2003). The mixedeffects model included fixed effects of irrigation, planting date, cultivar, N, and their interactions, and the random effects of irrigation*planting date and irrigation*planting date*block interactions. Because of the experimental design, irrigation effects have to be interpreted with caution-possible meaningful effects of irrigation and its interactions can be either caused by differences between the planned irrigation treatments, or by differences in the two sections of the field where the irrigation treatments were applied. Tukey's statistic was used to test differences (P ≤ 0.05) among means.

Correlations Between Yield and its Component Variables and Agronomic Treatments
The PCA comprising the first two principal components accounted for 51.5% of the total variance. In the loading plot (Fig. 4A) are the eigenvectors of characteristics for 31 variables related to growth and development of maize during the growing season. The maximum height (Hm), tassel length (TaLm), stem weight (SWm), leaf weight (LWm), tassel weight (TWm), and husk weight (HWm), and also maximum leaf area index (LAIm) shown in the loading plot belong to tasseling or silking stages (VT/R1).
Three informal clusters of correlated variables were made according to their distribution (similar eigenvector values and low angle between eigenvectors) in principal components 1 and 2. In the first cluster, GY is highly correlated to stover yield (SY), ear weight (EW), and RUE. The correlation of these variables is confirmed by the correlation analysis (supplementary 2, Table S2). For instance, the correlation between GY and RUE was 0.824 (Table S4). In the second cluster, number of grains row −1 (GR) is highly correlated to Hm and leaf greenness index (LGI) at silking (LGIs). The correlation between GR and LGIs was 0.508. The eigenvector for grains ear −1 (GE) points in the opposite direction; thus, the three clustered variables were negatively correlated to GE. For example, the correlation between GE and GR was − 0.652. In the third cluster, HI, and hundred-grain weight (HGW) were highly correlated to SWm, stem weight (SW), leaf weight (LW), ear diameter (ED), ear length (EL), and number of rows ear −1 (RE) at physiological maturity. As an example, the correlation between SWm and ED is 0.582.
The score plot (Fig. 4B) plots 2 years, two planting dates, two cultivars, two irrigation regimes, and two N rates, giving a total of 32 points, averaged across replicates. There was a noticeable pattern of different levels of irrigation and the years; so means in Fig. 4B were marked by these two factors, to highlight the effects. The first component approximately separates the high irrigation rate from the low irrigation rate-markers for the low irrigation rate, are often located in the left of the score plot, whereas markers for the high irrigation rate are often in the right. The second component clearly separates the first year of the experiment from the second year-markers for the first year (2018), are distributed toward the top of the score plot, whereas markers for the second year (2019) are located toward the bottom. The data points presented in score plots coincide with the direction of changes in variables in the loading plots. The variables from the first cluster and second cluster are associated most strongly with more frequent irrigation, and the year 2018. This means that plots with more frequent irrigation during 2018 are likely to have higher values for these measured variables. The variables from the third cluster are associated most strongly with more frequent irrigation, and the year 2019. The three clusters identified are all associated with more frequent irrigation. Also, LAI and GE are strongly associated with the 2019 season.

Treatment Effects on Grain Yield and its Components
The effects of treatment main effects on yield components are presented in Table 2 and significant interactions are plotted in Fig. 5 (unless no pairs of means were significantly different or the graph could not be clearly interpreted). P-values of all main effects and interactions for yield components are in supplementary material 3, Table S3. P-values and results of analysis of variance for vegetative characteristics are in supplementary material 4. In 2018, high N increased GY by one tonne ha −1 ( Table 2). In addition, KSC260 had higher GY than KSC704 by 0.96 tonnes ha −1 (Table 2). In 2019, there were interaction effects of irrigation and N rate (Fig. 5A) and planting date and cultivar (Fig. 5B) on GY. Irrigation rate had the greatest effect, and for either irrigation rate, the high N rate increased GY (Fig. 5A). KSC704 with an early planting date had the highest GY (10.84 tonnes ha −1 ), whereas for KSC260 there was no planting date effect.
In 2018, EW with an early planting date was higher than with a late planting date by 4.57 tonnes ha −1 (Table 2). In 2018, KSC260 had a higher EW than KSC704 by 1.78 tonnes ha −1 (Table 2). In 2019, there were interaction effects of irrigation and N rate (Fig. 5C), irrigation and cultivar (Fig. 5D), and N rate and cultivar (Fig. 5E) on EW. For any irrigation rate, the high N rate increased EW. In addition, EW was higher with the high irrigation rate, regardless of the N rate (Fig. 5C). For any cultivar, the high irrigation rate increased EW (Fig. 5D). For both cultivars, the high N rate increased EW (Fig. 5E). With high or low irrigation and N shows the eigenvector of each characteristic. The score plot B shows the means, grouped by year and irrigation treatments (I1, Irrigation at 12-day intervals; I2, Irrigation at 6 day-intervals). Each marker represents the mean of replicates in the field. Abbreviations for eigenvectors: BL barren length, CD cob diameter, CW cob weight, DG days to germination, DPhM days to physiological maturity, DS days to silking, DT days to tasseling, ED ear diameter, EW ear weight, EL Ear length, GE grains ear −1 , GFP grain filling period, GR grains row −1 , GY grain yield, HW husk weight, HWm maximum husk weight, HGW hundred-grain weight, HI harvest index, Hm maximum height, LAIm maximum leaf area index, LGId leaf greenness index at dough, LGIs leaf greenness index at silking, LW Leaf weight, LWm maximum leaf weight, RE rows ear −1 , RUE radiation use efficiency, SW stem weight, SWm maximum stem weight, SY stover yield, TaLm maximum tassel length, TWm maximum tassel weight, TW tassel weight. The variables including BL, CD, CW, ED, EW, EL, GE, GR, GY, HW, HW, HI, LW, SW, SB and TW were measured at physiological maturity rates, there was no significant difference between cultivars in terms of EW. In 2018, KSC260 had higher HI than KSC704 by 4.03% (Table 2). In 2019, high N increased HI by 2.27% (Table 2). In 2018, ED with an early planting date was higher than with a late planting date by 0.26 cm (Table 2). In 2019, there were main effects of irrigation, cultivar and N rate on ED (Table 2). High irrigation and N rate increased ear diameter by 0.35 and 0.15 cm, respectively. Cultivar KSC704 had higher ear diameter than KSC260 by 0.17 cm (Table 2). In 2018, there were interaction effects of irrigation and N rate (Fig. 5F) and cultivar and N rate (Fig. 5G) on ED. With the high irrigation rate, the high N rate increased ED (4.62 cm), whereas with low irrigation there was no N effect. In KSC704, the high N rate increased ED (4.65 cm), whereas in KSC260 there was no N effect. In 2019, high irrigation increased EL by 1.68 cm (Table 2). In 2018, there was an interactive effect of cultivar and N rate (Fig. 5H) on EL. In KSC704, the high N rate increased EL (4.65 cm), whereas in KSC260 there was no N effect. Cultivar KSC704 had higher EL than KSC260, regardless of the N rate.
In 2018, KSC704 had higher RE than KSC260 (Table 2). In addition in 2018, the high N rate and in 2019 the high irrigation rate increased RE (Table 2). In both years, high N and irrigation rate increased GR (Table 2). In 2018, there was an interactive effect of planting date and cultivar on GR (Fig. 5I). For KSC704, planting at early date increased GR, whereas for KSC260 there was no planting date effect. For KSC704, GR was higher than for KSC260, regardless of the planting date. High N rate increased HGW by 2.23 and 2.39 g in 2018 and 2019, respectively. In 2019, high irrigation rate increased HGW by 4.06 g. In 2019, KSC704 had higher HGW than KSC260 by 1.5 g. In 2018, there were interactive effects of cultivar and irrigation (Fig. 5J) and cultivar and planting date (Fig. 5K) on HGW. In KSC704, high irrigation increased HGW, whereas in KSC260 there was no irrigation effect. In both cultivars, HGW with an early planting date was higher than with a late planting date.

Physiological Basis of Correlations Between Yield and its Component Variables and AGRONOMIC Treatments
The results showed that grain yield (GY) is highly correlated with stover yield, ear weight and radiation use efficiency (RUE). This is in line with the established theory that high biomass is associated with better genetic gains to produce high grain yield (Reynolds et al., 2017). However, in the areas faced with drought stress, increasing biomass (especially due to mass and area of leaves) does not always related to achieving higher grain yield, because there is strongly negative correlation between leaf water loss, transpiration and leaf size across the species (Wang et al., 2019). RUE is highly correlated with biomass production (Soufizadeh et al., 2018) and chlorophyll in the leaves, an indicator of N content which improves assimilate synthesis, biomass production, grain filling, and yield (Liu et al., 2018).
The correlations between GY and RUE and grains row −1 (GR) and leaf greenness index (LGI) in the study indicate that results are in line with earlier studies. In the present study, LGI was measured in two stages (silking and dough), but LGI of ear leaves at silking had a higher correlation with grain yield and its components than LGI at dough stage. This Least squares means of maize grain yield (GY) in response to interaction effects of irrigation and nitrogen (A) and planting date and cultivar (B); ear weight in response to interaction effects of irrigation and nitrogen (C), cultivar and irrigation (D) and cultivar and nitrogen (E); ear diameter in response to interaction effects of irrigation and nitrogen (F) and cultivar and nitrogen (G); ear length in response to interaction effects of cultivar and nitrogen (H); number of grains row −1 in response to interaction effects of cultivar and planting date (I); 100-grain weight in response to interaction effects of cultivar and irrigation (J), cultivar and planting date (K). These graphs are only for significant interactions. I1 irrigation at 12-day intervals, I2 irrigation at 6 day-intervals, PD1 planting date June 21, PD2 planting date July 22, C1 cultivar KSC704, C2 cultivar KSC260, N1 zero-nitrogen, N2 184 kg N ha −1 . Least squares means labelled with the same letter do not differ significantly at P < 0.05, based on Tukey's test. Vertical bars represent the 95% confidence interval is because tasseling and silking are extremely critical stages for pollination, kernel formation, and increasing the number of filled grains, and N and chlorophyll status of the plant at this stage has a significant effect on grain yield (Gaju et al., 2016). Additionally, the carbohydrate and N in the leaves at this stage are remobilized towards developing grains (Ning et al., 2018) in the grain development phase. As the crop moves closer to the maturity, the supply of carbohydrate and N for grain development becomes limiting, because the crop starts shedding the lower leaves and supplies the resources to the grain (Ning et al., 2018). Hence, the higher correlation of greenness index (an indicator of chlorophyll and nitrogen content) at tasseling than at dough stage is in accordance with the physiological mechanisms. Similarly, the correlations between harvest index (HI), and grain yield components (hundred-grain weight (HGW), grains row −1 (GR) and ear diameter and length) with vegetative characteristics (maximum stem weight at flowering, and stem and leaves weight at physiological maturity) further indicates the source-sink balancing phenomenon by the crop. In mature grains, the biggest part of dry matter is remobilized from pre-anthesis accumulated reservoirs of organs, especially from the stems (Crevelari et al., 2018). Stem weight at flowering is important for increasing grain weight and yield; in other words, stronger sink activity (e.g. HGW and GR, and larger ear size) is one of the main reasons for achieving greater grain weight (Zhao et al., 2020).
The results of this study suggest that irrigation and year accounted for much of the data variation. Sah et al. (2020) reported recently that water stress explained up to 66% of changes in grain yield characteristics in non-drought tolerant maize lines. Two main potential sources can be considered for the differences among years herein: (a) conducting the first cultivation after 3 years fallow, and (b) differences of daily temperature during the reproductive stage among years. Keeping soil fallow for a season or more improves soil fertility, organic matter and physical properties to supply more elements needed for assimilation; whereas, continuous cultivation as performed in this study could make some elements deficient for grain growth and development. Zhu et al. (2019) reported that high temperature can affect maize grain yield greatly. More than half of the yield changes of maize in the U.S. Midwest was related to temperature variations among years during the grain filling period.

Coupling the Responses of Grain Yield and its Component with Agronomic Treatments
The results of the first year showed that high N rate increased GY. This result was not surprising because the application of N fertilizer increases the N content in leaves, which is strongly correlated with high greenness index or chlorophyll content and increasing the activities of chloroplast, thereby increasing the photosynthetic rate and grain yield . In the second year, for any irrigation rate, high N rate increased GY and EW, and they were higher in high irrigated treatments, regardless of the N rate. These results suggested that irrigation was more critical than N, and they are in line with the results obtained by Di Paolo and Rinaldi, (2008), who reported N application amplified the irrigation effect and vice versa.
Surprisingly in the first year, early maturity cultivar KSC260 had higher GY and EW than the late maturity cultivar KSC704. However, in the second year with the continued cultivation, KSC704 in early planting date had higher grain than KSC260, whereas it declined in late planting. Late maturing cultivars are high yielding because of the longer vegetative and reproductive stage, but have higher requirements for solar radiation, water, and nutrition (Huang et al., 2020). Besides the genotype factor (G), the agronomic management factor (M) of planting late maturity cultivars earlier in arid regions can take advantage of late spring rainfalls which facilitates its germination and establishment, and together with less chances to be exposed to the temperature drop in autumn at the end of the growing season, an environmental factor (E) (Wang et al., 2017), KSC704 might be exposed to conducive G × E × M in the second year, particularly for better grain development. Although higher grain yield was observed for KSC260 in the first year and also in late planting date of the second year compared to KSC704, there was no significant difference between GY of KSC704 and KSC260. Considering the results of the second year, it can be stated that KSC260 can provide more flexibility to crop managers to shift planting times and allows producers in certain production areas to harvest before the onset of inclement weather.
In addition, the reasons for higher grain yield of KSC260 in the first year can be linked with the individual genotypic traits. In KSC260, there was a high frequency of doubleear plants, whereas KSC704 plants were mostly single-ear. KSC260 had all the grains filled in the secondary ears, which was unusual as secondary ears are usually barren in maize or only can produce 5 to 10% of primary ear grains (Nafziger, 1996). However, some studies reported that the growth and development of secondary ears can be affected by environmental stress factors like temperature . This could be the potential reason for the increased frequency of barren secondary ears of KSC260 and reduced its GY, in 2019. In general, the inconsistency in GY among years and planting dates could suggest a large impact of exogenous factors like edaphic and climatic variables, in addition to the maturity group and genetic potentials on grain yield (Maresma et al., 2019), which can further be associated with the complexities of G x E x M.
In the first year, EW was higher in the early planting date than the late one. However, it seems that the increase in EW was due to the increase in the straw part, rather than the grain, because there was no significant difference between planting dates in terms of grain yield. In the second year, for any cultivars, the high irrigation and N rate increased EW, similar to Ercoli et al. (2009) who reported that higher irrigation and N regimes can increase organ biomass (e.g. ears) by increasing the rate of assimilate production.
The HI was higher in KSC260 than KSC704 in the first year, which can be associated with the higher GY and EW. However, the other yield components i.e. RE, GR, and HGW, which contribute to the yield and thus influence HI were lower than KSC704. This again suggests that the secondary ears of KSC260 and its filled grains were important factors for grain yield and HI differences among the cultivars. Considering the early maturing characteristics of KSC260 and its potential to escape from the high temperature stress during reproductive stages the responses of the cultivar are in accordance with the study of Kobata et al. (2018). In the second year, HI of KSC260 was higher than KSC704 but there was no significant difference between cultivars, which similar to the first year, can be attributed to EW of the cultivars. In addition, the HI increase with high N in the second year can be explained by the fact that N application enhances the physiological efficiency in the translocation of photosynthetic products to kernel number and weight which increases grain yield and subsequently HI (Asibi et al., 2019). In response to different irrigation and planting date treatments, HI during 2 years of cultivation was variable but not significantly different, which suggested that it was stable under those conditions. A similar response where HI was variable in different years and managements but was statistically stable, reported by Kumar et al. (2016).
The increased ED and EL in response to the main or interactive effect of high irrigation and N rate might explain water and nutrition utilization ability of the crop. Ear diameter as an important criterion to reflect crop status in terms of the interplant competition is expected to increase in response to high utilization of water and nutrition in the plant canopy, as a consequence of increasing accumulation of plant biomass, its translocation to ear biomass, and the efficiency of using this biomass for kernel set (Borrás & Vitantonio-Mazzini, 2018). However, the response of ED to N rates varied with maize cultivars. One possible reason for this discrepancy can be related to genotypic characteristics. It is reported that cultivar KSC704 potentially has a bigger ear size than earlier maturing cultivars (Rafiee & Kalhor, 2016).
For any year of cultivation, high irrigation and N rate increased GR and high N rate increased HGW. The high N rate increased RE in the first year, and high irrigation rate increased RE and HGW in the second year. These observations are expected because of high N and irrigation effects on forming and filling more ovules (potential grains) of maize ear and their effects on increasing grain cellular proliferation and growth (Echarte & Tollenaar, 2006).
In the first year, the response of GR to planting dates was variable between cultivars. In the second year, the response of HGW to irrigation and planting date was variable between cultivars. In general, KSC704 had the highest GR and HGW. Jahangirlou (2015) reported higher GY and yield components of KSC704 compared to other local cultivars. Plant maturity is usually an important factor in the development of yield in maize.
Prolonged maturity makes more days available for higher biomass and grain yield production, if the environmental conditions at the end of the growth season are favorable (Jahangirlou, 2015). The results of second year showed that HGW decreased with late planting, regardless of cultivar, which could be associated with low temperatures and high diurnal temperature range which decrease kernel growth rate and ultimately final kernel weight (Zhou et al., 2017).
The results from the analysis of variance were generally in line with the results of principal component analysis, however, some non-commonalities are expected, because the score plot only shows means grouped by irrigation and year of principle components 1 versus 2, whereas correlation matrix values take into account all dimensions (Parsons et al., 2006).

Conclusion
The findings provide insights into the phenotypic variables that are strongly correlated with grain yield and its components in accordance with the physiological mechanisms. Stover yield, RUE, height, LGI, stem weight at flowering, and ear size were strongly correlated with GY, yield components and HI. Irrigation and year had the greatest effect on maize crop characteristics. More irrigation (6-day intervals) or N rate (184 kg N ha −1 ), or their interaction significantly increased GY and yield components. KSC260 had a higher GY in the first year across the treatments and KSC704 had a higher GY only in the early planting date of the second year. The inconsistency of cultivars GY among years and planting dates can further be associated with the complexities of G × E × M. However, KSC260, the early maturity cultivar, is more recommended because of good performance, greater flexibility to be used at various planting times, a 15-20 days earlier maturing cycle, and 1000 m 3 ha −1 less water consumption than KSC704. The findings are useful for understanding the importance of nitrogen, and selecting cultivars and planning irrigation in semi-arid conditions, where the dry and short growing season are important challenges in dent maize cultivation. In addition, the reported responses of genotypic traits under varied management practices can