In the Middle Eocene period, the true Rhino family (Rhinocerotidae) was a very abundant and diverse group of mammals and widely distributed throughout Africa, Asia, Europe and North America, with 41 genera and 142 species    . By the end of Pleistocene, most rhino species became extinct in most areas they used to exist mainly due to hunting and over-exploitation, hence, today the family is now confined to the sub-Saharan Africa and south-east Asia. black rhino (Diceros bicornis) is a critically endangered mega-herbivore of sub-Sahara Africa  . As recognized  , the species experiences reduced population due to illegal hunting caused by the illicit trade in rhino horn. In Zimbabwe, the black rhino species is specially protected, and the rhino population is below the minimum viable population level required to provide long term genetic security  , however, illegal hunting remains a problem   .
Following local extirpation of black rhinos in Gonarezhou National Park (GNP), southeastern Zimbabwe around the 1940s, the first two rhinos appear to have been released during 1969  . The second introduction of 77 black rhinos was conducted in the present GNP, Zimbabwe during the 1970s, and this was after the first local extirpation in the park. Unfortunately, poaching of wildlife in GNP increased, e.g. the number of black rhinoceros killed peaked to 51 in 1984 and by 1994 the black rhinos were locally extinct in GNP for the second time   . Meanwhile, at the peak of rhino poaching in Zimbabwe in the late 1980s, the Zimbabwe Parks and Wildlife Management Authority (ZPWMA), a state agency with the overall responsibility for wildlife protection, decided to relocate all black rhinos from vulnerable areas across the country and confine them to four state owned Intensive Protection Zones (IPZs) in Zimbabwe (i.e., Sinamatella, Matopos, Matusadona National Park and Chipinge Safari Area/Intensive Protection Zone) and private conservancies for enhanced protection  .
However, to improve the diversity of wildlife in the GNP and enhance its attractiveness as a tourist destination, the ZPWMA proposes to establish an IPZ in northern GNP with the black rhino species being re-introduced by the 2020s and beyond  . The IPZ proposal is evaluated in this research article with specific reference to habitat variables, i.e., browse availability and suitability for black rhinos, within northern GNP. Wildlife translocation, introduction or reintroduction regulations require feasibility or an ecological assessment to inform the reintroduction planning  . Hence, the present study is aimed at contributing to the existing body of knowledge on northern GNP and similar habitat status prior to the black rhino re-introduction project following the last assessment conducted in 2005  . Specifically, the study objective was to assess the browse availability and suitability for black rhino re-introduction in northern GNP, based on selected woody vegetation attributes.
2. Materials and Methods
2.1. Study Area
GNP is situated in the southeastern region of Zimbabwe between latitude 21˚15'S and longitude 30˚32'E and covering 5053 km2 (Figure 1). GNP constitutes part of the Greater Limpopo Transfontier Conservation Area, together with other protected areas like the Kruger National Park in South Africa and the Limpopo National Park in Mozambique. GNP and its surrounding areas experience short dry winter season in June and July with temperatures below 30˚C and a hot wet summer season from November to April when temperatures can exceed 40˚C. The remaining months are hot and dry  . The mean annual rainfall is 466 mm but is highly variable   . Two severe droughts have occurred since 1961 and the 1991/1992 drought led to the death of large numbers of wildlife  . The climate of northern GNP and its surrounding Safari areas is semi-arid  , and vegetation that is predominantly dry deciduous savanna woodland of varying types dominated by Colophospermum mopane woodland
Figure 1. Location of the two study strata of north and south of Runde River in Gonarezhou National Park, southern Zimbabwe (Source,  ).
and riverine and alluvial woodlands   . In 2016, the estimated population numbers of the principal large herbivores in GNP were: African elephant (Loxodonta africana) 10,715 (upper and lower 95% confidence limits ± 17.5%); impala (Aepyceros melampus) 7421 (±27.3%); African buffalo (Syncerus caffer) 4797 (±50.7%), kudu (Tragelaphus strepsiceros) 1789 (±27.6%); zebra (Equus quagga) 1830 (±26.2%); wildebeest (Connochaetes taurinus) 929 (±43.4%); giraffe (Giraffa camelopardalis) 446 (±36.8%); eland (Taurotragus oryx) 241 (±81.6%); nyala (Tragelaphus angasii) 187 (±47.1%) and waterbuck (Kobus ellipsiprymnus) 326 (±60.3%)  .
The study area which is the Save-Runde River sub-region in northern GNP has a total spatial extent of approximately 3000 km2  .
2.2. Data Collection
We stratified the study area into two strata, both sides of the perennial Runde River in GNP. The basis for stratifying on the left and right side of the river was premised on the fact that the perennial river with its high levels of water flow would tend to be a physical barrier to the rhinos preventing them to easily access all the sides of the banks of Runde River   . The data reported in this study were also used in a previous study where soil type was assessed on how it influences the structure and composition of woody species within the same study area of northern GNP  . A total of 96 sample plots measuring 20 m × 30 m (i.e., 48 plots in each study stratum) were randomly pegged based on GNP topographical map using grid intercept method, within the defined study area, between March and April 2011. Woody plant species observed within a sample plot were identified using the field guide  and recorded. The Linnean system of binomial nomenclature was used to classify woody plant species. Woody species taxa procedure method by comparison of an unknown woody species with named specimens, photographs, illustrations or descriptions were done using the keys in a southern Africa flora manual, guide handbook by  . In this present study, trees were defined as rooted, woody, self-supporting plants ≥ 3 m high with one or a few definite trunks (≥6 cm in basal stem diameter) whereas shrubs were defined as rooted, woody, self-supporting, multi-stemmed or single stemmed plants < 3 m high and <6 cm in stem basal diameter  . Specifically, the following woody plant variables were recorded: shrub height, shrub canopy diameters, evidence of browsing on woody plants, number of shrubs, and woody species name following the methods outlined by  . A visual method of recording canopy cover and shrub height was followed by an assessment of the woody species contributing to the available browse for black rhino as outlined by  . When plant resources are used disproportionately to their availability, use is said to be selective. The availability of plant resources was considered as quantity which can be accessibly to a population of black rhinos during the study time.
2.3. Data Analysis
The Shannon-Weiner diversity index (H’) was used to calculate the woody species diversity for each of the study stratum  . Data were tested for normality using Shapiro-Wilk test in STATISTICA for Windows version 6 software  . As a requirement to satisfy the assumptions of normality and equality of variance, data were normalized using: log10 (y + 1) where y is the recorded plant variable. Measured variables were tested for differences between the study strata using the two tailed Student t-test. Black rhino’s principal browse frequency of occurrence and proportional quality of woody species was categorized and analyzed following a standard procedure by  . Browse quality ranking was conducted in combination with the knowledge on the woody species’ value to potential black rhino diet following  known forage data as reference. This was in accordance to the study by  , which used the backtracking method to track black rhinos and record their diet composition assessing woody species browsed by black rhinos. With reference to available forage data on the diet of black rhinos, we determined principal wood species that constitute the diet of black rhino as recorded across the study area.
In order to determine which woody species are preferred by black rhinos we followed the method of  , where the proportional value of each woody species in the black rhino’s diet was divided by the abundance of that woody species available in a sample plot. A unit value which serves as an indicator for the black rhino’s level of preference of particular woody species was calculated following the formula of  as cited in  : Ei = (ri – ni)/(ri + ni) where E is the browse suitability/electivity measure for woody species i, r is the percentage of woody species i in the diet, and n is the percentage of woody species i in the environment (sample plot). This equation yields values ranging from −1 to 1. A value close to −1 indicates an avoidance of a woody species, a value of 0 shows no preference, and a value close to 1 illustrates a high preference for a woody species. Values of neither −1 nor 1 can be reached. Calculation of the browse suitability/electivity index was based on a list of woody plant attributes, important for the black rhinos, combined in a simple equation of Ei. Range of variation of each parameter was divided into classes to which scores of browse suitability/electivity index ranging from −1 to 1 which were attributed to the standard measures by  . In order to estimate the browse suitability/electivity index scores, we considered the successional stage of woody plants in the plots (height), share of principal taxa/forage woody species preferred by black rhinos, canopy volume of the shrub layer, percentage cover of shrubs. The browse suitability/electivity index was calculated after  , from the equation Ei. Principal woody species were defined as those that comprised ≥ 5% of the total number of bites made by the black rhino as earlier reported elsewhere in Matusadonha National Park, northern Zimbabwe  and  . Preferred woody species were determined using the browse suitability/electivity index described above and in  . Following several studies    , those taxa with browse suitability/electivity indexes greater than 0.4 and comprising at least 1% of total measured diet are reported and discussed in this present study.
A total of 3201 individual woody plants were recorded within the two study strata. Woody species diversity and black rhino’s browse suitability/electivity value significantly differed (p = 0.001) across the study strata (Table 1). However, there were no significant differences in shrub density (p = 0.398), shrub canopy volume (p = 0.893) and browsed plant density (p = 0.219).
The number of woody species was more regular in the northern side of the Runde River study stratum in comparison to the southern side of the Runde River (Table 2). Woody species recorded in northern GNP as common principal forage which dominated the standard reference of black rhino diet, were: knob thorn Acacia (Acacia nigrescens), small-leaved sickle-bush (Dichrostachys cineria), wild rubber horn pod tree (Diplorhynchus condlycarpon), tamboti (Spirostachys africana), and white-berry bushweed (Flueggea virosa) (Table 2), with different frequency of occurrence and proportion of availability.
With reference to standard forage data of black rhino diet, Acacia spp. comprised 5.34% of the total diet and had a browse suitability/electivity index of 0.70. Euclea spp. comprised 20.53% of the total diet and had a browse suitability/electivity index of 0.99, while the remainder of the Combretum spp. taxa comprised an additional 6.96% of total diet and possessed a browse suitability/electivity index of 0.92. The Faidherbia spp. comprised 7.6% of the total
Table 1. Vegetation variables in both the northern and southern Runde River strata: results of a two-tailed Student t-test analysis.
Table 2. Mean frequency of occurrence and percentage proportional browse availability of five principal forage woody species.
diet and had a browse/electivity index that approached 1. The Grewia genus comprised 42.71% of total diet, contributing by far the largest amount to black rhinodiet and had a browse suitability/electivity index of −0.07. These five principal taxa recorded in the present study comprised a total of 83.15% of the total diet of the black rhino (Table 3).
The study results recorded no significant differences in shrub density, shrub canopy volume and browsed plant density across the study strata (Table 1). This suggested presence of highly browsed woody plants and dominance of shrubs in some sections of northern GNP as reported by earlier studies    . However, this can impact positively to black rhinos feeding habits and pattern which are known as selective for woody plant species and size  . The study highlighted a high shrub density across the study strata (Table 1). According to earlier studies  and  , this is a key parameter of woody vegetation structure of importance to black rhino forage since mostly consumed shrubs are of a height below 2 m. The importance of woody vegetation to black rhinos largely depends on the availability of preferred quality browse within the reach (plant height ≤ 2 m) of the rhinos  . The synergy between the commonly reported fires and intense herbivory in GNP    will likely continue to contribute to the general modification of habitats, for example, reduced woody plant height and increased number of stems per plant, thereby, replacing woodlands with shrublands patches. This woody vegetation structural change likely to proliferate and persist across northern GNP, given the recorded no significant difference in shrub density and shrub canopy volume across the study strata could continue to set this section of the park as an ideal area for black rhino re-introduction.
Spirostachys, Dichrostachys and Acacia genera which are known palatable quality browse suitable for black rhino’s feed   were recorded with high frequency occurrence within the northern side of the Runde River study stratum, suggesting it to be a more suitable sanctuary space for the proposed re-introduction of black rhinos in northern GNP. Earlier studies in GNP confirmed abundance of Spirostachys woodlands  and Acacia woodlands  within northern GNP. The study results highlighted that north of Runde River
Table 3. Reference forage data: diet components of black rhinos (adapted from  ).
had higher woody species diversity and black rhino’s browse suitability/electivity values compared to the south of Runde River study stratum in northern GNP, southeastern Zimbabwe (Table 1, Table 2). Consistent with other studies, we recorded woody species of the following genera, the African bush willow (Combretum), mallow-leaved cross-berry (Grewia), bush guarri (Euclea), small-fruited potato bush (Phyllanthus), apple-ring acacia (Faidherbia), and chaste tree (Vitex), which are important food sources/forage for black rhinos in southern Africa   . The woody plant genera that rank high in the diet of black rhinos in Zimbabwe are diverse, namely, Acacia, Aloe, Bauhinia, Combretum, Commiphora, Cordia, Croton, Diospyros, Diplorhynchus, Euphorbia, Grewia, Vitex and the fallen fruits of Kigelia africana   . It was a positive finding that most of these woody species were recorded in northern GNP during the present study, over and above other recorded habitat forage status suggesting north of the Runde River as the most appropriate sanctuary of the proposed IPZ for black rhino re-introduction.
The study concluded that black rhino’s browse suitability, preferred browse woody species, woody species diversity and proportional browse availability, as selected attributes are favorably more pronounced in the northern side of the Runde River in GNP, southeastern Zimbabwe, suggesting that it is potentially suitable as a sanctuary for the black rhino re-introduction programme. We recommend that the black rhinos could be re-introduced in the northern GNP area provided other important factors such as law enforcement and other IUCN reintroduction requirements are taken into consideration.
We thank the Zimbabwe Parks and Wildlife Management Authority and Frankfurt Zoological Society for supporting this study.
*Formerly, Scientific Services, GonarezhouNational Park, Parks and Wildlife Management Authority, Private Bag 7003, Chiredzi, Zimbabwe.
 Cerdeno, E. (1998) Diversity and Evolutionary Trends of the Family Rhinocerotidae (Perissodactyla). Palaeogeography, Palaeoclimatology, Palaeoecology, 141, 13-34.
 Harley, E.H., Cunningham, J. and O’Ryan, C. (2005) Genetic Variation and Population Structure in Remnant Populations of Black Rhinoceros, Diceros bicornis in Africa. Molecular Ecology, 14, 2981-2990.
 IUCN International Union of Conservation and Nature (2017) Red List of Threatened Species. Version 2017-3.
 Dunham, K.M and van der Westhuizen, H.F. (2016) Aerial Survey of Elephants and other Large Herbivores in Gonarezhou National Park (Zimbabwe) and Some Adjacent Areas: 2016. Frankfurt Zoological Society, Gonarezhou Conservation Project, Gonarezhou National Park, Chiredzi.
 Gandiwa, E. and Kativu, S. (2009) Influence of Fire Frequency on Colophospermum mopane and Combretum apiculatum Woodland Structure and Composition in Northern Gonarezhou National Park, Zimbabwe. Koedoe, 51, Article No. 685.
 Gandiwa, E., Zisadza-Gandiwa, Goza, P.D., Mashapa, C. and Muboko, N. (2014) Diversityand Structure of Woody Vegetation across Areas with Different Soils in Gonarezhou National Park, Zimbabwe. Southern Forests: A Journal of Forest Sciences, 76, 111-116.
 Zisadza-Gandiwa, P., Chiganze, S., Chirombe, P., Mashapa, C., Muboko, N. and Gandiwa, E. (2013) Abundance and Distribution of African Fish Eagles along Major Rivers in Gonarezhou National Park, Zimbabwe. Greener Journal Biological Science, 3, 36-341.
 Matipano, G. (2003) A Comparison of Woody Browse Selection by Hand-Raised, Boma Adapted and Black Rhinoceros, (Diceros bicornis) in Matusadona National Park, Zimbabwe. Koedoe, 46, 83-96.
 Clark, J.H. (2013) Habitat Use Analysis of a Reintroduced Black Rhino (Diceros bicornis) Population. Thesis, Honors College Capstone Experience/Thesis Project, Paper 402, Western Kentucky University, Bowling Green.
 Gandiwa, E., Magwati, T., Zisadza, P., Chinuwo, T. and Tafangenyasha, C. (2011) The Impact of African Elephants on Acacia tortilis Woodland in Northern Gonarezhou National Park, Zimbabwe. Journal Arid Environment, 75, 809-814.
 Mashapa, C., Nyabawa, P., Zisadza-Gandiwa, P., Muvengwi, J., Kativu, S. and Gandiwa, E. (2014) Status of African Baobab (Adansonia digitata) across Gonarezhou National Park, Zimbabwe. Journal of Applied Science & Environmental Management, 18, 139-143.
 Kotze, D.C. and Zacharias, P.J.K. (1993) Utilisation of Woody Browse and Habitat by the Black Rhino (Diceros bicornis) in Western Itala Game Reserve. African Journal for Science, 10, 36-40.
 Emslie, H.R. and Adcock, A. (1994) Feeding Ecology of the Black Rhinoceros. In: Penzhorn, B.L. and Kreik, N.P.L., Eds., Proceedings of a Symposium on Rhinos as Game Ranch Animal, Wildlife Group of the South African Veterinary Association, Onderstepoort, 9-10 September 1994, i-iv, 1-242.
 Muya, S.M. and Oguge, M. (2000) Effects of Browse Availability and Quality on Black Rhino (Diceros bicornis) Diet in Nairobi National Park, Kenya. African Journal of Ecology, 38, 62-71.
 Zisadza-Gandiwa, P., Mango, L., Gandiwa, E., Goza, D., Parakasingwa, C., Chinoitezvi, E., Shimbani, J. and Muvengwi, J. (2013) Variation in Woody Vegetation Structure and Composition in a Semi-Arid Savanna of Southern Zimbabwe. International Journal of Biodiversity and Conservation, 5, 71-77.
 Mashapa, C., Zisadza-Gandiwa, P., Gandiwa, E. and Kativu, S. (2014) Abundance and Structure of Baobab (Adansonia digitata) across Various Soil Group Substrate Strata in Gonarezhou National Park, Southeast Zimbabwe. International Journal of Biodiversity, 2013, Article ID: 874713.
 Gandiwa, E., Gandiwa, P. and Mxoza, T. (2012) Structure and Composition of Spirostachys africana Woodland Stands in Gonarezhou National Park, Southern Zimbabwe. International Journal of Environmental Science, 2, 234-241.
 Oloo, T.W., Brett, R. and Young, T.P. (1994) Seasonal Variation in the Feeding Ecology of Black Rhinos (Diceros bicornis L.) in Laikipia, Kenya. African Journal of Ecology, 32, 142-157.