The structure and composition of the woody plant communities of Majete Wildlife Reserve , Malawi

Results: Five distinct woody plant communities, two of which were subdivided into three sub-communities each, were recognised in MWR, and a total of 118 woody plant species within 31 families were identified. A description of the location, structure and species composition of each community is provided. Miombo was the most widespread community (covering 35.9% of the area), while the lower-altitude shrublands and woodlands were the richest floristically.


Introduction
Savanna ecosystems cover approximately 40% of Africa (Scholes & Walker 1993), and the most common variant in southern and eastern Africa, namely miombo woodland, covers approximately 2.7 million km 2 in seven countries (Frost 1996;Kutsch et al. 2011). Malawi falls entirely within the Miombo Ecoregion (~3.6 million km 2 ; Byers 2001), and the ecoregion was identified as one of five global wilderness areas that require conservation prioritisation (Mittermeier et al. 2003). The Miombo Ecoregion is an amalgamation of several smaller ecoregions, with the southern section of Malawi containing the Eastern Miombo Woodlands Ecoregion (EMWE), extending from southeastern Tanzania and northern Mozambique into parts of southeastern Malawi, and the Zambezian and Mopane Woodland Ecoregion (ZMWE) in southwestern Malawi (Malaisse & Parent 1985;Byers 2001). The EMWE is dominated by miombo sensu stricto, characterised by trees in the genus Brachystegia, occurring either alone or together with Julbernardia and Isoberlinia species (Geldenhuys & Golding 2008). The EMWE is further subdivided into drier miombo and wetter miombo by the 1 000 mm rainfall isohyet (White 1983;Byers 2001). The ZMWE consists mainly of broadleaf, deciduous savannas and woodlands interspersed with edaphic grasslands (White 1983;Chirwa, Syampungani & Geldenhuys 2014). Collectively, these two ecoregions were formerly classified by Huntley (1982) as Caesalpinoid (now Detarioideae) woodland, or broadleaved dystrophic savanna woodland, while White (1983) classified the area as the Zambezian Phytoregion, with Zambezian woodland being the most widespread woody plant community in the region.
True miombo woodland constitutes 92.4% of the natural vegetation of Malawi (Government of Malawi 2010). The vast majority of extant natural vegetation is, however, confined to small, isolated areas that are subjected to rapid deforestation rates of up to 0.9% annually (FAO 2015). Additionally, while projections of future drier climates in sub-Saharan Africa suggests the range of dry miombo may expand (Jinga & Palagi 2020), a sharp increase in the human population (Eastwood & Lipton 2011) will place undue pressure on remaining miombo woodlands, which already support the livelihoods of more than 100 million people either directly or indirectly (Syampungani et al. 2009). The role of protected areas (PAs) as sanctuaries for miombo vegetation, as well as other woodland vegetation types, will thus become even more significant in the future, and knowledge of their composition and structure is needed as a basis for their sustainable management. Majete Wildlife Reserve (MWR) in Malawi was proclaimed in 1955 and was taken over in 2003 by African Parks with the aim of restoring, developing and preserving the integrity of all facets of the area's biodiversity. The reserve is situated in southwestern Malawi, at the intersection of the EMWE and the ZMWE (Byers 2001;Munishi, Temu & Soka 2011), making this a unique study site to compare contrasting vegetative communities. Different plant communities may also require specific management approaches, adding to the need for an accurate vegetation description at a scale appropriate to management. Though useful, the coarse-scale classification of ecoregions fails to describe plant communities and their distributions at the scale of PAs. Classifying vegetative communities at scales finer that those of ecoregions can be extremely resource intensive and is therefore likely not realistic in most African countries. More cost-effective solutions are thus needed, and remote sensing may offer such a solution to delineate woody plant communities in Africa (Peel, Kruger & MacFadyen 2007;van Rooyen et al. 2008). The delineation of non-woody vegetative components, which may be crucial to certain aspects of environmental management, is however likely not achievable with remote sensing approaches.
The aims of this study were to (1) classify, describe and map the woody plant communities of MWR using a combination of remote sensing and on-the-ground surveys, and (2) to compile an inventory of the tree and shrub species present in MWR. This information is intended to provide a basis for improved management planning and policy development, including fire management, the placement of infrastructure, and the re-introduction of extirpated mammal species, as well as providing a baseline against which to monitor change (Kremen 1992;Peel et al. 2007).

Study site
Majete Wildlife Reserve (MWR) is situated in the Lower Shire Valley, part of the Great Rift Valley in the south of Malawi. The reserve covers roughly 700 km 2 and is situated in a tropical climate zone. The underlying climate can be categorised into three distinct seasons (Hall-Martin 1972), namely the wet season (December to March), the cool dry season (April to August), and the hot dry season (September to November). Mean monthly temperatures range from 20.5°C (June and July) to 27.8°C (November), with maximum temperatures often reaching ~45°C. A distinct N-S line divides MWR into the western upland area (altitude > 350 mamsl) and the eastern lowland area (altitude < 350 mamsl). The topography of the eastern lowlands is relatively flat with few rocky outcrops, with mean annual precipitation (MAP) ranging between 680-800 mm, while the western uplands consists of many rocky outcrops, and receiving slightly higher MAP (700-1 000 mm) (Hall-Martin 1972). A detailed description of the geology and soils of MWR is lacking, but most soils are shallow, stony lithosols of poor nutrient status, with fertile alluvial soil only occurring along a few river beds (Bell 1984;Sherry 1989). Two perennial rivers occur in the northeastern section of MWR, which, along with ten artificial water points and seven natural springs, provide the only source of water for animals in the dry season. The vegetation was previously described by Sherry (1989) as tropical dry deciduous woodland, miombo savanna woodland, riverine associations, and riparian thicket. Several browsing herbivore species that affect the structure of woody vegetation are present in MWR, notably elephant (Loxodonta africana), black rhinoceros (Diceros bicornis), giraffe (Giraffa camelopardalis), kudu (Tragelaphus strepsiceros) and impala (Aepyceros melampus).

Delineation of land cover types with remote sensing
The Sentinel-2 cloudless land cover product was obtained from the Copernicus global land service (Buchhorn et al. 2019), and used to map land cover at 100 m spatial resolution for MWR in QGIS v3.10. The product uses supervised classification of reflectance data and provides 23 distinct land cover classes globally. Five of these (open woodlands, closed woodlands, shrublands, herbaceous vegetation, and permanent inland water bodies) were present in MWR. The classification product was in sinusoidal UTM projected at Datum WGS84 and is updated at annual time steps. The latest available iteration (2018) was used. Additionally, the ASTER digital elevation model (DEM) version 3 (ASTGTM v003) was obtained from the U.S. Geological Survey and used to create an elevation profile for MWR at a spatial resolution of 1 arc second (~ 30 m). The inclusion of a DEM would likely allow for better separation of plant communities (Sedano, Gong & Ferrão 2005). The DEM was simplified in QGIS v3.10 to display three equal interval elevation categories in MWR, and used to divide land cover classes into separate elevation categories where relevant. The lowest point in MWR is where the Mwembezi River exits the reserve (90 mamsl) in the northeast, and the highest point is at Kapirimbewe in the northwest (880 mamsl). The final land cover classes were reprojected to 5 ha grid cells for subsequent analyses.

Delineation of plant communities
Following the delineation of land cover classes with remote sensing, 195 field sites (circular plots of 50 m radius, covering an area that was considered adequate to encapsulate the composition of the woody plant community) were selected across all classes for the final classification of woody plant communities based on the dominant tree species. Each site was chosen to be homogenous and representative of the respective land cover classes, and the sites were surveyed during May-June 2020. The number of sites surveyed in each land cover class was proportional to its extent (area) and ranged from two to 66 per land cover class. At each of the 195 field sites, a count of the number of individual trees or shrubs belonging to eight woody genera was made. Based on their prevalence in MWR, and woody plant community classifications in other savannas (White 1983;Byers 2001;Peel et al. 2007;Gandiwa & Kativu 2009;Munishi et al. 2011;Clegg & Connor 2012;Maquia et al. 2019), the chosen genera incorporated the most important characteristic species for classifying woody plant communities, and were ultimately used as a basis for describing the final plant communities present in MWR. The chosen genera were: Senegalia (including S. galpinii, S. nigrescens and S. burkei), Vachellia (including V. karroo, V. nilotica and V. torrei), Brachystegia (including B. boehmii, B. allenii, B. longifolia, B. spiciformis and B.utilis), Julbernardia (including J. globiflora), Sclerocarya (including S. birrea), Terminalia (including T. mollis, T. sericea and T. stenostachya), Combretum (including C. adenogonium, C. apiculatum, C. hereroense, C. imberbe, C. molle, C. mossambicense and C. zeyheri), and Colophospermum (including C. mopane).

Accuracy assessment
Woody vegetation structural data were collected at each of the 195 field sites to assess whether the remotely sensed land cover classes were matched by vegetation physiognomy on the ground. At each field site, the physiognomy and structure of the vegetation was classed according to the broad classification of Edwards (1983) by recording the dominant primary growth form (trees or shrubs), projected crown cover classes (closed, open, sparse or scattered), and the mean number of crown diameters by which plant crowns were separated. The accuracy assessment was done by comparing the classified pixel to the same site in the field through a confusion matrix that compared the Edwards classifications to the Sentinel land cover classes. Overall map accuracy was determined by summing the number of pixels classified correctly and dividing by the number of pixels (Mohd Hasmadi, Pakhriazad & Shahrin 2009).

Plant community richness and diversity
At 94 of the 195 field sites, a list of all tree and shrub species present was compiled by noting the species that were visible to an observer walking along a line transect of 50 m. The location of the transects was chosen to be relatively evenly distributed across MWR, while still encapsulating obvious landscape variations (such as topography and riverbeds) and to ensure that all woody vegetative variations were sampled. Prior to conducting the inventories, the species-area method (Smith, Meredith & Johns 1999) was used to determine the optimal length of transects by noting the increase in the number of species observed as the transect length increased. Transect length was deemed optimal when additional sampling did not identify a significant number of novel species. Species were identified on site where possible, but a sample was taken for later identification where uncertainty existed.
Sampling performance of species inventories made along the 94 transects was evaluated through the construction of species accumulation curves for each woody plant community and for MWR as a whole (Colwell & Coddington 1994). Species accumulation curves also provided a reliable estimate of species richness due to the variation in sample sizes, and provided an estimate of the predicted number of species that could be found if sampling continued indefinitely based on the rate of species discovery in this study (Colwell 2005). For this purpose, the performance of six non-parametric species richness estimators appropriate for incidence-based data were used, namely the Chao 2 estimator, the bias-corrected form of the Chao 2 estimator (Chao 2-bc), a newer version of the Chao 2 estimator (iChao2), the incidence-based coverage estimator (ICE), and two Jackknife species richness estimators (Jack 1 and Jack 2) (Burnham & Overton 1978;Chao 2005;Chiu et al. 2014). The sample order was randomised 100 times for each of the calculations to provide mean statistics at each sample order, thereby generating smooth accumulation curves (Colwell & Coddington 1994;Nieman, Leslie & Wilkinson 2019). Due to its relatively small extent in MWR, grasslands (Gr) were grouped with shrublands and woodlands (SW -see plant community descriptions) based on the clear similarity in geographical location and woody species similarity between these plant communities. In addition, the woody species diversity and evenness or equitability was quantified with the use of standard ecological diversity indices for MWR as a whole, as well as for each of the plant communities. Specifically, woody species diversity was expressed using the Shannon-Wiener diversity index (H'), the Simpson diversity index (1-D), and inverse Simpson diversity index (1/D). In all instances, higher values equate to a higher level of diversity in the sample, and thus higher entropy and less dominance by one or a few species. Evenness was expressed with the Shannon evenness index (J'), wherein values approaching zero represent low evenness in the sample, and thus dominance by one or a few species. The coefficient of variation (CV) was used to describe the degree of heterogeneity among species incidence, where zero represents a completely homogenous sample. Finally, the Jaccard coefficient was used to determine similarity between the species' assemblages compiled for each plant community, and a dendrogram was constructed using the k-means clustering algorithm with Euclidian distances. All species richness estimator and index values were calculated in R v4.0.1.

Woody plant communities recognised
The final classification of the land cover of MWR was done at 90.77% accuracy (i.e. 90.77% of pixels from the Sentinel map matched the corresponding woody plant community structure on the ground), and included five woody plant communities, two of which were further subdivided into three sub-communities each ( Table 1). The extent and distribution of these plant communities and sub-communities were subsequently mapped at the scale of the individual 5 ha pixels ( Figure  3). At the level of the individual plant communities and sub-communities, a clear distinction in terms of woody species assemblages was confirmed by the Jaccard similarity index ( Figure 1). The floristic composition of riparian woodland was the least similar to all other woody plant communities, while grassland was closely related to the shrubland and woodland sub-types. The floristic composition of all miombo sub-types were also closely related.

Riparian woodland
Riparian woodland (Rw) (Figure 2A) occurred on alluvium fringing the only two perennial rivers in MWR, namely the Mkulumadzi and Shire rivers. The rivers respectively form 9 km and 7 km of the outer reserve boundary in the northeastern section of the reserve. As expected (Smith, Meredith & Johns 1999), this plant community was characterised by a distinct species composition and density compared to the rest of MWR, but unlike the other plant communities, was not characterised by one or a few dominant genera or species. The species recorded most frequently (≥ 40%) were Philenoptera violacea, Allophylus africanus, Senegalia nigrescens, Terminalia sericea, Sclerocarya birrea, Cleistochlamys kirkii, Kigelia africana, Combretum imberbe, Pterocarpus rotundifolius, Cremaspora triflora, Gymnosporia senegalensis, Diospyros squarrosa, Croton megalobotrys and Grewia flavescens (Table 2). This woody plant community was made up by 52 pixels, thus covering 245.1 ha, or 0.34% of MWR. Tree cover was typically between 10 and 75%, and tree height often exceeded 15 m. The woody plant community occurred at low elevation (< 250 mamsl), and contained very little grass cover, while bare ground patches and forbs were prominent.

Grassland
Grassland (Gr) ( Figure 2B) was found in a few scattered locations in the low-altitude (< 250 mamsl) belt of the MWR, and its presence was thus not significant, as it only covered 221.6 ha, or 0.31% of MWR. The factors determining the occurrence and distribution of grasslands are not known, but likely pertain to edaphic factors. Additionally, the location of most grassland sites on the boundary of MWR suggests that these areas may have been illegally logged and cultivated by local communities prior to the erection of the reserve fence in 2004. Grasslands in MWR may thus represent a state of secondary vegetation succession, and it is possible that they may eventually revert back to shrubland and woodland states if left undisturbed. The dominant tree and shrub species were similar to those reported in shrublands and woodlands (SW), but woody species were mostly absent or their growth was stunted. Both tree and shrub cover were less than 0.1%, while grass height mainly exceeded 1 m, and grass cover was between 10 and 100%. Additionally, grassland sites were often accompanied by large expanses of bare ground, and are thus presumably exposed to higher erosion risk due to decreased soil stabilisation.

Shrublands and woodlands
Shrublands and woodlands (SW) ( Figure 2C-D) covered the largest portion of the low-altitudinal (< 250 mamsl) band in the east of MWR, and extended into some of the medium-altitude (250-450 mamsl) areas. This woody plant community was floristically rich and characterised by species in the genera Combretum, Terminalia, Vachellia, Senegalia and Sclerocarya. The most frequently recorded species (≥ 40%) were Combretum adenogonium, Senegalia nigrescens, Combretum imberbe, Diospyros mespiliformis, Terminalia sericea and Sclerocarya birrea (Table 2). Miombo dominants such as Brachystegia and Julbernardia are extremely gregarious, and rarely occur in other woody plant communities (Geldenhuys & Golding 2008). The shrublands and woodlands in MWR were thus distinguished by the absence of typical miombo species rather than by the presence of other species (White 1983). The woody plant community was further subdivided into three sub-communities on the basis of tree and shrub cover and height that were intermixed with each other, namely shrubland (SW1), open woodland (SW2) and closed woodland (SW3).
Open woodland (SW2) was the most widespread and prevalent of the three sub-communities, covering 15 207.4 ha, or 21.30%, of MWR. The sub-community was differentiated by its open canopy, with tree cover ranging between 1 and 10%. Closed woodland (SW3) was much less widespread, and confined to small, isolated patches covering 2 055.3 ha, or 2.88% of MWR. The sub-community was differentiated by its closed canopy, with tree cover ranging between 10 and 75%.
In both SW2 and SW3, tree height was generally greater than 10 m, and shrubs higher than 1 m in height covered less than 10% of the total area. Shrubland (SW1) was widely dispersed, but covered only 2 592.7 ha, or 3.63% of MWR. It was dominated by shrubs of 2-5 m in height, and stunted trees. Total tree cover was below 1%, while total shrub cover ranged between 10 and 100%. It is unclear why the woody component in these areas is stunted, although either short fire return intervals, high fire intensity or high herbivory pressure (particularly high elephant densities), or a combination of these, may have played a role (Trollope, 1998;Shannon et al. 2011).

Transitional woodland
Transitional woodland (Tw) ( Figure 2E) was the second most widespread plant community in MWR, and existed as a transitional state between shrublands and woodlands (SW), and miombo (M) at lower and higher elevations respectively. The dominant genera thus comprised of a mixture of the genera found in the neighbouring plant communities, as well as interspersed patches of Colophospermum mopane, and included predominantly Combretum, Terminalia, Senegalia, Colophospermum and Brachystegia. The most frequently recorded species (≥ 40%) were Combretum adenogonium, Diospyros mespiliformis, Diplorhynchus condylocarpon, Combretum hereroense, Colophospermum mopane and Terminalia sericea, while miombo species (e.g. Brachystegia longifolia and B. utilis) were also commonly recorded ( Table 2). The area occupied 5 392 pixels, thus covering 25 417.9 ha, or 35.59% of MWR. The plant community occurred at medium-altitude (250-450 mamsl). Tree cover ranged from 1 to 10%, and trees were generally higher than 10 m. Shrubs higher than 1 m covered less than 10% of the total area. The isolated pockets of mopane woodlands were dominated by small (< 5 m) Colophospermum mopane shrubs, either alone or together with Senegalia nigrescens, Vachellia nilotica, Cassia abbreviata, Dalbergia melanoxylon, Pterocarpus lucens, Pterocarpus rotundifolius and Xeroderris stuhlmannii. The extent of these patches was however too small to be included as a significant and separate woody plant community in the MWR.

Miombo
Miombo woodlands (M) ( Figure 2F) covered the largest area of MWR (35.9%), and were located at high altitudes (> 250 mamsl, but mostly > 450 mamsl) in the western upland escarpment, where steep slopes and high rock cover predominated. This thus substantiated the notion that slope and elevation are important determinants of miombo plant community patterns (Munishi, Temu & Soka 2011). As expected of typical miombo ecosystems, the miombo plant community was dominated by trees in the legume sub-family Detarioideae (formerly Caesalpinoideae), and the genera Brachystegia and Julbernardia. The widespread presence of Brachystegia boehmii, Brachystegia spiciformis, and to a lesser extent, Julbernardia globiflora (Table 2), coupled with typical MAP < 1 000 mm, placed the MWR's miombo woody plant communities in the dry miombo category (White 1983;Frost 1996;Kapinga et al. 2018;Moura et al. 2018). Additionally, characteristic species of wet miombo (e.g. Brachystegia floribunda, B. glaberrima, B. taxifolia and B. wangermeeana) (White 1983) were absent from MWR. Other commonly occurring species found in miombo in this study were Pseudolachnostylis maprouneifolia and Diplorhynchus condylocarpon, as previously found by Ribeiro et al. (2008) in northern Mozambique. A well-developed grass layer was also present throughout. The miombo woodland community was further subdivided into three sub-communities based on altitude and canopy cover, namely open miombo (M1), low-altitude, closed miombo (M2), and high-altitude, closed miombo (M3). In all three sub-communities, trees were generally taller than 10 m, and shrubs above 1 m in height covered less than 10% of the total area.

Woody species diversity
A total of 868 individuals of 118 woody plant species were identified along 94 transects in MWR, representing 16 orders and 31 families (Supplementary Material). The most species belonged to the families Fabaceae (37 species), Combretaceae (10 species), Phyllanthaceae (8 species) and Malvaceae (7 species). Together, these four families contributed more than half (52.5%) of all species recorded. Most species were recorded infrequently, with 29 species recorded only once, and a further 17 species only recorded twice.
The species accumulation curves for MWR (Figure 4) clearly indicated that the expected number of species had reached an asymptote, suggesting that sampling size was sufficient, and that an adequate proportion of the total woody species present in MWR were identified in our survey (Begossi 1996). In contrast, none of the individual woody plant communities ever reached an asymptote. However, given the rapid initial rates of discovery, the relatively large study area, and the high species richness present in all woody plant communities, it is unlikely that a true asymptote will ever be reached without excessive sampling (Williams, Witkowski & Balkwill 2007). riparian-, miombo-and transitional woodlands had comparatively lower diversity. Shannon-Wiener diversity (H') values ranged between 3.80 ± 0.50 (riparian woodland) and 4.02 ± 0.05 (shrublands and woodlands), while Simpson diversity values (1-λ) all ranged between 0.96 and 0.98, with the exception of riparian woodland (0.83). Evenness values were highest for miombo and riparian woodland (J' = 2.40 and 2.34, respectively), and comparatively lower for transitional woodlands, and shrubland and woodlands (J' = 2.15 and 2.12, respectively). Overall evenness in all samples was low (J' = 2.04). The species richness estimates and selected ecological diversity index values can be found in Table 3.

Discussion The plant communities of Majete Wildlife Reserve
Five distinct woody plant communities, along with six sub-communities, were identified and classified in this study (Table 1; Figure 3) based on their altitudinal and geographical location, vegetation physiognomy and structure, and woody plant species composition. Geographically, all of the woody plant communities can be broadly placed in the Miombo Ecoregion (Byers 2001), while floristically, all would fall within the Zambezian Phytoregion (White 1983;Byers 2001;Goyder et al. 2018). At the level of the entire MWR, it is clear from the species richness estimators that further sampling would not yield a significant number of novel, unrecorded tree or shrub species. However, at the level of individual plant communities, it was evident that further sampling could result in a far more complete species list. Similarly, the species richness estimators predicted that between 76 and 86% of all woody species in the MWR were recorded in our survey, which is regarded as adequate (Heck Jr, van Belle & Simberloff 1975). However, at the individual plant community level, it was predicted that up to 75 more woody species in addition to the 44 recorded, could be found in miombo alone. This is likely due to the comparatively low sampling intensity in miombo, as a result of difficulty in accessing the area. Mwase et al. (2007) and Missanjo et al. (2014) found 48 and 22 tree and shrub species in Malawian miombo ecosystems respectively, on plots ≤ 6.5 ha in size. Further sampling of the miombo community in MWR will thus likely be needed to provide a complete list of species present in MWR. Similarly, the riparian woodlands were slightly underrepresented in this study, and many more species are predicted to occur in this woody plant community. However, given the limited extent of riparian woodlands in MWR, it is unlikely that more sampling could occur without introducing spatial autocorrelation biases, and further sampling would thus have to occur outside the reserve perimeter. The reserve is also dominated by a vast network of seasonal streambeds that become active in the rainy season (December-March), and it is possible that sampling of these streambeds could reveal similar floristic compositions to those found along perennial rivers in this study. However, despite the apparent incompleteness of the species inventory, the genera used to differentiate plant communities (Combretum, Senegalia, Terminalia, Sclerocarya, Brachystegia, Julbernardia, Colophospermum and Vachellia) were all easily and frequently recorded. The failure to detect a potentially large number of species does therefore not affect the plant community classifications produced in this study. In fact, with the exception of the genus Julbernardia, and the addition of the genera Diplorhynchus and Diospyros, these were the most abundant genera recorded in MWR. All quadrats were thus dominated by a particular combination of these genera.

Remote sensing as an aid to vegetation surveys
Producing vegetation maps for large areas from groundbased measurements is extremely resource-intensive, and their accuracy would likely be compromised in the absence of intensive surveys. Remote sensing provides a relatively inexpensive and efficient way to overcome this logistical barrier through relatively fine-scale land cover classifications over a large area. In this study, remote sensing also allowed for the accurate placement of sampling plots that were used to add detail for the completion of the vegetation survey. Despite the overall congruence between remotely sensed and groundsurveyed information in this study, as determined by the accuracy assessment, remote sensing should be accompanied by ground-truthing methods in future studies, and in many instances (especially for fine-scale vegetation mapping), more than simple ground-truthing may be necessary (Smith, Meredith & Johns 1999). Remote sensing on its own also lacks the fine detail obtained from ground-based surveys, further substantiating the need for integrating both approaches to achieve optimal vegetation classifications (Kerr & Ostrovsky 2003;Clegg & Connor 2012). The combination of remote sensing and ground surveys was previously shown to be a successful method for delineating plant communities in southern Africa (Peel et al. 2007;van Rooyen et al. 2008) and elsewhere (Smith, Meredith & Johns 1999;Satyanarayana et al. 2011).
Because detailed information on the distribution and extent of ecosystems, and contractions or expansions of land-types, will become increasingly important in the face of global change (Jinga & Palagi 2020), the importance of rapid classification tools such as remote sensing will likely increase. Additionally, the relatively low cost of woody vegetation classifications through remote sensing, even when combined with field surveys, offers an attractive application for protected and near-natural areas in developing countries, as well as in regions that undergo rapid and constant change. Future management and conservation initiatives will thus likely rely on similar methods as those used in this study to identify and describe habitats or species that require prioritised attention. Some limitations of the method can however be expected in heavily degraded or non-natural areas, such as where natural areas have been converted for agricultural use. In these areas, poor resolution and misclassifications of remote sensing products may have a large impact on vegetation mapping and classification, and land cover classes alone will likely not be sufficient.

Implications for fire management
Most savanna trees are morphologically and physiologically adapted to persist in the presence of regular fires (Duvane et al. 2017), but fire effects may vary substantially depending on the frequency, intensity and season of the burn, as well as on woodland type (Smit et al. 2010) and interactions with other processes such as herbivory (Shannon et al. 2011). Woody vegetation maps thus provide a useful tool to plan fire management according to the distribution of dominant tree species. However, elements of fire regimes are also closely connected to non-woody components of vegetation, and information on non-woody vegetation structure and composition can thus improve approaches to fire management significantly (Nieman 2020).
In this study, the two dominant and contrasting woody plant communities, namely shrublands and woodlands, and miombo woodlands, differed substantially in terms of species composition and density, and the approaches to fire management in these areas may have to be adjusted accordingly (Nieman et al. 2021). The dominant genera of shrublands and woodlands are typically more tolerant of repeated fires of high-intensity than their miombo counterparts. In particular, Vachellia species display enhanced coppicing abilities, and fire is thus rarely able to prevent establishment and recovery, even after substantial damage to tree canopies and large-scale die-back are caused (Russell, Tedder & Demmer 2019), and even in drier systems where the effects of fire are more pronounced (van der Merwe et al. 2019). Similarly, fire has been found not to affect germination rates of Senegalia species (Walters, Midgley & Somers 2004), and Sclerocarya birrea was shown to be highly tolerant to fire (Luoga et al. 2004). For this reason, reversing undesirable woody encroachment in shrublands and woodlands would require the repeated application of high-intensity fires (Walters et al. 2004), which has been shown to be effective in reducing tree and shrub densities in Combretum, Senegalia and Sclerocarya woodlands (Sweet 1982;Enslin et al. 2000), as well as reduce tree heights and basal areas (Trapnell 1959;Shackleton & Scholes 2000;Nefabas & Gambiza 2007;Gandiwa & Kativu 2009), and hamper seedling establishment and growth (Jacobs & Biggs 2001). Ultimately, this may lead to the formation of extensively coppiced shrublands (or in extreme circumstances over a prolonged period, grasslands) (Mapaure 2001), by continually preventing trees from growing out of the fire trap and developing into tall adults (Bond & van Wilgen 1996), as seen in nearly 3 000 ha of MWR.
Increased fire frequency has also been shown to promote tree and shrub density in mopane woodlands (Gandiwa & Kativu 2009), and to increase the number of stems on individual plants (Jacobs & Biggs 2001;Kennedy & Potgieter 2001). In areas of MWR where the goal is to increase the incidence of taller, closed vegetation, fire should thus be excluded for a sufficient period, or applied as low-intensity burns in the wet or early dry season (Enslin et al. 2000;Kennedy & Potgieter 2001). Finally, Brachystegia-Julbernadia tree species have been extensively shown to be unable to adequately recover following annual or even biennial burns, and particularly from high-intensity fires in the late dry season (Thomson 1975;Ryan & Williams 2011), resulting in the conversion of miombo woodland to a grassland state (Furley et al. 2008). In contrast, the complete absence of fire may allow miombo woodland to develop into a closed canopy forest (Trapnell 1959). Fires in the miombo woodlands of MWR should therefore be limited to the early dry season, and on a rotation of at least two years (Nieman et al. 2021).
Additionally, riparian woodlands are arguably the most important woody plant community for ensuring the survival of large mammals due to their proximity to water, and provision of shade and highly palatable vegetation (Anderson & Walker 1974;Conybeare 2004), making appropriate conservation essential. Riparian woodlands do not burn frequently, and tend to persist in the absence of fire (van Wilgen et al. 2014), so that prescribed burning is not necessary. However, the only invasive alien plant species recorded in this study (i.e. Lantana camara and Eucalyptus camaldulensis) were found in or near the riparian woodland plant community. Presumably, the river systems act as a pathway for propagule transport originating outside MWR. These species may have the potential to facilitate fire spread (Dew et al. 2017), and fire has previously been applied by management to control L. camara spread in MWR (Nieman et al. 2021). Although not excessively encroached in the reserve at the time of the study, these invasive species pose a potential threat to the indigenous vegetation and should be removed.

Conclusions
The process of compiling information for the classification and identification of woody plant communities for MWR in this study was achieved at a relatively low cost and during a short time, thus providing an accurate and efficient method to elucidate fine-scale vegetation patterns in an area where little prior information existed. The short duration of the survey in a single season may translate to an underrepresented species inventory, and future studies may therefore benefit from identifying species in a variety of seasons. For example, a number of species identified by Sherry (1989) were not found in this study, including Vachellia tortilis, Breonadia microcephala, Sterculia appendiculata, Pterocarpus angolensis and Strychnos madagascariensis. Furthermore, though edaphic factors were not assessed, they are likely an important determinant of plant distributions. For example, miombo is typically found on freely drained (leached), acidic soils with a restricted rooting environment (White 1983;Jinga & Palagi 2020), typical of the western uplands of MWR. For this reason, and many other related to environmental understanding and management, a detailed description of the geology and soils of MWR is required. Nevertheless, the woody vegetation classification and mapping of MWR achieved in this study is comprehensive, and more than adequate to effectively guide management decision-making on a variety of environmental topics, as well as policy development and species-or habitat-conservation strategies. The spatial boundaries of woody plant communities produced in this study are predominantly determined by factors that only change over geological time (e.g. slope, elevation and soils), and should therefore remain unchanged in the long-term. In contrast, the sub-type vegetation classifications (as well as grasslands), based predominantly on vegetation composition, physiognomy and structure, may be temporally altered through the effects of, among others, elephants and fire, and will therefore need to be updated periodically.   Zanha golungensis (smooth-fruit zanha) 1 (1.06%) Tw Table S1. List of tree and shrub species found in Majete Wildlife Reserve in this study, as well as the number of transects and the woody plant communities in which species occurred. The woody plant communities are: Rw = Riparian woodland, Gr = Grassland, SW = Shrublands and woodlands, Tw = Transitional woodland, and M = Miombo (continued).