Assessing Soil Quality, Wheat Crop Yield, and Water Productivity under Condition of Deficit Irrigation

Wheat is one of the most important cereal crops in Egypt and all over the world. Its productivity is adversely affected by drought due to deficient irrigation to provide nutrients required for plant growth. In a field experiment, silicon foliar applications at concentrations of 0, 200, and 400 mg L−1 were performed at different irrigation rates ranging from 1000 to 4000 m3 ha−1 to assess water irrigation productivity and wheat crop yield in a calcareous soil under arid climate conditions. Increased irrigation rates led to a significant increase in soil nutrient dynamics, as well as in the number and weight of grains per spike, leaf area index, grain yield, straw yield, and biological yield, with the exception of the weight of 1000 grains. Spraying with sodium silicate had a significant impact on grain yield and harvest index but did not significantly impact the other traits. Furthermore, the interaction between irrigation and silicate application rates showed significance only for grain yield, the number of spikes/m2, and the harvest index. Applying three times irrigation could produce the highest nutrient retention, wheat yield, and water irrigation productivity. No significance was observed between 3000 m3 ha−1 and 4000 m3 ha−1 irrigation, indicating a saving of 25% of applied irrigation water. It can be concluded that applying irrigation at 3000 m3 ha−1 could be a supplemental irrigation strategy. High wheat grain yield can be achieved under deficit irrigation (3000 m3 ha−1) on the northwestern coast of Egypt with an arid climate by spraying crops with sodium silicate at a rate of 400 mg L−1.


Introduction
It is well known that Egypt is one of the top importers of basic crops, where 40% of its imports are grains, vegetable oils, and animal proteins due to the limited arable lands and water resources.Principally, it was reported as the top country for wheat imports in the world, with an expenditure of USD 2.69 million thousand in 2020, equivalent to 5.60% of the world's wheat imports, despite its production having reached 9.2 million tons of wheat in 2020 [1].
Consequently, there is an immediate necessity for effective management of agricultural production in Egypt to decrease land degradation and limited resources and to take a step forward toward self-sufficiency in agricultural productivity.Wheat (Triticum aestivum L.), family Poaceae, is known as Bread wheat [2] and is considered an excellent energy source due to its high contents of starch and protein [3,4].Wheat cultivation on a large scale is feasible, and its yield offers long-term food storage solutions.The growth of wheat occurs the soil C content was not enough to achieve the stoichiometric balance of soil nutrients.Moreover, it can be assumed that in all treatments, these soils can be classified as N-and P-limited according to their very low contents of organic C, N, and P causing an imbalance in soil nutrient dynamics.Nevertheless, a substantial enhancement in nutrient dynamics was observed, demonstrating high statistical significance (p < 0.05), across the irrigation periods, with the highest values being observed in those treatments receiving 4000 m 3 ha −1 compared with the other irrigation treatments.Therefore, we may assume that increasing the number of irrigation applications would result in enhanced soil nutrient cycling and consequently promote the stoichiometric balance of soil nutrients.Furthermore, an increase in the amount carbon dioxide (CO 2 ) released from the soil, referred to as its respiration capacity, was noted.The elevated CO 2 emissions from these soils with respect to their low contents of soil carbon pools can be attributed to the imbalance in soil nutrients that in turn increased the soil carbon loss assessed as C-CO 2 .No significant variability was observed in the CO 2 emitted from soil among the irrigation treatments at each sampling time, but high significance among the three sampling times was found for CO 2 (F = 45.46,p < 0.001), C-CO 2 (F = 45.46,p < 0.001), and C-CO 2 /SOC (F = 27.26,p < 0.001).This could imply that the significance of nutrient cycling found for molar C, N, and P did not contradict the carbon loss supplied as C/CO 2 , suggesting that nutritional cycling was essential despite the minor carbon loss.Accordingly, the nutrient contents varied significantly due to the changes in soil moisture regimes despite the fact that the applied irrigation rate was below the seasonal recommended water content of 5000 m 3 ha −1 (Table 1).
Soluble soil cations and anions were measured (Table 2).Moreover, Tukey's HSD test at α < 0.05 among the four irrigation treatments at the three sampling times showed highly significant variability in all soil cations but with no significant changes in soil anions.The cation exchange capacity (CEC) and the effective percentages of cations are presented in Table 3.The CEC of the soil and the effective percentages of calcium (Ca Effective ) and magnesium (Mg Effective ) presented a significant increase (p < 0.05) with the increase in the irrigation rate but no significant changes were observed.The correlation coefficients among all the studied soil parameters were measured (Figure 1).It can be found that soluble soil cations and anions decreased with the increase in soil organic reserve and nutrient limitation parameters.But the effective percentage of soil cations increased positively with the increase in contents of soil organic reserve parameters.Soil moisture was highly correlated with soil organic C and total N. Soil organic reserve parameters including SOC, TN, and DOC were positively correlated with the nutrient limitation parameters (C/P and N/P), soil moisture (SM), and soil carbon loss parameters (CO 2 and C-CO 2 ).Soil organic carbon was negatively correlated with DOC/SOC, indicating that increasing the irrigation rate may increase the transformation rate of dissolved or labile organic compounds into soil carbon pools (stable carbon forms).Additionally, soil carbon loss as a function of SOC (C-CO 2 /SOC) showed to be highly correlated with soil moisture and labile carbon fraction (DOC).This observation confirms that the dynamics of soil organic matter storage, as indicated by soil organics, exhibit nutrient-limiting increases with the increase in the irrigation rate.
CEC: cation exchange capacity; I1: one time irrigation; I2: two times irrigation; I3: three times irriga-tion; I4: four times irrigation.Different le ers indicate significant data variability; means with the same le er are not significantly different from each other at α < 0.05 as checked by Tukey's HSD test among all irrigation treatments.The same le ers indicate no significance.
The correlation coefficients among all the studied soil parameters were measured (Figure 1).It can be found that soluble soil cations and anions decreased with the increase in soil organic reserve and nutrient limitation parameters.But the effective percentage of soil cations increased positively with the increase in contents of soil organic reserve parameters.Soil moisture was highly correlated with soil organic C and total N. Soil organic reserve parameters including SOC, TN, and DOC were positively correlated with the nutrient limitation parameters (C/P and N/P), soil moisture (SM), and soil carbon loss parameters (CO2 and C-CO2).Soil organic carbon was negatively correlated with DOC/SOC, indicating that increasing the irrigation rate may increase the transformation rate of dissolved or labile organic compounds into soil carbon pools (stable carbon forms).Additionally, soil carbon loss as a function of SOC (C-CO2/SOC) showed to be highly correlated with soil moisture and labile carbon fraction (DOC).This observation confirms that the dynamics of soil organic ma er storage, as indicated by soil organics, exhibit nutrientlimiting increases with the increase in the irrigation rate.

Factor Analysis for Soil Response
Factor analysis was employed for all soil analyses (Table 4).The first three factors accounted for 85% of the total variance.The first factor, which explained 59% of the total variance, exhibited high positive loadings (>0.50) for electrical conductivity; soluble potassium, calcium, and magnesium; cation exchange capacity; and effective percentage of P and Mg.On the other hand, it displayed high negative loadings (>−0.50) for soil respiration (CO 2 ), carbon loss as C-CO 2 , C-CO 2 /SOC, and soluble and effective sodium.The structure of the second factor accounted for 17% of the total variance, displaying high positive loadings (>0.50) for soil pH, soil organic carbon, total nitrogen, and molar C/P and N/P ratios, as well as soluble and effective sodium.Conversely, high negative loadings (>−0.50) were observed for the dissolved proportion of SOC expressed as DOC/SOC, total P, bicarbonates, chlorine, SO 4 , and effective Ca.The third factor explained 9% of the total variance with high negative loadings (>−0.50) for soil moisture, soil organic carbon, dissolved organic carbon, total nitrogen, carbon loss, and C-CO 2 /SOC directly.
The distribution of factor loadings in the relationship between the first-two-factor structures is shown in Figure 2a.

Wheat Yield Assessment
The main parameters for assessing the wheat crop yield (Table 5), showed that grain number/spike, grain weight/spike, and the number of spikes/m 2 , as well as the LAI of wheat, grain yield, biological yield, straw yield, and harvest index, were significantly increased by increasing the number of irrigation applications.For example, grain/spike increased from 32.93 under I1 to 44.97 under the I4 treatment.Also, grain weight/spike increased from 1.04 g to 1.53 g under the I1 and I4 treatments, respectively.The LAI in-

Wheat Yield Assessment
The main parameters for assessing the wheat crop yield (Table 5), showed that grain number/spike, grain weight/spike, and the number of spikes/m 2 , as well as the LAI of wheat, grain yield, biological yield, straw yield, and harvest index, were significantly increased by increasing the number of irrigation applications.For example, grain/spike increased from 32.93 under I1 to 44.97 under the I4 treatment.Also, grain weight/spike increased from 1.04 g to 1.53 g under the I1 and I4 treatments, respectively.The LAI increased from 1.58 to 2.19 under the I1 and I4 treatments, respectively.Nonetheless, no significant differences at p < 0.05 were observed in the weight of 1000 grains.Grain yield, biological yield, and straw yield showed the highest significant increases (p < 0.05) under the two times irrigation (I2) treatment because of the soil's physical and chemical development at this irrigation rate.Under sodium silicate treatments, grain yield increased under the application of 200 mg L −1 , but HI it increased under 400 mg L −1 .The statistical interactions applied between irrigation and silicate treatments resulted in high significance only for spikes/m 2 , grain yield, and harvest index.Five effect size statistics were used to assess the relative importance of the studied factors.They were eta squared (η 2 ), partial eta squared (η 2 p ), omega squared (ω 2 ), partial omega squared (ω 2 p ), and epsilon squared (ε 2 ).

Grain Yield Assessment
Among the wheat yield assessment parameters, grain yield, harvest index, and number of spikes/m 2 were tested with Tukey's HSD test (Table 6) to check their significant response to irrigation and silicate treatments and their interactions.Highly significant differences at p < 0.05 were observed among all applied treatments for the three parameters.The data summarized that the lowest values of grain yield (0.95-ton ha −1 ), harvest index (0.11), and number of spikes/m 2 of the wheat crop (378.7) were for the interaction of one time irrigation with no Na-silicate.The four times irrigation treatment with 400 mg L −1 silicates (the highest application rate) yielded the highest values for both grain yield and harvest index.The highest values for the number of spikes were found under four times irrigation with no silicate.As the interaction between irrigation and silicate was significant for grain yield, harvest index, and number of spikes/m 2 , the response of these traits to irrigation depended on the silicate level, so high grain yield and harvest index were obtained under the four times irrigation treatment with 400 mg L −1 silicates.The correlation matrix depicted in Figure 3 illustrates the relationship between the number of spikes per square meter, the number of grains per spike, the weight of 1000 grains, and grain yield.The upper triangle represents the correlation coefficients of each pair of traits (parameters), while the lower triangle visualizes the trend relationship between the two variables.The diagonal histograms of the matrix show the distribution pattern of each trait.The matrix reveals that the traits were approximately normally distributed and that the number of grains/spikes was positively correlated with grain yield (r = 0.72, p < 0.01).The trend relationship of this correlation increased linearly along with the increase in the value of grain yield.A significant negative correlation was observed between the number of grains/spikes and the weight of 1000 grains (r = −0.37,p < 0.05).However, no significant correlation was observed for grain yield against the numbers of spikes/m 2 (r = 0.22) and the weight of 1000 grains (0.06).The trend relationship between the number of spikes/m 2 and grain yield was linear with a slight slope, which means a weak effect on grain yield.The trend relationship between the weight 1000 grains and grain yield formed a linear trend when the weight of 1000 grains increased by more than three grams.
Figure 4 illustrates the path diagram resulting from path analysis, revealing the direct effect of yield components on grain yield and their correlations.The analysis identified the number of grains per spike (0.87) as the most influential and crucial trait affecting grain yield.The number of spikes/m 2 was not adequate for grain yield, as its coefficient was (−0.05).Despite that, the number of grains/spike was positively correlated by the number of spikes/m 2 (0.31), while it was negatively affected by the weight of 1000 grains (−0.37), as there was no effect between the last two parameters (0.00).
Plants 2024, 13, 1462 9 of 20 between the number of grains/spikes and the weight of 1000 grains (r = −0.37,p < 0.05).However, no significant correlation was observed for grain yield against the numbers of spikes/m 2 (r = 0.22) and the weight of 1000 grains (0.06).The trend relationship between the number of spikes/m 2 and grain yield was linear with a slight slope, which means a weak effect on grain yield.The trend relationship between the weight 1000 grains and grain yield formed a linear trend when the weight of 1000 grains increased by more than three grams.Figure 4 illustrates the path diagram resulting from path analysis, revealing the direct effect of yield components on grain yield and their correlations.The analysis identified the number of grains per spike (0.87) as the most influential and crucial trait affecting grain yield.The number of spikes/m 2 was not adequate for grain yield, as its coefficient was (−0.05).Despite that, the number of grains/spike was positively correlated by the

Grain Yield Reduction
Table 7 presents the relative yield reduction affected by the number of irrigation applications at the silicon application rates.The data show the highest reduction under the one time irrigation treatment (an average of 53.7%) compared with the four times irrigation treatment.The highest reduction was under zero-silicate application in this treatment (70.9%).The lowest reduction was 6.4% under three times irrigation compared with the four times irrigation treatment, which also saved about 25% of the amount of irrigation.

Grain Yield Reduction
Table 7 presents the relative yield reduction affected by the number of irrigation applications at the silicon application rates.The data show the highest reduction under the one time irrigation treatment (an average of 53.7%) compared with the four times irrigation treatment.The highest reduction was under zero-silicate application in this treatment (70.9%).The lowest reduction was 6.4% under three times irrigation compared with the four times irrigation treatment, which also saved about 25% of the amount of irrigation.

Irrigation Water Productivity under Applied Treatments
The irrigation water productivity values are expressed as kg m −3 and concern the applied irrigation water (AIW) under irrigation and sodium silicate treatments (Table 8).Regarding the averages presented by the two observed seasons, the highest value of IWP was 2.1 kg m −3 under the one time irrigation treatment, while the lowest was 1.09 kg m −3 under the four times irrigation treatment.Within the treatments, the highest IWP was 2.98 kg m −3 under the one time irrigation treatment at the 400 mg L −1 silicate rate, while the lowest values was 0.815 kg m −3 under the four time irrigation treatment with zerosilicate application.

Discussion
The ideal ratios of C:N:P were previously presented as 106:16:1 as evaluated in the aquatic ecosystems [24].Soils with N/P ratios below 14 can be classified as nitrogen-limited, presenting very low contents of soil N forms [25].
Soil salinity and soil alkalinity were negatively correlated (r = −0.75,p < 0.01) across the irrigation treatments.Significant variation in soil nutrient dynamics was observed among the irrigation treatments, suggesting that the notable increase in soil water retention enhanced soil buffering capacity.This, in turn, led to heightened soil microbial activity and improved the cycling of soil nutrients in the upper soil horizons [26].
In the context of the factor analysis results, it becomes evident that changes in soil salinity and cation exchange capacity may exert a notable influence on soil respiration capacity and carbon loss.The observed negative loadings for soil respiration (CO 2 ), carbon loss as C-CO 2 , and C-CO 2 /SOC ratio in relation to factors associated with soil salinity and cation exchange capacity support this notion.
These findings suggest that increased soil salinity and cation exchange capacity may negatively impact soil respiration, leading to elevated carbon loss from soil.This could be attributed to several factors, including reduced microbial activity and organic matter decomposition in saline soils, which may inhibit the release of CO 2 through respiration processes.Additionally, the higher cation exchange capacity in these soils may result in increased adsorption of organic carbon, reducing its availability for microbial decomposition and subsequent CO 2 release.
The implications of these relationships are significant for soil management and carbon cycling in agricultural ecosystems.High soil salinity and cation exchange capacity could potentially compromise soil health and fertility by reducing microbial activity and organic matter decomposition rates.As a result, strategies aimed at mitigating soil salinity and optimizing cation exchange capacity may be crucial to maintaining soil respiration capacity and minimizing carbon loss [27].
The dataset shows a substantial correlation among the variables of soil pH, soil organic carbon (SOC), total nitrogen, molar C/P and N/P ratios, and soluble and effective sodium, as indicated by their high positive loadings (>0.50).This correlation could suggest common underlying mechanisms or environmental factors.For example, the presence of greater soil pH, SOC, and total nitrogen loadings may indicate that regions with elevated pH levels also exhibit increased levels of organic carbon and nitrogen.This correlation could be suggestive of more fertile or healthier soils.
In addition, the molar C/P and N/P ratios offer valuable information about nutrient cycling and stoichiometry in the soil system.Significant positive loadings for these ratios suggest a possible equilibrium or imbalance in the availability of carbon, phosphorus, and nitrogen, which could impact the dynamics of nutrients and the productivity of the ecosystem.
In contrast, variables such as the dissolved proportion of SOC represented as DOC/SOC, total phosphorus, bicarbonates, chlorine, SO4, and effective calcium showed significant negative loadings (>−0.50).The negative loadings indicate a negative correlation between these variables and the components with positive loadings.For instance, the presence of negative loading for total phosphorus and bicarbonates could suggest reduced quantities of these constituents in regions characterized by elevated soil pH and organic carbon content [28].
Furthermore, the adverse impact of effective calcium, bicarbonates, chlorine, and SO 4 loadings suggests the possibility of interactions or competition among these ions in the soil system.Gaining a comprehensive understanding of these connections is essential to deciphering the intricate dynamics of soil chemistry and its consequences for the functioning of ecosystems.
These findings indicate that irrigation can affect the processes involved in the movement of carbon, nitrogen, and phosphorus in soil by altering the levels of anions and alkalinity.These findings emphasize the significance of taking into account not only indi-vidual soil qualities but also their interconnections and potential drivers, such as irrigation techniques, in attempts to manage soil and conserve the environment.Additional investigation and analysis of these connections could offer significant knowledge for enhancing agricultural methods and reducing environmental consequences [29,30].
The third factor explained 9% of the total variance with high negative loadings (>−0.50)from soil moisture, soil organic carbon, dissolved organic carbon, total nitrogen, carbon loss, and C-CO 2 /SOC, directly indicating that the number of irrigation applications, when affecting soil moisture, directly influence the dynamic changes in soil carbon preservation and mineralization processes.
The highest communality values were for soil moisture, soil organic reserve parameters, soil carbon loss assessment parameters, and soluble cations, pointing out that these parameters were the most affected by the changes in irrigation rates.In dry areas, irrigation rates, which influence soil moisture regimes, tend to develop soil organic dynamics parameters and thus reduce soil carbon loss [31].
The number of irrigation applications plays a crucial role in the dynamics of soil organic matter for its storage and mineralization processes because of soil response to wetting/drying cycles.The highest contribution to the first two factors was from the four times irrigation treatment at each sampling time (Figure 2b).The increasing factor loadings of CEC and soil cations observed for the first factor, as well as soil organic C and N for the second factor, suggest the dynamic stability of soil organic matter contents (Figure 2a).
This enhancement in wheat crop yield due to increasing irrigation frequency aligns with the findings of the study in [32].
Grain yield, biological yield, and straw yield showed the highest significant increase (p < 0.05) under the three times irrigation (I3) treatment because of the soil's physical and chemical development at this irrigation rate.These results agree with those of [33], corroborating that silicon helps the plant to alleviate drought stress, thus enhancing crop productivity.Likewise, [34] revealed that foliar silicon application improved grain yield in drought-stressed wheat crops by further decreasing canopy temperature.It was also applied in priming seeds, which demonstrated an improvement in yield parameters in wheat crops, whose seeds were treated with silicon during seed preparation, which triggered water-deficit tolerance until plant maturity.
The effect size statistics revealed that the irrigation factor was more affected than the sodium silicate factor in all the studied parameters, even when the two factors were significant.This was obvious in grain yield and harvest index, where both irrigation and sodium silicate and their interaction were significantly affected.The most responsive parameter was grain yield, where the eta squared (η 2 ) values were 0.701 and 0.082, the partial eta squared (η 2 p ) values were 0.946 and 0.672, the omega squared (ω 2 ) values were 0.692 and 0.077, the partial omega squared (ω 2 p ) values were 0.885 and 0.461, and the epsilon squared (ε 2 ) values were 0.694 and 0.077 for irrigation and sodium silicate, respectively.
Additionally, the weight of 1000 grains (grain index) (0.38) showed significant influence on grain yield.These results are in agreement with the results of [35], who reported that a 20% reduction in irrigation requirements did not significantly affect plant height, spike length, the number of spikelets/spike, biological yield (tons/fed.), straw yield (tons/fed.), or harvest index.However, the number of grains per spike, grain weight per spike, and grain yield (tons per feddan) were significantly impacted.
Also, [36] reported that foliar potassium silicate rates improved drought tolerance in different growth stages by enhancing the morphological, physiological, and antioxidant potential of the wheat crop.
Silicon may be a low-cost way to reduce crop water stress.However, its effectiveness varies greatly, and the basic processes are unknown.Of concern is whether and how variability in plant silicon accumulation is related to its effects on plant growth during water stress.
Previous research has examined variations in silicon accumulation and their effects on plant growth and crop yield [37,38]. Regarding the relationship between different accumulation of silicon and the variations in stress tolerance, researchers investigated whether silicon affects osmotic and drought stresses differently in strains with higher silicon content compared with those with lower content.Some investigations indicated that silicon enhances wheat growth under osmotic and drought stress [2,10], but other studies did not find a significant increase in shoot dry weight biomass [20,39].Furthermore, in many other crops, such as soybean, growth has been reported to show no response to silicon under drought despite the reduced membrane damage and increased peroxidase activity.In barley, silicon reduced osmotic stress but did not affect biomass [18,19].These discrepancies could be partially clarified by methodological differences; for example, many research studies used sodium and potassium silicate as a silicon treatment but did not account for cation concentrations in the control treatment.As a result, the demonstrated silicon response was due to additional sodium and potassium.Another possibility of misunderstanding is the difference in interpretation; many researchers have reported that silicon positively affects tolerance to osmotic or drought stress, which can be explained when compared with control treatments with a higher response to the associated stress.Another significant component may be the genotype's silicon response.The authors of [40] realized that silicon positively affected growth in a poinsettia cultivar.Additionally, sugarcane response to silicon varied significantly in dry weight only in one of four cultivars under drought conditions [17].The present study states that foliar application of silicon may enhance plant resistance to water deficit.However, more detailed studies with more levels of both irrigation and silicon application should be conducted to achieve comprehensive conclusions about the role of silicon.

Study Area
The research site was situated in New Borg El-Arab City, within Alexandria Governorate, Egypt, approximately 10 km inland from the sea (see Figure 5).It is located at latitude 30.91 N, longitude 29.68 E, at an altitude of 54 m above sea level.Its location is representative of the arid climate of the northwestern coast of Egypt.

Study Area
The research site was situated in New Borg El-Arab City, within Alexandria Governorate, Egypt, approximately 10 km inland from the sea (see Figure 5).It is located at latitude 30.91 N, longitude 29.68 E, at an altitude of 54 m above sea level.Its location is representative of the arid climate of the northwestern coast of Egypt.
The arid climate of the area is typified by scorching, dry summers and mild winters, with annual precipitation ranging from 100 to 150 mm(see Figure 6).Short rainstorms sometimes occur, mainly in winter, most of which fall in January, February, and March.The Egyptian coastal zones receive noticeable amounts of rainfall, especially in winter; therefore, precipitation is considered the main source of groundwater recharge for aquifers in the area.Figure 5 represents the mean annual rainfall map during 2011-2020 according to the CRU TS (Climatic Research Unit gridded Time Series) dataset [41] and the World Imagery with Metadata web map (2023) as basemap Imagery with Metadata-Overview (arcgis.com).The field experiment spanned two wheat growing seasons, taking place from November to May in both the years 2021/2022 and 2022/2023.During the two seasons, the highest temperature was 37 ± 5 °C in May, while the lowest was 10 ± 4 °C in January.The average humidity varied between 60 and 70%.The average wind speed was around 14 km h −1 .The arid climate of the area is typified by scorching, dry summers and mild winters, with annual precipitation ranging from 100 to 150 mm (see Figure 6).Short rainstorms sometimes occur, mainly in winter, most of which fall in January, February, and March.The Egyptian coastal zones receive noticeable amounts of rainfall, especially in winter; therefore, precipitation is considered the main source of groundwater recharge for aquifers in the area.

Field Experiment
The field experiment was designed as a split-plot design, comprising three replications.Irrigation treatments were assigned to the main plots, whereas foliar silicon spraying was assigned to the sub-plots as subsequently explained.The irrigation treatments were applied as follows: (i) one time irrigation (I1) (1000 m 3 ha −1 ) at sowing in mid-November, (ii) two times irrigation (I2) (2000 m 3 ha −1 ) at planting and tillering, (iii) three times irrigation (I3) (3000 m 3 ha −1 ) at planting, tillering, and elongation, (iv) four times irrigation (I4) (4000 m 3 ha −1 ) at planting, tillering, elongation, and grain filling.Irrigation treatments were sequenced by accumulating the number of irrigation applications per 1000 m 3 ha −1 .Three different sodium silicate (Na2SiO3) (purchased from Al-Gomhoria company, Alexandria, Egypt) concentrations (0, 200, and 400 mg L −1 ) were applied one month after planting and repeated after two weeks.Each sub-plot area was delimited at 6 m 2 (2 × 3 m) and divided into six rows, separated by a 30 cm distance.Wheat grains cv.Giza 168 (pedigree: MIL/BUC//Seri CM93046-8M-0Y-0M-2Y-0B-0SH) were sourced from the Agriculture Research Center, Giza, Egypt, and were manually drilled into the rows in mid-November

Field Experiment
The field experiment was designed as a split-plot design, comprising three replications.Irrigation treatments were assigned to the main plots, whereas foliar silicon spraying was assigned to the sub-plots as subsequently explained.The irrigation treatments were applied as follows: (i) one time irrigation (I1) (1000 m 3 ha −1 ) at sowing in mid-November, (ii) two times irrigation (I2) (2000 m 3 ha −1 ) at planting and tillering, (iii) three times irrigation (I3) (3000 m 3 ha −1 ) at planting, tillering, and elongation, (iv) four times irrigation (I4) (4000 m 3 ha −1 ) at planting, tillering, elongation, and grain filling.Irrigation treatments were sequenced by accumulating the number of irrigation applications per 1000 m 3 ha −1 .
Three different sodium silicate (Na 2 SiO 3 ) (purchased from Al-Gomhoria company, Alexandria, Egypt) concentrations (0, 200, and 400 mg L −1 ) were applied one month after planting and repeated after two weeks.Each sub-plot area was delimited at 6 m 2 (2 × 3 m) and divided into six rows, separated by a 30 cm distance.Wheat grains cv.Giza 168 (pedigree: MIL/BUC//Seri CM93046-8M-0Y-0M-2Y-0B-0SH) were sourced from the Agriculture Research Center, Giza, Egypt, and were manually drilled into the rows in mid-November for both seasons at the rate of 143 kg ha −1 .The planting was carried out at 1 cm soil depth in mid-November, and harvest was in late April for each season.The experimental fields were fertilized with 357 kg ha −1 phosphorus (15.5% P 2 O 5 ), during seedbed preparation, mixed thoroughly into the soil before sowing, and then gently raked to a depth of 10-15 cm.Nitrogen fertilizer (ammonium nitrate of 33.5% N) was administered at a rate of 238 kg/ha in three equal doses at planting, tillering, elongation, and during rainfall.

Soil Quality Assessment
Soil samples were collected in triplicate from each irrigation plot, with sampling being conducted three times at two-month intervals.Fresh soil cores were collected from each plot, in triplicates, to determine soil moisture contents by using the gravimetric method.The collected samples were oven-dried at 105 • C for 24 h.Initial soil characteristics were analyzed before the experiment (Table 9).The soil samples were air-dried and then ground and sieved through a 2 mm mesh filter for subsequent soil analyses [42].By using Robinson's pipette method, soil particle analysis was performed (Eijkelkamp Agri research Equipment, Giesbeek, The Netherlands).Based on the soil particle size analysis, the soil texture was classified as sandy clay loam (62.43%), silt (13.02%), and clay (24.55%).The total CaCO 3 content was quantified by CO 2 release from fine soil upon acidification with 10% HCl.Soil pH and electrical conductivity (EC) were measured in 1:2.5 w/v and 1:1 w/v soil aqueous extracts, respectively.Soluble cations (Na + , K + , Ca 2+ , and Mg 2+ ) and anions (SO 4 2− , Cl − , and HCO 3 − ) were determined in the 1:1 soil aqueous extract.Soluble Na + and K + were estimated (flame photometer, PG Instruments M: FP902).The EDTA titration method measured soluble Ca 2+ and Mg 2+ in the extract.Soluble SO 4 2− was determined by barium sulfate precipitation.Soluble Cl − and HCO 3 − were determined by titration with silver nitrate and 0.01 N H 2 SO 4 , respectively.Soil organic carbon (SOC) was assessed via the Walkley-Black wet oxidation method.Dissolved organic carbon (DOC) was extracted from a 1:5 (w/v) suspension of deionized water at 25 • C. Following a 24 h shaking period, the suspension underwent filtration through a cellulose nitrate membrane filter with a pore size of 0.45 µm.DOC, representing the dissolved fraction of soil organic matter (SOM) in the resulting aqueous extract, was quantified by using a total organic carbon (TOC) analyzer (Torch Combustion TOC/TN Analyzer; Teledyne Tekmar, OHIO, USA) [43].The dissolved organic carbon fraction was expressed as a percent (DOC/SOC; %) of total SOC content.Total nitrogen (TN) was measured by using the Kjeldahl method.Available phosphorus was extracted with 0.5 M NaHCO 3 at pH 8.5 according to [44] and measured by the PG Instruments T80 UV/VIS Spectrophotometer.The stoichiometric balance of soil C:N:P ratio, known as the "Redfield ratio", was also calculated considering the molar ratios of C, N, and P [25].The cation exchange capacity (CEC) was also assessed, including the effective percentage of each cation concentration in cmol kg −1 , presented as Na effective , K effective , Ca effective , and Mg effective .
The soda lime absorption method was applied for soil respiration measurements (the metabolically emitted CO 2 from the soil) by using the dark cover box technique in the field.The detected CO 2 concentrations were calculated after 24 h of incubation of soda lime on the soil surface in the delimited area under the cover boxes installed on the soil surface.
The calculations were derived after the oven-drying of soda lime, followed by multiplication by a correction factor of 1.69, and presented in mg g −1 soil.The soil carbon loss emitted into the atmosphere as C-CO 2 was calculated by multiplying the CO 2 readings by 12/44 (the molecular weight of C/CO 2 ) and then expressed in mg g −1 .

Wheat Crop Yield Assessment
Plants were harvested and randomly chosen from each sub-plot to assess wheat yield measurements by measuring the LAI by using the disc method [45] the number of grains/spike and grain weight/spike (g) by taking the average of 10 random spikes from each plot; the number of spikes m −2 by counting the spikes in a random square meter in each plot; grain yield (ton ha −1 ) by weighing the grains threshed from the plot and calculating the yield for hectare; straw yield (ton ha −1 ) as the difference between biological yield and grain yield; biological yield (ton ha −1 ) by weighing all the plants from the harvested plot; 1000 grain wt.(gm); and the harvest index (HI) by dividing the grain yield by the biological yield.
Ref. [46] linked the relative yield decrease with the relative evapotranspiration deficit by using the following equation: where ET a represents the rate of actual evapotranspiration, quantifying the water stress in plants (actual water consumption; mm), and ET m denotes the rate of maximum evapotranspiration (maximum water consumption; mm).The relative yield, Yr, is calculated as Y a /Y m , where Y a is the actual yield and Y m is the maximum yield.When there is any reduction in irrigation, Y a equals Y m .K y represents the yield response factor.The recommended amount of applied irrigation water was calculated as described by [47] as follows: where Ea represents the irrigation efficiency (0.60 for surface irrigation).

Crop Evapotranspiration (ETc)
The FAO crop evapotranspiration method was used to calculate crop water requirements (Table 2) by using reference evapotranspiration (ETo; mm) and crop coefficient (Kc) as described by [48]: The ETo was calculated by using CROPWAT software based on the FAO Penman-Monteith equation [49].

Irrigation Water Productivity (IWP)
IWP (irrigation water productivity; kg m −3 ) calculations followed [50] method as follows: where Y a is the seed yield of various treatments (kg ha −1 ) and AIW is the seasonal amount of applied water (m 3 ha −1 ), which was calculated by using the following equation: where AIW denotes the applied irrigation water, Peff denotes adequate rainfall during the growing season, and S denotes groundwater contribution to crop water use (neglected because the water table was deeper than 2.0 m) [51].

Statistical Analysis
Following the assessment of error homogeneity through Bartlett's test, a combined analysis encompassing both seasons was undertaken.Soil quality parameters together with wheat yield assessment parameters were compared among all applied treatments by using Tukey's HSD (Honestly Significant Difference) tests at the p < 0.05 level, following ANOVA performed by using STATISTICA 10 by StatSoft, Inc. (Tulsa, OK, USA) [52].The statistical analysis of wheat and grain yield assessment parameters was conducted by using R Statistical Software (R_Core_Team, 2021).

Conclusions
The results confirm that three times irrigation application amounting to 3000 m 3 ha −1 for wheat in these soils as proposed can lead to the highest nutrient retention, wheat yield, and water productivity.This irrigation rate showed adequate soil carbon storage capacity against carbon loss that was highly dependent on the solubility of soil cations and anions, thus affecting soil microbial activity.The results also show that the changes in soil organic matter dynamics in preservation and mineralization across the different irrigation times in these calcareous soils varied due to soil salinity and alkalinity and soil cation exchange capacity.Further research is needed to identify the soil organic matter fractions more susceptible to the changes in soil moisture regimes in these calcareous soils with poor structure.
The three times irrigation rate represents 60% of the recommended total applied irrigation water (5000 m 3 ha −1 ) during the wheat season and produced a high grain yield and high water productivity on the northwest coast of Egypt.Moreover, foliar application of silicon can alleviate the adverse effect of water deficit in wheat production.Finally, a high yield of wheat grains can be achieved under deficit irrigation on the northwest coast of Egypt with a dry climate by spraying 400 mg L −1 sodium silicate.Deficit irrigation can be substituted by foliar silicon application, which can be, in addition, considered for climate-smart plants' long-term productivity under drought stress conditions.Further investigation is required to comprehend the dynamic mechanisms of silicate in different crops under diverse climatic conditions, which influence irrigation regimes.

Figure 1 .
Figure 1.Pearson correlation matrix among all the studied soil parameters.

Figure 1 .
Figure 1.Pearson correlation matrix among all the studied soil parameters.
Figure 2b also displays the distribution of factor scores in the relationship between the first two factors.Plants 2024, 13, x FOR PEER REVIEW 7 of 21

Figure 2 .
Figure 2. Relations between the first two factors: (a) factor loadings and (b) factor score values, red rectangle in (b) represents sampling times for the four-times irrigation.

Figure 2 .
Figure 2. Relations between the first two factors: (a) factor loadings and (b) factor score values, red rectangle in (b) represents sampling times for the four-times irrigation.

Figure 4 .
Figure 4. Path analysis of grain yield as affected by the number of spikes/m 2 , number of grains/spike, and weight of 1000 grains and their interaction effects.Green arrows indicated positive effect; red color indicates negative effect.Single headed arrows indicate direct effect, double headed arrows indicate correlation.

Figure 4 .
Figure 4. Path analysis of grain yield as affected by the number of spikes/m 2 , number of grains/spike, and weight of 1000 grains and their interaction effects.Green arrows indicated positive effect; red color indicates negative effect.Single headed arrows indicate direct effect, double headed arrows indicate correlation.

Figure 5 .
Figure 5. Annual rainfall (mm) in the study area.

Figure 5 .
Figure 5. Annual rainfall (mm) in the study area.
Figure 5 represents the mean annual rainfall map during 2011-2020 according to the CRU TS (Climatic Research Unit gridded Time Series) dataset [41] and the World Imagery with Metadata web map (2023) as basemap Imagery with Metadata-Overview (https://www.arcgis.com/).The field experiment spanned two wheat growing seasons, taking place from November to May in both the years 2021/2022 and 2022/2023.During the two seasons, the highest temperature was 37 ± 5 • C in May, while the lowest was 10 ± 4 • C in January.The average humidity varied between 60 and 70%.The average wind speed was around 14 km h −1 .Plants 2024, 13, x FOR PEER REVIEW 15 of 21

Table 1 .
General soil characteristics, as well as CO 2 emissions from soil and soil carbon loss.
EC: electrical conductivity; SOC: soil organic carbon; DOC: dissolved organic carbon; I1: one time irrigation; I2: two times irrigation; I3: three times irrigation; I4: four times irrigation.Different letters indicate significant data variability; means with the same letter are not significantly different from each other at α < 0.05 as checked by Tukey's HSD test among all irrigation treatments.The same letters indicate no significance.TN: total nitrogen; TP: total phosphorus; CO 2 : emitted CO 2 from the soil as soil respiration; C-CO 2 : soil carbon loss as C-CO 2 .

Table 2 .
Soluble cations and anions in soil under irrigation treatments.
I1: one time irrigation; I2: two times irrigation; I3: three times irrigation; I4: four times irrigation.Different letters indicate significant data variability; means with the same letter are not significantly different from each other at α < 0.05 as checked by Tukey's HSD test among all irrigation treatments.The same letters indicate no significance.

Table 4 .
First three-factor structures with communality values of factor analysis.

Table 5 .
Wheat yield assessment under irrigation and sodium silicate treatments and their interaction with the five effect size statistics tested by using R-2021software.

Table 6 .
Grain yield, harvest index, and number of spikes/m 2 of wheat under irrigation and silicate treatments.
I1: one time irrigation; I2: two times irrigation; I3: three times irrigation; I4: four times irrigation.Different letters (a-e) in columns indicate significant differences as tested by Tukey's HSD test at >0.05.

Table 7 .
Grain relative yield reduction (RYD) of application treatments.

Table 8 .
Irrigation water productivity (IWP) concerning applied irrigation water (AIW) under applied irrigation and sodium silicate treatments.

Table 9 .
The recommended amounts of applied irrigation water (AIW).) and applied irrigation water values are expressed in mm and m 3 ha −1 .