Evaluation of Agronomic Performance of Maize ( Zea mays L.) under a Fertilization Gradient in Transylvanian Plain

: In the last few years, Romania has become a top maize producer. Export potential is sustained by ensuring high-quantity and -quality maize. Success of maize crop is highly dependent on inputs. In this context, insight into the potential of different fertilizers to maximize crop performance could shed light on best practices to enhance yields and other traits of interest. The aim of this study was to assess the agronomic performance of maize under a fertilization gradient. Six fertilizer regimes were tested on three maize hybrids between 2018 and 2020, in conditions from the Transylvanian Plain. Results showed that fertilization had a signiﬁcant inﬂuence on yield, thousand kernels weight, grain quality (starch and protein content) and crop health. The experimental year also played a signiﬁcant role in the expression of productivity potential of maize genotypes. Different fertilizer regimes could be used for targeting desired outcomes, but top performance across all or multiple agronomic components remains a challenge and should receive further attention for optimization.


Introduction
Romania has been among the top 10 producers of maize worldwide since 2018 according to the Food and Agriculture Organization of the United Nations [1] and the most important maize grower within the European Union according to Eurostat [2]. Statistics from the Ministry of Agriculture and Rural Development from Romania indicate that, since the integration in the European Union, both the average yield-per-unit area as well as total production of maize increased at the national level [3], presenting attractive export potential. The possibility to increase exports is relying on the optimization of maize crop performance, in order to ensure high-quantity and -quality maize.
Yield of maize crop is to a high extent dependable on the application of fertilizers [4,5], and will continue to be so in the foreseeable future. Out of the nutrients required for crop growth, N (Nitrogen), as well as P (Phosphorus), and K (Potassium), are considered the most consequential for successful cultivation of maize [5][6][7][8][9]. It is standard knowledge that adequate fertilization is essential for realizing maize hybrid yield potential, but nutrient application has to be conducted both on sustainable environmental and financial considerations [10]. Yields can be obtained in a sustainable manner through balanced nutrition, which ensures adequate uptake of both macronutrients as well as micronutrients by the plants and that, in turn, ensure grains that are healthy, biofortified and nutritious [11]. Because optimal nutrient supply of the plants is related to health of the crop, it can be considered a component of Integrated Pest Management [12]. For maize crop, the pattern of nutrient uptake can vary in time, rate and duration, while also varying with each specific nutrient and its plant part where it is allocated [10]. Grain yield and nutrients accumulation in kernels are chiefly influenced by genotype and nutrients available in soil [8].
Rapid scientific advances in soil-plant interactions from recent years have generated a ripple effect for trends on fertilizers market. As a consequence, farmers today are presented with a variety of options that promise success for their crops. From mineral to organic components, besides various formulations that stimulate soil biota or plant performance, the list of options is expanding. However, clearly defined agronomic benefits of fertilization regimes following a comparative approach can highlight best practices in connection with desired outcomes, as well as areas that need further optimization. Understanding the behavior of widely used or new maize hybrids under novel fertilization regimes could help in identifying potential particularities that might affect nutrient recommendations in regards to timing and rate. Since sustainable stewardship of nutrient use in agriculture is at the nexus of successful crop-environment-economic efficiency [10,13], investigation into interaction between environment, genotypes and fertilization on agronomic outcomes becomes an indispensable approach for food security.
The aim of this research was to provide an evaluation of maize performance under a fertilization gradient, by screening a range fertilizer regime. The results are meant to provide insight into the performance potential of newly trended options to enhance productivity components, compared to the traditional basal complex fertilizer + calcium ammonium nitrate, a fertilization regime well established in agronomic practice for the past few decades. The objective of the research was to evaluate the influence of experimental factors and gradients on the following: plant growth, maize ear disease and pest's incidence, on plant quantitative and qualitative productivity components.

Location and Pedo-Climatic Conditions
The experiment was conducted at Agricultural Research and Development Station (ARDS) fields from Turda, Romania, situated at latitude 46.35, longitude 23.47 and altitude 345-493 m [14]. The ARDS Turda is situated in the Transylvanian Plain, characterized by hilly landscape [15,16]. Between hills, there are valleys oriented mainly on the direction east-west, with naturally poor drainage [17,18]. The characteristic landforms of this area are gentle and moderate slopes subject to erosion processes [14]. According to ARDS Turda, arable farming land presents a gentle sloping gradient (3-5%).
The climate is continental temperate with four seasons. According to Köppen-Geiger classification, Turda has a humid continental climate [19].
Weather data from the experimental years are presented in Figure 1, alongside with average values for the last 60 years, as registered at the weather station of ARDS Turda. The soil type is Phaeozem and having the following average values of agrochemical indices determined at depth of 0-30 cm: pH 7.00, humus 3.40%, Nitrogen 0.226%, Phosphorus 73 ppm, Potassium 295 ppm [20]. During the entire maize growing season (April-September), the average monthly temperatures in the year 2018 exceeded the 60-year average. Also, during the months of April until August, the average monthly temperatures were highest in 2018 compared to the other two experimental years.

Organization of the Experiment
Regarding precipitations, one can observe that the year 2019 experienced the highest precipitations during the spring months out of the three experimental years, as well as exceeding the multiannual average, but was lowest in precipitations during the months of June-July. By comparison, the year 2020 experienced lower precipitations than the multiannual average or the other two previous experimental years in the months of April-May, but afterwards, the precipitations were more abundant, reaching the highest levels in June-July.

Organization of the Experiment
The experiment was located in the field at ARDS Turda and had the following factors:

•
Factor Y (the experimental years) with three levels: The fertilization products were obtained from the companies Naturevo [21] and Timac Agro [22] from Romania. The description of the products according to the label is provided in Table 1. The products were chosen in order to cover a diverse range: both the classic, frequently used products as well as novel fertilization products. The rate of application of the products is mentioned in Table 1. The fertilization regime was comprised by basal complex fertilization and calcium ammonium nitrate application has been most commonly used in Romania by farmers and selected as control for the comparative analysis of productivity components. Out of the six fertilizer regimes tested, F6 was the only one that did not involve also phasic application.
Two of the maize hybrids studied were created at ARDS Turda (H1, H2) [24], while the third one (H3) was created by Pioneer [25]. All three hybrids were created to be cultivated for grains and belong to the same FAO maturity group (380). PR37N01 is the most cultivated maize hybrid from Pioneer in the Transylvanian Plain, while the other two maize hybrids (Turda 332; Turda 344) are among recent creations of the ARDS Turda.
Descriptions of the maize hybrids according to the breeders are the following: • "Turda 332" has good resistance to drought and shriveling of kernels. Plants are tall with 16-17 leaves. The ear is long with a cylindrical shape, having 18-22 rows and red rachis. Kernels are dentate with a dark yellow color. Thousand kernels weight is between 220 and 240 g [24]. • "Turda 344" is characterized by high productivity potential of 14 t/ha, good resistance to low temperatures during first growth stages. Plants are tall, reaching 3.30 m in height, and ear insertion is between 1.20 and 1.30 m. It is stay-green hybrid. The ear has 18-29 rows and thousand kernels weight of 300-320 g. Kernels are dentate [24]. • "PR37N01" is a hybrid with stable yields and resistant to drought, lodging and main maize crop pests and diseases. It has balanced height and medium-high ear insertion. The ear has 16-18 rows and 37-41 kernels per row. It is a stay-green maize hybrid [26].
The maize crop was cultivated under a conventional system. Sowing was performed in May, and maize was harvested in October. Plant density was 57,000/ha. No herbicides, pesticides or any disease control agents were used. Weed control was performed only mechanically.
The experimental variants were organized according to a randomized block design ( Figure 2) in field conditions. There were three repetitions in the field, and the experiment took place during the span of three years (2018-2020). Each variant in each repetition had four rows of plants, and all determinations were conducted on the plants from the two rows in the middle. The experimental variants were organized according to a randomized block design ( Figure 2) in field conditions. There were three repetitions in the field, and the experiment took place during the span of three years (2018-2020). Each variant in each repetition had four rows of plants, and all determinations were conducted on the plants from the two rows in the middle.
Plant characteristics were determined by measurements in the field, for ten plants from each variant in the month of August when plants were in full maturity, after the browning of the silk. After harvesting the maize in October, the hectoliter mass was determined using Granomat equipment [27], while grain quality parameters were determined with NIR Tango spectrophotometer [28]. Yield was determined following manual harvest and weighing, followed by correction at STAS humidity of 15%.
Disease incidence was calculated as frequency (1) for maize ears according to standard methodology [29]: incidence = infected or infested maize ears total number of ears analyzed × 100 (1)

Statistical Analysis
Data were analyzed for each of the 12 agronomic components by applying three-way ANOVA and using Origin software by OriginLab [30]. Significant influence of bipartite interaction between factors over a given agronomic component was further explored using POLIFACT (non-commercial software) by applying post hoc tests: either Duncan test
Plant characteristics were determined by measurements in the field, for ten plants from each variant in the month of August when plants were in full maturity, after the browning of the silk. After harvesting the maize in October, the hectoliter mass was determined using Granomat equipment [27], while grain quality parameters were determined with NIR Tango spectrophotometer [28]. Yield was determined following manual harvest and weighing, followed by correction at STAS humidity of 15%.
Disease incidence was calculated as frequency (1) for maize ears according to standard methodology [29]: incidence = infected or infested maize ears total number of ears analyzed × 100 (1)

Statistical Analysis
Data were analyzed for each of the 12 agronomic components by applying three-way ANOVA and using Origin software by OriginLab [30]. Significant influence of bipartite interaction between factors over a given agronomic component was further explored using POLIFACT (non-commercial software) by applying post hoc tests: either Duncan test or the Fishers' Least Significance Difference test. The influence of experimental gradients on each of the agronomic components are presented in Supplementary File 1.

Results
Analysis revealed that experimental factors, their interactions and their gradients had a different influence on the agronomic components of maize, ranging from highly significant to not significant.

Characteristics of Plants
All maize plant characteristics were highly significantly influenced by experimental years. The hybrid exercised a highly significant influence on three plant traits: plant height, ear insertion and leaf length ( Figure 3).

Results
Analysis revealed that experimental factors, their interactions and their gradients had a different influence on the agronomic components of maize, ranging from highly significant to not significant.

Characteristics of Plants
All maize plant characteristics were highly significantly influenced by experimental years. The hybrid exercised a highly significant influence on three plant traits: plant height, ear insertion and leaf length (Figure 3). The fertilization treatments and interaction between fertilization treatment and experimental years did not exercise a significant influence over the maize plant characteristics. Also, tripartite interaction between the factors did not significantly influence the plant morphometric traits (Figure 3).
There were registered differences in value under different fertilization treatments.  The fertilization treatments and interaction between fertilization treatment and experimental years did not exercise a significant influence over the maize plant characteristics. Also, tripartite interaction between the factors did not significantly influence the plant morphometric traits (Figure 3).
There were registered differences in value under different fertilization treatments. The tree-year average was highest for plant height (2.63 ± 0.21 m) and ear insertion height (1.12 ± 0.13 m) under fertilization with Eurofertil Top 51 NPK + CAN. Longest leaves (0.89 m) were obtained by applying Eurofertil Top 51 NPK + CAN as well as Basal + Synertec NK35 N-process. Highest number of leaves above ear (5.57 ± 0.38) was obtained by applying NovaTec Classic. Widest leaves (9.47 cm) were obtained by applying Basal + Synertec NK35 N-process as well as Basal + CAN + Rootip Basic + Energevo (Table S1).
Under the influence of experimental years, the plant morphometric characteristics of the hybrids displayed a variation. In the year 2018, as well as in 2019, the hybrid "Turda 344" had the tallest plants and highest ear insertion, while in 2020, the hybrid "Turda 332" presented the highest values for these two traits compared to values registered by the other two hybrids in the same year (Table 2). Plant height registered by the maize hybrids "Turda 332" and "Turda 344" in 2020, as well as by "Turda 344" in 2019 was significantly higher than values registered by these hybrids in 2018. Similarly, plant height of the hybrid "PR37N01" in 2018 was significantly lower compared to the one registered by the same hybrid in the next two years. Leaf length was significantly higher for the two Turda hybrids in 2020 compared to values registered in the other experimental years by either of the three hybrids tested. As for the leaf number, all three hybrids presented a significantly higher number in 2019 compared to 2018 (Table 2).

Quantitative and Qualitative Productivity Components of Maize
Two quantitative productivity traits of maize (TKW and yield) were both influenced highly significantly by the three experimental factors: years, hybrid and fertilization regime. Also, these two productivity components were also influenced significantly by bi-partite interaction years × hybrid and respectively years × fertilization. However, only TKW was also significantly influenced by tripartite interaction between these factors. Compared to these components, the hectoliter mass was influenced significantly only by the hybrid and not by the other factors or their interactions (Figure 4). This suggests that, while fertilization has the potential to maximize the yield, the hectoliter mass is less susceptible to variation under environmental conditions or inputs.
Quality of the grains was significantly influenced by all three factors, as well as by their bipartite and tripartite interaction. These indicated that the agronomic components could be enhanced through inputs (Figure 4), but fertilization can be optimized and given genotypes.
From the three maize hybrids tested, "Turda 332" led to the obtaining of the highest yield (10,607.98 ± 2657.89 kg/ha), followed by "PR37N01" (10,477.56 ± 2059.38 kg/ha) (Table S2).  Quality of the grains was significantly influenced by all three factors, as well as by their bipartite and tripartite interaction. These indicated that the agronomic components could be enhanced through inputs (Figure 4), but fertilization can be optimized and given genotypes. From the three maize hybrids tested, "Turda 332" led to the obtaining of the highest yield (10,607.98 ± 2657.89 kg/ha), followed by "PR37N01" (10,477.56 ± 2059.38 kg/ha) (Table S2).
A more detailed analysis on the influence of fertilization gradient across experimental years revealed that various fertilizers could enhance productivity components under certain climatic conditions. In the year 2018, all five fertilization regimes ensured significantly higher yields compared to the control, represented by basal fertilization at sowing followed by phasic application of calcium ammonium nitrate. In the year 2020, again, four fertilization regimes ensured significantly higher yields compared to the control ( Table 3).
On average, the highest yield was obtained in the year 2020 (12,432.48 ± 979.08 kg/ha) and the highest TKW was obtained in 2018 (320.80 ± 32.19 g). In the year 2019, was registered both the lowest yield (7482.28 ± 579.03 ka/ha) as well as lowest TKW (254.11 ± 33.38 g). The highest yield was obtained by the hybrid "Turda 332", but the highest TKW and HLM by hybrid "PR37N01" (Table S2). The highest starch content of grains was registered in 2020 (74.31%) and associated with the lowest protein content (6.85%). The highest protein content was registered in 2019 (8.39%) when the yield was lowest. The hybrid A more detailed analysis on the influence of fertilization gradient across experimental years revealed that various fertilizers could enhance productivity components under certain climatic conditions. In the year 2018, all five fertilization regimes ensured significantly higher yields compared to the control, represented by basal fertilization at sowing followed by phasic application of calcium ammonium nitrate. In the year 2020, again, four fertilization regimes ensured significantly higher yields compared to the control (Table 3). On average, the highest yield was obtained in the year 2020 (12,432.48 ± 979.08 kg/ha) and the highest TKW was obtained in 2018 (320.80 ± 32.19 g). In the year 2019, was registered both the lowest yield (7482.28 ± 579.03 ka/ha) as well as lowest TKW (254.11 ± 33.38 g). The highest yield was obtained by the hybrid "Turda 332", but the highest TKW and HLM by hybrid "PR37N01" (Table S2). The highest starch content of grains was registered in 2020 (74.31%) and associated with the lowest protein content (6.85%). The highest protein content was registered in 2019 (8.39%) when the yield was lowest. The hybrid "PR37N01" was characterized by higher average starch content, while the hybrid "Turda 344" by higher protein content compared to the other hybrids (Table S3).
Compared to the control, application of Eurofertil Top 51 NPK + CAN in 2018 ensured significantly higher starch levels compared to the control. In 2020, application of Basal complex fertilizer + Synertec NK35 N-process ensured significant increase in starch content, but at the same time, significant decrease in protein content compared to the control (Table 3).

Disease and Pest Incidence
All three studied factors exercised a significant influence over Fusarium sp., and Ostrinia nubilalis incidence on maize ears, while tripartite interaction between these did not exercise a significant influence. The interaction between experimental years and hybrid was highly significant for disease and pest incidence ( Figure 5). Compared to the three-year average, the Fusarium sp. incidence was significantly lower for the two Turda hybrids in 2020, and significantly higher for the same hybrids in 2018 (Table 4). Similarly, Ostrinia nubilalis incidence followed the same trend by displaying highly significant infestation for the two Turda hybrids in 2018 and significantly lower in 2020 compared to the control (Table 4).
Overall, both Fusarium sp. and Ostrinia nubilalis incidence on maize ears was higher in 2018 and lower in 2020. The lowest incidence was identified for "PR37N01" among the three genotypes studied (Table S4), which might indicate a lower susceptibility of this genotype. Compared to the three-year average, the Fusarium sp. incidence was significantly lower for the two Turda hybrids in 2020, and significantly higher for the same hybrids in 2018 (Table 4). Similarly, Ostrinia nubilalis incidence followed the same trend by displaying highly significant infestation for the two Turda hybrids in 2018 and significantly lower in 2020 compared to the control (Table 4).
Overall, both Fusarium sp. and Ostrinia nubilalis incidence on maize ears was higher in 2018 and lower in 2020. The lowest incidence was identified for "PR37N01" among the three genotypes studied (Table S4), which might indicate a lower susceptibility of this genotype.
Analysis revealed that Fusarium infection of maize ears was lowest under the Basal + CAN fertilization regime, while Ostrinia nubilalis infestation frequency was highly significantly lower than the control under fertilizer regime Basal + CAN + Rootip Basic + Energevo (Table 5).

Discussion
According to the results presented, fertilizer regimes options available today confer the possibility to farmers of enhancing traits of interest according to desired outcomes. Thus, increased biomass, or yields, low disease and pest incidence, high starch or protein content of grains could be obtained by selecting certain combination of fertilizers. The application of basal fertilization in conjunction of CAN, a classic and well-established fertilization regime was outperformed by the fertilization regimes based on the novel fertilizers, combination of several fertilizers or complementing the Basic + CAN, regarding yield and grain quality. This could be attributed to a more complete range of nutrients supplied and higher levels of macronutrients that can be associated with higher overall agronomic performance of maize, according to the current study performed in conditions from the Transylvanian Plain, that experiences continental climate with variations from one year to another.
It remains an important aspect to consider that experimental years had a highly significant influence over 11 out of 12 agronomic components of maize, particularly yield, grain quality and crop health. Low yield registered in 2019 could be attributed to the drought from the summer months of June and July, when very low precipitation values were reported concomitantly with high temperatures. According to the Ministry of Agriculture and Rural Development from Romania, the national average year for maize in 2019 was lower, only 6.502 kg/ha, compared to the average national yield from the previous year, of 7.644 kg/ha [3]. The higher yield of 2020 can also be associated with climatic conditions, since the months of June and July of the respective year were more abundant in precipitations, while also cooler compared to the two previous years enabling the expression of productivity potential of the hybrids tested. In this regard, the results seem to suggest that the fertilizer regime chosen can significantly enhance the productivity components of maize under favorable years for this crop.
Around the world, optimization of fertilizer use for increased production parameters of maize remains an important subject of research that enjoys attention [5,9,[31][32][33][34][35]. Besides application of mineral fertilizers, productivity components of maize were shown to be positively influenced by biostimulants [36] and foliar micronutrient application [12,37,38]. One way to optimize the nutrient use by the crop might be to ensure supply in quantity and at the moment when plants need it, in short, to synchronize the supply with plant demands [11,13]. Differences between genotypes in the uptake and nutrient efficiency use can be attributed to pattern of accumulation and partitioning across growing stages and plant parts [10].
Over a half century of maize crop research, yield per hectare increased continuously, with about half of the gain attributed to genetic improvement [10]. Thus, the role of genotypes and their adaptability to local conditions can prove essential performance criteria for maize genotypes, especially for the Transylvanian Plain, which experiences a continental climate [39]. The results of the study presented here have highlighted that genotype had a significant influence on all agronomic components, while genotype × experimental year interaction also influenced significantly several key agronomic components.
Historically, the use of NPK received the most attention as a means to enhance yields, while the use of micronutrients to increase productivity parameters of crops has been less considered, yet these are just as important for achieving balanced nutrient stewardship [13]. From macronutrients, nitrogen use efficiency has been a subject of attention for a long time in agronomic science, due to its foremost influence on maize yield. There is evidence that the benefits of delaying a part of nitrogen (N) application to late vegetative stages do not bring additional benefit for the crop [6]. Ca and Mg are macronutrients with a pivotal role in plant growth [11], while eight micronutrients are considered essential for healthy growth of plants in small doses and these are: B, Cl, Cu, Fe, Mn, Mo, Ni and Zn. Research shows that the form in which micronutrients are found when applied to the crop have influence on availability, and thus on their impact on the crop growth [13]. Farmers would benefit from increased awareness on the importance of micronutrient application for their crops.
Basic fertilizers containing either one-two components have certain trace elements and thus, mineral fertilizers can increase the level of trace elements in plants and soil, either essential or non-essential. Experimental results have shown that incremental potassium fertilization doses led to a linear increase of essential nutrients such as Ni, Zn and Mn, while higher N fertilization increased Mn and Fe. Potassium and Nitrogen were found to have antagonistic effect on Fe:Zn and Fe:Mn ratios in maize plants [40].
Recent research has demonstrated the lasting impact of P [41] and NP starter fertilization [42] on maize crop, results indicating as optimal localized fertilizer to be applied at 5 cm to the side and below the kernel at sowing [42]. Localized fertilization to ensure uniform plants and boost the crop start is perhaps one of the most interesting approaches for sustainable use of mineral fertilizers.
A comparative study conducted in the United States on maize hybrids created between 1960 and 2000 by DuPont Pioneer identified that newer maize hybrids had higher nutrient uptake capacity compared to older ones. The study also managed to put in evidence the differences in nutrient concentration across plant growth stages and plant-fraction levels. Results showed that macronutrients NP followed overall a similar trend of concentration in tissue that was contrasting with K pattern. Although overall accumulation of N in whole maize plants increased until dent stage, in plant fractions these had a specific pattern. Thus, at silking the N accumulation reached the peak in stalk, leaves and tassel. In vegetative tissue, accumulation of P reached the peak at the silking-blistering stage and decreased once with plant maturity. Similarly, K accumulation in tissues increased until blister stage followed by a decline. Notably, the K content was higher than N or P in cobs. For grains, nutrients accumulation and concentration increased with maturity [10]. Under increased NPK supply, maize shows delayed leaf senescence that ensure postsilking nutrient accumulation and remobilization that translates into higher yields. Modern genotypes that retain green foliage longer exhibit higher post-silking N accumulation and higher dry grain matter levels [8].
Research conducted in Brazil on four transgenic maize hybrids has put in evidence that macronutrient uptake displayed some differences when subject to high versus medium fertilization investment levels. Results showed that, although under high-fertilizer investment the total dry mater yield of plant shoot was higher, the harvest index was lower compared to medium investment. At the end of growing season, an average of 53% N found in the plant was allocated to grains [32]. A study conducted in Mexico on maize hybrid AS-757 cultivated for grain and silage showed that N-190 dose was optimal for delaying inflorescence and obtaining maximum grain yields as well as grain to forage and cobs to stover ratios [43]. An experiment conducted in West Bengal, India, on three maize hybrids subject to nine fertilizer regimes representing 50-150% of the recommended dose and omission of one of the macronutrients (either N, P or K), demonstrated that NPK levels exercised a significant influence on grain yield while incremental NPK levels caused increase in saturated fatty acids concomitantly with decrease of unsaturated fatty acids in maize oil [44]. In addition, omission of N was the most consequential on post-silking and total dry matter and nutrient accumulation followed by P and K omission [8], confirming the long-standing notion that maize is highly dependent on NPK fertilization.

Conclusions
This study assessed the agronomic performance of three maize hybrids in conditions of the Transylvanian Plain, under a fertilizer regime gradient meant to compare traditional fertilization approach with novel options of fertilization available to farmers.
Results showed that fertilization influenced significantly the yield, TKW, grain quality (starch and protein content) and crop health (Fusarium and Ostrinia incidence).
From the six fertilizer regimes tested, the basal fertilization at sowing (NP) coupled with calcium ammonium nitrate at 4-8 leaves was outperformed by all other five fertilizer regimes, that ensured yields of >10 t/ha. This proved the potential of novel fertilization products and approaches to maximize the productivity potential of the genotypes tested.
Experimental year influenced significantly eleven out of twelve agronomic components studied, indicating that fertilization can positively influence the expression of productivity potential of maize genotypes in favorable years for this crop.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/agriculture11090896/s1, Table S1: Influence of experimental gradients on plant characteristics; Table S2: Influence of experimental gradients on agronomic productivity components of maize; Table S3: Influence of experimental gradients on grain quality of maize; Table S4: Influence of experimental gradients on disease and pest incidence.