Simulation of air temperature and their influence on the potential distribution of Myracrodruon urundeuva, Copernicia prunifera and Cereus jamacaru in the Caatinga

Areas in the process of aridification in Caatinga phytogeographic domain in northeastern Brazil increase every year due to human intervention and increase in air temperature. The identification of long-term patterns and air temperature trends in the phytogeographic domain can express climate variability as well as a new phase of adaptation to some plant species. The temperature series from 1951 to 2018 obtained from the National Centers for Environmental Prediction data sets in four conservation areas with native vegetation, located in the North (A1), East (A2), South (A3) and West (A4) regions of this phytogeographic domain, show an increase in temperature between 0.5 and 1.4 °C over the 68-year period with the highest warming occurring in the months of March, April and May. The Maxent model is used to identify the influence of this increase on the presence potential of three species in the Caatinga, Myracrodruon urundeuva (aroeira), Copernicia prunifera (palmeira) and Cereus jamacaru DC (cactus) in the future time interval of 2041 to 2060, considering IPCC projected climate changes. The results show that climate change can lead to a reduction as well as redistribution of the potential areas of occurrence of the three species. Notable changes are: in the case of Carnauba, the high potential area reduces from 25.3% in the present state to 19.6% in 2050, and potential area for Aroeira diminishes in central Bahia and increases in Rio Grande do Norte. The projected changes for all three species are discussed.


Introduction
The Caatinga phytogeographic domain is the only type of plant formation that is exclusively Brazilian [4]. Earlier studies indicate that this phytogeographic domain has a great wealth of habitats and species, many of which are exclusive to the eastern portion of Brazil [39]. Many of the Caatinga species provide essential ecosystem services for the population, such as the development of phytopharmaceuticals [46] and biodiversity. This phytogeographic domain covers almost 10% of the Brazilian territory with a high rate of endemism. The Caatinga is among the most altered phytogeographic domains, almost 50% of the alteration is anthropogenic through habitat destruction and fragmentation.
According to Seddon et al. [52], the Caatinga showed high sensitivity to climate variability in the last 14 years, mainly associated with variations in water availability and increased air temperature. The increase in global air temperature is attributed to the increases in the concentrations of greenhouse gases (GHG) in the atmosphere. Although greenhouse effect is natural and essential for the maintenance of life on the planet, the increase in GHG concentrations elevates the mean temperature of the planet.
The impacts of new climate conditions due to changes in atmospheric composition were investigated through projections of future scenarios, including the so-called representative concentration pathways (RCP) of the International Panel for Climate Change (IPCC) AR5 derived from the global climate model HadGEM2-ES simulations [8]. The pathways are named according to the levels of radiative forcing [49,54]. The IPCC AR5 report is based on the RCP4.5 and RCP8.5 scenarios which consist of changes in greenhouse gas concentrations, aerosols and land use and land cover for the twenty-first century [57], suggesting a radiative forcing (W m −2 ) that would be caused by the increase in the retention of incident solar energy [55]. The RCP4.5 scenario assumes that the radiative forcing will stabilize shortly after 2100 at around 4.5 W m −2 and the RCP8.5 scenario assumes a high radiative forcing of around 8.5 W m −2 by the end of the twenty-first century [26,54]. Zanin et al. [63] observed that, despite the Caatinga and the Cerrado (another Brazilian phytogeographic domain) showing greater resilience against climatic variability and less vulnerability in all scenarios, they may suffer a reduction in their optimal climatic space.
The global temperature projections point to midtwenty-first century increases of up to 2 °C in air temperature (T_air), indicating the possibility of extinction of the phytogeographic domain due to a decrease in the system's water replacement capacity [53]. In the RCP4.5 scenario, the increase in T_air can be up to 2.1 °C and for RCP8.5 it can be higher, up to 4 °C [20]. The increase in T_air associated with climate changes due to global warming could increase evaporation from lakes, weirs, reservoirs and plant transpiration [36]. The increase in the output of water from the vegetation system, due to the increase in T_air, makes the soil drier, causing damage to the plants [37]. Thus, the increase in T_air can lead to a decrease in carbon sequestration by vegetation.
In addition to the changes that the Caatinga can undergo due to climate changes, there are other factors that lead to environmental degradation, such as deforestation. Data from the Instituto Chico Mendes de Conservação da Biodiversidade (ICMBIO) of Brazil show that the total territory of the Caatinga is over 844 thousand km 2 , with 46% of it being illegally deforested and exploited (ICMBIO, 2020).
The present study evaluates T_air trends over four conservation areas in northeastern Brazil (NEB), in the period 1951-2018, and their possible impacts on the species Myracrodruon urundeuva (aroeira), Copernicia prunifera (carnaúba) and Cereus jamacaru DC (mandacarú). These are three plant species present in the Caatinga with high socioeconomic value for the region [7].
The impacts of the projected future climate on the presence of the three species are obtained by the Maxent model [15]. This is an attempt to observe whether changes in the climatic conditions for the phytogeographic domain can lead to the decline of these species, which would lead not only to negative consequences for the vegetative composition of the phytogeographic domain, but also to economic losses for the region.

Study area
The study is carried out on four areas of native Caatinga vegetation in the conservation units located in the North (A1), East (A2), South (A3) and West (A4) of the phytogeographic domain shown in Fig. 1. A1 is the national park of Furna Feia which, with an area of 8518 ha, is located in the municipalities of Mossoró and Baraúna in the state of Rio Grande do Norte. A2 is the Raso da Catarina ecological station, the driest region in the State of Bahia that has very sandy, deep and infertile soils. In this location, the Caatinga is shrubby and very dense [58]. A3 is Serra do Barbado, with 63,652 ha of preserved Caatinga, located in the state of Bahia. A4 is called Dunas e Veredas of the Lower-Middle São Francisco basin, which, also in the state of Bahia, was created in 1997, and extends over 106 ha. Dunas e Veredas is located in the municipalities of Barra, Pilão Arcado and Xique-Xique.

Weather data
Surface air temperature (T_air) data are obtained from the NCEP/NCAR numerical prediction and data assimilation systems [27], and from meteorological stations of the National Institute of Meteorology (INMET) located close to the chosen conservation areas. In the first stage of the analysis of the meteorological data, the correlation coefficients (R 2 ) between the NCEP rainfall data and the historical data of the INMET conventional stations were calculated. The criterion for choosing the closest station to compare with the NCEP rainfall data was based on the

Species description
The aroeira (Myracrodruon urundeuva) is a slow growing tree species, late flowering and gives good quality wood. It is not endemic and occurs in the northeastern, centralwestern and southeastern regions [7] of Brazil. According to the ecological group, it is a species considered late secondary [16]. The Ministry of the Environment classified aroeira as one of the most important species for preservation, indicating its four most important uses: beekeeping, wood, medicinal, producing oil and waxes [38]. It is a species with wide geographic distribution and great ecological plasticity, easily seen in the coastal zone of the country, close to water bodies such as rivers, streams and wetlands, and it can also grow in dunes and dry and rocky terrains [43]. The number of recorded points of presence of the aroeira species for this study is 83. Its population density is around 6 trees per ha.
The carnauba (Copernicia prunifera) is a species of palm tree native to Brazil that usually occurs near rivers and marshy areas, preferring clayey and alluvial soils, being quite resistant to high levels of salinity. It also has a high capacity to adapt to heat, withstanding 3000 h of sunshine per year. According to Duque [13], the age of the palm trees, the type of soil, the climate and the proximity to the sea are factors that influence the production of wax. This native palm tree, admirably adapted to the climatic conditions of the semiarid tropics, has extraordinary resistance, both to the most severe droughts and to floods, two variables determined by the rainfall in the region. According to IBGE [24,22], this is one of the most important species in the Caatinga, both in terms of productivity and in financial terms, offering non-timber forest products. The number of recorded points of presence of the Carnauba species for this study is 64. Its population density is around 18 trees per ha. The mandacaru (Cereus jamacaru P.) is a native species of the Caatinga vegetation, belonging to the Cactaceae family. It grows in stony soils and together with other species of cactus, forms the typical landscape of the semiarid region of the northeast Brazil, found in the states of Piauí, Ceará, Rio Grande do Norte, Paraíba, Pernambuco, Alagoas, Sergipe, Bahia and northern Minas Gerais. Mandacaru plants develop in the drier areas of the northeastern semiarid, in shallow soils and/or on top of rocks, multiplying regularly and covering large areas, showing themselves adaptable to the edaphoclimatic conditions of these regions [19]. The number of recorded points of presence of the mandacaru species for this study is 63. Its population density is around 8 bushes per ha.

Modeling with Maxent
The Maxent model is a tool used for modeling the potential distribution of species, relating their places of occurrence with environmental characteristics [15]. Bioclimatic variable maps, elevation map and species presence records are used as input data. The model gives us the abundance of a given species in terms of percentage number of pixels (or cells) occupied by the species. It also gives us the performance statistic known as area under the curve (AUC).
Spatial data of bioclimatic variables with a resolution of 30 s of a degree of latitude and longitude (~ 1 km), cropped for the Caatinga phytogeographic domain, are obtained from the website http:// world clim. org (v.1.4) for the period from 1970 to 2000 (historical data). The terrain elevation ( Fig. 1) data are obtained from http:// eros. usgs. gov.
Before feeding the Maxent model with bioclimatic variable data, principal component analysis (PCA) is performed for each of the variables in combination with the presence distribution data for a given species. Principal component analysis (PCA) is a method whose basic purpose is to analyze the data used in order to reduce it, eliminate overlaps and choose the most representative forms of data based on linear combinations of the original variables. The goal is to reduce the number of variables. Correlation coefficients between each pair of principal components obtained for the different variables are calculated. (That is, the principal component obtained for a bioclimatic variable, say V1, is correlated with the principal components obtained for V2, V3, etc. Likewise, the principal component of V2 is correlated with those of V3, V4, etc., and so on.) If the correlation between the principal components of any two bioclimatic variables, in combination with the presence distribution, is 0.7 or greater, only one of the variables is considered to affect the species of plants, thus reducing the number of environmental variables that enter the model [60,11]. The understanding involved is that if two variables that affect plant species are strongly correlated, one of the two variables is sufficient to model the potential presence of the species.
Species presence data used as input to the model are in the form of decimal geographic coordinates. The location coordinates of each species are obtained from two online databases, SpeciesLink and Global Biodiversity Information Facility (GBIF). Each of the points obtained was individually visualized, in order to avoid duplicate records and to check if geographic position is wrongly inserted in water bodies or within urban areas. The three species selected for this study are chosen due to their socioeconomic importance for the Caatinga region. The present study is an attempt to project whether changes in the climatic conditions in the phytogeographic domain can lead to the decline of these species, which would lead not only to negative consequences for the vegetative composition of the phytogeographic domain, but also to economic losses for the region.
The validation of the model is based on the analysis of the area under the curve (AUC) which represents the overall accuracy or performance of the model, taking into account the influence of all the variables for which the species is sensitive. Maxent generates relative operating curves (ROC) from contingency tables of 'hits' and 'falsealarms' for different cutoff values of presence potential. The area under the ROC curve is the AUC and is calculated from the model projections and the observed species presence data for the present state. The AUC represents an independent index that assesses the discriminatory power of a model [47] and its calculation provides a unique measure of the model's performance, regardless of the previous choice of cutoff limits for obtaining ROC, as the curve is constructed from multiple cutoff limits (each cutoff limit being responsible for one point on the curve). The AUC value ranges from 0 to 1. An AUC value of 0.50 or less indicates that the model did not perform well, whereas an AUC value of 1.0 indicates total accuracy to discriminate the species [14]. Further details of the model and its performance evaluation are found in Elith et al. From the Maxent output data, maps are prepared in QGIS (Quantum Gis, V. 2.18) that show the potential distribution of species, through climate changes imposed by the scenarios of future projections, the RCPs.

Statistical treatments
The climatological characterization of a region is normally based on 30 years of data. Thus, the NCEP data series from 1951 to 2019 are divided into two parts of 34 years each (January 1951 through December 1984 and January 1985 through December 2019) generating two successive climatologies that could be compared in order to easily visualize whether the locations (conservation units) had T_air changes. Statistical treatment is performed using R software version 3.4.4.
Initially, a graphical and descriptive analysis of the series and the tests by Mann-Kendall (MK) and Pettitt are observed. The Mann-Kendall test is a nonparametric test [31,34], suggested by the Meteorological Organization (WMO) for trend assessment in environmental data time series. This test is suitable for detecting climatic changes in meteorological series. MK considers that, in the hypothesis of stability of a time series, the succession of values occurs independently and the probability distribution must always remain the same, characterizing a simple random series. It is best suited as a test to see whether the time series values tend to increase or decrease over time (monotonic change). The MK test, when compared to other parametric tests, is more robust in terms of deviations from normality and non-stationarity of historical series data, justifying its wide use in trend studies [3]. This test, together with Kendall's Tau correlation statistic (τ), is commonly used to assess trends in environmental, climatological and hydrological data [1].
In the MK test, a trend is considered negative or positive, indicating a decrease or increase in the attributes of the analyzed time series, if the Mann-Kendall score (Kendall's tau) is negative or positive, respectively. Furthermore, the trend indicated by this methodology is considered significant when the p value is lower than α = 0.05. As it is a nonparametric test, it has the advantage of not requiring the data series to be normally distributed [62]. However, the method requires the data to be independent and homogeneous [40]. The homogeneity of the data was tested according to the methodology proposed by Thom [56], verifying that the series is homogeneous. In the entire series, no flaws were found, a fact that did not require any special treatment to the data.
In turn, the nonparametric Pettitt statistical test [45] performs the verification of two entities, Y1, Y2, …, Yt and Yt + 1, Yt + 2, YT, coming from different in the same way. This methodology is able to locate the point where there was a sudden change (point of change) in the mean of the time series, in addition to providing information on the homogeneity of the data in the analyzed series.

Analysis of rainfall and air temperature (T_air)
The correlation coefficients between the NCEP data and the Apodi station (closest to A1) was R 2 = 0.8. Between NCEP and Paulo Afonso station (closest to A2) it was R 2 = 0.7, between NCEP and Itauçu (closest to A3) it was R 2 = 0.7, and finally for A4, Barra station was used, located in the state from Bahia, with R 2 = 0.65.
The rainy season lasts about four months for A1 and A2 conservation areas from January to April and about five months for A3 and A4 from November to March. The dry period lasts about seven months, from June to November for A1 and A2 and from April to October for A3 and A4.  Fig. 2. The place with the highest mean T_air is A1 despite having the highest amount of rainfall recorded over 68 years. This characteristic is due to its location, which is in an area of Northern Sertaneja Depression, with an average altitude of about 100 m. It is known that areas with lower altitudes tend to be warmer and with smaller thermal amplitudes. A2 has an average altitude of 276 m, and an average T_air of 23.8 °C. A3 is mountainous and has altitudes ranging between 480 and 2021 m, with the lowest accumulated annual rainfall, perhaps in a rain shadow region. The lowest mean value of T_ air (22.2 °C) of the NEB region occurs at A3. A4 has an average altitude of around 870 m but has warmer temperature and higher rainfall than A3. Figure 2 shows systematic and generalized increase in T_air in all the four locations.
The MK and Pettitt test results for the mean monthly T_air are given in Table 1. The MK test clearly indicates that air temperatures have increased in all areas. The average annual values showed positive trends (τ = + 0.025, + 0.171, + 0.181 and + 0.267) and all these trends are statistically significant (p value < α = 0.05). The Pettitt test showed positive T_air changes for A1 and A2 since 1996. For A3, the highest positive trend occurred from 1993 and for A4 from 1997.   Fig. 2. Table 2 shows the decadal mean T_air values, for the rainy and dry seasons separately, during the period studied. First, we note that the dry season temperatures are, in general, lower than the wet season temperatures in A2 and A3, while the wet season temperatures are lower in A1. In A1, the daytime convective cloudiness in the wet season masks solar radiation. In A2 and A3, the late evening rains and cloudiness do not affect the daytime temperature. A4, in general, is also warmer in the dry season.
In A1, the difference in the mean T_air value between the last three decades and the first four decades was 1.11 °C for the rainy season and 1.41 °C for the dry season. In A2, the rainy season showed higher increase in temperature, 1.41 °C, compared to the dry season warming of 0.95 °C. In A3, the increases were almost uniform in both the seasons, 1.31 °C for the rainy season and 1.22 °C for the dry season. In A4, the increase was 1.03 °C for the rainy season and 1.20 °C for the dry season. Table 2, thus, clarifies the seasonal differences in the warming trends and indicates more intense warming near the coast. Figure 3 shows maps of the differences between the two climatologies of T_air, one for the period 1951 to 1984 and the other for 1985 to 2019. Each panel represents a month and it is possible to see the predominance of positive values in the differences over the northeast Brazil (NEB) region, indicating increases from the first 34-year mean to the last 34-year monthly mean T_air. The annual cycle of spatial distribution of T_air differences over NEB shows migration and strengthening of higher positive values from the southern portion to the center and to the northern portion, during the first half of the year. In August and September months, the differences weaken and then strengthen in November and December. There is again a weakening from December to January near the north coast which is perhaps due to the dampening of T_air in the beginning of the rainy season in January. In austral autumn, the southward migration of the ITCZ heralds the beginning of the rainy season on the northern parts of NEB [48], which may explain why the increase in T_air was concentrated in the central part of the NEB during many months. The study of Santos et al. [50] Figures 4, 5 and 6 show the potential intensity of presence distribution simulated by the Maxent model for the species aroeira, Carnauba and mandacaru, respectively. Potential intensity of presence can be interpreted as the probability of occurrence of a given species in a given pixel or as the percentage of the number pixels with the presence of the species in a given area. Each figure has three panels, one for the current state, second for the future state in  Table 3) are shown in red. The gradation of potential intensities is defined as follows: pixels with values between 0 and 0.01 have null potential, between 0.01 and 0.25 have low potential, 0.25 and 0.50 have moderate potential, 0.50 and 0.75 have good potential and 0.75 to 1.0 have high potential.

Maxent model projections
Examining the areas with high potential of occurrence (> 0.75) for aroeira (Fig. 4), the spatial distributions in the future scenarios RCP4.5 and RCP8.5 in relation to the current state show some notable differences in the biogeographic arrangements. In the present state, high potential for the presence of the species is concentrated in Bahia state (left panel). In the future scenarios, the concentration in central Bahia is reduced and is increased in both the future scenarios (right column panels) in the state of Rio Grande do Norte (RN). That is, a shift in the presence potential of aroeira to north is projected. Table 3 shows the expected changes in different presence potential gradations. Combining the null and low potential areas, the area is projected to increase from 31.7% in the present state to 47.0% in the future states under the RCP4.5 and RCP8.5 scenarios. The table also shows a reduction of moderate to high potential areas to about 53.0% in 2050 from 68.0% in the present state (current scenario). This shows the sensitivity of aroeira species to temperature changes. The environmental variables that most contribute to the potential occurrence of aroeira, obtained from PCA, are the mean annual temperature (Bio1), isothermality (Bio3) and precipitation in the coldest quarter (Bio19). Surprisingly, the difference between the two RCPs is small, although the increase in air temperature in the Caatinga phytogeographic domain in the RCP8.5 scenario is 3 °C compared to 2 °C for the RCP 4.5 scenario. The areas with the greatest potential for the occurrence of Carnauba in the present state are situated mostly in the state of Ceará (CE), Paraíba (PB) and Piaui and are spotty (Fig. 5). They are places with good drainage. In the future scenarios (panels in the right column, Fig. 5), the high potential areas shrink or disappear in Piaui and Paraiba and appear in the coastal places of Ceará and in the extreme western portion of Pernambuco state. Table 4 shows the expected changes in the presence potential of the species in future scenarios compared to the present state in different gradations of potential. The combined good and high potential areas decrease slightly from 25.3% in the present state to 24.0% in the RCP4.5 scenario and to19.5% in the RCP8.5 scenario. A shift in the habitat suitable for the occurrence of Carnauba from the interior of the region to the coast of Ceará is notable.
The results of the PCA showed that the variables that most contribute to the presence of Carnauba are: minimum temperature in the coldest month (°C) (Bio6), precipitation in the driest quarter (mm) (Bio17) and precipitation in the warmest quarter (mm) (Bio18). This information suggests that the species adapts well to places with temporal water stress. Field data indicate that the carnauba species adapts well to intermittent wetlands, withstanding periods of flood and drought [12].
Mandacarú habitat with good and high potential is present in PB, RN and the eastern parts of BA in the present state (Fig. 6). Cavalcante et al. [5] mentioned that it is a species of low demand of water that can occur over a large part of the phytogeographic domain. However, in the present study, the high potential for its presence is seen more concentrated over the eastern belt of the Caatinga phytogeographic domain. The projections for future scenarios indicate that the areas with good potential will reduce in RN and will increase in western parts of the states of Alagoas and Sergipe. In the current scenario, there is 23.5% of the area with a good or high potential for the presence of the species. In the RCP4.5 scenario, this value decreases to 17.5%, and in RCP8.5, it jumps to 22.9% (Table 5). In the extreme RCP8.5 scenario, the rainfall changes may be favoring a rise in the distribution potential of mandacarú compared to RCP4.5. Areas with null and low potential for mandacarú species for the current scenario are 39.7%, increasing to 62.5% in the RCP4.5 scenario and decreasing to 44.0% in the RCP8.5 scenario (Table 5).  The variables that most contribute to the species distribution is precipitation in the driest month (mm) (Bio14), followed by precipitation in the wettest month (mm) (Bio13) and precipitation in the warmest quarter (mm) (Bio18). That is, it does not depend on temperature changes in the future scenarios as much as it depends on the precipitation. Mandacaru blooms in the rainy season [63] and presents seasonal reproductive events linked to the rainy season. The fruiting strategy during the beginning of the rainy season allows for greater efficiency in the germination process with greater use of water in the soil [32], 42].
Lastly, the average AUC values found in the present study for the three species were 0.84 for aroeira, 0.87 for carnaúba and 0.83 for mandacaru. AUC value over 0.8 indicates good accuracy of the model, and therefore, the potentials of presence presented above can serve as a guide for future estimates of evolution of the three species.

Discussion
We observe that A1 has higher T_air values than the other three areas, despite being the place with the highest rainfall among the four conservation areas studied due to its location in a low lying area with an average altitude of about 100 m. The lowest monthly mean values of T_air observed in A3, although it receives only 753 mm year −1 of rain, are due to its location in a well-drained region, with several bodies of water around the conservation area and with high levels of Improved Vegetation Index [59], in press), which may contribute to cooling the air temperature.
In general, the years with the lowest mean air temperature record for all locations were concentrated in the first climatology  and the highest mean air temperature values were recorded in the second climatology (1984-2019), which indicates warming. In the respective years in which there was maximum and minimum T_air for the analyzed regions, El Niño (possible cause of higher temperatures than normal) and La Niña (possible cause of lower temperatures) events also contributed. The A2 area, however, registered higher T_air in the La Nina year 2008 and lower T_air in El Niño year 1963, the reasons for which need investigation. Studies focusing on the El Niño and La Niña episodes have shown that these events can influence the rainfall pattern over NEB and, consequently, interfere with the air temperature [28,29]. El Niño aggravates the lack of rain in the NEB region, influencing the increase in air temperature. The changes observed in T_air were also reported in other works, such as Marengo [33] who presented projections for the twenty-first century pointing to warming throughout Brazil, reduction in precipitation in NEB, possibly as a cause of global warming and an increase in T_air at the end of the twenty-first century [35].
The increase in T_air, when associated with vegetation suppression, leads to aridification. Even in the so-called conservation units with a large extension of Caatinga vegetation, there was an increase in T_air. This can increase the risk of hydrological deficiency in the region by favoring an increase in the rate of evaporation. In addition, a greater number of hours of exposure to sun (approximately 3.000 h of sun per year) can further aggravate the hydrological deficit in this region for most of the year [2].
Salazar et al. [49] estimated a aridification situation in the semiarid northeast in the period 2090-2099, for an air temperature increase scenario of 2-6 °C in the A2 region, and 1-4 °C in the A1 region, indicating a direct influence on the characteristics and distribution of vegetation. Thus, changes in T_air values will have a direct influence on vegetation characteristics and distribution.
Kumar et al. [30] explained that part of the upward trend in T_air is not only due to global warming but also due to the higher frequency of El Niño events. It is important to emphasize that factors such as the increase in urban areas, reduction of vegetation cover, change in the atmospheric composition due to the increase in gases such as carbon dioxide (CO 2 ), all contribute to the change in T_air [50]. In the state of Paraíba, for example, Dantas et al. [10] showed an increase in air temperature, and they attributed this to the density of buildings, change in land use and urban growth.
The Maxent model utilized to project future situation of plant species directly depends on the existing species occurrence data, and thus, the output is sensitive to the quality of observed data. Errors in the latitude and longitude position in the occurrence data may cause the model make imprecise projections. For example, a species known to be endemic to the Amazon Basin may make its presence in the geographic area of Caatinga. Sometimes, in the place of the coordinates of the location of a species the coordinates of the nearby city are coded. It is also found that the resolutions of the numerical values at different places are different and these data are thoroughly scrutinized before feeding to the model. Precipitation, photoperiod and air temperature are related to flowering, fruiting, fall and sprouting of the mastic tree [44]. While precipitation positively affects the production of new leaves and budding, it is in the reduction of rainfall that flowering, fruiting and leaf fall occur, and it is only with increase in air temperature that the fruit ripens, showing a direct influence of rainfall and air temperature in the reproductive and vegetative phenophase [41].
In the annual average temperature map (Bio1, not shown) for the current climate, the occurrence of aroeira is concentrated in areas with air temperature between 20 and 30 °C, which is the same air temperature range considered optimal for germination of its seeds according to Medeiros et al. [34]. In the scenario RCP4.5, it is observed that there is an expansion of areas to the north of the phytogeographic domain and coastal areas with increased air temperature. The details of the relationship between this species and environmental variables are not very clear. It can be said, however, that this species needs a dry period, with an increase in T_air, for the maturation of its fruits. T_air is the only factor that affects the occurrence of ripe fruits [41].
A survey carried out to determine the relationship between climatic factors and aroeira seed germination and the trend of this germination for future climate scenarios showed that a reduction of the number of weeks with minimum rainfall directly affects seed germination and seedling recruitment, and that the quality of the seeds produced by aroeira is directly related to the maximum, average and minimum temperatures, average and minimum humidity, indicating that fluctuations in climate variables will lead to a reduction in the physiological quality of the seeds [43].
A study carried out by Francisco et al. on the detection of changes in the Caatinga vegetation showed that places with greater potential for recovery have alluvial soils with greater moisture available during the year and Litholic Neossolos in areas of steeper slope. In the state of Paraíba, for example, the central region showed low potential for the presence of Carnauba in the current scenario, worsening in scenarios RCP4.5 and RCP8.5. The central region, where the Taperoá river basin is located, is currently in the process of aridification as per Francisco et al. The potential for occurrence is lower in all places which are already susceptible to aridification.
Mandacaru is one of the main species, of the Cactaceae family, found in Brazil [6], it is a highly resistant species that adapts to places with little water resource. Despite its adaptation to the semiarid climate, mandacaru may suffer the effect of global climate changes that may favor the advance of aridification, and influence the distribution of this species [26]. In general, there is a spatial variation in the distribution of mandacaru in the two RCP scenarios, with a decline in the species abundance compared to the current scenario. The increased distribution of the species in the RCP8.5 scenario compared to RCP4.5 can be explained by its resistance to the dry period.
The increase in air temperature and the reduction of rainfall in the RCP scenarios can contribute to the aridification process, with negative consequences especially in the RCP8.5 scenario where there is spatial expansion of the mandacaru over a large part of the Caatinga with AUC values considered good, featuring an arid landscape. Williams et al. [61] explain the adaptation of cacti as a relationship between the physical structure of the species (morphology) and the climate conditions where the species will be established.
In an attempt to contain the advance of the degradation of native vegetation and aridification, the conservation areas created by the Federal Government in the Caatinga are of great importance. With the creation of such units, the protected area for conservation increased in the phytogeographic domain by about 7.5% of the total area of Caatinga. Even so, this will continue as one of the least protected phytogeographic domain in the country, as just over 1% of these units are fully protected [39].

Conclusion
Taking into account current and future climate conditions and biotic data for species distribution models, we observe that there is a possibility of contraction of area with high distribution potential for the three target species and expansion of areas with low potential by the end of this century.
The time series of T_air for the four areas of the Caatinga phytogeographic domain in northeast Brazil (NEB) studied here showed positive trends, indicating that important changes are present in the data. For all the four areas, the periods of increase in T_air were accompanied by a reduction in rainfall, which could lead to an increase in the aridification process. The records of El Nino and La Nina, in the same years as warmer and less warm years, may indicate the influence of the phenomenon on local climate conditions. The conservation areas studied here represent less than 12% of the region, and a large part of it has already been deforested, aggravating the increase in T_air. The deforestation, in addition to the global warming, can worsen the temperature rise.
The Maxent model results for the three important vegetation species in the phytogeographic domain indicate future reductions in potential areas for their presence. In the most catastrophic scenario, RCP8.5, the reduction in species presence is more accentuated for Carnauba, followed by aroeira and mandacaru. In the conservation | https://doi.org/10.1007/s42452-021-04886-w areas A1 and A4 (Fig. 1), high distribution potentials are observed for the aroeira and carnaúba species in the current scenario and in RCP4.5. These areas are the ones that register the highest values of annual precipitation.
In A3, despite being an area with little more than 700 mm year −1 of rain, there is a high potential for the distribution of aroeira and mandacaru, possibly as a consequence of the local hydrology. This may be the explanation for the high distribution potential that is maintained in this conservation unit in the future scenarios also.
A2 is an area with rarefied soils and dominated by shrub and herbaceous vegetation. In addition to these conditions found in the field, it is an area with very little annual rainfall. This may be the reason for the species occurrence potential to fluctuate between low and moderate for all the three species.
The possibility of the absence of these species in near future (approximately 30 years) may have consequences for the maintenance of animals and affect the local economy, since these species are highly productive and present high economic value. Its absence can pose a serious social and environmental problem. Thus, further investigations to monitor the increase in T_air and the consequences of this change for the vegetation of the Caatinga phytogeographic domain are needed.