ABSTRACT
In this work, the prediction of the distribution of M. urundeuva Fr. All. was performed based on the region of natural occurrence of the species. Its geographic coordinates were obtained from online databases CRIA and SpecialLinks, from scientific articles and fieldwork carried out by Universidade Estadual Paulista (UNESP) in Ilha Solteira, São Paulo, Brazil. M. urundeuva is a native tree species with great potential for commercial use in Brazil. For this purpose, ecological niche modeling was used, with current layers of climate variables and layers prepared for future climate scenarios, according to the 4th Report of the Intergovernmental Panel on Climate Change (AR4/IPCC), using Worldclim data on Brazil. With the Open Modeller and ArcGIS programs, maps were generated to predict its occurrence for the current period and future climate scenarios, made according to the projections of global climate changes. With the projection of increases in temperature and precipitation in the area where the species occurs, it tends to migrate to areas of Brazil where the climate is currently milder, in the south and southeast regions. Due to climatic changes, the species tends to undergo changes in distribution and area size until 2080. It was projected for Caatinga and Pantanal, in both periods, an increase in area, while for the Cerrado, in the first period, the area increased, and, for the second, it decreased. Therefore, according to the results of the maps of future projections for the next decades, it is concluded that there will be changes in the distribution of M. urundeuva, with a significant reduction of the potential area of occurrence in the region.
Keywords:
Brazilian biomes; Climatic changes; Ecological modeling niche
RESUMO
Neste trabalho, foi realizada a previsão da distribuição de M. urundeuva Fr. All., com base na região de ocorrência natural da espécie. Suas coordenadas geográficas foram obtidas nas bases de dados online CRIA e SpecialLinks, a partir de artigos científicos e trabalhos de campo realizados pela Universidade Estadual Paulista (UNESP) em Ilha Solteira, São Paulo, Brasil. M. urundeuva é uma espécie arbórea nativa com grande potencial para uso comercial no Brasil. Para tanto, foi utilizada a modelagem de nicho ecológico, com camadas atuais de variáveis climáticas e camadas preparadas para cenários climáticos futuros, de acordo com o 4º Relatório do Painel Intergovernamental sobre Mudanças Climáticas (AR4/IPCC), utilizando dados do Worldclim no Brasil. Com os programas Open Modeller e ArcGIS, foram gerados mapas para prever sua ocorrência para o período atual e para cenários climáticos futuros, feitos de acordo com as projeções de mudanças climáticas globais. Com a projeção de aumentos de temperatura e precipitação na região onde ocorre a espécie, ela tende a migrar para áreas do Brasil onde o clima é atualmente mais ameno, nas regiões Sul e Sudeste. Devido às mudanças climáticas, a espécie tende a sofrer mudanças na distribuição e tamanho de área até 2080. Foi projetado para a Caatinga e Pantanal, nos dois períodos, um aumento de área, enquanto para o Cerrado, no primeiro período, a área aumentou e para o segundo, diminuiu. Portanto, de acordo com os resultados dos mapas de projeções futuras para as próximas décadas, conclui-se que haverá mudanças na distribuição de M. urundeuva, com redução significativa da área potencial de ocorrência na região.
Palavras-Chave:
Biomas brasileiros; Mudanças climáticas; Nicho de modelagem ecológica
1. INTRODUCTION
Myracrodruon urundeuva Fr. All. (Anacardiaceae) is a tropical tree species, popularly known as aroeira, aroeira do sertão, aroeira preta, and is characterized by high basic density (1,19 g cm-3) and wood durability (Lorenzi, 2000Lorenzi H. Árvores brasileiras: manual de identificação e cultivo de plantas arbóreas nativas do Brasil. Nova Odessa: Plantarum; 2000. 3. ed. v1.). Due to its qualities, it was the subject of intense anthropic exploration, reaching significant importance (Monteiro et al., 2012Monteiro JM, Araújo EL, Amorim ELC, Albuquerque UP. Valuation of the Aroeira (Myracrodruon urundeuva Allemão): perspectives on conservation. Acta Botanica Brasilica. 2012;26(1):125-132. doi: 10.1590/S0102-33062012000100014
https://doi.org/10.1590/S0102-3306201200...
). As a result, it became a member of the Red List of Endangered Species and classified in the Vulnerable Category (Brazil, 2014).
M. urundeuva has a wide natural distribution area compared to most native forest species. It can be found in Argentina, Brazil, Bolivia, and Paraguay. In Brazil, it is widely found in the Midwest, Southeast, Northeast, and, to a lesser extent, in the South Region (Lorenzi, 2000Lorenzi H. Árvores brasileiras: manual de identificação e cultivo de plantas arbóreas nativas do Brasil. Nova Odessa: Plantarum; 2000. 3. ed. v1.; Rizzini, 1971Rizzini CT. Árvores e madeiras úteis do Brasil: manual de dendrologia brasileira. São Paulo: Edusp; 1971. v.1.). Its extensive geographic distribution indicates high genetic diversity thus facilitating its adaptation in different phytoecological regions (Kageyama et al., 2003Kageyama PY, Cunha GC, Barreto KD, Camargo FRA, Sebbenn AM. Genetic diversity and spatial autocorrelation in Ocotea odorífera (Lauraceae) populations. Scientia Forestalis. 2003;64:108-119.) as the Seasonal Forest, both Semi-deciduous and Decidual, Cerrado, Caatinga and Pantanal (Carvalho, 1994Carvalho PER. Espécies florestais brasileiras: recomendações silviculturais, potencialidades e uso da madeira. Colombo: Empresa Brasileira de Agropecuária - Centro Nacional Pesquisas Florestais; 1994. v.1.). In this way, the species presents a characteristic that favors studies related to genetic conservation and climatic changes. Monitoring in the main biomes of the species can be carried out in terms of adaptation, development, population density, and reproductive phenology. However, a priori, the projections of scenarios related to the distribution of the species need to be performed to identify the ecosystems that will require priority for conservation and studies.
Several changes in natural ecosystems in favor of economic development have negatively affected the natural habitat of numerous species in recent years, and consequently the genetic structure of their natural populations. This situation has generated a climate imbalance in several regions, and this has contributed significantly to the changes in the populations during the last years, having each species reacted differently to these climatic effects (Lorenzen et al., 2011Lorenzen ED, Nogués-Bravo D, Orlando L, Weinstock J, Binladen J, Marske, KA, Willerslev E. Species-specific responses of Late Quaternary megafauna to climate and humans. Nature. 2011;479:359-364.). Populations with lower genetic variability will be the most affected, thus characterizing the genetic diversity of populations will also assist in decision making on genetic conservation and use of genetic resources in the face of climate effects. As environmental impacts are evidenced, the demand for a reliable prediction of ecosystem change from these impacts increases (Wang et al., 2012Wang T, Campbell EM, O’neill GA, Aitken SN. Projecting future distributions of ecosystem climate niches: Uncertainties and management applications. Forest Ecology and Management. 2012;279:128– 140. doi: 10.1016/j.foreco.2012.05.034
https://doi.org/10.1016/j.foreco.2012.05...
).
The ecological niche modeling then appears as a promising tool, generating maps of probability of occurrence of a certain species (Terrible et al., 2012Terrible LV, Lima RMS, Nabout JC, Oliveira G, Rangel TF, Rabelo SG, et al. A coupled phylogeographic and species distribution modeling approach recovers the demographic history of a Neotropical seasonally dry forest tree species. Molecular Ecology. 2012;21(23):5845-5863. doi: 10.1111/mec.12071
https://doi.org/10.1111/mec.12071...
) and correlating it with possible environmental and climatic variables, predicting geographic distributions for any individual (Wang et al., 2012Wang T, Campbell EM, O’neill GA, Aitken SN. Projecting future distributions of ecosystem climate niches: Uncertainties and management applications. Forest Ecology and Management. 2012;279:128– 140. doi: 10.1016/j.foreco.2012.05.034
https://doi.org/10.1016/j.foreco.2012.05...
). From the modeling are realized projections, which can be for different ecosystems and be used in diverse ways, orienting in the assembly of methods of adaptation to the climatic changes for the conservation of habitats with species, vegetal or animal, threatened of extinction (Wang et al., 2012Wang T, Campbell EM, O’neill GA, Aitken SN. Projecting future distributions of ecosystem climate niches: Uncertainties and management applications. Forest Ecology and Management. 2012;279:128– 140. doi: 10.1016/j.foreco.2012.05.034
https://doi.org/10.1016/j.foreco.2012.05...
), caused by reduced population size, resulting in loss of genetic diversity (Lovato, 2010Lovato MB. Genetic diversity and structure of the critically endangered tree Dimorphandra wilsonii and of the widespread in the Brazilian Cerrado Dimorphandra mollis: Implications for conservation. Biochemical Systematics and Ecology. 2010;38(1):49-56. doi: 10.1016/j.bse.2009.12.038
https://doi.org/10.1016/j.bse.2009.12.03...
).
Work on climate change versus forest species is still incipient compared to other species. M. urundeuva. Due to its large occurrence and multiple uses, can be considered a reference species in these studies. Thus, the objective of this work was to predict future scenarios for M. urundeuva, using ecological niche models for this purpose, according to the 4th Assessment Report of the Intergovernmental Panel on Climate Change (IPCC, 2007)Intergovernmental Panel on Climate Change — IPCC. Climate change 2007: the physical science basis, summary for policy makers. 2007 [cited 2017 September 15]. Edition 1. Available from: https://www.ipcc.ch/report/ar4/wg1/
https://www.ipcc.ch/report/ar4/wg1/...
(Ar4 / IPCC). Based on this, will be proposed new sites for the conservation and planting of M. urundeuva.
2. MATERIAL AND METHODS
Primary and secondary data were obtained from M. urundeuva points of origin from four different databases, totaling 1,210 points of presence of the species distributed in Brazil.
2.1.Primary database
The primary databases contain georeferenced points of occurrence (latitude and longitude), with data from 19 natural populations, around 30 progenies per population, from the germplasm collection of UNESP, Ilha Solteira, Brazil.
2.2.Secondary database
The used secondary databases containing georeferenced points of occurrence (latitude and longitude) from information provided in the scientific literature. And data stored in Environmental Information Reference Center - CRIA (CRIA, 1999) and the Global Biodiversity Information Facility (GBIF) (GBIF, 2015) database. CRIA (1999) is a biological collections database freely accessible. It can be accessed through the SpeciesLink tool, a species information system (fauna, flora, and microbiota), which gathers historical information from several herbarium present in the whole country, associated to a system of prediction of geographic distribution of species, based on mathematical modeling. Likewise, the GBIF (2015) database collects information on species and their occurrence.
In a first analysis, the consistency and the errors of the data of the places of occurrence of M. urundeuva were verified. Outliers were isolated, which were isolated from the main groupings and did not represent the natural distribution zone of the species. In a second analysis, inverted coordinates (latitude with longitude and vice versa) were corrected and incorrect ones were excluded (occurrence on the sea, for example). There were excluded 295 occurrences, leaving 915 records of occurrence of M. urundeuva in the world. As the modeling of the occurrence prediction was made only for Brazil, the coordinates with the presence register of the species in other countries were removed, leaving 844 records remaining. For this purpose, a geographic information system - GIS was used.
2.3.Determination of bioclimatic variables
The climatic data of the base period (or the "current period") and the future scenarios were obtained from the Worldclim site and maps of the bioclimatic variables were generated with a spatial resolution of approximately 1 km2. The bioclimatic variables were organized in the monthly, seasonal and annual time scale, the main variables of importance for determining the distribution of the species, such as air temperature and rainfall, among others (Kumar and Stohlgren, 2009Kumar S, Stohlgren TJ. Maxent modeling for predicting suitable habitat for threatened and endangered tree Canacomyrica monticule in New Caledonia. Journal of Ecology and Natural Environment. 2009;1(4):94-98.). The current or present period is currently called the base period or reference period, since it deals with a period that, with climate change, has undergone major changes and has no meaning to be called current. Data from the base period consider the average of the climatic variables of the period between 1961 and 1990, according to the World Meteorological Organization (WMO). Future scenarios were projected for 2041-2060 and 2061-2080, according to the IPCC (2007). The Brazilian border was obtained from IBGE (2001).
The base period climate maps and future climate scenarios were elaborated using multiple linear regression, using latitude, longitude, and the digital elevation model (representing altitude), as predictor variables. The maps were made for Brazil, in a 1: 1,000,000 scale. The digital elevation model (DEM) used was the GTOPO30, developed by the United States Geological Survey (USGS, 1999)Survey National Mapping Division — USGS. 1999 [cited 2017 September 18]. Edition 1. Available from: http://edcwww.cr.usgs.gov/landdaac/gtopo30/gtopo30.html
http://edcwww.cr.usgs.gov/landdaac/gtopo...
, based on the Shuttle Radar Topography Mission (SRTM) (Farr and Kobrick, 2000Farr TG, Kobrick M. Shuttle Radar Topography Mission produces a wealth of data. American Geophysical Union. 2000;81(2):583-585.).
For the modeling of occurrence prediction, 19 bioclimatic variables were selected, including minimum and maximum temperatures (°C) and rainfall (mm) (Wrege et., 2017Wrege MS, Fritzons E, Soares MTS, Bognola IA, Sousa VA, Sousa LP, Ferrer R. S. Distribution and habitat of Brazilian-pine according to global climate change. Brazilian Journal of Forestry Research. 2017;37(91):331–346. doi: 10.4336/2017. pfb.37.91.1413
https://doi.org/10.4336/2017. pfb.37.91....
), described below: (Bio1) Annual mean temperature; (Bio2) Monthly mean of the daily temperature variation (maximum temp -minimum temp); (Bio3) Isothermality (Bio2/Bio7) (* 100); (Bio4) Seasonality of temperature (standard deviation *100); (Bio5) Maximum temperature in the hottest month; (Bio6) Minimum temperature in the coldest month; (Bio7) Annual temperature variation (Bio5-Bio6); (Bio8) Mean temperature in the wettest quarter; (Bio9) Mean temperature in the driest quarter; (Bio10) Mean temperature in the hottest quarter; (Bio11) Mean temperature in the coldest quarter; (Bio12) Cumulative rainfall in the year; (Bio13) Precipitation accumulated in the wettest month; (Bio14) Precipitation accumulated in the driest month; (Bio15) Seasonality of rainfall (coefficient of variation); (Bio16) Cumulative rainfall in the wettest quarter; (Bio17) Cumulative rainfall in the driest quarter; (Bio18) Cumulative rainfall in hottest quarter; (Bio19) Cumulative rainfall in the coldest quarter.
2.4.Prediction of occurrence of M. urundeuva
The climate data used in the work were organized, compiled, and consisted. The base period averages of the climate variables listed in Table 1 were calculated. Based on the trends in the 4th Assesment Report of the Intergovernmental Panel on Climate Change (AR4 / IPCC) (IPCC, 2007)Intergovernmental Panel on Climate Change — IPCC. Climate change 2007: the physical science basis, summary for policy makers. 2007 [cited 2017 September 15]. Edition 1. Available from: https://www.ipcc.ch/report/ar4/wg1/
https://www.ipcc.ch/report/ar4/wg1/...
, future climate scenarios were expected. The selected scenarios were A2 and B1, with A2 being the most pessimistic scenario, with the maintenance of current GHG emissions standards, and B1, a scenario of lower emissions or less pessimistic scenario (Nakicenovic et al., 2000Nakicenovic N, Alcamo J, Grubler A. Special report on emissions scenarios (SRES), a special report of Working Group III of the intergovernmental panel on climate change. Cambridge: University Press; 2000. ISBN 9291691135.). Climate projections for the coming decades were made for 2041-2060 and 2061-2080 periods.
The prediction of geographic distribution was made by six algorithm models: Bioclim, Climate Space Model, Envelope Score, Maximum Entropy, Niche Mosaic, and Enviromental Distance (Figure 1). The modeling of the future prediction of occurrence of the species occurred with the use of Open Modeller software (Muñoz et al., 2011Muñoz MES, Giovanni RS, Siqueira MF, Sutton T, Brewer P, Pereira RS, et al. OpenModeller: a generic approach to species potential distribution modelling. GeoInformatica. 2011;15(1):111-135.). This program works with geographic distribution data of species (latitude and longitude) and with maps or environmental layers (climate, soil and relief), composing a mathematical system of prediction of species distribution.
The environmental variables were the same for the base period and for the future. The output maps of the models were transformed into numerical values, varying between 0 and 1. Each pixel of the map came to represent a value, 0 or 1, 0 being related to the no possibility of occurrence of M. urundeuva and 1 the maximum possibility of occurrence. The maps formed in the Open Modeller in ASCII text format (American Standard Code for Information Interchange), containing these binary values, were imported into ArcGIS and transformed into 'raster' format. Classes were then created with a gradient varying from '0' to '1', representing from the zones with no possibility of occurrence until the zones with maximum possibility of occurrence, respectively, for the development of the species.
The evaluation of the quality of the adjusted models was done by calculating the area under the curve (AUC), obtained from the integration of the Receiver Operating Characteristic (ROC) curve. The maximum AUC value is theoretically 1.0 and indicates perfect discrimination, while values lower than 0.5 denote poor modeling performance. The Jackknife test was applied to diagnose the relative contribution of each bioclimatic variable.
Thus, in this study, through the prediction of occurrence of the species, sites were identified where conservation programs of natural populations of M. urundeuva should be prioritized and where there is aptitude for new plantings, considering the climate of the base period and that of the next decades, designed according to the analysis of global climate change.
3. RESULTS
From the used databases, it was possible to obtain 915 points of occurrence of M. urundeuva in the world, all concentrated in South America (Figure 2A) and 844 points within of the Brazilian territory (Figure 2B), after the elimination of discrepant points, with location error or located outside Brazil. The mapping of the current distribution of M. urundeuva (Figure 3), based on these points of presence of the species, was significant for all the models used (p <0.001). Among the generated models, the most representative of the distribution was selected. The most similarity to the species distribution was Environmental Distance, with AUC of 1.00. The model used allowed to express the potential occurrence of the species in an area that covers the domain of four Brazilian biomes, including the Pantanal, Cerrado, Mata Atlântica, and Caatinga limits. The projections for future climatic scenarios (Figure 4) point to the reduction of areas suitable for the development of the species in the two evaluated periods (2041-2060 and 2061-2080), with the most sensitive areas species are concentrated in the northern latitudinal limits of the Cerrado and Caatinga biomes.
The models generated indicate the fragility of M. urundeuva populations in the Caatinga, Pantanal, and Cerrado biomes, and the change for each period (2041-2060 and 2061-2080) was of a total increase in scenario A2, with Caatinga in both periods, an increase of 0.05% and 0.07%, while for the Cerrado an estimated loss of 5.6% for the period 2041-2060 and an increase of 1.2% for the second period, and in the biome Pantanal area increased by 0.2% in both periods. For scenario B1, loss of area in the Caatinga of 1.10% in the period of 2041-2060 and of 10.0% in the period of 2061-2080, while for the Pantanal the first period was an increase of 0.2% in the area and loss of 99% for the second period studied and, for the Cerrado biome, there was a loss of area in both periods of 4% and 38%, respectively (Figures 5A and 5B). Therefore, the natural populations located in these biomes should be prioritized for sampling, aiming the implantation of collections of germplasm, breeding, and commercial plantations.
4.DISCUSSION
The Caatinga is among the most sensitive ecosystems to climate change (Seddon et al., 2016Seddon AW, Macias FM, Long PR, Benz D, Willis KJ. Sensitivity of global terrestrial ecosystems to climate variability. Nature. 2016;531:229-232.). Vulnerable areas to population reduction have also been found to coincide with regions suffering from severe anthropogenic pressure from land use change. In the Cerrado biome, in its northern limit of distribution of the species, great loss of area with potential for development of the species in the next decades, coinciding with the Arch of Deforestation, could occur. In the state of Mato Grosso, most deforestation has occurred in these transitional areas between biomes (Fausto et al., 2016Fausto MA, Angelini LP, Marques HO, Filho AS, Machado NG, Biudes MS. Impact of land-use change in the net radiation of the Cerrado of the southern Mato Grosso. Revista Ambiente e Água. 2016;1(1):350-361. doi: 10.4136/ambi-agua.1843
https://doi.org/10.4136/ambi-agua.1843...
) and in non-flooded areas of the Pantanal, under pressure from agricultural occupation (Azevedo and Saito, 2013Azevedo AA, Saito CH. O perfil dos desmatamentos em Mato Grosso, após implementação do licenciamento ambiental em propriedades rurais. Cerne. 2013;19(1):11-122. doi: 10.1590/S0104-77602013000100014
https://doi.org/10.1590/S0104-7760201300...
).
The Pantanal Matogrossense is considered one of the most important Brazilian biomes, containing one of the largest continuous floodplains in the world, whose ecological functioning is dependent on the complex hydroclimatological dynamics of the region. Although there are uncertainties regarding the results obtained by climate models in relation to the behavior of the hydrological cycle in relation to future scenarios in the region (AR4/IPCC) (Marengo and Valverde, 2007Marengo JA, Valverde MC. Caracterização do clima no Século XX e Cenário de Mudanças de clima para o Brasil no Século XXI usando os modelos do IPCC-AR4. Revista Multiciência. 2007;8(1):5-28.), there are significant microclimatic changes in the conversion of forested areas to pasture in this biome, with potential changes in the regime of rainfall, temperature, and energy balance (Biudes et al., 2012Biudes MS, Nogueira JS, Dalmagro HJ, MAchado NG, Danelichen VHM, Souza MC. Change in the microclimate due to the conversion of a cambará forest to pasture in the northern Pantanal. Revista de Ciências Agro-Ambientais. 2012;10(1):61-68.). These anthropic alterations, caused by the change in the use of the soil, can potentiate the inherent risks of the climate changes on the native vegetation in the Brazilian biomes, as pointed out by Aleixo et al. (2010)Aleixo A, Albernaz AL, Grelle CEV, Vale MM, Rangel TF. Mudanças Climáticas e a Biodiversidade dos Biomas Brasileiros: Passado, Presente e Futuro. Natureza & Conservação. 2010;8(2):194-196. doi: 10.4322/natcon.00802016
https://doi.org/10.4322/natcon.00802016...
, for different terrestrial ecosystems.
Potential distribution modeling is an alternative tool for mapping potential areas for adaptation of species, using a smaller number of variables that would be required by the zoning method, making possible the extrapolation of species occurrence projections in future scenarios, according to global climate change (Bader et al., 2008Bader D, Covey C, Gutowski W, Held I, Kenneth K, Miller R, et al. [Internet]. 2008;3(1) [cited 2017 April 27]. Available from: https://dr.lib.iastate.edu/ entities/publication/40b61a66-b5e1-4823-b55a-7502db4c8c93
https://dr.lib.iastate.edu/ entities/pub...
; Melo et al., 2015Melo LC, Sanqueta CR, Corte APD, Virgens Filho JS. Cenários Climáticos Futuros Para O Paraná: Oportunidades Para O Setor Florestal. Revista Brasileira de Climatologia.2015;11(16):120-131. doi: 10.5380/abclima.v16i0.41149
https://doi.org/10.5380/abclima.v16i0.41...
). However, it has limitations and should be used in conjunction with other tools to aid decision-making (Garcia et al., 2014Garcia LG, Ferraz SFB, Alvares CA, Ferraz KMOM, Higa RCV. Modeling suitable climate for Eucalyptus grandis under future climates scenarios in Brazil. Scientia Forestalis. 2014;42(104):493-501.). An example is the various tests of provenances and/or progenies that present high values of genetic variability, which vary according to the place of origin. These also indicate a good adaptation to the climatic conditions, evidenced by the high survival rate and growing development of the M. urundeuva individuals (Guerra et al., 2009Guerra CRSB, Moraes MLT, Silva CLSP, Canuto DSO, Andrade JAC, Freitas ML M, Sebbenn AM. Selection strategies within progenies of two Myracrodruon urundeuva Fr. All. populations. Scientia Forestalis. 2009;37(1):79-87.; Tung et al., 2010Tung ESC, Freitas MLM, Florsheim SMB, Lima IL, Longui EL, Santos FW, et al. Genetic variation for silvicultural and anatomical traits of wood in progenies of Myracrodruon urundeuva (Engler) Fr. Allem. Scientia Forestalis. 2010;38(87):499-508.; Moraes et al., 2012Moraes MA, Filho WVVF, Resende MDV, Silva AM, Manoel ROM, Freitas MLM, et al. Productive, stability and adaptability in progenies of Myracrodruon urundeuva F.F. & M.F. Allemão – Anacardiaceae. Scientia Forestalis. 2012;40(93):69-76.; Bertonha et al., 2016Bertonha LJ, Freitas MLM, Cambuim J, Moraes MLT, Sebbenn AM. Selection of Myracrodruon urundeuva progenies based on phenological traits and growth for reconstitution of legal reserve areas. Scientia Forestalis. 2016;44(109):95-104. doi: 10.18671/scifor.v44n109.09
https://doi.org/10.18671/scifor.v44n109....
; Pupin et al., 2017Pupin S, Freitas MLM, Canuto DSO, Silva AM, Marin ALA, Moraes MLT. Genetic variability and genetic gains in progeniesof Myracrodruon urundeuva Fr. All. in multi-species system. Nativa. 2017;5(1):59-65. doi: 10.5935/2318-7670.v05n01a10
https://doi.org/10.5935/2318-7670.v05n01...
). In this way, the preservation of forest fragments and the conservation of ex situ materials are fundamental to maintain the genetic variability of this species over the years.
However, it is also important to consider that the presented ecological model was generated with a database of insufficient climatic data for the whole national territory, especially in remote areas of difficult access, where a large part of the individuals of this species are concentrated. A more detailed survey of the climate in the areas of occurrence of the M. urundeuva could contribute more efficiently to this study, as well as information on the genetic variability of the populations and, consequently, the potential of their adaptation in each environment.
The interactions, as edaphic conditions, not considered in this work, may also modify the projection of occurrence of the species (Garcia et al., 2014Garcia LG, Ferraz SFB, Alvares CA, Ferraz KMOM, Higa RCV. Modeling suitable climate for Eucalyptus grandis under future climates scenarios in Brazil. Scientia Forestalis. 2014;42(104):493-501.), which may further decrease the areas with potential for its occurrence in the coming decades, since the projections presented here did not consider such restrictions. The main biomes in which the M. urundeuva is found have distinct aspects, so it is of paramount importance a more detailed study regarding the characteristics of each one, in order to analyze the behavior of the species. In this way, it is possible to analyze its adaptation in the biomes, being more successful in the choice of the area for its conservation.
In the Caatinga biome, the annual temperature averages around 25°C and 30°C, with little difference between the colder and hotter months. The annual precipitation varies from 300 mm in the central area to 1,000 mm in the transition with other biomes. Although it has very shallow soils, which supply the need of the plant only for a few days, the greatest extension is of deep and well-drained soils, with good water retention (Garglio, 2010Garglio MA. Uso sustentável e conservação dos recursos florestais da caatinga. Brasília: Serviço Florestal Brasileiro; 2010. ISBN 978-85-63269-04-1), which is important, since the dry season in the Caatinga can last seven to eight months (Scariot et al., 2005Scariot A, Souza SJC, Felfili JM. Cerrado e Secas: ecologia, biodiversidade e conservação; Brasília: Ministério do Meio Ambiente; 2005. v.1.). The Cerrado biome, as well as the Caatinga, has several textures and soil depth. The accumulated rainfall in one year varies between 600 mm and 800 mm in the transition with the Caatinga and from 2,000 mm to 2,200 mm at the border with the Amazon. On average, it has accumulated rainfall in the year from 1,200 mm to 1,800 mm and a period of six months of drought (Scariot et al., 2005Scariot A, Souza SJC, Felfili JM. Cerrado e Secas: ecologia, biodiversidade e conservação; Brasília: Ministério do Meio Ambiente; 2005. v.1.). The average temperature varies between 22°C and 23°C (Klein, 2002Klein AL. Eugen Warming e o Cerrado brasileiro: um século depois. São Paulo: Unesp; 2002.). In the Pantanal, the average annual temperature is between 23°C and 25°C, with a rainfall index of 1.110 mm accumulated in the year and seven months of drought (Garcia, 1986Garcia EAC. Estudo técnico econômico da pecuária bovina de corte do Pantanal Mato-Grossense. Corumbá: Embrapa-CPAP;1986. v.1.; Narcuzzo et al., 2011Narcuzzo FFN, Rocha HM, Melo DCR. Rainfall mapping of the Pantanal biome in the state of Mato Grosso. Geoambiente. 2011;16(1):66-84.). In this biome, the soils are considered to be more fragile (Beirigo et al., 2011Beirigo RM, Torrado PV, Stape JL, Couto EG, Andrade GRP. Solos da reserva particular do patrimônio natural SESC Pantanal. Rio de Janeiro: SESC; 2011. v.7.; Batista et al., 2014Batista MA, Paiva DW, Marcolino A. Solos para todos: perguntas e respostas. Rio de Janeiro: Embrapa Solos; 2014. v.1.).
According to Costa et al. (2015)Costa AS, Freire ALO, Bakke IA, Pereira FHF. Respostas Fisiológicas e Bioquímicas de Plantas de Aroeira (Myracrodruon urundeuva Allemão) ao déficit hídrico e posterior recuperação. Irriga. 2015;20(4):705-717. doi: 10.15809/ irriga.2015v20n4p705
https://doi.org/10.15809/ irriga.2015v20...
M. urundeuva, which is tolerant to drought, has a reduced photosynthetic rate with water deficit, but it is a species that recovers rapidly when drought ceases, as it has great plasticity, adapted to the conditions from the Atlantic Forest to the Cerrado, which leads us to believe that the species will adapt, even if slowly, to climate change, so there is a decrease in the area of survival, but not extinction.
The fact that M. urundeuva shapes more efficiently the prediction conditions considered more pessimistic may be directly linked to global warming. According to Nunes (2008)Nunes YRF, Fagundes M, Almeida HS, Veloso MDM. Ecological aspects of Aroeira (Myracrodruon urundeuva Allemão - Anacardiceae): phenology and seed germination. Revista Árvore. 2008;32(2):233-243. doi: 10.1590/S0100-67622008000200006
https://doi.org/10.1590/S0100-6762200800...
, when there is a decrease in temperature, there is a greater fall of leaves in the individuals and the ripening of the fruit is affected, making it unfeasible, since it is carried out between the months of August and November, higher temperature.
Another point to be observed regarding the decrease of the area occupied by the species is the pollinator, which has its flight affected by the climatic conditions. Honey bee, the main pollinator of the M. urundeuva, makes its displacement in the morning, usually between 8 and 11:00 AM. With the increase of temperature, this time can be altered, harming the pollination of the individuals. The different responses of plants and insects to climate changes can bring changes in the time and space of activity of bees, respectively related to their phenology and distribution, bringing serious demographic consequences for tree species, including M. urundeuva.
Pollinating agents are the key to the planet's biodiversity, providing vital ecosystem services to tree species with the pollination of flowers and, consequently, the exchange of genomes between individuals and populations. About 75% of species on the planet are dependent on this service, including M. urundeuva.
Climate changes related to the increase in droughts and air temperature can also affect the flowers, increasing the fall and reducing the production of nectar and the protein content in the pollen. These changes may result in less insect visitation to flowers and reduced pollen deposition, reducing seed production and gene exchange between different populations of M. urundeuva. This species seeds quality is related to the maximum, average, and minimum temperature, average and minimum humidity, and precipitation (Oliveira et al., 2020Oliveira GM, Angelotti F, Marques EJN, Silva FFS, Pelacani CR, Dantas BF. Climate and the Myracrodruon urundeuva Allemão seed production. Revista Brasileira de Geografia Física. 2020;13(07):3689-3697.). However, these climate variables during the different phenological phases of M. urundeuva affect the physiological quality of the seeds, and, in climate change scenarios, there will be a reduction in the seed production of this species (Oliveira et al., 2020Oliveira GM, Angelotti F, Marques EJN, Silva FFS, Pelacani CR, Dantas BF. Climate and the Myracrodruon urundeuva Allemão seed production. Revista Brasileira de Geografia Física. 2020;13(07):3689-3697.).
The results indicate the need for population monitoring in the coming decades, especially in areas where climatic aptitude may change with climatic changes (Garcia et al., 2014Garcia LG, Ferraz SFB, Alvares CA, Ferraz KMOM, Higa RCV. Modeling suitable climate for Eucalyptus grandis under future climates scenarios in Brazil. Scientia Forestalis. 2014;42(104):493-501.), as well as the need for a better understanding of the direct and indirect effects of global climate change on the occurrence and development of the species, increased mortality and species vulnerability to pest and disease attack.
The first populations to be affected will be those located in the border areas of the species' distribution zone, where its genetic diversity is greater, due to climate change between the area where the species is registered and the area where it no longer occurs. The climate is the main factor related to the occurrence of species and vegetation represents the expression of the climate of a place.
Ecological niche modeling, together with field assessments, can greatly contribute to improving knowledge on species distribution and distribution trends in the future, based on global climate change, including collaborating to define areas where they can be prioritized field trips to improve population sampling, especially in the Midwest and Northeast regions, at the Northern border of the species (in the vicinity of the Deforestation Arc) (Moscoso et al., 2013Moscoso V, Albernaz AL, Salomão RP. Niche modelling for twelve plant species (six timber species and six palm trees) in the Amazon region, using collection and field survey data. Forest ecology and management. 2013;310:652-662. doi: 10.1016/j. foreco.2013.08.064
https://doi.org/10.1016/j. foreco.2013.0...
). This area is known as the region where there is great advance of the agricultural frontier and, thus, great deforestation. It comprises the regions between Maranhão and Rondônia (Ipam, 2015; Vieira et al., 2005Vieira ICG. Silva JMC, Toledo PM. Estratégias para evitar a perda de biodiversidade na Amazônia. Estudos Avançados. 2005;19(54):153–164.). Another factor responsible for the great loss of original vegetation cover in the region is cattle breeding (Carvalho et al., 2016Carvalho TS, Magalhães AS, Domingues EP. Desmatamento e a contribuição econômica da floresta na Amazônia. Estudos Econômicos. 2016;46(2):499–531. doi: 10.1590/0101-416146288tae
https://doi.org/10.1590/0101-416146288ta...
), which is responsible for soil erosion, silting of watercourses, and loss of drinking water quality (Capoane et al., 2016Capoane V, Santos DR, Kuplich TM. Efeito do uso da terra e relevo na distribuição espacial e vertical de carbono do solo em uma bacia hidrográfica agrícola. Revista Brasileira de Geografia Física. 2016;9(06):1629-1643.).
5. CONCLUSIONS
1. According to the results of the layers of future projections scenarios for the next decades, there are changes in the distribution of M. urundeuva in the climatic scenarios A2 and B1, with a significant reduction in the potential area of occurrence in the Northern latitudinal limits, mainly where the Arc of Deforestation;
3. The indication of strategic areas for the rescue of genetic material, such as the borders of the distribution area of M. urundeuva, may help to conserve the species.
4. With the increases in temperature and rainfall in the area where the species occurs, it tends to migrate to areas of Brazil where the climate is currently milder, in the south and southeast regions.
6. REFERENCES
- Aleixo A, Albernaz AL, Grelle CEV, Vale MM, Rangel TF. Mudanças Climáticas e a Biodiversidade dos Biomas Brasileiros: Passado, Presente e Futuro. Natureza & Conservação. 2010;8(2):194-196. doi: 10.4322/natcon.00802016
» https://doi.org/10.4322/natcon.00802016 - Azevedo AA, Saito CH. O perfil dos desmatamentos em Mato Grosso, após implementação do licenciamento ambiental em propriedades rurais. Cerne. 2013;19(1):11-122. doi: 10.1590/S0104-77602013000100014
» https://doi.org/10.1590/S0104-77602013000100014 - Bader D, Covey C, Gutowski W, Held I, Kenneth K, Miller R, et al. [Internet]. 2008;3(1) [cited 2017 April 27]. Available from: https://dr.lib.iastate.edu/ entities/publication/40b61a66-b5e1-4823-b55a-7502db4c8c93
» https://dr.lib.iastate.edu/ entities/publication/40b61a66-b5e1-4823-b55a-7502db4c8c93 - Batista MA, Paiva DW, Marcolino A. Solos para todos: perguntas e respostas. Rio de Janeiro: Embrapa Solos; 2014. v.1.
- Beirigo RM, Torrado PV, Stape JL, Couto EG, Andrade GRP. Solos da reserva particular do patrimônio natural SESC Pantanal. Rio de Janeiro: SESC; 2011. v.7.
- Bertonha LJ, Freitas MLM, Cambuim J, Moraes MLT, Sebbenn AM. Selection of Myracrodruon urundeuva progenies based on phenological traits and growth for reconstitution of legal reserve areas. Scientia Forestalis. 2016;44(109):95-104. doi: 10.18671/scifor.v44n109.09
» https://doi.org/10.18671/scifor.v44n109.09 - Biudes MS, Nogueira JS, Dalmagro HJ, MAchado NG, Danelichen VHM, Souza MC. Change in the microclimate due to the conversion of a cambará forest to pasture in the northern Pantanal. Revista de Ciências Agro-Ambientais. 2012;10(1):61-68.
- Centro Nacional de Conservação da Flora - BRASIL. Projeto Lista Vermelha: lista vermelha. 2014. [cited 2018 January 15]. Edition 1. Available from: http://cncflora.jbrj.gov.br/portal/pt-br/listavermelha
» http://cncflora.jbrj.gov.br/portal/pt-br/listavermelha - Capoane V, Santos DR, Kuplich TM. Efeito do uso da terra e relevo na distribuição espacial e vertical de carbono do solo em uma bacia hidrográfica agrícola. Revista Brasileira de Geografia Física. 2016;9(06):1629-1643.
- Carvalho PER. Espécies florestais brasileiras: recomendações silviculturais, potencialidades e uso da madeira. Colombo: Empresa Brasileira de Agropecuária - Centro Nacional Pesquisas Florestais; 1994. v.1.
- Carvalho TS, Magalhães AS, Domingues EP. Desmatamento e a contribuição econômica da floresta na Amazônia. Estudos Econômicos. 2016;46(2):499–531. doi: 10.1590/0101-416146288tae
» https://doi.org/10.1590/0101-416146288tae - Costa AS, Freire ALO, Bakke IA, Pereira FHF. Respostas Fisiológicas e Bioquímicas de Plantas de Aroeira (Myracrodruon urundeuva Allemão) ao déficit hídrico e posterior recuperação. Irriga. 2015;20(4):705-717. doi: 10.15809/ irriga.2015v20n4p705
» https://doi.org/10.15809/ irriga.2015v20n4p705 - Sistema de Informação Ambiental do Programa Biota — CRIA-BDT. Centro de Referência em Informação Ambiental, BDT e Instituto de Computação. 1999. [cited 2017 September 25]. Edition 1. Available from: http://www.biotasp.org.br/sia/
» http://www.biotasp.org.br/sia/ - Farr TG, Kobrick M. Shuttle Radar Topography Mission produces a wealth of data. American Geophysical Union. 2000;81(2):583-585.
- Fausto MA, Angelini LP, Marques HO, Filho AS, Machado NG, Biudes MS. Impact of land-use change in the net radiation of the Cerrado of the southern Mato Grosso. Revista Ambiente e Água. 2016;1(1):350-361. doi: 10.4136/ambi-agua.1843
» https://doi.org/10.4136/ambi-agua.1843 - Garcia LG, Ferraz SFB, Alvares CA, Ferraz KMOM, Higa RCV. Modeling suitable climate for Eucalyptus grandis under future climates scenarios in Brazil. Scientia Forestalis. 2014;42(104):493-501.
- Garcia EAC. Estudo técnico econômico da pecuária bovina de corte do Pantanal Mato-Grossense. Corumbá: Embrapa-CPAP;1986. v.1.
- Garglio MA. Uso sustentável e conservação dos recursos florestais da caatinga. Brasília: Serviço Florestal Brasileiro; 2010. ISBN 978-85-63269-04-1
- Global Biodiversity Information Facility — GBIF. GBIF Occurrence. 2015 [cited 2017 September 10]. Edition 1. Available from: http://doi.org/10.15468/dl.xpdycr
» http://doi.org/10.15468/dl.xpdycr - Guerra CRSB, Moraes MLT, Silva CLSP, Canuto DSO, Andrade JAC, Freitas ML M, Sebbenn AM. Selection strategies within progenies of two Myracrodruon urundeuva Fr. All. populations. Scientia Forestalis. 2009;37(1):79-87.
- Instituto Brasileiro de Geografia e Estatística — IBGE. Malha municipal digital do Brasil. 2001 [cited 2017 September 02]. Edition 1. Available from: https://www.ibge.gov.br/geociencias/organizacao-do-territorio/malhas-territoriais/15774-malhas.html
» https://www.ibge.gov.br/geociencias/organizacao-do-territorio/malhas-territoriais/15774-malhas.html - Instituto de Pesquisa Ambiental da Amazônia — IPAM. Arco do desmatamento. 2015 [cited 2017 September 02]. Edition 1. Available from: https://ipam.org.br/glossario/arco-do-desmatamento/
» https://ipam.org.br/glossario/arco-do-desmatamento/ - Intergovernmental Panel on Climate Change — IPCC. Climate change 2007: the physical science basis, summary for policy makers. 2007 [cited 2017 September 15]. Edition 1. Available from: https://www.ipcc.ch/report/ar4/wg1/
» https:// www.ipcc.ch/report/ar4/wg1/ - Kageyama PY, Cunha GC, Barreto KD, Camargo FRA, Sebbenn AM. Genetic diversity and spatial autocorrelation in Ocotea odorífera (Lauraceae) populations. Scientia Forestalis. 2003;64:108-119.
- Klein AL. Eugen Warming e o Cerrado brasileiro: um século depois. São Paulo: Unesp; 2002.
- Kumar S, Stohlgren TJ. Maxent modeling for predicting suitable habitat for threatened and endangered tree Canacomyrica monticule in New Caledonia. Journal of Ecology and Natural Environment. 2009;1(4):94-98.
- Lorenzen ED, Nogués-Bravo D, Orlando L, Weinstock J, Binladen J, Marske, KA, Willerslev E. Species-specific responses of Late Quaternary megafauna to climate and humans. Nature. 2011;479:359-364.
- Lorenzi H. Árvores brasileiras: manual de identificação e cultivo de plantas arbóreas nativas do Brasil. Nova Odessa: Plantarum; 2000. 3. ed. v1.
- Lovato MB. Genetic diversity and structure of the critically endangered tree Dimorphandra wilsonii and of the widespread in the Brazilian Cerrado Dimorphandra mollis: Implications for conservation. Biochemical Systematics and Ecology. 2010;38(1):49-56. doi: 10.1016/j.bse.2009.12.038
» https://doi.org/10.1016/j.bse.2009.12.038 - Marengo JA, Valverde MC. Caracterização do clima no Século XX e Cenário de Mudanças de clima para o Brasil no Século XXI usando os modelos do IPCC-AR4. Revista Multiciência. 2007;8(1):5-28.
- Melo LC, Sanqueta CR, Corte APD, Virgens Filho JS. Cenários Climáticos Futuros Para O Paraná: Oportunidades Para O Setor Florestal. Revista Brasileira de Climatologia.2015;11(16):120-131. doi: 10.5380/abclima.v16i0.41149
» https://doi.org/10.5380/abclima.v16i0.41149 - Monteiro JM, Araújo EL, Amorim ELC, Albuquerque UP. Valuation of the Aroeira (Myracrodruon urundeuva Allemão): perspectives on conservation. Acta Botanica Brasilica. 2012;26(1):125-132. doi: 10.1590/S0102-33062012000100014
» https://doi.org/10.1590/S0102-33062012000100014 - Moraes MA, Filho WVVF, Resende MDV, Silva AM, Manoel ROM, Freitas MLM, et al. Productive, stability and adaptability in progenies of Myracrodruon urundeuva F.F. & M.F. Allemão – Anacardiaceae. Scientia Forestalis. 2012;40(93):69-76.
- Moscoso V, Albernaz AL, Salomão RP. Niche modelling for twelve plant species (six timber species and six palm trees) in the Amazon region, using collection and field survey data. Forest ecology and management. 2013;310:652-662. doi: 10.1016/j. foreco.2013.08.064
» https://doi.org/10.1016/j. foreco.2013.08.064 - Muñoz MES, Giovanni RS, Siqueira MF, Sutton T, Brewer P, Pereira RS, et al. OpenModeller: a generic approach to species potential distribution modelling. GeoInformatica. 2011;15(1):111-135.
- Nakicenovic N, Alcamo J, Grubler A. Special report on emissions scenarios (SRES), a special report of Working Group III of the intergovernmental panel on climate change. Cambridge: University Press; 2000. ISBN 9291691135.
- Narcuzzo FFN, Rocha HM, Melo DCR. Rainfall mapping of the Pantanal biome in the state of Mato Grosso. Geoambiente. 2011;16(1):66-84.
- Nunes YRF, Fagundes M, Almeida HS, Veloso MDM. Ecological aspects of Aroeira (Myracrodruon urundeuva Allemão - Anacardiceae): phenology and seed germination. Revista Árvore. 2008;32(2):233-243. doi: 10.1590/S0100-67622008000200006
» https://doi.org/10.1590/S0100-67622008000200006 - Oliveira GM, Angelotti F, Marques EJN, Silva FFS, Pelacani CR, Dantas BF. Climate and the Myracrodruon urundeuva Allemão seed production. Revista Brasileira de Geografia Física. 2020;13(07):3689-3697.
- Pupin S, Freitas MLM, Canuto DSO, Silva AM, Marin ALA, Moraes MLT. Genetic variability and genetic gains in progeniesof Myracrodruon urundeuva Fr. All. in multi-species system. Nativa. 2017;5(1):59-65. doi: 10.5935/2318-7670.v05n01a10
» https://doi.org/10.5935/2318-7670.v05n01a10 - Rizzini CT. Árvores e madeiras úteis do Brasil: manual de dendrologia brasileira. São Paulo: Edusp; 1971. v.1.
- Scariot A, Souza SJC, Felfili JM. Cerrado e Secas: ecologia, biodiversidade e conservação; Brasília: Ministério do Meio Ambiente; 2005. v.1.
- Seddon AW, Macias FM, Long PR, Benz D, Willis KJ. Sensitivity of global terrestrial ecosystems to climate variability. Nature. 2016;531:229-232.
- Terrible LV, Lima RMS, Nabout JC, Oliveira G, Rangel TF, Rabelo SG, et al. A coupled phylogeographic and species distribution modeling approach recovers the demographic history of a Neotropical seasonally dry forest tree species. Molecular Ecology. 2012;21(23):5845-5863. doi: 10.1111/mec.12071
» https://doi.org/10.1111/mec.12071 - Tung ESC, Freitas MLM, Florsheim SMB, Lima IL, Longui EL, Santos FW, et al. Genetic variation for silvicultural and anatomical traits of wood in progenies of Myracrodruon urundeuva (Engler) Fr. Allem. Scientia Forestalis. 2010;38(87):499-508.
- Survey National Mapping Division — USGS. 1999 [cited 2017 September 18]. Edition 1. Available from: http://edcwww.cr.usgs.gov/landdaac/gtopo30/gtopo30.html
» http://edcwww.cr.usgs.gov/landdaac/gtopo30/gtopo30.html - Vieira ICG. Silva JMC, Toledo PM. Estratégias para evitar a perda de biodiversidade na Amazônia. Estudos Avançados. 2005;19(54):153–164.
- Wang T, Campbell EM, O’neill GA, Aitken SN. Projecting future distributions of ecosystem climate niches: Uncertainties and management applications. Forest Ecology and Management. 2012;279:128– 140. doi: 10.1016/j.foreco.2012.05.034
» https://doi.org/10.1016/j.foreco.2012.05.034 - Wrege MS, Fritzons E, Soares MTS, Bognola IA, Sousa VA, Sousa LP, Ferrer R. S. Distribution and habitat of Brazilian-pine according to global climate change. Brazilian Journal of Forestry Research. 2017;37(91):331–346. doi: 10.4336/2017. pfb.37.91.1413
» https://doi.org/10.4336/2017. pfb.37.91.1413
Publication Dates
-
Publication in this collection
21 Feb 2022 -
Date of issue
2022
History
-
Received
17 May 2021 -
Accepted
23 Dec 2021