Abstract Introduction Methods Results Ethical considerations Discussion Acknowledgements References About the Author(s) Keoikantse Sianga Okavango Research Institute, University of Botswana, Botswana Richard W.S. Fynn Okavango Research Institute, University of Botswana, Botswana Mpaphi C. Bonyongo Okavango Research Institute, University of Botswana, Botswana Citation Sianga, K., Fynn, R.W.S. & Bonyongo, M.C., 2017, ‘Seasonal habitat selection by African buffalo Syncerus caffer in the Savuti–Mababe–Linyanti ecosystem of northern Botswana’, Koedoe 59(2), a1382. https://doi.org/10.4102/koedoe.v59i2.1382 Original Research Seasonal habitat selection by African buffalo Syncerus caffer in the Savuti–Mababe–Linyanti ecosystem of northern Botswana Keoikantse Sianga, Richard W.S. Fynn, Mpaphi C. Bonyongo Received: 26 Jan. 2016; Accepted: 24 Nov. 2016; Published: 23 May 2017 Copyright: © 2017. The Author(s). Licensee: AOSIS. This is an Open Access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Abstract This study aimed to establish seasonal movement and habitat selection patterns of African buffalo Syncerus caffer in relation to a detailed habitat map and according to seasonal changes in forage quality and quantity in the Savuti–Mababe–Linyanti ecosystem (Botswana). Two buffalo were collared in November 2011 and another in October 2012. All three buffalo had greater activities in the mopane–sandveld woodland mosaic during the wet season, which provided high-quality leafy grasses and ephemeral water for drinking, but moved to permanent water and reliable forage of various wetlands (swamps and floodplains) and riverine woodlands during the dry season. Wetlands had higher grass greenness, height and biomass than woodlands during the dry season. Buffalo had similar wet season concentration areas in the 2011–2012 and 2012–2013 wet seasons and similar dry season concentration areas over the 2012 and 2013 dry seasons. However, their dry season location of collaring in 2011 differed dramatically from their 2012 and 2013 dry season concentration areas, possibly because of the exceptionally high flood levels in 2011, which reduced accessibility to their usual dry season concentration areas. The study demonstrates that extremely large and heterogeneous landscapes are needed to conserve buffalo in sandy, dystrophic ecosystems with variable rainfall. Conservation implications: This study emphasises the importance of large spatial scale available for movement, which enables adaptation to changing conditions between years and seasons. Introduction African buffalo (Syncerus caffer) have a large distributional range across the savannas of Africa. Their habitat selection and foraging ecology have been relatively well studied (Landman & Kerley 2001; Macandza, Owen-Smith & Cross 2004; Sinclair 1979; Taylor 1985). Because of their large body size and large groups, buffalo are able to fend off predators (Sinclair, Mduma & Brashares 2003), thereby enabling them to forage in wooded vegetation with relatively low visibility and high predation risk. The large absolute food demands of buffalo, which is a function of their large body size (Illius & Gordon 1987; Wilmshurst, Fryxell & Bergman 2000), combined with their inability to efficiently crop short grass (use of tongue to increase bite size cannot work on short grass) constrain them to foraging in vegetation with sufficient height and biomass of forage (Codron et al. 2008; Illius & Gordon 1987). Buffalo generally avoid heavily grazed regions of short grassland (Bhola et al. 2012; Jarman & Sinclair 1979; Traill & Bigalke 2007), preferring woodlands dominated by tufted and leafy perennial grasses during the wet season and often relying on more productive riverine vegetation during the dry season (Bell 1970; Bennitt, Bonyongo & Harris 2014; Fynn, Chase & Roder 2014; Macandza, Owen-Smith & Cain 2012; Sinclair 1979). While drylands often support less productive but higher quality forage than wetlands over the wet season (Fynn et al. 2014; Taylor 1985), their moisture-limited position in the landscape results in forage drying out and declining in forage quality and quantity over the dry season (Ellis & Swift 1988; Owen-Smith 2008; Taylor 1985). Thus, wetlands, where permanent water is easily accessible and perennial grasses continue to produce green forage over the dry season, would be expected to be favoured by many herbivores at this time with drylands being favoured over the wet season (Bell 1970; Fynn et al. 2014; Maddock 1979; Taylor 1985). In contrast to more mobile species such as zebra (Equus burchellii) and wildebeest (Connochaetes taurinus), which often migrate between distinct wet- and dry season ranges (Bartlam-Brooks, Bonyongo & Harris 2011; Maddock 1979), most buffalo studies have found strong overlap between wet- and dry season ranges (Macandza et al. 2012; Ryan 2006; Ryan, Knechtel & Getz 2006; Venter & Watson 2008). Many studies, however, were conducted in small reserves (Landman & Kerley 2001; Ryan et al. 2006; Tshabalala, Dube & Lent 2009) where opportunity for large seasonal movements are limited or in larger reserves where functional heterogeneity of resources may be well developed at the landscape (catena) scale (Bell 1970; Macandza et al. 2012; Perrin & Brereton-Stiles 1999; Sinclair 1979), thereby reducing the need for large seasonal movements (Hopcraft, Olff & Sinclair 2010). However, in large, relatively unfragmented ecosystems and where landscape scale functional heterogeneity may be poorly developed relative to regional-scale heterogeneity (Fynn et al. 2014), larger buffalo movements may be expected (Naidoo et al. 2012; Skarpe et al. 2004). The Savuti–Mababe–Linyanti ecosystem (SMLE) is part of the > 80 000 km2 northern conservation area of Botswana, one of the largest relatively unfragmented wildlife regions in Africa (Fynn & Bonyongo 2011). Despite several studies on buffalo in the SMLE (Fynn et al. 2014; Patterson 1972), detailed seasonal movement and habitat selection patterns of buffalo in the region have not been fully identified and established, especially seeing that no detailed habitat map has been available until this year (2016). Considering their non-specialised mouth anatomy and their large body size–mediated demands for absolute food intake, we expected buffalo in the SMLE to favour vegetation where tufted, leafy grasses provided sufficient height and biomass of grass to enable food intake requirements to be satisfied (Illius & Gordon 1987; Wilmshurst et al. 2000) (Hypothesis 1 – H1). The strongly developed regional-scale distribution (rather than landscape scale) of functional heterogeneity of wet- and dry season vegetation types in the SMLE (Fynn et al. 2014) is likely to promote extension of buffalo home ranges beyond the landscape scale (Hopcraft et al. 2010). Consequently, we expected seasonal movement patterns of buffalo in the SMLE to match the scale at which functional seasonal vegetation types are distributed in the ecosystem, which is strongly regional (Fynn et al. 2014) (Hypothesis 2 – H2). The objectives of this study were to (1) to examine seasonal movements and habitat selection of buffalo in the SMLE of northern Botswana and (2) to link seasonal movements to the quality and quantity (grass greenness, height and biomass) of vegetation in favoured seasonal regions of the landscape. Methods Study area This study was conducted in the SMLE (northern Botswana, Figure 1). Climate in the study area is described as semi-arid with mean annual rainfall ranging from around 500 mm in the Okavango region (western boundary of the SMLE) to over 600 mm in the Chobe Enclave (eastern boundary of the SMLE), most of which is received between November and April (Botswana Meteorological Services). Seasons may be functionally separated into a wet season (December–April), a cool early dry season (May–August) and a hot late dry season (September–November), where daily maximum temperatures are between 35 °C and 40 °C (Figure 2, adapted from Fynn et al. 2014). The SMLE is defined by a variety of vegetation types from swamps, floodplains and riverine in wetter areas to vast dryland woodland systems away from permanent water comprising a mosaic of mopane woodland on alluvial soils and sandveld woodland on Kalahari sands (Figure 1; Sianga & Fynn 2017; Wolski & Murray-Hudson 2006). In addition, very heavy clay soils formed under sedimentation in a paleolake system known as the Mababe Depression (Teter 2007) give rise to open grassland and sparse savanna in the eastern section of the SMLE (Figure 1). Another key feature are the dambo grasslands of the Chobe Enclave, which are seasonally flooded shallow, linear depressions (flooded by runoff from upslope regions) supporting tall grassland. FIGURE 1: Savuti–Mababe–Linyanti ecosystem vegetation. FIGURE 2: Rainfall and temperature patterns between 2010 and 2013. (a) 2010/2011, (b) 2011/2012 and (c) 2012/2013 rainfall season. Buffalo habitat use Single cows in three buffalo herds were fitted with satellite collars (Africa Wildlife Tracking, Pretoria, South Africa). Buffalo 1 (BH1) and Buffalo 2 (BH2) were collared in the Savuti Channel floodplains and near the Savuti Marsh, respectively, at the end of the late dry season in November 2011. Buffalo 3 (BH3) was collared in the Savuti Marsh in the late dry season (October 2012). Anaesthetic M3080 Xylazine was used to immobilise cows during capturing and later reversed using naltrexone after fitting collars. Collars were programmed to take a global positioning system (GPS) position of the animal every 4 hours. Every fix obtained by the satellite collars was downloaded and plotted in ArcGIS 10.1 (ESRI 2010) for various analyses. For analyses relating to seasonal differences in movement patterns and habitat selection, we defined six seasonal periods to which each GPS fix could be allocated. The seasonal periods were early wet season (mid-November to end of December), mid wet season (January and February), late wet season (March and April), early dry season (May and June), mid dry season (July and August) and late dry season (September to mid-November). For visual presentation of seasonal movements and locations for each buffalo, we plotted the home ranges of each buffalo (colour coded for each of the six seasonal periods) in a local convex hull (a-LoCoH, 95% isopleths) (Getz et al. 2007) in R (RCore-Team 2013) and later displayed as a shapefile on a habitat map of the region. The habitat map was developed from a detailed vegetation classification and mapping study of the SMLE funded by Southern African Science Service Centre for Climate Change and Adaptive Land Management (SASSCAL); a detailed account of the vegetation communities and habitat map can be seen in an accompanying paper in this special issue on Botswana (Sianga & Fynn 2017). LoCoH has been found to be an appropriate tool in GPS studies (Getz et al. 2007). Vegetation sampling GPS data from satellite collars were used to locate sites where buffalo have been during the wet and dry seasons. Sampling of buffalo’s wet season vegetation types was conducted in the 2012–2013 wet season, while their dry season vegetation types were sampled in the 2012 dry season. A total of 124 samples (grass greenness, height and biomass) were collected from sites where the buffalo have been. Seventy-five samples (20 and 55) were collected for BH1 during the early and late dry season of 2012, while 20 samples were collected for BH2 in the late dry season of 2012. At each site, five 0.25 m2 quadrats were set-up at the GPS position obtained from the collar. The first quadrat was set at the 0-m position on the hand-held GPS (Garmin GPS Map 62s) and the other four quadrats 5 m each side of the first (central) quadrat. Greenness was estimated visually as the percentage of green tissue of grasses and sedges rooted within the quadrat. Grass height was measured by lowering a brown paper sample bag at the centre of each quadrat, and the height above the soil surface was measured. All grasses rooted within the quadrats were clipped at ground surface level and air-dried during the field exercise. The air-dried grasses were oven-dried at 60 °C for 48 hours and weighed for biomass at the Okavango Research Institute Laboratory (Maun, Botswana). In addition, to determine the vegetation type, we noted grasses and trees common within approximately a 10-m radius around the central quadrat. Statistical analyses Home ranges and habitat selection For analysis of habitat selection, we determined seasonal minimum convex polygons (MCPs) in ArcGIS 10.1 (ESRI 2010) and local convex hulls in R (Getz et al. 2007; RCore-Team 2013) using seasonal location data of each collared buffalo. Seasonal habitat selection indices of each buffalo were determined by dividing the proportion of use of each habitat by the proportion of availability of each habitat following Jacobs Index (J.I) = (r ‑ p)/(r + p – 2rp), where r is the proportion of habitat used, and p is the proportion of habitat available (Jacobs 1974). J.I ranges between –1 (selected against) and +1 (selected for). To consider the effects of scale on resource availability (Gustine et al. 2006), we calculated habitat selection at second and third order. Second-order selection (Johnson 1980) was determined by comparing the availability of various habitat types in the individual buffalo MCPs against availability in the overall population MCP (Thomas & Taylor 1990), while third-order selection (Johnson 1980) was determined by comparing the availability of various habitat types in the individual buffalo local convex hulls to availability in the respective individual buffalo MCPs (Thomas & Taylor 1990). Owing to having only three collared individual buffalo, there was insufficient replication to effectively statistically test habitat selection. Thus, the J.I merely provides a guide to readers on potential habitat selection. To determine home-range overlaps between seasons, we used Intersect Tool (Universal Transverse Mercator, ESRI 2010) to estimate area percentage overlaps, where % overlaps = area overlap of two seasons/(area of season 1 + season 2) × 100. To examine how buffalo moved in relation to distance from perennial water sources differed seasonally, we determined the distance to the nearest perennial water source of every GPS location of all buffalo throughout the study period using Near Tool in ArcGIS 10.1 (Universal Transverse Mercator, ESRI 2010). Vegetation data Data of grass greenness, height and biomass were subjected to Shapiro–Wilk test (test of normality) and Levene statistic (test of homogeneity of variance) in R version 2.15.2 (RCore-Team 2013). Data sets that failed assumptions of normality or homogeneity of variance were natural log transformed, and failure to meet these assumptions after transformation led to the use of non-parametric tests (Kruskal–Wallis test). Multiple comparison tests (kruskalmc’s function) in ‘pgirmess’ package in R (RCore-Team 2013) was used to analyse data sets. Results All three buffalo had much greater activities in sandveld and mopane woodlands during the wet season and various wetland systems during the dry season (Figure 3). BH1 and BH2 had greater separation of their wet and dry season concentration areas than BH3 (quantified by percentage overlap of their wet and dry season 95% local convex hull isopleths – Table 1; also see Figure 4 vs. Figure 5). BH1 and BH2 demonstrated distinct wet and dry season concentration areas, returning to these respective seasonal locations in successive years (Table 1; Figure 6). FIGURE 3: Percentage of various vegetation types observed in the follow-ups of GPS collar locations during the wet season and late dry season. FIGURE 4: Movement patterns of buffalo 1 (BH1) and buffalo 2 (BH2) over the 2011/2012 annual cycle, in relation to vegetation. FIGURE 5: Movement patterns of buffalo 3 (BH3) over the 2012/2013 annual cycle, in relation to vegetation. FIGURE 6: Movement patterns of buffalo 1 (BH1) and 2 (BH2) in the mid wet and mid dry season of 2012 and 2013, in relation to vegetation. TABLE 1: Percentage overlaps of home ranges by buffalo in different seasons over the 2012 and 2013. For second-order habitat selection, the buffalo selected mopane and sandveld woodlands over the wet season, whereas riverine and wetland habitats were selected over the dry season (Table 2); for third-order habitat selection, the buffalo selected mopane woodlands over the wet season and acacia grasslands, whereas riverine and wetland habitats were selected over the dry season (Table 3). For more details of the plant species composition and classification of these plant communities and for details of the habitat map development, see a parallel paper in this special issue (Sianga & Fynn 2017). Over both dry seasons, BH1 concentration areas were in wetland and adjacent woodlands where it focused on the Linyanti Swamp in increasing proportions (and regularly across the international border with Namibia) as the dry season progressed (Figure 4). The concentrations of BH2 were in the Selinda Spillway and Kwando River wetland (floodplains) and adjacent woodlands in the 2012 and 2013 dry seasons (Figures 4 and 6). Interestingly, BH2 concentrated its activities in the Tsum Tsum wetlands (floodplains) of the Okavango Delta in the late wet season of April 2013 (which it did not do in 2012) before moving to the Selinda Spillway later in the dry season (Figure 6). BH3 had the wetland habitats of the Savuti Channel and Savuti Marsh at the core of its range, which it focused on during the dry season, and utilised adjacent woodlands on the peripheries of the home range during the wet season (Figure 5). TABLE 2: Second-order habitat selection (annual minimum convex polygons vs. wet or dry seasons minimum convex polygons) by three buffalo in the Savuti–Mababe–Linyanti ecosystem (northern Botswana). TABLE 3: Third-order habitat selection (annual minimum convex polygons vs. wet or dry season local convex hull) by three buffalo in the Savuti–Mababe–Linyanti ecosystem (northern Botswana). During the early and late dry season of 2012, wetland habitats (floodplain not sampled in the early dry season) had higher grass greenness, height and biomass than mopane, sandveld and riverine habitats, respectively (p < 0.05; Figure 7), while mopane, sandveld and riverine habitats were not significantly different in grass greenness, height and biomass, respectively (p > 0.05; Figure 7). FIGURE 7: (a) Grass greenness, (b) height and (c) biomass of mopane woodland, sandveld, riverine and wetland used by buffalo herds during the early and late dry season. Ethical considerations A veterinarian registered with the government of Botswana conducted all darting operations in principle with the research permit (EWT 8/36/4 XVII [31]) and the supplementary collaring permit (EWT 8/36/4 XVII [41]) issued by the Ministry of Environment, Wildlife and Tourism and the Department of Wildlife and National Parks (Gaborone, Botswana). Darting operations or removing collars from cows were done from a vehicle and a helicopter, respectively. Collars were removed successfully after two successive years. Discussion Buffalo shifted their seasonal concentration areas between wetlands and woodlands (Figures 3 and 4) as an adaptive strategy to seasonal forage dynamics and water availability over the annual cycle, with BH1 and BH2 having a migratory strategy with little overlap between wet and dry season concentration areas. However, the distance moved from perennial water sources over the wet season varied significantly across the buffalo, with BH3 being more sedentary. Bennitt et al. (2014) and Taolo (2003) reported similar results, noting that buffalo selected contrasting seasonal vegetation types, using woodlands far from permanent water during the rainy season and seasonally flooded vegetation close to permanent water during the early and late flood seasons. As with our findings, several other studies in the region have also found a mix of strategies (migratory vs. sedentary) between buffalo (Bennitt et al. 2014; Bennitt, Bonyongo & Harris 2015; Naidoo et al. 2012, 2014). Cornélis et al. (2011) also noted similar findings where buffalo were observed to range between woodlands and riverine habitats during the wet and dry season, respectively. A key factor that may influence buffalo concentration patterns during the wet season appears to be the availability of higher quality, soft-leaved grasses (such as Digitaria eriantha) in back country woodlands than in wetlands (floodplains), which in contrast support tougher-leaved grasses adapted to shading and litter accumulation under the more productive conditions of wetlands (Fynn & Bonyongo 2011), and where greater biomass dilutes nutrient concentrations (Jarrell & Beverly 1981), while greater silica contents of wetland versus dryland grasses further reduce digestibility (Mosimane 2015). One of the most favoured grasses by buffalo, D. eriantha, peaks in abundance beyond 10 km from the permanent water of the Okavango Delta and Linyanti Swamps, whereas shorter grazing lawn–forming grasses such as Urochloa trichopus peaked in abundance with 5 km of these extensive wetland systems (Sianga et al. 2017). The medium to tall, leafy high-quality D. eriantha would be expected to be optimal in height and digestibility for maximising protein and energy intake rates in buffalo, which explains why buffalo focused their wet season foraging beyond 10 km from wetland systems (Figure 4). By contrast, wildebeest and impala, which favour shorter grasses, focused within 5 km of these wetland systems (Fynn et al. 2014). This suggests that short and tall grass grazers in the SMLE have distinct spatial niche separation during the wet season driven by distance to wetland systems. Similar observations have been made in the Maasai-Mara, where buffalo favour the taller grass areas within the reserve and wildebeest and gazelles the shorter grass areas outside the reserve (Bhola et al. 2012). Thus, a high conception rate of buffalo during the wet season is probably because of the availability of high-quality green forage of optimal height and biomass (Ryan, Knechtel & Getz 2007). Thus, concentration during the wet season by the three buffalo in vegetation types (sandveld and mopane) where tufted, soft-leaved grasses are more abundant supports our first hypothesis (H1). Drying up of ephemeral pans over the dry season forced buffalo to move away from favoured wet season ranges towards permanent water sources (Figures 3 and 4), which included floodplains and swamps (wetlands) and adjacent woodlands within the study area. It is this forced movement away from favoured wet season ranges towards permanent water that likely contributes to reduced grazing pressure in favoured wet season ranges and, therefore, higher abundances of favoured grasses such as D. eriantha far from permanent water. This indicates that artificial water provision in back country woodlands far from permanent water may be an unwise management option that could reduce niche diversity in the system leading to declines in taller grass grazers such as buffalo, as well as roan and sable antelope. These floodplain and swamp (wetlands) systems provided buffalo not only with reliable drinking water during the dry season but also with forage significantly greener than that found in woodlands, as observed in other studies in the region (Bennitt et al. 2014; Fynn et al. 2014; Taolo 2003). More productive wetland areas are likely to provide a critical reserve and critical buffer resources (Owen-Smith 2002) for the late dry season, or elsewhere referred to as key resources (Illius & O’Connor 2000), especially during drought years when more preferred resources have been depleted (Owen-Smith 2002). The ability of wetlands to provide taller green forage with higher energy and protein levels relative to drylands during the late dry season has been shown to greatly elevate buffalo population productivity (Taylor 1985). The taller grass of wetlands is likely to be especially critical for buffalo, which are highly vulnerable to competition for forage (Bhola et al. 2012), owing to their tongue sweep foraging strategy, which cannot deal with short grass. Thus, our findings of increasing use of swamp, marsh and floodplains (wetlands) as the dry season progressed by our three buffalo together with similar observations in other studies (Prins & Beekman 1989; Taylor 1985; Tinley 1977; Vesey-FitzGerald 1960; Western 1973), demonstrate the importance of wetland systems as dry season key-resource habitat for buffalo, as is the case for cattle and other wild herbivores across Africa (Fynn et al. 2015). Our collared buffalo migrated at an intermediate scale from wetlands to woodlands rather than from wetlands to the high-quality Mababe Depression grasslands, favoured by a migratory zebra population over the wet season (Sianga 2014) and, thus, our second hypothesis (H2) that buffalo seasonal movements would match the regional-scale distribution of functional seasonal vegetation types was not supported. Though this finding does not support our H2, it does not rule out the possibility of buffalo migrating into the depression during the wet season because not all buffalo herds in the ecosystem were fitted with collars to determine their movements. Also, aerial surveys conducted by Fynn et al. (2014) did not observe any buffalo in the Mababe Depression at that time, suggesting that buffalo probably do not favour the Mababe Depression during the wet season, unlike the zebra population (Sianga 2014). However, the buffalo appear to show adaptive variation in seasonal presence to changing environmental conditions in different years. For instance, in 2011, BH1 and BH2 avoided floodplains and swamps (wetlands) in their usual dry season concentration areas of the Linyanti Swamps and Selinda Spillway (2012 and 2013) because of the exceptionally high floods of 2011, being forced to use the more elevated floodplains of the Savuti Channel. This shift in location of different dry season concentration areas emphasises the importance of having a large spatial scale available for movement (Fynn et al. 2014), which enables adaptation to changing conditions between years. Finally, the finding of cross-border movements by BH1, together with those of Naidoo et al. (2014) and Patterson (1972), demonstrates the importance of Transfrontier Conservation Areas such as this one (Kavango Zambezi Transfrontier Conservation Area, KAZA), which transcends the borders of Angola, Namibia, Botswana, Zimbabwe and Zambia allowing adaptive foraging to a variety of different functional seasonal resources. Acknowledgements We thank the German Ministry of Education and Research (BMBF) who funded this work through ‘The Future Okavango project’ (buffalo study) and through Southern African Science Service Centre for Climate Change and Adaptive Land Management (SASSCAL) (vegetation study and buffalo habitat selection analyses). Dr Andy Lyons is acknowledged for his guidance with TLoCoH. Dr Rob Jackson is appreciated for buffalo darting, collar fitting and removal. Finally, we thank the Ministry of Environment, Wildlife and Tourism and the Department of Wildlife and National Parks (Botswana) for permitting us to conduct research in northern Botswana. Competing interests The authors declare that they have no financial or personal relationships that may have inappropriately influenced them in writing this article. Authors’ contributions All authors of this research article have directly participated in the planning, execution and analysis of this study. References Bartlam-Brooks, H., Bonyongo, M. & Harris, S., 2011, ‘Will reconnecting ecosystems allow long-distance mammal migrations to resume? A case study of a zebra Equus burchelli migration in Botswana’, Oryx 45, 210–216. https://doi.org/10.1017/S0030605310000414 Bell, R., 1970, ‘The use of the herb layer by grazing ungulates in the Serengeti’, in A. Watson (ed.), Animal populations in relation to their food resources, pp. 111–124, Blackwell, Oxford. Bennitt, E., Bonyongo, M.C. & Harris, S., 2014, ‘Habitat selection by African buffalo (Syncerus caffer) in response to landscape-level fluctuations in water availability on two temporal scales’, PLoS One 9, e101346. https://doi.org/10.1371/journal.pone.0101346 Bennitt, E., Bonyongo, M.C. & Harris, S., 2015, ‘Behaviour-related scalar habitat use by Cape buffalo (Syncerus caffer)’, PLoS One 10, e0145145. https://doi.org/10.1371/journal.pone.0145145 Bhola, N., Ogutu, J.O., Piepho, H.P., Said, M.Y., Reid, R.S., Hobbs, N.T. et al., 2012, ‘Comparative changes in density and demography of large herbivores in the Masai Mara Reserve and its surrounding human-dominated pastoral ranches in Kenya’, Biodiversity & Conservation 21, 1509–1530. https://doi.org/10.1007/s10531-012-0261-y Codron, D., Brink, J.S., Rossouw, L., Clauss, M., Codron, J., Lee-Thorp, J.A. et al., 2008, ‘Functional differentiation of African grazing ruminants: An example of specialized adaptations to very small changes in diet’, Biological Journal of the Linnean Society 94, 755–764. https://doi.org/10.1111/j.1095-8312.2008.01028.x Cornélis, D., Benhamou, S., Janeau, G., Morellet, N., Ouedraogo, M. & Visscher, M.N., 2011, ‘Spatiotemporal dynamics of forage and water resources shape space use of West African savanna buffaloes’, Journal of Mammalogy 92(6), 1287–1297. https://doi.org/10.1644/10-MAMM-A-397.1 Ellis, J.E. & Swift, D.M., 1988, ‘Stability of African pastoral ecosystems: Alternate paradigms and implications for development’, Journal of Range Management Archives 41, 450–459. https://doi.org/10.2307/3899515 ESRI, 2010, ArcGIS desktop, Environmental Systems Research Institute, Redlands, CA. Fynn, R.W.S. & Bonyongo, M.C., 2011, ‘Functional conservation areas and the future of Africa’s wildlife’, African Journal of Ecology 49, 175–188. https://doi.org/10.1111/j.1365-2028.2010.01245.x Fynn, R.W.S., Chase, M. & Roder, A., 2014, ‘Functional habitat heterogeneity and large herbivore seasonal habitat selection in northern Botswana’, South African Journal of Wildlife Research 44, 1–15. https://doi.org/10.3957/056.044.0103 Fynn, R.W.S., Murray-Hudson, M., Dhliwayo, M. & Scholte, P., 2015, ‘African wetlands and their seasonal use by wild and domestic herbivores’, Wetlands Ecology and Management 23, 559–581. https://doi.org/10.1007/s11273-015-9430-6 Getz, W.M., Fortmann-Roe, S., Cross, P.C., Lyons, A.J., Ryan, S.J. & Wilmers, C.C., 2007, ‘LoCoH: Nonparametric kernel methods for constructing home ranges and utilization distributions’, PLoS One 2, e207. https://doi.org/10.1371/journal.pone.0000207 Gustine, D.D., Parker, K.L., Lay, R.J., Gillingham, M.P. & Heard, D.C., 2006, ‘Interpreting resource selection at different scales for woodland caribou in winter’, Journal of Wildlife Management 70, 1601–1614. https://doi.org/10.2193/0022-541X(2006)70[1601:IRSADS]2.0.CO;2 Hopcraft, J.G.C., Olff, H. & Sinclair, A., 2010, ‘Herbivores, resources and risks: Alternating regulation along primary environmental gradients in savannas’, Trends in Ecology & Evolution 25, 119–128. https://doi.org/10.1016/j.tree.2009.08.001 Illius, A.W. & Gordon, I.J., 1987, ‘The allometry of food intake in grazing ruminants’, The Journal of Animal Ecology 56, 989–999. https://doi.org/10.2307/4961 Illius, A.W. & O’connor, T.G., 2000, ‘Resource heterogeneity and ungulate population dynamics’, Oikos 89, 283–294. https://doi.org/10.1034/j.1600-0706.2000.890209.x Jacobs, J., 1974, ‘Quantitative measurement of food selection’, Oecologia 14, 413–417. https://doi.org/10.1007/BF00384581 Jarman, P.J. & Sinclair, A.R.E., 1979, ‘Feeding strategy and the pattern of resource partitioning in African ungulates’, in A.R.E. Sinclair & M. Norton-Griffiths (eds.), Serengeti: Dynamics of an ecosystem, pp. 130–163, University of Chicago Press, Chicago, IL. Jarrell, W. & Beverly, R., 1981, ‘The dilution effect in plant nutrition studies’, Advances in Agronomy 34, 197–224. https://doi.org/10.1016/S0065-2113(08)60887-1 Johnson, D.H, 1980, ‘The comparison of usage and availability measurements for evaluating resource preference’, Ecology 61, 65–71. https://doi.org/10.2307/1937156 Landman, M. & Kerley, G., 2001, ‘Dietary shifts: Do grazers become browsers in the Thicket Biome?’, Koedoe 44, 31–36. https://doi.org/10.4102/koedoe.v44i1.183 Macandza, V.A., Owen-Smith, N. & Cain, J.W., 2012, ‘Habitat and resource partitioning between abundant and relatively rare grazing ungulates’, Journal of Zoology 287, 175–185. https://doi.org/10.1111/j.1469-7998.2012.00900.x Macandza, V.A., Owen-Smith, N. & Cross, P.C, 2004, ‘Forage selection by African buffalo in the late dry season in two landscapes’, South African Journal of Wildlife Research 34, 113–121. Maddock, L., 1979, ‘The migration and grazing succession’, in A.R.E. Sinclair & M. Norton-Griffiths (eds.), Serengeti: Dynamics of an Ecosystem, pp. 104–129, University of Chicago Press, Chicago, IL. Mosimane, K., 2015, ‘A first quantification of biogenic silica pools in an evaporation-dominated wetland system: The Okavango Delta’, MPhil thesis, University of Botswana. Naidoo, R., Du Preez, P., Stuart-Hill, G., Beytell, P. & Taylor, R., 2014, ‘Long-range migrations and dispersals of African buffalo (Syncerus caffer) in the Kavango–Zambezi Transfrontier Conservation area’, African Journal of Ecology 52, 581–584. https://doi.org/10.1111/aje.12163 Naidoo, R., Du Preez, P., Stuart-Hill, G., Jago, M. & Wegmann, M., 2012, ‘Home on the range: Factors explaining partial migration of African buffalo in a tropical environment’, PLoS one 7, e36527. https://doi.org/10.1371/journal.pone.0036527 Owen-Smith, N., 2002, Adaptive herbivore ecology: From resources to populations in variable environments, Cambridge University Press, Cambridge. Owen-Smith, N., 2008, ‘Changing vulnerability to predation related to season and sex in an African ungulate assemblage’, Oikos 117, 602–610. https://doi.org/10.1111/j.0030-1299.2008.16309.x Patterson, L., 1972, Studies on the population biology of the African buffalo Syncerus caffer in Botswana, Department of Wildlife and National Parks, Gaborone, Botswana. Perrin, M. & Brereton-Stiles, R., 1999, ‘Habitat use and feeding behaviour of the buffalo and the white rhinoceros in the Hluhluwe-Umfolozi Game Reserve’, South African Journal of Wildlife Research 29, 72–80. Prins, H.H.T. & Beekman, J.H., 1989, ‘A balanced diet as a goal for grazing: The food of the Manyara buffalo’, African Journal of Ecology 27, 241–259. https://doi.org/10.1111/j.1365-2028.1989.tb01017.x Rcore-Team, 2013, R: A language and environment for statistical computing, R Foundation for Statistical Computing, Vienna, Austria, viewed 11 April 2015, from http://www.R-project.org Ryan, S.J., 2006, ‘Spatial ecology of African buffalo and their resources in a savanna Ecosystem’, PhD thesis, University of California, USA. Ryan, S.J., Knechtel, C.U. & Getz, W.M., 2006, ‘Range and habitat selection of African buffalo in South Africa’, Journal of Wildlife Management 70, 764–776. https://doi.org/10.2193/0022-541X(2006)70[764:RAHSOA]2.0.CO;2 Ryan, S.J., Knechtel, C.U. & Getz, W.M., 2007, ‘Ecological cues, gestation length, and birth timing in African buffalo (Syncerus caffer)’, Behavioral Ecology 18, 635–644. https://doi.org/10.1093/beheco/arm028 Sianga, K., 2014, ‘Habitat use by buffalo and zebra in relation to spatial and temporal variability of resources in the Savuti-Mababe-Linyanti ecosystem of northern Botswana’, Mphil thesis, University of Botswana, Botswana. Sianga, K. & Fynn, R., 2017, ‘The vegetation and wildlife habitats of the Savuti-Mababe-Linyanti ecosystem, northern Botswana’, Koedoe 59(2), a1406. https://doi.org/10.4102/koedoe.v59i2.1406 Sianga, K., Van Telgen, M., Vrooman, J., Fynn, R.W.S. & Van Langevelde, F., 2017, ‘Spatial refuges buffer landscapes against homogenisation and degradation by large herbivore populations and facilitate vegetation heterogeneity’, Koedoe 59(2), a1434. https://doi.org/10.4102/koedoe.v59i2.1434 Sinclair, A.R.E., 1979, ‘The eruption of the ruminants’, in A.R.E. Sinclair & M. Norton-Griffiths (eds.), Serengeti – Dynamics of an ecosystem, pp. 82–103. Chicago University Press, Chicago, IL. Sinclair, A.R.E., Mduma, S. & Brashares, J.S., 2003, ‘Patterns of predation in a diverse predator–prey system’, Nature 425, 288–290. https://doi.org/10.1038/nature01934 Skarpe, C., Aarrestad, P.A., Andreassen, H.P., Dhillion, S.S., Dimakatso, T., Du Toit, J.T. et al., 2004, ‘The return of the giants: Ecological effects of an increasing elephant population’, Ambio: A Journal of the Human Environment 33, 276–282. Taolo, C., 2003, ‘Population ecology, seasonal movement and habitat use of the African buffalo (Syncerus caffer) in Chobe National Park, Botswana’, PhD thesis, Norwegian University of Science and Technology, Trondheim, Norway. Taylor, R.D. 1985, ‘The response of buffalo syncerus caffer to the Kariba lakeshore grasslands (Panicum repens) in Matusadona National Park’, PhD thesis, University of Zimbabwe, Zimbabwe. Teter, K.L., 2007, ‘Paleoenvironmental reconstruction of Paleolake Mababe, northwestern Botswana from sediment chemistry and biological productivity data’, MSc thesis, Oklahoma State University, Stillwater, OK. Thomas, D.L. & Taylor, E.J., 1990, ‘Study designs and tests for comparing resource use and availability’, The Journal of wildlife management 54, 322–330. https://doi.org/10.2307/3809050 Tinley, K.L., 1977, ‘Framework of the Gorongosa Ecosystem’, PhD dissertation, University of Pretoria, South Africa. Traill, L.W. & Bigalke, R.C., 2007, ‘A presence-only habitat suitability model for large grazing African ungulates and its utility for wildlife management’, African Journal of Ecology 45, 347–354. https://doi.org/10.1111/j.1365-2028.2006.00717.x Tshabalala, T., Dube, S. & Lent, P.C., 2009, ‘Seasonal variation in forages utilized by the African buffalo (Syncerus caffer) in the succulent thicket of South Africa’, African Journal of Ecology 48, 438–445. https://doi.org/10.1111/j.1365-2028.2009.01132.x Venter, J.A. & Watson, L.H., 2008, ‘Feeding and habitat use of buffalo (Syncerus caffer) in the Nama-Karoo, South Africa’, South African Journal of Wildlife Research 38, 42–51. https://doi.org/10.3957/0379-4369-38.1.42 Vesey-Fitzgerald, D.F., 1960, ‘Grazing succession among east African game animals’, Journal of Mammalogy 41, 161–172. https://doi.org/10.2307/1376351 Western, D., 1973, ‘The structure, dynamics and changes of the Amboseli ecosystem’, PhD dissertation, University of Nairobi, Kenya. Wilmshurst, J.F., Fryxell, J.M. & Bergman, C.M., 2000, ‘The allometry of patch selection in ruminants’, Proceedings of the Royal Society of London. Series B: Biological Sciences 267, 345–349. https://doi.org/10.1098/rspb.2000.1007 Wolski, P. & Murray-Hudson, M., 2006, ‘Flooding dynamics in a large low-gradient alluvial fan, the Okavango Delta, Botswana, from analysis and interpretation of a 30-year hydrometric record’, Hydrology and Earth System Sciences Discussions 2, 1865–1892. https://doi.org/10.5194/hessd-2-1865-2005