Bulb Size Interacts with Lifting Term in Determining the Quality of Narcissus Poeticus L. Propagation Material

: Geophytes are an increasingly widely used group of plants. Its high biodiversity, however, signiﬁcantly impacts the cultivation process. Therefore, every new insight into the behavioral aspects of each single species is crucial for production. This study, for the ﬁrst time, aims at contributing to a better understanding of the interplay between the bulb size potential of Narcissus poeticus and the pattern of bulb lifting time to ﬁnd the optimal combination of these factors in obtaining the best quality propagation material. Twelve bulb size groups were examined. Plants were planted in the open air in two locations of Central East Europe and lifted in three terms: immediately and two and four weeks after ﬂowering. Bulb fresh and dry weight, as well as leaf, root and total plant dry weight, increased with growing bulb size and delayed the time of lifting. The bulb weight increase indicator was twice as high for bulbs lifted in the third date compared to those lifted earlier. Its value was the highest for the smallest bulbs. Bulb shape ratio gradually decreased together with the increase of bulb weight and later lifting term. Principal components analysis allowed singling out two principal components that accounted for 86.58% of total variability.


Introduction
Species belonging to the genus Narcissus (Amaryllidaceae) are known and appreciated ornamental plants worldwide. The most popular are cultivars of the Pseudonarcissus section (e.g., N. pseudonarcissus L.) and Hermione section (N. tazetta L.) [1,2]. Narcissus poeticus cultivars have been cultivated as ornamentals since the beginning of the 20th century and nowadays are still clearly desirable because of their unique decorative and sensory values. The species is included in the Narcissus section [1]. It is characterized by the latest flowering date within the genus, pure white perianth segments and a very short light-yellow corona with a distinct reddish edge. N. poeticus is a spring-flowering plant, producing usually one fragrant flower (uncommon feature) per stem [3]. The species shows a high horticultural potential. Its use as an ornamental bulbous plant for growing in the ground, containers or cut flowers is not the only one. In plant breeding programs, N. poeticus is a parent form in obtaining multi-flowered hybrids and large-cupped narcissi [4]. It found an application in the pharmaceutical industry as a source of biologically active compounds, such as alkaloids, including galantamine [5], which alleviates the symptoms of Alzheimer's disease [6]. Bulbs of Narcissus poeticus are also a source of narpoetan-a growth-stimulating substance recommended for pre-sowing treatment of spiked cereals and cotton in order to stimulate seedlings growth [7]. The invigorating scent of the absolute extracted from N. poeticus Table 1. Characteristics of bulb size groups and bulb spacing in the field experiment. 5.0 < r ≤ 5.5 5.3 20.5 12 S11

Field Experiment
The study was conducted in the 2017-2018 growing season under field conditions in two locations-experimental stations of Mendel University in Lednice in Czech Republic (48 • 47 59,70" N 16 • 48 12,22" E) and University of Agriculture in Krakow in Poland (50 • 03 41" N 19 • 56 18" E). In the last week of October 2017, the bulbs were planted in the experimental plots in both locations. Spacing among rows was 0.5 m, and bulbs' spacing depended on the bulb size (Table 1). Based on a general rule, planting depth comprised twice of bulb height. During the experiment, plants were cultivated in accordance with the typical practice [21,24]. From the beginning of the vegetative season 2018, observations of aboveground development of plants had been recorded. In both locations, only bulbs representing S9-S12 size groups flowered and hence were adopted for all bulb size groups as the determinants for defining the flowering period in order to assess the bulbs' lifting terms. Flowering began in Lednice in the last week of April (27 weeks from planting time) and in Krakow two weeks later. It lasted for three weeks until the last bulb ended flowering, which indicated the end of the flowering season. Bulbs were lifted in 2-week intervals in three terms: immediately after the end of flowering (LT1), two weeks after the end of flowering (LT2) and four weeks after the end of flowering (LT3).
The experiment consisted of 72 combinations (3 lifting terms × 12 bulbs size groups × 2 localizations) with 3 replications, each of which encompassed 23-37 bulbs (depending on the bulb size and spacing), 7300 in total. To observe the morphological parameters from each replication, 5 plants were randomly selected. The layout in particular field localizations was of split-block design.

Climatic and Soil Condition
Central Europe has a continental climate with a warm summer but also with a cold season (Dfb according to Köppen classification) [25], and it covers both experimental locations. Lednice is located around 263 km in a straight line southwest from Krakow. Detailed data concerning climatic conditions such as an average daily temperature (ADT), minimal (MinDT) and maximal (MaxDT) daily average temperature and the sum of precipitation (SP) are gathered as weekly averages or sum (respectively) and presented in Figure 1 (Lednice) and Figure 2 (Krakow). Data were collected for Czech Republic from synoptic station of Mendel University and for Poland from the IMGW (Polish Institute of Meteorology and Water Management) synoptic station in Krakow-Balice and synoptic station of University of Agriculture. In general, the weather conditions in both locations during the experiment were similar and moderate with slightly higher daily temperatures occurring in Lednice. Autumn and winter were warm and humid, while spring was warm with two dry periods observed in April and in the beginning of June.

Plant's Parameters Measurement
In each term of bulb lifting, immediately after harvest, final bulb fresh weight (bulb FW) was estimated. Straight after, the following morphological parameters were measured: bulbs height matter content form 1.5-2.2%. Based on the soil analysis, the following fertilization in both locations was applied: N 2 × 50 kg/ha, P2O5: 45 kg/ha, K2O: 50 kg/ha.  ADT-average daily temperature; MinDT-minimal daily average temperature; MaxDT-maximal daily average temperature; SP-sum of precipitation; Week no.-week number since planting.

Plant's Parameters Measurement
In each term of bulb lifting, immediately after harvest, final bulb fresh weight (bulb FW) was estimated. Straight after, the following morphological parameters were measured: bulbs height In Krakow and in Lednice, the experiment was located on soil classified as brown type (with loess as the basement complex in Krakow). The soil was deeply cultivated and fertile, with pH from slightly acidic to neutral (pH 6.1-6.8), electrical conductivity (EC) below 0.14 dS m −1 and organic matter content form 1.5-2.2%. Based on the soil analysis, the following fertilization in both locations was applied: N 2 × 50 kg/ha, P 2 O 5 : 45 kg/ha, K 2 O: 50 kg/ha.

Plant's Parameters Measurement
In each term of bulb lifting, immediately after harvest, final bulb fresh weight (bulb FW) was estimated. Straight after, the following morphological parameters were measured: bulbs height (bulb H) and biggest diameter (bulb BD) (for the purpose of further calculations), leaf width (leaf W), length (leaf L) and total area (leaf TA), plants total length (plant TL). Leaf W was measured as the widest leaf width, leaf L as the longest leaf length including etiolated part and leaf TA as the sum of the green parts of leaf area of all leaves of the measured plant. Leaf TA was measured with WinDIAS 3.2. Leaf Image Analysis System (Delta-T Devices, Cambridge, UK). The system comprised a video camera pointing vertically down to an illuminated table and connected to a computer. After calibration to accuracy 0.01 cm 2 , WinDIAS-3 analyzed the colored threshold regions to measure leaf parameters. Plant TL was calculated as a sum of bulb H and leaf L values.
Afterward, plant material was washed and dried in a forced ventilation heater at 50 • C to constant weight. Dry weight of organs-root (root DW), bulbs (bulb DW) and leaves (leaf DW)-was measured on a digital scale with measurement accuracy 0.01 g. Plant total dry weight (plant TDW) was calculated by summing up root DW, bulb DW and leaf DW values.
On the basis of the obtained results, the bulb shape ratio (bulb S ratio) was estimated as bulb H-to-bulb BD ratio, and bulb weight increase indicator (bulb WI indicator) was calculated in accordance to Equation (1)

Statistical Analysis
Statistical analyses were performed using Statisica 13.3 software (TIBCO Software Inc., Palo Alto, CA, USA). In order to verify normality of data distribution, a Shapiro-Wilk test was applied. A combined analysis of variance (ANOVA) followed by Duncan's honest significance test was used to determine the main effects of and interactions between the size group, lifting term and localization at the p ≤ 0.05 significance. To determine the relation between all studied traits, the Pearson's correlation coefficient (r) was calculated. Simple regression analysis was performed for establishing the relationship between the initial bulb FW and the final bulb FW in particular lifting time. To evaluate the relationships between original variables on the basis of morphological and ratios traits, the principal component analysis (PCA) was performed (using standardized data). The first two components (PC1 and PC2) explaining the maximum variance in the datasets were used to generate plots.

Fresh and Dry Weight Organs Parameters
The effect of bulb size and lifting term on the quality of Narcissus poeticus propagation material have to our knowledge not been reported before. In this study, the analysis of the main effects proved that lifting term, bulb size, location and interactions between them significantly influenced dry and fresh weight of plant organs (Table 2).
In Lednice, the highest mean fresh weight (FW) was noticed for the bulbs from S11 and S12 size groups and lifted as the latest (S11-LT3 13.6 g; S12-LT3 14.0 g) ( Figure 3). The lowest mean FW (1.3 g) was observed for the smallest bulbs lifted immediately after flowering (S1-LT1); however, it was not significantly different from FW of S1 bulbs group lifted two and four weeks after flowering. Within each size group S4-S12, the FW statistically increased with each other lifting term. Generally, along with increasing bulb size group and delaying lifting term, the DW of plant organs increased, but these differences were not always of statistical significance. The highest mean values for bulb DW (4.8 g), leaf DW (0.7 g) and plant TDW (6.1 g) were noticed for the combination S12-LT3, which, however, did not statistically differ from S11-LT3 combination. All mentioned parameters were the lowest for Agronomy 2020, 10, 975 6 of 18 S1-LT1 bulbs (0.3 g; 0.08 g; 0.5 g), but in comparison to some other combinations, those values were statistically alike. Similar trends were observed in Krakow conditions ( Figure 4). Most often, with the size of planted bulbs their FW increased after the bulb lifting. The highest mean value of bulb FW was obtained in combination S12-LT3 (10.7 g), while the lowest in S1-LT1/LT2/LT3 (1.4-1.7 g); however, it was significantly alike with the S2-LT1/LT2 combination. Within each bulb size group (except S1-S3), delaying the term of lifting resulted in statistically higher bulb FW. Similarly, the later the lifting term and the bigger the bulbs, the higher the dry weight (DW) of bulb and plant (TDW). The highest mean values of these parameters (3.85 g and 4.45 g, respectively) were indicated for S12-LT3 combinations, wherein plant TDW was not significantly different from S12-LT2. Bulb DW and plant TDW of S12-LT3 were 7-8 times higher than those noticed for S1-LT1/LT2/LT3. According to Rees [26], the decrease or increase of Narcissus bulb weight is related to the carbohydrate status of the plant. Moreover, the weight of harvested bulbs is higher in plants with more leaves (e.g., in Lilium) [27] and more leaves per plant are usually produced by larger bulbs than by smaller ones (e.g., in Lachenalia) [13]. In our study, an increase in leaf DW was observed along with an increase in the size of planted bulbs. However, this tendency was not of a model nature. Table 2. Results of ANOVA for the analyzed parameters affected by lifting time, bulb size, localization and interactions.
As reported [19], very early narcissus bulb lifting term or leaves defoliation conducted two to four weeks after flowering have a significant and detrimental effect on bulbs' yields. Based on the results obtained, it can be hypothesized (stated) that in the case of the smallest bulbs (S1-S2), the date of lifting does not affect the value of the dry matter of the organs (which determines its quality), and therefore they can be lifted at any convenient time within 4 weeks from flowering. In the case of larger bulbs, the time of lifting is an important intervention-it is worth leaving the bulbs longer in the ground to get better quality organs. However, in determining the time of bulb lifting, weather conditions in the growing season should be taken into account. Prolonged moist periods favor the development of fungal diseases (e.g., Fusarium, Botrytis, Penicillium) which may infect narcissus bulbs staying too long in the ground [4]. Late bulb lifting may also increase viruses' infections even by three-fold [28]. In view of the mentioned difficulties, the priority is the health of the material obtained, not its size.
Rees [29] indicated that during the spring development of N. pseudonarcissus "Fortune" and "Golden Harvest", and N. poeticus "Actaea", the curves of growth for bulb dry weight and total plant dry weight are close to parallel for all varieties. In the presented study, the noticed relation between those curves in the corresponding period of cultivation is similar, but not so obvious. The mentioned author, however, studied larger, round and double-nosed bulbs material, which suggests that in the case of smaller bulbs, the curves of growth for plant organs may be slightly different.          Root DW grew regularly as the size of the planted bulbs increased, reaching eight times higher values in the S12 group compared to the smallest bulbs of S1 ( Figure 5). In the case of tulips, a similar tendency was observed-better quality of root system (roots length, weight and number) was typical for bigger bulbs [30]. In this study, the mean dry weight of roots developed by the bulbs lifted as the earliest was the lowest (0.21 g), whereas that produced by the plants lifted later (LT2-LT3) was about 20% higher compared to LT1. Such an increase of N. poeticus root DW may be explained by a large increase in fatty acids that was observed in the roots of N. pseudonarcissus towards the end of the season [31]. Strong and positive correlation between root DW and bulb weight parameters (FW, DW) was observed (r = 0.899 and r = 0.873, respectively, at p ≤ 0.01) ( Table 3). Root DW grew regularly as the size of the planted bulbs increased, reaching eight times higher values in the S12 group compared to the smallest bulbs of S1 ( Figure 5). In the case of tulips, a similar tendency was observed-better quality of root system (roots length, weight and number) was typical for bigger bulbs [30]. In this study, the mean dry weight of roots developed by the bulbs lifted as the earliest was the lowest (0.21 g), whereas that produced by the plants lifted later (LT2-LT3) was about 20% higher compared to LT1. Such an increase of N. poeticus root DW may be explained by a large increase in fatty acids that was observed in the roots of N. pseudonarcissus towards the end of the season [31]. Strong and positive correlation between root DW and bulb weight parameters (FW, DW) was observed (r = 0.899 and r = 0.873, respectively, at p ≤ 0.01) ( Table  3).

Morphological Parameters of Plants and Leaves
The impact of the studied factors was not as significant for plant TL and leaf L as it was for DW and FW parameters. In both locations, plant TL is strongly correlated with leaf L (Table 3, Figure 6, Figure 7). In Lednice, the tendency to increase leaf and plant TL with delaying the lifting term was most often observed. Results show that leaf W widens with increasing bulb size both in Lednice and Krakow, generally in accordance with the rule that the bigger the mother bulb, the wider the leaf. The interactions between leaf W and lifting term were noticed (

Morphological Parameters of Plants and Leaves
The impact of the studied factors was not as significant for plant TL and leaf L as it was for DW and FW parameters. In both locations, plant TL is strongly correlated with leaf L (Table 3, Figure 6, Figure 7). In Lednice, the tendency to increase leaf and plant TL with delaying the lifting term was most often observed. Results show that leaf W widens with increasing bulb size both in Lednice and Krakow, generally in accordance with the rule that the bigger the mother bulb, the wider the leaf. The interactions between leaf W and lifting term were noticed ( Table 2), but no clear tendency can be defined.  Taking into consideration the average results of leaf TA measurements from two locations, it was found that initially, larger bulbs gave rise to plants whose leaf TA increased more than in plants originating from initially smaller bulbs ( Figure 8). This is in agreement with Rees et al. [32], who report that the bulb size determines the leaf area. Comparison of the results in the context of lifting term shows that the leaf TA value had been increasing up to two weeks after flowering (LT2) and decreasing afterwards, reaching the lowest level in LT3. In the period preceding the second term of bulb lifting (LT2), the highest values of leaf TA (LT2) (Figure 8) could have coincided with the peak of incoming solar radiation. Correlation coefficient indicated that leaf TA was positively associated with the most of attributes studied, except bulb WI indicator (Table 3). Moreover, that strong positive correlation with leaf TA concerned parameters connected with the quality of obtained bulbs (bulb FW, r = 0.819; bulb DW, r = 0.781; at p ≤ 0.01). These data agree with the results reported by Rees [17] and Weraduwage et al. [33] and indicate that leaf area growth is favorable for dry matter production and the light interception and may be used as the criterion of plant productivity. Our experiment shows that Narcissus poeticus leaves (TA) and roots (DW) exhibit active growth until two weeks from flowering, while growth of bulbs (DW, FW) still continues (Figures 3-8). Taking into consideration the average results of leaf TA measurements from two locations, it was found that initially, larger bulbs gave rise to plants whose leaf TA increased more than in plants originating from initially smaller bulbs (Figure 8). This is in agreement with Rees et al. [32], who report that the bulb size determines the leaf area. Comparison of the results in the context of lifting term shows that the leaf TA value had been increasing up to two weeks after flowering (LT2) and decreasing afterwards, reaching the lowest level in LT3. In the period preceding the second term of bulb lifting (LT2), the highest values of leaf TA (LT2) (Figure 8) could have coincided with the peak of incoming solar radiation. Correlation coefficient indicated that leaf TA was positively associated with the most of attributes studied, except bulb WI indicator (Table 3). Moreover, that strong positive correlation with leaf TA concerned parameters connected with the quality of obtained bulbs (bulb FW, r = 0.819; bulb DW, r = 0.781; at p ≤ 0.01). These data agree with the results reported by Rees [17] and Weraduwage et al. [33] and indicate that leaf area growth is favorable for dry matter production and the light interception and may be used as the criterion of plant productivity. Our experiment shows that Narcissus poeticus leaves (TA) and roots (DW) exhibit active growth until two weeks from flowering, while growth of bulbs (DW, FW) still continues (Figures 3-8).

Shape Ratio and Weight Indicator as Bulbs Quality Determinants
The bulb S ratio provides information about the shape of a bulb. Its value closer to 1 indicates an approximately round (globular) shape, whereas values receding from unity show bulb elliptical shape. Predominantly in both locations bulb, S ratio gradually decreases together with the increase of bulbs weight (higher size group) and later lifting term (Figure 9, Figure 10). Bigger bulbs have lower bulb S ratio values in comparison to the smaller bulbs. This trend is confirmed by correlation analysis, presented in Table 3, revealing a strong, negative relationship between the bulb S ratio and bulb FW (r = −0.776, p ≤ 0.01). Similarly, strong negative correlations occur in relations of bulb S ratio LE S1 S2 S3 S4 S5 S6 S7 S8 S9 S10 S11 S12 0 10 20 30 40 50 KR S1 S2 S3 S4 S5 S6 S7 S8 S9 S10 S11 S12 LT3

24.3B a ab a-c c-e c-e f-j h-j h-k i-l i-l i-m k-m a-c a-c a-d d-f e-h f-j h-k i-l j-m l-m m m ab a-d a-e b-e d-g g-j h-j i-l l-m k-m lm lm a a-c b-d d-f i-l j-m k-o o-s tu m-p o-s f-i ab a-c d-g j-m g-j l-p m-p r-t o-s o-s u tu ab c-e e-h g-j g-j g-j h-k j-n k-o p-t k-o n-r
Size group

Shape Ratio and Weight Indicator as Bulbs Quality Determinants
The bulb S ratio provides information about the shape of a bulb. Its value closer to 1 indicates an approximately round (globular) shape, whereas values receding from unity show bulb elliptical shape. Predominantly in both locations bulb, S ratio gradually decreases together with the increase of bulbs weight (higher size group) and later lifting term (Figure 9, Figure 10). Bigger bulbs have lower bulb S ratio values in comparison to the smaller bulbs. This trend is confirmed by correlation analysis, presented in Table 3, revealing a strong, negative relationship between the bulb S ratio and bulb FW (r = −0.776, p ≤ 0.01). Similarly, strong negative correlations occur in relations of bulb S ratio with root DW, bulb DW and plant TDW. The results discussed above correspond to the findings of Hanks [21], who claims that large round bulbs are obtained in one year from larger offsets, while small offsets need two or even three years to achieve a large round shape. According to Marques and Draper [34], bulb size (which in the case of narcissus is combined with the shape) indicates the degree of bulb maturity, and this feature is genotype-dependent. Moreover, Pütz [35] claims that illumination may play a role in determining bulb shape. The illumination of the basal parts of the leaf of N. tazzeta resulted in obtaining the bulbs with oval or elliptical shape, whereas in not illuminated plants, more globular bulb shape was observed. Additionally, in the present study, the initial bulbs of elliptical shape were able to flower (see Section 2) however only those of adequate high weight (4.8-6.3 g).  The course of the bulb WI indicator values distribution for Lednice and Krakow is presented in Figure 11. In Krakow conditions, the values of this parameter obtained in the third lifting term were different from those determined in LT2 and LT1, however not always relevant. Results for Lednice revealed a similar tendency, and the observed differences were statistically significant. Bulb WI Indicator successively increased with delaying the lifting term. In both locations, bulb WI indicator parameters were twice as high in the third lifting term as those recorded in the LT1. Bulb WI indicator was significantly the highest (115.3%) for the smallest bulbs (S1) and the lowest (65%) for the biggest ones (S12). The latter mentioned value did not differ significantly from the ones recorded for the S11 (75%) and S10 (74.9%) bulbs. In practice, these results mean that the smallest bulbs g   The course of the bulb WI indicator values distribution for Lednice and Krakow is presented in Figure 11. In Krakow conditions, the values of this parameter obtained in the third lifting term were different from those determined in LT2 and LT1, however not always relevant. Results for Lednice revealed a similar tendency, and the observed differences were statistically significant. Bulb WI Indicator successively increased with delaying the lifting term. In both locations, bulb WI indicator parameters were twice as high in the third lifting term as those recorded in the LT1. Bulb WI indicator was significantly the highest (115.3%) for the smallest bulbs (S1) and the lowest (65%) for the biggest ones (S12). The latter mentioned value did not differ significantly from the ones recorded for the S11 (75%) and S10 (74.9%) bulbs. In practice, these results mean that the smallest bulbs increase their volume more than once, while larger by only more than half. Thus, the important S1 LT1 LT2 LT3   The course of the bulb WI indicator values distribution for Lednice and Krakow is presented in Figure 11. In Krakow conditions, the values of this parameter obtained in the third lifting term were different from those determined in LT2 and LT1, however not always relevant. Results for Lednice revealed a similar tendency, and the observed differences were statistically significant. Bulb WI Indicator successively increased with delaying the lifting term. In both locations, bulb WI indicator parameters were twice as high in the third lifting term as those recorded in the LT1. Bulb WI indicator was significantly the highest (115.3%) for the smallest bulbs (S1) and the lowest (65%) for the biggest ones (S12). The latter mentioned value did not differ significantly from the ones recorded for the S11 (75%) and S10 (74.9%) bulbs. In practice, these results mean that the smallest bulbs increase their volume more than once, while larger by only more than half. Thus, the important practice is to wait with lifting the largest bulbs as their increase during the vegetative season is not too spectacular.
Further results indicate a positive association of bulb WI indicator with leaf and total plant length as well as with leaf DW (r = 0.486, r = 0.452, r = 0.225 respectively at p ≤ 0.01) ( Table 3). It has previously been proposed for some ornamental geophytes that the size of planting material is important for determining the increase in weight of bulbs-greater weight gain during the season was observed for smaller plant materials than for bigger ones [17,36,37]. A similar observation concerning higher values of bulb weight increase shown by smaller bulbs of Narcissus 'Fortune' compared to larger ones was noticed by Fenlon et al. [38]. This relation turned out to be independent of method of propagation from which bulbs were obtained. This can be explained by the fact that larger bulbs are closer to achieve maximum final size than smaller ones [16]. The general remark that arises after analysis of all parameters presented in sections 3.1-3.3. is that in both locations, similar tendencies were observed, although the values of particular parameters sometimes slightly differ between locations and were higher in Lednice. Subtle differences in the temperature and precipitation pattern could have affected the individual values of the studied parameters ( Figure 1, Figure 2, Table S1).

Linear Response of the Bulb Fresh Weight (Bulb FW) under Different Lifting Terms (LT)
Linear response of the final bulb FW to the initial bulb fresh weight and the three different lifting terms are shown in Figure 12. When bulbs are lifted just after flowering (LT1) the final bulb FW is expected to increase by 1.36 of the initial bulb fresh weight, while lifted two (LT2) and four (LT3) weeks later, by 1.64 and 1.94, respectively. The R 2 value showed that 97, 88 and 85% of variation in final bulb FW for subsequent lifting terms could be explained by a linear regression model.

S1
S2 S3 S4 S5 S6 S7 S8 S9 S10 S11 S12 Size group Location Mean LE-KR S1-S12  The general remark that arises after analysis of all parameters presented in Sections 3.1-3.3 is that in both locations, similar tendencies were observed, although the values of particular parameters sometimes slightly differ between locations and were higher in Lednice. Subtle differences in the temperature and precipitation pattern could have affected the individual values of the studied parameters ( Figure 1, Figure 2, Table S1).

Linear Response of the Bulb Fresh Weight (Bulb FW) under Different Lifting Terms (LT)
Linear response of the final bulb FW to the initial bulb fresh weight and the three different lifting terms are shown in Figure 12. When bulbs are lifted just after flowering (LT1) the final bulb FW is expected to increase by 1.36 of the initial bulb fresh weight, while lifted two (LT2) and four (LT3) weeks later, by 1.64 and 1.94, respectively. The R 2 value showed that 97, 88 and 85% of variation in final bulb FW for subsequent lifting terms could be explained by a linear regression model. In the presented experiment, flowering bulbs weighted from 4.8 g (S9). Theoretically, considering the cultivation of small, for example 2 g, bulbs, based on the presented linear regression equations (Figure 12), delaying lifting term for 4 weeks enabled reaching the target weight for flowering bulb earlier. Thus, this approach may improve the production process, and from the point of view of producers, it is of great importance. However, as other authors report, each species, or even cultivar, may have specific critical bulb parameters for flowering [4,13,29,39,40].

Principal Components Analysis
Principal components analysis provides a new artificial variable (PCs) calculated on the basis of linear combinations of original variables [41]. In this study, the first two principal components evaluated for 11 tested traits identified 86.58% of total variability (Table 4, Figure 13). The other PCs, that had a minor effect on the model, were discarded.   In the presented experiment, flowering bulbs weighted from 4.8 g (S9). Theoretically, considering the cultivation of small, for example 2 g, bulbs, based on the presented linear regression equations (Figure 12), delaying lifting term for 4 weeks enabled reaching the target weight for flowering bulb earlier. Thus, this approach may improve the production process, and from the point of view of producers, it is of great importance. However, as other authors report, each species, or even cultivar, may have specific critical bulb parameters for flowering [4,13,29,39,40].

Principal Components Analysis
Principal components analysis provides a new artificial variable (PCs) calculated on the basis of linear combinations of original variables [41]. In this study, the first two principal components evaluated for 11 tested traits identified 86.58% of total variability (Table 4, Figure 13). The other PCs, that had a minor effect on the model, were discarded.  Figures 14 and 15 show the location on the plane of the first two main components, presented in relation to bulb size group (S1-S12) and lifting term (LT1-LT3). Comparison of Figures 13 and 14 and      Although positive interaction of almost all studied features of plants with location factor was observed (Table 2), the similar tendency in growth and development of bulbs both in Lednice and Scatter plot (Total -86.58%) S1 S1 S1 S1 S1 S1 S1 S1 S1 S1 S1 S1 S1 S1 S1 S1 S1 S1 S10 S10 S10 S10 S10 S10 S10 S10 S10 S10 S10 S10 S10 S10 S10 S10 S10 S10 S11 S11 S11 S11 S11 S11 S11 S11 S11 S11 S11 S11 S11 S11 S11 S11 S11 S11 S12 S12 S12 S12 S12 S12 S12 S12 S12 S12 S12 S12 S12 S12 S12 S12 S12 S12 S2S2   Although positive interaction of almost all studied features of plants with location factor was observed (Table 2), the similar tendency in growth and development of bulbs both in Lednice and Scatter plot (Total -86.58%) S1 S1 S1 S1 S1 S1 S1 S1 S1 S1 S1 S1 S1 S1 S1 S1 S1 S1 S10 S10 S10 S10 S10 S10 S10 S10 S10 S10 S10 S10 S10 S10 S10 S10 S10 S10 S11 S11 S11 S11 S11 S11 S11 S11 S11 S11 S11 S11 S11 S11 S11 S11 S11 S11 S12 S12 S12 S12 S12 S12 S12 S12 S12 S12 S12 S12 S12 S12 S12 S12 S12 S12 S2S2  Although positive interaction of almost all studied features of plants with location factor was observed (Table 2), the similar tendency in growth and development of bulbs both in Lednice and Krakow may suggest that probably other factors besides environmental factors influenced plant behavior. According to Rees [29], the size of a bulb has a greater impact on the way of plants behavior during development process rather than other attributes such as a cultivar type. Seasonal changes in the primary and secondary metabolite profiles observed in the bulbs, leaves and roots which are involved in organogenesis processes, stress response, hormone transport, plant architecture and in many other phases of plant development [31] might have also influenced the obtained results.

Conclusions
In the presented study we proved that harvesting bulbs of Narcissus poeticus at an early time-point is possible; however, bulbs of a higher quality are obtained with the delay of the harvest date. Thus, taking into account the needs of the producers and consumers who expect to obtain plants of the highest quality, we recommend leaving the bulbs in the ground four weeks after flowering because they achieve the higher quality parameters. Derogation from this recommendation may apply to the bulbs with the smallest initial weight, on which the lifting term does not have as much impact as on larger bulbs. Of all studied factors, bulbs size appeared to determine the dynamics of plants growth and influence the final plant parameters the most. We found that root dry weight and leaf total area can be adopted as the criterion of plant productivity as they were positively correlated with bulb weight parameters. Based on the results presented here, a rule can be formulated-the smaller the mother bulb, the greater weight gain during the season. The obtained results relate to a moderate climate, which diverges from the optimal climatic conditions for N. poeticus. Such a climatic area of research broadens the knowledge about this species in an ecological and geographical context and may form part of new horticultural practices or the specialized protocols used for pharmaceutical or other industrial targets.