The Pongola Bush: Tree diversity assessment in a KwaZulu-Natal forest patch

recruitment rate of small stemmed individuals and few large individuals. This pattern is consistent with the disturbance history of the site, with limited recruitment of certain species, mainly limited to early successional species. The Pongola Bush forest is particularly diverse in terms of tree species: 41 species were recorded, and it is estimated from the species area curve that there may be 70 tree species present. Conservation implications: This survey is the first formal tree diversity survey conducted on the Pongola Bush which may assist future research and conservation strategies.


Introduction
Roughly half of the world species are found in forests and they continue to be exploited for their socio-economic value, including the need for industrial timber, medicinal products, and use of productive agricultural land (Grainger 2013). This exploitation often leads to the transformation of systems, negatively influencing diversity and ecosystem functionality (Falk, Palmer & Zedler 2013;Hooper et al. 2005;Johnson et al. 1996). Subject to the extent of transformation, ecosystems can regenerate after disturbance events (Vieira & Scariot 2006). Regeneration is dependent on a variety of factors including seed dispersal and seedling establishment, structure and richness of existing species, and resprouting ability (Vieira & Scariot 2006). Identifying natural ecosystem restoration processes in forests is challenging and requires systematic assessment of structural composition (Chazdon et al. 2016).
Forests are scattered along eastern and southern margins of South Africa, from the Soutpansberg in the north and Maputaland in the east to the Cape Peninsula in the West (Mucina & Rutherford 2006a). With their patchy distribution within more open systems, limited to fire protected refugia and steep south facing slopes, indigenous forests account for less than 0.5% of the land cover in South Africa (Berliner 2009;FAO 2010;Geldenhuys 1994;Huntley 1984;Rutherford, Mucina & Powrie 2006). Mucina and Rutherford (2006a), using a floristic biogeographic classification of the indigenous forests of South Africa, recognised 26 forest types, grouped into eight zonal groups and one azonal group. Afrotemperate forests stretch from the Cape Peninsula via the Southern Cape Coast, the Eastern Cape and KwaZulu-Natal (KZN) Midlands and Drakensberg to the North-eastern Escarpment and Soutpansberg Mountains. Each of these forest types can have both high and scrub forest depending on climatic conditions and landscape. Afromontane forests are generally found on steep south-facing slopes with lowland forests on the coast having a lower canopy.
The Pongola Bush is placed in the Northern Afrotemperate zone by Mucina and Rutherford (2006b), but Lötter, Mucina and Witkowski (2014) provided an updated classification of the The Pongola Bush Nature Reserve lies in a narrow band along the escarpment between Mpumalanga and KwaZulu-Natal (KZN). Although referred to as 'bush', this vegetation type may be considered to fit into the intersection of two types of scarp forest: the Northern Afrotemperate Forest and the Southern Mistbelt Forest. The area was heavily logged in the early 1900s with the need for timber for the emerging gold mining industry in Barberton. The forest fragment is now protected but lacks formal tree diversity assessments, and this study sets out to establish its community composition and richness using a range of standard techniques. Distributed throughout the forest, 127 circular plots (each 80 m 2 ) were surveyed. In these, 1152 stems were measured for species, height, and diameter at breast height (DBH). Four species contributed roughly 70% of the 70 Mg ha −1 above ground woody biomass. Stem size class distributions showed a high recruitment rate of small stemmed individuals and few large individuals. This pattern is consistent with the disturbance history of the site, with limited recruitment of certain species, mainly limited to early successional species. The Pongola Bush forest is particularly diverse in terms of tree species: 41 species were recorded, and it is estimated from the species area curve that there may be 70 tree species present.  Lötter et al. (2014). The Pongola bush is believed to be a remnant of ancient forests (perhaps from the Tertiary), as evidenced by the presence of Onychophora (velvet worms or peripatus), for instance, that flourished in the area (Ruhberg & Hamer 2005). It does not have a formal inventory of tree species or community structure. This pilot study determined the tree diversity and community composition which would provide a foundation for a more extensive study on changes in the floristic-structural composition along an altitudinal gradient. The diversity and composition was achieved through the documentation of species and species metrics to quantify the richness and relative abundances as well as the total above ground biomass of species within the forest. This study also investigated whether there would be a high dominance of early successional species.

Study site
The Pongola Bush Nature Reserve covering 8.82 km 2 (882 ha) is managed by Ezemvelo Wildlife (KwaZulu-Natal Nature Conservation Board). It is situated immediately south of the border of Mpumalanga and KZN (Figure 1; 27.32349S, 30.47642N). The Pongola Bush is the largest remaining forest fragment in the area, which historically provided timber for the development of both the Barberton and Witwatersrand goldfields. Sawpits and wagon roads are present in the upper reaches of the forest, suggesting an extensive logging history. Much of the remaining extent of the Pongola Bush is situated within the reserve boundaries: an additional 3 km 2 lies outside the borders of the reserve. The forest extends less than 1 km wide in a north-south band, and several kilometres in an east-west direction. The forest lies along the southern edge of an escarpment cliff, extending down a steep slope to grassland on a gentle gradient, marking the southern forest edge. The vegetation is tall (> 10 m) and has a closed, interlocking canopy with multiple layers. The climate in the area is temperate, with an average annual temperature of 14.1 °C, and is characterised by wet summers (November-February) and dry winters (May-August) with an average annual precipitation of 790 mm.

Species data collection
Transects were walked from the southern edge of the Pongola Bush patch towards the northern cliff base. Approximately every 50 m along each transect that was on average 600 m long (450 m -800 m), a circular plot with a 5 m radius was surveyed. A total of 127 circular plots were sampled through the forest. All tree stems greater than 20 cm in circumference (diameter ≥ 6.37 cm) at breast height (1.3 m) where the centroid of the stem cross-section is within 5 m of the plot centre were identified, and the height of the top of the crown of that individual was estimated using a ranging rod and clinometer with a precision of about 1 m.

Statistical analyses
Species biomass values were estimated using the above ground biomass (AGB) model of Chave et al. (2014): Source: Base-layer map created using ArcGIS® software by Esri® where ρ is wood density in Mg m −3 , D is stem diameter at breast height (DBH) in m and H is height in m. Wood densities of the species encountered were taken from the African Wood Density Database (Carsan et al. 2012), shown in Appendix 1 (Table 1-A1). Species were ranked based on their total above ground biomass (in Mg), and statistical significance for biomass measures among and between species was determined using an ANOVA and a post-hoc Tukey test at a 95% confidence interval (α = 0.05). Stem density (stems ha −1 ) was determined for all species using count data within the known plot area. Stem size-class distributions (in terms of DBH) were determined for the entire sampled population as well as the four most abundant species, plus an economically valuable species (Afrocarpus falcatus Thunb.). Species area relationship analyses were conducted using the 'sars' (species-area relationships) software developed by Matthews et al. (2019). The speciesarea model (the increase in species richness recorded as the sampled area increases) was computed using eight random bootstrapped samples from the plot data, which were then averaged.
Biodiversity and community structure estimates were conducted using Paleontological Statistical Software Package (PAST) 1.81 (Hammer, Harper & Ryan 2001). Diversity indices were calculated in two ways: firstly, for number of individuals, which is most appropriate for population dynamics interpretations, and secondly, based on contribution to biomass, which is more important for functionality (Mason et al. 2013;Schleuter et al. 2010). Three diversity indices were calculated: dominance (1-Simpson), Simpson 1-D and Shannon (entropy). Community composition analyses using ordination techniques were conducted using a presence-absence matrix to run a non-metric multidimensional scaling (NMDS). The similarity index input for the NMDS was the Bray-Curtis index, as it is the most valuable metric for showing the compositional dissimilarity of sites using abundance data (Gauch 1973).

Ethical considerations
This article followed all ethical standards for research without direct contact with human or animal subjects.

Results
A total of 1152 stems with a circumference greater than 20 cm were sampled, in an aggregated area comprising one ha, from a forest with a total area of about 882 ha. The sampled stems represented 41 tree or tall shrub species (Table 1).  Figure 2). Eleven of the species contributed 95% of the AGB with over 30 other species accounting for the remaining 5% of the Pongola Bush tree AGB ( Figure 2).
There seems to be no clear dominance in terms of stem density within the Pongola Bush ( Figure 3). The four species which constituted the greatest biomass are not the four species with the greatest stem densities, but they do appear in the quartile of species with the greatest stem densities. Gymnosporia mossambicensis (Koltzsch), Afrocarpus falcatus, Cussonia spicata (Thunb.) and Diospyros whyteana (Hiern) contributed the largest number of stems per unit ground area ( Figure 3). Gymnosporia mossambicensis has the greatest stem density yet it does not fall into the list of species which contribute to the greatest AGB, because it is a multi-stemmed tall shrub, rather than a canopy-forming tree. The stem size distributions (in DBH) of all individuals of all species, and the stem size distributions for species with the highest biomass and a previously economically important species (A. falcatus) are shown in Figure 4. Stem diameters, both in aggregate and for some individual species populations, follow a classical 'reverse J' distribution -a more-or-less smooth exponential decrease in numbers with increasing stem diameter. This indicates steady recruitment following a period of past disturbance, and a constant stem mortality rate. The populations of S. chirindensis, C. kraussii and K. africana all show a hump-shaped distribution, with few small (< 20 cm) and very large (> 80 cm) stemmed individuals. The bulk of the population have medium stem diameters (20 cm -60 cm), indicating an episode of strong recruitment a century ago, followed by reduced recruitment. This is contrasted by the size class distribution of A. falcatus and C. africana, which both have populations dominated by smaller stemmed individuals and fewer larger individuals, suggesting that the demographic is in a recruitment and regeneration phase, following exploitation in the past. For all five of the chosen species, there are individuals within each population which have very large DBHs compared to the rest of the population, suggesting these are the remnants of large individuals that were established before the timber exploitation era, and were missed in harvesting operations.
The significant relationship between species richness and area, expressed as a log-shaped curve, shows the asymptotic species-area boundary of this curve has not been reached at 1 ha (10 000 m 2 ), the upper limit of the aggregate sample ( Figure 5; r 2 = 0.84 , p < 0.01). At this point, the graph still has an upward trend. Extrapolating the relationship to the full extent of the Pongola Bush (~882 ha) suggests that there may be 70 woody plant species in the forest as a whole, thus about 30 undiscovered by this survey.
Diversity indices, calculated on two fundamentals, are shown in Table 2. The dominance index ranges from 0 to 1 with 1 being total dominance by an individual species, in this case dominance was lower for species diversity (0.065) compared to functional diversity (0.1492), which is supported by the biomass dominance by a small number of species shown in Figure 3, and little dominance in Figure 4, in terms of the gradual decrease of stem densities with no clear dominance. The Simpsons index for individual counts (0.94) shows a high diversity and evenness (the scale is 0-1, with 1 high richness and evenness). When expressed as functional diversity, there is still a high degree of richness and evenness (0.85). The upper limit of the Shannon index is unbounded, and shows the species diversity (2.94) to be higher than that of the functional diversity (2.25).
One invasive species was identified in this study, Solanum mauritianum Scoop., but few individuals existed and were confined to areas where the canopy had been disturbed and establishment could take place, suggesting that the potential of rapid colonisation is unlikely given the interlocking canopy of the Pongola Bush. Upon calculating community composition using multi-dimensional scaling based on the Bray-Curtis similarity coefficient, there appeared to be no clear tree compositional changes along the forest gradient. In the field there were apparent compositional changes of canopy floor species (forbs) and species which fell outside of the threshold stem diameter; thus these were not explicitly addressed in this study. This means that there are no predictable communities of certain species based on depth into the forest (or equivalently, altitude, because the slope is always upward with distance into the forest), and that the spread of species is relatively even throughout the forest.

Discussion
The biomass is higher than those of fragmented dryer Afromontane forests found in Kenya and other regions (~50 Mg ha −1 ), but lower than tropical African forests (150 Mg ha −1 -300 Mg ha −1 ) in the Albertine Rift region (Adhikari et al. 2017;Baccini et al. 2008;Glenday 2008;Imani et al. 2017). The Pongola Bush is similar in biomass compared to African scarp subtypes at roughly 70 Mg ha −1 (Glenday 2008), but, as expected, only roughly a third of the Note: If the letters above the bars differ then the total above ground biomass is significantly different across the species (p < 0.05). AGB, above ground biomass; Mg, megagram.  biomass of pristine mature Afrotemperate forests (Adie, Rushworth & Lawes 2013). Biomass in the Pongola Bush is not spread evenly among species but rather is dominated by four species; this dominance mirrors some of the other Afrotemperate studies (Adie et al. 2013). Species that are commonly found in Northern Afrotemperate forests include many that were found in this study including, Celtis, Afrocarpus, Halleria and Scolopia; those missing were Olinia, Pittosporum and Rothmania. When extrapolated, the species area relationship of roughly 70 woody species in the forest patch as a whole is consistent with other Afromontane forest patches of similar size (50-110 species) (Gemeda, Lemenih & Gole 2018;Giliba et al. 2011;Mensah et al. 2020). These Northern Afrotemperate forest patches are generally dominated by late successional, canopy forming Podocarpus but due to relatively recent disturbance events common genera such as Combretum and Celtis have taken their place (Asrat et al. 2020) resulting in different types of forest.
Diversity indices coincide with forests of similar composition in areas around Africa, particularly Afromontane forests. In terms of dominance, Afromontane forests in Tanzania have  Note: Log-log relationship assumes a form S = bAz, where S is richness, A is area, and b and z are constants.   (Giliba et al. 2011). In the cases where there was higher dominance, these were due to disturbances experienced by the forest subtypes (Mafupa 2006). Species richness is comparatively higher per ha in the Pongola Bush compared to patches seen both in specific studies in Ethiopia and Tanzania (Gemeda et al. 2018;Giliba et al. 2011), yet lower than southern African forests of like composition (Adie et al. 2013). Comparisons with other South African forests have also been made (Geldenhuys 1992;Mucina & Rutherford 2006). Shannon indices were generally lower in Ethiopian Afromontane forests (2.67-2.87). This may be due to various factors including disturbance and limited rainfall (500 mm per annum) compared to the Pongola Bush (Gemeda et al. 2018;Mekonnen 2006). Geldenhuys (1992) stated that the number of tree species in southern Africa is lower than that in the tropics. It was unfortunate that information on other plant forms and species, for example, ferns was not collected as this would have provided an additional element to this biodiversity study.
There are cases, in other scarp forests, of unsustainable harvesting to the point of local extinction of medicinal species such as Ocotea bullata. This species was not encountered in this survey, suggesting it may too have extirpated from this particular forest but has been known to occur in the forest in the past. Medicinal plant exploitation may reduce the species diversity of the forest which would be of concern if there was extensive extraction of all medicinal species.
Regeneration potential using stand density in the Pongola Bush forest fragment is similar to some other studies in Southern African Afromontane patches, but these studies show highly variable results based on a range of biotic and abiotic factors. The Pongola Bush shows a high recruitment of small stems to a decrease in larger stems in a 'reverse J' distribution, as demonstrated in other studies (Adie et al. 2013;Asrat et al. 2020;Eriksson, Teketay & Granström 2003;Mensah et al. 2020). Smaller-stemmed plants such as S. mossambicensis and many smaller individuals of potentially large trees, such as A. falcatus dominate the forest. Populations of C. kraussii have few small stems which contrasts with previous stem diameter studies of forests in the region (Mensah et al. 2020). The limited small stem recruits could potentially affect the future composition of the forest as C. kraussii is a significant contributor to AGB and stem density. Further studies on this species in this location should be included in an integrated monitoring plan. This trend of small stemmed recruits is also occurring for S. chirindensis and K. africana. Combretum erythrophyllum has been recorded in Northern Afrotemperate forests but the occurrence of C. kraussi and K. africana are interesting. Given that these species form a significant part of the forest composition there is likely to be a large shift in dominant species many decades into the future. Species such as A. falcatus and C. africana show the classical 'J' shaped curve which shows that there is a stable regeneration potential for the future of their populations (Botzat, Fischer & Farwig 2015).
These stem diameter differences between species may be due to previous exploitation which has put the forest in earlier stages of succession as seen with other South African forests as well as South American forests of similar composition (Adie et al. 2013;Echeverría et al. 2007).
Emergence of earlier successional features is correspondingly correlated to a decreased patch size or fragmentation (Echeverría et al. 2007), but this forest extent is regulated by various external factors such as historic grassland fire regimes and specific abiotic factors; hence, fragmentation or decreased patch size may not have induced the early successional traits.
Despite there being limited current forest fragmentation of this patch due to its conservation, other patches of like composition may be under threat from fragmentation, particularly from fire, and this patch may have been part of a much larger constellation of related forests in the distant past (e.g. in the Pliocene, when the climate was wetter).
In conclusion, the Pongola Bush is a highly diverse, relatively dense, and stable habitat at present due to current conservation efforts. At the time of the study, four species dominated in terms of biomass, but this was only a preliminary study. A much more detailed study over time is needed in order to contribute to the management and conservation of these types of forests. There remains an information gap on non-tree biodiversity elements within the forest, such as forest floor species, forbs and animal diversity. Having a comprehensive study on the full diversity of the forest will boost the knowledge of the forest and may aid in future conservation. Additionally, the dynamic between forest patches in the region and dispersal mechanisms of both seeds and pollen as genetic material between patches could be better understood. This forest is an understudied remnant of perhaps an early evolutionary stage in the Northern Afrotemperate forest system and its conservation should be encouraged to ensure further studies on the flora and fauna.