Conservation of Tree Species Richness in a Traditional Agroforestry Landscape in the Vhembe Biosphere Reserve, South Africa

: Tree species richness is a critical element concerning trees on farms, on communal land and in protected areas to support biodiversity and socio-economic livelihoods in traditional agroforestry landscapes. Tree species richness is directly linked to the use of provisioning ecosystem services and to management practices in traditional agroforestry landscapes. The study aimed to investigate the link between socio-ecological and conservation strategies regarding tree species richness in traditional agroforestry landscapes. The study was conducted in the Damani, Thenzheni, Tshiombo and Tshipako villages located in Thulamela Municipality of the Vhembe Biosphere Reserve, South Africa. The data were collected using a mixed method approach combining forestry inventory and focus group discussion. The study recorded a total number of 126 tree species: 83 communal-land-hosted species, 68 species of trees on farms and 81 species in the protected areas. The indigenous species Englerophytum magalismontanum (Sond.) T.D.Penn. was the most cited (62%) by interviewees, with a primary use for wild fruits, followed by Pteleopsis myrtifolia (M.A. Lawson) Engl. & Diels. (57%) for fuelwood, Combretum molle R.Br. ex G.Don (36%) for traditional medicine and Albizia adianthifolia (Schumach.) W.F.Wight (12%) for fodder. Species richness was found to be commonly driven by provisioning ecosystem services with trees on farms and on communal land. Distance was found to be major driving factor of species richness in protected areas. This study found that the local people have no conservation strategy and practices targeting the enhancement of tree species richness in the traditional agroforestry landscape. This study advocates for the establishment of a conservation strategic framework for restoring tree species richness by targeting traditional agroforestry landscapes. trees on farms, live fences, protected areas, communal land use areas and forests [23]. The local people of VBR conserve indigenous tree species within and outside the homesteads. Such indigenous trees are commonly used for timber, medicine, wild fruits, fuelwood, fodder, wild food, timber, soil improvement, fencing and building materials [22].


Introduction
The majority of people in Limpopo Province in South Africa depend on forest ecosystems for their daily livelihoods [1,2]. This close dependency on the forest ecosystem is largely influenced by socio-economic and cultural factors [1]. The Vhembe Biosphere Reserve (VBR) is one of the tropical climatic zones in Limpopo Province and is known for its few remaining forest patches and its richness of floristic diversity [3]. A larger portion of the population in the biosphere live in rural areas [2] and sustain their livelihoods by relying on provisioning ecosystem services (PES) such as the harvesting of natural resources for in South Africa is exerting more pressure on TALs, indigenous forests and woodlands [18]. However, the extent of resources used and the impact thereof on biodiversity within TALs are not well documented. The main objective of the study was to investigate the link between socio-ecological and conservation strategies on tree species richness in traditional agroforestry landscapes. The leading questions were: (i) Are there differences in tree species richness between different land uses within TALs? (ii) Is there a link between conservation strategies and tree species richness in TALs? (iii) Do socio-ecological factors have an impact on tree species richness in TALs?

Study Area
The study area was conducted in the Thulamela local municipality in the south-eastern part of the Soutpansberg Mountain range in the Vhembe Biosphere Reserve (VBR), South Africa. The VBR was recognised and designated in 2009 by the United Nations Educational, Scientific and Cultural Organisation (UNESCO) [19]. The Soutpansberg Mountain range is classified as the centre for botanical endemism in southern Africa and the biodiversity hotspot of South Africa [20,21]. The study was conducted in four villages (Damani (22 • • 28 59 E)) located in Thulamela Municipality, Vhembe District, Limpopo province, which is part of the VBR (Figure 1). The study area is characterised by tropical climatic conditions with high annual precipitation, which favours both agriculture and forestry activities [22]. The traditional agroforestry landscape of the VBR consists of different agroforestry practices and land uses including alley cropping, home-gardens, windbreaks, silvopasture, trees on farms, live fences, protected areas, communal land use areas and forests [23]. The local people of VBR conserve indigenous tree species within and outside the homesteads. Such indigenous trees are commonly used for timber, medicine, wild fruits, fuelwood, fodder, wild food, timber, soil improvement, fencing and building materials [22]. South Africa [3]. Furthermore, a huge dependency on forest resources by rural communities in South Africa is exerting more pressure on TALs, indigenous forests and woodlands [18]. However, the extent of resources used and the impact thereof on biodiversity within TALs are not well documented. The main objective of the study was to investigate the link between socio-ecological and conservation strategies on tree species richness in traditional agroforestry landscapes. The leading questions were: (i) Are there differences in tree species richness between different land uses within TALs? (ii) Is there a link between conservation strategies and tree species richness in TALs? (iii) Do socio-ecological factors have an impact on tree species richness in TALs?

Study Area
The study area was conducted in the Thulamela local municipality in the south-eastern part of the Soutpansberg Mountain range in the Vhembe Biosphere Reserve (VBR), South Africa. The VBR was recognised and designated in 2009 by the United Nations Educational, Scientific and Cultural Organisation (UNESCO) [19]. The Soutpansberg Mountain range is classified as the centre for botanical endemism in southern Africa and the biodiversity hotspot of South Africa [20,21]. The study was conducted in four villages (Damani (22°50′45 S, 30°31′38 E), Thenzheni (22°49′54 S, 30°28′57 E), Tshiombo (22°48′30 S, 30°30′53 E) and Tshipako (22°51′14 S, 30°28′59 E)) located in Thulamela Municipality, Vhembe District, Limpopo province, which is part of the VBR (Figure 1). The study area is characterised by tropical climatic conditions with high annual precipitation, which favours both agriculture and forestry activities [22]. The traditional agroforestry landscape of the VBR consists of different agroforestry practices and land uses including alley cropping, home-gardens, windbreaks, silvopasture, trees on farms, live fences, protected areas, communal land use areas and forests [23]. The local people of VBR conserve indigenous tree species within and outside the homesteads. Such indigenous trees are commonly used for timber, medicine, wild fruits, fuelwood, fodder, wild food, timber, soil improvement, fencing and building materials [22].

Data Collection-Fieldwork
This study used mixed methods (quantitative and qualitative); this approach helped to understand the socioecological problems under investigation in depth [24,25]. The quantitative method was used to collect vegetation data through inventory [25]. The study sites were selected based on the premise that the study communities practice land management systems that fit with traditional agroforestry landscapes and were in proximity to the biodiversity hot spot of the Soutpansberg Mountain [26,27] in the VBR. Data collection was conducted in all four study communities during August 2021. In each community, three land-uses were identified: (i) trees located within farmed landscapes (hereafter, "trees on farms"), (ii) trees within communal lands and (iii) tribal protected areas (hereafter, "protected areas"). Three transects, each being 150 m long, were established 200 m apart along a linear transect in each land use. Five rectangular plots of 20 m × 10 m (200 m 2 ) in each transect were established for inventory, with a sampling intensity of 10% [28]. The trees were categorised into three categories based on growth form and size: (i) a woody plant with a diameter at breast height (DBH, at 1.3 m) of greater than 5 cm and a total height greater than 2 m; (ii) a shrub with a DBH of less than 5 cm and a height of less than 2 m; (iii) seedlings which were counted in subplots. A calliper and diameter tape (for trees with a girth, >50 cm) were used to measure the DBH, and an electronic clinometer (Haglof, Sweden) was used to measure tree height. The first plot was marked randomly, and then the following plots were marked systematically along the transect line, 25 m apart, corresponding to the first plot. Within each rectangular plot, five 2 m × 2 m sub-plots were established for sampling; only the number of individuals, local names and scientific names were recorded.

Data Collection-Focus Group Discussion
The qualitative method was used to collect data through focus group discussions [29]. Focus Group Discussions (FGDs) were conducted to obtain more data that would help to understand the community members' perceptions concerning tree species richness and conservation strategies. This allowed participants to provide a wide range of information. The FGD participants were a combined group of community leaders and community members. Four FGDs were conducted with four communities, Damani (n = 15), Tshipako (n = 23), Thenzheni (n = 16) and Tshiombo (n = 13). Each FGD session lasted for about 1 h. The leading questions were structured as follows: (i) What are the tree species used for fodder, fuelwood, traditional medicine and food? (ii) What are the conservation strategies for traditional agroforestry and other tree-based systems? Each participant was given an opportunity to discuss and raise questions under the facilitation of the researcher. The discussions from the FGDs were recorded for content analysis.

Vegetation Analysis
Plymouth Routines in Multivariate Ecological Research (PRIMER-e) software: (Version 7.0.21, Auckland, New Zealand) with add-on PERMANOVA+ was used to analyse the biological data. The permutation number to estimate the probabilities was set at 999 [30].

Species Richness in Different Land-Uses
The sampling effectiveness effort was tested to evaluate the sufficiency of sampling [31]; the sufficient sampling effort captures about ≥80% of estimated species richness [6]. The species richness estimate was performed using a bootstrapping technique [30]. Then, the species accumulation curve was computed to demonstrate the species richness; these indices compared the variety of species richness in different land-uses. To determine the significant difference in species richness among the land-uses, the Margalef index (d) for overall species richness was determined using the DIVERSE function in PRIMER-e [30]. The Margalef index was calculated as follows in Equation (1): where S = the species count, and N = the total number of individuals [32]. In a following step, the Euclidean distance sample plots were developed using the resemblance matrix of the Margalef index [6]. Based on the species richness data, significant differences in species richness among different land-uses (trees on farms, communal land and protected areas) were tested with permutation multivariate analysis of variance (PERMANOVA) [30]. Pairwise comparison tests with 999 PERMANOVA permutations and t statistics were conducted to assess the significance of difference between the land-use pairs "trees on farms × communal land, communal land × protected areas, trees on farms × protected areas" for species richness data. Additionally, species richness analysis and assemblage illustration among the land-uses were analysed using non-parametric multidimensional scaling (nMDS) in PRIMER-e using the Bray-Curtis dis/similarity index [33].

The Variance of Species Richness among Different Land-Uses
To assess the significance of variance among the land-uses, the homogeneity of dispersion (PERMDISP) was tested using the Bray-Curtis coefficient matrix with 999 permutations [30]. A pairwise comparison of land-uses was conducted to show variance between land-use groups. Then, a dissimilarity within-group means square was generated using the Bray-Curtis dis/similarity index in PERMDISP [30]. The principal component analysis (PCO) was used to test the homogeneity of variance among land-uses [30]. In this analysis, the distance among land-uses was obtained using the Euclidean distance.

Drivers of Species Richness in Different Land-Uses
The inter-correlations among the change drivers were tested using the Multicollinearity test of the correlation matrix in PRIMER-e Draftsman plot. The pairs of all variables' correlation were below the threshold (|r| ≥ 0.95); therefore, they contain effectively different information and are not redundant [30]. Then, the link between the change drivers and land-uses was tested using distance-based linear modelling (DISTLM) and distancebased redundancy analysis (dbRDA) [6]. A dbRDA plot obtained from PERMANOVA using DISTLM was used to visualise the patterns of the DISTLM results.

Species Richness Distribution in Different Land-Uses
Significant species richness distribution differences among the land-uses were tested using analysis of similarities (ANOSIM), and similarities percentages analysis (SIMPER) was used to identify dominant tree species (cut of 70%) and the abundance contribution to each land-use. The tests were conducted using the Bray-Curtis dis/similarity index [32]. Then, the dominance of each tree species (top 30) was visualised in shade plots.

Social-Ecological Data Analysis-Focus Group Discussion
To understand the perceptions of the local people and the community leadership and the use and management of the tree species richness in TALs, a comprehensive qualitative analysis of focus group discussion data was performed. The focus group discussion participants and the informants' consensus on utility categories were analysed using fidelity level [34] and were calculated as follows in Equation (2): where FL = fidelity level, I p = number of participants that cited the principal use of the species and I u = total number of participants that cited the species for any purpose. Content analysis was performed for the themes or subjects discussed with participants [35]. The categorising of themes or subjects was carried out through manifest analysis [36]. The variability and reliability were assessed through a crossing evaluation of the participant's responses.

Species Richness in Different Land-Uses
A total of 2716 individual trees were recorded, of which 636 were trees on farms, 853 were individual trees on communal land and 1227 were in protected areas. Significant differences among the study villages were tested using analysis of similarities (ANOSIM); the results were statistically not significant (p > 0.5). Therefore, the different land uses were compared. The bootstrapped species accumulation curve demonstrated that the species sampling captured the full range of species within the defined land-use categories ( Figure 2). The variability and reliability were assessed through a crossing evaluation of the participant's responses.

Species Richness in Different Land-Uses
A total of 2716 individual trees were recorded, of which 636 were trees on farms, 853 were individual trees on communal land and 1227 were in protected areas. Significant differences among the study villages were tested using analysis of similarities (ANOSIM); the results were statistically not significant (p > 0.5). Therefore, the different land uses were compared. The bootstrapped species accumulation curve demonstrated that the species sampling captured the full range of species within the defined land-use categories ( Figure  2). This study recorded a total number of 126 tree species, of which communal land, farms and the protected areas hosted 83, 68 and 81 species, respectively ( Table 1). The sampling effectiveness of the whole study was calculated as 92% of the effort, and in each land-use, trees on farms, communal land and protected areas were individually 84%, 87% and 96%. This study recorded a total number of 126 tree species, of which communal land, farms and the protected areas hosted 83, 68 and 81 species, respectively ( Table 1). The sampling effectiveness of the whole study was calculated as 92% of the effort, and in each land-use, trees on farms, communal land and protected areas were individually 84%, 87% and 96%.
The PERMANOVA test detected a significant difference in species richness among land-uses (F 18.015 ; p ≤ 0.001), while there were no significant differences among the villages (F 1.9945 ; p = 0.016) ( Table 2). Then, a pairwise comparison was conducted to compare the species richness significance difference between the land-uses. The species richness significantly differed across all the land-uses (Trees on farms and Communal land, t = 4.5, p ≤ 0.05; Trees on farms and Protected areas, t = 5.5, p ≤ 0.05; and Communal land and Protected areas, t = 2.3, p ≤ 0.05). The average similarity showed that there was a high similarity between the communal land and protected areas (25%), followed by a medium-level similarity (21%) between trees on farms and communal land and a low average similarity between trees on farms and protected areas (17%) ( Table 3).

The Variance of Species Richness among Different Land-Uses
The homogeneity of multivariate dispersion (PERMDISP) results showed a significant difference in species variation among the land-uses (F = 15.79, p ≤ 0.001). The PERMDISP pairwise comparisons support the significant difference in land-uses ( Table  4). The trees on farms category differed significantly from communal land (t = 4.4022, p ≤ 0.001). Likewise, the trees on farms differed significantly from protected areas (t = 4.2802, p ≤ 0.001). The species variation of communal land and a protected area was not statisti-

The Variance of Species Richness among Different Land-Uses
The homogeneity of multivariate dispersion (PERMDISP) results showed a significant difference in species variation among the land-uses (F = 15.79, p ≤ 0.001). The PERMDISP pairwise comparisons support the significant difference in land-uses (Table 4). The trees on farms category differed significantly from communal land (t = 4.4022, p ≤ 0.001). Likewise, the trees on farms differed significantly from protected areas (t = 4.2802, p ≤ 0.001). The species variation of communal land and a protected area was not statistically distinct (t = 0.45084, p = 0.673). The species richness variation average dissimilarity within the land-uses was 45% at trees on farms, 52% on communal land and 51% in protected areas ( Table 4). The PCO explained 27% of the variation among the land-uses-9.5% and 17.5%, respectively, between axes ( Figure 4).

Drivers of Species Richness in Different Land-Uses
The DISTLM results showed that there was a significant correlation between species most of the socio-ecological factors examined (R 2 = 0.194, p ≤ 0.001). The marginal te significant relationship (p ≤ 0.001) between species richness (all land-uses) and most soc

Drivers of Species Richness in Different Land-Uses
The DISTLM results showed that there was a significant correlation between species richness and most of the socio-ecological factors examined (R 2 = 0.194, p ≤ 0.001). The marginal test showed a significant relationship (p ≤ 0.001) between species richness (all land-uses) and most socio-ecological factors (elevation, gradient, PES harvesting, grazing and distance), except for fire occurrence (p = 0.384). The variations explained by each variable are as follows: distance (12.2%), elevation (2.6%), gradient (5.5%), fire occurrence (0.6%), grazing (1.6%) and PES harvesting (1.5%) and as shown in Table 5. The first two axes of dbRDA captured 83.2% of the variability in the fitted model and 16.2% of the total variation in the data cloud ( Figure 5). Concerning trees on farms, species richness is driven by fire occurrence, grazing and PES harvesting. On communal land, species richness is driven by PES harvesting, gradient, distance, elevation and grazing. The tree species richness in the protected area is driven by distance.

Species Abundance Distribution in Different Land-Uses
The SIMPER results showed the dominating species and contribution of each species in different land-use abundances ( Table 6). A full species list is given in Table S1 in the Supplementary Materials. Non-native tree species have largely dominated and contrib-

Species Abundance Distribution in Different Land-Uses
The SIMPER results showed the dominating species and contribution of each species in different land-use abundances (Table 6). A full species list is given in Table S1 in   An average dissimilarity of 97% species richness was observed between the trees on farms and on communal land (Table 7). A total of 27 of the 126 shared species between the trees on farms and communal land contributed above 70% to the average dissimilarity between the two land-uses. The most dominating species in terms of dissimilarity were An average dissimilarity of 89% species richness was observed between the communal land and protected areas. A total of 30 of the 126 shared species between the communal land and protected areas contributed above 70% to the average dissimilarity between the two land-uses. The most dominating species in terms of dissimilarity were Pteleopsis myrtifolia (6%), Tabernaemontana elegans (5%), Parinari curatellifolia (4%), Pseudolachnostylis maprouneifolia Pax. (3%) and Combretum molle (3%). The indigenous species Faurea saligna was completely absent in the communal land. The visual impact of the greyscale intensities in a shade plot ( Figure 6) gives a strong idea of which species are highly dominated in each land-use. White denotes the absence of that species in that land-use, and full black represents the maximum abundance in the land-use. The shade plot shows the top 30 dominating tree species; the trees on farms were dominated by non-native tree species, while both the communal land and the protected areas were dominated by indigenous species. Notably, only one endemic species was found in the top 30 and was observed only in protected areas.

Tree Species Use and Conservation Strategies
In the FGDs, the participants listed the important tree species based on the utility categories such as fuelwood, fodder, wild food and traditional medicine. Based on species nomination within the FGDs, the most important tree species in each category were evaluated using the fidelity level formula. A full species list is given in Table S1 in the supplementary material. The overall number of citations of the most important tree species varied between 4% and 62% in the use categories (Figure 7). The indigenous species Englerophytum magalismontanum was the most reported tree species (62%) for wild fruits in the study sites, while it was also highly reported for fuelwood utilisation (43%). Pteleopsis myrtifolia was the most described (57%) tree species for fuelwood, while Combretum molle was the most cited (36%) tree species for traditional medicine. Albizia adianthifolia (Schumach.) W.F.Wight was most frequently given (12%) as a species for fodder. From 29 tree species, 8 tree species were cited as multiple purpose species, such as Englerophytum magalismontanum and Parinari curatellifolia (both appearing in the wild food and fuelwood utility categories). Peltophorum africanum Sond. and Combretum molle were suggested in both the traditional medicine and fuelwood utility categories. Artabotrys brachypetalus Benth. and Annona senegalensis Pers. were suggested in both the wild food and traditional medicine utility categories, Diplorhynchus condylocarpon Müll.Arg. was suggested in both the fodder and wild food utility categories and Albizia adianthifolia was recorded in both the fuelwood and fodder utility categories.

REVIEW
14 of 21 Figure 6. Shade plot showing tree species distribution between land-uses, further subdivided by native or non-native status in South Africa.

Tree Species Use and Conservation Strategies
In the FGDs, the participants listed the important tree species based on the utility categories such as fuelwood, fodder, wild food and traditional medicine. Based on species nomination within the FGDs, the most important tree species in each category were evaluated using the fidelity level formula. A full species list is given in Table S1 in the supple-  Overall, in the FGDs sessions, three themes (national conservation policy, community conservation policy and conservation practices) were selected for an in-depth understanding of the conservation strategies.  Overall, in the FGDs sessions, three themes (national conservation nity conservation policy and conservation practices) were selected for an standing of the conservation strategies.

Theme 1: National Conservation Policy
In all study sites, the participants are not aware of conservation po they knew that the government protects natural forests and trees. a partic village said: "A few decades ago, there were forest guards in our communitie government to protect forests and trees. But

Species
Wildfood Traditional medicine Fodder Fuelwood Figure 7. Nominated important tree species per defined utility category, as cited by participants in the focus group discussions.

Theme 1: National Conservation Policy
In all study sites, the participants are not aware of conservation policies. However, they knew that the government protects natural forests and trees. a participant at Damani village said: "A few decades ago, there were forest guards in our communities employed by the government to protect forests and trees. But, they are no longer there and our forests are overexploited". The challenge of lacking knowledge of conservation policies was further emphasised by the chief of one of the communities, who said: "We just hear that are some of the trees which are protected by government policies, but we don't know these trees".

Theme 2: Community Conservation Policy
In all the study sites, the Chiefs or community leadership reported that it is their responsibility to protect the forest trees within their jurisdiction. Although there are no documented policies or rules provided by the traditional leaders, the message is conveyed through word-of-mouth during the communities' gatherings. An FGD participant said: "Our leadership always reminds us to protect and report a person who fells live trees" Furthermore, "People from outside our communities are not allowed to harvest or collect forest re-sources in our forests". The leadership said: "We have people who harvest fuelwood for selling, we always advised them not to harvest in bulk and they should frequently change harvesting sites to allow regeneration".

Theme 3: Conservation Practices
In all the study sites, there is no formal project implementing conservation practices. The participants indicated: "We are protecting the existing trees on our homesteads and communal farms, and plant non-native trees such Mangifera indica". The leadership said: "We manage indigenous trees in our farms and homesteads by doing minimal harvesting instead of felling the whole tree". Furthermore, "We are interested in planting indigenous trees in our communal land. However, we do not know how we can plant indigenous trees". Trees within the community are highly managed by the community members because of the benefits provided by the trees. A participant said: "Trees help us during natural disasters floods, storms and drought".

Discussion
This study utilised both quantitative and qualitative methods, employing a mixed method approach. The quantitative method was used to collect vegetation data through inventory for an assessment of biodiversity and species richness within the study areas. The qualitative method was used to collect social data through FGDs. FGDs were critical in this study to gather more detailed information of how the local people conserve and understand the conservation of trees in their traditional agroforestry landscape. Behaviours that can be explained by social science methodology are integral when evaluating agroforestry systems and can complement the empirical data collection rooted in natural science practices. The combination of inventory and FGDs helped to show a link between the species richness and conservation strategies being applied by the local people. Such an approach also helped in understanding changes in species richness, in combination with the vegetation survey data concerning the three categories: trees on farms, communal land and protected areas.

Species Richness in Different Land-Uses
The tree species richness was higher on communal land than it was in protected areas and on farms. These findings coincide with the results of the studies conducted by Tezbew and Asfaw [37] and Eyasu et al. [38], where they reported that human-modified landscapes harbour over 90% of the tropical species richness of protected areas. However, due to anthropogenic activities, the tropics are experiencing deforestation and forest degradation [38]. The high number of tree species in communal land in comparison to protected areas and trees on farms was associated with the disturbances in communal land. The intermediate disturbance hypothesis suggests that the species diversity maximises when a disturbance occurs at an intermediate frequency or with an intermediate intensity [39,40]. This may show that communal land could harbour a large number of species and disturbances due to different activities; therefore, such land could serve to increase species richness. Nevertheless, the high disturbance is associated with the loss of the tree species richness in natural forests [40]; therefore, the magnitude of disturbance has to be monitored to help to conserve the tree species richness in communal land. Other studies conducted in different parts of the world also reported similar findings; higher species numbers are found on communal land than within protected areas or associated with trees on farms [41]. The lower level of tree species richness concerning trees on farms was associated with land users and managers planting more non-native species and retaining few indigenous tree species. Local farmers tend to keep more species that provide multiple important economic benefits [42]. The value of tree species' heterogeneity in trees on farms is regarded as a long-term benefit [43,44]. The limitation of planting more indigenous tree species is a delayed-return-investment from indigenous trees [45]. However, Borsari [46] argued that agriculture production is heavily dependent on the diversity and abundance of existing and planted tree species on a farm. The integration of diverse native tree species enhances the landscape and habitat heterogeneity [37]. Modern agriculture homogenised most of the heterogonous landscape by maximising the production at the expense of the tree species heterogeneity [43]. A traditional agroforestry practice such as trees on farms has been appreciated in terms of enhancing tree species heterogeneity on farms and in promoting conservation [41]. Trees on farms provide similar resources as native forests; for example, the local people of Nepal collect 43% of the fuelwood and fodder from trees on farms [41]. Overexploitation has a negative impact on the richness of trees on farms, resulting in poor production, thus affecting local livelihoods [47]. In protected areas, the lower level of tree species richness is highly likely to be associated with succession, because the protected areas are intact and not disturbed due to restrictions [23]. The results on the tree species richness similarity of the land-use indicated that the communal land and protected area are similar in terms of floristic richness and harbour similar tree species. This fits well with the notion that trees on farms harbour most non-native tree species; hence, the similarity is very low within communal land and protected areas.

Drivers of Species Richness Variance in Different Land-Uses
The distance between households and land-uses such as communal land, trees on farms and protected areas has a positive impact on species richness, which is in line with the findings of Bhandari et al. [41], who reported that the distance between forests and people's homes influences the species richness in agroforestry landscapes. When there is a large distance between the forest and the people, people save the time and costs of the regular collection of forest products by collecting in nearby sites such as from trees on farms and on communal land [39,40], and they sometimes tamper with protected areas [23]. Hence, the current study observed human disturbances in protected areas. It is more likely that the human disturbances in protected areas are unauthorised. Furthermore, the distance compels the local people to harvest the forest products from the trees on farms and on communal land [41]. The close dependence on communal land and trees on farms reduces the pressure on native forests and eventually contributes to conservation [41]. A great species richness in communal land and on farms is not only contributing to conservation but also contributes to the stand productivity and diversity of the forest products [38]. The tree species richness in communal land and on farms provides a diversity of forest products such as fruits, fodder, construction materials, food and firewood [37].

Tree Species Richness, Uses and Conservation Strategies
This study found that the TALs of the VBR host up to 29 of the most appreciated tree species, of which 15 are used for fuelwood, 14 are used for wild food, 7 are used for traditional medicine and 2 are used for fodder. On the contrary, studies conducted in traditional agroforestry systems in India found that most tree species are used for fuelwood, followed by traditional medicine uses [48]. Unsurprisingly, most of the tree species in TALs have a high value to the local people [14], determined by human needs such as food, energy and shelter. Concerning trees on farms, the non-native species contributed the greatest number of tree species, which were represented by two species, Mangifera indica and Persea Americana. These findings coincide with the results of the study conducted by Yashmita-Ulman [49], who reported that farmers in India plant Mangifera Indica along the farm boundaries or within the farms. A study conducted by Borsari [46] reported that Mangifera indica is one of the iconic tree species in agroforestry practices such as silvopasture and trees on farms. Both the communal land and protected area were dominated by several indigenous tree species such as Pteleopsis myrtifolia, Albizia adianthifolia, Tabernaemontana elegans, Parinari curatellifolia, Englerophytum magalismontanum, Bridelia micrantha, Pseudolachnostylis maprouneifolia and Combretum molle.
This study found that the local people have no conservation strategy and conservation practices to enhance the tree species richness for the sustainable functionality of the traditional agroforestry landscape. The only strategy that seems to be effective is the set of traditional rules of protecting trees by restricting people to felling live trees. The study conducted by Duriaux-Chavarria et al. [47] highlighted that, in Niger and Ethiopia (north and central), rural communities and farmers restored tree species richness through the planting of native tree species and subsequently improved tree species richness and ecosystem services. Considerable restoration practices such as planting native tree species on traditional agroforestry land reduce disservices by human beings.

Implications for Conservation
The findings revealed a lower tree species richness among trees on farms compared to the communal land and protected areas. The participants acknowledged and highlighted the lack of awareness of conservation policies; others highlighted the lack of knowledge of planting native trees in the landscape. This may have an overall negative impact on rural forest tree conservation and the enhancement of tree species richness in traditional agroforestry landscapes. Therefore, identifying native tree species that can be planted to enhance forest tree species richness and diversity is paramount in this traditional agroforestry landscape. The Climate Smart Agriculture Strategic Framework for Agriculture, Forestry and Fisheries in South Africa developed actions to increase the forest cover or tree species in agroforestry systems by identifying agroforestry species for different agroecological zones and supporting farmers in increasing the tree cover on their land [15]. However, practically, nothing has been achieved, as local people are still only depending on the existing tree species. An additional constraint is that the developed Agroforestry Strategic Framework [15] largely promotes modern agroforestry systems over traditional agroforestry. Interdepartmental collaboration in enhancing traditional agroforestry tree species diversity would play an important role in conserving and improving the local people's livelihoods. Moreover, implementing the actions and strategic frameworks would fast-track the mitigation of biodiversity loss in South Africa.

Conclusions
The study found that the TALs were characterised by multipurpose tree species such as Englerophytum magalismontanum, Parinari curatellifolia, Peltophorum africanum, Combretum molle, Artabotrys brachypetalus, Annona senegalensi, Diplorhynchus condylocarpon and Albizia adianthifolia for different uses such as fuelwood, fodder, traditional medicine and wildfruits. Species richness was found to be higher on communal land commonly driven by PES. A reduced overharvesting of PES harvesting, in combination with the effort to conserve trees, would, therefore, restore and augment species richness in both trees on farms and communal land uses. The physical distance between resources and village settlements was the major factor driving species richness within protected areas. Therefore, there is a greater need for protected areas in close proximity to communities with strict traditional rules to conserve tree species richness. The results of this study confirmed that local people are promoting the conservation of native tree species within traditional agroforestry landscapes and are also conserving the existing native tree species governed by local and traditional rules. In addition, there were no forest restoration or conservation practices or projects implemented in the communities to support biodiversity and species richness. Lastly, an inadequate legislative framework on conservation in rural communities is a critical challenge considering the unprecedented rise of biodiversity loss and forest degradation. This study advocates for the adoption of native tree species planting and conservation in traditional agroforestry landscapes to enhance the landscape productivity and functionality. There is, therefore, a need to establish a conservation strategic framework for restoring tree species richness, targeting traditional agroforestry in the VBR region and beyond.