Nano Urea Fertilization for more Wheat Yield under North Western Plains Zone for Irrigated Conditions Assessed by AMMI and BLUP Analysis

Highly significant variations due to treatments, locations and TxL interactions were observed for wheat yield by analysis of variance during field evaluation of nano urea formulations at number of locations in the north western plans zone of the country. First component of AMMI analysis shared of about 51.3% while second accounted for 29.9%. Significant variations due to locations, TxL interactions and treatments effects were observed for thousands grains weight. Around 82.9% of the interaction effects accounted by two significant components while total of significant components were up to 96.9%. Maximum yield was observed at Karnal followed by Hisar more over maximum yield for 56.1 was expressed by T3 treatment. ASV1 & ASV measures based on 81.3% of interaction effects selected T5, T8 and T8, T5 treatments. Average of thousands grains weight found T3, T8, T9 treatments as suitable for maximum realization and as per MASV1 and MASV settled T6, T2 treatments for thousands grains weight. Adaptability measures corresponding to BLUP estimates of yield i.e. PRVG, PRVG*Gu, HMPRVG*Gu, HMPRVG measures found T3, T4, T2 treatments for yield. Superiority index measures considering average thousands grains weight and stability in 65 and 35 ratios for weighted average settled for T3, T7, T8 treatments formulation of nano urea in the study. Biplot analysis observed Pantnagar center had expressed strong bondage with superiority index measures while considering mean, GAI and HM of treatments based on fixed and random effects of treatments. WAASB, W3, W2, W5 had maintained direct association with MASV, MASV1 on right hand side and with ASV1, W1, IPC4 on left side for yield. Treatments T13, T5 and T4 would express unstable yield as compared to T11, T6, T7 placed near to origin of biplot analysis


Introduction
The UN has proclaimed the 17 sustainable development goals as an effective method of global mobilization to achieve social priorities around the world, such as zero hunger. 1 Sustainability includes agricultural practices without adverse environmental impacts, ensuring the production and quality of fruits and vegetables.Fertilizers have taken axial role with respect to boosting crops yield and nutritional quality especially after the development of fertilizer responsive crop varieties. 2The crop growth need essential nutrients to be appended in the soil are nitrogen, phosphorous, and potassium apart from other micro elements. 3Nitrogen required by the plants is usually supplied in the form of nitric acid, ammonium nitrate, synthetic ammonia, urea, or sodium nitrate.In our country unfortunately, large quantities of urea are applied to fertilize the soils. 4lthough fertilizers are essential for agriculture to feed the growing population, the excessive use of large amounts of chemical fertilizers leads to environmental pollution. 5Besides, only 20-50% of the applied fertilizers are used efficiently; the other 50-80% is lost through leaching, emissions, or incorporation into the soil by microorganisms in the long term, generating ecological problems such as reduced soil fertility and economic losses. 6owever, the key macronutrient elements applied to the soil reported the loss to the tune of 40-70, 80-90, and 50-90%, respectively, causing a considerable loss of resources.Moreover, nitrogen volatilization results in the release of nitrous oxides and thus being the greenhouse gases, contribute to the global warming. 7Loss of mineral nutrients through leaching and runoff to surface and ground water along with abundant volatilization constitute growing concerns owing to economic losses and environmental pollution.Nanofertilizers have been advocated owing to higher NUE as plants cell walls have small pore sizes (up to 20 nm) which result in higher nutrient uptake. 8Plant roots which act as the gateways for nutrients, have been reported to be significantly porous to nano materials compared to conventional manuring materials.Nano fertilizers are more efficacious in terms of nutrients absorption and utilization owing to considerably lesser losses in the form of leaching and volatilization. 9The research findings of a field investigations had proved nano nitrogen fertilizers improved the productivity of crops. 10It was inferred that nano nitrogen fertilizer hold potential to be used in place of mineral urea and it can also reduce environmental pollution caused by leaching, de-nitrification and volatilization of chemical fertilizers.The present study was planned to evaluate the performance of nano urea formulations on wheat crop and association among the commonly exploited measures in the multi locations trails for number of crops.

Materials and Methods
The treatments were evaluated at eight locations (Delhi, Gurdaspur, Gwalior, Hisar, Jammu, Karnal, Ludhiana and Pantnagar centers during 2021-22 cropping season with objective to maximize wheat productivity by optimizing the nitrogen dose and nano urea under irrigated conditions.The recommended agronomical interventions were followed after thorough ploughing and field layering as three replications were maintained.One third of nitrogen along with full dosage of phosphorus and potash as basal one and the remaining 1/3rd of nitrogen at first irrigation and 1/3rd at second irrigations to the plots.The experimental plots were of size 1.80 m x 8 m = 14.40 meter 2 to accommodate the 9 rows with 20 cm spacing among them.Total quantity of spray solution was 400 litre of water/ha with @ of 4 ml of nano urea/litre.The harvested produce of plot size 1.40 m x 7 m = 9.80 meter 2 ( as 7 inner rows x 7 m long were recorded to overcome the effects of border rows) were analysed statistically and homogeneity of locations mean squares were compared by Bartlett where WAASB i was the weighted average of absolute scores of the ith genotype; IPCA ik was the score of the ith genotype (or environment) in the kth IPCA, and EP k was the amount of the variance explained by the kth IPCA.Superiority index has been devised that allowed weights between yield and WAASB as index SI = where rG i and rW i were the rescaled values for yield and, respectively.The superiority index had weighted between yield and stable performance of treatments to be of 65% and 35% respectively.

Analysis of Variance Yield
Highly significant variations due to treatments, locations and TxL interactions were observed (Vaezi et al., 2019) 13 for wheat yield by analysis of variance with 48.5%, 22.2% and 17.6% respectively (Table 2).First component of AMMI analysis shared share of about 51.3% while second accounted for 29.9% and third and fourth were of for 9.3%, 4.3% respectively of interactions sum of squares.Nearly 81.3% of the contributions were of first two significant components while total of 94.9% by significant components.The sums of squares for T×L signal were 79.7% and 20.3% by noise factor towards total of interaction effects.The less share by the residual terms in the analysis suggested further calculations based on interaction principal components for yield.

Thousands Grains Weight
Analysis of variance in table found highly significant variations due to locations, TxL interactions and treatments effects observed for thousands grains weight with respective share of 73.4%, 8.8% and 1.9% respectively (Table 2).First interaction component of AMMI analysis had contributed 51.3% while second accounted for 31.6% and third and fourth were of 9.7%, 4.2% respectively of interactions sum of squares.First two significant components had accounted for 82.9% of the contributions while total of 96.9% by significant components.

Treatments Performance Assessed by Ammi Analysis Measures Yield
Maximum yield values were observed at Karnal location followed by Hisar more over maximum yield for 56.1 was expressed by T3 treatment (Fig 1).IPC1 measure selected T5, T8 treatments and T3, T12 would be desirable as per IPC2 values (Table 3).Measure IPC3 settled for T8, T4 while IPC4 favoured T1, T3 treatment formulation of the current study.IPC5 exhibited suitability of T11, T1 and IPC6 pointed for T13, T8.AMMI analysis based measures (ASV1 & ASV) considering the first two interaction principal components (81.3%%) selected T5, T8 and T8, T5 treatments, while MASV1 along with MASV settled for T8, T2 treatments. 14Average values found T3, T4, T2 treatments as suitable for maximum realization other measures GAI and HM also selected the same treatments.Adaptability measures PRVG, PRVG*Gm selected the treatments T3, T4, T2 as suitability of these treatments backed by values of HMPRVG, HMPRVG*Gm measures.
The recent analytic measure while considering the weighted average of mean value and stability in the ratios of 65 to 35 i.e.Superiority Index pointed for T3, T8, T2 treatments as observed by other SIMe, SIGe, SIHe measures.

Thousands Grains Weight
Gwalior, Hisar and Karnal locations of the zone were able to express the large values for thousands grains weight whereas the maximum thousands grains weight obtained by T3 treatment (Fig 4).Treatments T13, T10 were selected by IPC1 measure and T11, T3 would be desirable as per IPC2 values (Table 6).Measure IPC3 settled for T8, T3 while IPC4 favoured T4, T10 treatment formulation of the current study.IPC5 exhibited suitability of T9, T6 and IPC6 pointed for T9, T1.ASV1 & ASV measures utilized 82.9% of interaction effects selected T7, T5, T6 treatments while MASV1 along with MASV settled for T6, T2 treatments.Average values for thousands grains weight found T3, T8, T9 treatments as suitable for maximum realization other GAI and HM measures values also selected the same treatments.Adaptability measures PRVG, PRVG*Gm selected the treatments T3, T8, T9 as suitability of these treatments were also backed by HMPRVG, HMPRVG*Gm measures.Superiority Index values had pointed for T3, T8, T7 along with values of SIMe, SIGe, SIHe measures.

Superiority Index Measures based on Blup of Treatments Yield
The least value of W1 measure favoured T5, T8 for stable performance whereas W2 selected T5, T8 for the present study (Table 4).T5, T8 by value of W3 whereas W4 favoured T8, T5 while as per W5 treatments T8, T5 would express stable performance. 15Measure WAASB pointed for T8, T5 treatments also.Average value as per BLUP of treatments Gu pointed for T3, T4, T2 and measure GAIu found the higher values of T3, T4, T2 and Hmu measure settled for T3, T4, T2 nano urea formulations.Adaptability measures PRVG, PRVG*Gu found T3, T4, T2 and HMPRVG*Gu, HMPRVG measures had pointed towards T3, T4, T2.Superiority index values had settled for T3, T8, T2 nano urea treatments formulation of in the study.7).Value of W3 and W4 favoured the T11, T6 while W5 selected T11, T6 and WAASB pointed for T11, T6 treatments.Average value as per BLUP estimates of treatments Gu pointed for T3, T4, T7 and GAIu measure found T3, T4, T7 with higher values and Hmu measure settled for T3, T4, T7 treatments for nano fertilization formulations.Adaptability measures corresponding to BLUP estimates found T3, T4, T8 and T3, T4, T7 by PRVG, PRVG*Gu and HMPRVG*Gu, HMPRVG measures respectively.Treatments T13, T5 and T4 would express unstable yield as compared to T11, T6, T7 placed near to origin of biplot analysis.Unstable performance of T3,T4,T5,T8 treatment combinations would be observed as compared to T11, T6, T9, T10 treatments of the present study as placed near to origin in biplot analysis.The measures that brought the lowest variability of the interaction effects were IPC6, IPC5, IPC1 which was manifested by the their shortest vectors among all the studied measures.

Fig. 1 :
Fig. 1: Treatments performance at evaluated centers for wheat yield

Fig. 2 :
Fig. 2: Biplot analysis of Nano treatments and adaptability measures based on PC1 vs PC2 for yield Thousands Grains Weight T13, T10 whereas T11, T6 by least W1 and W2 values for the present study (Table7).Value of W3 and W4 favoured the T11, T6 while W5 selected T11, T6 and WAASB pointed for T11, T6 treatments.Average value as per BLUP estimates of treatments Gu pointed for T3, T4, T7 and GAIu measure found Pantnagar had expressed strong bondage with superiority index measures while considering mean, GAI and HM of treatments based on fixed and random effects of treatments (Fig 2).WAASB, W3, W2, W5 had maintained direct association with MASV, MASV1 on right hand side and with ASV1, W1, IPC4 on left side.IPC3 with Hisar, Ludhiana and Gurdaspur showed direct relationship while Delhi, Karnal, Gwalior with adaptability measures values as per BLUE and BLUP of treatments effects.W1 showed right angles with Hisar, Ludhiana, Gurdaspur values, IPC3 with MASV, MASV1 values, IPC2 with ASV1, WAASB measures, adaptability measures with W1 values in biplot analysis.Right angle of IPC2 observed with Hisar, Ludhiana, Gurdasur centres.

Fig. 3 :
Fig. 3: Clustering pattern of considered measures for multi locations evaluations of treatments for yield

Fig. 4 :
Fig. 4: Variation in thousands grains weight as per nano urea formulations at evaluated centers Thousands Grains Weight First two principal components of the studied measures and treatments combinations explained 75% of total variation.First component accounted for 50.4% while only 24.6% contributed by second one.Measures Hmu, HMPRVGu, HMPRVG*Gu, GAIu, PRVGu, PRVG*Gu were major contributors for first component whereas WAASB, W4, W5, W3, W2, ASV, ASV1 accounted more for second component.Gurdaspur and Jammu were major contributors for first whereas Delhi, Ludhiana was for the second one.Nono urea treatments T13, T3, T7 had contributed more in first component whereas T9, T4, T6 shared more for second principal component (Table 5).T9, T4, T3,T8 treatment combinations would be of unstable performance as far as thousands grains weight was concerned as compared to T2, T1, T5 treatments of the study as placed near to origin in biplot analysis.The measures that brought the lowest variability of the interaction effects were IPC6, IPC4, IPC5, IPC1 which was manifested by the their shortest vectors among all the studied measures.Straight line angle of IPC3 observed with ASV values, Karnal with Ludhiana, Hisar with Jammu, IPC2 with PRVG, MHRVG values (Fig 5).Right angle of W1 was observed with superiority index measures.ASV, ASV1 had ninety degree angles with adaptability measures based on

Fig. 5 :Fig. 6 :
Fig. 5: Biplot analysis of Nano treatments and adaptability measures based on PC1 vs PC2 for thousands grains weight

Table 1 : Details of Nano fertilization treatments and locations of the north western plains zone Code Treatment details of the experiment Code Locations
's test and further analysis was carried out by AMMI soft and SAS 9.3 version software's.A number of AMMI and BLUP measures (Anuradha et al., 2022) 11 mentioned below for ready reference and details about treatments and locations in table 1.
12e stability measure as weighted Average of Absolute Scores has been defined (Olivoto et al., 2019)12as WAASB =

Association Among the Measures as per Biplot Analysis Yield An
advantage of association analysis to club together the treatments or measures into homogenous groups those would assist to select genotypes or measures from different groups as per the objectives of the study.First two principal components of the studied measures and treatments combinations explained 84.7% of total variation.First component accounted for 68.8% while only 15.9% contributed by second one (Table5).Measures SIMu, SIGu, SIHu, HMPRVG, HMPRVG*Gm, HMPRVGu, HMPRVG*Gu, PRVGu, PRVG*Gu were major contributors for first component whereas W1, ASV, ASV1, W4, W5 WAASB IPC3 accounted more for second component.Centres Jammu and Karnal had contributed more for the first while Hisar, Ludhiana for second one.T13, T3, T8 treatments had contributed more in first component whereas T4, T5, T8 shared more for second principal component for yield.
Karnal, Hisar joined with IPV2 in first cluster while Gurdaspur with Gwalior joined hands with superiority index measures in second cluster of biplot analysis (Fig 6).Next cluster was of adaptability measures PRVG, HMPRVG, PRVG*Gu, HMPRVG*Gu, mean, GAI, HM with Jammu, Delhi and Ludhiana centres.Next cluster of IPC4, IPC6 with Pantnagar observed near to large cluster of adaptability measures PRVG, HMPRVG, PRVG*Gu, HMPRVG*Gu.Values of W3 with W2, W4 formed cluster with WAASB values and last cluster of MASV and MASV1 was showed in fourth quadrant.ConclusionsAMMI analysis observed highly significant variations due to treatments, locations and TxL interactions during evaluation of nano urea formulations at number of locations in the north western plans zone of the country.Maximum yield was observed at Karnal followed by Hisar more over maximum yield was expressed by T3 treatment.ASV1 & ASV measures selected T8, T5 treatments.Average of thousands grains weight found T3, T8, T9 treatments as suitable for maximum realization and as per MASV1 and MASV settled T6, T2 treatments for thousands grains weight.Adaptability measures corresponding to BLUP estimates of yield i.e.PRVG, PRVG*Gu, HMPRVG*Gu, HMPRVG measures found T3, T4, T2 treatments for yield.Biplot analysis observed WAASB, W3, W2, W5 had maintained direct association with MASV, MASV1 on right hand side and with ASV1, W1, IPC4 on left side for yield.Treatments T13, T5 and T4 would express unstable yield as compared to T11, T6, T7 placed near to origin of biplot analysis.Thousands grains weight found that Karnal, Hisar centers with IPV2 formed the first cluster while Gurdaspur with Gwalior joined hands with superiority index measures in second cluster of biplot analysis.Next cluster was of adaptability measures PRVG, HMPRVG, PRVG*Gu, HMPRVG*Gu, mean, GAI, HM with Jammu, Delhi and Ludhiana centres.Next cluster of IPC4, IPC6 with Pantnagar observed near to large cluster of adaptability measures PRVG, HMPRVG, PRVG*Gu, HMPRVG*Gu.