Effect of Foliar Application of Micronutrients and Fungicides on the Nitrogen Use Efﬁciency in Winter Wheat

: The increase in the grain density (GD) is the key factor for effectively controlling Nitrogen-Use Efﬁciency (NUE) in winter wheat. Winter-wheat protection with fungicides and/or foliar fertilization with micronutrients during the critical stages of yield formation affects the grain yield by increasing GD. This hypothesis was veriﬁed in a two-factor ﬁeld experiment, carried out in the 2013/2014, 2014/2015, 2015/2016 growing seasons. A ﬁeld experiment in a two-factor split-plot design, included three systems of wheat foliar protection (FP): (i) N + Mi (N + macronutrients), (ii) N + P (N + fungicides), (iii) N + Mi + P, and 0, 40, 80, 120, 160, 200, 240 kg N ha − 1 . The grain yield, despite the signiﬁcant effect of years, depended on the interaction of FP and N doses. The maximum yield of 8.1, 10.7, 11.1 t ha − 1 for the optimal N dose of 79.8, 227.4, and 245.2 kg ha − 1 , was achieved, respectively, for N + Mi, N + N, N + Mi + P. The wheat-grain yield depended signiﬁcantly on GD (R 2 = 0.98, p ≤ 0.001). Wheat protected with fungicides increased GD gradually with the N dose, reaching over 3000 grains per m 2 when fertilized with 200 kg N ha − 1 . The lack of fungicide protection, reducing GD, consequently resulted in a greater accumulation of N in straw at the expense of grain. Fungicide protection of wheat through the positive effect on N management by wheat should be considered as an agronomic measure that supports N-Use Efﬁciency.


Introduction
Wheat (Triticum aestivum L.) is one of the three main cereals responsible for bridging the food gap in the world [1]. It has been well-documented that differences in yields between regions of the world are due to biophysical factors such as temperature, rainfall amount and their distribution during the growing season [2,3]. The importance I of soil fertility cannot be ignored in food production [4]. In wheat, nutrient deficiency reduces nitrogen (N) uptake and, consequently, grain yield and its quality [5]. Modern wheat varieties have achieved a high yield potential in recent decades. However, their sensitivity to environmental factors has increased. This negative trend is especially visible in Europe [6]. The main factors, significantly affecting the realization of the yield potential of each crop are water and N, which are not complementary to each other [7,8].
The accumulation of N by wheat during the growing season occurs according to the so-called sigmoid pattern [9]. The main feature of this general model is the mathematically described pattern of N-accumulation during the successive stages of wheat growth. On this wheat) and Puccina triticina (leaf rust). On the opposite side are plants deficient in N due to (i) N shortage, and/or (ii) late stage of development. Plants deficient in N are attacked by necrotrophs, for example Pyrenophora tritici-repentis (tan spot) [35,36]. Wheat, in its life cycle, and during the growing season undergoes various nutritional states, which makes it constantly susceptible to pathogens [37]. The effective use of the yield potential of modern wheat varieties requires measures and actions that allow for an effective control of (i) the rate of the applied fertilizer N, (ii) the nutritional status of the plant in critical stages of grain yield formation, and (iii) the pressure of pathogens. The implementation of all these measures together is the only possible way to increase the Nitrogen-Use Efficiency (NUE).
The present state of knowledge assumes that the increase in grain density is a key factor that effectively control the NUE of winter wheat. The objective of the study was to evaluate the effect foliar protection of winter wheat with micronutrients and/or fungicides on the development of yield components and, finally, on grain yield as a prerequisite for effective NUE control.

Experimental Setup
The study on the grain yield of winter wheat, N accumulation, N use efficiency (NUE) indices in response to the management system (foliar fertilization with micronutrients, plant protection-fungicides), was carried out in 2013/14, 2014/2015, and 2015/2016 seasons in Smolice (52 • 42 N; 17 • 10 E), Poland. The field experiment was carried out on soil formed from loamy sand over loamy sand, classified as Albic Luvisol. The organic carbon (C org ) content and pH values were variable, peaking in the 2014/2015 growing season. The content of available nutrients, measured before wheat sowing, i.e., before the application of fertilizers, was in the medium class at least for P, K, Fe, and Zn. This means very favorable conditions for wheat growth. More variable levels of soil fertility were found for other nutrients, such as Mg, Ca, Cu, and Mn, which in the 2013/2014 growing season were in low or very low classes. The amount of mineral N (N min ), measured just before the spring regrowth of winter wheat in a 0.0-0.6 m layer, was generally high or very high as in 2015 (Table 1). The local climate of the study area, classified as intermediate between Atlantic and Continental, is seasonal, especially in Summer. The meteorological data are presented in Figure 1. The growth conditions of winter wheat assessed on the basis of the Sielianinov hydrothermal indices in the subsequent years of study were very diverse ( Table 2). The early vegetation phase of winter wheat in spring 2014 was very good, but May was wet. The second part of the season was less favorable as June was very dry and July was dry. The 2015 growing season was generally unfavorable for wheat growth, as indicated by the predominance of dry conditions. The beginning of the 2016 growing season was wet.
May was semi-dry, June was dry. These two months are critical for the development of the number of grains in an ear of winter wheat. Most of the wheat-grain filling period was wet again.

Experimental Design
The field experiment, arranged in a two-factor split-block design, was replicated four times, included:

1.
Three systems of wheat foliar protection in the spring part of the growing season (FP): 1.1. N + foliar fertilization with macronutrients (N + Mi, fungicide control − FC) 1.2.
The total area of one plot was 22.5 m 2 (1.5 × 15 m). The winter wheat cv. Wydma was sown annually in the fourth week of September in the rate of 300 grain m −2 . The forecrop was winter oilseed rape. The wheat was harvested the following year at the end of July from an area of 19.5 m −2 . Nitrogen was applied in the form of ammonia nitrate (34:0:0) in accordance with the experimental schedule: (1) 80 kg N ha −1 -the late winter, before beginning of winter wheat vegetation in spring; (2) 160 kg N ha −1 -at the end of tillering/beginning of shoot elongation (BBCH 29/30); (3) 240 kg N ha −1 -at the stage of a flag leaf becoming visible (BBCH 39).
Phosphorus was applied at the rate of 17.2 kg P ha −1 in the form of triple superphosphate (46% P 2 O 5 ). Potassium was applied at the rate of 100 kg K ha −1 as Korn-Kali (K-MgO-Na 2 O-SO 3 → 40-6-3-12.5). Both fertilizers were applied two weeks before wheat sowing. The foliar application of micronutrients and fungicides was carried out in accordance with the experimental schedule, as shown in Table 3.

Plant Sampling
The plant material for the determination of dry matter and the N content was collected at harvest. A single sample was first partitioned into grain and vegetative biomass (harvest residues = straw + chaffs). The N content was determined in both parts of the plant, using the standard macro-Kjeldahl procedure [42]. The content of N was expressed on a dry matter basis.

Regression Models
The general pattern of wheat grain yield (GY) and nitrogen accumulation (GN, TN) in response to increasing N doses (N fi ) was described by two regression models [43]:

2.
Quadratic: The main estimated parameters of these function are: The computational procedures for determining the models and parameters of GN and TN were the same. where: GY-grain yield, t ha −1 ; GN-N accumulated in grain, kg ha −1 ; TN-N accumulated in total wheat biomass, kg ha −1 ; N fi -"i" dose of fertilizer N, kg ha −1 ; N op -optimal N rate, kg ha −1 ; GY max -maximum GY for the N op , t ha −1 ; A, B, C, a, b, c-constants of the respective regression model.

Statistical Analysis
The effects of individual research factors (year, foliar protection, N doses) and their interactions on the grain yield and indices of nitrogen use efficiency were assessed by means of a two-way ANOVA. Means were separated by honest significant difference (HSD) using Tukey's method, when the F-test indicated significant factorial effects at the level of p < 0.05. The relationships between the traits were analyzed using the Pearson correlation and linear regression. The stepwise regression was applied to define the optimal set of variables for a given plant characteristics. The best regression model was selected based on the highest F-value for the entire model. STATISTICA 12 software was used for all statistical analyses (StatSoft Inc., Krakow, Poland, 2013).

Grain Yield and Grain Yield Structure
The experimental factors and the years had a significant impact on the wheat-grain yield ( Table 4). The highest GY recorded in dry 2015 was 66% and 22% higher, respectively, than in 2016 and 2014. The effect of foliar protection (FP) of wheat was significant. The yield increased in the following order: N + Mi < N + P < N + Mi + P. The increasing dose of N caused a progressive increase in the grain yield, averaged over the years and FP. The highest value was obtained on the plot fertilized with 240 kg ha −1 . Similar letters means a lack of significant differences using Tukey's' test; ***, **, * indicate significant differences at p < 0.001, p < 0.01, and p < 0.05, respectively; ns-non significant; 1 GY-grain yield; SY-straw yield; TB-total biomass; HI-harvest index; ED-ears density; GR-grain per ear; GD-grain density; TGW-thousand grain weight.
The wheat grain yield, despite the significant effect of experimental factors and years, depended on the interaction of FP and N doses, which did not interact with years ( Figure 2). The yield on the N + Mi object increased significantly to the N dose of 80 kg ha −1 , and then continued its growth, but at a much slower rate. The resulting N dose-response pattern can be described by a double linear model, indicating the existence of two sub-phases: (1) progressive: (2) smooth: Response of winter wheat grain yield to interaction of foliar protection systems and N rates. * Similar letters means a lack of significant differences using Tukeys' test. Legend: N + Mi-N and micronutrients; N + P-N + fungicides; N + Mi + P-N plus micronutrients and fungicides.
The intersection of both lines provided a peak, indicating a change in the actual rate of increase in the grain yield in response to the dose of N. As indicated by the value of the direction coefficient, each dose of N fertilizer above 78.9 kg ha −1 resulted in a 10-fold reduction in the rate of increase in the yield compared to its rate below this value. For the N + P and N + P + Mi objects, the grain yield response to gradually increasing N doses can be best described by the quadratic regression model: (1) N + P: GY = 0.000082N 2 + 0.037N + 6.432 for n = 7, R 2 = 0.99, P ≤ 0.01 (2) N + P + Mi: GY = 0.000074N 2 + 0.036N + 6.611 for n = 7, R 2 = 0.99, P ≤ 0.01 Both models are only seemingly the same, as indicated by the parameters of the equations obtained. The optimal N dose of 227.4 and 245.2 kg N ha −1 resulted in the maximum grain yield of 10.674 and 11.063 t ha −1 for N + P and N + P + Mi, respectively.
The straw yield, despite significant variability from year to year and for N doses, did not show a significant response to interactions of the studied factors (Table 4). The effect of the years was the same as that reported for the grain yield, but the differences were much smaller. The straw yield in 2015 was 21% higher than 2016 and 33.8% higher than to the lowest, yield recorded in 2014. The significant difference between the years resulted from the specific tendency of the harvest index (HI). It was significantly higher in 2014 and on the plots treated with fungicides. HI in 2014 was more than 10-percentage points higher than in 2016. This index showed a positive and significant response to plant protection and its combination with the Mi application. The gradual increase in N doses caused a progressive increase in HI up to 120 kg N ha −1 , where it reached a stable level of 48%. The overall response of total wheat biomass to the experimental factors and the years was consistent with the patterns observed for the grain yield.
The wheat grain yield significantly depended on the number of grains per unit area (Grain density, GD) (Table A1). GD variability explains 98% of the GY variability: GY = 0.52GD − 0.69 for n = 21, R 2 = 0.98 and P ≤ 0.01 (5) This component of the yield-grain structure also showed a significant dependence on the interaction of both experimental factors. At the same time, it showed no interaction with the years ( Table 4). The primary GD components, i.e., ear density (ED) and the number of grains per ear (GR), did not show a significant response to the interaction of experimental factors and years. ED changed significantly over the study years, showing the same trend as that observed for the grain yield, but the differences between the years were much smaller. This yield component progressively responded to the gradual increase in the N dose, achieving stability on the plot fertilized with 200 kg N ha −1 . The seasonal variability in the GR trend was almost at the same level as that for the grain yield. In 2015, it was 62% and 12.5% higher compared to 2016 and 2014. The combined action of micronutrients and fungicides (Mi + P) resulted in a significant increase in GR compared to the fungicide control. The effect of N doses was very specific, because GR increased progressively only up to the N dose of 80 kg ha −1 , and then stabilized. The only thousand grain weight (TGW) responded significantly to experimental factors and years. Nevertheless, its variability was low, as indicated by the means of 48.1 ± 2.54 g.

Nitrogen Accumulation
The amount of N accumulated in wheat grain (NG) showed the same tendency as for the grain yield in response to the experimental factors and the years (Tables 4 and 5). The main difference between the trends obtained resulted from the effect of N doses. As shown in Figure 3, the increase in NG almost linearly corresponded to the increasing doses of N. The highest value of NG was achieved on the plot fertilized with 240 kg N ha −1 . Its average value for other objects and years was 2.3-fold higher than in the N control. The key differences, as shown in Figure 3, resulted from the value of the direction coefficient. The FP effect was significant, and NG increased at a different rate in relation to the gradually increasing N doses. Averaged over the years and N doses, the order of the objects is as follows: N + Mi < N + P ≤ N + Mi + P.
The accumulation of N in the harvest residues (N r ) of wheat did not show a significant response to the interaction of all factors. The effect of the years was recorded only for 2014, when N r was 1 / 3 lower compared to other years. The effect of N doses, averaged over other factors, was only slightly smaller than that reported for NG. The highest dose of N, equal to 240 kg ha −1 resulted in double the amount of N in the crop residues compared to the N control. The total N accumulated in wheat (TN) was consistent with the general trend described for total wheat biomass (Tables 4 and 5). The Nitrogen Harvest Index (NHI) responded significantly to FP systems and N doses, but not to the interaction of experimental factors and years.
The overall effect of the experimental factors and years on the N utilization by wheat was assessed using three indicators. The first is Nitrogen Unit Accumulation (NUA). This index showed a significant response to all factors, analyzed separately, but not to the interaction of FP systems with years. The highest value of this index, obtained in 2016, was 30% and 40.6% higher than in 2015 and 2014. The FP effect was significant, but small. Higher NUA values were recorded for the N + Mi object. The effect of N doses was significant, but significantly lower than that observed for NG or N r . The highest value of NUA was recorded on the plot fertilized with 240 kg N ha −1 , which was 33% higher than in the N control plot.
The value of the next index, i.e., Total Unit Productivity of N (NUP-T) was the lowest in 2016 and the highest in 2014. As expected, the FP effect was the highest for the N + Mi + P object. Increased doses of N fertilizer caused a corresponding decrease in the NUP-T, which was the lowest for the plot fertilized with 240 kg N ha −1 . The same trend was observed for the NUP-G index, which expresses the productivity of the N unit accumulated in wheat grain. The effect of the years was the same as that recorded for the NUP-T, but at a higher level. No significant impact of FP on this index was observed. The effect of the N doses was significant, showing an interaction with years. As presented in Figure A1, NUP-G decreased linearly each year, which logically corresponds to the effect of the increased N dose. The main difference between the years is in the value of the direction coefficient. In 2016, this coefficient was 43% lower than in 2015, when it was the highest. The highest values of NUP-G recorded in the N control were, respectively, 58.   Similar letters means a lack of significant differences using Tukeys' test; ***, **, * indicate significant differences at p < 0.001, p < 0.01, and p < 0.05, respectively; ns-non significant; Legend: 1 NG, Nr, NT-nitrogen accumulation in grain harvest residues, total uptake; NHI-nitrogen harvest index; NUA-nitrogen unit accumulation; NJP-G-nitrogen unit productivity-grain; nitrogen unit productivity-total biomass.
The analysis of the relationships between the NUE indices and the yield components clearly revealed the dominant role of ED, which significantly influenced the NG (Table A2). The effect of NT and N r on the other yield components was significant, but much weaker. GR, as the basic component of GD, showed the highest positive relationship with NHI and at the same time negative with NUA. GD, as an aggregate component of the yield, showed the greatest correlation with NG, but not with N r . The second major grain yield component, i.e., TGW, was significantly, but weakly, correlated with the NUE indices, but mostly with N r .

Nitrogen Efficiency Indices
Partial Factor Productivity of fertilizer N (PFP-N) was the only index that responded significantly to the interaction of experimental factors, but not to their interaction with years ( Table 6). The values of this index showed high variability, which was mainly affected by the N dose ( Figure 4). The PFP-N response to increasing doses of N was the best fit to the power regression model. The key difference between the FP systems was the value of the constant, related to the theoretical efficiency of 1 kg of applied N fertilizer. This value was the highest for the N + Mi object and significantly lower for the other two. In the experimental reality, the PFP-N values for the plot with 40 kg N ha −1 ranged from 182 (100%) to 201 (104%) and to 206 (113%) kg grain kg −1 N for N + Mi, N + P and N + Mi + P, respectively. The lowest values of this index were typically found for the N treatment with the highest N dose. The differences between the main objects relatively increased on the fertilized plot with 240 kg N ha −1 , and amounted to 37 (100%), 45 (121.6%), and 47 (127%) kg grain kg −1 N, respectively. Similar letters means a lack of significant differences using Tukeys' test; ***, * indicate significant differences at p < 0.001, p < 0.05, respectively; ns-non significant; NAE-net nitrogen fertilizer productivity; R-N-nitrogen fertilizer recovery N-PhE-physiological efficiency of fertilizer taken up by the crop. Nitrogen Agronomic Efficiency (NAE) responded significantly to experimental factors and years. However, significant interactions were noted only for Y × P and Y × N (Table 6). This index was more sensitive to weather variability in the studied years than grain yield and N accumulation by wheat. Its highest value, recorded in 2015, was 3.35-fold and 1.5-fold higher than in 2016 and 2014, respectively. It also increased in response to FP systems, and was significantly higher in objects treated with fungicides. The NAE response to gradually increasing N doses was consistent with the classical rule, i.e., it decreased gradually in response to increasing N doses. A progressive decrease in the index was recorded up to an N dose of 160 kg ha −1 , and then it stabilized. A key indicator of the activity of N is its recovery from the applied fertilizer (R-N). This index followed the same patterns as those for NAE, but the differences were much weaker. The R-N index in 2015, as shown by the analysis of the Y × N interaction for all N doses exceeded 100%. Its values ranged from 181% to 101% on plots fertilized with 40 and 240 kg N ha −1 , respectively. In 2014, R-N ranged from 128% to 66%, and in 2016 only from 71% to 60%, respectively. The Physiological efficiency of N fertilizer (PhE-N) showed slightly different patterns of response to the examined factors, mainly due to the impact of years. In 2016, PhE-N values were 37% lower compared to the first two years, but followed the same trend in response to increasing N doses.
The examined NUE indices showed a negative correlation with the components of yield (Table A3). The highest values of the correlation coefficients were found for ED, and then for GD. It should be emphasized that the negative sign for the analyzed dependencies is misleading. In fact, it reports that these two groups of wheat characteristics are opposite to one other. A classic example is the relationship between R-N and ED: ED = −1.64R-N + 831.5 for n = 21, R 2 = 0.70 and P ≤ 0.05 (6) R-N = 0.0042ED 2 − 6.077ED + 2277.8 for n = 21, R 2 = 0.70 and P ≤ 0.01 (7) The first equation clearly shows that the lowest ED was associated with the highest R-N and vice versa. The second equation states that a minimum R-N of 77.8% was achieved for ED of 723 ears per m 2 . Each decrease in ED resulted in a significant increase in the value of the R-N index. No significant correlations were found between the NUE indices and TGW.

Grain Yield and Components of Yield Structure
The average wheat yields in the studied years were variable and ranged from 6.7 t ha −1 in 2016 to 11 t ha −1 in 2015. These two figures unequivocally confirm the hypothesis that modern wheat varieties are very sensitive to harsh conditions in the growing season, regardless of the European region [6,46]. Nevertheless, the interactive effects of the applied agronomic measures such as micronutrients and fungicides applied to wheat foliage, were the key driving forces in stabilizing the grain yield of wheat at a high level (Figure 2). The observed types of wheat responses to field protection systems resulted mainly from their interaction with a wide range of the applied N doses, which increased from 40 to 240 kg ha −1 . Due to the lack of fungicide protection, despite of the use of micronutrients, the average yield was 8.4 t ha −1 . This yield was achieved by applying only 80 kg N ha −1 , which can be considered very high in Poland [5,47]. The main reason for such high yield was the forecrop, i.e., winter oilseed rape. This crop creates good conditions for wheat as a succeeding crop [48]. A progressive increase in N doses resulted in a gradual increase in the yield, but up to only 9.1 t ha −1 . The yield growth rate on the plot fertilized with 80 to 240 kg N ha −1 was 10-folds slower compared to the rate obtained, applying less than 80 kg N ha −1 . The obtained quasi stagnation in the grain yield suggests the dysfunctionality of the sink strength of wheat, i.e., grain number per unit area (GD, grain density). This conclusion was directly confirmed by the lower number of grain per ear (GR) but not by the number of ears per m 2 (ED, Table 4). The results obtained confirm the common opinion that the stages immediately preceding flowering are of key importance for the yield of wheat grain [20,49].
Wheat grain yields showed a linear, and at the same time, a strong dependence (R 2 = 0.98) on the grain density (GD). It has been well-documented that this trait of wheat is characteristic of modern wheat varieties. The size of this component of the yield depends on the nitrogen supply to plants over a period from the stem elongation to heading [16,50,51]. The lowest yield recorded in 2016, despite the good amount and distribution of rainfall, was the result of the lower number of grains per ear (GR). This component of the yield was 28% lower than in other years. Such a large decrease in GR, due to its significant impact on both the GD and grain yield, necessitates an identification of the cause. At the same time, no compensation in TGW was found. Slafer et al., [52] defined this cause as "yield regulator". The main reasons for the GR reduction was the excess of N taken up by wheat in the vegetative phase of development. This explanation is supported by the value of Nitrogen Unit Accumulation (NUA), which was 1 / 3 higher compared to its values in the other years of study. An excess of N during the booting phase, which is critical for GR formation, leads to increased competition between the growing ear and the vegetative parts of wheat for the assimilates. Excessive and prolonged growth of the vegetative tissues of wheat, as a consequence, reduces the number of fertile flowers [53,54].
Globally, yield losses due to wheat diseases are in the range of 10-28% [55]. The yield gap due to pathogen pressure can be bridged by an increase in N doses, but only with the simultaneous use of fungicides. Kulig [56] found that under balanced conditions of the use of N, with an appropriate level of fungicide protection, the wheat grain yield increased by 35%. These two complementary facts provide a basis for an investigation of the most controversial problem in cereal production, i.e., the assessment of the yieldforming role of plant protection. In the Rabbinge [57] concept, plant protection belongs to the group of yield-reducing factors. This means that plant protection does not directly increase the yield by affecting the yield components. According to this view, the role of plant protection is to maintain the maximum yield produced by plants due to its optimal supply with water and N. Our study has provided evidence that contradicts this view. The conducted study showed that the yield of wheat treated with fungicides increased in line with N doses in comparison with the fungicide control, i.e., the object treated only with N and with micronutrients. The relative increase in GD compared to the fungicide control was approximately 18%. The GD increment level clearly emphasizes the yield-forming effect of fungicidal protection on yield components of wheat. This result can be explained by the effect of fungicides on N management in wheat during the formation of fertile flowers in the wheat ear. It has been well documented that the high dose of the applied N creates a favorable environment for the growth of biotrophs [35]. The positive effect of plant protection was indirectly confirmed by the yield stagnation on a plot fertilized with 80 kg N ha −1 in the fungicide control. At the same time, the amount of N accumulated by wheat increased linearly with the applied N doses (Figure 3). Pathogen pressure in combination with excess of N leads to a decrease in the green surface of the physiologically active parts of wheat [58,59]. In our case, this unexpected interaction resulted in a reduction in GD. The yield increase in wheat grain yield in response to fungicide application was 2.951 t ha −1 compared to the fungicide control. Thus, the hypothesis on the yield-forming effect of plant protection was fully confirmed. In light of the obtained results, the yieldforming effect of the applied micronutrients to wheat was found to be insignificant. This does not mean that the applied micronutrients did not fulfill their yield-forming function, taking into account GD and NUE. Their effect was detected, but wheat was protected with fungicides. The applied doses of N were exploited up to 227 kg ha −1 on the object with fungicides and theoretically up to 245 kg ha −1 on the object with both fungicides and micronutrients. The simultaneous combination of fungicide and micronutrients resulted in an additional increase in the yield by 0.39 t ha −1 .

Agronomic Measures and Nitrogen Use Efficiency (NUE)
Intensive wheat fertilization with N is a common strategy to achieve two goals, i.e., high yield and high N content in the grain [60,61]. Therefore, the wheat management strategy, based on the strong protection of plants against the pressure of pathogens, requires an in-depth assessment of the effectiveness of the applied rates of N fertilizers [62].
The effect of the tested winter wheat foliar protection systems on NUE showed the existence of two different regression models ( Figure 5). In the N + Mi system, i.e., fungicide control, the optimal total N accumulation (TN op ), amounting to 298 kg N ha −1 , resulted in a grain yield of 9.067 t ha −1 . In contrast, in both systems with fungicides, the grain yield responded linearly to the increase in total N accumulated by wheat. The differences between the wheat protection systems were assessed on the basis of the Nitrogen Unit Productivity (NUP). This index directly relates to Nitrogen Utilization Efficiency, describing the production value of N accumulated by a crop. The values of this index calculated for TN op of 298 kg N ha −1 reached 30.4, 34.1, and 34.5 kg grain per kg of NT, respectively for N + Mi, N + P, N + Mi + P. This set of figures clearly stresses the greater effectiveness of N accumulated in wheat on the fungicide-treated objects. This study confirms the trends observed in other studies [62,63]. However, these reports did not explain the reason behind higher NUE. In our study, the ineffectiveness of N in wheat in the fungicide control resulted from inadequate levels of GD due to a reduction in GR. This condition is synonymous with the underdevelopment of the physiological sink [54,55]. The reasons for the ineffectiveness of N management in the fungicide control can be explained on the basis of the Nitrogen Unit Accumulation index (NUA). The increase in NUA in response to gradually increasing N doses showed a negative correlation with NHI, NUP-G, NUP-T, as well as with GR. This set of relationship indicates an unsustainable redistribution of N taken up by wheat between the grain and its vegetative tissues. The high share of N in the vegetative wheat tissues was a consequence of the reduced grain density on the object not protected with fungicides. This type of response suggest that these processes were revealed in the early stages of the ear growth [64,65]. The grain yield was not compensated for by the increase in grain weight (TGW), as is often assumed [66,67]. This response explains the reason behind the relatively low accumulation of N in wheat grains during the grain-filling period [68]. However, on objects treated with fungicides, the maximum grain yield in relation to the N supply could be achieved with TN op of 370 and 424 kg N ha −1 , for the N + P and N + Mi + P systems, respectively. Therefore, it can be concluded that, in these two systems, the increase in the wheat grain yield was limited due to size of the source, i.e., the deficiency of NT in the wheat biomass at harvest.
The productivity of the applied N fertilizer doses was also assessed with the use of four classical NUE indices. Among the examined indices, PFP-N allows for a quick assessment of the effectiveness of the applied N [69]. In general, the course of the obtained regression model best fit the power function ( Figure 4). The effect of the wheat foliar protection systems was revealed by a slightly smoother slope of the curve in response to gradually increasing N doses. In the P + Mi object, the PFP-N value of 55 kg grain kg −1 N f was obtained on the plot fertilized with 160 kg N ha −1 . In the objects with fungicides, this value was recorded on the plot fertilized with 200 kg N ha −1 . The increase in the productivity of the N dose by 40 kg ha −1 indicates a longer activity of wheat. In the studied case, it was directly related to the noticeably higher GD, and thus the increased wheat demand for N ( Table 5).
The key index, summarizing the yielding-forming effect of the applied N is its recovery (R-N) by the currently grown crop. The R-N values showed a high variability between years ( Figure A2). A phenomenon known as the priming effect as observed for all N rates in 2015 an on the plot with 40 kg N ha −1 in 2014. The observed excessive N uptake from the soil can only be explained by the extremely high N demand of the growing ear wheat. Two factors support this explanation, i.e., high GD and stable TGW. The observed type of wheat response is in line opinion provided by Kong et al. [17]. These authors found that under certain growing conditions, i.e., with increased grain demand for N, it is also taken from the soil during the grain filling period.

Conclusions
A comparison of the three foliar protection systems for winter-wheat production clearly shows that the yield gap was reduced by the use of fungicides. It can be further reduced by the simultaneous use of fungicides and micronutrients. This protection system is specifically recommended in wheat production systems with intensive use of N fertilizer. The strong dependence of the increase in the grain yield on the number of grain-per-unitarea clearly confirms that fungicides in wheat should be treated as the yield-forming factor. The ineffectiveness of the applied N in the fungicide control was fully confirmed by the Nitrogen Use Efficiency indices. The Nitrogen Unit Accumulation index, which expresses the efficiency of N uptake was not correlated with the amount of N accumulated in wheat grain, but with the amount in harvest residues. These non-productive dependencies were the result of disturbance of the N partition between the wheat grain and its vegetative tissues on the object without fungicide protection. The key reason of yield reduction was the lower number of grain per unit area, which consequently reduced the N demand of the enlarging grains for N during the grain-filling period. The observed disturbance in the N-flow was confirmed by both the significantly lower Nitrogen Harvest Index and its negative response to increased NUA. The increase in the indices of the Partial Factor Productivity of fertilizer on objects treated with fungicides indirectly indicates a higher productivity of the applied N fertilizer. On the one hand, the obtained results confirm the dominance of the physiological sink over the source. On the other hand, they justify the basic goal of the study, i.e., the yield forming role of plant protection. The foliar application of microelements supports the action of the fungicidal protection of wheat, increasing the productivity of high N rates. Funding: This research received no external funding.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the private property of Witold Szczepaniak.