Effects of geoclimatic factors on the variability in Pinus pinea cone, seed, and seedling traits in Turkey native habitats

Stone pine (Pinus pinea L.) is harvested for its edible kernels. It is an important forest tree species, and also plays an important role in afforestation. Successful afforestation activities strongly depend on the quality of seeds and seedlings. Five eastern native populations of stone pine were analyzed, in terms of their geographical and climatic parameters, in order to identify their cone, seed, pine nut, and seedling traits. With this regard, we aimed to contribute to the reforestation, yield, and breeding of this species. All native populations occur in the sub-humid bioclimatic zone, according to Emberger’s bioclimatic classification, ranging 74.7–63.2 of Emberger’s Q Index value (Q), which synthesizes humidity and temperature, with lower values indicating more arid conditions. With regard to geographical factors, altitude was found to have the greatest effect on cone diameter and seed weight, while there was a weak correlation of cone, seed, and pine nut dimension with longitude and latitude. In terms of climatic factors, Q and mean of the maximum temperature of June, July, and August (ME) were found to be the main drivers behind significant effects in cone, seed, and pine nut morphometric traits. The most negative effect on seed number cone−1 (r: − 0.897), pine nut number cone−1 (r: − 0.923), pine nut yield (r: − 0.903), and pine nut weight (r: − 0.878) was Q, whereas cone diameter, seed length, and seed weight were positively correlated with the ME (r: 0.939, r: 0.889, and r: 0.866, respectively). Cone and seed diameter and seed weight increased with decreasing summer drought index (r: − 0.806, r: − 0.846, and r: − 0.866, respectively). Rainfall in the June of the third year prior to harvest time (RJ) had a negative effect on the percentage of sound seeds per cone. Annual rainfall in the third year prior to harvest time (AR) significantly affected pine nut weight (r: 0.889), although rainfall during summer had a negative effect. The empty seed ratio of the different populations ranged from 20 to 68%, while considering all populations, the empty seed ratio was 46%. The study suggests that the Q Index and ME should be taken into account to determine cone to pine nut yield. With this regard, plantations of stone pine should be established on sub-humid areas and Q values should not exceed 75.


Background
Broad range of ecological conditions could create a variety of ecotypes (Huang 1993). This could contribute to the differences in characteristics of species such as cone, seeds, and shape which are important parameters related to improving seed yield. A better understanding of cone-yield variation patterns would facilitate many aspects of management (Mutke et al. 2005). Many studies have assessed cone and seed morphology of conifer populations regarding geographic and climatic variation (Şefik 1964;Kaya and Işık 1997;Mutke et al. 2005;Liu et al. 2013;Loewe-Muñoz et al. 2020). Climatic variables such as annual and autumn thermal oscillation, spring hydric index, and spring and summer precipitation are directly linked with cone production Loewe-Muñoz et al. 2016).
Stone pine is one of the commercially and ecologically important forest trees. Nuts constitute the most important non-wood product obtained from Mediterranean forests (Calama et al. 2016); it is one of the most characteristic species in the Mediterranean basin, with its umbrella shape and its use in horticulture, being cultivated for ornamental purposes. The cones are harvested for their edible kernels (pine nuts), which have been used as a food item since Paleolithic times (Mutke et al. 2005). Stone pine nuts are among the most expensive nuts; this high price makes them attractive as a crop. The species performs well in poor or eroded soils; its resistance to pests and diseases reduces the need for cultivation practices, and its high resistance to drought makes it a good candidate for conversion to a horticultural crop. Around the Mediterranean Sea, there are about 700,000 ha of stone pine-dominated forests, stretching from the Atlantic coast in Portugal to the shores of the Black Sea and Mount Lebanon (Mutke et al. 2012). Its distribution is mainly concentrated in Spain (450,000 ha), Portugal (90, 000 ha), Turkey (50,000 ha), and Italy (40,000 ha) (Pereira et al. 2015).
Stone pine is characterized by high phenotypic plasticity and adaptability, but low genetic variability (Loewe-Muñoz et al. 2016). Genetic variation is generally considered to be an important adaptation to environmental conditions. Genetically depauperate, but geographically widespread, species are rare, and no other abundant and widespread plant species has as little genetic diversity as the stone pine. However, the species does have a considerable amount of variation in its adaptive traits (Vendramin et al. 2008).
Besides, seed production of trees is key to regeneration of forests and to successful conservation management (Gonçalves and Pommerening 2012). Successful restoration and rehabilitation in Mediterranean areas greatly depends on the quality of the seeds and seedlings used (Boydak and Çalışkan 2014). Significant degradation of the vegetation across large areas of the Mediterranean Basin has called for rehabilitation using native coniferous and broadleaved species (Boydak et al. 2006). Stone pine has been used in a variety of remedial ways, such as for ecological restoration, watershed and soil protection, the stabilization of dunes, and afforestation in urban areas (Mutke et al. 2012;Çalışkan 2014, 2015;Çalışkan and Boydak 2017). The first dune afforestation project using Pinus pinea is thought to have been by the Romans along the Turkey-Antalya (Köprüçay) coastline, in an effort to stop mass movement and to prevent invasions of their agricultural land (Boydak and Çalışkan 2014).
On the other hand, Mutke et al. (2005) stated that further efforts will be required in order to accurately estimate the impacts of temperature rise and changes in rainfall patterns on cone development yield in stone pines. In the last few decades, efforts have been made in Spain, Turkey, and Portugal to domesticate the species, including intensive planting, reproductive propagation by grafting selected clones for cone production (Mutke, cited in Calama et al. 2016), and by employing irrigation and fertilization techniques Loewe-Muñoz et al. 2017). Despite this, the main crop of cones collected and marketed are from either natural forests or afforestations (Mutke et al. 2012).
There is limited literature concerning variations among native populations of the species, in terms of cone, seed, germination, and seedling traits, and their correlation with climatic variables (Baskin and Baskin 1998;Ganatsas et al. 2008). The aims of the present study were to determine (a) the variation in cone, seed, and seedling morphological traits between, and among, the eastern stone pine populations (specifically, in Turkey); and (b) the relationships among climatic and geographical variables with cone and seed traits. With this regard, the present study will contribute to the reforestation, yield, and breeding of this species.

Seed material collection
Natural stands of stone pine were selected from five populations across Turkey (Fig. 1, Table 1). Mature cones were collected in December 2016 from 15 to 25 randomly selected, mature trees in each population, which were growing at a distance of at least 50-100 m from each other, in order to capture a range of genetic variation (Ürgenç 1982). Ten cones were collected per tree, each placed in separate, labeled cloth bags (to prevent mildew). The cones were transported to the laboratory where they were stored at room condition until analyzed (Buru et al. 2016). The seeds were then manually extracted from 50 randomly selected cones from each of the five populations.
In Table 1, the coefficient 'S' explains the severity of the summer drought, with higher values indicating weak summer drought conditions. The coefficient 'Q' refers to climatic humidity and temperature, with lower values indicating more arid conditions. The coefficient 'm' represents the duration of the frost cycle, with higher values indicative of a shorter cycle (Dirik 2008).
We used Emberger's equation, which is commonly used for the Mediterranean region, to categorize the bioclimate (Daget et al. 1988), the pluviothermic quotient (Q) calculated as: where P is the mean annual rainfall (in mm), M is the mean of the maximum temperatures in the hottest month (in°C), and m is the mean of the minimum temperatures in the coldest month (in°C). The summer drought index (S) was calculated as S = PE/ME, where PE is the sum of the rainfall (in mm) in June, July, and August, and ME is mean of the maximum temperatures in June, July, and August.

Cone and seed measurements
Fifty randomly selected cones, collected from 16 to 25 families in each population, were completely dissected to  determine the seed yields of the five populations. The seeds extracted from each cone were separately immersed in water to distinguish the sound (those that sank) seeds (Table 2). All the seed measurements were performed on a random sample of 10 seeds per cone (Ganatsas et al. 2008). The 1000-seed weight (1000-SW) was determined using the ISTA (1999) procedure.

Seed germination
Germination tests were performed in a growth chamber (model VT 3 4034; Vötsch Industrietechnik, Germany), maintained at 15, 20, and 25°C until no additional seed germinated was observed. Randomly selected, sound seeds, with three replications (3 × 50) for each population, were used for germination. These were placed in 11-cm-diameter glass Petri dishes that were lined with two layers of filter paper, moistened with deionized water. The germinations were assessed daily, with germinated seeds being removed from the Petri dishes when their radicles had protruded (Farelli et al. 1997;Saracino et al. 2017). At the end of the test, the ungerminated seeds were subjected to a cut test to determine the proportion of sound to empty seeds in each replication of each treatment. Stone pine seeds were also sown in the Forest Nursery Directorate in Bahçeköy/Sarıyer/Istanbul, Turkey at the beginning of March of 2017. The nursery is located at an altitude of 126 m above sea level at 41°10′ 56′′ N, 28°59′ 14′′ E, and the climate around the nursery is humid, mesothermal, and maritime, with a moderate deficit of water in the summer months, according to Thornthwaite's classification. The mean annual precipitation is about 1111.4 mm, and the mean annual temperature is 12.8°C (Akburak et al. 2018).
Plastic containers (18 cm tall and 190 cm 3 ), filled with conventional media (soil, peat, and river sand in a ratio of 2:1:1) were used for sowing the seeds. For the trial, a randomized block design, with five replications, was used (135 seeds × 5 replications × 5 populations = 3375 sound seeds used). The germinations in the nursery were recorded two times in a week until no further germinated was observed. The germination percentage (GP) and mean germination time (MGT) for laboratory and field conditions (FGP and FMGT) were calculated thus: GP or FGP (%) = G/T × 100 where G is the total number of germinated seeds and T is the total number of seeds; and MGT or FMGT = ∑ (t × n) / ∑ n where t is the number of days from the beginning of the test and n is the number of germinated seeds on day t (Bewley and Black 1994).

Statistical analysis
All statistical analyses were carried out using SPSS (2010) software. The GP and FGP data were arcsine transformed. The GP, FGP, MGT, and FMGT data of the five populations were analyzed using ANOVA. The following ANOVA models were used: where ϒ i is the observed value of the cone, seed, and seedling traits; μ is the overall mean; φ i represents the effect of population; and ε i is the effect caused by error. The populations were excepted as fixed factors; and where ϒ ik is the observed value of the germination traits; μ, φ i , k , and ε ik are the overall mean, population effect, temperature effect, and error, respectively. The populations and temperatures were excepted as fixed factors. Duncan's post hoc test was used to determine homogeneous subgroups. Pearson's correlation coefficients between the cone, seed and seedling traits, germination values, and geographical and climatic variables were determined.

Results
The ANOVA results for the cone and seed traits showed significant differences among the five native populations from Turkey. The range of variation in the cone diameter (CD), cone length (CL), and cone weight (CW) ranged from 48-87 mm, 58-150 mm, and 82-435 g, respectively. The seed diameter (SD), seed length (SL), and 1000-SW varied between 6.5 and 10.9 mm, 11.8 and 22.4 mm, and 589 and 903 g, respectively. The pine nut yield (PY) was low in all the populations, ranging from 1.20 to 1.55%, except for the Muğla-Katrancı (3.53%) and Aydın-Koçarlı (4.47%) populations (Table 3).
The germinations in the growth chamber started at days 10, 12, and 15, at 15, 20, and 25°C, respectively, for the five stone pine seed lots. The effects of temperature were statistically significant for GP and MGT. Interaction between provenances and temperatures did not reveal any significant effect. The GP and MGT at 15°C in the growth chamber showed significant differences among the five populations. The highest GP (93%) and lowest MGT (12 days) were from the Aydın-Koçarlı population at 20°C. The optimum seed germination temperature was 20°C for all populations, except Mugla-Katrancı, which germinated 10% better at 15°C than at 20°C. The MGT values were the lowest for all populations at 20°C, and the mean MGT for the five populations was 13 days. The germinations in the field emerged at around 50 days (in May) after sowing (in March). The field germination percentage (FGP) and field mean germination time (FMGT) values were 50% and 69 days, respectively; they showed no significant differences within the populations (Fig. 3, Table 4).
The number of cotyledons (NC) varied between 8 and 16. The highest NC (16) belonged to the K. Maraş-Önsen population. This population also had the tallest seedling height in the first year (SdL1). In the second year, the seedlings from the İzmir-Kozak population were the tallest (38 cm), with the K. Maraş-Önsen being the second tallest (34 cm). The Muğla-Katrancı population had the highest root-collar diameter in the first (SdD1) and second year (SdD2) (3.3 and 6.1 mm, respectively). There was a relatively large variation in the number of lateral branches in the first (LB1) and second-year (LB2) seedlings. No terminal buds were observed in the populations in the first year. In the second year, no terminal buds were observed in the İzmir-Kozak population, while the other populations' bud sets (BS) ranged from 8 to 58%. In the K. Maraş-Önsen population, more than half the seedlings formed terminal buds (Table 5).
According to the results of the correlation analysis, there were highly significant correlations between cone size and seed size. CD and CL had a significant (the strongest) positive effect on CW (r: 0.885 and r: 0.849, respectively). SD and SL significantly affected the SW (r: 0.921 and r: 0.866, respectively). SL had a significant positive effect on CD and CW (r: 0.840 and r: 0.765, respectively). SN increased with increasing cone size. PY was negatively related to ES. In addition, PY had a significantly positive correlation with cone size, including CW and SL (Table 6).
The correlation analysis revealed that ME had the most significant effect on CD and SL (r: 0.939 and r: 0.889, respectively), and also had a positive effect (r: 0.866) on the 1000-SW, while summer drought index (S) had a negative effect on this (r: − 0.860). Although the sum of rainfall in June, July, and August (PE) and S factors had a negative effect on the 1000-SW, the mean of the maxima of the hottest month (M) had a positive effect (r: 0.723). Q had substantial negative effects on SN (r: − 0.897), CW (r: − 0.820), CL (r: − 0.823), and CD (r: − 0.25). S had substantial negative effects on CD (r: − 0.806) and SD (r: − 0.846). Q had a significant negative effect on SN, pine nut per cone (PN), and PY (r: − 0.897, r: − 0.923, and r: − 0.903, respectively). In addition, Q had a negative correlation with all cone, seed, and pine nut parameters, except for ES, while ME positively correlated with all the parameters, except for ES. RJ had a negative influence on ES (r: 0.846), while AR is significantly related to pine nut weight (PW) (r: 0.889) ( Table 7).
The SdL1 affected the SdL2 (r: 0.760). LB2 had a substantial negative effect on the SdL1 and SdL2 (r: − 0.915 and r: − 0.808, respectively). The correlation coefficient between the number of cotyledons (NC) and the BS at the end of second year was r: 0.876. The longitudes of the five native stone pine populations positively correlated with BS and NC (r: 0.976 and r: 0.892, respectively). The increasing altitude of the populations resulted in increased SdD1 and NC (r: 0.713 and r: 0.758, respectively) ( Table 8).

Cone traits
The mean CW of 251 g is consistent with reported ranges of 227-304 g, based on 3-year surveys (Gonçalves and Pommerening 2012) In the present study, the average CL (9.7 cm) and CD (7.5 cm) are smaller than those reported by Loewe-Muñoz et al. (2018), which were 11.3 cm (CL) and 9.4 cm (CD), but were bigger than those from southern Greece (CL 8.3 cm and CD 6.8 cm) (Ganatsas et al. 2008).
Stone pine cones have a 3-year reproductive development cycle, which begins when the female strobili (conelets) burst in the spring of the first year, with pollination occurring a few weeks after that. Strobilus fertilization, however, takes place only in the spring/summer of the third and final year of cone maturation, when the nutritious, edible seeds are formed (Valdiviesso et al. 2017). It is expected that a reduction in rainfall will have a negative effect on cone production.  indicated that spring and summer rainfall has a positive impact on cone production, meaning cone number per tree, CW, and PY. A decrease in stone pine seed production could also result from heavy rains during pollination (Parlak et al. 2013). In our study, rainfall in the June of the third year prior to the cone-collecting year (RJ) (Mutke et al. 2005), during pollination, is a climatic factor that might cause a decrease in the sound seed ratio. This supports the findings of Parlak et al. (2013), who highlighted the importance of the pollination period on cone production. They also pointed out  that non-productive areas located at altitudes below 500 m were identifiable by a thermal oscillation of 22°C and a higher frequency of days with temperatures below 10°C , early and late frosts, high relative humidity, and fog. Annual rainfall in the third year prior to harvest time (AR) had a positive impact on CW, as also reported by Loewe-Muñoz et al. (2016). AR showing inter-annual variation has also been positively related to cone yield (Loewe-Muñoz et al. 2020).
The climatic factors Q and ME were found to be the main drivers behind the stone pine cone, seed, and pine nut morphometric traits. The most significant limiting factors to annual cone yield appeared to be water stress during different stages of cone development and the negative effect of hot midsummers on cone setting (Mutke et al. 2005). Loewe-Muñoz et al. (2016) suggested that the selection of sites for fruiting should ensure minimum average spring temperatures above 7°C, annual and autumn thermal oscillations below 12°C, and a high spring rainfall, except during the male flowering period. In the final year, losses in seed yield (with up to 50% of the seeds being empty) have been reported from Portugal, Spain, Italy, and Turkey (Mutke et al. 2014), but in Chilean plantations, the percentage of damaged/empty seeds was relatively very low, at around 9% (Loewe-Muñoz et al. 2019). We found that the overall sound seed ratio was 54%, ranging from 32 to 80%. This ratio could change in the future because of stand  properties. The number of sound seeds is an important factor in afforestation purposes. Seed production is a limiting factor for regeneration, and it is necessary to improve seedling survival in stone pine stands .

Seed traits
We found that the average SW (0.8 g) differed among the populations, and was higher than the 0.6 g determined by Mutke et al. (2011) and the 0.7 g of Ganatsas et al. (2008), while it was lower than the 0.9 g found by Loewe-Muñoz et al. (2018). The average SD and SL determined by Loewe-Muñoz et al. (2018) were 1 and 2 cm, respectively, while in the present study, they were 0.9 and 1.8 cm, respectively; these latter measurements are similar to those reported by Ganatsas et al. (2008) and Carcel et al. (2012). Loewe-Muñoz et al. (2019) report that the mean 6year measurement of SN in plantations in Chile was 113 which is more than double that found in the present study (50), with the Çanakkale-Kirazlı population having the lowest (25). A positive correlation was found between SN and CW; this finding supports several previous studies that have indicated that heavier cones contain more seeds Boutheina et al. 2013;Loewe-Muñoz et al. 2019). In our findings, the average cone to seed yield (SY) ranged between 8.1 and 18.5%, while Loewe-Muñoz et al. (2019) reported 6-year-mean values ranging from 17.7 to 22.6%.
The 6-year measurement of PY varied from 3.6 to 5.0% in plantations in Chile (Loewe-Muñoz et al. 2019), which is approximately 1.5 times higher than was measured in our populations. We also observed that populations with larger and heavier cones had the highest PW, consistent with the findings of Loewe-Muñoz et al. (2019), who highlighted that PY depended first on PN, then on PW.
Loewe-Muñoz et al. (2019) stated that cone morphometry is not a good indicator of PY, due to the inverse relationship between CW, CD, CL, and PY, so selecting cones for size/weight would not improve PY. However, our results did not support this idea--increased CW, Table 5 Seedling traits of the five natural stone pine (P. pinea) populations in Turkey

Parameters
Populations Mean 1 K.Maraş-Önsen 2 Muğla-Katrancı 3 Aydın-Koçarlı 4 İzmir-Kozak 5 Çanakkale-Kirazlı NC 13 (9-16) c 11 (9-13) ab 12 (9-13) b 11 (8-13) a 11 (9-13) a 12 (9-14) Means within columns following by same letters are not statistically different at 0.05 significance level in Duncan post hoc test. For abbreviations, please see Table  2. Minimum and maximum values are presented in parenthesis * and ** stating significant differences at 0.05 and 0.01 levels, respectively. For abbreviations, please see Table 2 CD, and CL resulted in increased PY. The PY is particularly interesting because of an increase of empty and damaged seeds (Loewe-Muñoz et al. 2019), which have been associated with increased drought or phenological shifts due to climatic change (Mutke et al. 2014) and infestation with Leptoglossus occidentalis Heidemann (Farinha et al. 2018). Many studies have found that the seed size in many species is nonlinearly related to altitude and latitude, being dependent on low temperatures and radiation and a shorter growing season at higher latitudes or altitudes (Moles and Westoby 2003;Murray et al. 2004;Moles et al. 2007). Şefik (1964) determined that the SW of Pinus brutia (Turkish red pine) decreases from low to high altitudes, and from south to north. Also, Kaya and Işık (1997) found that the SW in P. brutia decreases as it ascends above sea level. Some seeds, such as those of Picea mariana and Pinus banksiana, negatively correlate with latitude and altitude, while those of Picea glauca have shown no correlations with latitude or altitude in  the Boreal Shield Ecozone. In the same study, in the Atlantic Maritime Ecozone, seed size in Picea mariana and Picea glauca had no relationship with latitude or altitude, while seed size in P. banksiana did show positive correlations with latitude and altitude (Liu et al. 2013).
In the present study, we have confirmed that this nonlinear relationship is not consistent within populations of stone pine, which indicate an increase in seed and cone size with increasing altitude and decreasing latitude. Furthermore, in the present study, longitude had no influence on cone or seed size. We can conclude that mean annual rainfall was not a determinant of cone or seed size in stone pine populations.
Considering that the SW decreases as elevation increases, the altitude factor should be taken into consideration for stone pine plantations being used for seed production. It should be noted that, where the soil properties are suitable, plantations above 400 m might be more productive. The impact of drought is likely to be lessened in such locations, compared to areas close to sea level, and the nutritional environment might be more suitable.
Seed germination and seedling traits Ganatsas and Tsakaldimi (2007) found 87.2% to 100% stone pine seed germination, with a MGT of 13.4-15.6 days at 20°C. According to Escudero et al. (2002), GP in stone pine decreased at 25°C, and differences in the GPs among populations was more than 50%. Population effects have been found to be significant in P. pinea and P. pinaster, and in their relationship to temperature. In addition, significant interaction has been reported among population, temperature, and germination rate in P. pinea. However, it has also been pointed out that, unlike mountain pines (Pinus nigra, Pinus sylvestris), lowland pines (Pinus halepensis, P. pinaster, P. pinea) and Pinus canariensis have exhibited significant effects of temperature on germination response, and germination is higher between 15 and 20°C than at warmer and alternating temperatures.
It has also been reported that successful seedling establishment depends on rapid germination and growth, such as in Pinus strobus (Parker et al. 2006) and Pinus densata . As the germination speed of a seed lot increases, the height of the seedlings from the seed lot may also increase. In other words, the lower the MGT, the taller the seedling height (Çalışkan 2006). The present study concluded that there was no significant intraspecific variation in the germinations between the field and laboratory (20°C), although all populations showed variety in seed size, with bigger seeds germinating faster (Moles and Westoby 2004). Our results support certain studies in which it was indicated that seed weight had little effect on germination in Picea sitchensis (Chaisurisri et al. 1992) and P. pinaster (Wahid and Bounoua 2013). Additionally, bigger seeds are prone to producing bigger seedlings in Pinus yunnanensis ).
In the present study, NC exhibited significant variety across the native stone pine habitats, while Ganatsas et al. (2008) reported that this was lower in over-mature stands. Moreover, it was determined that NC positively correlated with the longitude of the origins of the stone pines. Also, there was a strong correlation between NC and BS. It could be that BS relates better to mean annual rainfall in our study, while Meng et al. (2015) indicated that annual mean temperature and mean temperature of the coldest month have a stronger relationship with BS in P. yunnanensis over 2 years.
A limitation of this study was that it depended on 1year results, whereas inter-annual variations in seed production have been documented by several authors (Mutke et al. 2005;Calama et al. 2008Calama et al. , 2011Calama et al. , 2016Loewe-Muñoz et al. 2019, 2020. Such variations could be the result of climatic factors affecting the 3-year reproductive development cycle of cones in stone pine.

Conclusions
The findings of this study contribute to a better understanding of native populations (with an eastern distribution) of stone pine, and should be taken into consideration in forest management especially relating to afforestation, yield, and breeding activities. Stone pine plays an important role in the rehabilitation and restoration of the Mediterranean Basin, the success of which depends on the quality of seeds and seedlings. Considering intraspecific variation, Çanakkale-Kirazlı population should be assessed within a framework of in-situ conservation due to having distinctive value in all cone and seed morphometric traits among the studied populations. The overall empty seed ratio was 46%, ranging from 20 to 68% among the populations. This dramatic decrease in seed production has been reported in previous works, and is understood to be associated with increasing drought, or phenological shifts due to climatic change and infestation with Leptoglossus occidentalis Heidemann. With seed weight decreasing with elevation, the altitude factor should be taken into consideration in stone pine plantations for seed production. It should be noted that, where the soil properties are suitable, plantations above 400 m may be more appropriate. In terms of climatic factors, Q and ME were found to be the main drivers behind significant effects in stone pine cone, seed, and pine nut morphometric traits. We suggest that the Q Index should be used to determine cone to pine nut yield. Moreover, stone pine plantations should be focused on sub-humid areas with Q values of less than 75. Further studies are needed to interpret these results in relation to the stand properties, and long-term research is needed to investigate how the stone pine will react to the novel environmental conditions of the future.