Next Article in Journal
Use of Gamma Radiation for the Genetic Improvement of Underutilized Plant Varieties
Previous Article in Journal
Distribution Patterns of Essential Oil Terpenes in Native and Invasive Solidago Species and Their Comparative Assessment
Previous Article in Special Issue
Survey of the Solidago canadensis L. Morphological Traits and Essential Oil Production: Aboveground Biomass Growth and Abundance of the Invasive Goldenrod Appears to Be Reciprocally Enhanced within the Invaded Stands
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 11
Issue 9
10.3390/plants11091160
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Distribution Patterns of Invasive Buffelgrass (Cenchrus ciliaris) in Mexico Estimated with Climate Niche Models under the Current and Future Climate

by
Pablo Siller-Clavel
1,
Ernesto I. Badano
2,
Federico Villarreal-Guerrero
1,
Jesús A. Prieto-Amparán
1,
Alfredo Pinedo-Alvarez
1,
Raúl Corrales-Lerma
1,
Alan Álvarez-Holguín
3 and
Nathalie S. Hernández-Quiroz
1,*
1
Facultad de Zootecnia y Ecología, Universidad Autónoma de Chihuahua, Periférico Francisco R. Almada km. 1, Chihuahua 31453, Mexico
2
IPICYT División de Ciencias Ambientales, Instituto Potosino de Investigación Científica y Tecnológica, Camino a la Presa San José 2055, Colonia Lomas 4a Sección, San Luis Potosí 78216, SLP, Mexico
3
Instituto Nacional de Investigaciones Forestales, Agrícolas y Pecuarias (INIFAP), Campo Experimental La Campana, Carretera Chihuahua-Ojinaga km. 33.3, Aldama 32190, Mexico
*
Author to whom correspondence should be addressed.
Plants 2022, 11(9), 1160; https://doi.org/10.3390/plants11091160
Submission received: 7 March 2022 / Revised: 1 April 2022 / Accepted: 18 April 2022 / Published: 26 April 2022
(This article belongs to the Special Issue Ecology of Invasive Plants)
Browse Figures
Versions Notes

Abstract

:
In Mexico, buffelgrass (Cenchrus ciliaris) was introduced in the middle of the 20th century. Currently, buffelgrass has become an invasive species and has colonized various ecosystems in the country. In addition to its invasive capacity, climate change is a factor that has to be taken into account when considering how to effectively manage and control this species. The climatic niche models (CNM) and their projections for climate change scenarios allow for estimating the extent of biological invasions. Our study aimed to calibrate a CNM for buffelgrass in Mexico under the current climatic conditions and to project the extent of its biological invasion under climate change scenarios. For that, we used MaxEnt to generate the current CNM and to detect if climate change could cause future changes, we then evaluated the distribution patterns over the periods of 2041–2060, 2061–2080, and 2081–2100 for all the shared socioeconomic pathways (SSPs). Linear regressions were used to compare the outputs between current and future scenarios. Under the current climate, the CNM estimated that 42.2% of the continental surface of Mexico is highly suitable for buffelgrass. The regression analyses indicated no effects from climate change on the distribution of buffelgrass. Moreover, when the projected period is further in the future, and when the SSPs intensify, the surface of suitable areas for the species increases. These analyses clearly suggest Mexico is facing a biological invasion from buffelgrass, which may represent a threat to native biodiversity.

Graphical Abstract

1. Introduction

Buffelgrass (Cenchrus ciliaris L., Poaceae) is a perennial species native to arid and semi-arid ecosystems of Africa. It can tolerate high temperatures that approach 50 °C [1,2]. This species was introduced to southern USA and northern Mexico for cattle foraging in the middle of the 20th century due to its resistance to drought conditions [1,3]. In Mexico, this grass spreads across semi-arid, subtropical, and tropical climates at elevations ranging from 500 to 2000 m. It is also found in marginal habitats with extreme temperatures (5.8–40.6 °C) and nutrient-poor soils, which are not usually suitable for the development of native plants [4,5]. The invasive status of this exotic grass in America, however, is a matter of debate. On the one hand, some authors postulated that buffelgrass is not an invasive plant because it neither disperses away from its introduction sites nor affects native species [2,6]. On the other hand, studies indicated buffelgrass is an aggressive invader because it quickly colonizes disturbed habitats, outcompeting native plants and altering ecosystem functioning by increasing the intensity and frequency of fires [7,8]. This affects the resilience capacity of habitats [9,10], thereby reducing the provision of ecosystems services and restricting the capacity of ecosystems to deal with a pest, among other events, which might affect the ecosystem’s resilience and make it prone to invasive species [11,12,13]. In Mexico, buffelgrass is listed in the “Agreement to determine the list of invasive exotic species for Mexico” due to the associated impacts on native biodiversity, the economy, and public health that is causes in the country [14,15]. For these reasons, and given the elevated diversity of ecosystems, which buffelgrass has colonized in Mexico, there is a need for developing tools to appropriately manage and control this species. Before the development of these tools, knowing about the potential extension (distribution) that buffelgrass could have in the country is necessary.
As climate is the primary driver of the distribution of species, modelling the climatic niches of exotic plants (i.e., the fraction of the fundamental niche that embraces those climatic factors that determine the survival of species) allows for estimating the extent that biological invasions can reach [16,17,18]. These climate niche models (hereafter, CNM) result from correlating the spatial distribution of a species (i.e., georeferenced occurrence data) with the spatial distribution of bioclimatic variables (i.e., climatic variables derived from temperature and rainfall that explain the distribution of living organisms). The geographical expression of the CNM (i.e., their projection on the physical space using geographic information systems) are maps indicating the probability of finding the target species in a given region as a function of the values these variables take across different habitats [19,20,21]. In this way, the potential distribution range of an invasive species can be estimated by determining what habitats in a region have an elevated likelihood of containing the environmental conditions for meeting its climatic niche [18].
All CNM are calibrated under the current climatic conditions because species occurrence data constitute historical records [22]. However, the management and control of invasive plants now require future projections due to the advance of climate change. The latest general circulation models of the World Climate Research Programme 6 (CMIP6 models, https://www.worldclim.org/data/cmip6/cmip6climate.html (accessed on 31 July 2021)) predict that global temperature will progressively rise during this century. This will be caused by the increased radiative forcing (i.e., the difference between in-going and outgoing energy in the Earth’s system in terms of W/m2), which will result from growing atmospheric concentrations of greenhouse gasses. This will also alter rainfall regimens and intensify drought in several regions [23,24]. The magnitude of these climatic changes, however, will depend on whether human societies choose to reduce carbon-dependent energy sources or continue their use. As these decisions are influenced by local and global environmental policies, the advance of carbon-free technologies, and human population growth, among other factors, four shared socioeconomic pathways (SSPs) were proposed. These SSPs describe different trends of greenhouse gases emissions, [25] and, consequently, each SSP has an associated radiative forcing level (2.6, 4.5, 7.0 and 8.5 W/m2) [26]. Therefore, the future distribution ranges of invasive plants can be predicted for different periods and SSPs projecting their CNM on the climatic envelopes are generated by general circulation models [27].
Considering these issues, our study was aimed (1) to calibrate a CNM for buffelgrass in Mexico to assess its invasibility under the current climatic conditions (i.e., their probability to invade; Lonsdale [28]) and (2) to evaluate the effect of climate change scenarios on the projection of the extent of this biological invasion. These projections can generate three different results, including expansions, contractions, or no changes in the distribution range of the target compared to what it has under the current climate [29,30]. As we are unable to predict the exact climatic conditions that will occur during this century, the future distribution of buffelgrass was predicted at different periods under contrasting SSPs. This information will allow for the development of adaptative control programs for this exotic plant.

2. Results

Most of the occurrence records used to calibrate the CNM of buffelgrass fell in Mexico (271 points = 93%), while just a few of them fell in other countries of North and Central America (17 in USA, 2 in Guatemala, and 2 in Honduras). The average AUC of the integrative CNM that resulted from the 100 bootstrap runs was 0.949 (±0.003 SE). The bioclimatic variable with the highest explanatory power in the CNM was precipitation seasonality, which accounted for 31.3% of the variance in the distribution of occurrence probabilities of buffelgrass. These values increased with higher precipitation seasonality (Figure 1a). The variable with the second highest explanatory power was the temperature annual range (explaining 23.0% of the variance; occurrence probabilities peaked around 25 °C; Figure 1). The third most important explanatory variable was the average temperature of the driest quarter of the year (explaining 19.3% variance; occurrence probabilities peaked between 15 and 25 °C; Figure 1). The fourth explanatory variable was the precipitation during the coldest quarter of the year (explaining 16.8% of variance; occurrence probabilities peaked continuously and decreased as precipitation increased; Figure 1). The bioclimatic variables with the lowest contributions to explain CNM variance were the average temperature during the warmest quarter of the year and the precipitation during the wettest quarter of the year. The former explained 5.7% of the variance and its occurrence probabilities continuously increased until temperature reached 30 °C (Figure 1), while the latter explained 4.0% of the variance and its occurrence probabilities continuously peaked around 500 mm and later decreased (Figure 1). The Jackknife test corroborated these results, indicating that precipitation seasonality and temperature annual range are the most important variables, which explain the distribution of buffelgrass. These tests also indicated that the explanatory power of the CNM is reduced if any of the bioclimatic variables used for calibration are excluded—i.e., the gain of the CNM that includes all the bioclimatic variables is lower than the gain of the models that exclude each variable (Figure 1b).
The geographical expression of the CNM under the current climate (Figure 2) indicated that 42.2% of the continental surface of Mexico comprises habitats that provide highly suitable climatic conditions for buffelgrass (i.e., map pixels with occurrence probabilities above 0.5 for the target species). Meanwhile, 45.6% of the country’s surface offers habitats with moderately suitable climates for the species (i.e., map pixels with occurrence probabilities between 0.1 and 0.5 for the target species). Only 12.3% of the surface of Mexico comprised climatically unsuitable habitats for the exotic species (i.e., map pixels with occurrence probabilities below 0.1 for the target species). An interactive map illustrating the current distribution of occurrence probabilities of the species estimated with the CNM is available in Supplementary Material 02 (Zenodo repository https://doi.org/10.5281/zenodo.6323654). The overlapping of this map with the climate units of Mexico indicated that warm habitats (highly warm, warm, or semi-warm climates in Table 1) with low rainfall (highly dry, dry, or semi-dry climates in Table 1) are highly susceptible to being invaded by buffelgrass. Conversely, this analysis indicated that the invasibility of habitats declines as temperature decreases and rainfall increases (temperate, semi-cold, cold climates, sub-humid, humid, or wet climates in Table 1).
The projection of the CNM under climate change scenarios predicted that the cover of highly suitable habitats for buffelgrass (i.e., pixels with occurrence probabilities above 0.5) will rise as the projections are further in the future and radiative forcing intensifies, as compared to the cover under the current climate (Figure 2; Table 2). Conversely, the availability of moderately suitable habitats (i.e., pixels with occurrence probabilities between 0.1 and 0.5) and unsuitable habitats for the species (i.e., pixels with occurrence probabilities below 0.1) was predicted to decrease for projections further in the future and for increased radiative forcing (Figure 2; Table 2). When these predictions were superimposed on their respective MESS maps (Figure 3), the overlapping between suitable and unsuitable habitats for buffelgrass was below 10% in most climate change scenarios (Table 2). The exceptions occurred in the period of 2081–2100 with the SSPs that had higher radiative forcing (7.0 and 8.5 W/m2), wherein the overlapping between the predicted suitable habitats and climatically unsuitable areas for buffelgrass were above 20% (Table 2). This indicates that the CNM generates highly reliable predictions for the future distribution of buffelgrass on most climate change scenarios expected in Mexico. Meanwhile, predictions performed by the end of this century at elevated radiative forcing levels must be taken with caution, since climate change may induce extremely unfavorable environmental conditions for this exotic species. Interactive maps illustrating the future distribution of buffelgrass occurrence probabilities estimated with the CNM are available in Supplementary Material 02 (Zenodo repository https://doi.org/10.5281/zenodo.6323654).
The future occurrence probabilities of buffelgrass estimated with the CNM at the ten-thousand random coordinates plotted on Mexico positively related with their current occurrence probabilities for all climate change scenarios (Table 3). The intercepts and slopes of all these linear regression functions significantly differed from the theoretical values (β0 = 0 and β1 = 1). The theoretical values represented the expected values from no effects of climate change on the distribution of buffelgrass (Table 3). In all climate change scenarios, the empirical regression curves dropped above the theoretical curve (Figure 4).

3. Discussion

The elevated AUC of the CNM for buffelgrass, calibrated with data from North and Central America, shows the estimates of the model with elevated accuracy, thereby indicating the distribution of occurrence probabilities of this exotic species in Mexico. This has been indicated for most species’ distribution models with AUC values above 0.9 [17,18,31,32,33]. The geographical projection of the CNM suggested that warm habitats under the current climate with limited precipitation are more susceptible to become invaded by buffelgrass (i.e., they have a high invasibility) than temperate, cold, and wet habitats of this country. Our findings concur with the generalized suggestion that after its introduction in new regions, buffelgrass mainly colonizes arid and semi-arid habitats [1]. Similar results have been reported with the application of a different software to model the distribution of buffelgrass [4]. A recently developed distribution model, focusing on southern USA and northern Mexico, indicated that the invasibility of habitats increases with rising annual mean temperature and decreasing annual rainfall [34]. However, these environmental variables were not included in our model, since they were spatially autocorrelated with the bioclimatic variables used to calibrate the CNM of buffelgrass. This exclusion of annual mean temperature and annual precipitation from our model could be a consequence of considering a larger geographical extent than that used in studies focused on particular fractions of the buffelgrass distribution range. The increases or reductions in the spatial extent of CNM can change the nature of the bioclimatic variables, thereby explaining the distribution of species [35].
In our study, the bioclimatic variables with higher explanatory power on the distribution of buffelgrass in Mexico were precipitation seasonality and temperature annual range. High probabilities of occurrence were estimated for precipitations from 120 mm to 160 mm and temperatures between 15 °C and 30 °C. These results suggest that buffelgrass would prefer habitats where rainfall concentrates in a particular period of the year. Since Mexico has tropical cyclones that generate summer storms frequently [36], it is feasible to propose that the elevated explanatory power of this bioclimatic variable in the CNM occurs because the species recruits in this rainy period. In addition, these results suggest buffelgrass would mainly invade habitats where temperature widely varies across the year, which is a typical climatic feature of arid and semi-arid Mexican ecosystems [37]. These suggestions concur with other findings, which have indicated that this exotic grass can get established at habitats where annual rainfall is below 300 mm and annual temperature ranges from −5 to 45 °C [1,2,38].
The projection of the CNM on climate change scenarios suggested that the distribution range of buffelgrass will continually expand during this century. This may be due to the increasing availability of suitable habitats for the species across Mexico, with larger expansions of the distribution range at a higher radiative forcing. In other words, the invasibility of habitats is predicted to increase with larger climatic changes. In contrast, de Albuquerque [27] reported that buffelgrass will decrease in the tropical zones of the Sonoran Desert and the lowest part of the Colorado River. Our predictions are supported by the results of the regression analyses conducted between future and current occurrence probabilities, which were estimated with the CNM. Such results indicated that the occurrence probabilities of this exotic plant will intensify across habitats of Mexico in all climate change scenarios (i.e., empirical regression curves drop above the theoretical curve with β0 = 0 and β1 = 1). Thus, the ecological niche of buffelgrass seems to coincide with the environmental conditions that predominate in warm drylands [1]. The predicted expansion of its distribution range is supported by the application of regional circulation models in Mexico. This indicates increases in temperature and decreases in rainfall due to climate change, which will increase aridity across this country [39,40].
An important caveat about the current and future invasibility of habitats reported in this study is that the CNM considers that buffelgrass can colonize all climatically suitable habitats available in Mexico. This is a common assumption for models estimating the distribution of species relying on climatic variables only [27,41,42]. Thus, the geographical expression of our CNM illustrates the potential distribution range that buffelgrass can reach in the absence of natural enemies and dispersal limitations. However, as these factors might also influence the establishment of species [43,44,45], they must be discussed to assess whether the current and future distribution ranges of buffelgrass estimated with the CNM are reliable.
In Mexico, buffelgrass is promoted as cattle forage, and in production fields these plants are commonly attacked by the pathogenic fungus Magnaporthe oryzae (Hebert) Barr (Magnaporthaceae), thereby reducing seed harvesting [46]. Given the extensive meadows this species has developed in human-disturbed sites of northern Mexico [47], it is feasible to propose that parasites and predators have negligible effects on naturally established populations of buffelgrass. It is important to note that buffelgrass seeds are dispersed by wind and water run-off, which are common dispersal mechanisms in successful invasive plants [48]. In addition, buffelgrass seeds are often dispersed in mud attached to vehicles, while they can also adhere to cattle fur and reach remote areas that are inaccessible to vehicles [49]. This variety of highly efficient dispersal mechanisms suggests that buffelgrass seeds can reach habitats located far away from parental populations, making it unfeasible to postulate that the species is affected by dispersal limitations. These lines of evidence suggest that both the current and future distribution ranges that were estimated with the CNM of buffelgrass in Mexico are reliable enough to be used for developing regional management and control programs for this exotic species.

4. Materials and Methods

Although our study focused on continental Mexico only, the occurrence data required for calibrating the CNM of buffelgrass were gathered across North and Central America to cover an extensive range of environmental conditions where the species is present, as this improves the predictive accuracy of the model. This information was obtained from the Global Biodiversity Information Facility (https://doi.org/10.15468/dl.acn5rp (accessed on 9 December 2019)). This online platform compiles the largest collection of corroborated plant records, even those located in regions where species are exotic [50]. We visualized the occurrence data in QGIS 3.12.0 (available at https://qgis.org (accessed on 15 August 2021)) and removed those points located within cities and towns to avoid presences due to human interventions (e.g., unintentional water supply) that might locate the species at sites where climate is not necessarily favorable for its development [17,18]. Climate data associated with these occurrence points were obtained from the WorldClim geodatabases (available at http://worldclim.org (accessed on 31 July 2021)), which interpolate temperature and rainfall information gathered between 1970 and 2000 to estimate the current worldwide values of 19 bioclimatic variables commonly used to explain species distribution [51]. WorldClim also downscales the predictions of CMIP6 general circulation models to estimate the values of these bioclimatic variables in climate change scenarios, but since the finest resolution at which these data are provided is 2.5 arc-minutes (~21 km2/pixel), the buffelgrass CNM was calibrated at this spatial scale.
Once the spatial scale of the CNM was defined, we performed a new depuration of buffelgrass presences to avoid including duplicated information within the same spatial unit, as this may lead to model overfitting [52,53]. For this, we drew a circular buffer of 2.5-km radius around each occurrence point and retained a single, randomly selected record when two or more buffers overlapped. This resulted in 292 occurrence points for the modeling procedure. Given CNM may overpredict species occurrence probabilities when redundant environmental variables are included in their calibration [22,54], we conducted Spearman rank-order tests among all pairwise combinations of bioclimatic variables (i.e., bioclimatic data associated with occurrence points). We ran these analyses in R 4.0 software (available at http://cran.r-project.org (accessed on 15 August 2021)) and, when correlation coefficients above 0.60 were found, we only retained the variable that had the highest number of correlations with others. Besides minimizing environmental redundancy, this allows for an explanation of the distribution of species with relevant bioclimatic variables only [55]. With these analyses we selected six bioclimatic variables, including temperature annual range, average temperature of the driest quarter of the year, average temperature of the warmest quarter of the year, precipitation seasonality, precipitation of the wettest quarter of the year, and precipitation of the coldest quarter of the year. The dataset with the occurrence points of buffelgrass and their associated values of bioclimatic variables, as well as the results of the correlation tests used to select bioclimatic variables for CNM, are available in Supplementary Material 01 (deposited in the Zenodo repository https://doi.org/10.5281/zenodo.6323654).
We calibrated the CNM of buffelgrass with MaxEnt 3.4 software [56] because its maximum entropy algorithm performs better than other computer routines (e.g., GARP) for estimating the occurrence probabilities of species when presence-only data are available, as was our case [22,54,56]. We instructed MaxEnt to randomly select 75% of the occurrence data to calibrate the CMN (training points), while the remaining 25% of the buffelgrass data were used to validate the resulting model (test points) with receiver operating characteristic curves. These curves indicate the fraction of test points correctly classified by the CNM, while the area under the curves (AUC) estimates the accuracy of the model to estimate the occurrence probabilities of the target species as a function of the bioclimatic variables used to calibrate it. AUC varies between one (1) and zero (0), where values close to 1 indicate the CNM properly estimates the occurrence probabilities of the species in the physical space and values below 0.5 indicate that the CNM cannot differentiate among suitable and unsuitable habitats for the species [56]. To prevent potential biases due to the random selection of training points, the CNM was calibrated 100 times using the bootstrap resampling procedure of MaxEnt. After that, we integrated the 100 outputs in a single CNM and averaged their AUC to assess the accuracy of the integrative model [20]. Furthermore, we used the jackknife test of MaxEnt to identify what bioclimatic variables have important individual effects on the distribution of buffelgrass. This procedure generates a single-variable model for each environmental variable used to calibrate the CNM to identify which of them make larger contributions to the model gain (goodness of fit of the full model). Meanwhile, the Jackknife test generates models that exclude each of these variables, while including the rest of them, and compare the gain of these models against the gain of the full model. In these comparisons, an environmental variable is assumed as relevant to explain the distribution of the target species if the model that excludes it has a lower gain than the full model [57].
We plotted the CNM of buffelgrass in QGIS using an environmental background with current values for the bioclimatic variables, which were used to calibrate it, and this resulted in a map whose pixels have occurrence probabilities ranging from zero (0) to one (1). For identifying which sites have an elevated likelihood of containing the species, we reclassified the pixels of this map in probability deciles and superimposed on it the occurrence points of buffelgrass for counting the number of true presences at each decile. Since only 1.3% of buffelgrass presences dropped in pixels with occurrence probabilities below 0.1, we assumed the pixels in the lowest probability decile (0.0–0.1) constituted habitats climatically unsuitable for the species. The accumulated fraction of true presences increased to 18.9% in pixels with occurrence probabilities between 0.1 and 0.5; however, 79.8% of buffelgrass presences were concentrated in pixels with occurrence probabilities above 0.5. The number of buffelgrass occurrences that dropped at each probability decile can be consulted in Supplementary Material 01 (deposited in the Zenodo repository https://doi.org/10.5281/zenodo.6323655). Considering these values, we assumed that pixels with occurrence probabilities between 0.1 and 0.5 are habitats that offer moderately suitable climatic conditions for buffelgrass. Meanwhile, pixels with occurrence probabilities above 0.5 are highly suitable habitats for the development of the species. After that, we constrained this occurrence probability map to the continental surface of Mexico and superimposed it on a map containing the climate units of this country (available at https://www.inegi.org.mx (accessed on 9 December 2019)). According to the latter, Mexico can be split into 21 climate units (Table 1). Within each climate unit, we computed the degree of overlapping among highly suitable, moderately suitable, and unsuitable habitats to identify the most susceptible to become invaded by buffelgrass under the current climatic conditions (i.e., we assessed the invasibility of these habitats in accordance to Lonsdale [28]).
The future distribution of buffelgrass was estimated using the predictions of the Canadian Earth System Model version 5 (CanESM5) as a climatic background. For that, we chose CanESM5 because it has a higher equilibrium climate sensitivity and transient climate response than other CMIP6 models included in WorldClim [58]. As these parameters estimate the thermal sensitivity of the Earth’s system in the face of rising concentrations of greenhouse gases, the elevated values of CanESM5 make its predictions highly realistic [59,60]. Given the climate predictions from the CanESM5 cover for 20-year intervals until the end of this century (2041–2060, 2061–2080, and 2081–2100), we estimated the future distribution of buffelgrass for all of these periods considering the radiative forcing associated with the four SSP (2.6, 4.5, 7.0, and 8.5 W/m2). Then, we projected the CNM calibrated at each bootstrap run on all climate change scenarios (three periods × 4 SSPs = 12 climate change scenarios) and then integrated the output of these projections in a single prediction for each climate change scenario. These predictions were visualized as maps of occurrence probabilities whose pixels were classified using the aforementioned criteria. To assess the reliability of these predictions, we ran a multivariate environmental similarity surface (MESS) analysis each time the CNM was projected on a climate change scenario. The MESS analyses compare the current and future values of the bioclimatic variables used to calibrate the CNM across the target region. It differentiates among habitats that will become unsuitable for the species due to climate change (pixels with negative MESS values) and habitats that will concur with the climatic requirements of the species (pixels with positive MESS values) [27]. For each climate change scenario, we integrated the outputs of the MESS analyses that resulted from the 100 bootstrap runs in a single MESS output and expressed it as a map with pixels categorized as positive or negative. On these MESS maps, we superimposed the occurrence probabilities of buffelgrass predicted at the respective climate change scenario and calculated the percent overlapping of moderately and highly suitable habitats (i.e., pixels with occurrence probabilities above 0.1) with areas classified as unsuitable for the species in the MESS maps (i.e., areas with negative MESS values). For this procedure, we assumed that the reliability of the CNM to predict the future distribution of the exotic species increases as this percentage decreases [18].
To determine whether climate change will substantially affect the distribution of buffelgrass in Mexico, we compared the future distribution of occurrence probabilities for the species against the current distribution of occurrence probabilities. For this, we generated ten-thousand random geographic coordinates across the continental Mexico and plotted them on the current and future occurrence probabilities maps estimated with the CNM. After that, we extracted the probability values associated to each of these points and performed simple linear regression analyses using the future occurrence probabilities as response variables and the current occurrence probabilities as predictive variables. These analyses were conducted separately for each climate change scenario in R 4.0. software. In all the cases, we obtained a linear regression function of the type Pf = β0 + β1 Pc, where Pf is the future occurrence probability of the species at a given point of the geographical space and Pc is the current occurrence probability of the species at the same geographical coordinates. Meanwhile, β0 and β1 represent the intercept and slope of the regression function, respectively, which were estimated with the method of least squares [61]. This approach allows for proposing the following null hypothesis: if climate change has no effects on the distribution of buffelgrass, its current and future occurrence probabilities should be similar across the geographical space, and, therefore, the parameters of the empirical regression function should not differ from their theoretical values (β0 = 0 and β1 = 1). Alternatively, if these regression parameters differ from their theoretical values (i.e., β0 ≠ 0 and/or β1 ≠ 1), the former null hypothesis must be rejected. We used Student’s t-test to determine whether empirical and theoretical values of β0 and β1 statistically differ.

5. Conclusions

The CNM calibrated in this study indicates that the warm drylands of Mexico constitute habitats that are highly susceptible to being invaded by buffelgrass. In addition, as the advance of climate change will increase aridity across this country, the model predicts that buffelgrass populations will continually expand during this century. These analyses clearly suggest that Mexico is facing a biological invasion by buffelgrass and given the aggressive invasiveness of the species, it may constitute a threat for native biodiversity. Therefore, regional management and control programs should be implemented to prevent the expansion of buffelgrass populations. These programs must also consider that habitat invasibility will increase in the future, which could make it harder to control this biological invasion.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants11091160/s1, Supplementary material 01, excel file that contains the information that supports the analyses performed in the study, Supplementary material 02, Keyhole Markup Language file with maps that can be visualized in Google Earth.

Author Contributions

Conceptualization, P.S.-C., E.I.B. and N.S.H.-Q.; Data curation, P.S.-C., R.C.-L. and A.Á.-H.; Formal analysis, E.I.B. and F.V.-G.; Methodology, E.I.B. and N.S.H.-Q.; Project administration, N.S.H.-Q. and E.I.B.; Software, J.A.P.-A. and A.P.-A.; Supervision, A.P.-A. and R.C.-L.; Validation, E.I.B. and A.Á.-H.; Writing–original draft, P.S.-C., N.S.H.-Q. and E.I.B.; Writing–review & editing, F.V.-G., R.C.-L., A.Á.-H., A.P.-A. and J.A.P.-A. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by Programa de Mejoramiento al Profesorado (PROMEP), grant number UACH-EXB-250 to N.S.H.-Q. and Fondo Sectorial de Investigación Ambiental SEMARNAT-CONACYT, grant FSSEMARNAT01-C-2018-1-A3-S-80837 to E.I.B.

Data Availability Statement

All the data presented in this study are available in the article and in Supplementary Materials.

Acknowledgments

P.S.-C. Thank the scholarship provided by Consejo Nacional de Ciencia y Tecnología (CONA-CYT-México).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Marshall, V.M.; Lewis, M.M.; Ostendorf, B. Buffel grass (Cenchrus ciliaris) as an invader and threat to biodiversity in arid environments: A review. J. Arid Environ. 2012, 78, 1–12. [Google Scholar] [CrossRef]
  2. de la Barrera, E. Recent invasion of buffel grass (Cenchrus ciliaris) of a natural protected area from the southern Sonoran Desert. Rev. Mex. Biodivers. 2008, 79, 385–392. [Google Scholar]
  3. Ibarra-Flores, F.A.; Marín-Rivera, M.H.; Quero-Carrillo, A.R. Pasto buffel Cenchrus ciliaris L. Link. (Pennisetum ciliare Lam.). In Gramíneas Introducidas: Importancia e Impacto En Ecosistemas Ganaderos, 1st ed.; Quero-Carrillo, A.R., Ed.; Colegio de Postgraduados: Texcoco, México, 2013; Volume 1, pp. 73–102. [Google Scholar]
  4. Arriaga, L.; Castellanos, V.A.E.; Moreno, E.; Alarcón, J. Potential ecological distribution of alien invasive species and risk assessment: A case study of buffel grass in arid regions of Mexico. Conserv. Biol. 2004, 18, 1504–1514. [Google Scholar] [CrossRef]
  5. Brenner, J.C.; Kanda, L.L. Buffelgrass (Pennisetum ciliare) invades lands surrounding cultivated pastures in Sonora, Mexico. Invasive Plant. Sci. Manag. 2013, 6, 187–195. [Google Scholar] [CrossRef]
  6. Blanco, L.J.; Ferrando, C.A.; Biurrun, F.N.; Orionte, E.L.; Namur, P.; Recalde, D.J.; Berone, G.D. Vegetation responses to roller chopping and buffelgrass seeding in Argentina. Rangel. Ecol. Manag. 2005, 58, 219–224. [Google Scholar] [CrossRef]
  7. Saucedo-Monarque, E.; García-Moya, E.; Castellanos-Villegas, A.E.; Flores-Flores, J.L. La riqueza, una variable de respuesta de la vegetación a la introducción del zacate buffel. Agrociencia 1997, 31, 83–90. [Google Scholar]
  8. Williams, D.G.; Baruch, Z. African grass invasion in the Americas: Ecosystem consequences and the role of ecophysiology. Biol. Invasions 2000, 2, 123–140. [Google Scholar] [CrossRef]
  9. Chaffin, B.C.; Garmestani, A.S.; Angeler, D.G.; Herrmann, D.L.; Stow, C.A.; Nyström, M.; Sendzimir, J.; Hopton, M.E.; Kolasa, J.; Allen, C.R. Biological invasions, ecological resilience and adaptive governance. Environ. Manag. 2016, 183, 399–407. [Google Scholar] [CrossRef] [Green Version]
  10. Seipel, T.; Rew, L.J.; Taylor, K.T.; Maxwell, B.D. Disturbance type influences plant community resilience and resistance to Bromus tectorum invasion in the sagebrush steppe. Appl. Veg. Sci. 2018, 21, 385–394. [Google Scholar] [CrossRef] [Green Version]
  11. Lien, A.M.; Baldwin, E.; Franklin, K. Collective action and invasive species governance in Southern Arizona. Range Ecol. Manag. 2021, 74, 151–164. [Google Scholar] [CrossRef]
  12. Kumar, R.P.; Singh, J.S. Invasive alien plant species: Their impact on environmental, ecosystem services and human health. Ecol. Indic. 2020, 111, 106020. [Google Scholar] [CrossRef]
  13. Brooks, M.L.; Chambers, J.C. Resistance to invasion and resilience to fire in desert shrublands of North America. Rangeland Ecol. Manag. 2011, 64, 431–438. [Google Scholar] [CrossRef]
  14. Diario Oficial de la Federación. Available online: https://www.dof.gob.mx/nota_detalle.php?codigo=5464456&fecha=07/12/2016 (accessed on 31 March 2022).
  15. Método de Evaluación Rápida de Invasibidad (MERI) para Especies Exóticas en México. Cenchrus ciliaris L. Available online: http://sivicoff.cnf.gob.mx/ContenidoPublico/MenuPrincipal/07Fichas%20tecnicas_OK/02Fichas%20tecnicas/Fichas%20t%C3%A9cnicas%20SEMARNAT_especies%20ex%C3%B3ticas/Cenchrus%20ciliaris.pdf (accessed on 31 March 2022).
  16. Hui, C.; Krug, R.M.; Richardson, D.M. Modelling spread in invasion ecology: A synthesis. In Fifty Years of Invasion Ecology: The Legacy of Charles Elton, 1st ed.; Richardson, D.M., Ed.; Wiley: Oxford, UK, 2010; Volume 1, pp. 329–343. [Google Scholar] [CrossRef]
  17. Ramírez-Albores, J.E.; Bustamante, R.O.; Badano, E.I. Improved predictions of the geographic distribution of invasive plants using climatic niche models. PLoS ONE 2016, 11, e0156029. [Google Scholar] [CrossRef]
  18. Guerra-Coss, F.A.; Badano, E.I.; Cedillo-Rodríguez, I.E.; Ramírez-Albores, J.E.; Flores, J.; Barragán-Torres, F.; Flores-Cano, J.A. Modelling and validation of the spatial distribution of suitable habitats for the recruitment of invasive plants on climate change scenarios: An approach from the regeneration niche. Sci. Total Environ. 2021, 777, 146007. [Google Scholar] [CrossRef]
  19. Guisan, A.; Thuiller, W. Predicting species distribution: Offering more than simple habitat models. Ecol. Lett. 2005, 8, 993–1009. [Google Scholar] [CrossRef]
  20. Phillips, S.J.; Dudík, M. Modeling of species distributions with Maxent: New extensions and a comprehensive evaluation. Ecography 2008, 31, 161–175. [Google Scholar] [CrossRef]
  21. Peterson, A.T.; Soberón, J.; Pearson, R.G.; Anderson, R.P.; Martínez-Meyer, E.; Nakamura, M.; Araújo, M.B. Ecological Niches and Geographic Distributions, 1st ed.; Princeton University Press: Princeton, NJ, USA, 2011; Volume 328, pp. 97–137. [Google Scholar]
  22. Elith, J.; Phillips, S.J.; Hastie, T.; Dudík, M.; Chee, Y.E.; Yates, C.J. A statistical explanation of MaxEnt for ecologists. Divers. Distrib. 2011, 17, 43–57. [Google Scholar] [CrossRef]
  23. Forster, P.M.; Maycock, A.C.; McKenna, C.M.; Smith, C.J. Latest climate models confirm need for urgent mitigation. Nat. Clim. Chang. 2019, 10, 7–10. [Google Scholar] [CrossRef]
  24. Cook, B.I.; Mankin, J.S.; Marvel, K.; Williams, A.P.; Smerdon, J.E.; Anchukaitis, K.J. Twenty-first century drought projections in the CMIP6 forcing scenarios. Earths Future 2020, 8, e2019EF001461. [Google Scholar] [CrossRef] [Green Version]
  25. O’Neill, B.C.; Kriegler, E.; Riahi, K.; Ebi, K.L.; Hallegatte, S.; Carter, T.R.; Mathur, R.; van Vuuren, D.P. A new scenario framework for climate change research: The concept of shared socioeconomic pathways. Clim. Chang. 2014, 122, 387–400. [Google Scholar] [CrossRef] [Green Version]
  26. Riahi, K.; van Vuuren, D.P.; Kriegler, E.; Edmonds, J.; O’Neill, B.C.; Fujimori, S.; Bauer, N.; Calvin, K.; Dellink, R.; Fricko, O.; et al. The Shared Socioeconomic Pathways and their energy, land use, and greenhouse gas emissions implications: An overview. Glob. Environ. Chang. 2017, 42, 153–168. [Google Scholar] [CrossRef] [Green Version]
  27. Elith, J.; Kearney, M.; Phillips, S. The art of modelling range-shifting species. Methods Ecol. Evol. 2010, 1, 330–342. [Google Scholar] [CrossRef]
  28. Lonsdale, W.M. Global patterns of plant invasions and the concept of invasibility. Ecology 1999, 80, 1522–1536. [Google Scholar] [CrossRef]
  29. Bradley, B.A.; Blumenthal, D.M.; Wilcove, D.S.; Ziska, L.H. Predicting plant invasions in an era of global change. Trends Ecol. Evol. 2010, 25, 310–318. [Google Scholar] [CrossRef] [PubMed]
  30. Petitpierre, B.; Kueffer, C.; Broennimann, O.; Randin, C.; Daehler, C.; Guisan, A. Climatic niche shifts are rare among terrestrial plant invaders. Science 2012, 335, 344–1348. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Araújo, M.B.; Pearson, R.G.; Thuiller, W.; Erhard, M. Validation of species–climate impact models under climate change. Glob. Chang. Biol. 2005, 11, 1504–1513. [Google Scholar] [CrossRef] [Green Version]
  32. Peña-Gómez, F.T.; Guerrero, P.C.; Bizama, G.; Duarte, M.; Bustamante, R.O. Climatic niche conservatism and biogeographical non-equilibrium in Eschscholzia californica (Papaveraceae), an invasive plant in the Chilean Mediterranean region. PLoS ONE 2014, 9, e105025. [Google Scholar] [CrossRef] [Green Version]
  33. Wan, J.Z.; Wang, C.J.; Yu, F.H. Effects of occurrence record number, environmental variable number, and spatial scales on MaxEnt distribution modelling for invasive plants. Biology 2019, 74, 757–766. [Google Scholar] [CrossRef]
  34. Suzart de Albuquerque, F.; Macías-Rodríguez, M.Á.; Búrquez, A.; Astudillo-Scalia, Y. Climate change and the potential expansion of buffelgrass (Cenchrus ciliaris L., Poaceae) in biotic communities of Southwest United States and northern Mexico. Biol. Invasions 2019, 21, 3335–3347. [Google Scholar] [CrossRef]
  35. Stankowski, P.A.; Parker, W.H. Species distribution modelling: Does one size fit all? A phytogeographic analysis of Salix in Ontario. Ecol. Modell. 2010, 221, 1655–1664. [Google Scholar] [CrossRef]
  36. Agustín Breña-Naranjo, J.; Pedrozo-Acuña, A.; Pozos-Estrada, O.; Jiménez-López, S.A.; López-López, M.R. The contribution of tropical cyclones to rainfall in Mexico. Phys. Chem. Earth Parts A/B/C 2015, 83–84, 111–122. [Google Scholar] [CrossRef]
  37. González-Medrano, F. Las Zonas Áridas y Semiáridas de México y Su Vegetación, 1st ed.; Instituto Nacional de Ecología: México City, Mexico, 2012; pp. 9–159. [Google Scholar]
  38. De la Barrera, E.; Castellanos, A.E. High temperature effects on gas exchange for the invasive buffel grass (Pennisetum ciliare [L.] Link). Weed Biol. Manag. 2007, 7, 128–131. [Google Scholar] [CrossRef]
  39. Karmalkar, A.V.; Bradley, R.S.; Diaz, H.F. Climate change in Central America and Mexico: Regional climate model validation and climate change projections. Clim. Dyn. 2011, 37, 605. [Google Scholar] [CrossRef] [Green Version]
  40. Colorado-Ruiz, G.; Cavazos, T.; Salinas, J.A.; De Grau, P.; Ayala, R. Climate change projections from Coupled Model Intercomparison Project phase 5 multi-model weighted ensembles for Mexico, the North American monsoon, and the mid-summer drought region. Int. J. Climatol. 2018, 38, 5699–5716. [Google Scholar] [CrossRef]
  41. Araújo, M.B.; Pearson, R.G. Equilibrium of species’ distributions with climate. Ecography 2005, 28, 693–695. [Google Scholar] [CrossRef]
  42. Václavík, T.; Meentemeyer, R.K. Equilibrium or not? Modelling potential distribution of invasive species in different stages of invasion. Divers. Distrib. 2012, 18, 73–83. [Google Scholar] [CrossRef]
  43. Sol, D. Do successful invaders exist? Pre-adaptations to novel environments in terrestrial vertebrates. In Biological Invasions, 1st ed.; Nentwig, W., Ed.; Springer: Berlin, Germany, 2007; Volume 193, pp. 127–144. [Google Scholar] [CrossRef]
  44. Coutts, S.R.; van Klinken, R.D.; Yokomizo, H.; Buckley, Y.M. What are the key drivers of spread in invasive plants: Dispersal, demography or landscape: And how can we use this knowledge to aid management? Biol. Invasions 2011, 13, 1649–1661. [Google Scholar] [CrossRef]
  45. Meijer, K.; Schilthuizen, M.; Beukeboom, L.; Smit, C. A review and meta-analysis of the enemy release hypothesis in plant-herbivorous insect systems. PeerJ 2016, 4, e2778. [Google Scholar] [CrossRef] [Green Version]
  46. Ibarra-Flores, F.A.; Martin-Rivera, M.H.; Moreno-Medina, S.; Ibarra-Martín, F.A.; Retes-López, R. Impacto económico asociado con los daños del tizón foliar en la producción de semilla del zacate buffel en el centro de Sonora, México. Rev. Mex. Agronegocios 2019, 45, 313–324. [Google Scholar]
  47. Morales-Romero, D.; Lopez-Garcia, H.; Martinez-Rodriguez, J.; Molina-Freaner, F. Documenting a plant invasion: The influence of land use on buffelgrass invasion along roadsides in Sonora, Mexico. J. Arid Environ. 2019, 164, 53–59. [Google Scholar] [CrossRef]
  48. Pyšek, P.; Richardson, D.M. Traits associated with invasiveness in alien plants: Where do we stand? In Biological Invasions, 1st ed.; Nentwig, W., Ed.; Springer: Berlin, Germany, 2007; Volume 193, pp. 97–125. [Google Scholar] [CrossRef]
  49. Stevens, J.; Falk, D.A. Can buffelgrass invasions be controlled in the American Southwest? Using invasion ecology theory to understand buffelgrass success and develop comprehensive restoration and management. Ecol. Restor. 2009, 27, 417–427. [Google Scholar] [CrossRef]
  50. James, S.A.; Soltis, P.S.; Belbin, L.; Chapman, A.D.; Nelson, G.; Paul, D.L.; Collins, M. Herbarium data: Global biodiversity and societal botanical needs for novel research. Appl. Plant. Sci. 2018, 6, e1024. [Google Scholar] [CrossRef]
  51. Fick, S.E.; Hijmans, R.J. WorldClim 2: New 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 2017, 37, 4302–4315. [Google Scholar] [CrossRef]
  52. Barry, S.; Elith, J.; Barry, S.; Elith, J. Error and uncertainty in habitat models. J. Appl. Ecol. 2006, 43, 413–423. [Google Scholar] [CrossRef]
  53. Elith, J.; Graham, C.H.; Anderson, R.P.; Dudík, M.; Ferrier, S.; Guisan, A.; Hijmans, R.J.; Huettmann, F.; Leathwick, J.R.; Lehmann, A.; et al. Novel methods improve prediction of species’ distributions from occurrence data. Ecography 2006, 29, 129–151. [Google Scholar] [CrossRef] [Green Version]
  54. Merow, C.; Smith, M.J.; Silander, J.A. A practical guide to MaxEnt for modeling species’ distributions: What it does, and why inputs and settings matter. Ecography 2013, 36, 1058–1069. [Google Scholar] [CrossRef]
  55. Beaumont, L.J.; Hughes, L.; Poulsen, M. Predicting species distributions: Use of climatic parameters in BIOCLIM and its impact on predictions of species’ current and future distributions. Ecol. Modell. 2005, 186, 251–270. [Google Scholar] [CrossRef]
  56. 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]
  57. 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]
  58. Swart, N.C.; Cole, J.N.S.; Kharin, V.V.; Lazare, M.; Scinocca, J.F.; Gillett, N.P.; Anstey, J.; Arora, V.; Christian, J.R.; Hanna, S.; et al. The Canadian Earth System Model version 5 (CanESM5.0.3). Geosci. Model. Dev. 2019, 12, 4823–4873. [Google Scholar] [CrossRef] [Green Version]
  59. Meehl, G.A.; Senior, C.A.; Eyring, V.; Flato, G.; Lamarque, J.F.; Stouffer, R.J.; Taylor, K.E.; Schlund, M. Context for interpreting equilibrium climate sensitivity and transient climate response from the CMIP6 Earth system models. Sci. Adv. 2020, 6, eaba1981. [Google Scholar] [CrossRef]
  60. Zelinka, M.D.; Myers, T.A.; McCoy, D.T.; Po-Chedley, S.; Caldwell, P.M.; Ceppi, P.; Klein, S.A.; Taylor, K.E. Causes of higher climate sensitivity in CMIP6 models. Geophys. Res. Lett. 2020, 47, e2019GL085782. [Google Scholar] [CrossRef] [Green Version]
  61. Kutner, M.H.; Nachtsheim, C.J.; Neter, J.; Li, W. Applied Linear Statistical Models, 5th ed.; McGraw-Hill Irwin: New York, NY, USA, 2005; pp. 2–39. [Google Scholar]
Plants 11 01160 g001 550
Figure 1. Response curves for each variable and Jackknife test. (a) The curves show how the predicted probability of occurrence changes as each environmental variable varies. (b) The Jackknife test indicates in blue the gain of the model from keeping only that variable, in red it shows the gain without the variable, and in green it shows the gain of the model with all the variables.
Figure 1. Response curves for each variable and Jackknife test. (a) The curves show how the predicted probability of occurrence changes as each environmental variable varies. (b) The Jackknife test indicates in blue the gain of the model from keeping only that variable, in red it shows the gain without the variable, and in green it shows the gain of the model with all the variables.
Plants 11 01160 g001
Plants 11 01160 g002 550
Figure 2. Distribution of suitable and unsuitable habitats for buffelgrass under the current and climate change scenarios. Occurrence probabilities below 0.1 are considered unsuitable, values between 0.1 to 0.5 offer moderate conditions, and above 0.5 are highly suitable habitats for the development of the species.
Figure 2. Distribution of suitable and unsuitable habitats for buffelgrass under the current and climate change scenarios. Occurrence probabilities below 0.1 are considered unsuitable, values between 0.1 to 0.5 offer moderate conditions, and above 0.5 are highly suitable habitats for the development of the species.
Plants 11 01160 g002
Plants 11 01160 g003 550
Figure 3. Analysis of multivariate environmental similarity surface (MESS) for each climate change scenario. Negative values indicate unsuitable habitat for the species because of climate change and positive values specify habitats that will concur with the climatic requirements of the species.
Figure 3. Analysis of multivariate environmental similarity surface (MESS) for each climate change scenario. Negative values indicate unsuitable habitat for the species because of climate change and positive values specify habitats that will concur with the climatic requirements of the species.
Plants 11 01160 g003
Plants 11 01160 g004 550
Figure 4. Graphs of the linear regressions analysis between current climate and climate change scenarios for the periods of 2041–2060, 2061–2080, and 2081–2100 for all the shared socioeconomic pathways. The color blue indicates the theoretical curve, and the color red shows the linear regression functions for each scenario.
Figure 4. Graphs of the linear regressions analysis between current climate and climate change scenarios for the periods of 2041–2060, 2061–2080, and 2081–2100 for all the shared socioeconomic pathways. The color blue indicates the theoretical curve, and the color red shows the linear regression functions for each scenario.
Plants 11 01160 g004
Table
Table 1. Results of the overlapping between the occurrence probability map of buffelgrass under the current climate and the map of the climate units of Mexico. For each climate unit, the table indicates the percent cover of habitats with highly suitable climatic conditions (occurrence probabilities above 0.5), moderately suitable climatic conditions (occurrence probabilities between 0.1 and 0.5), and unsuitable climatic conditions (occurrence probabilities below 0.1) for buffelgrass.
Table 1. Results of the overlapping between the occurrence probability map of buffelgrass under the current climate and the map of the climate units of Mexico. For each climate unit, the table indicates the percent cover of habitats with highly suitable climatic conditions (occurrence probabilities above 0.5), moderately suitable climatic conditions (occurrence probabilities between 0.1 and 0.5), and unsuitable climatic conditions (occurrence probabilities below 0.1) for buffelgrass.
Climate UnitHighly Suitable
Habitats		Moderately
Suitable Habitats		Unsuitable
Habitats
Cold0%8%92%
Semi-cold semi-dry 0%48%52%
Semi-cold sub-humid0%31%69%
Temperate dry 30%65%5%
Temperate highly dry3%72%26%
Temperate humid 1%54%45%
Temperate semi-dry 45%53%3%
Temperate sub-humid 37%54%9%
Semi-warm dry 66%34%0%
Semi-warm highly dry51%47%2%
Semi-warm humid 10%43%46%
Semi-warm semi-dry 62%38%0%
Semi-warm sub-humid 73%27%0%
Warm dry89%11%0%
Warm highly dry84%15%2%
Warm semi-dry 94%6%0%
Warm sub-humid37%60%2%
Warm wet0%14%86%
Highly warm highly dry 100%0%0%
Highly warm semi-dry 85%15%0%
Highly warm dry 91%9%0%
Table
Table 2. Percent cover of highly suitable habitats (pixels with occurrence probabilities above 0.5), moderately suitable habitats (pixels with occurrence probabilities between 0.1 and 0.5), and unsuitable habitats (pixels with occurrence probabilities below 0.1) for buffelgrass estimated with the CNM under the current climate and the climate change scenarios predicted by CanESM5 on 20-year intervals considering the radiative forcing associated with the different SSPs. The last column of the table shows the percent overlapping of habitats predicted as moderately and highly suitable for buffelgrass (i.e., occurrence probabilities above 0.1) with areas containing negative MESS values at each climate change scenario.
Table 2. Percent cover of highly suitable habitats (pixels with occurrence probabilities above 0.5), moderately suitable habitats (pixels with occurrence probabilities between 0.1 and 0.5), and unsuitable habitats (pixels with occurrence probabilities below 0.1) for buffelgrass estimated with the CNM under the current climate and the climate change scenarios predicted by CanESM5 on 20-year intervals considering the radiative forcing associated with the different SSPs. The last column of the table shows the percent overlapping of habitats predicted as moderately and highly suitable for buffelgrass (i.e., occurrence probabilities above 0.1) with areas containing negative MESS values at each climate change scenario.
Climate ScenarioHighly Suitable
Habitats		Moderately
Suitable Habitats		Unsuitable
Habitats		Overlapping with
Negative MESS Areas
Current climate42.2%45.6%12.3%–
2041–2060SSP2.6SSP2.652.4%39.2%8.4%
SSP4.5SSP4.554.5%37.6%7.9%
SSP7.0SSP7.056.8%35.8%7.4%
SSP8.5SSP8.558.1%34.6%7.3%
2061–2080SSP2.6SSP2.652.1%39.1%8.7%
SSP4.5SSP4.557.2%35.4%7.5%
SSP7.0SSP7.062.2%31.0%6.9%
SSP8.5SSP8.564.2%29.3%6.6%
2081–2100SSP2.6SSP2.651.4%39.8%8.8%
SSP4.5SSP4.557.9%34.6%7.5%
SSP7.0SSP7.068.1%25.9%6.1%
SSP8.5SSP8.572.0%22.7%5.2%
Table
Table 3. Result of the regression analyses conducted between current and future occurrence probabilities of buffelgrass across Mexico. Future occurrence probabilities were estimated by projecting the CNM on the current climate and the climate change scenarios predicted by CanESM5 on 20-year intervals considering the radiative forcing associated with the different SSPs. The table also provides the linear regression function that resulted from these analyses, where Pf is the future occurrence probability of the species and Pc is the current occurrence probability of the species.
Table 3. Result of the regression analyses conducted between current and future occurrence probabilities of buffelgrass across Mexico. Future occurrence probabilities were estimated by projecting the CNM on the current climate and the climate change scenarios predicted by CanESM5 on 20-year intervals considering the radiative forcing associated with the different SSPs. The table also provides the linear regression function that resulted from these analyses, where Pf is the future occurrence probability of the species and Pc is the current occurrence probability of the species.
Climate ScenarioResults of
Regression Analyses		Empirical Linear
Regression Function
2041–2060SSP2.6F(1, 9998) = 80,688.383, p < 0.001, R2 = 0.890 P f = 0.069 * + 1.032 * Pc
SSP4.5F(1, 9998) = 59,827.316, p < 0.001, R2 = 0.857Pf = 0.090 * + 1.025 * Pc
SSP7.0F(1, 9998) = 43,507.023, p < 0.001, R2 = 0.813Pf = 0.119 * + 1.013 * Pc
SSP8.5F(1, 9998) = 40,987.125, p < 0.001, R2 = 0.804Pf = 0.132 * + 0.988 * Pc
2061–2080SSP2.6F(1, 9998) = 80,439.742, p < 0.001, R2 = 0.889Pf = 0.063 * + 1.036 * Pc
SSP4.5F(1, 9998) = 41,775.102, p < 0.001, R2 = 0.807Pf = 0.121 * + 1.010 * Pc
SSP7.0F(1, 9998) = 24,058.211, p < 0.001, R2 = 0.706Pf = 0.198 * + 0.931 * Pc
SSP8.5F(1, 9998) = 19,815.479, p < 0.001, R2 = 0.665Pf = 0.234 * + 0.884 * Pc
2081–2100SSP2.6F(1, 9998) = 84,131.766, p < 0.001, R2 = 0.894Pf = 0.058 * + 1.038 * Pc
SSP4.5F(1, 9998) = 36,641.332, p < 0.001, R2 = 0.786Pf = 0.133 * + 0.981 * Pc
SSP7.0F(1, 9998) = 13,323.927, p < 0.001, R2 = 0.571Pf = 0.293 * + 0.804 * Pc
SSP8.5F(1, 9998) = 10,316.182, p < 0.001, R2 = 0.508Pf = 0.344 * + 0.738 * Pc
* The empirical value statistically differs from the theoretical value (t-tests critical α = 0.05).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Siller-Clavel, P.; Badano, E.I.; Villarreal-Guerrero, F.; Prieto-Amparán, J.A.; Pinedo-Alvarez, A.; Corrales-Lerma, R.; Álvarez-Holguín, A.; Hernández-Quiroz, N.S. Distribution Patterns of Invasive Buffelgrass (Cenchrus ciliaris) in Mexico Estimated with Climate Niche Models under the Current and Future Climate. Plants 2022, 11, 1160. https://doi.org/10.3390/plants11091160

AMA Style

Siller-Clavel P, Badano EI, Villarreal-Guerrero F, Prieto-Amparán JA, Pinedo-Alvarez A, Corrales-Lerma R, Álvarez-Holguín A, Hernández-Quiroz NS. Distribution Patterns of Invasive Buffelgrass (Cenchrus ciliaris) in Mexico Estimated with Climate Niche Models under the Current and Future Climate. Plants. 2022; 11(9):1160. https://doi.org/10.3390/plants11091160

Chicago/Turabian Style

Siller-Clavel, Pablo, Ernesto I. Badano, Federico Villarreal-Guerrero, Jesús A. Prieto-Amparán, Alfredo Pinedo-Alvarez, Raúl Corrales-Lerma, Alan Álvarez-Holguín, and Nathalie S. Hernández-Quiroz. 2022. "Distribution Patterns of Invasive Buffelgrass (Cenchrus ciliaris) in Mexico Estimated with Climate Niche Models under the Current and Future Climate" Plants 11, no. 9: 1160. https://doi.org/10.3390/plants11091160

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