1. Introduction
Ornamental sunflower (
Helianthus annuus L., Asteraceae) includes genotypes for using as garden plants, cut flower, and potted plant [
1]. Its beauty and profitability potential has attracted the interest of producers and entrepreneurs seeking new opportunities for business because of several properties such as short developmental cycle, ease of propagation, plant hardiness, and its attractive inflorescence, widely used in floral arrangements [
2,
3]. These properties are aligned with the species of flowers in the “Flowers for All” project [
4,
5], a nationwide project that aims to promote floriculture among small and medium producers as a source of income for rural communities, considering the pillars of economic, environmental, and social sustainability. In this context, cut sunflower is used as one of the main species cultivated by the “Flowers for All” project teams across Brazil [
3].
The ornamental quality of cut sunflower is influenced by various factors, among which plant density (the number of plants per unit area) and location where it is cultivated. Plant density affects intraspecific competition for resources such as soil nutrients, soil water, and solar radiation so that an inadequate plant density can affect agronomic performance of the crop and flower quality [
6,
7,
8,
9]. Location is responsible for providing different environments that affect plant development and growth, which may result in consequences on the quality of cut sunflower. Previous studies have shown the effects on agronomic performance of sunflowers due to the cultivation location, influencing the morphophysiological characteristics of the plant and crop developmental cycle because of different meteorological variables among locations such as precipitation, solar radiation, and above ground and soil temperature [
10,
11,
12]. Therefore, the hypothesis in this study is that there is a plant density x location interaction and that the optimum agronomic plant density may differ among locations in cut sunflower.
Most of the flowers in Brazil are produced in the state of São Paulo [
13]. In a large country like Brazil, the concentration of production in a single state or region has advantages e disadvantages. Among the disadvantages is that the chain from producers to consumers is long, which impacts directly on flowers price for consumers. Therefore, it is important to identify other regions of Brazil, which has subtropical, tropical and equatorial climates, with potential to produce high quality flowers and therefore shorten the chain. The objective in this study was to evaluate the effect of plant density and cultivation location on the development, growth, and quality of cut sunflower in a tropical and in a subtropical environment.
3. Discussion
Comparing the agronomic performance of cut sunflower in the two locations, which represent the tropics (Dourados/MS) and subtropics (Santa Maris/RS) of Brazil, that the environmental conditions in Dourados/MS resulted in higher PH, SD, LA, LAI, and PHY. This trend suggests that the tropical environment, characterized by higher temperatures and greater solar radiation, promotes faster growth and more efficient radiation use by plants.
Plant density played a crucial role in the results. In general, plants cultivated at a lower density (D10) showed greater CD, SD, FLN, and LA, while the highest density (D50) favored LAI and PHY. This reflects the complex dynamics of intraspecific competition, where plants at a lower plant density can better use the available resources, leading to greater growth.
The cut sunflower cultivated in Santa Maria /RS stood out in terms of CD and FLN. The higher meteorological variability and cooler temperatures in the Subtropics of Brazil seem to have positively influenced these variables, resulting in a larger capitulum diameter and number of leaves, with a longer developmental cycle.
The location x plant density interaction shown in this study revealed interesting results for CD and FLN. In Dourados/MS, sunflower plants cultivated at D10 showed the highest values of CD, while in Santa Maria/RS, this density and D20 also stood out for FLN. This result indicates that the choice of plant density should be made based on the specific cultivation location, considering the environmental and climatic conditions.
The cultivation location of the cut sunflower hybrid Vincent´s Choice significantly influenced the growth variables, confirming that the crop response is linked to the cultivation environment and plant density. Compared with previous studies, such as those by [
14] in Santa Catarina and [
15] in Rio de Janeiro, Brazil, with sunflower genotypes for grain and oil where the sunflower plant heigth was 178 and 170 cm, respectively, the values found in our study were lower, highlighting the variability of the culture's development in different conditions and cultivation locations.
Plant height was one of the variables most influenced by the cultivation location, being greater in plants when cultivated in Dourados/MS. This location had a shorter photoperiod throughout the sunflower developmental cycle, suggesting that favorable conditions, such as temperature accumulated, radiation, and precipitation contributed to accelerated plant growth.
Other growth variables CD, SD, FLN, LA, LAI, are also important for crop development and growth. In this study, these variables were influenced by location, with CD and FLN being greater in Santa Maria/RS, which can be attributed to the greater oscillations in solar radiation and cooler temperatures, demonstrating the influence of these meteorological variables on development and growth variables of sunflower plants.
The SD, along with CD and PH, are variables that indicate the commercial value of cut sunflower. The values found in this study, especially at lower plant densities, align with the quality standards of cut sunflower floral stems of the Veiling Holambra Cooperative, indicating that the produced stems meet the market demands. The phyllochron is an important variable of plant development during the vegetative phase as is a measure of the rate that new leaves appear on the stem, which are essential for the interception and absorption of solar radiation. A higher phyllochron indicates a longer time interval for a new leaf to appear and vice-versa [
16]. In this study, it was observed a higher phyllochron in Dourados/MS and, consequently, a longer time interval for each leaf to appear. This explains why the final number of leaves was lower in Dourados than in Santa Maria.
The temperature difference between the two locations, along with plant density, can significantly influence the crop developmental cycle. The increase in plant density can lead to greater shading, directing more photoassimilates towards plant height growth up to a limit, after which resource scarcity can compromise development [
17,
18]. Air temperature directly influences the phenology of agricultural crops, especially as it is associated with the rate of photosynthesis, likewise affecting the phyllochron, influencing the crop developmental cycle [
19]. In Dourados/MS, the higher temperatures and reduced photoperiod during the crop developmental cycle may have contributed to accelerating its development by about 5 days. However, when the duration of the developmental cycle is in thermal time (
oC day), the plants cultivated in Dourados/MS presented a longer cycle. This demonstrates the importance of defining which time (calendar days or
oC day) is used for defining the duration of developmental phases.
In both locations, the duration of developmental cycle of the cut sunflower was shorter than for sunflower genotypes grown for grain and oil purposes. While our results varied from 49 to 52 days, [
2], working with the cultivars BRS Oasis, BRS Refúgio M, and BRS Paixão M in Chapecó - Santa Catarina, reported a developmental cycle varying from 59 to 81 days. [
20] noted an inverse relationship between air temperature and developmental cycle duration, that is, higher average air temperatures during the cycle tend to reduce its duration.
Plant density affects intraspecific competition, determining variables of economic yield [
7]. Or results indicate that a plant density of 50 plants/m² results in a smaller capitulum diameter and stem, emphasizing the need to consider the growth and production potential of each genotype in different edaphoclimatic conditions to determine the optimum agronomic plant density [
21,
22]. [
11] also observed that the plant density in sunflower varies according to regional environmental conditions, influencing crucial characteristics for the commercialization of cut sunflower, such as the capitulum diameter. Higher planting density intensifies intraspecific competition, potentially reducing the availability of nutrients for growth [
17]. [
23] noted a reduction in the number of capitulum within the standards for the sunflower variety IAC Uruguay with the increase in planting density, attributing this to greater growth in height and the number of normal capitulum, intensifying intraspecific competition.
For the final leaf number, [
24] obtained an average of 23 leaves/plant, a value close to that found in this study in both locations. However, these numbers were lower than the 28 leaves reported by [
25] in Campinas, São Paulo.
The results obtained in this study highlight the complex interaction between edaphoclimatic factors and plant density on growth variables in the cultivation of cut sunflower. It has become evident that the cultivation environment plays a crucial role in the development of the crop, with variables such as temperature, radiation, and photoperiod directly affecting the developmental cycle. Moreover, plant density had a considerable impact on commercial characteristics of the flower stem, such as capitulum and stem diameter.
Intraspecific competition at higher plant densities can result in a limited allocation of resources, hindering growth, while lower densities may favor more robust plants at the end of the developmental cycle. These findings are of high practical relevance for cut sunflower producers, offering valuable insights for choosing the optimum plant density and effectively managing edaphoclimatic conditions, aiming at optimizing production. Furthermore, the results emphasize the importance of taking into account regional and environmental variations when planning the production of cut sunflower, adapting cultivation to each location. This strategy can not only enhance the quality of the floral stems but also contribute to the stability of supply throughout the year, positively impacting market prices.
Finally, it is important to emphasize that the focus of this study is on optimizing the management of cut sunflower in tropical and subtropical regions of Brazil, allowing family farmers to cultivate and market ornamental sunflower locally, fostering income generation for these small producers. Therefore, all the stems produced, regardless of plant density and cultivation location, can be marketed without significant losses in the quality of the floral stems. With the data presented in this study, the use of a plant density of 30 plants/m² is recommended, which is pretty close to the plant density of 32 plants/m² used in the “Flowers for All” project [
4,
5], as it allows a better use of space, increasing the number of flower stems per unit area without significant losses in the quality of the floral stems.
4. Materials and Methods
4.1. Study Locations
The study was conducted in two locations: in the experimental area of Horticulture at the Faculty of Agricultural Sciences (FCA) of the Federal University of Grande Dourados (UFGD), in Dourados/MS, in the tropics of Brazil) (22°11’45” S, 54°55’18” W, altitude 446 m) in Central-West Brazil, and in a farm in Santa Maria/RS in the subtropics of Brazil (29°41’29” S, 53°48’ 3” W, altitude 96 m) in Southern Brazil. According to the Köppen classification, the climate of Dourados/MS is Aw, tropical with a rainy season in Summer and a dry season in Winter, and in Santa Maria/RS the climate is Cfa, subtropical humid without dry season and with hot Summer [
26].
4.2. Meteorological Data Collection
Daily meteorological data of minimum and maximum air temperature, precipitation, and solar radiation during the experimental period at Dourados/MS were from a weather station located 9 km from the experiment and missing data were filled from the weather station of the Brazilian National Weather Service (INMET). In Santa Maria/RS the meteorological data were from a weather station of INMET located about 6 km from the farm.
4.3. Plant Material and Experimental Design
The cut sunflower hybrid Vincent´s Choice was used. The experimental design was a 2 x 5 factorial in randomized blocks with four replications. Factor A was location (Dourados/MS and Santa Maria/RS) and Factor B was plant density (D10 = 10 plants/m², D20 = 20 plants/m², D30 = 30 plants/m², D40 = 40 plants/m², and D50 = 50 plants/m²). Each plot (replication) was as 1 m wide x 2 m long (2 m2) raised bed (0.25 m height), with 4 rows of plants, in a plant spacing of 0.20 m among rows and 0.40 m within plants in each row (D10), 0.20 m within plants in each row (D20), 0.13 m within plants in each row (D30), 0.10 m within plants in each row (D40), and 0.08 m within plants in each row (D50).
4.4. Seedling Preparation and Greenhouse Conditions
In both locations, seeds were sown on 15 October 2021 in plastic trays with 128 cells filled with the same commercial substrate, sowing one seed per cell at 1cm depth. The trays were kept inside a low-density polyethylene plastic greenhouse during 10 days in Dourados/MS during and 11 days in Santa Maria/RS to protect the seeds and seedlings from rain and birds until transplanting. During the period inside the plastic greenhouse, the emergence of the seedlings was counted daily until complete emergence (final stand).
4.5. Soil Preparation and Characteristics
The experimental areas were prepared with plowing and harrowing, followed by raising beds 1 m wide and 0.25 m height with a bed-forming rototiller. Soil samplings were taken for chemical and physical tests (
Table 1) at a depth of 0-10 cm, and based on the results of the tests, soil correction was performed. The soil in Dourados is classified as Dystroferric Red Latosol (SiBCS), with a clay texture (60%), while the soil in Santa Maria is a type of transition between typical Aluminic Brown-Ashy Argisol and Arenic Red Dystrophic Argisol.
Table 3.
Soil chemical properties of the sunflower cultivation areas cv. Vincent's Choice in Dourados/MS and Santa Maria/RS.
Table 3.
Soil chemical properties of the sunflower cultivation areas cv. Vincent's Choice in Dourados/MS and Santa Maria/RS.
Soil Chemical Characteristics |
pH in water pH in CaCl2 P (mg/dm3) Ca (cmolc/dm3) Mg (cmolc/dm3) K (cmolc/dm3) Sum of bases (cmolc/dm3) Base saturation (cmolc/dm3) Organic matter (%) |
Dourados – MS* 5.97 5.30 5.90 2.72 1.57 0.50 4.79 55.38 3.57 |
Santa Maria – RS** 4.55 -- 88.30 3.70 1.70 0.52 5.92 49.80 2.10 |
*Methodology Embrapa and IAC; ** Methodology of the Official Network of Soil and Plant Tissue Analysis Laboratories of the States of Rio Grande do Sul and Santa Catarina. |
4.6. Transplanting and Fertilization
The transplanting point was defined when the cotyledonary leaves were open at a 180° angle, the first pair of true leaves had a blade length between 1 and 2 cm, and the clod was well-formed before transplanting, 1 L m-2 of semi-decomposed chicken manure-based organic fertilizer was broadcast in both experimental areas. After incorporating this organic fertilizer, the beds (1 m x 2 m) for each replication were delimited, followed by chemical fertilization of the soil. In Santa Maria, 50 g m² of NPK 05-20-20 was applied on the day of planting, while in Dourados an equivalent fertilization was used, applying 66.7 g m-2 of NPK 10-15-15 on the day of planting. The transplanting was done on 25 October 2021 in Dourados/MS and on 26 October 2021 in Santa Maria/RS, with plants arranged according to the plant density treatments described previously.
4.7. Irrigation and Crop Management
Right after transplanting, mulching with straw was added to the surface of the soil on the beds, and three plants in each of the two central lines of each plot were randomly tagged with colored wires. Ten days after transplanting, side dressing fertilization was performed with 25 g m-2 of urea and 25 g m-2 potassium chloride (KCl). During the experimental period, plants were irrigated as needed with drip tubes for approximately 45 minutes, with a flow rate of 0.003 L min-1. Management practices included manual control of weeds, staking with bamboo stakes and twine at 30, 60, and 90 cm heights, and pest control as needed, using chemical insecticides based on fipronil and biological with Beauveria bassiana fungus, following the manufacturer instructions.
4.8. Growth Measurements and Developmental Stages
Starting four days after transplanting, the number of leaves on the tagged plants was counted twice a week until the final leaf (final leaf number). A leaf was counted when it blades length was greater than 2 cm. Plant height and leaf area were measured once a week on two tagged plants per plot up to the beginning of the reproductive phase. The greatest length (L) and the greatest width (W) of each leaf were measured, and the area of each leaf (LA) was calculated using the equation for sunflower by [
27] as LA = 0.733 (L * W). The plant leaf area was calculated as the sum of the area of each leaf and then the leaf area index (LAI, cm
2 leaf/cm
2 soil) was calculated areas as the plant leaf area divided by the area of soil occupied by each plant [
28]. On the tagged plants, the data of the following developmental stages during the reproductive phase (R) were determined according to the PhenoGlad developmental scale for cut sunflower [
3]: R1 (visible floral bud), R4 (color of the petals of the outer ligulate flowers of the capitulum is visible), and R5 (harvest point – ligulate flowers at 90° with the capitulum disk) (
Figure 4).
Figure 4.
Reproductive stages R1, R4, and R5 in cut sunflower according to the PhenoGlad developmental scale used in the study.
Figure 4.
Reproductive stages R1, R4, and R5 in cut sunflower according to the PhenoGlad developmental scale used in the study.
4.9. Harvest and Quality Evaluation @
Daily photoperiod, considering the duration of civil twilight of 6° below the horizon was calculated using the algorithm by [
29]. The phyllochron (°C day leaf
-1) as an indicator of leaf appearance rate was calculated for each plot by the inverse of the angular coefficient of the linear regression between the average number of leaves of the six plants in the plot and the accumulated thermal time (ATt, °C day) [
30]. The accumulated thermal time (ATt, °C day) was calculated as ATt = ∑TTd where TTd is the thermal time of the day, calculated as TTd = Tmed - Bt, with Tmed being the daily mean air temperature and Bt is the basal temperature of the crop assumed 7.2 °C [
31]. The duration of the developmental phases from sowing to R1 (SO-R1), from sowing to R4 (SO-R4), and from sowing to R5 (SO-R5) was calculated in days and in °C day.
When plants reached the harvest point (R5), the tagged plants in each plot were harvested by cutting them at the soil surface and the stems were cut to standardize floral stems with 0.7 m length (Veiling Holambra standard), leaves were dropped off and floral stems were placed in containers with tap water. The floral stems harvested in Dourados were taken to the Laboratory of in vitro Cultivation of Flowers and Ornamental Plants of UFGD and in Santa Maria to the Crop Science Department of UFSM.
At harvest, the total height of the plant, from the insertion of the stem in the soil to the inflorescence, was the diameter of the stem at the height of the cut and the diameter of the floral capitulum were measured. The floral stems were evaluated for quality, including stem tortuosity, stem diameter, and capitulum diameter. For classification, the quantitative standards of [
32] were used, classifying the floral stems in the standard 70: stem length of 0.70 m, minimum stem diameter of 0.80 cm, and the capitulum diameter (open flower) of at least 6.0 cm. Stems that did not meet one of the criteria for minimum diameter or stem length were classified as non-marketable.
4.10. Statistical Analysis
The data were subjected to analysis of variance, and where significant differences were found, the Bonferroni t-test (p ≤ 0.05) was applied for the Factor A (location) and regression analysis for Factor B (plant density) using the SISVAR software (Statistical Analysis Program v.5.6. Federal University of Lavras - MG) [
33].