Appropriate Application Methods for Salicylic Acid and Plant Nutrients Combinations to Promote Morpho-Physiological Traits, Production, and Water Use Efficiency of Wheat under Normal and Deficit Irrigation in an Arid Climate

Freshwater shortage and inadequate nutrient management are the two major challenges for sustainable wheat production in arid agro-ecosystems. Relatively little is known about the positive roles of the application methods for the combination of salicylic acid (SA) and plant nutrients in sustaining wheat production under arid climatic conditions. A two-year field study was undertaken to assess the impact of seven treatments for the integrated application of SA, macronutrients, and micronutrients on the morpho-physiological traits, yield, and irrigation water use efficiency (IWUE) of wheat subjected to full (FL) and limited (LM) irrigation regimes. The results showed that the LM regime caused a significant reduction in different plant growth traits, relative water content, chlorophyll pigments, yield components, and yield, while a significant increase was observed in IWUE. The sole application of SA or co-application with micronutrients through soil did not significantly affect the studied traits under the FL regime, while they achieved some improvement over untreated plants under the LM regime. Based on the different multivariate analyses, the soil and foliar applications for the combinations of SA and micronutrients, as well as a foliar application for the combinations of SA, macronutrients, and micronutrients were identified as an efficient option for mitigating the negative impacts of water deficit stress and enhancing the growth and production of wheat under normal conditions. In conclusion, the results obtained herein indicated that the co-application of SA and macro- and micronutrients is an effective option to greatly enhance and improve the growth and production of wheat crops in water-scarce countries of arid regions, such as Saudi Arabia, while an appropriate application method for this combination was required for positive effects.


Introduction
Water scarcity is one of the significant environmental issues that restrict food crop production and affect all factors related to ensuring global food security. Unfortunately, several studies have predicted that water scarcity will increase due to continued human activity and increased emissions leading to global warming and climate change, particularly in arid and semiarid regions [1,2]. It is predicted that water scarcity will affect more than 50% of global croplands by 2050. Therefore, food crop production is still far below the expected quantities to meet the increase in food demand.
Bread wheat (Triticum aestivum L.) is one of the most important cereal crops produced and consumed worldwide. It provides approximately 20% of the protein and daily calories for 4.5 billion people around the world. It is grown on a total area of approximately 219 million ha, with a total production of over 760 million tons in 2020 [3]. Demand for wheat is estimated to increase by 60% by 2050, posing a serious concern for meeting that

Growth Parameters
The analysis of variance showed that the F-values, due to irrigation regimes (IR) and different combination treatments between SA, macronutrients and micronutrients (T), were significant at p ≤ 0.05 or 0.01 for all growth parameters in each growing season and the combined two seasons, except for plant height (PH) in the second season and the combined two seasons and tiller number (TN) in the two seasons (Table 1). All growth parameters in both growing seasons and the combined two seasons, except PH and TN, were significantly affected at p ≤ 0.01 by the interaction IR × T. Seasons (S) had an insignificant effect on all growth parameters, except shoot dry weight (SDW). The interaction IR × S, T × S, and IR × T × S showed insignificant effects on all growth parameters at p ≤ 0.05 or 0.01, except TN, green leaf number (GLN), leaf area index (LAI), and shoot fresh weight (SFW), which were significantly affected by IR × S interaction (Table 1). Averaged over the two seasons, the limited irrigation regime (LM) hampered PH, TN, GLN, LAI, SFW, and SDW by 12.2%, 21.5%, 43.6%, 47.9%, 35.7%, and 32.1%, respectively, when compared with the full irrigation regime (FL) ( Table 2). However, the exogenous applications of different T substantially improved all the parameters mentioned above as compared to the non-treated control treatment (T1). In general, T4 (the soil and foliar application of SA and micronutrients) exhibited the highest values for all growth parameters followed by T7 (foliar application of SA, macronutrients, and micronutrients), while the T1 exhibited the lowest ones (Table 2). Over the two seasons, the T2 (soil application of SA), T3 (soil application of SA and micronutrients), T4, T5 (soil application of SA and foliar application of macronutrients and micronutrients), T6 (foliar application of SA), and T7 increased GLN by 4.3-41.3%, LAI by 10.5-53.3%, SFW by 5.1-26.0%, and SDW by 6.3-27.1%, as compared to T1 (Table 2). Table 2. Effects of irrigation regimes and combination treatments (T) on growth parameters of wheat in first season (S1), second season (S2), and across both seasons (comb).
T S1 S2 Comb. S1 S2 Comb. S1 S2 Comb. The responses of growth parameters to different T under each IR regime over both growing seasons are shown in Figure 1. The different T (T2−T7) improved the LAI, GLN, SFW, and SDW under FL by 1.9-41.1%, 1.6-43.2%, 2.0-26.3, and 2.6-23.8%, respectively, in comparison to T1. The maximum increase in the parameters mentioned above was observed with T7 and T4, which increased by 22.5-40.7% in comparison to their corresponding T1. However, the T2 and T3 did not improve the aforementioned growth parameters under the FL regime, and the values of these parameters for both treatments were statistically on par with the T1 (Figure 1). Contrarily, all T improved the parameters mentioned above under the LM regime, and T4 showed a maximum increase in these parameters by 36.8-74.2%, as compared with the T1 (Figure 1).

Physiological Parameters
The IR regime, T, and S had a significant main effect on all physiological parameters, namely relative water content (RWC), chlorophyll a (Chla), chlorophyll b (Chlb), and total chlorophyll content (Chlt) in each growing season as well as the combined two seasons (Table 3). Additionally, the four physiological parameters were significantly affected by the interaction IR × T, except Chlb in the first season, IR × S, and T × S, except RWC. The interaction effect of IR × T × S had a significant effect on Chla at the 5% probability level and a non-significant effect on other physiological parameters (Table 3).  5  T1  T2  T3  T4   T5  T6  T7 Irrigation treatments Response of different growth parameters to the interaction between irrigation regimes and treatments over the two growing seasons. The full name of the abbreviations of different combination treatments (T1-T7) is listed in the footer of Table 2. The different letters indicate statistical significance at the 0.05 level.

Physiological Parameters
The IR regime, T, and S had a significant main effect on all physiological parameters, namely relative water content (RWC), chlorophyll a (Chla), chlorophyll b (Chlb), and total chlorophyll content (Chlt) in each growing season as well as the combined two seasons (Table 3). Additionally, the four physiological parameters were significantly affected by the interaction IR × T, except Chlb in the first season, IR × S, and T × S, except RWC. The interaction effect of IR × T × S had a significant effect on Chla at the 5% probability level and a non-significant effect on other physiological parameters (Table 3).
Irrespective of T, averaged over both seasons, the LM regime significantly reduced RWC, Cha, Chb, and Cht by 22.4 % , 39.8%, 40.8%, and 41.1%, respectively, as compared to the FL regime (Table 4). Regarding the effect of different T, the T4 showed the highest values for all physiological parameters, followed by T7. Averaged over both seasons, the T4 and T7 enhanced 12.0% and 7.4% more RWC, 32.7% and 30.6% more Chla, 33.3% and 29.0% more Chlb, and 32.7% and 29.9% more Chlt as compared to the T1, respectively (Table 4). Similarly, the other treatments (T2, T3, T5, and T6) were accompanied by in-    (Table 4). Regarding the effect of different T, the T4 showed the highest values for all physiological parameters, followed by T7. Averaged over both seasons, the T4 and T7 enhanced 12.0% and 7.4% more RWC, 32.7% and 30.6% more Chla, 33.3% and 29.0% more Chlb, and 32.7% and 29.9% more Chlt as compared to the T1, respectively (Table 4). Similarly, the other treatments (T2, T3, T5, and T6) were accompanied by increases in RWC by 5.9-8.7%, Chla by 15.9-24.4%, Chlb by 17.0-21.9%, and Chlt by 16.0-23.2% compared to the T1 (Table 4). Table 4. Effects of irrigation regimes and combination treatments (T) on physiological parameters of wheat in the first season (S1), second season (S2), and across both seasons (comb). The full name of the abbreviations of different physiological parameters is listed in the footer of Table 3. The full name of the abbreviations of the irrigation regimes and different combination treatments is listed in the footer of The responses of physiological parameters to different T under each irrigation regime are presented in Figure 2. Under the FL regime, plants treated with foliar application (T6 and T7 treatments) or combined foliar and soil applications (T4 and T5 treatments) showed significant superiority in RWC over those plants treated with soil application only (T2 and T3 treatments) or untreated plants (T1). However, under the LM regime, the T2-T7 treatments significantly outperformed the T1 in RWC by 7.9-25.4% ( Figure 2). Regarding chlorophyll pigments under the FL regime, the highest values for these parameters were achieved with T7 and T4, followed by T3, T5, and T6, while the T2 were statistically on par with the T1. However, under the LM regime, the T2-T7 exhibited higher chlorophyll pigments than those of the T1, with the treatments of the foliar application or a combination of soil and foliar applications were better than those treatments of sole soil application in enhancing the chlorophyll pigments ( Figure 2). a e d a bc ab cd ab these parameters were achieved with T7 and T4, followed by T3, T5, and T6, while the T2 were statistically on par with the T1. However, under the LM regime, the T2-T7 exhibited higher chlorophyll pigments than those of the T1, with the treatments of the foliar application or a combination of soil and foliar applications were better than those treatments of sole soil application in enhancing the chlorophyll pigments ( Figure 2).  Table 2. The different letters indicate statistical significance at the 0.05 level, according to Duncan's test.

Yield and Yield Attributes Parameters
The F-values of the ANOVA analysis showed that all yield components and irrigation water use efficiency (IWUE) were significantly affected by the IR regime, T, and their interaction, with the exception of spike length (SL), in which the T and IR × T interaction had non-significant effects on this parameter ( Table 5). The S showed a significant effect on only the grain weight per spike (GWS), number of grains per spike (NGPS), and harvest index (HI). The IR × S interaction showed non-significant effects on any yield parameters and IWUE, and the T × S and IR × T × S interaction showed significant effects on only HI (Table 5).  Table 2. The different letters indicate statistical significance at the 0.05 level, according to Duncan's test.

Yield and Yield Attributes Parameters
The F-values of the ANOVA analysis showed that all yield components and irrigation water use efficiency (IWUE) were significantly affected by the IR regime, T, and their interaction, with the exception of spike length (SL), in which the T and IR × T interaction had non-significant effects on this parameter ( Table 5). The S showed a significant effect on only the grain weight per spike (GWS), number of grains per spike (NGPS), and harvest index (HI). The IR × S interaction showed non-significant effects on any yield parameters and IWUE, and the T × S and IR × T × S interaction showed significant effects on only HI (Table 5). Table 5. Analysis of variance (F-values) for the effects of irrigation regimes (IR), combination treatments (T), and their interaction on yield components, yield, and irrigation water use efficiency for each season and combined across the two seasons. Averaged over both seasons, the LM regime significantly reduced SL, GWS, GNPS, thousand grain weight (TGW), grain yield per ha (GY), biological yield per ha (BY), and HI by 16.2%, 38.5%, 26.5%, 16.5%, 41.8%, 29.1%, and 18.1%, respectively, while it increased irrigation water use efficiency (IWUE) by 14.1%, when compared with the FL regime ( Table 6). Regardless of the IR, the T2−T7 treatments considerably enhanced GWS, GNPS, TGW, GY, BY, HI, and IWUE over the two seasons by 6.4-30.5%, 1.4-14.0%, 6.3-17.1%, 12.7-39.0%, 9.0-26.8%, 5.0-11.8%, and 16.5-45.3%, respectively, as compared with the T1 ( Table 6). The maximum increases in yield components, yield, and IWUE were recorded with T4, followed by T7. Averaged over two seasons, the T4 and T7 increased the different parameters of yield and IWUE by 17.1-45.3% and 7.2-34.5%, respectively, in comparison to the T1, with the Duncan's test in decreasing order of these parameters were T4 > T7 > T6 > T5 > T3 > T2 > T1 ( Table 6).

Sources
The responses of yield components, yield, and IWUE to the exogenous applications of combination treatments under each irrigation regime are shown in Figure 3. Under the FL regime, the maximum values for yield parameters and IWUE were obtained with the treatments of foliar application (T6 and T7) and a combination of foliar and soil application (T4 and T5), with the T4 and T7 always exhibiting the highest values for these parameters. The treatments of soil application only (T2 and T3) are less efficient than the other combination treatments in improving the yield parameters and IWUE, and the values of these parameters for both treatments were statistically on par with the T1 under the FL regime ( Figure 3). However, under the LM regime, all combination treatments improved yield parameters and IWUE compared to the T1. Compared with the T1, the increases in GWS, NGPS, TGW, GY, BY, HI, and IWUE due to treated plants with different combination treatments (T2-T7) was 4.2-2.7%, 13.9-30.5%, 28.1-64.6%, 16.1-35.5%, 10.6-22.1%, and 28.1-64.6%, respectively ( Figure 3). Overall, the treatment of a combination of soil and foliar applications (T4) recorded the highest values for all yield parameters and IWUE and statistically outperformed the other treatments, except for TGW and BY (Figure 3). Table 6. The effects of irrigation regimes and combination treatments (T) on yield components, yield, and irrigation water use efficiency of wheat in first season (S1), second season (S2), and across both seasons (Comb). T S1 S2 Comb S1 S2 Comb S1 S2 Comb S1 S2 Comb

Pearson's Correlation Coefficient among Different Parameters under the Full and Limited Irrigation Regimes
In Table 7, which shows the Pearson's correlations values between growth, physiological, yield components, yield, and IWUE parameters under each IR regime, all the studied parameters showed a significant positive correlation (r = 0.77-0.99) with each other, except PH and SL under both irrigation regimes and TN under the FL regime, which displayed a weak and occasionally even a negative correlation with all parameters. Furthermore, there were positive and significant correlations of GY and IWUE with all studied parameters (r = 0.78-0.99) except PH, TN, and SL under the FL regime. Similarly, GY and IWUE were strongly positively correlated with all studied parameters (r = 0.77-0.99) except SL under the LM regime. The RWC and Chlt showed significant positive correlation with almost all growth and yield parameters under both irrigation regimes (Table 7).  The full name of the abbreviations of different growth, physiological, and yield and irrigation water use efficiency parameters is listed in the footers of Tables 1, 3 Table 2. The different letters indicate statistical significance at the 0.05 level.

Principal Component Analysis and Heatmap Analysis for the Responses of Plant parameters to Different Combination Treatments
The associations between the 16 studied parameters and 7 combination treatments under each irrigation regime were evaluated via a principal component analysis (PCA) scatter plot (Figure 4) and heatmap clustering ( Figure 5). The first two principal components (PC1 and PC2) explained 85.94% and 92.10 of the total variations derived from the 16 parameters under the FL and LM regimes, respectively (Figure 4). The PC1 explained 76.59% and 85.90% of the total variance, while the PC2 explained 9.35% and 6.19% of the total variability between all parameters under the FL and LM regimes, respectively. Both PCs grouped all studied parameters together in one group in a positive direction of PC1 under both irrigation regimes, with the angle between the vector of these parameters being acute, which indicates a positive association between these parameters ( Figure 4). Additionally, all parameters were strongly positively associated with T4 (first quarter), T7 and T6 (second quarter) under the FL regime ( Figure 4) and T4 (first quarter), T7, T6, and T5 (second quarter) under the LM regime ( Figure 4). The combination treatments in the first and second quarters of PCA were leaned toward the highest PC1 and lowest or highest PC2. All studied parameters did not show any associations with T1, T2, and T3 under both irrigation regimes, in addition to T5 under the FL regime ( Figure 4).   Figure 5 shows that the heatmap dendrogram distributed the different combination treatments into three and four groups under the FL and LM regimes, respectively. The T4 and T7 under the FL regime and T4 alone under the LM regime were clustered into one group and displayed the highest values for all studied parameters under FL, except the PH and TN, and the LM regimes ( Figure 5). The T1 was clustered alone in one group under the LM regime, while it clustered with T2 and T3 (treatments of soil application only) in one group under the FL regime. The T5 and T6 were clustered together in one  Table 2.  Tables 1, 3, and 5, respectively. The full name of the abbreviations of irrigation regimes and different combination treatments is listed in the footer of Table 2.

Discussion
In this study, all studied parameters were significantly affected by irrigation regimes (Tables 1, 3, and 5) and the LM regime considerably reduced the growth, physiological, and yield parameters compared with the FL regime (Tables 2, 4, and 6). These reductions in the different studied parameters under the LM regime may be due to the exposure of plants to water deficit stress, resulting in substantial inhibition in several morphological, physiological, and biochemical attributes, including limiting RWC, photosynthetic rate, stomatal conductiance, activity of photosynthetic enzymes, chlorophyll biosynthesis, cell division, cell expansion, and biomass accumulation. Additionally, water deficit stress indirectly affects plant growth and production by hindering the access of nutrients to the plant and reducing their mobility and absorption even with their availability in the soil [39][40][41]. A shortage of nutrients, or at least some of them, leads to physiological and biochemical disturbances in the plant as a result of the role of these elements in important vital processes, such as cell osmotic relations and turgor-related processes, photosynthesis reactions, enzymatic activity, building nucleic acids, and plant reproduction [12,15,42]. Furthermore, water deficit stress induces oxidative damage that may affect plant water relations, leaf pigments, and cell membrane integrity as well as leading to the generation of ROS [43,44]. Therefore, previous studies have reported that the synthesis  Figure 5 shows that the heatmap dendrogram distributed the different combination treatments into three and four groups under the FL and LM regimes, respectively. The T4 and T7 under the FL regime and T4 alone under the LM regime were clustered into one group and displayed the highest values for all studied parameters under FL, except the PH and TN, and the LM regimes ( Figure 5). The T1 was clustered alone in one group under the LM regime, while it clustered with T2 and T3 (treatments of soil application only) in one group under the FL regime. The T5 and T6 were clustered together in one group under the FL regime, while they clustered with the T7 treatment in one group under the LM regime ( Figure 5).

Discussion
In this study, all studied parameters were significantly affected by irrigation regimes (Tables 1, 3 and 5) and the LM regime considerably reduced the growth, physiological, and yield parameters compared with the FL regime (Tables 2, 4 and 6). These reductions in the different studied parameters under the LM regime may be due to the exposure of plants to water deficit stress, resulting in substantial inhibition in several morphological, physiological, and biochemical attributes, including limiting RWC, photosynthetic rate, stomatal conductiance, activity of photosynthetic enzymes, chlorophyll biosynthesis, cell division, cell expansion, and biomass accumulation. Additionally, water deficit stress indirectly affects plant growth and production by hindering the access of nutrients to the plant and reducing their mobility and absorption even with their availability in the soil [39][40][41]. A shortage of nutrients, or at least some of them, leads to physiological and biochemical disturbances in the plant as a result of the role of these elements in important vital processes, such as cell osmotic relations and turgor-related processes, photosynthesis reactions, enzymatic activity, building nucleic acids, and plant reproduction [12,15,42]. Furthermore, water deficit stress induces oxidative damage that may affect plant water relations, leaf pigments, and cell membrane integrity as well as leading to the generation of ROS [43,44]. Therefore, previous studies have reported that the synthesis and accumulation of compatible solutes are considered one of the most common methods of maintaining plant water relations and removing excess levels of ROS under water deficit stress [31][32][33][34][35]. SA is one of the most promising of these compatible solutes in enhancing the growth and production of crops under water deficit stress by enhancing osmotic regulation, reserving water in plant cells, promising water potential gradients, enhancing the activities of antioxidant enzymes, concealing the ROS, and regulating the uptake and assimilation of many nutrients [34,35,45]. Therefore, we hypothesized in this study that the exogenous co-application of SA, macronutrients, and micronutrients through soil and/or foliar spraying methods could be considered as a cost-effective and easy strategy for not only mitigating the water deficit stress on wheat growth and production but also enhancing the performance of wheat crop under normal irrigation conditions in arid agro-ecosystems.
The current study showed that the combined application of SA, macronutrients, and micronutrients through soil and foliar spray methods improved the performance of wheat under the FL regime as well as having the ability to ameliorate the negative impacts of water deficit stress on morpho-physiological traits and the production of wheat under LM regimes (Tables 2, 4 and 6). Several previous studies indicated the efficacy application of SA and/or nutrients, such as N, P, K, Zn, and Mn, to enhance the growth and yield of many crops such as wheat, rice, soybeans, mung beans, and maize under water deficit stress and reduce the adverse effects of water shortage on the production of these crops [33,[46][47][48][49]. In this regard, El Sherbiny et al. [33] concluded that the foliar application of SA at 700 µM could be used to improve the growth, production, and WUE of rice as well as to mitigate the negative impacts of limited water irrigation. Razmi et al. [46] also found that spraying soybean plants exposed to water deficit stress with SA improved the parameters of LAI, RWC, and chlorophyll a and b. Similarly, Ahmad et al. [47] indicated that the exogenous application of SA in wheat enhanced the leaf water potential, osmotic potential, chlorophyll contents, photosynthetic rate, RWC, and antioxidant enzymatic activities under both normal and stress conditions. On the other hand, Waraich et al. [50] reported that application of some macronutrients and micronutrients, such as N, P, K, Zn, Si, and Mg, can reduce the negative impacts of drought stress on plants. There are several mechanisms by which these elements reduce the negative impacts of water deficit stress as well as enhance the performance of plant crops under normal conditions. For example, foliar-applied fertilizer of N enhanced the process of photosynthesis and regulated several metabolic processes in plants as well as enhancing the leaf Chl contents and the amount of photosynthetically active radiation intercepted by plants, which eventually leads to further improving plant growth and yield parameters [24]. K is involved in a wide range of plant processes, including photosynthesis, stomatal regulation, enzyme activity, osmoregulation and membrane stability, carbohydrate synthesis, and the transportation of assimilates from source to sink, which all are important to enhance the growth and production of crops under both normal and stress conditions [14,15,22]. P is essential for numerous cellular functions, including energy transfer, enzyme activation, membrane structure maintenance, biomolecule synthesis, and the creation of high-energy molecules (adenosine triphosphate). Moreover, it enhances root development, resulting in increased water uptake and enhancing the growth and productivity of crop plants in water-scarce environments [13,[51][52][53].
Additionally, foliar-applied micronutrients, such as Zn and Mn, enhanced the growth and production of several field crops under both normal and water stress conditions because both ions are involved in a wide range of developmental and physiological processes, including enhancing auxin formation, cell expansion, stem elongation, source-sink relationship, membrane stability, hormone and chlorophyll synthesis, enzyme activation, and antioxidant defense [18,20,54,55]. Furthermore, the exogenous application of SA, which also is involved in several physiological, biochemical, and development processes in plants, also play an important role in regulating plant growth and development under both normal and stress conditions. It improves the growth and production of crops by stimulating the cell division, photosynthetic activity, and water potential gradients, regulating some physiological responses related to carbon uptake and/or fixation, increasing the creation of photosynthetic pigments, raising stomatal conductance, delaying the senescence of plant organs, and regulating the source-to-sink relationship [33][34][35].
In this study, the response of different morpho-physiological and yield parameters of wheat to the co-application of SA, macronutrients, and micronutrients varied according to the irrigation regimes (Figures 1-3). Under the FL regime, the co-application of SA, macronutrients, and micronutrients through foliar application (T7) or foliar and soil applications (T4) were more efficient for enhancing the different studied parameters than those of sole application of SA (T1 and T6), sole application of SA through soil with the sole application of nutrients through foliar (T5), or the co-application of SA and micronutrients through soil application only (T3), while the T2 and T3 did not show any significant variation from T1 (Figures 1-3). However, under the LM regime, the T4 achieved the highest values for all studied parameters, followed by T5, T6, and T7. Additionally, the T2 and T3 achieved higher values for the most studied parameters than the T1 under the LM regime (Figures 1-3). These results indicate that the co-application of SA with micronutrients through both soil and foliar spray methods (T4) is a useful and simple approach to enhance the growth and production of wheat crops under both normal and water deficit stress conditions. This means that applying a small amount of SA and micronutrients through foliar sprays directly on the foliage, in addition to applying a low dose of these materials through soil application, are efficient approaches to improving the growth, physiological, production, and IWUE of the wheat crop when compared with using only one mode of application for these materials. As the roots are the first organ to experience water deficit stress, in addition to the leaves being the key organ in which the majority of biochemical and physiological processes occur, this may explain why the combination of soil and foliar applications methods for the combinations of SA and micronutrients seems to be an effective approach for enhancing the growth and production of wheat under normal and water deficit stress conditions. In addition, the foliar application of SA and micronutrients helps their entry directly into the leaves, which raises their concentration in them sufficiently. This helps protect various biochemical and physiological processes from the negative effects of water deficit stress, leading to rapid relief of physiological stress, and also promoting photosynthesis and regulating other metabolic processes in the plant. While the application of SA and micronutrients through soil significantly raises their concentration in the roots, which enhances root growth under both normal and stress conditions as well as improving the osmotic adjustment in roots under stress conditions. This enables the roots to extract water and nutrients from deeper soil layers and the presence of SA as compatible solutes in the roots and micronutrients, such as Zn and Mn, enhance water uptake and promising water potential gradients under water deficit stress conditions. These findings are in line with the previous study conducted by Karim and Rahman [21], who reported that the soil application of micronutrients (Zn, B, and Mn) at an early stage of cereal crops in combination with foliar application at a late stage is a promising approach to alleviating the negative impacts of drought stress, which often occur during the late stages of cereal production in Asia's least developed countries. They also found that the foliar application of micronutrients (Zn, B, and Mn) along with the soil application of macronutrients (NPK) significantly enhanced cereal production under drought stress conditions. The results of this study also found that the co-application of SA, macronutrients,-and micronutrients through foliar sprays only (T7) is a more efficient strategy for enhancing the performance of wheat crops under the FL regime than under the LM regime, as shown with the heatmap clustering for different combination treatments under the FL and LM conditions ( Figure 5). This finding indicates that the foliar application of macronutrients and micronutrients can be considered an essential strategy for achieving a higher productivity of wheat under normal conditions. This result can be attributed to the fact that the foliar application of nutrients can increase the efficiency of nutrient absorption and delay leaf senescence; then, photosynthesis can be enhanced and/or prolonged. Similarly, improved productivity, quality, and profitability in soybean crops were reported by Dass et al. [24] and Gheshlaghi et al. [56] using the foliar application of essential nutrients under normal conditions. However, the soil application of micronutrients has been noticed to be very effective in alleviating the adverse effects of water deficit stress on wheat performance under the LM regime. Similarly, Ma et al. [57] found that the soil application of Zn led to an increase in wheat grain yield under stress conditions as well as under adequate water supply. This increase was directly proportional to the intensity of the stress and was 28.2%, 22.6%, and 10.5% in severe drought, moderate drought, and adequate water supply, respectively. The same was reported by Dimkpa et al. [58], who found that the soil application of micronutrients was more effective in reducing the negative effects of water deficit stress on soybeans.
It is observed from this study that yield, yield components, and IWUE were significantly and strongly correlated with each other and with the different morpho-physiological parameters under both irrigation regimes (Table 7). This result indicates that the coapplication of SA, macronutrients, and micronutrients positively affected several morphophysiological and biochemical aspects of wheat plants under both water and non-water deficit conditions. This may be because the integration between SA and plant nutrients may play a defensive role in mitigating the negative impacts of water deficit stress on wheat performance under stress conditions while playing a synergistic role in enhancing wheat performance under normal conditions. Therefore, a statistically significant correlation was observed between most studied traits under the FL and LM regimes (Table 7). Similarly, Dass et al. [24] and El Sherbiny et al. [33] reported that strong and positive correlations were found between most studied parameters of soybean and rice crops under both water and non-water deficit conditions due to the exogenous application of SA, macronutrients, and micronutrients. This is because the SA and/or macronutrients, and micronutrients are involved in a wide range of developmental and physiological processes, not only under stress conditions but also under normal conditions. For instance, plants of many species treated with SA increase their concentrations of essential elements such as N, P, K, Ca, S, Fe, Mn, and Mg [59][60][61]. This mineral nutrient ultimately results in a significant enhancement in the chlorophyll a, b, and carotenoids content, photosynthetic rate, enzymatic and nonenzymatic antioxidant activity, and the production of carbohydrates, thus improving the growth and production of plants treated with SA under both normal and environmental stress conditions as compared to untreated plants under stress conditions. Zn is a constituent of various enzymes that facilitate the diffusion of CO 2 from stomatal apertures to the carboxylation sites by Rubisco and the export of carbohydrates from source to sink [24,62]. It also plays a special role in cell division, nitrogen metabolism, photosynthesis, and the synthesis of protein auxin, RNA, and DNA [63]. The macronutrients present in the combination treatments, particularly NPK, when applied through the foliar spray could have boosted the photosynthesis process and regulated other metabolic processes in the plant [24,64]. K plays a positive role in the accumulation of carbohydrates in source (leaves) and its translocation to sink (grains), thus improving the yield. It is also plays a key role in the maintenance of turgor pressure and osmotic adjustment under stress conditions [14,15,22]. All the abovementioned advantages of treated plants with SA, macronutrients, and micronutrients under both normal and stress conditions might explain why strong and positive correlations between most studied parameters were observed under both the FL and LM regimes.

Experimental Site Describtion and Cultivation Conditions
Two field experiments were conducted from December to April in two successive seasons of 2020/2021 and 2021/2022 at the Agricultural Research Farm of the Plant Production Department, College of Food and Agriculture Sciences, King Saud University, Riyadh, Saudi Arabia (24 • 25 03" N 46 • 39 17" E, 570 m above sea level) to investigate the effect of the exogenous co-application of SA, macronutrients, and micronutrients through soil and/or foliar spray methods on the morpho-physiological traits, yield attributes, and IWUE of wheat under FL and LM irrigation regimes. The soil texture of the research farm is classified as a sandy loam (57.92% sand, 28.65 % silt, and 13.42% clay), and it is climatically categorized as typical arid climatic conditions. The monthly averages of the meteorological data at the research farm during both growing winter seasons are shown in Table 8. Table 8. Monthly averages of the climatic data at the research farm during the growing period of wheat in the first (S1) and second (S2) seasons.

Precipitation
(mm) S1 S2 S1 S2 S1 S2 S1 S2 The soil was ploughed twice, leveled, and divided into plots (4 m × 2 m each). Phosphorus fertilizer was applied during the seedbed preparation at the rate of 90 kg P 2 O 5 ha −1 in the form of calcium superphosphate (17% P 2 O 5 ). Then, seeds of the spring wheat cultivar Summit were planted manually in ten rows, each 4 m long and 20 cm apart from the adjacent row, at a seeding rate of 150 kg ha −1 . The seeds were planted on December 8th and 1st in the first and second seasons, respectively. Other nutrients, potassium (K) and nitrogen (N), were applied at the rate of 60 kg K 2 O and 180 kg N ha −1 in the form of potassium sulfate (50% K 2 O) and urea (46% N), respectively. Potassium and nitrogen fertilizers were applied at two and three equal doses, respectively. The first and second doses of both fertilizers were applied at the seedling and stem elongation stages, while the third dose of N was applied at booting stage. Other agronomic practices, such as the control of weeds and disease, were done in a timely manner.

Experimental Design and Treatments
The experiments were accomplished in a randomized complete block design with a split plot arrangement and replicated three times. The two irrigation regimes (FL and LM regimes) were assigned to the main plot, while the different combination treatments were randomly distributed in subplots. The first irrigation regime represents non-waterstressed plants (100% of the estimated crop evapotranspiration; ETc), while the second one represents water-stressed plants (50% ETc). The different combination treatments included a possible combination between the application methods for the co-application of SA, macronutrients (NPK), and micronutrients (Zn and Mn). These treatments included control (T1), the soil application of SA (T2), soil application of SA and micronutrients (T3), soil and foliar application of SA and micronutrients (T4), soil application of SA and foliar application of macronutrients and micronutrients (T5), foliar application of SA (T6), and foliar application of SA, macronutrients, and micronutrients (T7).
Based on the reference evapotranspiration rate (ETo) and crop coefficient (Kc), the amount of irrigation water required (ETc) for the FL regime was calculated using the following Equation: ETc = ETo × Kc (1) The ETo was determined based on the modified Penman-Monteith equation [65], while the Kc values of spring wheat provided in the FAO-56 were used after adjustment based on the climatic conditions of the research farm. Based on this calculation, the cumulative irrigation volume for the FL regime was 6470 and 6500 m 3 ha −1 for the first and second seasons, respectively. The LM regime received 50% of the irrigation volume of the FL regime. The irrigation water was applied using the modified surface irrigation method, as outlined by El-Hendawy et al. [66]. The irrigation treatments were applied 20 days after sowing.
The soil application of SA and micronutrients were applied 30 days after sowing. The SA, Zn, and Mn were applied at a rate of 3, 20, and 15 kg ha −1 , respectively. SA, Zn and Mn were applied in the form of HOC 6 H 4 COOH, ZnSO 4 ·7H 2 O, and MnSO 4 ·3H 2 O, respectively. The foliar application treatments were applied twice at 40 and 60 days after sowing at the concentrations of 0.5% for Mn, 1% for N, P, K, and Zn, and 2.0 mM for SA. Salicylic acid was dissolved in absolute ethanol, while the other compounds were dissolved in distilled water to make a stock. Then, the appropriate quantities of this stock were gently added to the distilled water to make spray solutions with the required concentrations. The various combinations of SA solution and nutrient elements containing 0.1% Tween-20 were sprayed directly onto the leaves of the plants to the point of runoff using a back-mounted pressurized sprayer (16 L) with a T-jet nozzle that was calibrated to deliver 15 mL s-1 at a pressure of 207 kPa.

Growth Parameters
Ten plants were randomly harvested from each plot 100 days after sowing to record SFW, PH, TN, and GLN. Then, the surface area of all green leaf blades was measured using an area meter (LI 3100; LI-COR Inc., Lincoln, NE, USA) to record the leaf area per plant, which was used to calculate the leaf area index (LAI) according to Equation (2). All parts of the ten plants were oven-dried at 80 • C to a constant weight and then weighed to obtain the SDW.

LAI =
Leaf area per plant cm 2 Ground area (cm 2 )

Physiological Parameters
A second leaf of five randomly selected plants was taken to record the RWC using the following equation [67]: where FW is the fresh weight of an area of 7-10 cm 2 from each leaf, TW is the turgid weight of leaf samples after being rehydrated in distilled water for 24 h in the darkroom at 25 • C, and DW is the dry weight of leaf samples after being dried at 80 • C for 72 h. The chlorophyll pigments, namely Chl a, Chl b, and Chlt, were determined following the methods of Arnon [68] and Lichtenthaler and Wellburn [69] and calculated using the following Equations: where A is the absorbance at specific wavelengths, V is the final volume of extract (ml), and FW is the fresh weight of tissue extracted (g).

Yield and Yield Components
After the plants reached maturity on 14th April in both seasons, fifty randomly selected spikes from each plot were collected to determine SL, GWS, NGPS, and TGW. To determine BY and GY, five inner rows of 3 m (3.0 m 2 ) from each plot were harvested manually, sun-dried, and then weighed to determine BY in kg per 3.0 m 2 before being converted to ton ha −1 . All spikes of the harvested area were threshed and the grains were collected, cleaned, dried, and weighed to determine GY in kg per 3.0 m 2 and then converted to ton ha −1 . The HI, which is the ratio of GY and BY, was calculated using Equation (7), while the IWUE, which is the ratio of GY and ETc, was calculated using Equation (8) [70]: IWUE (kg ha mm −1 ) = GY ETc (8)

Data Analysis
An analysis of variance (ANOVA) was performed according to the split plot in a randomized complete block design by using the SAS statistical software 9.3 for Windows (SAS Institute, Cary, NC, USA). The mean values of treatments were compared using Duncan's test at a p ≤ 0.05 significance level. For a better understanding of the relationship between the traits across studied factors, Pearson's correlation coefficient and principle components analysis (PCA) were applied and performed using a computer software program XLSTAT statistical package (vers. 2019.1, Excel add-ins soft SARL, New York, NY, USA). A heatmap was performed using heatmap packages in R statistical software (version 4.2.2). The Pearson's correlation, PCA, and heatmap were made based on the pooled data of the two seasons. All figures were plotted using SigmaPlot 14 software program (v. 14.0; SPSS, Chicago, IL, USA).

Conclusions
The results of this study confirmed that it is difficult to grow wheat under deficit irrigation in arid climatic conditions without a significant decrease in its growth and production. Although the LM regime improved the IWUE, it markedly reduced the different morphophysiological traits, which ultimately lead to a wheat yield and yield components reduction as compared to the FL regime. Nevertheless, the co-application of SA, macro nutrients, and micronutrients through soil and foliar spray methods effectively mitigated the negative impacts of the LM regime, mainly by improving plant growth and physiological attributes. It also enhanced the performance of wheat under the FL regime. Therefore, strong and positive correlations among most studied parameters were observed under both irrigation regimes. In conclusion, combinations of osmolytes compounds, macronutrients, and micronutrients that are exogenously applied through foliar and soil methods could be recommended as an effective strategy to augment the growth and production of wheat crops under both full and deficit irrigation regimes in arid agro-ecosystems.