Next Article in Journal
Optimizing Antioxidant Activity and Phytochemical Properties of Peppermint (Mentha piperita L.) by Integrative Application of Biofertilizer and Stress-Modulating Nanoparticles under Drought Stress Conditions
Next Article in Special Issue
Variation in Seed Dormancy of Chaco Seasonally Dry Forest Species: Effects of Seed Traits and Population Environmental Conditions
Previous Article in Journal
Wheat Crop under Waterlogging: Potential Soil and Plant Effects
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Potential Distribution of Cedrela odorata L. in Mexico according to Its Optimal Thermal Range for Seed Germination under Different Climate Change Scenarios

by
Salvador Sampayo-Maldonado
1,
Cesar A. Ordoñez-Salanueva
1,
Efisio Mattana
2,
Michael Way
2,
Elena Castillo-Lorenzo
2,
Patricia D. Dávila-Aranda
3,
Rafael Lira-Saade
3,
Oswaldo Téllez-Valdés
3,
Norma I. Rodríguez-Arévalo
3,
Cesar M. Flores-Ortiz
1,4,* and
Tiziana Ulian
2
1
Plant Physiology Laboratory, Biotechnology and Prototypes Unit (UBIPRO), FES Iztacala, Universidad Nacional Autónoma de Mexico (UNAM), Tlalnepantla 54090, Estado de Mexico, Mexico
2
Royal Botanic Gardens, Kew, Wakehurst, Ardingly, Haywards Heath, West Sussex RH17 6TN, UK
3
Natural Resources Laboratory, Biotechnology and Prototypes Unit (UBIPRO), FES Iztacala, Universidad Nacional Autónoma de Mexico (UNAM), Tlalnepantla 54090, Estado de Mexico, Mexico
4
National Laboratory in Health, FES Iztacala, Universidad Nacional Autónoma de Mexico (UNAM), Tlalnepantla 54090, Estado de Mexico, Mexico
*
Author to whom correspondence should be addressed.
Plants 2023, 12(1), 150; https://doi.org/10.3390/plants12010150
Submission received: 6 December 2022 / Revised: 23 December 2022 / Accepted: 24 December 2022 / Published: 28 December 2022
(This article belongs to the Special Issue Advances in Seed Physiology)

Abstract

:
Cedrela odorata is a native tree of economic importance, as its wood is highly demanded in the international market. In this work, the current and future distributions of C. odorata in Mexico under climate change scenarios were analyzed according to their optimal temperature ranges for seed germination. For the present distribution, 256 localities of the species’ presence were obtained from the Global Biodiversity Information Facility (GBIF) database and modelled with MaxEnt. For the potential distribution, the National Center for Atmospheric Research model (CCSM4) was used under conservative and drastic scenarios (RCP2.6 and RCP8.5 Watts/m2, respectively) for the intermediate future (2050) and far future (2070). Potential distribution models were built from occurrence data within the optimum germination temperature range of the species. The potential distribution expanded by 5 and 7.8% in the intermediate and far future, respectively, compared with the current distribution. With the increase in temperature, adequate environmental conditions for the species distribution should be met in the central Mexican state of Guanajuato. The states of Chihuahua, Mexico, Morelos, Guerrero, and Durango presented a negative trend in potential distribution. Additionally, in the far future, the state of Chihuahua it is likely to not have adequate conditions for the presence of the species. For the prediction of the models, the precipitation variable during the driest month presented the greatest contribution. When the humidity is not limiting, the thermal climatic variables are the most important ones. Models based on its thermal niche for seed germination allowed for the identification of areas where temperature will positively affect seed germination, which will help maximize the establishment of plant populations and adaptation to different climate change scenarios.

1. Introduction

Climate change is going to alter the patterns of rainfall and temperature in various regions of the country. As a result, a higher frequency of extreme meteorological events will affect subtropical and tropical forests, affecting species phenology (e.g., flowering and germination time) [1,2,3,4,5,6,7], and modifying the structure and composition of the tree communities, which will alter migration and its geographical distribution [8,9,10,11,12,13,14,15,16]. Among the factors that shape species distribution, the dispersal of its seeds and the environmental conditions suitable for germination play important roles [17,18,19]. Therefore, it is necessary to consider species optimum germination temperature range in order to identify the new areas in which it will be possible to find the species as well as where the environmental conditions will not allow its germination and establishment in the future.
Cedrela odorata L. provides several ecosystem services, among which CO2 capture and water balance stand out [20]. The wood of C. odorata L., due to its color and texture, is sold as sawn wood and plywood for the manufacture of furniture, for which there is a high demand in the national and international markets. Considering the C. odorata wood reserves in natural forests, according to the information generated by the National Forest Inventory, Romo-Lozano et al. [20] estimated that the species represents a financial value of USD 288,328,788 for Mexico. For these reasons, it is a species of major importance in the Mexican forestry industry [21]. However, the intensive extraction of this wood has decreased its populations and areas of distribution. Therefore, the species has been overexploited and threats to this species are exacerbated by its vulnerability to pests [7].
Due to its vulnerable conservation status in natural forests, C. odorata is included in Appendix II of the Convention on International Trade in Endangered Species [22]. In addition, it is listed as an endangered species subject to special protection in the Official Mexican Standard 059 (NOM-059-SEMARNAT) [23]. Finally, climate change has already been reported to alter its phenology, interactions, distribution, morphology, and net primary productivity in tropical forests [24].
As an effect of climate change, an increase in temperatures and a decrease in precipitation are forecasted for Mexico [7]. Thus, studies on potential species distribution are needed in order to understand how a warming climate will impact the development of a single species [25]. Several studies have already been carried out to model and predict the future distribution of C. odorata. Hernández-Ramos et al. [26] generated an ecological niche model using 19 climate and terrain variables as well as altitude, orientation, slope, and soil type layers, which suggested a reduction in the ecological niche of the species in the country. In addition, Gómez-Díaz et al. [25] estimated the potential climate areas of the distribution of C. odorata in the state of Hidalgo, and they obtained conflicting results according to the two applied climate change scenarios. For the American model (GFDL-R30) [27], an expansion in the surface of 3.1 and 4.4% is expected for the years 2020 and 2050, respectively, while for the English model (HadCM3) [27], a contraction in the surface of 0.9% for 2020 and 0.2% for 2050 is expected.
Species distribution models (SDM) have been used to analyze the impact of climate change on the distribution of species, among which the Maximum Entropy model (MaxEnt) stands out for its precision [28]. It feeds on species presence records as well as environmental and climatic variables. However, when physiological variables are incorporated into the model, it becomes stronger and increases in precision. The model performs a maximum likelihood estimation of the relative probability of presence and analyzes the relationship between species distribution dynamics and climate change [17,19,28].
For the study of plant distribution, bioclimatic indices are usually used, which are based on temperature and precipitation data as climatic variables. These indices define the physiological limitations of the species, where the minimum temperature is a factor that can limit the growth and development of plants, especially for those from warm regions such as the species under study [29]. Funes et al. [30] identified temperature as the main regulating factor that triggers germination and maximizes seedling establishment, survival, and regeneration in the understory of forests. Durán-Puga et al. [31] mentioned that faster germination rate confers a competitive advantage in colonizing fragmented habitats. Temperature is an important element of the climate, suitable for analyzing the response of seed germination under different climate change scenarios [32].
González [33] stated that temperature influences the physiological processes of seed germination since it determines both the speed and maximum percent of germination. García-Huidobro et al. [34] concluded that the maximum and minimum temperatures determine the limits within which germination is possible, while the optimum temperature, in agreement with Gilbertson et al. [35], is the temperature at which germination is fastest and seedlings are acquired in the amount least time.
Calzada-López et al. [36] support the use of the cardinal temperatures—in order to find the optimum temperature where the highest amount of germination occurs in the shortest time—as a criterion for the distribution of the species as a direct effect of adaptation to different climate change scenarios [37]. C. odorata seeds are reported to be non-dormant [7,38]. Seeds of this species have a wide range of temperatures for germination, ranging from −0.5 ± 0.09 to 53.3 ± 2.1 °C, while the thermal time was 132.74 ± 2.60 °C, which accumulates in a shorter period due to climate change, thus accelerating the germination of seeds in the soil. [7]. It is important to mention that there are no previous studies that take into account optimum germination temperature range in the potential distribution of this species.
According to Adam et al. [39], cardinal temperatures for germination tend to be similar to those of normal vegetative growth of plants, although, in some cases, it may need higher temperatures for flowering. Thus, the environmental conditions where the species is distributed (suitable climate habitat), especially where the temperature is in the range of temperature thresholds for germination, are of paramount importance. In this sense, to consider the optimum temperature range when modelling the potential distribution of C. odorata lays the basis for locating areas where the temperature will have a positive effect on seed germination, thus establishing the species in changing climate conditions [40]. The objective of this study was to model the current and future distributions of C. odorata according to its thermal niche for seed germination under different climate change scenarios using environmental variables with two representative concentration pathways (RCP) of CO2 emissions.

2. Results

2.1. Current Distribution of C. odorata

Using the georeferenced data from the Global Biodiversity Information Facility (GBIF) and the BioClim database from the period of 1970 to 2000, a map of the current distribution of optimal climate habitat was generated for C. odorata. The model was built from species presence data that were within the optimum germination temperature range. The current surface of C. odorata with optimal climatic habitat for germination was 228,780.66 km2. According to the MaxEnt model, it is distributed from 10 to 28.7 °C annual average temperature, which falls into the suboptimal range for seed germination (Figure 1). According to the jackknife, the variables that contributed most to the model were Bio14 (31.4%), Bio12 (29.6%), and Bio6 (23%).

2.2. Potential Distribution of C. odorata

According to the NCAR model (CCSM4) for the intermediate future (2050) under the conservative scenario (RCP2.6 Watts/m2), it is predicted that the potential distribution will be of 243,890.85 km2 with optimal climate habitat (Figure 2 and Table 1). This represents a 6.6% expansion compared with the current distribution. According to the jackknife, the variables that contributed most to the model were Bio14 (36.8%), Bio6 (28.3%), Bio12 (11.6%), and Bio7 (9.3%).
The states in which the distribution surface with optimal temperatures for germination will expand compared with the current model are Michoacan (219%), Hidalgo (179%), Queretaro (92%), Colima (70%), Jalisco (70%), Sinaloa (50%), Oaxaca (27%), and Veracruz (26%). In contrast, the states whose surface with temperatures for germination of C. odorata will contract are Chihuahua (93%), Morelos (84%), Mexico (84%), Guerrero (64%), Puebla (36%), Sonora (30%), San Luis Potosi (28%), Durango (25%), Tamaulipas (18%), and Quintana Roo (14%), compared with the current model (Table 1). The state of Guanajuato, which does not have C. odorata populations in the current scenario, in the conservative scenario will present 15.38 km2, which will be located in the northeast part of the state in the municipalities of Xichu and Victoria.
For the drastic scenario (RCP8.5 Watts/m2), it is predicted that the potential distribution will be of 241,533.21 km2 with optimal climate habitat (Figure 2 and Table 1). This represents a 5.57% expansion compared with the current distribution. However, it represents a 1.03% contraction compared with the conservative scenario in the intermediate future. According to the jackknife, the variables that contributed most to the model were Bio14 (43.3%), Bio6 (28.7%), Bio12 (9.7%), and soils (4.3%).
The states in which the distribution surface with optimal temperatures for germination will increase compared with the current model are Hidalgo (134%), Sonora (98%), Jalisco (75%), Sinaloa (47%), Queretaro (34%), Veracruz (28%), Oaxaca (27%), and Chiapas (22%). In the opposite case, the states whose surface with temperatures for germination of C. odorata will decrease are Chihuahua (99%), Morelos (91%), Mexico (81%), Guerrero (62%), Durango (48%), Colima (37%), San Luis Potosí (25%), Tamaulipas (25%), and Quintana Roo (18%) compared with the current model (Table 1). In the intermediate future, Guanajuato will present 6.64 km2, which will be located in the northeast part of the state in the municipalities of Xichu and Victoria.
For the far future (2070) in the NCAR model (CCSM4) under the conservative scenario (RCP2.6 Watts/m2), it is predicted that the potential distribution will be of 239,938.67 km2 with optimal climate habitat (Figure 3 and Table 1). This represents a 4.87% expansion compared with the current distribution. According to the jackknife, the variables that contributed most to the model were Bio14 (37.9%), Bio6 (24.5%), Bio12 (18.5%), and Bio7 (5.2%).
The states in which the distribution surface with optimal temperatures for germination will increase compared with the current model are Michoacan (199%), Hidalgo (173%), Queretaro (120%), Jalisco (81%), Sinaloa (48%), Colima (41%), Tabasco (29%), Veracruz (28%), and Chiapas (22%). In the opposite case, the states whose surface with temperatures for germination of C. odorata will decrease are Chihuahua (96%), Morelos (92%), Mexico (89%), Sonora (77%), Guerrero (64%), San Luis Potosi (35%), Durango (26%), Quintana Roo (22%), and Tamaulipas (18%) compared with the current model (Table 1). In the intermediate future, Guanajuato will present 65.1 km2, which will be located in the northeast part of the state in the municipalities of Xichu and Victoria.
For the drastic scenario (RCP8.5 Watts/m2), it is predicted that the potential distribution will be of 246,681.45 km2 with optimal climate habitat (Figure 3). This represents a 7.82% expansion compared with the current distribution, representing a 2.81% expansion compared with the conservative scenario for the far future. According to the jackknife, the variables that contributed most to the model were Bio14 (43.8%), Bio6 (25.8%), Bio7 (11.8%), and soils (4.3%).
The states in which the distribution surface with optimal temperatures for germination will expand compared with the current model are Colima (177%), Hidalgo (162%), Sonora (94%), Michoacan (83%), Jalisco (78%), Sinaloa (49%), Puebla (38%), Veracruz (31%), Chiapas (28%), and Oaxaca (19%). In the opposite case, the states whose surface with temperatures for germination of C. odorata will decrease are Chihuahua (100%), Mexico (98%), Morelos (91%), Guerrero (55%), Durango (49%), San Luis Potosi (24%), and Tamaulipas (23%) compared with the current model (Table 1). The state of Guanajuato in the intermediate future, will present 2.7 km2, which will be located in the northeast part of the state in the municipalities of Xichu and Victoria.
From the current to the potential distribution of C. odorata, the states of Chihuahua, Mexico, Morelos, Guerrero, and Durango presented a negative trend in land surface with optimal temperatures for germination. The most extreme case was the state of Chihuahua, which, in the current scenario, shows 433.45 km2, but in the intermediate future, it loses more than 93%, and for the far future, no surface is predicted. In the opposite case, the states of Michoacan and Hidalgo will present positive trends in surface with optimal temperature for germination. On the other hand, the models generated predict the distribution in optimal germination temperature range for the state of Guanajuato, which is indicative of the impact of climate change in that region. According to the above, the state of Guanajuato will present favorable conditions for C. odorata. Therefore, climate change can have a positive impact on the state’s forest production, but at the same time, it shows the lack of policies and programs to mitigate the negative impacts on the country’s fragile ecosystems.
According to the jackknife method, the variables that alone explain, to a greater extent, the models of the distribution of C. odorata were precipitation during the driest period, minimum temperature of the coldest month, and the annual temperature oscillation variable. The soil type variable has information that is not contained in any other variable; its presence is relevant since eliminating it would result in a decrease in the precision for the models. The generated prediction models presented area under the curve (AUC) values of 0.974 ± 0.007, which indicates a good adjustment of the models to the prediction of the areas most likely to find the species based on its temperature thresholds. For each scenario, more than 1500 candidate models were generated; the selected model was the one with an omission rate of 0.05%, a receiver operating characteristic (ROC) of 1.3, and a model complexity (AIC) of 2907.

3. Discussion

3.1. Current Distribution of C. odorata

The modeling of current and future distributions was based on the statistical principle of maximum entropy, according to Philips et al. [41], allowing predictions considering only records of species presence. For modeling the optimal distribution on the basis of its temperature thresholds, the environmental variables obtained from BioClim were used. The climate variables of BioClim have been used previously in the modeling of potential optimal distribution of forest species from Mexico, including Pinus chiapensis (Martínez) Andresen [42], Abies religiosa (Kunth) Schltdl and Cham [43], Pinus leiophylla Schiede ex Schltdl and Cham [44], Swietenia macrophylla King [45,46], and Cedrela odorata L. [26]. The current distribution model of C. odorata coincides with the distribution of the humid tropical forest biome and part of the seasonally dry tropical forest in Chihuahua, Sonora, Durango, and Sinaloa [47]. They present annual precipitations of more than 1200 mm, with maximum temperatures of 34 to 40 °C during spring and summer (end of March to end of August) and without frost.
The current distribution models obtained in this work coincide with those reported by Hernández-Ramos et al. [26]. The variables that contributed highly to the prediction of the model of current distribution of C. odorata by cardinal temperatures were precipitation during the driest period, annual precipitation, and minimum temperature of the coldest month, with 84%. The variables for the ecological niche modeling in this species are in agreement with Hernández-Ramos et al. [26], and they were annual precipitation, precipitation during the driest period, diurnal temperature oscillation, and altitude, with an 88.1% contribution. At the same time, the results are comparable to the potential distribution in Mexico of the optimal climate habitat for S. macrophylla (classified in the same family, Meliaceae), in which the variables with the largest contribution were precipitation during the driest month, minimum temperature of the coldest month, temperature seasonality, and annual temperature range; with a contribution of 78.1% [46], while in another study of S. macrophylla [45], the variables with the largest contribution were mean temperature during the coldest month, mean temperature, minimum temperature, and precipitation during the growing season, with a contribution of 66.7%. According to Harmann et al. [48], the distribution of forest species is determined mainly by temperature. A recent investigation [49] confirms that the temperature thresholds on the coldest days is the most important environmental variable to determine the distribution of the types of vegetation that exist in the world. This observation confirms the main patterns observed in this and previous studies [29] and highlights the importance of temperature in predicting the distribution of tree species in Mexico.

3.2. Potential Distribution of C. odorata

The variables that contributed the most to the models’ prediction of the potential distribution of C. odorata, based on its temperature thresholds for seed germination, were precipitation during the driest period, minimum temperature of the coldest month, annual precipitation, and annual temperature oscillation, with 86% in the conservative scenario for the future 2050 and 2070. However, in the model of the future distribution of the species proposed by Hernández-Ramos et al. [26] for 2050, the mean temperature of the warmest quarter, altitude, minimum temperature of the coldest month, and temperature seasonality contributed, with a 92% in the conservative scenario. In turn, Sampayo-Maldonado et al. [46] built a model to predict the distribution of S. macrophylla in the conservative scenario for the future 2050 and 2070. The variables that contributed the most were precipitation during the driest month, minimum temperature of the coldest month, temperature seasonality, annual temperature range, and altitude, with 85.5%. In another model for the species [45], for the 2030 in the intermediate scenario (RCP6), the variables that contributed the most were average temperature of the coldest month, minimum average temperature, and growing season precipitation, with 66.7%.
For the drastic scenario in the intermediate future (2050) in this work, the variables precipitation during the driest period, minimum temperature of the coldest month, annual precipitation, and soil type contributed, with 86%, while for the far future (2070), they contributed with 85.7%, which is very similar for the distribution of S. macrophylla [46]. Since the thermal variables are the same, they contributed from 80.4 to 83.2% for the construction of the distribution models of this species. Although the annual precipitation variable did not contribute, the annual temperature does an oscillation variable. Therefore, the temperature is a variable that should always be considered [29,49]. According to González [33], if the humidity is not limiting, the temperature is the most determining factor for the germination and distribution of the species. Each species has its own place in space and time that can be represented on a map [50].
According to the NCAR model (CCSM4) with a probability greater than 20%, under the intermediate future (2050) in the conservative and drastic scenarios, the potential distribution will expand 6.6 and 5.57%, respectively, compared with the current distribution. On the other hand, for the far future (2070), the potential distribution will expand 4.87 and 7.82% for the conservative and drastic scenarios, respectively, compared with the current distribution. This is without considering external factors, such as anthropogenic activities that lead to changes in land use [51]. Using the same model (CCSM4), Hernández-Ramos et al. [26] reported a contraction in the distribution of the species by 2050, with a model probability greater than 50%; however, they did not mention the percent of contraction. Another study used the same model (CCSM4) for the potential distribution of S. macrophylla [46], and they reported for 2050 an expansion of 18.3% in the conservative scenario and 14.4% in the drastic scenario. For 2070, they reported expansions of 17% and 12.3%, respectively, all with respect to the current distribution, with a probability greater than 20% in the construction of the models. Other authors [45] reported a 60% contraction in the distribution of S. macrophylla in the Yucatan Peninsula by 2030. They used a model based on averaging 18 climate models with intermediate greenhouse gas concentrations (RCP6.0 Watts/m2), which is a moderate scenario, although the probability for the construction of the model is not mentioned.
According to the species’ distribution models, an expansion in the area is expected in the future with respect to the current distribution, likely related to the rise in minimum temperature. However, when comparing the scenarios in the intermediate future (2050), in the drastic scenario (RCP8.5 Watts/m2), there was a contraction of 1.03% compared with the conservative scenario (RCP2.6 Watts/m2), which can be explained by the 1.9 °C temperature increase between the scenarios. For the far future (2070), in the drastic scenario (RCP8.5 Watts/m2), there was a 2.95% expansion compared with the conservative scenario (RCP2.6 Watts/m2), which can be explained by the 2.0 °C temperature increase between the scenarios. Therefore, according to the above, the increase in temperature benefits the distribution of the species, but this is valid up to a specific temperature in the near future. Parmoon et al. [52] calls it the optimum temperature, since a small increase affects the germination of the species, contracts the area of distribution, and reduces the ecological niche. According to Gómez-Días et al. [25], an increase above the optimum temperature affects the hormonal mechanisms involved in flowering, fruiting, and germination of the species, which is reflected in the decrease of surface suitable for germination and/or growth. Changes in distribution and phenology are expected with increasing temperature [7]. Temperature affects the hormonal mechanisms involved in the flowering and fruit set, which impacts the dispersal and distribution of species [46].
Particularly, in relation to flowering, it is known that this is determined by temperature and photoperiod [2]. It is a stage that is sensitive to photothermal changes; therefore, a change in the average temperature can modify the synchrony with pollinators and rainy periods. In this sense, climate change could promote earlier flowering, as temperatures continue to rise [3]. This will impact seed production and all associated phenological processes [1]. Additionally, the production, dispersal, and germination of seeds not only respond to instantaneous thermal signals, but also incorporate information from the thermal history of the plant’s ancestors (thermal memory) [19].
Since seeds are the dispersal unit of plants, germination determines the regeneration site, which affects the assembly of species in plant communities. Therefore, thermal niche breadth of germination is related to species distribution across different ecosystems, delimiting the regeneration site [53]. According to Fernández-Pascual et al. [19], models of potential distribution should consider how the future climate will affect germination percentage and velocity, as these two traits are important for the individuals’ dispersibility of populations as well as the vigor of their seedlings.
In general, there was an expansion of the potential distribution in most of the states studied for both the conservative and drastic scenario for 2050 and 2070. These results are in agreement with those found in Hernández-Ramos et al. [26], showing that C. odorata has a high probability of being distributed in the southern Yucatan Peninsula, north and south of Chiapas, and in the Gulf coastal plain, which is reflected in our study. On the other hand, the states of Chihuahua, Mexico, Morelos, Guerrero, and Durango suffer a contraction of the potential distribution for both the conservative and drastic scenarios. The most extreme case will be the state of Chihuahua, which will lose environmental conditions for the distribution of the species in the distant future. Therefore, it is recommended that the populations of the species in the state are monitored and the collection of seeds to conserve local genotypes is carried out. Similarly, other studies predict a 0.2% decrease in the potential distribution for the state of Hidalgo under the English model (HadCM3 phase 3 of the IPCC) by 2050 [25]. At the same time, Garza-López et al. [45] predicted a decrease in the distribution of S. macrophylla in Campeche, Chiapas, and Yucatán as well as complete disappearance in Quintana Roo by 2030.
On the other hand, the models generated indicate that the state of Guanajuato in the future will present favorable climatic conditions for the species, of which there is no record in previous studies. Therefore, it is recommended that provenance genetic tests are carried out and the genotypes that show the greatest adaptation are selected. With the increase in temperature, Michoacan and Hidalgo’s states will present a positive trend in surface for the distribution of species. This can be explained by the temperature rise from 1.9 to 4.7 °C, which the NCAR model (CCSM4) predicts in the far and intermediate future, which will favor the ideal conditions for the potential distribution of C. odorata. This coincides with Ortiz-Solorio [54], who mentioned that many physiological processes of trees are influenced by temperature. Therefore, increases in this bioclimatic element will accelerate the hormonal mechanisms of flowering and fruiting, which favors seed dispersal and impacts the distribution of the species.
According to the American model (CCSM4), an increase in temperature, a decrease in rainfall, and marked drought events are expected, altering the climatic conditions of forest regions, making them warmer and drier. Therefore, changes are expected in the five most important biomes for Mexico [47]. Tropical rain forests, mountain rain forests, and temperate forests are expected to contract dramatically. In contrast, the seasonal dry tropical forest and xerophilous scrub will expand its distribution considerably. Contrary to what would be expected, these climatic changes will favor C. odorata, according to the distribution model. Based on the range of optimal germination temperatures, conditions with suitable climatic habitat for the species are expected until 2070, which coincides with that which was reported for another tropical species S. macrophylla [46]. However, Fremout et al. [55] mentioned that tropical forests face greater threats from climate change, such as change in land use, overexploitation of forests, and overgrazing—activities that have reduced two-thirds of the original coverage in the Americas. Another serious threat is pests, which, with the increase in temperature, are expected to shorten their life cycle and can cause further damage. Therefore, the species under study is highly vulnerable to all these threats.
This study is the first to report the potential distribution of C. odorata based on its thermal niche for seed germination. The approach developed here constitutes an alternative methodology for estimating the potential distribution of species, considering the optimum germination temperature range, which is the most critical stage in the establishment and development of plant species. The results of this study complement, in ecophysiological terms, the conventional methods of studies of potential distribution and allow for the visualization of the perspectives including physiological, ecophysiological, and phenological components in the bioclimatic variables that are conventionally used to define the distribution niches of the species.

4. Materials and Methods

4.1. Study Species

Cedrela odorata L. (Meliaceae) is a tree native to tropical America (Figure 4) [56], distributed from northern Mexico to northern Argentina, up to 1200 m a.s.l. [25]. It is found in the humid and subhumid tropics and grows in association with tropical deciduous, sub-deciduous, and rain forests, as well as tropical montane cloud forest [57]. In Mexico, it is distributed along the Gulf slope, from southern Tamaulipas and southeastern San Luis Potosí to the Yucatan Peninsula. On the Pacific slope, it is distributed from southern Sonora to southern Guerrero, the Central Depression, and the coast of Chiapas [58], where it can be found in warm and semi-warm climates, and it requires soil that is fertile with good drainage [57].

4.2. Study Area

The study area was defined for Mexico according to the database of the National Inventory of Forests and Soils 2009–2013, in which the distribution of C. odorata was delimited with the humid and subhumid tropics and grows in association with tropical deciduous forest, sub-deciduous, and evergreen or rain forests, as well as tropical montane cloud forest (Figure 5) [20,58].

4.3. Data Acquisition

A georeferenced database of C. odorata was constructed from the presence data of the species collected in the country using the information contained in the Global Biodiversity Information Facility (GBIF) platform of copies deposited in herbaria around the world (https://www.gbif.org/ accesed on 18 November 2021) [59]. The database was supplemented with records obtained in field collections by personnel from the Seed Bank of the Iztacala Faculty of Higher Studies, UNAM. The database was filtered to eliminate data that were incomplete, as well as when located in urban areas, roadways, rivers, and agricultural areas, for which the land use map was used [60]. For the construction of the distribution models of C. odorata, a special filtering of the records of the presence of the species was carried out, which consisted of selecting only the records that met the condition of being in the optimal range of cardinal temperatures for germination, determined for the species in a previous work [7]. According to the nearest neighbor method in ArcMap 10.3 and at the spatial resolution of the layers in vector format, the minimum spatial distance between occurrence records was 5 km [61].

4.4. Climate Variables

Nineteen environmental variables were used at a spatial resolution of 30 arc seconds, of which eleven were temperature variables. These data were obtained from the BioClim database (Table 2) from the period of 1970 to 2000 for the current distribution [62], at a resolution of 1 km2 per pixel [61]. The layers, in vector format (INEGI, VI series layer) of land use and vegetation, were taken from CONABIO [63].

4.5. Distribution Model and Optimum Germination Temperature Range

To model the distribution of C. odorata according to the cardinal temperatures for seed germination, we used the values previously identified in a previous study by Sampayo-Maldonado et al. [7]. These authors determined a minimum temperature (Tb) of −0.5 ± 0.09 °C, an optimum temperature (To) of 38 ± 1.6 °C, and a maximum temperature (Tc) of 53.3 ± 2.1 °C for this species [7]. To define the optimum germination temperature range, the normal behavior assumption of germination values was presumed. In this sense, considering the mean value of To ± 2 standard deviations, the optimal range of 34.8 to 41.2 °C was defined, which included 95% of the population in this parameter. To model the distribution of C. odorata, the collection sites that were located in the zones that meet the optimum germination temperature range were selected.

4.6. Modeling Current Distributions of C. odorata

Species distribution models (SDM) were built with MaxEnt (software version 3.4.1®) (Steven J. Phillips, Miroslav Dudik and Robert E. Schapire. Nueva York, NY, USA) to carry out the analyses according to the collection sites and distribution of C. odorata based on the maximum entropy model. The algorithm estimates the probable distribution of species and uses a set of environmental variables to determine the maximum entropy distribution. High predictive precision is achieved, in a logistical format, by improving the calibration of the model, which provides greater representativeness of suitability [64]. In the analysis, 1373 georeferenced sites with the presence of populations of the species were used. Of the sites with the presence of the species, 256 were used and were those that were distributed in the optimum germination temperature range (34.8 to 41.2 °C). These numbers of sites was considered a sufficient sample size [65], so the prediction of the geographic distribution for the species was acceptable.
To generate the current distribution of optimal climate habitat for C. odorata, BioClim version 2.0 (Richard G. Pearson. MD, USA) [62] http://www.worldclim.org accessed on 18 November 2021 was used, with data from the period 1970 to 2000. The distribution of optimal climate habitat for the species under the current climate was modeled using the MaxEnt algorithm. This process is known as environmental niche modeling [41,66] or climate niche modeling [45]. Additionally, the jackknife method was used, implemented in MaxEnt, to indicate the relative contribution of each environment al variable in the model. With the environmental variables that contributed the most to the model, the program was executed again to generate the maps. The spatial distribution of the optimal climate habitat was obtained using the tool ArcMap 10.3®, obtaining the number of pixels and transforming to km2. Finally, the maps of the distribution of the species were generated for the current climate, with a probability greater than 20% of identifying areas where the temperature will have positive effects on seed germination in the understory.
To overcome the limitations of the distribution models of the species in this study, the Pearson correlation coefficient was used to test the multicollinearity among the 21 environmental variables. The selection of predictive variables was based on the following criteria [64,65,66,67]: those with the greatest contribution according to the jackknife test, uncorrelated (r ≤ 0.8), with normal distribution, and by the type of response of the environmental variable, in addition to taking into account previous work on the distribution of the species, eliminating areas taking into account barriers that make the dispersal of the species impossible, restricting the results by biogeographical regions, and using physiological data on the germination of the species.
To validate the model, the recommendations of Peterson and Soberón [67] were followed. The test was carried out using 30% of the data, which were separated randomly from the total localities with presence of the species. The goodness-of-fit of the predictions of the model were evaluated using the area under the curve (AUC) and the receiver operating characteristics (ROC). The results were expressed as a ratio of the ROC curve observed to a random curve, where both were truncated to the area delimited by the error threshold. The model obtained, using 70% of the data, was reclassified to values between 0 and 1, where 0 is less ideal (lower probability), and 1 is the most ideal (higher probability) for the presence of the species. The jackknife procedure was used to reduce the number of environmental variables so that only the variables that had the highest contribution to the model prediction are shown. In addition, the generated models were evaluated through the omission rate, Aikaike information criteria, or model complexity (AIC) [67]. Following the parsimony principle, we discard the most complex models and choose the simplest model, the one with the lowest (AIC), according to the Aikaike information criterion.

4.7. Model of Potential Distribution of C. odorata

Climate change scenarios were estimated on the basis of the difference between current and future climate layers according to the United States National Center for Atmospheric Research (NCAR) model (CCSM4). For this, the General Circulation Model (GCM) climate layers for CCSM4 were downloaded, which were generated on the basis of Regional Models from the Coupled Model Intercomparison Project, Phase 5 (CMIP 5), of the Intergovernmental Panel on Climate Change (IPCC), projected for a horizon of intermediate future 2050 (average for 2041–2060), and for 2070, as a horizon of far future (average of 2061–2080), with two contrasting representative concentration pathways of Watts/m2 (RCP2.6 and RCP8.5) representing two extreme scenarios, with the lowest and the highest CO2 emissions classified as conservative and drastic scenarios, respectively [68]. For the intermediate future (2050), it predicts average annual temperature increase for RCP2.6 of 1.9 and of 3.8 °C for the RCP8.5, while for the far future horizon (2070), it predicts 2.7 and 4.7 °C for the RCP2.6 and RCP8.5, respectively, of the average annual temperature increase [69].
The potential future distribution was projected using the MaxEnt algorithm (SDM), using the geographic records and the intermediate (2050) and far (2070) future temperature grids. Firstly, this tool models optimal climate habitat under the current climate. Once that model has been constructed, the projection is repeated, but this time using the model with the grids of climate variables for the decades centered on 2050 and 2070, estimated using the NCAR model (CCSM4) with an RCP2.6 and RCP8.5 Watts/m2 (lowest and the highest CO2 emissions), classified as conservative and drastic models, respectively. Finally, maps of the distribution of the species in future climate were generated, with a probability greater than 20%.

5. Conclusions

In general, due to climate change, the potential distribution of C. odorata will expand in most of the states of the country with some exceptions (Chihuahua, Mexico, Morelos, Guerrero, and Durango). The most extreme case will be the state of Chihuahua, which will lose optimal conditions for the distribution of the species. With the increase in temperature, the appearance of adequate environmental conditions for the species distribution in the state of Guanajuato is expected. From the thermal variables, the coldest month and the annual temperature oscillation are those that contributed the most to the prediction models. Therefore, on the basis or our results, it was possible to predict the potential distribution of C. odorata based on its thermal niche for seed germination. This will allow actions to be undertaken to monitor populations that are vulnerable to climate change, as well as activities such as seed collection to conserve local genotypes whose populations will disappear, and to carry out genetic studies of provenances in the new areas that will have adequate environmental conditions to the potential distribution to the species.

Author Contributions

C.M.F.-O., C.A.O.-S. and E.M. conceived the idea and designed the experiment; N.I.R.-A. collected the samples; S.S.-M. and C.A.O.-S. conducted the experiment; T.U., M.W. and E.C.-L. verified the analytical methods; P.D.D.-A., R.L.-S. and O.T.-V. supervised the development of this work and analyzed the data; S.S.-M. and C.M.F.-O. wrote the article. All authors interpreted the results and contributed to the final manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This study has been funded by the Garfield Weston Foundation as part of the Global Tree Seed Bank Project, and by the British Council under the Newton Fund Institutional Links programme with the aim to conserve and sustainably use the native flora of Mexico. Both projects are lead by the Royal Botanic Gardens, Kew.

Data Availability Statement

Geographic coordinate data of collected species and environmental variables were described in Section 4.2, Section 4.3 and Section 4.4 of the article.

Acknowledgments

We thank the Postdoctoral Fellowship Program at the UNAM, of the General Direction of Academic Personnel Affairs (DGAPA). This work was supported by UNAM-PAPIIT IA210020.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Franks, S.J.; Sim, S.; Weis, A.E. Rapid evolution of flowering time by an annual plant in response to a climate fluctuation. Proc. Natl. Acad. Sci. USA 2007, 104, 1278–1282. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Craufurd, P.Q.; Wheeler, T.R. Climate change and flowering time in annual crops. J. Exptl. Bot. 2009, 60, 2529–2539. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Anderson, J.T.; Inouye, D.W.; McKinney, A.M.; Colauti, R.I.; Mitchell-Olds, T. Pheotypic plasticity and adaptive evolution contribute to advancing flowering phenology in response to climate change. Proc. R. Soc. B 2012, 279, 3843–3852. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Iler, A.M.; Høye, T.K.; Inouye, D.W.; Schmidt, N.M. Nonlinear flowering responses to climate: Are species approaching their limits of phonological change? Phil. Trans. R. Soc. B 2013, 368, 20120489. [Google Scholar] [CrossRef] [PubMed]
  5. Gaira, K.S.; Rawal, R.S.; Rawat, B.; Bhatt, I.D. Impact of climate change on the flowering of Rhododendron arboretum in central Himalaya, India. Curr. Sci. 2014, 106, 1735–1738. [Google Scholar]
  6. Wang, B.; Liu, D.L.; Asseng, S.; Macadam, I.; Yu, Q. Impact of climate change on wheat flowering time in eastern Australia. Agric. Forest. Meteorol. 2015, 209, 11–21. [Google Scholar] [CrossRef]
  7. Sampayo-Maldonado, S.; Ordoñez-Salanueva, C.A.; Mattana, E.; Ulian, T.; Way, M.; Castillo, E.; Davila-Aranda, P.D.; Lira-Sadde, R.; Téllez-Valdéz, O.; Rodríguez-Arevalo, N.I.; et al. Thermal time and cardinal temperatures for germination of Cedrela odorata L. Forests 2019, 10, 841. [Google Scholar] [CrossRef] [Green Version]
  8. Parmesan, C.; Yohe, G. A globally coherent fingerprint of climate change impacts across natural systems. Nature 2003, 421, 37–42. [Google Scholar] [CrossRef]
  9. Root, T.L.; MacMynowski, D.P.; Mastrandrea, M.D.; Schneider, S.H. Human modified temperature induce species changes: Joint attribution. Proc. Natl. Acad. Sci. USA 2005, 102, 7465–7469. [Google Scholar] [CrossRef] [Green Version]
  10. Iverson, L.R.; Prasad, A.M.; Schwartz, M.W. Predicting potential changes in suitable habitat and distribution by 2100 for tree species in the eastern United States. J. Agric. Meteorol. 2005, 61, 29–37. [Google Scholar] [CrossRef] [Green Version]
  11. Iverson, L.R.; Prasad, A.M.; Matthews, S.N.; Peters, M. Estimating potential habitat for 134 Eastern US tree species under six climate scenarios. For. Ecol. Manage 2008, 254, 390–406. [Google Scholar] [CrossRef]
  12. McKenney, D.W.; Pedlar, J.H.; Lawrence, K.; Campbell, K.; Hutchison, M.F. Potential impacts of climate change on the distribution of North American trees. Bioscience 2007, 57, 939–948. [Google Scholar] [CrossRef]
  13. Kelly, A.E.; Goulden, M.L. Rapid shifts in plant distribution with recent climate change. Proc. Natl. Acad. Sci. USA 2007, 105, 11823–11826. [Google Scholar] [CrossRef] [Green Version]
  14. Thuiller, W.; Albert, C.; Araújo, M.B. Predicting global change impacts on plant species distributions: Future challenges. Perspect. Plant. Ecol. Evol. Syst. 2008, 9, 137–152. [Google Scholar] [CrossRef]
  15. Pautasso, M. Forest ecosystems and global change: The case study of insurbia. Ann. Bot. 2013, 3, 1–29. [Google Scholar]
  16. Iverson, L.R.; McKenzie, D. Tree-species range shifts in the changing climate: Detecting, modeling, assisting. Landsc. Ecol. 2013, 28, 879–889. [Google Scholar] [CrossRef]
  17. Kurpis, J.; Serrato-Cruz, M.A.; Arroyo, T.P.F. Modeling the effects of climate change on the distribution of Tagetes lucida Cav. (Asteraceae). Glob. Ecol. Conserv. 2019, 20, e00747. [Google Scholar] [CrossRef]
  18. Walk, J.L.; Hidayati, S.; Dixon, K.W.; Thompson, K.; Poschlod, P. Climate change and plant regeneration from seed. Glob. Chang. Biol. 2011, 17, 2145–2161. [Google Scholar] [CrossRef]
  19. Fernández-Pascual, E.; Mattana, E.; Pritchard, H.W. Seeds for future past: Climate change and the thermal memory of plant reproductive traits. Biol. Rev. 2019, 94, 439–456. [Google Scholar] [CrossRef]
  20. Romo-Lozano, J.L.; Vargas-Hernández, J.; López-Upton, J.; Ávila-Angulo, M.L. Estimate of the financial value of timber stocks of red cedar (Cedrela odorata L) in Mexico. Madera y Bosques 2017, 23, 111–120. [Google Scholar] [CrossRef] [Green Version]
  21. Sampayo-Maldonado, S.; Jiménez-Casas, M.; López-Upton, J.; Sánchez-Monsalvo, V.; Jasso-Mata, J.; Equihua-Martínez, A.; Castillo-Martínez, C.R. Cedrela odorata L. mini-cutting rooting. Agrociencias 2016, 50, 919–929. [Google Scholar]
  22. Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES). In Proceedings of the CoP14 Prop. 33, Fourteenth Meeting of the Conference of the Parties, The Hague, The Netherlands, 3–15 June 2007.
  23. Secretaria del Medio Ambiente y Recursos Naturales (SEMARNAT). NORMA Oficial Mexicana NOM-059-SEMARNAT-2010, Protección Ambiental -Especies Nativas de México de Flora & Fauna Silvestres- Categorías de Riesgo & Especificaciones Para su Inclusión, Exclusión o cambio-Lista de Especies en Riesgo. Available online: http://www.cec.org/wp-content/uploads/wpallimport/files/18-2-det1412_es.pdf (accessed on 12 September 2022).
  24. Chun-Jing, W.; Zhi-Xiang, Z.; Ji-Zhong, W. Vulnerability of global Forest ecoregions to future climate. Glob. Ecol. Conserv. 2019, 20, 1–10. [Google Scholar] [CrossRef]
  25. Gómez-Díaz, J.D.; Monterroso-Rivas, A.I.; Tinoco-Rueda, J.A. Distribution of red cedar (Cedrela odorata L.) in the Hidalgo state, under current conditions and scenarios of climate change. Madera y Bosques 2007, 13, 29–49. [Google Scholar] [CrossRef]
  26. Hernández-Ramos, J.; Reynoso-Santos, R.; Hernández-Ramos, A.; García-Cuevas, X.; Hernández-Máximo, E.; Cob-Uicab, J.V.; Sumano-López, D. Historical, current and future distribution of Cedrela odorata in Mexico. Acta Bot. Mex. 2018, 124, 117–134. [Google Scholar] [CrossRef]
  27. CICS. Pridicting Species Distribution Use of Climatic Parametrs. Canadian Institute for Climate Studies. 2006. Available online: http://www.cics.uvic.ca/ (accessed on 13 December 2022).
  28. Duan, X.; Li, J.; Wu, S. MaxEnt Modeling to Estimate the Impact of Climate Factors on Distribution of Pinus densiflora. Forests 2022, 13, 402. [Google Scholar] [CrossRef]
  29. Duval, V.S.; Benedetti, G.M.; Campo, A.M. Relación clima-vegetación: Adaptaciones de la comunidad del jarillal al clima semiárido, Parque Nacional Lihué Calel, provincia de la Pampa, Argentina. Investig. Geog. 2015, 88, 33–44. [Google Scholar] [CrossRef] [Green Version]
  30. Funes, G.; Díaz, S.; Venier, P. La temperatura como principal determinante de la germinación en especies del Chaco seco de Argentina. Ecol. Austral 2009, 19, 129–138. [Google Scholar]
  31. Durán-Puga, N.; Ruíz-Corral, J.A.; González-Eguiarte, D.R.; Núñez-Hernández, G.; Padilla-Ramírez, F.J.; Contreras-Rodríguez, S.H. Development cardinal temperatures of the planting-emergence stage for 11 forage grasses. Rev. Mex. Cienc. Pecu. 2011, 2, 347–357. [Google Scholar]
  32. Nakao, E.A.; Cardoso, V.J.M. Analysis of thermal dependence on the germination of braquiarão seeds using the thermal time model. Braz. J. Biol. 2016, 76, 162–168. [Google Scholar] [CrossRef] [Green Version]
  33. González, J.E. Seed collection and germination of 26 tree species from the tropical rain forest. Rev. Biol. Trop. 1991, 39, 47–51. [Google Scholar] [CrossRef]
  34. Garcia-Huidobro, J.; Monteith, J.L.; Squire, G.R. Time, temperature and germination of pearl millet (Pennisetum typhoides S. & H.). I. Constant temperature. J. Exp. Bot. 1982, 33, 288–296. [Google Scholar]
  35. Gilbertson, P.K.; Berti, M.T.; Johnson, B.L. Borage cardinal germination temperatures and seed development. Ind. Crops Prod. 2014, 59, 202–209. [Google Scholar] [CrossRef]
  36. Calzada-López, S.G.; Kohashi-Shibata, J.; Uscanga-Mortera, E.; García-Esteva, A.; Yáñez-Jiménez, P. Cardinal temperatures and germination rate in husk tomato cultivars. Rev. Mex. Cienc. Agrí. 2014, 8, 1451–1458. [Google Scholar]
  37. Ruíz-Corral, J.A.; Flores-López, H.E.; Ramírez-Díaz, J.L.; González-Eguiarte, D.R. Cardinal temperatures and length of maturation cycle of maize hybrid H-311 under rainfed conditions. Agrociencia 2002, 36, 569–577. [Google Scholar]
  38. Pernús-Alvarez, M.; Sánchez-Rendón, J.A. Efecto de la procedencia sobre rasgos seminales y requerimientos germinativos de Cedrela odorata (Meliaceae). Revista Jard. Bot. Nac. 2016, 37, 203–211. [Google Scholar]
  39. Adam, N.R.; Dierig, D.A.; Coffelt, T.A.; Wintermeyer, M.J.; Mackey, B.E.; Wall, G.W. Cardinal temperatures for germination and early growth of two Lesquerella species. Ind. Crops Prod. 2007, 25, 24–33. [Google Scholar] [CrossRef]
  40. Hardegree, S.P. Predicting germination response to temperature. I. Cardinal temperature models and subpopulation-specific regression. Ann. Bot. 2006, 97, 1115–1125. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Phillips, S.J.; Anderson, R.P.; Schapire, R.E. Maximum entropy modeling of species geographic distributions. Ecol. Modell. 2006, 190, 231–259. [Google Scholar] [CrossRef] [Green Version]
  42. Sáenz-Romero, C.; Rehfeldt, G.E.; Crookston, N.L.; Duval, P.; St-Amant, R.; Beaulieu, J.; Richardson, B.A. Spline models of contemporary, 2030, 2060 and 2090 climates for Mexico and their use in understanding climate-plant impacts on vegetation. Clim. Chang. 2010, 102, 595–623. [Google Scholar] [CrossRef]
  43. Sáenz-Romero, C.; Rehfeldt, G.E.; Duval, P.; Lindig-Cisneros, R.A. Abies religiosa habitat prediction in climatic change scenarios and implications for monarch butterfly conservation in Mexico. For. Ecol. Manag. 2012, 275, 98–106. [Google Scholar] [CrossRef]
  44. Sáenz-Romero, C.; Rehfeldt, G.E.; Ortega-Rodríguez, J.M.; Marín-Togo, M.C.; Madrigal-Sánchez, X. Pinus leiophylla suitable habitat for 1961-1990 and future climate. Bot. Sci. 2015, 93, 709–718. [Google Scholar] [CrossRef] [Green Version]
  45. Garza-López, M.; Ortega-Rodríguez, L.M.; Zamudio-Sánchez, F.J.; López-Toledo, J.F.; Domínguez-Álvarez, F.A.; Sáenz-Romero, C. Calakmul como refugio de Swietenia macrophylla King ente el cambio climático. Bot. Sci. 2016, 94, 43–50. [Google Scholar] [CrossRef] [Green Version]
  46. Sampayo-Maldonado, S.; Ordoñez-Salanueva, C.A.; Mattana, E.; Way, M.; Castillo-Lorenzo, E.; Dávila-Aranda, P.D.; Lira-Sadde, R.; Téllez-Valdéz, O.; Rodríguez-Arevalo, N.I.; Ulian, T.; et al. Thermal Niche for Seed Germination and Species Distribution Modelling of Swietenia macrophylla King (Mahogany) under Climate Change Scenarios. Plants 2021, 10, 2377. [Google Scholar] [CrossRef] [PubMed]
  47. Villaseñor, J.L.; Ortiz, E. Biodiversity of flowering plants (Division Magnoliophyta) in Mexico. Rev. Mex. Biodivers. 2014, 85, 134–142. [Google Scholar] [CrossRef] [Green Version]
  48. Hartmann, H.T.; Kester, D.E.; Davies, F.T.; Geneve, R.L. Plant Propagation: Principles and Practices, 8th ed.; Prentice-Hall: Hoboken, NJ, USA, 2013; 913p. [Google Scholar]
  49. Beigaité, R.; Tang, H.; Bryn, A.; Skarpaas, O.; Stordal, F.; Bjerke, J.W.; Zliobaite, I. Identifying climate thresholds for dominant natural vegetation types at the global scale using machine learning: Average climate versus extremes. Glob. Chang. Biol. 2022, 28, 3557–3579. [Google Scholar] [CrossRef]
  50. Leal-Narez, O.; Mendoza, M.E.; Pérez-Salicrup, D.; Geneletti, D.; López-Granados, E.; Carranza, E. Distribución potencial de Pinus Martinezii: Un modelo espacial basado en conocimiento ecológico y análisis multicriterio. Rev. Mex. Biodivers. 2012, 83, 1152–1170. [Google Scholar] [CrossRef]
  51. Gutiérrez, E.; Trejo, I. Effect of climatic change on the potential distribution of five species of temperate forest trees in Mexico. Rev. Mex. Biodivers. 2014, 85, 179–188. [Google Scholar] [CrossRef] [Green Version]
  52. Parmoon, G.; Moosavi, S.A.; Akbari, H.; Ebadi, A. Quantifying cardinal temperatures and thermal time required for germination of Silybum marianum seed. Crop J. 2015, 3, 145–151. [Google Scholar] [CrossRef] [Green Version]
  53. Daibes, L.F.; Amoêdo, S.C.; Nascimento-Moraes, J.; Fenelon, N.; da Silva, D.R.; de Melo-Lopes, M.J.; Vargas, L.A.; Monteiro, E.F.; Frigeri, R.B.C. Thermal requirements of seed germination of ten tree species occurring in the western Brazilian Amazon. Seed Sci. Res. 2019, 29, 115–123. [Google Scholar] [CrossRef]
  54. Ortiz-Solorio, C.A. Elementos de Agrometeorología Cuantitativa Con Aplicaciones en la República Mexicana; Universidad Autónoma Chapingo: Texcoco, Mexico, 1987; 327p. [Google Scholar]
  55. Fremout, T.; Thomas, E.; Gaisberger, H.; Van-Meerbeek, K.; Muenchow, J.; Briers, S.; Gutierrez-Miranda, C.E.; Marcelo-Peña, J.L.; Kindt, R.; Atkinson, R.; et al. Mapping tree species vulnerability to multiple threats as a guide to restoration and conservation of tropical dry forests. Glob. Chang. Biol. 2020, 26, 3552–3568. [Google Scholar] [CrossRef]
  56. Pennington, T.D.; Sarukhán, J. Tropical trees of Mexico: Manual for the Identification of the Main Species, 3rd ed.; University Scientific Text; UNAM: Mexico City, Mexico, 2005; pp. 294–296. [Google Scholar]
  57. Mendizábal-Hernández, L.C.; Alba-Landa, J.; Márquez-Ramírez, J.; Cruz-Jiménez, H.; Ramírez-García, E.O. Carbon sequestration by Cedrela odorata L. in a genetic trial. Rev. Mex. Cienc. Forest. 2011, 2, 107–113. [Google Scholar]
  58. Instituto Nacional de Estadística y Geografía (INEGI). Conjunto de Datos Vectoriales Edafológicos, Escala 1: 250,000 Serie IV; Continuo Nacional, Escala: 1: 250,000 Ed. 2; Instituto Nacional de Estadística y Geografía: Aguascalientes, Mexico, 2014. [Google Scholar]
  59. Global Biodiversity Information Facility. GBIF Ocurrence Cedrela odorata L. Available online: https://api.gbif.org/v1/occurrence/download/request/0216476-220831081235567.zip (accessed on 18 November 2021).
  60. Instituto Nacional de Estadística y Geografía (INEGI). Conjunto De Datos Vectoriales de Uso Del Suelo Y Vegetación. Escala 1: 250,000 Serie VI; Continuo Nacional, Escala: 1: 250,000 Ed. 2; Instituto Nacional de Estadística y Geografía: Aguascalientes, Mexico, 2018. [Google Scholar]
  61. Hijmans, R.j.; Cameron, S.E.; Parra, J.L.; Jones, P.G.; Jarvis, A. Very higt resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 2005, 25, 1965–1978. [Google Scholar] [CrossRef]
  62. WorldClim. Free Climate Data for Ecological Modeling and GIS. WorldClim-Global Climate Data. 2006. Available online: http://www.worldclim.org/ (accessed on 18 November 2021).
  63. Comisión Nacional para el Conocimiento y Uso de la Biodiversidad (CONABIO). Red mundial de información sobre biodiversidad (REMIB). Comisión Nacional para el Conocimiento & Uso de la Biodiversidad. Mexico City, México. 2015. Available online: https://www.academia.edu/23597234/Comisi%C3%B3n_Nacional_para_el_Conocimiento_y_Uso_de_la_Biodiversidad_Comisi%C3%B3n_Nacional_de_%C3%81reas_Naturales_Protegidas (accessed on 25 November 2020).
  64. Phillips, S.J.; Anderson, R.P.; Dudík, M.; Schapire, R.E.; Blair, M.E. Opening the black box: An open-source release of MaxEnt. Ecography 2017, 40, 887–893. [Google Scholar] [CrossRef]
  65. Baldwin, R.A. Use of maximum entropy modeling in wildlife research. Entropy 2009, 11, 854–866. [Google Scholar] [CrossRef]
  66. Elith, J.; Phillips, S.J.; Hastie, T.; Dudík, M.; Chee, Y.E.; Yates, C.J. A statistic explanation of MaxEnt for ecologists. Divers. Distrib. 2011, 17, 43–57. [Google Scholar] [CrossRef]
  67. Peterson, A.T.; Soberón, J. Rethinking receiver operating characteristic analysis applications in ecological niche modelling. Ecol. Modell. 2008, 213, 63–72. [Google Scholar] [CrossRef]
  68. Fernández-Eguiarte, A.; Zavala-Hidalgo, J.; Romero-Centeno, R.; Conde-Álvarez, A.C.; Trejo-Vázquez, R.I. Actualización de los escenarios de Cambio Climático Para Estudios de Impactos, Vulnerabilidad y Adaptación en México y Centroamérica; Centro de Ciencias de la Atmósfera, Universidad Nacional Autónoma de México. Instituto Nacional de Ecología y Cambio Climático: Mexico City, Mexico, 2015; 22p. [Google Scholar]
  69. Weyant, J.; Azar, C.; Kainuma, M.; Kejun, J.; Nakicenovic, N.; Shukla, P.R.; La Rovere, E.; Yohe, G. Future Intergovernmental Panel on Climate Change activities- New scenarios. Report of 2.6 Versus 2.9 Watts/m2 Representative Concentration Pathways (RCP´s) Evaluation Panel. In Proceedings of the Intergovernmental Panel on Climate Change (IPCC)-XXX/INF.6, Antalya, Turkey, 21–23 April 2009. [Google Scholar]
Figure 1. Current distribution of C. odorata. In green, the optimal distribution according to the optimum germination temperature range.
Figure 1. Current distribution of C. odorata. In green, the optimal distribution according to the optimum germination temperature range.
Plants 12 00150 g001
Figure 2. Potential distribution for the intermediate future (2050) of C. odorata in the NCAR model (CCSM4); for the conservative (RCP2.6) and drastic (RCP8.5) scenarios according to the optimum germination temperature range.
Figure 2. Potential distribution for the intermediate future (2050) of C. odorata in the NCAR model (CCSM4); for the conservative (RCP2.6) and drastic (RCP8.5) scenarios according to the optimum germination temperature range.
Plants 12 00150 g002
Figure 3. Potential distribution for the far future (2070) of C. odorata, in the NCAR model (CCSM4); for the conservative (RCP2.6) and drastic (RCP8.5) scenarios according to the optimal germination temperature range.
Figure 3. Potential distribution for the far future (2070) of C. odorata, in the NCAR model (CCSM4); for the conservative (RCP2.6) and drastic (RCP8.5) scenarios according to the optimal germination temperature range.
Plants 12 00150 g003aPlants 12 00150 g003b
Figure 4. Seeds, seedlings, and young trees of C. odorata.
Figure 4. Seeds, seedlings, and young trees of C. odorata.
Plants 12 00150 g004
Figure 5. Distribution of the tropical forest in Mexico, as area of study for the distribution of C. odorata.
Figure 5. Distribution of the tropical forest in Mexico, as area of study for the distribution of C. odorata.
Plants 12 00150 g005
Table 1. Surface with potential distribution of C. odorata according to climate change scenarios by state in the model (CCSM4).
Table 1. Surface with potential distribution of C. odorata according to climate change scenarios by state in the model (CCSM4).
StateCurrent (km2)2050 (km2)2070 (km2)
RCP2.6RCP8.5RCP2.6RCP8.5
Campeche49,206.3746,696.0346,596.3345,717.0546,299.80
Chiapas26,474.2731,882.0232,529.2032,424.9433,982.65
Chihuahua433.4528.612.2316.59-
Colima335.29572.07208.92475.87929.28
Durango1508.061128.22773.411101.02759.43
Guanajuato-15.386.6465.102.70
Guerrero6033.252138.622268.992150.802684.70
Hidalgo668.781867.851566.701831.651755.27
Jalisco5044.108603.378863.379131.178990.25
Mexico433.5467.7680.2246.474.41
Michoacan435.001391.38409.871,301.65798.17
Morelos802.23122.1570.7162.8272.15
Nayarit8899.308895.379016.189022.268872.49
Oaxaca16,410.2020,896.9220,852.3218,317.9519,579.65
Puebla1341.35848.881368.411309.981862.81
Queretaro71.61137.7796.53158.1265.26
Quintana Roo20,136.1917,307.7516,330.7815,596.0118,351.02
San Luis Potosi6848.314888.475112.184445.605200.63
Sinaloa6586.409883.499728.929754.789823.81
Sonora16.4511.4232.583.7831.98
Tabasco14,217.6617,776.3916,326.2518,410.8516,609.55
Tamaulipas4623.803757.833467.143779.013526.20
Veracruz24,716.4631,145.7531,821.3831,705.6532,537.99
Yucatan33,538.9633,827.3534,003.9533,109.5533,941.25
Suitability area (km2)228,780.66243,890.85241,533.21239,938.67246,681.45
Increment area (km2)N/A15,110.1912,752.5511,158.0117,900.79
Percentage of increment (%) N/A6.605.574.877.82
Table 2. Environmental variables used in this study.
Table 2. Environmental variables used in this study.
AbbreviationEnvironmental VariableUnitsUsed for Modeling
Bio1Annual average temperature(°C)No
Bio2Diurnal temperature oscillation(°C)Yes
Bio3Isothermality (Bio2/Bio7) × 100(°C)No
Bio4Temperature seasonality
(standard deviation × 100)
(°C)Yes
Bio5Maximum average temperature of the warmest period(°C)No
Bio6Minimum temperature of the coldest month(°C)Yes
Bio7Annual temperature oscillation (Bio5–Bio6)(°C)Yes
Bio8Average temperature of the wettest month(°C)No
Bio9Average temperature of the driest month(°C)Yes
Bio10Average temperature of the warmest quarter(°C)No
Bio11Average temperature of the coldest quarter(°C)No
Bio12Annual precipitation(mm)Yes
Bio13Precipitation during the wettest period(mm)No
Bio14Precipitation during the driest period(mm)Yes
Bio15Precipitation seasonality
(coefficient of variation)
CVNo
Bio16Precipitation during the wettest trimester(mm)No
Bio17Precipitation during the driest trimester(mm)No
Bio18Precipitation during the warmest trimester(mm)No
Bio19Precipitation during the coldest trimester(mm)No
AltitudeElevation(m)No
SoilsSoil type-Yes
Source: WorldClim [62]; CONABIO [63].
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Sampayo-Maldonado, S.; Ordoñez-Salanueva, C.A.; Mattana, E.; Way, M.; Castillo-Lorenzo, E.; Dávila-Aranda, P.D.; Lira-Saade, R.; Téllez-Valdés, O.; Rodríguez-Arévalo, N.I.; Flores-Ortiz, C.M.; et al. Potential Distribution of Cedrela odorata L. in Mexico according to Its Optimal Thermal Range for Seed Germination under Different Climate Change Scenarios. Plants 2023, 12, 150. https://doi.org/10.3390/plants12010150

AMA Style

Sampayo-Maldonado S, Ordoñez-Salanueva CA, Mattana E, Way M, Castillo-Lorenzo E, Dávila-Aranda PD, Lira-Saade R, Téllez-Valdés O, Rodríguez-Arévalo NI, Flores-Ortiz CM, et al. Potential Distribution of Cedrela odorata L. in Mexico according to Its Optimal Thermal Range for Seed Germination under Different Climate Change Scenarios. Plants. 2023; 12(1):150. https://doi.org/10.3390/plants12010150

Chicago/Turabian Style

Sampayo-Maldonado, Salvador, Cesar A. Ordoñez-Salanueva, Efisio Mattana, Michael Way, Elena Castillo-Lorenzo, Patricia D. Dávila-Aranda, Rafael Lira-Saade, Oswaldo Téllez-Valdés, Norma I. Rodríguez-Arévalo, Cesar M. Flores-Ortiz, and et al. 2023. "Potential Distribution of Cedrela odorata L. in Mexico according to Its Optimal Thermal Range for Seed Germination under Different Climate Change Scenarios" Plants 12, no. 1: 150. https://doi.org/10.3390/plants12010150

APA Style

Sampayo-Maldonado, S., Ordoñez-Salanueva, C. A., Mattana, E., Way, M., Castillo-Lorenzo, E., Dávila-Aranda, P. D., Lira-Saade, R., Téllez-Valdés, O., Rodríguez-Arévalo, N. I., Flores-Ortiz, C. M., & Ulian, T. (2023). Potential Distribution of Cedrela odorata L. in Mexico according to Its Optimal Thermal Range for Seed Germination under Different Climate Change Scenarios. Plants, 12(1), 150. https://doi.org/10.3390/plants12010150

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop