1. Introduction
Campo rupestre is a plant physiognomy of mountain tops in eastern Brazil, characterized by a rich mosaic of herbaceous vegetation and sparse shrubs [
1]. The region’s soil is compact, deficient in water and nutrients, and high in heavy metals [
2,
3]. The rocky outcrops in the rupestrian fields, particularly the iron ones (ferruginous outcrops), are rich in endemic plants and substantially contribute to the biodiversity of the domains in which they are inserted, Cerrado and Atlantic Forest [
4,
5]. With the region’s climate being considered as high altitude subtropical (dry winters and wet summers), according to the Köppen climate classification [
6], the rocky fields present daily high thermal amplitudes (minimum of 14.1°C and maximum of 14.1°C). 23.7ºC) and low relative humidity [
7,
8]. Thus, soil, relief, and climate conditions are the main factors for establishing its unique and dominant vegetation cover [
9].
The Cerrado and the Atlantic Forest are considered hotspots for preserving biodiversity as they harbor great biological diversity, which is highly threatened by human exploitation [
10]. This constant threat to rupestrian grasslands, due to the loss and fragmentation of these ecosystems [
11], and their high rates of endemism, makes their conservation extremely important, urgent, and necessary. In addition to the threats mentioned above, we find in the rupestrian fields one of Brazil’s largest iron ore extraction areas, the Quadrilátero Ferrífero (Q.F.), placing mining as the main factor of destruction in this region [
12]. In the years 1960 to 2012, the Q.F. lost about 40% of its total area of cangas (about 100 ferruginous outcrops) to the extraction of iron ore [
5,
12]. Today, nearly 225 remaining ferruginous outcrops and 83% of these show signs of loss or degradation by mining [
5].
The biodiversity of the Q.F. is currently under intense threat, as they make up one of the most important Brazilian mineral provinces in the world in iron ore [
1,
13], and the loss of species from mining activities is marked by a high irreversible environmental impact [
14]. The intense mining activities practiced for centuries on a large scale continue to impact the flora of the rupestrian fields, as they resulted in areas with exposed subsoil or gravel, which, even after decades, no natural recovery was recorded [
15]. Plants from these locations are susceptible to soil removal. They have a low rate of dispersal and slow growth [
16], hindering the natural regeneration of these areas of altered rupestrian fields.
The recurrent mineral exploration in the Q.F. has caused severe impacts on the environment, biodiversity and ecosystem services [
17,
18]. About 89 rare and threatened species in the Q.F.. are plants with high conservation value as they occur exclusively in cangas in the region [
12]. One of these species impacted on the mountain tops of rocky fields by mineral exploration is the microendemic and endangered cactus
Arthrocereus glaziovii, which has suffered from reducing its area of occurrence [
17,
19].
The rupestrian fields make up an ecosystem that is clearly at risk, as they are highly endemic; studies are needed to elucidate the distribution of exclusive, threatened, and rare species for their conservation. The geographic distribution modeling, or niche model, is the most indicated tool to predict the potential distribution of species [
20,
21,
22]. Thus, this work aims to investigate i) the potential distribution of
A. glaziovii, based on abiotic variables and soil elements; and ii) characterize the distribution of mineral titles that may restrict the areas of occurrence of
A. glaziovii.
3. Results
3.1. Selected climatic and edaphic variables
After performing Pearson’s correlation, six of the 20 initial climatic variables were selected (
Figures S1 and S2) 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 VSURF importance values (
Table 1), and solar radiation was not significant (SRAD).
Of the 18 initial soil variables, ten were selected by Pearson’s correlation (
Figure S3) and also 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 others were classified in order of influence (
Table 1).
3.2. Current potential environmental suitability for Arthrocereus glaziovii
All models had the area under the curve (AUC) > 0.98, showing its 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. (
Figure 3a, d;
Figure 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 mountainous areas of the studied site (
Figures 3b, e;
Figure 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 (
Figure 3c;
Figure 4c).
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 most essential variable for the niche was annual precipitation, followed by annual temperature variation, seasonal precipitation, average annual temperature, and temperature of the hottest month.
By Bioclim’s edaphic model, sand and clay showed high levels, from 41 to 53 and 25 to 27, respectively, then silt with 15 to 20, pH at 5, H from 3 to 4, and Mg 0 to 1. By the Domain model, the soil presented similar or close values, sand between 30 and 60, clay between 15 and 38, silt with 10 to 15, pH also in the value of 5, H with 1 to 6, and Mg again 0 a 1. Nitrogen and K had values equal to 0. For MaxEnt, the most critical edaphic variable for the niche was Mg (50%), followed by K, N, clay, and silt. As for R.F., the essential variable for the niche was the concentration of N, followed by silt, clay, K, and Mg. Sand, H, and water pH did not show significant importance in both models.
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 Q.F.) and northern regions. South of the study area (Serra da Mantiqueira) and north (
Figure 3g).
3.3. Mining titles
Overlapping all layers of mineral titles in the Q.F., we observe few remaining areas in this region (
Figure 5). It is clear that the distribution of
A. glaziovii, both actual and potential, is an area of important mining extraction. The distribution of future mines largely overlaps the occurrence and the suitable regions of the potential distribution of the species. Iron ore extraction is the most frequent, followed by gold ore.
4. Discussion
Our results demonstrate that the area of the potential occurrence of A. glaziovii considering both climatic and edaphic variables is more significant and covers areas outside the Q.F.. However, currently licensed and under exploration mining titles threaten the entire identified area, jeopardizing the species’ survival.
Seasonal and annual rainfall, warmest month temperature, 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 cacti’s reproductive period [
42].
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 [
43,
44]. The temperature has a positive influence on the production of flowers and fruits, but rainfall does not influence this factor [
45]. The fruiting of the species takes place in the middle of the rainy season, and both temperature and rainfall positively influence this phenophase [
45]. The increase and decrease in temperature influence the germination of the species, with an optimal temperature between 25 and 30ºC, even in the rainy season to establish seedlings [
1,
46].
The species A. glaziovii currently occurs in the Q.F. region in the center-southeast of the state of M.G., 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 smaller occurrence than the current. The GLM indicated a wide distribution throughout the state of M.G. The cut of the averages of the models stated an occurrence area located to the south (climate model) and to the north of the state of M.G. 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 state of M.G., with a concentration along the Espinhaço Range, to the south and north of the state of M.G.
The area south of M.G. state in the climate models is called Serra da Mantiqueira, and the indication in the model is probably due to 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 [
47]. In Serra da Mantiqueira, the vegetation found is called campo de altitude and grows under granitic outcrops of igneous or metamorphic rocks, which differs from the physiognomies occurring on ferruginous outcrops of the Q.F. and quartzite fields of the Espinhaço Range [
47]. 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 [
47].
Only in the MaxEnt model of climate variables a small area north of the state of M.G. was indicated as a possible area of environmental suitability for the occurrence of
A. glaziovii. The north of M.G. is a transition region of Cerrado and Caatinga with hematite in the soil, with a humid tropical climate of savannah. The driest season coincides with winter, and rainfall is less than 60 mm [
48]. The soil of this indicated region also presents high concentrations of Fe and 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 influence directly on the plant, indicating that the species has some mechanism of accumulation or tolerance to it [
35].
The final indication of the area to the north in the edaphic models is due to other elements of its composition. For the species, Mg, K, N, silt, and clay were the most important elements in determining the distribution model of edaphic variables. K has a positive correlation in 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 [
35]. Sand and clay contents are directly related to the water retention capacity of canga soils [
49]. The northern region soil has high Mg levels, oscillating between high and low K, which would go against the preference found in our results. However, it has high values of silt, with 50% (high) of its composition being clay and 20% (low) of 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 [
50].
The overlapping of existing mining titles in the Q.F. indicates that the remaining areas of canga are entirely compromised by mining, and mining tends always to increase and surpass its past productions. Mining has already been identified as the leading cause of the irreversible loss of 40% of
canga areas in the last 40 years [
12]. For the next 50 years, forecasts indicate a catastrophic scenario for cangas [
51]. Climate change would result in losses of up to 82% and mining impacts in losses of up to 60% of the remaining areas of rupestrian fields.
Although abiotic factors indicate an expansion in the distribution of the species, mining activities compromise this possibility with habitat loss and degradation. Fragmented environments are more vulnerable to the adverse effects of impacts, and a large part of the current area of
canga is constituted by canga fragmented by losses already suffered [
51]. The areas of occurrence and potential distribution of
A. glaziovii are of important mining extraction, around 102
cangas have mining activities close to them, and the indirect impacts of mining reach up to 5 km from the mines [
12,
52,
53]. The expansion of the mining sector in rupestrian fields will continue to generate a domino effect of unprecedented impacts to this ecosystem and its endemic species [
18,
53].
Canga species threatened with extinction by the IUCN, such as
A. glaziovii, could quickly disappear and become extinct if action is not taken to preserve them. The transmission of knowledge and awareness of the importance of these areas and their endemic species, and the magnitude of the losses already suffered by their misuse, is necessary for their effective conservation.