Impacts of N-P-K-Mg Fertilizer Combinations on Tree Parameters and Fungal Disease Incidences in Apple Cultivars with Varying Disease Susceptibility

Adequate mineral fertilization helps to ensure optimal tree growth, fruit development, and predictable yield of apple trees. This 7-year study (2016–2022) aims to investigate the effect of nitrogen (N), phosphorus (P), potassium (K), and magnesium (Mg) fertilizer combinations (NP, NPK, NPKMg, and control) on eight parameters (trunk cross-sectional area—TCSA; fruit yield—FY; number of fruit per tree—FNT; crop load—CL; fruit diameter—FD; fruit weight—FW; fruit scab incidence—FSI; and powdery mildew incidence on shoot—PMIS) on the cultivars (cvs) ‘Golden Reinders’ (disease susceptible) and ‘Pinova’ (scab and mildew tolerant). In the 7-year period, TCSA values continuously increased for both cultivars over the years. Fertilizer treatments showed significant differences on TCSA but the effect varied greatly annually among fertilizer treatments. Fertilizer treatments had increasing effects on FY and FNT in 2018 and 2022, on CL in 2018, on FD in 2018 and 2019, and on FW in 2016 and 2018 in both cultivars compared to the control treatment. FSI values were the lowest in the NPKMg treatment for cv. ‘Golden Reinders’ in 2016, 2017, and 2022; for cv. ‘Pinova’ in 2016; PMIS values for cv. ‘Golden Reinders’ in 2017, 2018, 2021, and 2022; and for cv. ‘Pinova’ in 2018. Correlation and regression analyses revealed strong and significant (p = 0.05) relationships between FNT versus (vs.) TCSA, FNT vs. FY, FW vs. TCSA, CL vs. FY, FW vs. FD, and FSI vs. FW. In conclusion, our study showed that multiyear application of fertilizer combinations can successfully increase TCSA and yield parameters as well as reduce fungal disease incidences, especially on the disease-susceptible cultivar in sandy soil with moderate fertility, under Central-European continental climate conditions.


Introduction
Apples are one of the largest fruit crops in temperate climates, with worldwide production showing a growing trend that exceeded 90 million tons in 2022 [1].A key component of apple production technology is the provision of adequate nutrients to the trees through mineral fertilization, which ensures optimal growth, development, and predictable productivity of apple trees [2][3][4].
Flowering and fruit set in apple trees are highly dependent on phosphorus and potassium, which are crucial for the initiation and development of flower buds.These nutrients play significant roles in the processes of flowering, pollination, and subsequent fruit set [19][20][21][22][23]. Potassium, in particular, is essential for fruit size and quality [24,25].The demand for potassium in apple trees peaks during ripening [24,25].Deficiencies in phosphorus and potassium can result in reduced flower formation, poor pollination, and low fruit set, ultimately leading to a decrease in overall fruit quantity [14,15,[19][20][21][22][23][24][25].
Adequate nutrient levels ensure optimal fruit development (size and weight) and reduce premature fruit drop [14,18,[26][27][28] as well as support the accumulation of carbohydrates and other compounds in the fruit [15,17,21].
Optimal uptake of NPK and Mg nutrients by apple trees increases tolerance to diseases and pests, while imbalances or deficiencies in nutrients can lead to reduced photosynthesis and increased susceptibility to diseases and pests [26,29,30].Certain nutrients, such as potassium and calcium, help strengthen cell walls and improve the tree's resistance to diseases.However, N excess often increases the susceptibility of plants to fungal diseases [31][32][33][34][35][36][37][38].
Increasing N fertilization was shown to enhance powdery mildew development on strawberries [36] and dogwoods [34].In the case of apples, an excess of soil-applied nitrogen was demonstrated to increase disease development caused by Neonectria ditissima [37,38].Nitrogen excess induces vigorous shoot growth in apple trees, which elevates the susceptibility of shoots to Podosphaera leucotricha [35,38] and reduces leaf resistance to Venturia inaequalis [31,32].
Although numerous studies have evaluated nutrient supply for apple vegetative and generative parameters, the effect of nutrient supply on apple tree vegetative and generative parameters, along with fungal diseases in tolerant/resistant and susceptible cultivars, has rarely been investigated in long-term multiyear studies.Additionally, correlations among tree parameters and disease incidences have not been explored, but such analysis may help us better understand the influence of nutrient supply on fungal diseases in relation to tree parameters.This 7-year study aims to investigate the effect of four fertilizer treatments including nitrogen, phosphorus, potassium, and magnesium fertilizer combinations (NP, NPK, NPKMg, and control) on six tree parameters (trunk cross-sectional area, fruit yield, fruit number, crop load, fruit size, fruit weight) and on the incidences of apple scab on fruit and powdery mildew on shoots.The experiment was conducted on the disease-susceptible cultivar (cv.) 'Golden Reinders' and on the scab and powdery mildew-tolerant cv.'Pinova'.

Study Area, Plant Material, Orchard Management, and Meteorological Assessment
A seven-year study from 2016 to 2022 was performed in an experimental apple orchard of the University of Debrecen, Pallag Experimental Station, Debrecen-Pallag, Hungary (47 • 35 ′ 31.5 ′′ N, 21 • 38 ′ 19.3 ′′ E).Soil properties of the experimental site are summarized in Table 1.The soil type is light and sandy with low humus content.The soil is slightly alkaline (pH 7. 5-7.6).The nitrogen of the soil was lower than the optimum, while phosphorus, potassium, and magnesium in the topsoil were optimal for the standard sandy soil type (Table 1), according to the Agricultural Technical Guidelines [39].
The experimental orchard was established in the spring of 2006.Trees were grafted on M.26 rootstock and trained to a slender spindle canopy (4.0 × 1.5 m = 1667 trees ha −1 ).Before the experiment, annual NPK fertilization was applied at rates of 30-40-50 kg ha −1 from 2006 onwards.At the commencement of the first year of experimentation in 2016, the trees were 10 years old with an average height of 3.5 m.Two apple cultivars were evaluated: the disease-susceptible cultivar 'Golden Reinders' [40], and the apple scab and powdery mildew-tolerant cultivar 'Pinova' [41].In the soil samples, the AL-soluble content of P, K, Ca, and Mg was determined according to Egnér et al. [42].The measurement of NO 3 − + NO 2 − − N content was performed according to Skalar [43].CaCO 3 , content was determined by the methods of Filep [44], using a Scheibler-type calcimeter.
Orchard management was carried out according to the European Integrated Fruit Production (IFP) guidelines [45].The plantation was equipped with a drip irrigation system.Winter pruning was performed once a year, in February.Mechanical fruit thinning was conducted annually in mid-June for cv.'Pinova', whereas it was only carried out in 2016, 2018, 2020, and 2022 for cv.'Golden Reinders', owing to reduced fruit set caused by spring frosts in 2017 and 2021.Concurrently, chemical thinning was applied in mid-May using 6-benzyladenine (Globaryll-100, Globachem Nv, Sint-Truiden, Belgium) during these years.
Meteorological parameters (mean temperature, minimum temperature, and precipitation) were recorded daily from January 2016 to December 2022 by the Metos Agrometeorological Station located in the experimental station.
For each fertilizer treatment and cultivar, seven trees were selected as assessment plots and replicated four times.Subsequently, the middle five trees were chosen for assessments within each plot.The first and seventh trees in the plots served as buffer trees without assessment.The application timing of the fertilizers is summarized in Table 3. Nitrogen fertilizer was always applied in spring due to the risk of nutrient leaching during winter.Fertilizers were applied as soil application at a distance of 1 m from the tree trunk on both sides of the trees on an annual basis.Fertilizers were incorporated into the soil at a shallow depth.The soil in the control treatment was also disturbed to the same shallow depth, without fertilizer application.

Assessment of Six Tree Parameters
Six tree parameters were assessed: trunk cross-sectional area (TCSA), fruit yield (FY), number of fruit per tree (FNT), crop load (CL), fruit diameter (FD), and fruit weight (FW).For each cultivar, 5-tree replicates per assessment plot were selected for assessments of all parameters.
Trunk thickness was measured after leaf fall in November each year with a Vernier caliper at the trunk halfway between the graft union and the main scaffold branches.From trunk thickness, TSCA was calculated in cm 2 .For TSCA, the difference in each year was also calculated, with 2016 as the reference point (0), and then Y2-Y1, Y3-Y2, Y4-Y3, Y5-Y4, Y6-Y5, and Y7-Y6.
Fruit picking was performed at the biological maturity stage of both cultivars, determined by the starch-iodine test.Harvest dates were 18, 15, 13, 4, 23, and 21 September in 2016, 2017, 2018, 2019, 2020, and 2022, respectively, for cv.'Golden Reinders', as well as 3 October 2016, 27 September 2017, 29 September 2018, 1 October 2019, 7 October 2020, 29 September 2021, and 28 September 2022 for cv.'Pinova'.Fruit yield (kg tree −1 ) and number of fruit per tree were measured by weighing/counting all the fruit on each selected tree.Crop load was calculated as fruit yield divided by trunk cross-sectional area (kg cm −2 ).Fruit diameter (mm) was assessed with a Vernier caliper based on 20 fruits per tree (100 fruits per treatment).Fruit weight (g) was measured using the same 100 fruit with a digital scale, with a precision of two decimal places.

Assessment of Incidences of Apple Scab and Apple Powdery Mildew
Apple scab [Venturia inaequalis (Cooke) G. Winter] assessments were conducted on fruit for each fertilizer treatment, cultivar, and year.For fruit assessment, 25 fruits with typical cultivar characteristics were observed on each quadrant of a selected tree at harvest in each year (4 × 25 fruits per tree and 5 × 4 × 25 fruits per assessment plot).A fruit was considered diseased if at least one visible scab lesion was observed on the fruit (Figure 1a).
Fruit scab incidences (FSI) from the four quadrants per tree were averaged to obtain the percentage of diseased fruits per tree.For each cultivar, the same trees were selected for fruit scab assessment as for the other tree parameters, in the years 2016, 2017, 2018, 2019, 2021, and 2022.1a).Fruit scab incidences (FSI) from the four quadrants per tree were averaged to obtain the percentage of diseased fruits per tree.For each cultivar, the same trees were selected for fruit scab assessment as for the other tree parameters, in the years 2016, 2017, 2018, 2019, 2021, and 2022.
Apple powdery mildew [Podosphaera leucotricha (Ellis & Everh.)E. S. Salmon] assessments were conducted at harvest on shoots of the same trees used for apple scab and other tree assessments in the same years.All shoots from each quadrant of a tree were examined for disease symptoms.Shoots were considered diseased if at least one leaf was covered with mycelium and/or spores (Figure 1b).Powdery mildew incidences on shoot (PMISs) from the four quadrants per tree were averaged to obtain the percentage of diseased shoots per tree.Data for each parameter were averaged to obtain a single value per tree for each cultivar, fertilizer treatment, and year.The data were subjected to analyses of variance (ANOVA) in order to determine the effect of cultivars, fertilizer treatments, and years, and their interactions were also assessed on the eight parameters.Subsequently, significant F-tests (p < 0.05) were followed by a least significant difference (LSD) test for each year, fertilizer treatment, and cultivar in order to compare the means of all parameters using LSD0.05 values.Genstat Release 9.1 (Lawes Agricultural Trust, IACR, Rothamsted, UK) was used for the analyses.Apple powdery mildew [Podosphaera leucotricha (Ellis & Everh.)E. S. Salmon] assessments were conducted at harvest on shoots of the same trees used for apple scab and other tree assessments in the same years.All shoots from each quadrant of a tree were examined for disease symptoms.Shoots were considered diseased if at least one leaf was covered with mycelium and/or spores (Figure 1b).Powdery mildew incidences on shoot (PMISs) from the four quadrants per tree were averaged to obtain the percentage of diseased shoots per tree.
2.5.Data Analyses 2.5.1.Analysis of Variance (ANOVA) Data for each parameter were averaged to obtain a single value per tree for each cultivar, fertilizer treatment, and year.The data were subjected to analyses of variance (ANOVA) in order to determine the effect of cultivars, fertilizer treatments, and years, and their interactions were also assessed on the eight parameters.Subsequently, significant F-tests (p < 0.05) were followed by a least significant difference (LSD) test for each year, fertilizer treatment, and cultivar in order to compare the means of all parameters using LSD 0.05 values.Genstat Release 9.1 (Lawes Agricultural Trust, IACR, Rothamsted, UK) was used for the analyses.

Correlation and Linear Regression Analyses
Relationships among all parameters were analyzed by Pearson correlation coefficients (r), and their associated significance levels (p = 0.05) were separately determined for the two cultivars.Subsequently, the strongest significantly correlated pairs were identified for the two cultivars and the overall data set.Additionally, the strongest significantly correlated pairs were plotted against each other, and linear regression functions were fitted for the four fertilizer treatments.A t-test was then used to determine whether the regression slopes were significantly (p = 0.05) different among the four fertilizer treatments.For the analyses, Genstat Release 9.1 (Lawes Agricultural Trust, IACR, Rothamsted, UK) was utilized.

Analysis of Variance
Analysis of variance on the parameters of trunk cross-sectional area, fruit yield, number of fruit per tree, crop load, fruit diameter, fruit weight, and incidences of apple scab and powdery mildew showed significant (p < 0.05) effects for fertilizer treatments, cultivars, and years (Table 4).The only exception was the parameters of trunk cross-sectional area, where the effects for cultivars were nonsignificant at the p < 0.05 level.Interactions for the treatment effects were nonsignificant for all treatments (Table 4).

Weather Conditions
The monthly mean temperature ranged from −6.0 to 25 • C during the assessed period from January 2016 to December 2022 (Table S1).The lowest minimum temperature values indicated that late spring frost in April (below −2 • C) occurred in 2020, 2021, and 2022.In 2021, cv.'Golden Reinders' suffered total frost damage; accordingly, it was not possible to evaluate all parameters in that year.The annual rainfall ranged from 397 to 720 mm between 2016 and 2022 (Table S1).

Trunk Cross-Sectional Area
In the 7-year period, trunk cross-sectional area (TCSA) values continuously increased for both cultivars with years (Table 5).For cv. 'Golden Reinders', TCSA was between 31.4 and 35.6 cm 2 in 2016, while after seven years, values ranged from 72.9 to 92.9 cm 2 in 2022.In cv.'Pinova', TCSA was between 32.0 and 48.0 cm 2 in 2016, while the values ranged from 72.0 to 102.5 cm 2 in 2022.Although 'overall (treatments)' TCSA values were higher for cv.'Pinova' compared to cv. 'Golden Reinders', the TCSA values of the two cultivars did not differ significantly at p = 0.05.indicates that there are no significant differences between the treatment values.
Significant differences among the four fertilizer treatments were found for cv.'Golden Reinders' in 2019, 2021, and 2022; and for cv.'Pinova' in 2016, 2017, 2018, 2020, 2021, 2022, and 'overall (year)' (Table 5).The values of the NPKMg treatment were significantly higher compared to the values of the NP treatment for cv.'Golden Reinders' in 2019, 2020, and 2022, while the two treatments were not significantly different from the other two treatments (control and NPK).In the case of cv.'Pinova', the values of the NP treatment were significantly higher compared to the values of the NPK treatment in all years (with the exception of 2019), while the two treatments were not significantly different from the other two treatments (control and NPKMg).
The overall values for the four fertilizer treatments, 'overall (cultivars)', showed that the overall values of the NPKMg treatment were significantly higher compared to the values of the NPK treatment in 2020, 2021, and 2022, while the two treatments were not significantly different from the other two treatments (control and NP) (Table 5).

Fruit Yield and Number of Fruit per Tree
The cultivar 'Golden Reinders' suffered from 100% frost damage of fruit in 2021, and partial fruit loss in 2017 and 2019, while cv.'Pinova' provided a more homogenous yield over the assessed period (Table 6).The highest fruit yield was recorded in the NPK treatment in 2022 on cv.'Golden Reinders' (69.1 kg tree −1 corresponding to 115.3 t ha −1 ).The 'overall (treatments)' fruit yield data showed significant variability in the case of cv.'Golden Reinders' (0-61.3kg tree −1 ) compared to cv. 'Pinova' (29.1-53.1 kg tree −1 ); the fruit yield of the two cultivars was significantly different in 2017 and 2019 at p = 0.05 (Table 6).Significant differences in fruit yield values among the four fertilizer treatments were found for cv.'Golden Reinders' in 2017, 2018, and 2022, and for cv.'Pinova' in 2016, 2018, 2019, and 2022 (Table 6).
The overall values of fruit yield for the four fertilizer treatments, 'overall (cultivars)', showed that the overall values of NPK and NP treatment were significantly higher compared to values of the control treatment in 2018 and 2022, respectively, while the two treatments were not significantly different from the other two treatments (NP and NPKMg, and NKP and NKPMg, respectively) (Table 6).
The cultivar 'Golden Reinders' had a wide range of number of fruit per tree (between 14 and 612), while cv.'Pinova' provided a more homogenous number of fruit per tree (between 143 and 387) over the assessed period (Table 7).The highest number of fruit per tree was recorded in the NPK treatment in 2020 on cv.'Golden Reinders' (612), with the lowest one in the NPK treatment in 2019 on cv.'Golden Reinders' (14).The 'overall (treatments)' data showed that the number of fruit per tree of the two cultivars was significantly different in 2017, 2019, 2020, and 2022 at p = 0.05 (Table 7).
Significant differences in the number of fruit per tree were found among the four fertilizer treatments for cv.'Golden Reinders' in all years, and for cv.'Pinova' in 2016, 2018, 2019, and 2022 (Table 7).For example, the values of the NP treatment were significantly different from the NPK treatment for cv.'Golden Reinders' in 2017, 2018, 2019, and 2020.In 2019, the values of the NP and NKPMg treatments for cv.'Golden Reinders' were significantly higher compared to the NKP treatment.In 2022, the values of the NPKMg treatment were significantly higher than the control.
In the case of cv.'Pinova', the values of NP, NKP, and NKPMg treatments were significantly higher compared to the values of the control treatment in 2016; the values of the NP treatment were significantly higher compared to the values of the control treatment Plants 2024, 13, 1217 9 of 21 in 2018 and 2022; and the values of the NPK treatment were significantly higher compared to the values of the control treatment in 2019 (Table 7).The overall values of the number of fruit per tree for the four fertilizer treatments, 'overall (cultivars)', showed that the overall values of the NPK treatment were significantly higher compared to values of the control treatment in 2016 and 2018, while the two treatments were not significantly different from the other two treatments (NP and NPKMg) (Table 7).In 2022, the overall values of the NPKMg treatment were significantly higher compared to the values of the control treatment, while the two treatments were not significantly different from the other two treatments (NP and NPK) (Table 7).

Crop Load
Similarly to other assessed parameters, cv.'Golden Reinders' had a wide range of crop load values (between 0.05 and 1.25 kg cm −2 ), while cv.'Pinova' provided a more homogenous crop load (between 0.42 and 1.00 kg cm −2 ) over the 7-year period (Table 8).The highest crop load was recorded in 2016 in the NPK treatment for both cultivars: 1.25 kg cm −2 for cv.'Golden Reinders' and 1.00 kg cm −2 for cv.'Pinova'.The lowest values were also found in the NPK treatment in 2019 on cv.'Golden Reinders' (0.05 kg cm −2 ).The 'overall (treatments)' data showed that the crop load of the two cultivars was significantly different only in two years (2019 and 2020) at p = 0.05 (Table 8).
Values of crop load among the four fertilizer treatments were significantly different for cv.'Golden Reinders' in 2018 and 2019, and for cv.'Pinova' in 2016, 2018, 2019, and 'overall (year)' (Table 8).The values of the NPK treatment were significantly higher than the NP and control treatments for cv.'Golden Reinders' in 2018.However, in 2019, the values of the NPK treatment for cv.'Golden Reinders' were significantly lower than the values of the NP the treatment, while the two treatments were not significantly different from the other two treatments (control and NPKMg).In the case of cv.'Pinova', the values of NKP treatment were significantly higher compared to values of the control treatment in 2016, 2018, 2019, and 'overall (year)', while the two treatments were not significantly different from the other two treatments (NP and NPKMg) (Table 8).
The overall values of crop load for the four fertilizer treatments 'overall (cultivars)' showed significant differences only in 2018.In this year, the values of the NPK treatment were significantly higher compared to the values of the control treatment, while the two treatments were not significantly different from the other two treatments (NP and NPKMg) (Table 8).

Fruit Diameter and Fruit Weight
During the 7-year period, the mean fruit diameter of both cultivars reached the market requirements of 70 mm, with the exception of cv.'Pinova' in 2021 in the control and NPK treatments (Table 9).The largest fruit (83.0 mm) was harvested in 2016 in the NPKMg treatment for cv.'Golden Reinders', while the smallest one (69.1 mm) in 2021, in the control treatment for cv.'Pinova' (Table 9).
Significant differences in the fruit diameter were found among the four fertilizer treatments for cv.'Golden Reinders' in 2016-2019, and for cv.'Pinova' in 2018, 2019, 2021, and 2022 (Table 9).In the case of cv.'Golden Reinders', the values of the NPKMg treatment were significantly higher than those of the NP treatment in 2016.In 2017, the values of the NPK treatment were significantly higher compared to the NP treatment.In 2018, the values of the NP and NKPMg treatments were significantly higher compared to the control treatment.In 2019, the values of the NPKMg treatment were significantly higher than the other three treatments (Table 9).In the case of cv.'Pinova', fruit diameter values of the NKPMg treatment were significantly higher compared to the values of the control treatment in 2018; the values of NP and NPKMg treatments were significantly higher compared to the values of the control treatment in 2019; the values of the NP treatment were significantly higher compared to the values of the control treatment in 2020; and the values of NPK and NPKMg treatments significantly differed from the control treatment (Table 9).
The overall values for the four fertilizer treatments 'overall (cultivars)' showed that the overall fruit diameter values of NPKMg treatment were significantly higher compared to the values of the control treatment in 2016, 2018, and 2019, while the two treatments were not significantly different from the other two treatments (NP and NPK) (Table 9).In 2021, the overall values of the NP treatment were significantly higher compared to the values of the control treatment, while the two treatments were not significantly different from the other two treatments (NPK and NPKMg).
Similarly to fruit diameter values, the fruit weight was the lowest in 2021 in the control treatment (133 g) for cv.'Pinova', while the highest incidence was observed in 2016 in the NPKMg treatment (238 g) for cv.'Golden Reinders' (Table 10).The 'overall treatment' data showed that the fruit weight of the two cultivars was significantly different in 2016, 2017, 2018, and 2020 at p = 0.05 (Table 10).
Significant differences in fruit weight were found among the four fertilizer treatments for cv.'Golden Reinders' in 2016-2019, and for cv.'Pinova' in 2016-2022 (Table 10).The values of the NPKMg treatment were significantly higher compared to the other three treatments in 2016.In 2017, the values of the NP treatment were significantly lower com-pared to the other three treatments.In 2018, the values of the NP and NKPMg treatments were significantly higher compared to the control and NPK treatments.In 2019, the values of the NPK and control treatments were significantly different from the NP and NPKMg treatments (Table 10).In the case of cv.'Pinova', fruit weight values of the NKPMg treatments were significantly higher compared to values of the control treatment in 2016, 2018, 2019, 2020, and 2021.In 2017, the values of the NPKMg treatment were significantly higher compared to the values of the NP treatment.In 2022, the values of the NPK and NKPMg treatments were significantly different from the control treatment (Table 10).
The overall values for the four fertilizer treatments 'overall (cultivars)' showed that the overall values of NPKMg treatment were significantly higher compared to values of the control treatment in 2016, 2018, and 2021 (Table 10).In 2017, the overall values of the NP treatment were significantly lower compared to the values of the NPK and NPKMg treatments.In 2019, the overall values of the NPK treatment were significantly lower compared to the values of the NP and NPKMg treatments.

Incidences of Apple Scab and Apple Powdery Mildew
Apple scab incidence on fruit was the smallest in 2018, 2019, and 2021 on cv.'Pinova' in all four fertilizer treatments (0%), while the largest one was observed in 2016 on cv.'Golden Reinders' in the NP treatment (14.0%) (Table 11).The 'overall (treatments)' data showed that the fruit scab incidence of cv.'Golden Reinders' was significantly higher compared to the values of cv.'Pinova' in 2016-2019 and 2022 at p = 0.05 (Table 11).
Significant differences in fruit scab incidence were found among the four fertilizer treatments for cv.'Golden Reinders' in 2016, 2017, and 2022, and for cv.'Pinova' in 2016 (Table 11).The values of the NPKMg treatment were significantly lower compared to the control treatment in 2016 for both cultivars.In 2017, the values of the NPKMg treatment were significantly lower compared to the control and NP treatments for cv.'Golden Reinders'.In 2022, the values of the NPKMg treatment were significantly lower compared to the control treatment for cv.'Golden Reinders' (Table 11).The overall values of fruit scab incidence 'overall (cultivars)' showed that the value of the NPKMg treatment was significantly lower compared to the values of the control and NP treatments in 2016 and to the value of the control treatment in 2017 (Table 11).
Powdery mildew incidence on shoots was the lowest (1.3%) in 2016 on cv.'Pinova' in the NKPMg treatment, while the highest incidence (7.3%) was observed in 2018 on cv.'Golden Reinders' in the control treatment (Table 12).The 'overall (treatments)' data showed that the powdery mildew incidence of cv.Golden Reinders' was significantly higher compared to the values of cv.'Pinova' in 2016-2019, 2021, and 2022 at p = 0.05 (Table 12).
Significant differences in powdery mildew incidence were found among the four fertilizer treatments for cv.'Golden Reinders' in 2017, 2018, 2021, and 2022, and for cv.'Pinova' in 2018 (Table 12).The disease incidence of the NPKMg treatment was significantly lower compared to the control treatment in 2017, 2018, 2021, and 2022 for cv.'Golden Reinders' and for cv.'Pinova' in 2018 (Table 12).
The values of powdery mildew incidence in 'overall (cultivars)' for the four fertilizer treatments showed that the overall value of the NPKMg treatment was significantly lower compared to the value of the control treatment in 2018 (Table 12).
a Differences among treatments are represented by LSD 0.05 values at p = 0.05.b Values coupled with different letters are significantly different at p = 0.05 according to LSD t-tests.If there are no letters beside the values, it indicates that there are no significant differences between the treatment values.c '-' missing data in 2020.

Correlation among Parameters and Linear Regression
Pearson correlation coefficient (r) was the highest (r = 0.95) between FNT and FY in the overall data analyses including both cultivars (Table 13).Six correlation pairs were significant (p = 0.05) in the overall data set and for each cultivar; five parameter pairs correlated positively (FNT vs. TCSA, FNT vs. FY, CL vs. FY, FW vs. FD, and FSI vs. FW), and one negatively (FW vs. TCSA) (Table 13).In the case of cv.'Pinova', three additional correlation pairs (FY vs. TSCA, FD vs. FNT, and FW vs. FNT) also showed strong and significant (p = 0.05) relationships (Table 13).In the case of cv.'Golden Reinders', two additional correlation pairs (CL vs. FNT and FD vs. TCSA) also showed strong and significant (p = 0.05) relationships (Table 13).
The linear regression analysis showed significant relationships for all six pair variables with r = 0.714-0.922and p = 0.045-0.001for the four fertilizer treatments.However, no differences were observed among the slope parameters for the six variable pairs among the control, NP, NPK, and NPKMg treatments as t-tests showed p-values ranging from 0.799 to 0.118.Slope parameters showed no significant differences between the two cultivars.FNT vs. TCSA relationships showed that most of the increasing values for the number of fruit per tree corresponded with increases in TCSA values (Figure 2A).In the case of TCSA vs. FW relationships, increasing TCSA values resulted in slight decreases in fruit weight for all four fertilizer treatments (Figure 2B).In the case of FNT vs. FY and CL vs. FY relationships, the majority of the number of fruit per tree values and crop load values were directly proportional to the increase in fruit yield for all four fertilizer treatments (Figure 2C,D).In the case of FD vs. FW relationships, increasing fruit diameter values resulted in sharp increases in fruit weight for all four fertilizer treatments (Figure 2E).In the case of FSI vs. FW relationships, the increase in the fruit scab incidence was not directly proportional to the increase in fruit weight (Figure 2F).Here, two clusters were separated.The cluster with values of zero or close to zero scab incidences includes the fruit of cv.'Pinova'.The linear regression analysis showed significant relationships for all six pair variables with r = 0.714-0.922and p = 0.045-0.001for the four fertilizer treatments.However, no differences were observed among the slope parameters for the six variable pairs among the control, NP, NPK, and NPKMg treatments as t-tests showed p-values ranging from 0.799 to 0.118.Slope parameters showed no significant differences between the two cultivars.FNT vs. TCSA relationships showed that most of the increasing values for the number of fruit per tree corresponded with increases in TCSA values (Figure 2A).In the case of TCSA vs. FW relationships, increasing TCSA values resulted in slight decreases in fruit weight for all four fertilizer treatments (Figure 2B).In the case of FNT vs. FY and CL vs. FY relationships, the majority of the number of fruit per tree values and crop load values were directly proportional to the increase in fruit yield for all four fertilizer treatments (Figure 2CD).In the case of FD vs. FW relationships, increasing fruit diameter values resulted in sharp increases in fruit weight for all four fertilizer treatments (Figure 2E).In the case of FSI vs. FW relationships, the increase in the fruit scab incidence was not directly proportional to the increase in fruit weight (Figure 2F).Here, two clusters were separated.The cluster with values of zero or close to zero scab incidences includes the fruit of cv.'Pinova'.2.

Discussion
In this 7-year study, we evaluated the effects of four fertilizer treatments (control, NP, NPK, NPKMg) on vegetative, generative, and disease incidence parameters of two apple cultivars ('Golden Reinders' and 'Pinova').The effects of fertilizer treatments on the observed parameters were dependent on year and cultivars.
The fruit yield of the two cultivars was highly dependent on the differences in the yearly abiotic environmental factors, especially for cv.'Golden Reinders'.Severe fluctuations in yield from year to year are reported for several apple cultivars [47,48], and alternate bearing is a common feature of cv.'Golden Reinders' in many growing areas [49,50].Previous studies reported that alternate fruit bearing was not observed for cv.'Pinova' [41], which was also confirmed by this 7-year study, as cv.'Pinova' produced a more  2.

Discussion
In this 7-year study, we evaluated the effects of four fertilizer treatments (control, NP, NPK, NPKMg) on vegetative, generative, and disease incidence parameters of two apple cultivars ('Golden Reinders' and 'Pinova').The effects of fertilizer treatments on the observed parameters were dependent on year and cultivars.
The fruit yield of the two cultivars was highly dependent on the differences in the yearly abiotic environmental factors, especially for cv.'Golden Reinders'.Severe fluctuations in yield from year to year are reported for several apple cultivars [47,48], and alternate bearing is a common feature of cv.'Golden Reinders' in many growing areas [49,50].Previous studies reported that alternate fruit bearing was not observed for cv.'Pinova' [41], which was also confirmed by this 7-year study, as cv.'Pinova' produced a more homoge-nous annual yield than cv.'Golden Reinders' (Tables 6-8).However, our findings clearly demonstrate a biennial pattern of fruit production for cv.'Golden Reinders', with reduced yields observed in 2017, 2019, and 2021, and improved yields in 2016, 2018, 2020, and 2022.This alternation coincides with varying levels of frost damage observed in our study.Our results align with previous research by Monselise and Goldschmidt [51], Atay et al. [51], and Netsawang et al. [48], which suggests that alternate bearing is influenced by multiple factors, including spring frost events for deciduous trees and drought stress during fruit set.
In previous studies on apples, fertilization has been reported to increase the TSCA values of the trees [52] or have no effect on TSCA values [53,54].In our study, TCSA values increased annually in both cultivars (Table 5).However, the fertilization effect was predominantly observed in the latter half of the 7-year experimental period, particularly for cv.'Pinova', which exhibits more balanced fruit-bearing characteristics.Pole et al. [52] demonstrated that nitrogen fertilization increased TSCA in apples, while our long-term results revealed a significant increase in TCSA with the combined NPKMg treatment in some years compared to NP or NPK fertilizer treatments.Conversely, previous studies have shown that fertilization did not affect the TCSA values of cultivars 'Sampion' and 'Golden Delicious Reinders' [53,54].Differences in the effects of fertilization on TSCA are likely attributed to variations in plant density, tree age, ecological conditions, and fertilization dosages in our study.Our fertilization study revealed significant correlations between TSCA and the number of fruit per tree, and between TSCA and fruit weight, which were not explored in previous studies (Table 13 and Figure 2A,C).
The response of fruit yield to nutrient supply has often been inconsistent over the years or has yielded controversial results [52][53][54][55][56].For instance, experiments involving nitrogen fertilization of apple orchards on fertile soil have frequently shown a lack of effectiveness of the fertilizer treatments [52,55,56], a trend that was also observed in this study in certain years when fertilization had no significant effect for either cultivar (Table 6).In our study, we observed optimum values of the nutrients P, K, and Mg (Table 1) in our sandy soil characterized by moderate fertility.This observation may elucidate the absence of detectable fertilization effects in certain years.Conversely, it has been noted that the optimal levels of P, K, and Mg in sandy soil are usually lower than those in 'standard soils' with higher fertility [39].This finding potentially explains why the effects of fertilization became evident in other years of our study (Table 6).These latter results were consistent with findings by Stefanelli et al. [57], who reported positive effects of fertilization within moderate doses.In another trial, the response in yield was observed only after 10 years of fertilization when nitrogen, phosphorus, and potassium fertilizers were applied singly and in combination [58].The authors demonstrated that while a single application of nitrogen reduced yield, combined NPK fertilizers significantly increased fruit yield [58].These results align with our findings on TSCA values (Table 5) and fruit yield values in 2018 and 2022, where the highest yield response was observed in the NPK treatments for cv.'Golden Reinders' (Table 6).Although our study did not involve the single application of nitrogen fertilizers, previous studies have indicated that the application of various doses of nitrogen had either a low or no effect on apple yield [59,60], but it improved shoot extension growth in cv.'Golden Reinders' apple trees [59].In addition, to our knowledge, previous studies have not demonstrated correlations among fruit yield parameters under fertilization treatments.Overall, our results on fruit yield indicate that the effects of combined NP, K, and Mg nutrients on fruit yield vary depending on several factors, including seasonal features, cultivar characteristics, and nutrient type combinations.This suggests that a reliable yield response to fertilization can be determined after long-term (min.10 years) experimentation and under specific orchard ecological conditions.
In many cases, external features such as fruit size and fruit weight hold greater importance to consumers than the internal traits [61].The cultivar 'Pinova' typically exhibits an average fruit weight of 175 g and a fruit diameter of 75 mm [41,62].In our study, under various fertilization treatments, cv.'Pinova' produced fruits of 133-203 g, with diameters ranging from 69.1 to 77.9 mm.These values align with findings from a study utilizing nanotechnology-based foliar fertilizers [63].For cv. 'Golden Reinders', the fruit diameter consistently reached 70 mm each year with fruit weights ranging from 145 to 238 g, consistent with results reported in other studies [20,64].Bielicki and Pasko [65] reported a yield of 42.2 kg tree −1 for cv.'Golden Reinders' on M.26 rootstock.Notably, our fertilization trial yielded higher values for cv.'Golden Reinders', such as in 2022, when fruit yield reached 64.0-69.1 kg tree −1 (Table 6).Our study indicates that fertilization treatments were effective in improving fruit weight and size in both cultivars, as evidenced by generally smaller fruit diameters and weights observed on control trees in most years.This was also revealed with a positive correlation between fruit weight and fruit diameter (Table 13, Figure 2E).Our findings are consistent with the results of Zhao et al. [66] and Zijian et al. [67], who demonstrated that fertilization increased fruit quality, including mean fruit weight and diameter.Although we employed identical fertilization regimes for both cultivars in our study, owing to their comparable growth characteristics, it is conceivable that cultivar-specific fertilization strategies could offer greater benefits to growers, as evidenced by prior research [68,69].
Disease-tolerant/resistant and -susceptible cultivars exhibited differing levels of susceptibility in terms of FSI and PMIS in this study (Tables 11 and 12), consistent with findings from previous studies [70][71][72].Our fertilizer treatments clearly demonstrated that complete nutrient supply (NPKMg) consistently reduced disease incidence for both diseases across most cases and years, especially on the disease-susceptible cv.'Golden Reinders' (Tables 11 and 12).This observation is consistent with prior research, indicating that optimal uptake of NPK and Mg nutrients by apple trees enhances tolerance to diseases, while nutrient imbalances may increase susceptibility to diseases and pests [26,29,30].In this study, we adhered to a standard fungicide spray regimen in accordance with integrated fruit production guidelines.It is plausible that the disparities would have been more pronounced in the absence or reduction of efficient fungicidal sprays.We also need to note that some nutrients, such as N, can increase apple tree canker development caused by N. ditissima [37,38].In addition, N excess induces vigorous shoot growth in apple trees, which elevates the susceptibility of shoots to P. leucotricha [35,38] and reduces leaf resistance to V. inaequalis [31,32].However, for instance, N fertilization in autumn can enhance the decomposition of leaf litter, thereby reducing the overwintering inoculum of V. inaequalis and consequently decreasing scab incidence in the following spring [35,[73][74][75][76][77].
Our study also revealed that the effect of fertilization on disease incidence was notably lower in the disease-tolerant cultivar 'Pinova' (Tables 11 and 12).This can be attributed to the inherently lower disease incidences of the disease-resistant cultivar in the control treatment, a phenomenon also observed in previous studies investigating cultivar susceptibilities to diseases [40,41,70,78].However, while fertilization treatments for FSI and PMIS showed little significant correlation with other measured parameters and between each other (Table 13), we found that only fruit weight exhibited a significant positive correlation with fruit scab incidence (Figure 2F).This finding is consistent with practical experience, as scabbed fruits tend to be smaller in size and weight compared to healthy fruits [78,79].

Conclusions
Our study demonstrated that the multiyear application of fertilizer combinations could effectively increase TCSA and yield parameters in two apple cultivars while also reducing fungal disease incidences in the disease-susceptible apple cultivar grown in sandy soil with moderate fertility, under Central-European continental climate conditions.However, the effect could be largely influenced by environmental and meteorological conditions of the year, as well as by the alternate bearing features of the cultivar.Furthermore, the outcome of unbalanced nutrient addition often did not differ significantly from that of balanced fertilizer combinations.Our study suggests a recommendation for orchardspecific practices, which necessitates thorough local investigations for several years.This duration allows for a comprehensive understanding of each tree's specific requirements within the given environmental circumstances.

Figure 1 .
Figure 1.Representative symptoms of apple scab on fruit (a) and representative symptoms of powdery mildew on shoot (b).Photos by I.J.Holb.

Figure 1 .
Figure 1.Representative symptoms of apple scab on fruit (a) and representative symptoms of powdery mildew on shoot (b).Photos by I.J.Holb.

Table 1 .
Seven soil parameters in the experimental site (Debrecen-Pallag, Hungary, 2016).Optimal values are given for sandy soil according to the Agricultural Technical Guidelines [39].

Table 2 .
Four fertilizer treatments applied and doses of fertilizer active ingredients (nitrogen-N; phosphorus-P 2 O 5 ; potassium-K 2 O; and magnesium-MgO) in the experimental apple orchard at Debrecen-Pallag, Hungary, from 2016 to 2022.

Table 3 .
Application time (dd mm yyyy) of the fertilizers (nitrogen-N; phosphorus-P; potassium-K; and magnesium-Mg) in the experimental apple orchard at Debrecen-Pallag, Hungary, from 2016 to 2022.
a Differences among treatments are represented by LSD 0.05 values at p = 0.05.b Values coupled with different letters are significantly different at p = 0.05 according to LSD t-tests.If there are no letters beside the values, it
a Differences among treatments are represented by LSD 0.05 values at p = 0.05.b Values coupled with different letters are significantly different at p = 0.05 according to LSD t-tests.If there are no letters beside the values, it indicates that there are no significant differences between the treatment values.c '-' No data available for cv.'Golden Reinders' in 2021 due to frost damage.
a Differences among treatments are represented by LSD 0.05 values at p = 0.05.b Values coupled with different letters are significantly different at p = 0.05 according to LSD t-tests.If there are no letters beside the values, it indicates that there are no significant differences between the treatment values.c '-' No data available for cv.'Golden Reinders' in 2021 due to frost damage.
a Differences among treatments are represented by LSD 0.05 values at p = 0.05.b Values coupled with different letters are significantly different at p = 0.05 according to LSD t-tests.If there are no letters beside the values, it indicates that there are no significant differences between the treatment values.c '-' No data available for cv.'Golden Reinders' in 2021 due to frost damage.
a Differences among treatments are represented by LSD 0.05 values at p = 0.05.b Values coupled with different letters are significantly different at p = 0.05 according to LSD t-tests.If there are no letters beside the values, it indicates that there are no significant differences between the treatment values.c '-' No data available for cv.'Golden Reinders' in 2021 due to frost damage.
a Differences among treatments are represented by LSD 0.05 values at p = 0.05.b Values coupled with different letters are significantly different at p = 0.05 according to LSD t-tests.If there are no letters beside the values, it indicates that there are no significant differences between the treatment values.c '-' No data available for cv.'Golden Reinders' in 2021 due to frost damage.
a Differences among treatments are represented by LSD 0.05 values at p = 0.05.b Values coupled with different letters are significantly different at p = 0.05 according to LSD t-tests.If there are no letters beside the values, it indicates that there are no significant differences between the treatment values.c '-' missing data in 2020, and no data available for cv.'Golden Reinders' in 2021 due to frost damage.