Diversity and Structure of an Arid Woodland in Southwest Angola, with Comparison to the Wider Miombo Ecoregion

: Seasonally dry woodlands are the dominant land cover across southern Africa. They are biodiverse, structurally complex, and important for ecosystem service provision. Species composition and structure vary across the region producing a diverse array of woodland types. The woodlands of the Huíla plateau in southwest Angola represent the extreme southwestern extent of the miombo ecoregion and are markedly drier than other woodlands within this ecoregion. They remain understudied, however, compared to woodlands further east in the miombo ecoregion. We aimed to elucidate further the tree diversity found within southwestern Angolan woodlands by conducting a plot-based study in Bicuar National Park, comparing tree species composition and woodland structure with similar plots in Tanzania, Mozambique, and the Democratic Republic of Congo. We found Bicuar National Park had comparatively low tree species diversity, but contained 27 tree species not found in other plots. Plots in Bicuar had low basal area, excepting plots dominated by Baikiaea plurijuga . In a comparison of plots in intact vegetation with areas previously disturbed by shifting-cultivation agriculture, we found species diversity was marginally higher in disturbed plots. Bicuar National Park remains an important woodland refuge in Angola, with an uncommon mosaic of woodland types within a small area. While we highlight wide variation in species composition and woodland structure across the miombo ecoregion, plot-based studies with more dense sampling across the ecoregion are clearly needed to more broadly understand regional variation in vegetation diversity, composition and structure.


Introduction
Tropical woodlands extend over 12 countries in central and southern Africa, with an estimated area of~3.7 million km 2 [1][2][3]. Within this, miombo woodlands are the dominant vegetation type, characterised by trees of the Brachystegia, Julbernardia and Isoberlinia genera, all within the Fabaceae family, subfamily Detaroideae [4][5][6]. These genera are seldom found as dominant species outside miombo woodlands, and while their contribution to the biomass of miombo woodlands is substantial, it varies throughout the region [5]. Across the range of southern African woodlands, variation in climate, edaphic factors, disturbance regimes and biogeography maintain a diverse array of woodland It is important to focus not only on the biodiversity of undisturbed woodland areas but also previously disturbed land in order to properly assess the biodiversity and woodland structure of the Park. Woodland disturbance through shifting cultivation practices produces novel habitats which are not necessarily of lower conservation value [29,30]. Since Bicuar National Park's rejuvenation following the reinforcement of park boundaries after the civil war, many areas of woodland that were previously heavily grazed, farmed via shifting cultivation techniques, and used for timber extraction have been allowed to re-establish and are now protected from further human resource extraction. This presents a unique opportunity to compare the species composition of these disturbed areas with areas of nearby woodland that have not been farmed in living memory.
In this study, we present results of the tree diversity and woodland structure of miombo woodlands found at the far western extent of miombo woodlands in Bicuar National Park, Huíla province, Angola. Our study utilised recently installed biodiversity monitoring plots set up within the Park in 2018 and 2019. We compare the tree diversity and woodland structure of Bicuar National Park with biodiversity monitoring plots previously established in other areas of miombo woodland across the miombo ecoregion which use a common plot biodiversity census methodology. In addition, we take advantage of a unique opportunity to compare the tree species composition of areas of abandoned and now protected farmland that have begun to re-establish as woodland. Specifically, this study aims to: 1.
Describe the tree species diversity and structure of woodlands in Bicuar National Park, and compare this composition with other woodlands across the miombo eco-region.

2.
Explore the role of environmental factors in driving changes in tree species composition across the miombo ecoregion. 3.
Describe variation in tree species composition and woodland structure between disturbed and undisturbed woodland patches within Bicuar National Park.

Study Area
We chose three areas of miombo woodland across the miombo ecoregion to compare with those in Bicuar National Park, Angola (S15.1 • , E14.8 • ). The three sites were Gorongosa District in central Mozambique (S19.0 • , E34.2 • ) [31], Kilwa District in southern Tanzania (S9.0 • , E39.0 • ) [32], and the Mikembo Natural Reserve in Katanga, southern Democratic Republic of Congo (DRC) (S11.5 • , E27.7 • ) [33]. Within each of these woodland sites, multiple one hectare square plots had been installed previously to monitor biodiversity and biomass dynamics. In Katanga, a larger 10 ha plot was subdivided into ten 1 ha plots for this study. We used these previous censuses, collected between 2010 and 2019, to estimate tree biodiversity and woodland structure. Sites range in Mean Annual Precipitation (MAP) from 864 mm y −1 in Bicuar to 1115 mm y −1 in Katanga. Mean Annual Temperature ranges from~20.5 • C in Bicuar and Katanga to~25.8 • C in Kilwa (Figure 1b, Table 1). q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q (a) q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q q showing the plot locations compared to the climate space of the miombo ecoregion estimated using the WorldClim dataset over the miombo woodland vegetation extent with a pixel size of 30 arc seconds (0.86 km 2 at the equator) [34]. Please note that the density colour scale is log-transformed for visual clarity. Bicuar National Park covers an area of~7900 km 2 , established as a hunting reserve in 1938, and later as a national park in 1964 ( Figure 2). While fauna populations in the Park were severely damaged by the Angolan civil war, the interior of the Park remains as a largely intact mosaic of miombo woodland, Baikiaea-Burkea woodland, shrub/thicket vegetation and seasonally flooded grassland. Encroachment of agriculture and grazing, particularly along the northwest and western boundaries of the Park, has led to a fragmented park boundary with patches of diminished thicket and woodland in areas of previously farmed land that have been protected since park boundaries were re-established following the end of the civil war.
Plots in Tanzania were located predominantly within or near the Mtarure Forest Reserve, administrated by the Tanzania Forest Service and protected from human incursion since their installation. Plots were established between 2010 and 2011 in grassy savanna/woodland areas, with plots located along the road network with a 1 km buffer from the road. Plots in Mozambique were established in 2004, in areas of miombo woodland that had been previously used for agriculture but since left fallow, and areas of undisturbed miombo woodland, located along the road network, with all plots >250 m from the road. Plots in DRC were established in 2009 and located within a larger 800 ha miombo woodland reserve, which consists of undisturbed miombo woodlands. All plots were located quasi-randomly, with consideration to accessibility for future woodland censuses.

Plot Data Collection
We sampled 15 one hectare plots in Bicuar National Park and collated data from a total of 64 one hectare plots across the miombo ecoregion within four sites. Figure 1a and Table 1 show the locations and general description of each site, respectively. Plots in Bicuar were situated at least 500 m from the edge of a woodland patch to prevent edge effects which may have altered tree species composition.
Within each plot, every tree stem ≥5 cm stem diameter was recorded, except in the DRC plots, where only stems ≥10 cm stem diameter were recorded. For each tree stem the species and stem diameter were recorded. Tree species were identified using local botanists at each site and taxonomy was later checked against the African Plant Database [35]. At all sites, we used Palgrave [36], along with other texts, to identify tree species. Specimens that could not be identified in the field, or subsequently at herbaria, were described as morphospecies. All tree species within the Bicuar National Park plots were identified. Tree coppicing due to fire, herbivory, and human actions is common in miombo woodlands, therefore, for trees with multiple stems, each stem ≥5 cm stem diameter was recorded, while the parent tree was also recorded for diversity analyses described below.
Stem diameter was recorded at 1.3 m from the ground along the stem (diameter at breast height, DBH) as per convention using a diameter tape measure [37]. Where stem abnormalities were present at 1.3 m from the ground, which precluded the accurate estimation of stem diameter at 1.3 m, the stem diameter was recorded at the nearest 10 cm increment above 1.3 m without significant stem abnormalities [37]. To ensure consistency among stem diameter values recorded at different heights, when the stem diameter was recorded at a height other than 1.3 m the stem diameter at 1.3 m was estimated from the recorded stem diameter using a cubic polynomial equation which adjusts for tree stem taper. This equation was calibrated on 100 stems measured at multiple heights in Niassa Province, Mozambique (Appendix A). Stems below 10 cm stem diameter were not measured in the DRC plots. We therefore estimated the number of 5-10 cm stems in each these plots by extrapolating a linear regression of log stem abundance across the available stem diameter classes.
In addition to the one hectare plots across the miombo ecoregion, we compared the tree biodiversity of undisturbed areas of miombo woodland in Bicuar National Park with areas of disturbed woodland around the edge of the Park that had been previously farmed via shifting cultivation methods, and had since been abandoned and reclaimed within the Park boundaries ( Figure 2). We identified areas previously farmed with the help of park rangers and local residents who identified these areas from memory. We conducted 20 plot surveys of woodland diversity and structure in these areas with 20 × 50 m (0.1 ha) plots, and compared their diversity and structure with 20 × 50 m subsamples of the 15 one hectare plots within the Park interior. Like the one hectare plots, within these smaller 20 × 50 m plots we recorded the species and stem diameter of every tree stem ≥5 cm stem diameter.

Climatic Data
The WorldClim dataset [34] was used to gather data on plot-level climatic conditions. We estimated Mean Annual Precipitation (MAP) as the mean of total annual precipitation values between 1970 and 2000, and Mean Annual Temperature (MAT) as the mean of mean annual temperatures between 1970 and 2000. The seasonality of temperature (MAT SD) was calculated as the standard deviation of monthly temperature per year, respectively. We estimated Climatic Water Deficit (CWD) for each plot according to Chave et al. [40], as the sum of the difference between monthly rainfall and monthly evapotranspiration when the difference is negative, using the dataset available at http://ups-tlse.fr/pantropical_allometry.htm, which uses data from the WorldClim dataset 1970-2000.

Data Analysis
We calculated the basal area of each stem (g i ) using: where d i is the estimated stem diameter of stem i at 1.3 m having accounted for tree taper. We then calculated the total basal area of each plot as the sum of each stem's basal area. For the DRC plots which lacked 5-10 cm stems, we estimated basal area in this stem diameter class from our extrapolation of stem abundance in the 5-10 cm diameter class, assuming a mean stem diameter of 7.5 cm.
All diversity measures were calculated on individual tree-level data, rather than stem-level data, to avoid artificial inflation of abundance for those species which readily coppice. We calculated the alpha diversity of each plot using both the tree species richness of trees with stems ≥5 cm diameter, and the Shannon-Wiener index (H ) (Equation (2)), using the vegan package in R [41]: where S is the total number of species in the plot, p i is the proportional abundance of the ith species and ln is the natural logarithm. We calculated the pairwise beta diversity among sites using the Sørensen coefficient (S S ) (Equation (3)) [42]: where a is the number of species shared between two sites, b is the number of species unique to site 1 and c is the number of species unique to site 2. We calculated S S for each pairwise combination of sites using aggregated species composition data from all plots in each site. The value of S S , which ranges between zero and one, was multiplied by 100 to give a "percentage similarity" between communities in species composition. We estimated abundance evenness for each plot using the Shannon equitability index (E H ) [43] which is the ratio of H to the log transformed species richness.
We analysed the difference in alpha diversity measures and woodland structural variables among groups of plots using Analysis of Variance (ANOVA) statistical models, with a null hypothesis that there was no difference among the mean values of groups of plots. Post-hoc Tukey's HSD tests were used to investigate the degree to which pairwise combinations of plot groups differed in each case.
We used Non-metric Multidimensional Scaling (NMDS) to assess the variation in species composition among one hectare plots, and also between disturbed and undisturbed 20 × 50 m plots within Bicuar National Park, using the vegan R package. The number of dimensions for NMDS was minimised while ensuring the stress value of the NMDS fit was ≤0.1. NMDS analyses were run with 500 random restarts to ensure a global solution was reached. We used Bray-Curtis dissimilarity as the optimal measure of ecological distance [44]. We fit plot-level estimates of MAP, MAT, the seasonality of MAT and CWD to the first two axes of the resulting ordination using the envfit function in the vegan R package to investigate how these environmental factors influenced the grouping of species composition among plots. All analyses were conducted in R version 3.6.1 [45].

Alpha Diversity
In Bicuar National Park we measured a total of 6565 trees within the one hectare plots, and across the four sites, a total of 25,525 trees were sampled. Trees in Bicuar National Park belonged to 48 species within 18 families. Across all four sites we recorded 468 species from 43 families. The most diverse family within each site and among all plots was Fabaceae with 61 species. We encountered 27 tree species in Bicuar National Park which were not found in the other miombo woodland plots ( Table 2). The most common of these unique species were Brachystegia tamarindoides (n = 576), Baikiaea plurijuga (n = 331) and Baphia massaiensis (n = 303). Four species unique to Bicuar National Park within this dataset only had one individual recorded: Elachyptera parvifolia, Entandrophragma spicatum, Oldfieldia dactylophylla, Peltophorum africanum. Alpha diversity in Bicuar National Park was low compared to other sites ( Figure 3). Mean H across plots in Bicuar National Park was 1.6 ± 0.13. An ANOVA showed a significant difference in H among sites (F(3,60) = 7.54, p < 0.01, Table 3), and a post-hoc Tukey's test showed that H in plots in Bicuar National Park was significantly different from those in DRC (H = 2.7 ± 0.19, p < 0.01), Mozambique (H = 2.4 ± 0.2, p < 0.01) and Tanzania (H = 2.2 ± 0.11, p < 0.05). Variation in H is large within Bicuar National Park, with H ranging from 0.85 to 2.56, but this was a similar range to other sites. In contrast, the range of species richness within Bicuar National Park was much lower than other sites, suggesting that the wide range in H was caused by variation in abundance evenness.

Beta Diversity
The NMDS of plot species composition among one hectare plots was run with four dimensions. The stress value was 0.10. Plot diversity in Bicuar National Park formed three distinct groups within axes 1 and 2 of the NMDS ordination. Bicuar plots 9, 13 and 15 were characterised by high abundances of Baikiaea plurijuga, Baphia massaiensis and Croton gratissimus, according to species scores from the NMDS. Bicuar plots 4, 11 and 12 were characterised by Brachystegia tamarindoides and Ochna pulchra. The third group consisting of the remaining seven plots surprisingly had a species composition most similar to that of plots in the DRC group according to the NMDS, sharing the core miombo species of Julbernardia paniculata and Pterocarpus angolensis. This group of plots in Bicuar National Park was further characterised by the abundance of Pterocarpus lucens, Strychnos pungens and Bridelia mollis however, which were not present in the DRC plots. All environmental factors fitted to the NMDS ordination correlated significantly with the grouping of plots (Figure 4a). MAT explained the most variation in plot position on the first two NMDS axes (R 2 = 0.75, p < 0.01), followed by CWD (R 2 = 0.54, p < 0.01), the seasonality of MAT (R 2 = 0.46, p < 0.01) and MAP (R 2 = 0.4, p < 0.01). Variation in MAP explained much of the difference among plots in Bicuar National Park versus those in Tanzania    The pairwise Sørensen coefficient of percentage similarity (S S ) showed that the species composition of plots in Bicuar National Park had low similarity with other sites in the study, sharing a few species with other sites (Table 4). Similar to the NMDS, these results show that plots in Bicuar National Park are most similar to those found in DRC.
The stem diameter abundance distribution in Bicuar National Park was comparable with other sites (

Effect of Disturbance via Shifting Cultivation on Diversity within Bicuar National Park
There was a clear difference in the species composition of previously farmed disturbed woodland plots and undisturbed woodland plots, but with some overlap (Figure 6). Notably, Plots 4 and 7 in putatively undisturbed woodland have a species composition more resembling the disturbed plots.

Comparison of Bicuar National Park with Other Woodlands across the Miombo Ecoregion
We compared the tree species diversity and woodland structure of arid woodlands in Bicuar National Park in southwest Angola with three other woodland sites across the miombo ecoregion. Our results show that Bicuar National Park is distinct in both woodland structure and species composition from these other woodlands. Notably, plots in Bicuar National Park contained 27 tree species which did not occur at other sites. This lends support for the Huíla Plateau as an important area for conservation of southern African woodland landscapes. The woodlands in Bicuar National Park were of low tree basal area, with few large trees except in plots dominated by Baikiaea plurijuga. Many other studies have drawn a relationship between water availability and basal area [46,47], and our study supports this, with Bicuar National Park being the most arid of the four sites considered in our study. The NMDS of species composition also suggests that plots in Bicuar National Park are influenced by aridity. While there are more arid woodlands within southern Africa, with Mopane woodlands for example often being particularly dry, these plots in Bicuar National park represent particularly dry miombo woodlands.

Delineation of Woodland Types within Bicuar National Park
Within Bicuar National Park, three distinct woodland types were identified. The first, dominated by Baikiaea plurijuga and Baphia massaiensis represents the Baikiaea woodland type commonly found to the south of the miombo ecoregion [48]. This is supported by Chisingui et al. [23] who also found Baikiaea woodlands as a distinct woodland type in the Park. B. plurijuga has been identified as an important species for conservation, being attractive for selective logging due to its large stature [49,50]. The woodlands created by B. plurijuga are also an important habitat for elephants (Loxodonta africana) [51,52], with Bicuar National Park and Mupa National Park being key refugia for this animal in the Huíla plateau region. The second woodland type, dominated by Brachystegia tamarindoides and Ochna pulchra represents a form of small stature woodland with a shrubby understorey and sparse canopy trees, which commonly occurs as a result of repeated disturbance by fire, or poor soil structure [53]. The remaining plots resemble the more archetypical miombo woodland with Julbernardia paniculata, though with several species not seen in plots further to the east in the miombo ecoregion such as Strychnos pungens. This mosaic of woodland types makes Bicuar National Park a valuable reservoir of diversity and strengthens the case for the Park being a key conservation asset within the Huíla plateau and the larger southern African region. While there are regional boundaries between Baikiaea and miombo woodlands [1], within Bicuar National Park it is likely that the mosaic of woodland types has been created by a combination of soil water capacity and disturbance history. Bicuar has a distinct landscape of wide shallow grassy valleys surrounded by woodland on higher ground ( Figure 2). On some of these high points the soil is particularly sandy, resembling the Kalahari sand soils found further east and south [20], and these areas coincide with the presence of Baikiaea woodlands [5]. High levels of disturbance by fire in these Baikiaea patches may additionally prevent a transition to an alternative woodland type via the control of sapling growth.

Comparison of Disturbed and Undisturbed Woodland Plots
Previously disturbed woodlands around the edge of Bicuar National Park were found to share many species with undisturbed plots in the Park, but with some additional species which did not occur in the undisturbed plots. They also lacked notable archetypical miombo species which tend to form larger canopy trees such as Brachystegia spiciformis and contained very few Julbernardia paniculata, leading to a distinct woodland composition. The species diversity of these disturbed patches was higher on average than was found in the undisturbed plots, a result which has been corroborated by other studies in miombo woodlands [54][55][56]. Other studies have shown a peak in species richness during woodland regrowth as pioneer species take advantage of a low competition environment, while some later stage woodland species remain as residuals that survived the original disturbance [30,57]. Gonçalves et al. [30] particularly, notes the dominance of Pericopsis angolensis and Combretum spp. as light-demanding pioneer species, which were found to be abundant in the disturbed plots here. This suggests that reclamation of previously farmed and abandoned land for landscape conservation in this ecological context is a valuable management strategy.
In disturbed plots near the edge of the Park, there was a lack of species which tend to grow to large canopy trees, possibly due to them being repeatedly felled for timber prior to reclamation by the Park, or due to them being unable to recruit into a more open, shrubby woodland. Despite this lack of canopy forming tree species, some disturbed plots had a greater basal area than undisturbed plots, possibly due to high levels of coppicing in these plots or a divergent fire history. Indeed, mean stem density was higher in undisturbed plots. This can lead to species that would otherwise remain small producing a much larger basal area as they grow multiple stems under high disturbance conditions [58]. The most common species in the disturbed plots were Combretum psidioides, Combretum collinum and Terminalia sericea, members of the Combretaceae family, all of which more commonly remain as smaller multi-stemmed trees in disturbed woodlands, rather than growing to larger canopy trees [59]. This result could be considered at odds with other studies which report lower woody biomass in plots that have experienced harvesting (e.g., Muvengwi et al. [60]). It is important to consider however that our study took place in plots that were measured after farming had been abandoned for at least 7 years, with time for regeneration to occur. It is possible that over time tree basal area will decrease as coppiced shrubby trees are replaced by core miombo species in the transition back to miombo woodland [30]. Indeed, other studies in miombo woodlands across the ecoregion have reported substantial recovery within seven years, with high levels of biomass accumulation in previously disturbed plots [30,61]. Bicuar National Park offers a valuable case study to track woodland regeneration in real time over the next decade in these previously farmed and now protected woodland plots, which could improve our understanding of this potential post-disturbance peak in basal area.
In conclusion, the woodlands of Bicuar National Park represent an important woodland refuge at the far western extent of the miombo ecoregion. These woodlands, both those disturbed by previous farming activity and those which remain undisturbed, possess several species not found commonly in other miombo woodland plots around the region. They may also house important genetic variation for widespread species, representing populations adapted to more arid conditions. Our study highlights the variation in species composition across the miombo ecoregion and underlines the need for studies which incorporate plot data from multiple locations to reach generalisable conclusions about the region as a whole. Additionally, the installation of 15 one hectare woodland monitoring plots and a further twenty 20 × 50 m plots in previously farmed and now protected land offer a valuable natural laboratory to further explore the dynamics of dry miombo woodlands of the Huíla plateau. Bicuar National Park should be considered a key conservation asset within the Huíla plateau and within the miombo ecoregion as a whole, as a successfully protected example of an arid woodland mosaic.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Abbreviations
The following abbreviations are used in this manuscript: