An Announced Extinction: The Impacts of Mining on the Persistence of Arthrocereus glaziovii , a Microendemic Species of Campos Rupestres

: The mountaintops of eastern Brazil harbor the highest rates of plant endemism in South America. However, local biodiversity faces constant threats due to habitat loss and mining activities. About 89 rare and endangered species are exclusive to this region, including the threatened species Arthrocereus glaziovii . This study aims to evaluate the potential distribution of A. glaziovii based on abiotic variables and soil elements and to characterize the distribution of mineral titles that may restrict the species’ occurrence areas. We used the Bioclim, Domain, MaxEnt, GLM, and Random Forest algorithms to model this ecological niche under future climatic scenarios, in addition to modeling the layers of mineral titles corresponding to areas already mined and those slated for future mining projects. Our predictions indicate an expansion in the future distribution of A. glaziovii . Nevertheless, the future predicted occurrence areas of the species are already compromised due to mining. According to our findings, we emphasize the looming threat of the predicted extinction of this species. Therefore, implementing conservation strategies to ensure the survival of A. glaziovii is imperative.


Introduction
Brazil is the second-largest mineral producer in the world [1].Thus, mineral extraction is one of Brazil's most important economic activities, mainly in the eastern Brazil mountaintops, which encompasses the threatened rupestrian grassland (hereafter, campo rupestre) [2].The campo rupestre is characterized by a rich mosaic of herbaceous vegetation and sparse shrubs [3].The region's soil is compact, deficient in water and nutrients, and high in heavy metals [4,5].
The rocky outcrops in the campo rupestre, particularly the iron ones (hereafter, canga), are rich in endemic plant species and substantially contribute to the biodiversity of the domains in which they are inserted, namely the Cerrado and Atlantic Forest [6,7].These plants evolved on extremely nutrient-poor soils [8] and are adapted to distinctly seasonal climates [9].Plants from these locations are very susceptible to soil removal [10].They have a limited dispersal distance and slow growth [11], hindering the natural regeneration of altered campos rupestres [12].It has recently been observed that the combination of this edaphic endemism and limited seed dispersal threatens the persistence of these mountaintop species [13].
Vegetation in the campo rupestre is highly threatened due to its restricted area [14] and the climate changes forecasted for the future [15], making their conservation extremely important, urgent, and necessary.The climate change scenarios predicted suggest catastrophic consequences for this ecosystem [15], and models project a reduction of up to 95% in the campo rupestre's current suitable area [16].
In addition to the threats mentioned above, there is a high mining pressure in this ecosystem [11] since we find in the campo rupestre one of Brazil's largest iron ore extraction areas, the Iron Quadrangle (local name Quadrilátero Ferrífero (Q.F.)), placing mining as the main factor of destruction in this region [17].From 1960 to 2012, the Q.F.lost about 40% of its total area of cangas (about 100 ferruginous outcrops) to iron ore extraction [7,17].Today, there are nearly 225 remaining ferruginous outcrops, and 83% of these show signs of loss or degradation by mining [7].
The Q.F. is the central mineral-producing province in the world in terms of iron ore [18,19], and the loss of species from mining activities is marked by a high, irreversible environmental impact [20].The intense mining activities practiced in the area since the early 17th century on a large scale [21] continue to impact the flora of the campo rupestre, resulting in areas with exposed subsoil or gravel, which, even after decades, no natural recovery has been recorded [22].In addition, about 89 rare and threatened species in the Q.F. are plants with high conservation value, occurring exclusively in cangas in the region [17].One of these species on the mountaintops of rocky fields impacted by mineral exploration is the microendemic and endangered cactus Arthrocereus glaziovii, which has suffered from a reduced area of occurrence [14,23].
In this study, we predicted the potential and future geographical distribution of A. glaziovii, an exclusive, threatened, and rare species of the campo rupestre.Specifically, we (i) study the potential distribution of A. glaziovii, based on abiotic variables and soil elements, to understand the impacts of these climatic variables in the future distribution of the species and its ecosystem and (ii) characterize the distribution of mineral titles that may restrict the areas of occurrence of A. glaziovii in order to understand the impact of mining on the potential distribution of the species.

Study Species
The Cactaceae family is neotropical, distributed almost exclusively on the American continent, and its main center of diversity is East Brazil [24,25].Cacti have been the subject of diverse studies since they are used as food, fodder, medicinal, and ornamental plants.Thus, the Cactaceae family has always sparked interest because of their morphological peculiarities and chemical properties [25,26].The highest biodiversity is observed in arid and semi-arid regions, where cacti serve as crucial resource reservoirs for numerous species of vertebrates and invertebrates, particularly during droughts characterized by resource scarcity [26].
The genus Arthrocereus stands out among the Cactaceae family, comprising six species, all endemic to Brazil.They occur in shrub and sub-shrub forms and have particular habitats.Consequently, they are restricted to their localities [27].Arthrocereus glaziovii (K.Schum.)N. P. Taylor & Zappi is among the three taxa of this genus present in Southeastern Brazil and stands out as a rare columnar-looking cactus.The species is considered microendemic, as it has few populations restricted to the ferruginous outcrops of the campos rupestres [28].This cactus is considered threatened with extinction, according to the International Union for the Conservation of Nature (IUCN) [29], due to anthropogenic activities which cause the suppression of its habitat by removing soil for mining [30].The species is considered a priority in the recovery of the ferruginous rocky fields of the Q.F., as it is unique to the region [31].

Study Area
Considering that the microendemic distribution of the species manly occurs in places rich in iron ore, the study region is concentrated in the Espinhaço Range, a mountain range characterized by the dominance of diverse and abundant rupestrian grasslands, where the Q.F. is located.The Espinhaço Range is situated within the Brazilian Central Plateau and stands as one of the most significant mountain ranges in the country, boasting one of the highest levels of floristic diversity in South America [4,9,32,33].
In the large iron ore deposits of the mountain ranges that delimit the Q.F., this is where the ferruginous outcrops, known as cangas, are located [14].The cangas are composed of interspersed quartzitic, granitic, and hematitic outcrops, forming different microhabitats where A. glaziovii occurs [30] (Figure 1).With the region's climate being considered as high-altitude subtropical (dry winters and wet summers), according to the Köppen climate classification [34], the campo rupestre presents daily high thermal amplitudes (minimum of 14.1 • C and maximum of 23.7 • C) and low relative humidity [30,35].Thus, soil, relief, and climate conditions are the main factors for establishing its unique and dominant vegetation cover [36].
characterized by the dominance of diverse and abundant rupestrian grasslands, where the Q.F. is located.The Espinhaço Range is situated within the Brazilian Central Plateau and stands as one of the most significant mountain ranges in the country, boasting one of the highest levels of floristic diversity in South America [4,9,32,33].
In the large iron ore deposits of the mountain ranges that delimit the Q.F., this is where the ferruginous outcrops, known as cangas, are located [14].The cangas are composed of interspersed quartzitic, granitic, and hematitic outcrops, forming different microhabitats where A. glaziovii occurs [30] (Figure 1).With the region's climate being considered as high-altitude subtropical (dry winters and wet summers), according to the Köppen climate classification [34], the campo rupestre presents daily high thermal amplitudes (minimum of 14.1 °C and maximum of 23.7 °C) and low relative humidity [30,35].Thus, soil, relief, and climate conditions are the main factors for establishing its unique and dominant vegetation cover [36].This region supplies the demand for iron ore in Brazil and abroad.The annual Brazilian iron ore production is approximately 400 million tons, and more than half of this production is concentrated in the Q.F.[37].In 2019, the iron ore production of the largest Brazilian mining company, and third in the world, located in Minas Gerais state (M.G.), was 301.9 million tons [38].

Species Distribution Modeling
To determine the ecological niche of A. glaziovii, we collected georeferenced specimen records from the Global Biodiversity Information Facility (GBIF) [39], from speciesLink databases [40], from the herbarium from the Federal University of Minas Gerais (UFMG), and from literature data [41].Records containing wrong geographic coordinates, no coordinates, or incomplete data were excluded from our analysis.
We used ecological niche modeling (ENM) to model the potential distribution of A. glaziovii in the Minas Gerais state region, which comprises a large part of the Espinhaço Range, including the Q.F.For the modeling, we use five types of ENM algorithms.It is This region supplies the demand for iron ore in Brazil and abroad.The annual Brazilian iron ore production is approximately 400 million tons, and more than half of this production is concentrated in the Q.F.[37].In 2019, the iron ore production of the largest Brazilian mining company, and third in the world, located in Minas Gerais state (M.G.), was 301.9 million tons [38].

Species Distribution Modeling
To determine the ecological niche of A. glaziovii, we collected georeferenced specimen records from the Global Biodiversity Information Facility (GBIF) [39], from speciesLink databases [40], from the herbarium from the Federal University of Minas Gerais (UFMG), and from literature data [41].Records containing wrong geographic coordinates, no coordinates, or incomplete data were excluded from our analysis.
We used ecological niche modeling (ENM) to model the potential distribution of A. glaziovii in the Minas Gerais state region, which comprises a large part of the Espinhaço Range, including the Q.F.For the modeling, we use five types of ENM algorithms.It is essential to use several algorithms to ensure the prediction is as accurate as possible.Each method has a different predictive ability, degrees of complexity, and stability, gener-ating different responses for each variable.The algorithms were applied to all compiled georeferenced record points (Figure 2; Table S1).
Conservation 2024, 4, FOR PEER REVIEW essential to use several algorithms to ensure the prediction is as accurate as possible.E method has a different predictive ability, degrees of complexity, and stability, generat different responses for each variable.The algorithms were applied to all comp georeferenced record points (Figure 2; Table S1).
To avoid redundancy in the variables, we performed a pair-by-pair Pear correlation for each dataset (climate and edaphic), available in the Raster package [45 R software 4.3.2[46].Only one variable was selected with 70% or more correlation to av
To avoid redundancy in the variables, we performed a pair-by-pair Pearson correlation for each dataset (climate and edaphic), available in the Raster package [45] of R software 4.3.2[46].Only one variable was selected with 70% or more correlation to avoid an excess of autocorrelated variables, following [47].These chosen variables were submitted to a Random Forest selection, eliminating variables irrelevant to the data and ordering the variables in order of importance.This selection was made using the VSURF package [48].

Modeling and Evaluation Procedures
To carry out the modeling, we chose five algorithms: Bioclim, Domain Generating Algorithm (DGA), General Linear Model (GLM), Maximum Entropy (MaxEnt), and Random Forest (R.F.), commonly used in ENM analysis.
The Bioclim algorithm [49] is a climate envelope model that calculates similarity in amplitude based on presence-absence.This algorithm compares the values of variables in any location with the values that appear in the grid of species occurrence points to identify locations with better environmental suitability [49].The Domain algorithm [50] calculates the Gower distance to find the minimal distance between environmental variables at a location and the known locations of occurrence.The General Linear Model (GLM) is a logistic regression model that analyzes cause and consequence, indicating which variables are responsible for the presence or absence of the species.The Maximum Entropy (MaxEnt) model [51] is a machine learning model based on maximum entropy.This method creates random absence points in the background, which, together with the presence points, will determine the possible area of occurrence and restrict the distribution according to the variables [52].The Random Forest (RF) algorithm [53] is a regression and machine learning model that creates "trees".It groups them to obtain a prediction with greater accuracy and stability.The advantage is that it adds extra randomness to the model without overfitting the data.
As we do not have actual absence points, we created one thousand pseudo-absence points for the studied species to be used in the GLM, MaxEnt, and R.F. models.After analyzing each current model generated by the algorithms, we created a mean between all models in each set (climate and edaphic).We used 80% of the points to construct the models and the remaining 20% of the points for model validation.We used the Dismo package [54] for the Bioclim, Domain, and MaxEnt models and the randomForest package [55] for the RF model.All analyses were performed using R version 3.5.3software [46].

Mining Titles
The layers of the mineral titles were obtained in vector format from Instituto Pristino's Geoenvironmental Digital Atlas database [56].We use 52 layers corresponding to areas that are already mined and those that will be mined in the future.

Selected Climatic and Edaphic Variables
After performing Pearson's correlation, 6 of the 20 initial climatic variables were selected (Tables S4 and S5) and submitted to VSURF.Among these six input climate variables, seasonal precipitation (BIO15), mean annual temperature (BIO1), annual precipitation (BIO12), temperature seasonality (BIO4), and annual temperature variation (BIO7) were the most influential according to the VSURF importance values (Table 1), and solar radiation was not significant (SRAD).Ten of the eighteen initial soil variables were selected by Pearson's correlation (Table S4) and submitted to VSURF.Among these ten selected edaphic variables, silt, potassium (K), nitrogen (N), and hydrogen (H) were the four most influential, and only aluminum (Al) did not influence the model.After excluding the irrelevant variables, the other variables were classified in order of influence (Table 1).

Current Potential Environmental Suitability for Arthrocereus glaziovii
All models had an area under the curve (AUC) > 0.98, showing their effectiveness.When considering each of the models individually for both climatic and edaphic variables, together, the most restrictive algorithms for the area were Bioclim and R.F., which identified a small space with environmental suitability located in the Q.F.(Figures 3A,D and 4A,D).Domain and MaxEnt were the algorithms that presented a moderately wide distribution, identifying restricted distributions in the center, south, and southeast of Minas Gerais state and in the mountainous areas of the studied site (Figures 3B,E and 4B,E).The algorithm with the broadest distribution was the GLM, which identified almost every area analyzed with a probability of environmental suitability for both variables (Figures 3C and 4C).According to the Bioclim climate model, the average annual temperature allowing for the occurrence of A. glaziovii is very restricted (between 16 and 20 • C), with annual precipitation ranging between 1500 and 1600 mm, seasonal precipitation ranging between 84 and 86 mm, the temperature in the hottest month ranging between 18.2 and 19.1 • C, and a variation in annual temperature between 1.7 and 1.8 • C. According to the Domain climate model, the restriction of species occurrence was similar, with a yearly average temperature between 16 and 21 • C, a variation annual precipitation between 1450 and 1650 mm, seasonal precipitation between 82 and 89 mm, the temperature in the hottest month being between 17.5 and 20.1 • C, and a variation in annual temperature between 1.6 and 1.8 • C. For MaxEnt, the most critical climate variable for the niche was the yearly mean temperature, followed by seasonal precipitation, annual precipitation, and annual temperature variation.For this model, the temperature of the warmest month was not significant.As for the R.F., the most essential variable for the niche was annual precipitation, followed by yearly temperature variation, seasonal precipitation, average annual temperature, and temperature of the hottest month.According to the Bioclim climate model, the average annual temperature for the occurrence of A. glaziovii is very restricted (between 16 and 20˚C), annual precipitation ranging between 1500 and 1600 mm, seasonal precipitation between 84 and 86 mm, the temperature in the hottest month between 18.2 and 19.1 ˚C and the annual temperature variation between 1.7 and 1.8˚C.According to the Domain climate model, the restriction of species occurrence was similar, with the average annual temperature between 16 and 21˚C, the annual precipitation with a variation between 1450 and 1650 mm, the seasonal precipitation between 82 and 89 mm, the temperature in the hottest month between 17.5 and 20.1 ˚C and the annual temperature variation between 1.6 and 1.8˚C.For MaxEnt, the most critical climate variable for the niche was the mean annual temperature, followed by seasonal precipitation, annual precipitation, and annual temperature variation.For this model, the temperature of the warmest month was not significant.As for the R.F., the Using Bioclim's edaphic model, sand and clay showed high levels, from 41 to 53 and 25 to 27, respectively, followed by silt with 15 to 20, pH at 5, H from 3 to 4, and Mg 0 to 1. Using the Domain model, the soil presented similar or close values, sand presented values between 30 and 60, clay between 15 and 38, silt at 10 to 15, pH at 5, H at 1 to 6, and Mg again at 0 a 1. Nitrogen and K had values equal to 0. For the MaxEnt model, the most critical edaphic variable for the niche was Mg (50%), followed by K, N, clay, and silt.As for R.F., the most essential variable for the niche was the concentration of N, followed by silt, clay, K, and Mg.Sand, H, and water pH were not important in any model.
The average of all models showed a wide potential distribution of the species (Figure 4F).According to the average cut-off of all environmental niche models, the potential distribution of A. glaziovii is currently located in the central-southeast and south of M.G., with a particular concentration in the Espinhaço Range (and the Q.F.) and the northern regions, as well as south (Serra da Mantiqueira) and north of the study area (Figure 3G).

Mining Titles
Overlapping all layers of mineral titles in the Q.F., we observe few remaining miningfree areas in this region (Figure 5).Iron ore extraction is the most frequent mining activity, followed by gold ore extraction.In fact, the Q.F. has numerous deposits of iron and gold ore, which makes the state of Minas Gerais one of the largest producers of these minerals in the world [18].So, it becomes clear that the distribution of A. glaziovii, both actual and potential, is over an essential area for mining extraction.The distribution of future mines largely overlaps current occurrence areas and suitable regions for potential distribution of the species.
Using the Domain model, the soil presented similar or close values, sand presented va between 30 and 60, clay between 15 and 38, silt at 10 to 15, pH at 5, H at 1 to 6, and again at 0 a 1. Nitrogen and K had values equal to 0. For the MaxEnt model, the m critical edaphic variable for the niche was Mg (50%), followed by K, N, clay, and sil for R.F., the most essential variable for the niche was the concentration of N, followe silt, clay, K, and Mg.Sand, H, and water pH were not important in any model.
The average of all models showed a wide potential distribution of the species (Fi 4F).According to the average cut-off of all environmental niche models, the pote distribution of A. glaziovii is currently located in the central-southeast and south of M with a particular concentration in the Espinhaço Range (and the Q.F.) and the nort regions, as well as south (Serra da Mantiqueira) and north of the study area (Figure 3

Mining Titles
Overlapping all layers of mineral titles in the Q.F., we observe few remaining min free areas in this region (Figure 5).Iron ore extraction is the most frequent mining acti followed by gold ore extraction.In fact, the Q.F. has numerous deposits of iron and ore, which makes the state of Minas Gerais one of the largest producers of these min in the world [18].So, it becomes clear that the distribution of A. glaziovii, both actual potential, is over an essential area for mining extraction.The distribution of future m largely overlaps current occurrence areas and suitable regions for potential distribu of the species.

Discussion
Our results demonstrate that the potential occurrence area of A. glaziovii, conside both climatic and edaphic variables, is larger and covers areas outside the Q.F.How currently licensed and under-exploration mining titles and areas licensed for fu mining activities threaten the entire identified area, jeopardizing the species' surviva

Discussion
Our results demonstrate that the potential occurrence area of A. glaziovii, considering both climatic and edaphic variables, is larger and covers areas outside the Q.F.However, currently licensed and under-exploration mining titles and areas licensed for future mining activities threaten the entire identified area, jeopardizing the species' survival.
Seasonal and annual rainfall, temperature in the warmest month, and annual temperature were the abiotic factors determining the geographic distribution of A. glaziovii.Temperature and rainfall are important, as they influence the length and intensity of the cactus's reproductive period [57].Precipitation is the primary factor influencing seedling emergence and the survival of cacti.Studies have indicated that spatial and temporal fluctuations in precipitation significantly shape the establishment patterns observed in this family [25].A. glaziovii has an annual flowering period lasting from 1 to 5 months, starting between the end of the dry period and the beginning of the rainy period [58,59].The temperature positively influences the production of flowers and fruits, but rainfall does not influence this factor [60].The fruiting of the species takes place in the middle of the rainy season, and both temperature and rainfall positively influence this phenophase [60].
Increases and decreases in temperature influence the germination of the species, with an optimal temperature between 25 and 30 • C, even in the rainy season, needed to establish seedlings [3,61].
The species A. glaziovii currently occurs in the Q.F.region in the center-southeast of the Minas Gerais state, and the Domain and MaxEnt models indicated a potentially larger distribution area for this species in this same region.The Bioclim and R.F. models showed a minor occurrence compared to the current model.The GLM indicated a wide distribution throughout the Minas Gerais state.The cut of the averages of the models stated an occurrence area located to the south (climate model) and north of the state of Minas Gerais.The potential distribution models of A. glaziovii showed, therefore, that this species has a potential geographic distribution restricted to the central-southeast region of the Minas Gerais state, with a concentration along the Espinhaço Range, to the south and north of this state.
The south area of Minas Gerais state in the climate models is called Serra da Mantiqueira, and its indication in the model is probably due to it being a region with high elevations.According to Koppen's classification, the climate in this region is subtropical in altitude, equal to that of the species' current occurrence area [62].In Serra da Mantiqueira, the vegetation found is called campo de altitude.It grows under granitic outcrops of igneous or metamorphic rocks, which differ from the physiognomies occurring on the cangas of the Q.F. and the quartzite fields of the Espinhaço Range [62].Although both occur on rocks and have similar climates, they are of different lithotypes, explaining why the edaphic model did not point to this region [62].
Only the MaxEnt model of climate variables indicated a small area in the north of Minas Gerais state as a possible area of environmental suitability for the occurrence of A. glaziovii.This region is a transition area of Cerrado and Caatinga domains with hematite in the soil and a humid tropical savannah climate.The driest season coincides with winter, and rainfall is less than 60 mm [63].The soil of this indicated region also presents high concentrations of Fe, as seen in the areas of occurrence of the species in the Q.F.However, this element was not analyzed in the models, as no studies prove its direct influence on the plant, indicating that the species has some mechanism of accumulation or tolerance to it [41].
Furthermore, the indication of the area in the north in the edaphic models is due to other elements of its composition.For the studied species, Mg, K, N, silt, and clay were the most important elements in determining the distribution model of its edaphic variables.K has a positive correlation with the development and growth of A. glaziovii and can promote its germination.At the same time, N in large concentrations decreases the germination rate of the species [41].Sand and clay contents are directly related to the water retention capacity of canga soils [64].The soil of the northern region has high Mg levels, oscillating between high and low K, which would contradict the preference found in our results.However, it has high values of silt, with 50% (high) of its composition being clay and 20% (low) being sand, and the pH of the area is slightly acidic, between 5.8 and 6.6, which explains the indication of this area in the model [65].
The overlapping of existing mining titles in the Q.F.indicates that the remaining areas of the canga are entirely compromised by mining, and mining activity tends to always increase and surpass past production levels.In addition, the mining lobby has considerable influence in Brazil country [66], and the Brazilian Government expects an increase in mineral production by 2030 [67].This will increase the pressure to exploit mineral reserves further, which increases threats to the existence of A. glaziovii.Mining has already been identified as the leading cause behind the irreversible loss of 40% of canga areas in the last 40 years [17].For the next 50 years, forecasts indicate a catastrophic scenario for cangas [2].Moreover, climate change could result in losses of up to 82%, and mining impacts could result in losses of up to 60% of the remaining areas of the campos rupestres.Therefore, although the forecasts obtained in our study indicate a suitable environment for A. glaziovii's expansion, the speed of loss of the current occurrence areas of the species will hinder its survival until this natural expansion occurs.

Conclusions
The conservation status of the endangered microendemic species Arthrocereus glaziovii could change, as our models indicate an expansion in its future distribution.Unfortunately, this may not be the fate of the species.
Firstly, expansion is hindered because the species' occurrence areas are already under immense pressure from mining activities, and the future predicted occurrence areas of the species are already compromised due to mining.Secondly, climate change predictions specific to the campos rupestres indicate a catastrophic scenario for the cangas, resulting in losses of up to 82% of their areas [2].The risk of species extinction is further exacerbated when considering that the loss and fragmentation of suitable areas could reduce genetic diversity within the species, making survival more difficult in the face of detected environmental changes.Additionally, the limited seed dispersal distance observed in endemic species from the campos rupestres could worsen the survival of A. glaziovii, which, for this species, could be more severe considering that germination rates are consistently below 50% [41,60].
The preservation of both the species and the ecosystems is urgent and imperative.We recommend community involvement and awareness for effective conservation by disseminating existing studies on the campos rupestres and its cacti.Additionally, the establishment of conservation units capable of safeguarding the areas and individuals of A. glaziovii from existing threats is paramount.Finally, further studies to aid the management and propagation of the species are urgently required, such as studies focused on the genetic diversity of its remaining populations, identification of seed propagation matrices, and restoration of degraded habitats to prevent the imminent extinction of this cactus.

Figure 1 .
Figure 1.Regions of occurrence of cangas in the Minas Gerais state.Source: Pristino Institue.

Figure 1 .
Figure 1.Regions of occurrence of cangas in the Minas Gerais state.Source: Pristino Institue.

Figure 2 .
Figure 2. The location of the studied species, Arthrocereus glaziovii (red circles indicate the occurre locations used for ecological niche modeling).

Figure 2 .
Figure 2. The location of the studied species, Arthrocereus glaziovii (red circles indicate the occurrence locations used for ecological niche modeling).

Figure 3 .
Figure 3. Arthrocereus glaziovii prediction models for the climatic variables of the algorithms: (A) Bioclim, (B) Domain, (C) GLM, (D) R.F., (E) MaxEnt, (F) average between all algorithms, (G) potential geographic distribution of the species generated from the average of the algorithms.

Figure 4 .
Figure 4. Arthrocereus glaziovii prediction models for the edaphic variables of the algorithms: (a) Bioclim, (b) Domain, (c) GLM, (d) R.F., (e) MaxEnt, (f) average between all algorithms, (g) potential geographic distribution of the species generated from the average of the algorithms.

Figure 4 .
Figure 4. Arthrocereus glaziovii prediction models for the edaphic variables of the algorithms: (A) Bioclim, (B) Domain, (C) GLM, (D) R.F., (E) MaxEnt, (F) average between all algorithms, (G) potential geographic distribution of the species generated from the average of the algorithms.

Figure 5 .
Figure 5. Mining titles currently licensed in the Iron Quadrangle (Q.F.) region.Source: Pr Institute.

Author Contributions:
Conceptualization, P.d.A.M. and M.d.V.B.; Data curation, A.P.; Formal analysis, A.P. and M.d.V.B.; All authors were involved in the Investigation, Methodology, Validation, Writing-original draft; Visualization, A.P.; Supervision, P.d.A.M. and M.d.V.B.; Writing-review & editing, P.d.A.M. and M.d.V.B.All authors have read and agreed to the published version of the manuscript.

Table 1 .
Climatic and edaphic variables used as predictors in ENM, in order of importance.