1. Introduction
Vegetables from the
Allium genus, such as onion (
A. cepa), garlic (
A. sativum), and leek (
A. ampeloprasum var.
porrum), are among the main sources of bioactive compounds in the human diet all over the world. Onion (
Allium cepa L.) is the most widely cultivated species of the genus
Allium. Both onion parts—leaves, but especially the bulb—are rich sources of bioactive compounds, mainly polyphenols, among which flavonoids dominate, characterized by elevated antioxidant and antibacterial activities. The qualitative health attributes of onion are widely exploited in human protection against cancer, high blood pressure, and infection [
1,
2].
On a worldwide scale, onions are among the most commonly consumed vegetables after carrots, cabbages, tomatoes, and cucumbers. Onion has gained in popularity during the last decades of the 20th century. Onion consumption has reached high levels in the United States and Europe [
3]. The largest onion producers are China, India, the USA, Pakistan, and Turkey. In Europe, the leading producers are Russia, the Netherlands, Spain, Poland, and Germany [
4,
5,
6]. The highest yields, at levels of 48–52 t ha
−1, are reached in the United States and the Netherlands, while the world average is much lower at around 19 t ha
−1. In Poland, onion yields are low, ranging in the last decade around 20–25 t ha
−1 [
4]. In spite of this, onion currently accounts for almost 40% of domestic exports of unprocessed vegetables [
5].
Onion is a vegetable crop with a high requirement for nutrients. The unit nutrient accumulation of a moderately yielding onion is as follows: N: 4.2; P: 0.8; K: 3.9; and Ca: 2.73 kg t
−1 of fresh weight (FW) of total onion biomass at harvest [
7]. These values are much higher compared to data for high yielding onion: N: 2.2; P: 0.5; K: 2.2; Ca: 0.6; Mg: 0.2; and S: 0.6 kg 1.0 t
−1 FW [
8]. Applying fertilizers is therefore an important agronomic practice to increase the yields of vegetables including onion, but its response to applied amounts mostly depends on the fertility status of the soil. Continuous work on improving soil fertility conditions for onions creates a basis for the effective use of fertilizer N [
9,
10]. Applied N doses fluctuate from about 50 to above 200 kg ha
−1, depending on the region, type of growing onion, and other production practices, like irrigation [
11,
12]. Too high an N dose or late application leads to prolonged growth of leaves and delayed bulb maturity. In addition, prolonged leaf growth delays or prevents bulbing, consequently resulting in a yield decline [
13].
The role of nutrients in the efficiency of N fertilizers is frequently limited to phosphorus (P) and potassium (K) [
12,
14]. However, there is growing interest among farmers in the application of S fertilizers to onion. There are two key reasons for S use in onion cultivation: yield and, even more important, the expected impact of S fertilizers on bulb pungency and the production of bioactive compounds [
15,
16]. The uptake of S by a high yielding onion can reach as high as 100 kg ha
−1. Therefore, onion is frequently referred to as a sulfur-loving crop. It is for these reasons that S is the essential fertilizer for proper development, high bulb yield, and high levels of bioactive compounds.
In spite of the above presented facts, onion response to S fertilizer is a matter of controversy [
17], based on data from pot experiments showing that increasing N rates significantly affected biomass, N and S uptake, and pungency of onion cultivars [
17,
18]. Some onion covarieties showed a stronger response to N than to S supply. However, extremely high rates of N can reduce S uptake by onion plants and decrease bulb pungency. In most cases field studies have shown a positive impact of applied S on onion yield. The optimum range of applied S ranges from 30 to 60 kg ha
−1 as recorded in different countries [
18,
19,
20]. In contrast, no impact of S application to field-grown onion on either yield or bulb pungency was found in [
21]. One of the key reasons for the very broad S ranges is the type of applied fertilizer.
The dominant S carrier is gypsum [
18,
22,
23], and to a much lower degree elemental sulfur (S
0) [
24,
25]. A comparative study on the impact of gypsum and S
0 on onion yield revealed an advantage of gypsum [
20,
26]. An analysis of the interactional impact of N and S on onion growth, yield, and bioactive compounds highlighted the necessity to optimize the doses of both nutrients [
22].
In spite of numerous experiments devoted to the response of onion to N or S application and their interactive impact on yield, the in-season trend of dry matter, N, and S accumulation is not known. This is the main weakness of any study on the yield-forming function of these two nutrients. The key objective of the present study was to define in-season models of dry matter, N, and S accumulation by onion and the parameters of the functions obtained. The next objective of this study was to define a set of parameters responsible for onion yield.
3. Results
3.1. Yield of Onion Bulbs
The key factor affecting the yield of onion bulbs is the course of the weather in a particular growing season, i.e., year (
F2.54 = 297.8,
p < 0.001). The highest yield was recorded in 2009, followed by 2011, and the lowest in 2010 (
Figure 2a). Averaged over experimental years, harvested yield was at the level of 2338 g m
−2 of bulb FW. The effect of the S rate on onion yield was independent of year and N rate. Each year, a linear yield increase up to an S rate of 60 kg ha
−1 was recorded. As a result, S application significantly increased the mean yield compared to the control (
F2,54 = 10.1,
p < 0.001). The relative yield increase averaged over N rates and years was 8.9% and 12.9% for an S rate of 30 and 60 kg ha
−1 (
Figure 2a).
In the case of N, its effect was variable in consecutive years of the study (
F4.54 = 9.9,
p < 0.001). In 2009, the yield significantly responded to N rates, reaching the highest value with fertilizer treatment of 60 kg ha
−1 (+32.8% compared to the N control). Plants fertilized with an N rate of 120 kg ha
−1 had much lower yield (+19.6% with respect to the N control). In 2010 and 2011, the effect of N doses was not significant, but a positive trend was observed on plots fertilized with the highest rate of N (
Figure 2b).
3.2. Dry Matter Accumulation Dynamics
Weather was the key factor, significantly differentiating the total dry weight (TDW) of onion plants (leaves + bulb) in consecutive samplings during the season (
Figure 3). At harvest, the highest TDW of onion was recorded in 2009, and the lowest in 2010. The difference between these two years was double, reaching 122%, but between 2010 and 2011 was only 50.2%. Following the 50th DAE, the onion plant began a phase of bulbing and a partition of total onion biomass between leaves (chive) and bulb was recorded. In 2009, the bulb biomass doubled in the period from DAE 70 to 80. A significant but smaller increase was recorded following the 80th DAE, although the biomass of leaves increased at a much slower rate, but progressively. A slightly different pattern was revealed in 2010. The progressive bulb increase lasted up to the 80th DAE, and then stopped. At the same time, leaf biomass did not show any increase. A completely different pattern of onion growth was observed in 2011. The bulb biomass increased step-by-step up to harvest, but growth of leaves stopped on the 70th DAE.
As shown in
Table 2, the bulb dry weight index (BDWI) showed a linear trend, reaching values ranging from 72% in 2010 to 77% at harvest in 2009. At the onset of bulbing (50 DAE), the lowest BDWI index of 23% was recorded in 2011, whereas in other years it was 2.5-fold higher. The key reason for the low value in this particular year was a high content of soil mineral N at sowing (
Table 1). This pattern of onion growth fully corroborates the opinion in [
13] that a high N supply during the early stages of development extends the phase of leaf growth. Nevertheless, the rate of bulb growth in later stages was sufficiently high to reach a BDWI of 74% at harvest.
The in-season trend of TDW compared to onion yield was, in general, significantly affected by the weather. On consecutive sampling dates, however, onion biomass was differentiated by the interaction of applied doses of N and S. The degree of impact of a particular factor on onion TDW was DAE-specific. Experimental factors did not affect TDW on the 20th and 60th DAE (
Table 3).
In general, beginning on the 40th DAE, the lowest onion TDW was recorded each year on the absolute control plot (S0 + N0 treatment, acronym AC).
On the 30th DAE, the highest biomass was recorded at the S30 + N0 plot. In the next 10 days, the highest biomass was recorded at two plots, S30 + N60 and S60 + N120. On the 70th DAE, the highest TDW was again found at the S30 + N60 plot, and 20 days later at the S60 + N120 plot. A significant difference with respect to the absolute control (AC) at DAE 90 was also recorded on the S0 + N120 and S60 + N60 treatments.
The significant impact of fertilization on BDWI values was recorded on two sampling dates, on the 50th and 90th DAE (
Table 3). The highest BDWI value of 49% at 50 DAE was found at the AC plot. Significantly lower indices were recorded on plots fertilized with S. This phenomenon can be explained by the higher supply of N to the growing plants, subsequently resulting in delayed onion bulbing [
13].
The parameters of the expolinear model characterizing the dynamics of onion TDW during the growing season are shown in
Table 4. In 2009, plants showed a higher growth rate, as expressed by C
m, during the linear phase of growth. The value of this parameter was ninefold and fivefold higher as recorded in 2010 and 2011, respectively. The main reason for these differences was the delay in the day of growth pattern change from exponential to linear (T
b).
In 2009, this particular day occurred on the 86th DAE, whereas in 2010 and 2011, it was on the 29th and 39th DAE, respectively. However, the predicted onion biomass (TDW
Tb) at this T
b was several-fold higher in 2009 compared to other years, clearly corroborating the close relationship with C
m (r = 0.99 ***). The relative rate of DW growth during the linear phase (R
L) displayed the same trend as the above-described parameters. A reverse trend was shown by R
m, a parameter describing the relative rate of onion growth during the exponential phase of development. This parameter was negatively correlated with R
L, T
b, TDW
Tb, and finally onion yield (
Table S1).
The interactional effect of S0 and S60 and increasing N rates significantly and positively affected four of five parameters of the developed model: Cm, RL, Tb, and TDWTb. The opposite relationship was recorded for Rm. On the S60 mai plot, the highest Cm and RL values were recorded for plants fertilized with an N rate of 60 kg ha−1, consequently resulting in t he highest TDWTb. Plants fertilized with an N rate of 120 kg ha−1, with the exception of Rm, had much lower values of all parameters of the model. Among the studied treatments, the highest Cm was recorded for plants grown on the S60 + N120 plot. Its value ranged from 3.5 to 9-fold higher compared to the values obtained on the S0 + N0 and S0 + N120 plots, respectively. It is necessary to stress that the predicted TDWTb, in spite of being the highest, had only a theoretical value. It was obtained on the 115th DAE, whereas the onions had been harvested 25 days earlier (DAE 90).
3.3. Nitrogen Uptake Dynamics
The course of the weather during the study was the key factor affecting total N uptake (TNU) during consecutive stages of onion growth (
Table 5). In general, the highest TNU, except on the 70th DAE, was observed in 2009. The greatest differences between years were found at the beginning of onion growth. A significantly lower N uptake was recorded in 2011 in comparison to 2009 and 2010. The same relationships were again revealed on the 50th DAE. At harvest, the TNU course was in accordance with onion yield (
Table 2).
The bulb N index (BNI) increased in a linear mode from about 40% at the onset of bulbing to about 70% at the 80th DAE in 2009 and 2010, when it reached its maximum. A slight drop was recorded between the 80th and 90th DAE. In 2011, this index showed a very steep increase from 40% on the 50th DAE to 85% on the 90th DAE (
Table 5).
During the growing season, except for the 60th DAE, both TNU and BNI significantly depended on the N × S interaction (
Table 6). Sulfur application was a factor that considerably increased N uptake compared to the absolute control (AC) plot. The differences between experimental treatments depended on the sampling date. At harvest, a significantly higher TNU compared to the AC plot was recorded for the S
30 + N
60, S
60 + N
60 and S
60 + N
120 treatments.
The effect of fertilization on BNI values was also DAE dependent (
Table 6). The highest value of this index was recorded on the 50th DAE on the AC plot, and the lowest on plots fertilized with S alone and plots with the highest N dose (S
30 + N
120 and S
60 + N
120). At harvest, the highest BNI was found for the S
30 + N
120 plot. A significant drop in the index value in the period extending from the 80th and 90th DAE was recorded for the AC plot.
As with TDS, the highest dynamics of N uptake (NC
m) were recorded in 2009 (
Table 7). This index was 3.25-fold and 1.63-fold higher compared to 2010 and 2011, respectively. The same trend was observed for N
Tb, i.e., the day when the N uptake trend changed from exponential to linear and predicted N uptake (TNU
Tb). The differences in TNU
Tb were extremely pronounced, reaching 3.03 g m
−2 in 2009, but only 0.11 g m
−2 in 2010. All of these parameters were significantly correlated with each other. In addition, they were positively but not significantly correlated with yield (
Table S1). It is necessary to stress that NC
m and TNU
Tb were considerably correlated with C
m (r = 0.95 ***). A completely reverse annual trend was observed for NR
m, a parameter that describes the maximum rate of N uptake during the exponential phase of onion growth. The highest value of this parameter was recorded in 2010, i.e., the year with the lowest yield. It was significantly but negatively correlated with parameters of the linear phase of the developed model, but its impact on onion yield was negligible.
3.4. Sulfur Uptake Dynamics
Total sulfur uptake (TSU) by onion in consecutive stages of plant growth was significantly affected by the year. The observed variability was, however, DAE specific (
Figure 4). As a rule, TSU, with the exception of the 70th DAE, was the highest in 2009. The lowest values, except for the 30th DAE, were recorded in 2010. The maximal differences between years, depending on DAE, ranged from 38.8% (DAE 70) to 368% (DAE 90). The differences in TSU observed at harvest were in accordance with the differences in bulb yield. Moreover, the observed patterns of TSU were very similar to those described for TDW (
Figure 3). A key difference refers to trends in S accumulation in leaves. Its decrease was recorded on the 60th DAE in 2009, and on the 70th DAE in 2010. In spite of this, S accumulation showed a tremendous increase following the 70th DAE in 2009. In 2010, only a slight increase in S in a bulb was observed, which progressed up to the 80th DAE and then decreased. A significant decrease of S in leaves followed the 70th DAE. A completely different pattern was observed in 2011, because up to the 70th DAE, the amount of S in the leaves was higher than in the bulb. During the next stages, it suddenly decreased, subsequently resulting in a progressive increase in bulb S.
The bulb S index (BSI), which describes the percentage of TSU in onion bulbs, was the lowest at the onset of onion bulbing (DAE 50), increasing progressively in accordance with crop development (
Table 8). The greatest differences between years were observed in the early stages of onion bulb growth. The slowest increase was observed in 2011, which was significantly lower, especially when compared with 2009. In that year, 89% of the S taken up by onion finished in the bulb.
Values of TSU were significantly modified by the interaction of N and S (
Table 9). However, a significant impact of fertilization treatment was recorded in the S
60 + N
60 treatment on the 60th DAE. At harvest, the highest TSU values, significantly higher compared to AC, were recorded for the same set of treatments as on the 60th DAE and for the S
60 + N
120 treatment. A significant impact of fertilization on the BSI was only observed on the 60th DAE. Its lowest value was recorded on the S
30 + N
60 plot, and the highest on the S
30 + N
0 plot. The BSI values at harvest were non-significant and ranged from 74% to 83%. The highest value was recorded on the S
30 + N
0 plot.
The in-season TSU trend was completely different from TDW and TNU, fitting the quadratic regression model best (
Table 10). TSU parameters, analogous to the dynamic parameters described for onion DW and TNU, were variable year to year. The absolute S uptake (SC
m) and predicted S accumulation (TSU
t2) reached their highest values in 2009. In that year, the two parameters were 4-fold and 1.7-fold higher compared to 2010 and 2011, respectively. The opposite trend was obtained for the relative S uptake rate (SR
m).
The impact of N doses on SCm was independent of S doses, but there some trends were observed. On the S30 plot, the highest SCm was recorded on the S30 + N60 plot. The N rate of 120 kg ha−1 resulted in a slightly lower SCm increase with respect to N0. On the other S plots, the highest SCm values were recorded on the plot fertilized with an N rate of 120 kg ha−1. This parameter displayed the highest impact as well as a positive impact on TDW (r = 0.88 ***). In addition, it was significantly correlated with NCm (r = 0.82 ***).
TSUT2 increased in accordance with increasing S rates. The effect of N rates was, in general, only related to the respective N controls. This parameter of the TSU model showed the same level of correlation with yield (r = 0.87 ***) as SCm. The maximum relative S uptake parameter showed an opposite trend to the previously described parameters of the model. The impact of fertilization was weakly differentiated, ranging from 0.110 to 0.138 mg mg−1 d−1 for S30 + N0 and S30 + N120, respectively. This parameter did not show a significant correlation with any of the studied parameters of the developed model.