Future uncertainties for the distribution and conservation of Paubrasilia echinata under climate change

ABSTRACT

Paubrasilia echinata is a widely cultivated endangered tree species with small populations restricted to a narrow strip of habitats along the Brazilian coast. The potential impacts of climate change on the distribution of P. echinata have yet to be investigated, and so it remains unknown whether protected areas will ensure the persistence of the species in the future. Here, we estimate the impacts of climate change on the distribution of P. echinata inside and outside of protected areas considering different climate change scenarios and two different sets of presence records: natural distribution and cultivated records. Future scenarios showed a gradual reduction in climatically suitable area both inside and outside of protected areas. Projections indicate a trend for a shift to the highlands of Southeast Brazil, and the loss of several areas throughout the entire distribution of the species. Predicted climatic conditions will be unsuitable for P. echinata inside most protected areas. Information provided here will be relevant in planning future national actions for this species, which is a must to properly protect this long-exploited tree.

Keywords: biodiversity conservation; ecological niche modeling; Fabaceae; species distribution modeling; tree species

Introduction

Climate change is a major threat to biodiversity and natural ecosystems on a global scale (Parmesan & Yohe 2003; Root et al. 2003). Shifts in species phenology and distribution in response to ongoing climate change have been reported for marine, freshwater, and terrestrial groups (Parmesan 2006; Franklin et al. 2016). Several species of vertebrates and plants are expected to lose suitable climate conditions both inside and outside of protected areas (Araújo et al. 2011; Li et al. 2015). Conserving future biodiversity in a changing climate is a complex issue, as some protected areas will become less effective in the future, while others will maintain their conservation value (Thomas & Gillingham 2015; Berteaux et al. 2018). Therefore, predicting future climate scenarios and their effects on the species ranges is crucial to effective conservation strategies and to mitigate the impacts of climate change on biodiversity.

When planning for conservation strategies, it is crucial to take climate change and species redistribution into account, if not, jeopardizing species persistence. Ecological niche modeling is the main tool to predict the impacts of climate change in species distribution (Peterson et al. 2000; 2002; Peterson 2001). With the development of new algorithms (Guisan & Thuiller 2005) and more accurate future climate predictions with multiple scenarios (Moss et al. 2010), we can reduce uncertainties to better support conservation decisions (Elith et al. 2006).

Paubrasilia echinata (Fabaceae), locally known as ‘Pau-Brasil’ is an endemic and endangered tree species with small populations restricted to narrow habitats along the Brazilian coast. Pau-Brasil is estimated to have been widely distributed before the arrival of Europeans, as documented by Américo Vespúcio (Fontana et al. 1994) when he wrote to the king of Portugal saying the Brazilian coast had no profitable resource, but an infinity of Pau-Brasil. Commercial exploitation until 1875 drastically reduced populations (Souza 1939), despite conservation attempts such as the Pau-Brasil statute (Brasil 1605). Since the development of synthetic dyes, the main use of Pau-Brasil shifted to luthiery, as the building material of high-end violin and cello bows (Skeaping 1955; Longui et al. 2010). In 1978, Pau-Brasil was declared as Brazil’s national tree (Lei N° 6.607). Currently it is widely distributed due to conservation efforts. It is cultivated in plant nurseries for urban landscaping and educational purposes. Although the species is found in 11 protected areas, it is unknown whether these areas will maintain viable populations under global climate change, which is the main driver of species redistribution in the XXI century (Thuiller 2004; Diniz-Filho et al. 2009; Thuiller et al. 2011; Loyola et al. 2014; Sales et al. 2017). So far no studies have estimated the potential impacts of climate change on the distribution of P. echinata. This is particularly alarming given the low levels of gene flow and genetic diversity of its populations in nature (Lira et al. 2003; Cardoso et al. 2005).

Our goal is to estimate the impacts of climate changes on the distribution of P. echinata, searching for climatic stable areas where a long-term conservation strategy could be applied, as well as to assess its potential distribution inside of protected areas in different climate change scenarios. We expect a decline of climatically suitable areas for P. echinata both inside and outside of protected areas, as the Atlantic forest habitats where it occurs are highly vulnerable to climate changes (Bellard et al. 2014; Esser et al. 2019). A secondary goal is to estimate the contribution of cultivated records to the models. We expect those records to help find sites where the species will have a long-term persistence and in which the species would never reach despite human mediated migration. It is important to point out that the use of non-natural records is not recommended in ecological niche modeling, since they may include biased information (i.e. environmental information where the species is not actually suitable without human management), and further inferences over this results should be taken with caution.

Materials and methods

Study species

Paubrasilia echinata (Lam.) Gagnon, H.C.Lima & G.P.Lewis is a semideciduous or deciduous medium sized tree (up to 20m tall), armed with prickles, chestnut brown to almost black bark; red heartwood and when injured the trunk exudates a red sap (Gagnon et al. 2016). Paubrasilia is a monospecific genus endemic to Eastern Brazil, with occurrence from the Brazilian states of Rio Grande do Norte to Rio de Janeiro within the Brazilian Atlantic Forest hotspot. Its main habitats are well-drained soils in coastal white-sand woodlands (tall restinga), seasonal forest (semideciduous forest) and rain forest (Gagnon et al. 2016; Lima 2019). This species is listed in CITES Appendix II and is listed as endangered at the national (CNCFlora 2012) and global (IUCN 2019) level due to a population size reduction of 50 % in three generations.

Datasets

Species occurrence data was obtained from SpeciesLink, GBIF, JABOT and Reflora, and were subject to a data cleaning approach, where we deleted miswritten coordinates, as well as records with zero degrees. We also deleted duplicated records within the same cell of a raster with 5-arcmin resolution. Data were then divided in a natural distribution and a cultivated group. The first group was build considering Castro (2002), which after an extensive literature review, based on 23 studies, proposed that the natural distribution of P. echinata is the brazilian coast, extending from 5º39'S in Rio Grande do Norte to 23^oS in Rio de Janeiro. Records selection was made manually using QGIS. The second group comprised all records, including both natural distribution and cultivated records (Fig. 1). Presence records summed a total of 237 sites, of which 180 were considered from its natural distribution (Tabs. S1, S2 in supplementary material). We considered this second approach due to the potential resistance to be unveiled by environmental conditions in those sites, expanding the niche breath of the species. The following routine was then applied for both datasets.

Figure 1
Records of Paubrasilia echinata, segregated in natural occurrences and cultivated occurrences.

Variable selection

Variable selection was conducted in two steps. We firstly pre-selected WorldClim 1.4 (Hijmans et al. 2006) 5-arcmin bioclimatic variables that matched species ecology and habitats where it is found (Esser et al. 2019). Rainforest variables were mean temperature of wettest quarter (BIO 8) and precipitation of warmest quarter (BIO 18), which are both proxy of tropical climate, and annual precipitation (BIO 12), which, together with BIO 18, tracks the lack of dry season. Restinga forest was represented by mean diurnal range, temperature seasonality and precipitation of driest month (BIO 2, 4 and 14), which are respectively proxy of maritimity, altitude and dry. BIO 14 was also representative of seasonal forest, since it is also a dry formation. To represent this last habitat, were selected also mean temperature of driest quarter (BIO 9), which is proxy of dry and heat, and precipitation seasonality (BIO 15), which is proxy of seasonality. The summarized pre-selection is presented in Table 1. Secondly, we conducted a statistical selection using a Variance Inflation Factor (VIF) approach with usdm package from R environment (Naimi et al. 2014; R Development Core Team 2018), selecting variables with correlation between 0.5 and -0.5; and VIF lower than 3. This grants a lower dependency from variables and high predictability. This statistical approach is dependent on the background data, so we masked bioclimatic variables with a one-degree width buffer for each presence coordinate, which better describes the distribution of the species and their surroundings.

Table 1
Pre-selection of variables according to each habitat. Variables marked with a * are those that most appeared in the selection routine, thus used to build models.

Habitat	Variables
Rainforest	*BIO8 = Mean Temperature of Wettest Quarter
	*BIO12 = Annual Precipitation
	BIO18 = Precipitation of Warmest Quarter
Restinga	BIO2 = Mean Diurnal Temperature Range
	*BIO4 = Temperature Seasonality
	BIO14 = Precipitation of Driest Month
Seasonal Forest	BIO5 = Max Temperature of Warmest Month
	BIO9 = Mean Temperature of Driest Quarter
	BIO14 = Precipitation of Driest Month

Ecological niche modeling

Models were generated using 10 algorithms (Artificial Neural Network, ANN, Classification Tree Analysis, CTA, Flexible Discriminant Analysis, FDA, Generalized Additive Model, GAM, Generalized Boosting Model, GBM, Generalized Linear Model, GLM, Surface Range Envelope, SRE, Multiple Adaptive Regression Splines, MARS, Random Forest, RF, and MaxEnt) in three independent runs, selecting in each run 1000 pseudo-absences randomly distributed in the study area, of 5-fold cross-validation, keeping 20 % of records for model evaluation. We calculated True Skill Statistic (TSS) that minimized the error for each model and made a weighted ensemble model summing those with optimum TSS greater than 0.7. Then, models were projected to current and four future scenarios (RCP-2.6, RCP-4.5, RCP-6.0, RCP-8.5; IPCC 2013) to the years of 2050 and 2070. To reduce uncertainties, we used 11 General Circulation Models (BCC-CSM1-1, CCSM4, GISS-E2-R, HadGEM2-AO, HadGEM2-ES, IPSL-CM5A-LR, MIROC-ESM-CHEM, MIROC-ESM, MIROC5, MRI-CGCM3 and NorESM1-M). Finally, we calculated the committee average of projections, which returns both a prediction and a measure of uncertainty, since it is the simple average of binary predictions (i.e. cells with values close to one are those that most models agree with a presence, cells with values close to zero are those that most models agree with an absence and cells with value of 0.5 are those that half of models predict a presence, while the other half predict an absence). Ecological niche models were made using biomod2 package in R.

Area calculation

Area value was obtained with a weighted method, where presence probability is multiplied by cell’s area, then summing all raster values (Esser et al. 2019). This approach results in a conservative area that considers that occupancy is not equal between cells. When we first obtained area values we realized the Amazônia biome had a great bias in results. Given the fact that its soil properties are a key factor to species persistence and we have not considered it in these models, we decided to exclude suitability values in Amazônia in this step. Despite this fact, we present here the complete projections.

Results

Models TSS’ values ranged in natural distribution from 0.186 to 0.744 with mean value of 0.562 and standard deviation of 0.112. When considering complete data, the values ranged from 0.151 to 0.787 with mean value of 0.569 and standard deviation of 0.116. Ecological niche model unveiled that currently P. echinata has a potential suitable distribution ranging from 1,137,920 to 1,794,434 km² (Tab. 2). Of this total area, 7.10 % to 7.39 % are considered to be within protected areas (Tab. 3). Future scenarios showed a gradual reduction in species suitable area (with increasing CO₂ concentration), ranging from 1,648,148 km² for RCP 2.6 in the year 2070, considering both natural and cultivated records, to 376,479 km² for RCP 8.5 in the year 2070, considering just records in the natural occurrence area (Tab. 2). In the same way, potential distribution inside protected areas in future scenarios ranges from 144,155 km² for RCP 4.5 in the year 2070, considering both natural and cultivated records, to 45,645 km² for RCP 8.5 in the year 2070, considering records in the natural occurrence area (Tab. 3).

Table 2
Total summed area in current and future scenarios using a weighted method, which results in conservative values.

Scenario	Natural (km2)	Cultivated (km2)	Difference (km2)
Current	1,137,920	1,794,434	656,513
RCP-2.6/2050	850,905	1,648,148	797,243
RCP-2.6/2070	844,532	1,691,666	847,135
RCP-4.5/2050	718,347	1,532,665	814,318
RCP-4.5/2070	608,911	1,528,558	919,647
RCP-6.0/2050	780,776	1,586,449	805,673
RCP-6.0/2070	609,799	1,487,229	877,430
RCP-8.5/2050	605,881	1,533,199	927,318
RCP-8.5/2070	376,479	1,409,438	1,032,959

Table 3
Protected area in current and future scenarios considering current Brazilian protected areas network. Percentage values are in terms of the area presented in Table 2.

Scenario	Natural (km2)	Natural (%)	Cultivated (km2)	Cultivated (%)	Difference (km2)
Current	84,048	7.39	127,347	7.10	43,299
RCP-2.6/2050	84,293	9.91	137,040	8.31	52,747
RCP-2.6/2070	87,183	10.32	141,573	8.37	54,390
RCP-4.5/2050	81,693	11.37	143,090	9.34	61,397
RCP-4.5/2070	71,646	11.77	144,155	9.43	72,509
RCP-6.0/2050	88,431	11.33	143,465	9.04	55,034
RCP-6.0/2070	70,951	11.64	138,180	9.29	67,229
RCP-8.5/2050	68,561	11.32	141,270	9.21	72,708
RCP-8.5/2070	45,645	12.12	132,435	9.40	86,791

Projections from natural occurrence models revealed a considerably larger distribution than the previously known for P. echinata (Fig. 2). The states of Mato Grosso, Goiás, Minas Gerais and the continental portion of Bahia seem to have the same environmental suitability as its natural distribution. Furthermore, the state of São Paulo, as well as scattered sites in Mato Grosso do Sul and even Paraná may also be considered good places for that species to grow.

Figure 2
Committee average of projections for current scenario.

Future projections indicate a shift trend to Southeast Brazil and loss of several areas throughout its distribution range. It is expected to remain, in the worst cases, in few places along the Brazilian coast (Espírito Santo and Alagoas) and in the state of Minas Gerais and São Paulo (Fig. 3). Considering cultivated records, projections showed a similar trend along its distribution, but with higher overall suitability, remaining, in the worst case, with suitable areas in the states of São Paulo and Rio de Janeiro in addition to projections considering just natural occurrences (Fig. 4). Most of all, the coastal region of Pernambuco and Alagoas, as well as the South of Espírito Santo and São Paulo are key regions to the preservation of P. echinata, since they will have high suitability even in the worst scenarios. We identified three main climatic stable areas: the Northeast Brazil (Fig. 5A, a), the highlands of Southeast Brazil (mountains of Espírito Santo and Minas Gerais, Fig. 5A, b) and highlands of São Paulo (Fig. 5A, c). Stable areas presented a low number of future scenarios agreeing with the presence of P. echinata inside most protected areas.

Figure 3
Future projections (committee average) of Paubrasilia echinata considering records in its natural occurrence only.

Figure 4
Future projections (committee average) of Paubrasilia echinata considering records in its natural occurrence, as well as cultivated records.

Figure 5
Climatic stable areas to Paubrasilia echinata along future scenarios considering natural presence records. Red cells indicate high-stability sites. Blue cells indicate low-stability sites. In zoomed panels, green polygons indicate protected areas’ boundaries. Lowercase letters identify the three main climatic stable areas: the Northeast Brazil (a), the highlands of Southeast Brazil (mountains of Espírito Santo and Minas Gerais, b) and highlands of São Paulo (c).

Discussion

Climatically suitable area for P. echinata is larger than previously considered (Gagnon et al. 2016), although it may not fully colonize area within the continent since the coast of Brazil has a major altitudinal shift. In this way, sierras of the Brazilian coast are possibly acting as a barrier to migration, even though there is a suitable area for the species to grow. The Mean Diurnal Temperature Range (BIO2) variable was firstly selected to generate models, since it is a proxy of continentality, but it had high correlation with other variables, thus it was discarded after our selection routine. P. echinata thrives in a relatively wide variety of habitats, from coastal outcrops and white-sand forests to seasonal and rain forests, but always in well-drained soils (Gagnon et al. 2016). To consider soil attributes in ecological niche models may be particularly difficult, especially when projecting models into future scenarios, since soil will not change as quickly as the climate, resulting in a higher importance to this variable. The lack of studies on the species distribution limits our inferences over which environmental filters are working to prevent P. echinata’s expansion into areas where it could potentially occur. Its relation with well-drained soils may give us a clue on why it does not colonize inner continent areas and other biomes such as Amazônia and Pantanal, despite the evidenced climatic suitability. Regardless of that, projections show an altitudinal shift in suitability through future scenarios.

As we didn’t consider landscape metrics in this study, such as structural and functional connectivity, edge effect or patch size, our inferences over populations dynamics are weak, and it is essential in further studies. The increase of relatively protected area within future scenarios, followed by a decrease in absolute protected and total area, indicates that habitat loss will be more intense outside conservation units. This leads us to a problematic scenario, since a major driver of species extinction is habitat fragmentation, which is particularly severe in P. echinata populations due to exploitation and selective logging (Cardoso et al. 2005). Despite that, it makes necessary to point out that, due to this long-term relative isolation of the natural populations, there is a high genetic diversity, with markedly morphological differences that may be recognized in future studies as new taxa (Lima 2019), and which may grant species quick adaptation, resilience and persistence through scenarios, although genetic diversity within populations is very low due to the discontinuity of the distribution.

Despite the fact that none of the protected areas that houses Pau-brasil nowadays may sustain it in future scenarios, the species will not be unprovided of protection. Parque Nacional da Serra de Itabaiana (Sergipe), Refúgio da Vida Silvestre Mata do Urucu (Pernambuco), Reserva Biológica da Pedra Talhada (Alagoas), Reserva Biológica Augusto Ruschi (Espírito Santo), Área de Proteção Ambiental Bacia do Paraíba do Sul (São Paulo) and Parque Estadual da Serra do Mar (São Paulo), as well as three different private reserves (Reservas Particulares do Patrimônio Natural - RPPNs) are predicted to harbor P. echinata in the next 50 years.

Climate stability maps (Fig. 5) can also be interpreted as potential areas for reforestation and restoration, which is particularly interesting for P. echinata, since light intensity in open areas of the forest positively influences growth rates (Zani et al. 2012). This species can act as a colonizer and thus as a promoter of the restoration process in degraded sections on forest edges (Sattler et al. 2018), facilitating as well the species migration.

Paubrasilia echinata will have a major obstacle in future scenarios, which is to adapt to new environments where it is not growing naturally. It may be capable to do that, since cultivated organisms show a high adaptation capacity to a myriad of habitats, developing relatively well. The significant difference in area between models (natural vs. cultivated) unveiled a huge potential for species adaptation. To consider cultivated records can drive future researches with ecological niche models to other cultivated species as well. Although some cultivated organisms may not bloom or fruit (e.g. due to ecological needs) or generate viable seeds (e.g. due to pollinators absence), these records amplify the niche breath of the species.

Even if P. echinata suffers greatly with climate changes, it is not likely to be extinct since it is an important symbol of our country, having not just a crucial role in Brazil’s history, but also in the development of Brazilian culture. P. echinata may be a potential flagship species for the conservation of Brazilian flora, since it has public appeal and is not restricted to one habitat. This research may also be key to a future national action plan for this species, which is a must to properly protect this long-exploited tree.

Acknowledgements

This study was financed in part by the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - Brasil (CAPES) - Finance Code 001. LFE was supported by CAPES PhD scholarship and DDS was supported by a postdoctoral fellowship (PNPD-CAPES). We thank three helpful anonymous reviewers for their contribution to improve this manuscript.

References

• Araújo MB, Alagador D, Cabeza M, Nogués-Bravo D, Thuiller W. 2011. Climate change threatens European conservation areas. Ecology Letters 14: 484-492.
• Bellard C, Leclerc C, Leroy B, et al 2014. Vulnerability of biodiversity hotspots to global change. Global Ecology and Biogeography 23: 1376-1386.
• Berteaux D, Ricard M, St-Laurent MH, et al 2018. Northern protected areas will become important refuges for biodiversity tracking suitable climates. Scientific Reports 8: 4623. doi: 10.1038/s41598-018-23050-w
  » https://doi.org/10.1038/s41598-018-23050-w
• Brasil. 1605. Regimento do Pau-Brasil. https://www.historiadobrasil.net/documentos/pau_brasil.htm
  » https://www.historiadobrasil.net/documentos/pau_brasil.htm
• Cardoso SRS, Provan J, Lira CDF, Pereira LDOR, Ferreira PCG, Cardoso MA. 2005. High levels of genetic structuring as a result of population fragmentation in the tropical tree species Caesalpinia echinata Lam. Biodiversity and Conservation 14: 1047-1057.
• Castro CFA. 2002. Gestão florestal no Brasil colônia. PhD Thesis, Universidade de Brasília, Brasília.
• CNCFlora. 2012. Caesalpinia echinata In: Lista Vermelha da flora brasileira versão 2012.2. Rio de Janeiro, Centro Nacional de Conservação da Flora.
• Diniz-Filho JAF, Bini LM, Rangel TF, et al 2009. Partitioning and mapping uncertainties in ensembles of forecasts of species turnover under climate change. Ecography 32: 897-906.
• Elith J, Graham H, Anderson R, et al 2006. Novel methods improve prediction of species’ distributions from occurrence data. Ecography 29: 129-151.
• Esser LF, Neves DM, Jarenkow JA. 2019. Habitat‐specific impacts of climate change in the Mata Atlântica biodiversity hotspot. Diversity and Distributions 2019;00: 1-11. doi: 10.1111/ddi.12984
  » https://doi.org/10.1111/ddi.12984
• Fontana R, Cunha EA, Mendes JP. 1994. O Brasil de Americo Vespúcio. Brasilia, Editora Universidade de Brasilia.
• Franklin J, Serra-Diaz JM, Syphard AD & Regan HM. 2016. Global change and terrestrial plant community dynamics. Proceedings of the National Academy of Sciences of the United States of America 113: 3725-3734.
• Gagnon E, Bruneau A, Hughes CE, Paganucci L, Lewis GP. 2016. A new generic system for the pantropical Caesalpinia group (Leguminosae). PhytoKeys 71: 1-160.
• Guisan A, Thuiller W. 2005. Predicting species distribution: Offering more than simple habitat models. Ecology Letters 8: 993-1009.
• Hijmans RJ, Cameron SE, Parra JL. 2006. Worldclim global climate layers. Version 1.4. WorldClim database. http://www.worldclim.org
  » http://www.worldclim.org
• IPCC. 2013. IPCC Fifth Assessment Report (AR5). IPCC, s. 10-12.
• IUCN 2019. The IUCN Red List of Threatened Species. Version 2019-1. http://www.iucnredlist.org
  » http://www.iucnredlist.org
• Li X, Clinton N, Si Y, Liao J, Liang L, Gong P. 2015. Projected impacts of climate change on protected birds and nature reserves in China. Science Bulletin 60: 1644-1653.
• Lima HC. 2019. Paubrasilia In: Flora do Brasil 2020 em construção. Rio de Janeiro, Jardim Botânico do Rio de Janeiro.
• Lira CF, Cardoso SRS, Ferreira PCG, Cardoso MA, Provan J. 2003. Long-term population isolation in the endangered tropical tree species Caesalpinia echinata Lam. revealed by chloroplast microsatellites. Molecular Ecology 12: 3219-3225.
• Longui EL, Yojo T, Lombardi DR, Alves ES. 2010. The potential of ipê (Handroanthus spp.) and Macarandura (Manilkara spp.) woods in the manufacture of bows for string instruments. IAWA Jounal 31:149-160.
• Loyola RD, Lemes P, Brum FT, Provete DB, Duarte LD. 2014. Clade-specific consequences of climate change to amphibians in Atlantic Forest protected areas. Ecography 37: 65-72.
• Moss RH, Edmonds JA, Hibbard KA. 2010. The next generation of scenarios for climate change research and assessment. Nature 463: 747-756.
• Naimi B, Hamm Na, Groen TA, Skidmore AK, Toxopeus AG. 2014. Where is positional uncertainty a problem for species distribution modelling. Ecography 37: 191-203.
• Parmesan C. 2006. Ecological and evolutionary responses to recent climate change. Annual Review of Ecology, Evolution, and Systematics 37: 637-669.
• Parmesan C, Yohe G. 2003. A globally coherent fingerprint of climate change impacts across natural systems. Nature 421: 37-42.
• Peterson AT. 2001. Predicting species’ geographic distributions based on ecological niche modeling. The Condor 103: 599-605.
• Peterson AT, Egbert SL, Sánchez-Cordero V, Price KP. 2000. Geographic analysis of conservation priority: Endemic birds and mammals in Veracruz, Mexico. Biological Conservation 93: 85-94.
• Peterson AT, Ortega-Huerta MA, Bartley J, et al 2002. Future projections for Mexican faunas under global climate change scenarios. Nature 416: 626-629.
• R Development Core Team. 2018. R: A language and environment for statistical computing. Vienna, R Foundation for Statistical Computing. https://www.R-project.org/
  » https://www.R-project.org/
• Root TL, Price JT, Hall KR, Schneider SH, Rosenzweig C, Pounds JA. 2003. Fingerprints of global warming on wild animals and plants. Nature 421: 57-60.
• Sales LP, Neves OV, Marco Jr. P, Loyola R. 2017. Model uncertainties do not affect observed patterns of species richness in the Amazon. PLOS ONE 12: e0183785. doi: 10.1371/journal.pone.0183785
  » https://doi.org/10.1371/journal.pone.0183785
• Sattler D, Raedig C, Hebner A, Wesenberg J. 2018. Use of native plant species for ecological restoration and rehabilitation measures in Southeast Brazil. In: Nehren U, Schlüter S, Raedig C, Sattler D, Hissa H. (eds.) Strategies and tools for a sustainable rural Rio de Janeiro. Rio de Janeiro, Springer International Publishing.
• Skeaping K. 1955. Some speculations on a crisis in the history of the violin. Galpin Society Journal 8: 3-12.
• Souza BJ. 1939. O Pau-brasil na História Nacional. São Paulo, Companhia Editora Nacional e MEC.
• Thomas CD, Gillingham PK. 2015. The performance of protected areas for biodiversity under climate change. Biological Journal of the Linnean Society 115: 718-730.
• Thuiller W. 2004. Patterns and uncertainties of species’ range shifts under climate change. Global Change Biology 10: 2020-2027.
• Thuiller W, Lavergne S, Roquet C, Boulangeat I, Lafourcade B, Araujo MB. 2011. Consequences of climate change on the tree of life in Europe. Nature 470: 531-534.
• Zani LB, Sarnaglia Junior VB, Gomes JML, Thomaz LD. 2012. Estrutura de um fragmento de Floresta Atlântica em regeneração com ocorrência de Caesalpinia echinata Lam. (pau-brasil). Biotemas 25: 75-89

Publication Dates

• Publication in this collection
  25 Nov 2019
• Date of issue
  Oct-Dec 2019

History

• Received
  15 May 2019
• Accepted
  26 Sept 2019