*Article* **Predicting Suitable Habitats of the African Cherry (***Prunus africana***) under Climate Change in Tanzania**

#### **Richard A. Giliba 1,\* and Genesis Tambang Yengoh <sup>2</sup>**


Received: 5 August 2020; Accepted: 25 August 2020; Published: 15 September 2020

**Abstract:** *Prunus africana* is a fast-growing, evergreen canopy tree with several medicinal, household, and agroforestry uses, as well as ecological value for over 22 countries in sub-Saharan Africa. This species is under immense pressure from human activity, compounding its vulnerability to the effects of climate change. Predicting suitable habitats for *P. africana* under changing climate is essential for conservation monitoring and planning. This study intends to predict the impact of climate change on the suitable habitats for the vulnerable *P. africana* in Tanzania. We used maximum entropy modeling to predict future habitat distribution based on the representative concentration pathways scenario 4.5 and 8.5 for the mid-century 2050 and late-century 2070. Species occurrence records and environmental variables were used as a dependent variable and predictor variables respectively. The model performance was excellent with the area under curve (AUC) and true skill statistics (TSS) values of 0.96 and 0.85 respectively. The mean annual temperature (51.7%) and terrain ruggedness. index (31.6%) are the most important variables in predicting the current and future habitat distribution for *P. africana*. Our results show a decrease in suitable habitats for *P. africana* under all future representative concentration pathways scenario when compared with current distributions. These results have policy implications for over 22 countries of sub-Saharan Africa that are facing problems associated with the sustainability of this species. Institutional, policy, and conservation management approaches are proposed to support sustainable practices in favor of *P. africana*.

**Keywords:** habitat suitability; species distribution; climate change; conservation; *P. africana*

#### **1. Introduction**

Climate change is a serious threat to floral biodiversity conservation [1–3]. Over the next century, the global temperature is projected to rise by 0.3–4.8 ◦C [4]. The projected increase in temperature is anticipated to influence both species and habitat distribution in several different ways. Some of the climate change impacts on flora include changes in species distribution, the increased extinction rate of species, changes in the length of growing season, and reproduction timings for plants [5]. With a rise in temperature, plant species are likely to shift their distribution patterns depending on resource availability [6]. This may lead to species range expansion or range contraction, or range shift when respond to changing climate [7]. Predicting suitable habitat of species under climate change using species distribution models (SDM) is one of the key important steps to undertake conservation planning and management [6,8]. Greenhouse gases (GHGs) from various human activities are the primary agents answerable for climate change and the current emission rates are at the highest level in the recorded history [9].

SDM have been widely used to monitor the impacts of climate change on the floral distribution [10] and identification of suitable habitats of species [11]. SDM uses species occurrence and environmental variables [12] to predict the distribution of a species in a geographical or an environmental space [13]. SDM engage a variety of methods for predicting species distribution and mapping habitat suitability [14]. These include: MAXENT, BIOCLIM, DOMAIN artificial neural networks, generalized linear models, and generalized additive models. Maxent, Bioclim, and Domain use presence-only data [15] while others require both presence absent data [14]. Presence-only methods include bioclimatic envelope algorithm BIOCLIM, DOMAIN.

Maxent has been widely used and gives better results when compared to other different modeling methods that use presence-only data are used [14,16,17]. Several studies have demonstrated Maxent's ability to accurately predict species distribution in a wide range of ecological and geographical regions [10,18,19]. Subsequently, conservation practitioners have been increasingly using habitat suitability models from Maxent to make conservation management decisions [20,21]. We used Maxent: first to identify the most important variables that govern the current distribution of *P. africana*; second to predict current suitable habitats; third to predict future suitable habitats based on two representative concentration pathways (RCP 4.5 and RCP 8.5) for the mid-century 2050 and late-century 2070 in Tanzania. We selected *P. africana* due to local and international economic importance for medicinal purposes, contributing to its overexploitation. Due to the changing climate and overexploitation and, there is a need to use species occurrence data for identification of key conservation sites so as to develop a countrywide conservation strategy.

#### *1.1. P. africana–Its Value, Demand and Conservation Pressures*

*P. africana* is a fast-growing, evergreen canopy tree about 30–40 m in height. It has a wide range, spanning several countries in central, western, southern, and eastern Africa (including Madagascar). It occupies habitats in upland rain-forest, montane and riverine forests, moist evergreen forest, and the edges of dry gallery forests [22]. The uses of *P. africana* are many, and explain the huge demand for this tree and its products throughout the world:

Medicinal uses: The tree is widely used in both traditional and modern medicine to treat a variety of ailments [23–25]. *P. africana* contains several medically active compounds including the cyanogenic glycoside amygdalin, which is found in the bark, leaf, and fruit; phytosterols such as. β-sitosterol 15–18%, and its 3-O-glycoside, β-sitostenone, campesterol, and aucosterol; pentacyclic triterpenoids [26,27]. The bark is highly valued for its medicinal properties, particularly as a treatment for benign prostatic hyperplasia and prostate gland hypertrophy, diseases that commonly affect older men in Europe and North America [28]. In traditional medicinal practices throughout sub-Saharan Africa, the bark is used in traditional medicine as a purgative and as a remedy for stomachache, while the leaves are used as an inhalant for fever or are drunk as an infusion to improve appetite. Demand for *P. africana* bark for medicinal uses has been high and growing, putting immense pressure on the species throughout sub-Saharan Africa. In Tanzania, this has raised issues of sustainability for the species [29]. While data on current demand and supply of *P. africana* barks in the international and local markets is scarce, it was expected in 2000 that this demand would triple or quadruple to 7000 to 11,000 tons/year in export and about 500 tons/year for use in Africa the years after the report [30].

Household uses: Domesticated trees serve as shade in compounds, as windbreaks, and as ornamental trees. The tree yields a high-quality fuel, and so is a favorite for the production of charcoal or for use as fuel wood in many communities. Regarding household socioeconomics, *P. africana* supports revenue streams of communities in a wide variety of ways. As an input to furniture production, the seasoned wood saws easily and cleanly; works well with hand and machine tools; and polishes and finishes well. It also serves as highly desirable timber for flooring and heavy construction where durability is not required.

Agroforestry uses: Where the tree has been domesticated for integration into agricultural systems, the tree is used for erosion control, while the leaves are incorporated into the system of organic manure

complex of crop production. There have been efforts to more tightly incorporate *P. africana* into the agroforestry mix of agricultural systems in some parts of the continent. For example, in Cameroon where the harvesting and commercialization have come under more targeted scrutiny given the large volumes involved, the government established a "National Management Plan for *P. africana*" that supported greater efforts towards reducing the pressures on wild tree species with locally cultivated trees [31].

Ecological value: One of the most cited ecological value of *P. africana* is its preference by the black and white colobus monkeys (*Colobus guereza*) as their top food species [32]. In forests, agricultural fields, as in homes, *P. africana* also serves as shelter for a variety of bird species. It is a valuable species for beekeeping, and therefore an important contributor in supporting pollination services that are relied on by forest species, and for the success of agriculture.

The heavy pressures to which *P. africana* has recently come under in most African countries because of wild harvesting for the medicinal plant trade have not gone unnoticed in the international and biodiversity community. In 1995, *P. africana* was added to Appendix A of the Convention on International Trade in the Endangered Species of Wild Fauna and Flora's (CITES) list of endangered species, for the regulation of its trade from wild harvesting [33]. Currently, all exports of *P. africana* should therefore subject to a CITES export permit to protect the tree from depletion in Africa. In response, a European Union (EU) ban on imports of *P. africana* bark came into force in 2007 to help stocks to recover [34]. The International Union for Conservation of Nature (IUCN) recognize *P. africana* as a vulnerable species on its red list (https://www.iucnredlist.org) and it was categorized as a species "of urgent concerns" by CITES [22].

#### *1.2. Forests and the Circular Economy*

The circular economy is defined as "the concept can, in principle, be applied to all kinds of natural resources, including biotic and abiotic materials, water and land. Eco-design, repair, reuse, refurbishment, remanufacture, product sharing, waste prevention and waste recycling are all important in a circular economy" [35]. The concept of the circular economy seeks to achieve a shift from the linear economy, which is characterized by less than optimal recycling and reuse of materials and resources in human societies. The overall goal of the circular economic model is to reduce the undesirable impacts of the linear economy by achieving a systemic transition into a more sustainable approach to natural resource exploitation and use built on long-term sustainability. One of the main objectives of the circular economy is to reduce the impact of human activities on the planet's ecosystems by reducing the excessive exploitation of natural resources and minimizing the pressure of human actions on the functioning of these ecosystems.

The role of forests in human societies and development, as well as the nature of forestry and forestry-based industries makes them a prime candidate for contributing to the global drive towards achieving the goals of a circular economy [36,37]. By striving towards achieving the objectives of the circular economy, positive contributions can be made towards achieving several sustainable development goals [38], many of which are relevant for Tanzania's development.

#### **2. Material and Methods**

#### *2.1. Study Area*

The United Republic of Tanzania is located in East Africa between longitude 29◦ and 41◦ East and latitude 1◦ and 12◦ South (Figure 1). Tanzania is endowed with a wide range of natural resources as well as ecological and cultural diversity including extensive areas of arable land, wildlife reserves and parks, mountains, forest reserves, rivers, and lakes. The mean annual rainfall varies from below 500 mm to over 2000 mm per annum while the mean temperature ranges from −4.9 ◦C to 27.9 ◦C per annum. The central and western plateau is relatively dry while, the northern and southern highland are cool. Rainfall for large parts of the country is bimodal with short rains from October-December

and long rains from March to May [39]. The country has 7 hotspots including forest reserves, nature reserves game reserves, and national parks that are recognized by The United Nations Educational, Scientific and Cultural Organization as World Heritage sites (https://whc.unesco.org). A significant number of world endemic and threatened species are reported from Tanzania (https://www.cbd.int). However, the country has lost its forest cover from land use change and it is threatened by changing climate [40].

**Figure 1.** Location map of the study area Tanzania. Black points show species occurrence. Note that the data for *P. africana* is superimposed on top of elevation layer.

#### *2.2. Species Presence Records*

We obtained the present locations of *P. africana* from a 5-year field survey done across the country by National Forest Resources Monitoring and Assessment Project, and different online sources, TROPICOS an online botanical database containing taxonomic information on plants (http://www.tropicos.org) and Global Biodiversity Information Facility database (http://www.gbif.org). A total of 187 records were collected, and after screening, 57 duplicate records were removed, and finally 120 records were used to run the model (Figure 1). To model potential attribution of *P. africana* across the country.

#### *2.3. Environmental Variables*

We collected 19 bioclimatic variables from WorldClim dataset (https://www.worldclim.org). To derive elevation and terrain ruggedness index, we downloaded a digital elevation model from Shuttle Radar Topography Mission dataset (http://srtm.csi.cgiar.org). Soil type (Table 1) were obtained from the International Soil Reference and Information Centre database (https://www.isric.org). We resampled both soils and topographic layers to the resolution of bioclimatic variables (~1 km) using ArcGIS 10.6. To reduce multi-collinearity of climate variables, the two variables that found to have a high correlation coefficient (|r| > 0.7), as suggested by [41], we selected "one variable for modeling due to its ecological importance for the survival of *P. africana*" [42]. This resulted in the inclusion of eight variables for modeling. Table 2 lists the general statistics of the major environmental variables used in this study. We used Climate Community Climate System version four (CCSM4) bioclimatic variables to predict

the future distribution of *P. africana* under future climate scenarios, namely representative concetration pathway (RCP) 4.5 a moderate greenhouse gas emission scenario and RCP 8.5 a extreme greenhouse gas emission scenario for mid-century and late-century. The CCSM4 is among the most commonly used bioclimatic variables to predict the impact of climate change on plant distribution [1,3].




**Table 2.** Environmental variables used to model distribution of *P. africana.*

#### *2.4. Species Distribution Modeling*

We used Maxent version 3.3.3; [17] to model the distribution of *P. africana* in this study due to the unavailability of absence records. Maxent uses presence records in combination with environmental conditions the species is present to model the spatial distribution based on the theory of maximum entropy [14]. During modeling, we selected 75% of presence records to training the model and 25% for testing the model [1,3,15], while changing Maxent setting. We tried to set various values for the regularization multiplier and the number of iterations and changed feature types. We obtained the good results with the following settings; cross-validate with iterations set to 5000, regularization multiplier set to 1, and feature type set to quadratic, hinge, and linear. Further, the maximum number of background points was set to 10,000, and replicates were set to 30. Afterward, we imported current and future predicted maps for *P. africana* from Maxent models into ArcGIS 10.6 and reclassified into five classes of potential habitats according to [43]: unsuitable habitat (0–0.2); barely suitable habitat (0.2–0.4); suitable habitat (0.4–0.6); highly suitable habitat (0.6–0.8); very highly suitable habitat (0.8–1). Finally, current distribution maps were subtracted from future maps to compute the relative changes in species range (decreasing or increasing) [2].

#### *2.5. Model Evaluation and Validation*

We used the area under receiver operating characteristic (AUC) and true skill statistic (TSS) to assess the performance of model. The AUC values range between 0–1; higher AUC values suggest the better and higher performance of a model [14,17]. "TSS values range between +1 to −1; a values > 0.8 suggest excellent, 0.4–0.8 useful, and <0.4 poor model performance" [2]. Finally, we selected the model with highest AUC and TSS. Besides, we used jackknife test to identify important variable governing the distribution of *P. africana*. Further, we use response curves to show how the predicted probability of presence changes as each environmental variable is varied.

#### **3. Results**

#### *3.1. Model Validation and Influencing Bioclimatic Variables*

Model for *P. africana* provided satisfactory results, with AUC and TSS values of 0.957 and 0.845 respectively. These suggest that the model for *P. africana* produced good results. Annual mean temperature (bio1) contributed most to the model, followed by terrain ruggedness index (tri) (Table 3). The cumulative contribution of these two variables is 83.30%. On the other hand, the variable with the highest gain when used in isolation is bio1, this implies this variable has the most useful information by itself. The variable that decreases the gain most when it is omitted is tri, this implies that the tri has the most information that is not present in the other variables (Figure 2). These results signify that bio1 and tri a proxy measure of topographic heterogeneity are the master variables governing the current and future distribution of *P. africana* in Tanzania.


**Table 3.** Environmental variables used in the study and their percentage contributions, and the maps that show the spatial distribution of the important variables are presented in Figure A1.

**Figure 2.** Jackknife test results indicating variable with highest gain when used in isolation and variable that decreases the gain the most when it is omitted. The test results indicate annual mean temperature (bio1) has the most useful information by itself while terrain ruggedness index (tri) has the most information that is not present in the other variables. Jackknife of regularized training gain for *P. africana*.

Figure 3, below, shows how the predictions depend on the variables, as mean annual temperature (bio1), increases habitat suitability for *P. africana* decreases while its habitat suitability increases with annual precipitation (bio12). On the other hand, as elevation (eleva) and terrain ruggedness index (tri), increases habitat suitability for *P. africana* increases indicating that it prefers undulating upland areas. Further, *P. africana* prefers to reside on nitisols, histosols, leptosols, and acrisols soils, which is widely distributed in undulating upland areas.

#### *3.2. Current and Future Distribution of P. africana*

Predicted distributions under current conditions revealed that highly and moderately suitable areas for *P. africana* covers only 6.69% (62,388.75 sq. km), while low and very low suitable areas cover the large portion, 93.31% (869,937.04 sq. km), of the study area (Table 4; Figure 4). The highly suitable areas to a large extent are identified in the southern highlands, western, and northern zones of the study area (Tanzania, Figure 4). Predicted distributions under future conditions indicates decline in suitable areas and increase in suitable areas under RCP 4.5 and RCP 8.5 scenarios for mid-century 2050 and late-century 2070 (Table 4; Figure 4). Climatically highly and moderately suitable areas will decline by 2.29% and 3.07% under RCP 8.5 for 2050 and 2070, respectively while under RCP 4.5 for 2050 and 2070 highly and moderately suitable areas will decline by 2.10% and 2.20% respectively. Climatically very low areas will increase by 6.85% and 8.59% under RCP 8.5 for 2050 and 2070, respectively while under RCP 4.5 for 2050 and 2070 suitable conditions will decline by 6.25% and 6.59% respectively (Table 4; Figure 4). Southern highlands, western, and northern zones are anticipated to lose large portion of suitable areas in the future for all climate scenarios under mid-century and late-century (Table 5; Figure 4).

**Figure 3.** Relationships between selected environmental variables and probability of species suitability of *P. africana* (**A**) mean annual temperature (bio1) (**B**) annual precipitation (bio12), (**C**) elevation (eleva), (**D**) terrain ruggedness index (tri), (**E**) precipitation of warmest quarter (bio18), (**F**) soil types (soils). The y-axis represents the probability of presence (cloglog output). Red curves show the average response and blue margins are ± SD calculated by 30 replicate runs. For the interpretation of soil type legend, refer Table 1.


**Table 4.** Change in suitable areas of *P. africana* country-wide for mid-century and late-century under representative concertation pathway (RCP) 4.5 and RCP 8.5 scenarios.

**Figure 4.** Climate Community Climate System version four (CCSM4) climate model based predicted future suitability of *P. africana* species: (**A**) current potential distribution; (**B**,**C**) RCP 4.5, and 8.5 emission scenario for 2050; (**E**,**D**) RCP 4.5, and 8.5 scenario for 2070.


**Table 5.**Change in suitable areas of*P. africana*zone wide for mid-century and late-century under RCP 4.5 and RCP 8.5 scenarios.

#### **4. Discussion**

The findings of the study show that the suitability distribution of the *P. africana* is largely controlled by annual temperature, terrain ruggedness index, elevation, and soil type (contributed more than 90%). *P. africana* showed sharp decline in response to an increase in annual temperature beyond signifying that the probability of occurrence of the species may be affected with higher temperature. This is in line with the results from a study in eastern arc mountains Tanzania, where the studied tree species showed a declining trend in response to an increase in temperature [44]. With warming trends, plant species are expected to track the changing climate and shift their distributions to the extent that resource availability allows [6]. *P. africana* showed an increasing trend in response to an increase of the terrain ruggedness index suggesting that the species prefers rugged or undulating areas [45]. "Terrain ruggedness index as a measure of terrain heterogeneity is an important variable for predicting which habitats are used by a species and the density at which species occur in a variety of environments" [46]. *P. africana* showed an increasing trend in response to the increase of elevation indicating that the species prefers high elevation areas. Higher temperature is stated to cause shifts in plant distribution along the elevation gradients [47]. *P. africana* distribution appears to associate with Nitisols, Histosols, and Leptosols soils that are found in undulating upland areas. Soil type plays a major role in the heterogeneity of habitats, thus determining the distribution of plant species [48].

#### *4.1. Management Implications*

Our results indicate that climate change will pose a severe impact on the future distribution of *P. africana* in Tanzania. Research institutions and public universities can take an interest in both in situ and ex situ long-term monitoring trends of *P. africana* distribution in a country. In-situ interventions should focus on the "recruitment and regeneration of the species while ex-situ interventions should target to promote tree retention on farms, or advocate further planting, collect specimens, and to establish gene banks and botanical gardens to ensure the survival of the species" [42].

#### *4.2. Institutional and Policy Context for Addressing Challenges Associated with P. africana*

Addressing challenges associated with the vulnerability of *P. africana* in Tanzania can benefit from the institutional context already in place in the country. The government has put in place institutional frameworks to manage natural resources and environment-related initiatives and challenges countrywide. The President's Office-Regional Administration and Local Government (PORALG) works closely with Local Government Authorities through their various departments in collaboration with the respective sector ministries to implement the strategic interventions at the local level (municipalities, districts, wards, villages, and sub-villages). This is important for addressing the immediate management decisions that directly affect the health and survival of *P. africana*, as the trees are directly impacted by community demand for livelihoods at the local level. Successful implementation of policies, laws, and plans also requires enhanced engagement with Civil Society Organizations, development partners, the private sector, and academic and research institutions.

Addressing the more systemic and long-term environmental challenges facing *P. africana* would be best addressed within the institutional arrangements for environmental management functions in Tanzania. These are two basic types of such functions: (i) Sectoral Environmental Management Functions (also known as Type A functions) that are concerned with the management of specific natural resources or environmental services, such as forestry, agriculture, fisheries, wildlife, mining, and waste management. These functions are to a large extent operational and guided by sector-specific policies and acts such as the Forest Act (2002), which should be directly relevant for addressing challenges of *P. africana*. (ii) Coordinating and Supporting Environmental Management Functions (commonly referred to as Type B functions) involve the task of providing central support functions by coordinating and supporting the different and sometimes conflicting Type A activities and integrating them into an overall sustainable system. Specific tools within this Type B functions relevant for addressing challenges

of *P. africana* are the National Environmental Policy (1997) and the Environmental Management Act (2004), which provide policy and legislative framework for the coordination of the implementation of policies and laws on environmental and natural resources management.

#### *4.3. Conservation and Management Approaches to Support Sustainable Practices in Favor of P. africana*

Three main conservation and management approaches can offer possibilities to address the vulnerability of *P. africana* in Tanzania. These draw from experiences in other parts of the developing world facing conservation challenges of their own and the lessons learned from their conservation and management approaches.

#### 4.3.1. Supporting Inclusive Conservation Approaches

While climate changing is increasingly representing a challenge to the distribution and health of *P. africana*, this challenge comes to compound existing pressures imposed by anthropogenic pressures of the species. Given the high degree of anthropic intervention contributing to the vulnerability of the species, conservation and management strategies need to be closely linked to the needs of indigenous peoples and local communities. There is therefore a serious need to consider the genuine and effective participation of indigenous peoples and local communities in the definition and application of resource management options when addressing challenges to the sustainable management of *P. africana*. Given the importance of this specie to the health of local populations, socio-economic and environmental welfare, emphasis on supporting its sustainable and adaptive use, supported by initiatives that reduce local reliance on the harvesting of wild resources (such as agroforestry) may prove to be more locally acceptable. This principle of "conservation through use" is an example of the application of community-based natural resources management models to address issues of resource degradation. These types of co-management models have been applied to the successful management of protected resources in other parts of the world [41]. Such an approach will be in alignment with current governmental efforts towards a more inclusive management of forests and natural resources. A program of Participatory Forest Management has been introduced and operationalized through the Joint Forest Management (JFM) as well as a Community Based Forest Management (CBFM) processes across the country. Under JFM, agreements between community groups and the Government have been developed with a view to promoting the participation of communities in the management and utilization of forest resources. The Community Based Forest Management program encourages communities to set up forest reserves from the general lands for economic and conservation activities.

#### 4.3.2. Collaboration to Streamline and Align Regional and International Efforts

Given the regional and international character of challenges of *P. africana* vulnerability, there is a need for collaboration to improving forest management, share best practices, and support effective conservation as well as the production and trade of forest products. The Collaborative Partnership on Forests (CPF) is an informal, voluntary arrangement among 14 international organizations to share experiences and build on them to produce new benefits for their respective national stakeholders in the forest resources sector [41]. Addressing the challenges of *P. africana* has the potential of benefiting from the CPF whose mission is to promote sustainable management of all types of forests and to strengthen long-term political commitment to this end.

#### 4.3.3. Leveraging the Potential of Payments for Ecosystem Services (PES)

Payments for ecosystem services (PES) refer to voluntary transactions between users and suppliers of environmental services, such that suppliers are subject to natural resource management and handling rules within and outside of service provision areas [49]. Under a PES scheme, users of land upstream may agree to voluntary limitation or diversification of their activities in return for an economic benefit. In many parts of the world, the positive impact on forest cover and species diversity of the implementation of PES) schemes have been documented [50,51]. Developing PES schemes that specifically target compounding factors contributing to the vulnerability of *P. africana* can contribute to reducing these vulnerabilities. Such schemes would also concomitantly contribute to combatting the degradation or loss of essential ecosystems and ecosystem services without sacrificing the well-being of people—an essential element in the portfolio of Tanzania's sustainable development goals. It must be noted however that the design and implementation of PES programs must be carefully done, adopting best practices as well as the best and most recent scientific guidance on the subject matter. This is because PES per se is not a panacea for addressing underlying deficiencies in natural resources governance policies and practices where they exist. For example, Tuanmu Viña [50] observed that the effectiveness of a PES program depends on who receives the payment and on whether the payment provides sufficient incentives.

#### 4.3.4. Incorporating Forest Management into the Circular Economy

In striving to incorporate forests and forestry into the circular economy, it is important to understand some of the strategies required to achieve the transition from a linear to a circular economy. In a study aimed at understanding strategies for an effective transition into sustainable forest-based bioeconomy in Italy, Falcon, Tani [52] identified that four strategies are most effective. These include defining viable methods of circular management to improve environmental and forest planning; investing in forest infrastructure; supporting entrepreneurship programs for professionals in the forestry sector; and enhancing the development and application of innovative forest-based value chains.

One of the key challenges to transitioning the exploitation and use of *P. africana* in Tanzania from a linear to a circular economic model is that of adding value to the main product, as well as to bio-residuals. Much of the chain for its value-addition for its many uses (medicinal, furniture, fuelwood, etc.) is minimal to non-existent. This is in line with the recognition of that weak market pull, needs for big investments, and the adoption of risk-averse approaches among the few incumbent firms in the sector are reducing the potential for the forestry industry to invest in technological and market capabilities for valorizing residuals [53].

#### **5. Conclusions**

Bioclimatic predictors mainly mean annual temperature presented high contribution and important information in predicting distribution and mapping habitat suitability for *P. africana* in Tanzania. Suitable habitats for *P. africana* will decline in mid- and late-century for both RCP 4.5 and RCP 8.5 scenarios when compared with baseline conditions. For instance, southern highlands and northern zones will constantly lose much more suitable habitats for *P. africana* in the future. The areas mapped in this study as suitable habitats for the tree species could be advantageous for conservation planning and afforestation interventions.

**Author Contributions:** R.A.G. and G.T.Y. contributed equally in the following tasks: conceptualization of the study; methodology development; investigation, data analysis; writing—original draft preparation; as well as the review and editing of the final draft. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was not funded by any external party.

**Acknowledgments:** The authors acknowledge the occurrence data providers Tanzania Forest Service Agency, Global Biodiversity Information Facility and TROPICOS. Further, the authors are grateful to Mathew Mpanda for his valuable suggestions and comments on the manuscripts. Finally, the authors are grateful to Elikana John for his guidance to access the soils data.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **Appendix A**

**Figure A1.** Distribution of the selected environmental variable in the study area.

#### **References**


© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
