Hostname: page-component-78c5997874-mlc7c Total loading time: 0 Render date: 2024-11-13T02:24:51.415Z Has data issue: false hasContentIssue false

A newly discovered wildlife migration in Namibia and Botswana is the longest in Africa

Published online by Cambridge University Press:  27 May 2014

R. Naidoo*
Affiliation:
WWF–US, 1250 24th Street NW, Washington, DC, USA.
M. J. Chase
Affiliation:
Elephants Without Borders, Kasane, Botswana
P. Beytell
Affiliation:
Directorate of Natural Resource Management, Ministry of Environment & Tourism, Windhoek, Namibia
P. Du Preez
Affiliation:
Directorate of Natural Resource Management, Ministry of Environment & Tourism, Windhoek, Namibia
K. Landen
Affiliation:
Elephants Without Borders, Kasane, Botswana
G. Stuart-Hill
Affiliation:
WWF in Namibia, Windhoek, Namibia
R. Taylor
Affiliation:
WWF in Namibia, Windhoek, Namibia
*
(Corresponding author) E-mail robin.naidoo@wwfus.org
Rights & Permissions [Opens in a new window]

Abstract

Migrations of most animal taxa are declining as a result of anthropogenic pressures and land-use transformation. Here, we document and characterize a previously unknown multi-country migration of Burchell's zebra Equus quagga that is the longest of all recorded large mammal migrations in Africa. Our data from eight adult female zebras collared on the border of Namibia and Botswana show that in December 2012 all individuals crossed the Chobe River and moved due south to Nxai Pan National Park in Botswana, where they spent a mean duration of 10 weeks before returning, less directly, to their dry season floodplain habitat. The same southward movements were also observed in December 2013. Nxai Pan appeared to have similar environmental conditions to several possible alternative wet season destinations that were closer to the dry season habitat on the Chobe River, and water availability, but not habitat or vegetation biomass, was associated with higher-use areas along the migratory pathway. These results suggest a genetic and/or cultural basis for the choice of migration destination, rather than an environmental one. Regardless of the cause, the round-trip, straight-line migration distance of 500 km is greater than that covered by wildebeest Connochaetes taurinus during their well-known seasonal journey in the Serengeti ecosystem. It merits conservation attention, given the decline of large-scale ecological processes such as animal migrations.

Type
Papers
Copyright
Copyright © Fauna & Flora International 2014 

This paper contains supplementary material that can be found online at http://journals.cambridge.org

Introduction

Seasonal migration, where individuals move to and from geographically separated home ranges to exploit changes in environmental conditions, is common among many taxonomic groups (Dingle & Drake, Reference Dingle and Drake2007). The long-distance seasonal migrations of monarch butterflies Danaus plexippus in North America, wildebeest Connochaetes taurinus in the Serengeti, songbirds across the Americas, and grey Eschrichtius robustus and humpback whales Megaptera novaeangliae in the Pacific Ocean are among the most spectacular natural phenomena (Wilcove, Reference Wilcove2009). However, migrations across many different taxa and in many locations have been extinguished or are under threat from anthropogenic pressures such as habitat destruction, construction of barriers, overexploitation and climate change (Wilcove & Wikelski, Reference Wilcove and Wikelski2008; Dobson et al., Reference Dobson, Borner, Sinclair, Hudson, Anderson and Bigurube2010). In addition, migrations are often poorly understood and many may remain uncharacterized, as evidenced by the rediscovery of a migration of large mammals in Sudan that rivals that of the Serengeti in abundance, which remained largely undocumented until 2007 (Harris et al., Reference Harris, Thirgood, Hopcraft, Cromsigt and Berger2009; Wilcove, Reference Wilcove2009). Yet understanding migrations and the movement ecology of wildlife is an important consideration in the design of functionally connected landscapes, in the management of species of conservation concern and for the protection of threatened natural phenomena (Brower & Malcolm, Reference Brower and Malcolm1991; Briers, Reference Briers2002; Crooks & Sanjayan, Reference Crooks and Sanjayan2006; Fynn & Bonyongo, Reference Fynn and Bonyongo2011; Epps et al., Reference Epps, Mutayoba, Gwin and Brashares2011).

Here we document for the first time the existence of a cross-border, long-distance migration of Burchell's zebra Equus quagga that ranks as the longest of all known mammal migrations in Africa. This previously undiscovered phenomenon involves several thousand zebras covering a straight-line distance of 500 km over several months, spans parts of Namibia and Botswana and was observed in consecutive years. We explore the influence of the external environment in shaping the migratory route and destination, and hypothesize that this migration may represent a conserved, ancient route rather than a response to current environmental conditions.

Methods

Fieldwork was conducted during the dry season in 2012, in mid September on the Botswana side of the Chobe River and mid October on the Namibia side (Fig. 1). Adult female zebras (n = 8) were darted from the air or ground and immobilized using a mixture of etorphine hydrochloride, azaperone and hyaluronidase. The age and family group size were estimated for each individual, and a satellite-tracking collar attached. Tracking collars (Africa Wildlife Tracking, Pretoria, South Africa) were programmed to record global positioning system (GPS) locations at either 4- or 5-hour intervals. The total number of data points per individual was 562–1,403 (mean 1,008 ± SD 215). Animals were captured and treated according to the protocols approved under research permit 1,537/2010 from the Ministry of Environment and Tourism in Namibia, and EWT 8/36/4 XVII (57) from the Department of Wildlife and National Parks in Botswana.

Fig. 1 Movement trajectories of eight female zebra Equus quagga collared on the Chobe floodplains in Botswana and Namibia from the late dry season (September–October) 2012 until June 2013. Grey polygons indicate actual (Nxai Pan National Park) and potential alternative (Savuti Marsh and Seloko Plains) migration destinations. The rectangle on the inset indicates the location of the main map in southern Africa.

GPS locations were used to map and animate movement trajectories of all individuals. We calculated net squared displacement distances from the collaring site over time to characterize the scale and timing of migration (Bunnefeld et al., Reference Bunnefeld, Börger, van Moorter, Rolandsen, Dettki, Solberg and Ericsson2011).

We assessed whether migration to Nxai Pan could be driven by unique environmental conditions at the site, first using spatial data layers on soil fertility (a proxy for the nutrition content of forage, and therefore a critical determinant of grazing potential for herbivores such as zebras; Fynn & Bonyongo, Reference Fynn and Bonyongo2011) in Botswana. Nxai Pan is an important wet season habitat for wildlife because of the filling of waterholes and the growth of new grass on the pan (Department of Wildlife and National Parks, 1995). Nxai Pan soils are categorized as fertile to very fertile and are therefore attractive to grazers because of their high nutritional value (Lindsay et al., Reference Lindsay, Muir and Robson1998). There are only two other locations in northern Botswana that have similarly fertile soils; both are closer to the Chobe River and zebras could reach them without having to navigate around fences or other potential anthropogenic barriers (Fig. 1). The Seloko Plains are c. 60 km from zebra dry season ranges along the Chobe River and have fertile to very fertile soils (Lindsay et al., Reference Lindsay, Muir and Robson1998). A small section of the eastern part of these grasslands has been converted to commercial farming but the remainder is unconverted and zebras and other wildlife are present (Chase, Reference Chase2013). The Savuti Marsh/Matabe Depression (Savuti Marsh hereafter) is c. 100 km from zebra dry season ranges and has moderately fertile to fertile soils. Zebras and other ungulates are known to migrate to Savuti from the Chobe and Linyanti (the Linyanti is the name of the Chobe River further westwards) floodplains (Vandewalle, Reference Vandewalle2000), and the grasslands of this marsh support large herds of wildlife (Fynn & Bonyongo, Reference Fynn and Bonyongo2011; Sianga et al., Reference Sianga, Bonyongo and Fynn2013). These three sites (Nxai Pan, Seloko Plains and Savuti Marsh) are also broadly classified as tree and grass savannah (Lindsay et al., Reference Lindsay, Muir and Robson1998). Remaining habitats accessible to zebras moving southwards from the Chobe River all occur on moderate- or low-fertility soils that would have limited appeal to migrating grazers.

We characterized environmental conditions at Nxai Pan and at the two potential alternative destinations (Seloko Plains and Savuti Marsh); our hypothesis was that Nxai Pan should appear significantly different to the other two sites if the environment was driving zebra migration there. We used available data on vegetation biomass, tree cover and rainfall, and averaged the values of these variables over all cells in each defined site area (Fig. 1). We assessed whether temporal changes in mean vegetation biomass varied across these three potential wet season migration destinations, using the MODIS 250-m resolution enhanced vegetation index, which is produced at 16-day intervals (Pettorelli et al., Reference Pettorelli, Vik, Mysterud, Gaillard, Tucker and Stenseth2005). We used remotely sensed data from the Tropical Rainfall Monitoring Mission (half-degree resolution) to characterize daily precipitation levels and cumulative precipitation at Nxai Pan, Seloko Plains and Savuti Marsh during September 2012–March 2013. Finally, we characterized tree cover at each of these three sites by averaging the percentage of each cell covered by woody vegetation, using 250-m MODIS continuous field data from 2010 (Hansen et al., Reference Hansen, DeFries, Townshend, Carroll, Dimiceli and Sohlberg2003).

We conducted an aerial census of zebras at Nxai Pan, using a fixed-wing aircraft, on 21 February 2013. Strip transects were flown, using standard methodology (Norton-Griffiths, Reference Norton-Griffiths1978). All zebras within 200 m of either side of a transect were counted, with 27 parallel transects flown during 07.30–10.00. This design ensured 100% coverage of the pan and resulted in a complete census of individuals in the area. High-resolution photographs were taken for more accurate counting post-flight.

To assess the relationship between environmental conditions and zebra movements along the main migratory pathway we calculated the minimum convex polygon across all GPS locations of all individuals moving from the Chobe River to Nxai Pan and back. We then overlaid a 5 × 5 km grid across this polygon and calculated the number of GPS points in each cell separately for two time periods: the southward (4 December 2012–10 January 2013) and the northward (18 February 2013–13 June 2013) migratory phases. We classified grid cells as stopover areas when at least one zebra spent at least 48 consecutive hours in an encamped state in the cell (Morales et al., Reference Morales, Haydon, Frair, Holsiner and Fryxell2004; Sawyer & Kauffman, Reference Sawyer and Kauffman2011).

Using high-resolution imagery from Google Earth we assessed the availability of standing water to zebras in each grid cell as absent (no visual evidence of standing water), present (1–3 waterholes) or abundant (> 3 waterholes and/or the presence of a stream or river). All cells were categorized by the same observer to maintain consistency. We were unable to assess whether waterholes contained standing water that was accessible to zebras at the time they were encountered, only whether cells contained waterholes or watercourses that had the potential to be filled. Similar to the minimum convex polygon environmental analyses above we also calculated mean cell fractions covered by woody vegetation and mean enhanced vegetation index values for each 5 × 5 km cell, separately for the southward and northward movement phases.

We conducted four separate statistical analyses, using generalized linear models and multi-model inference, to investigate how environmental characteristics of cells along the migratory pathway affected zebra movements.

Southward migration

(1) $${\rm Prob}(Y = y_i ) = \displaystyle{{e^{ - u_i} [u_i ]y_i} \over {y_i !}},y_i = 0,1,2,{\rm \ldots} $$

where

$$\eqalign{{\rm In}(u_i ) =\ & a + \beta _1 *{\rm Water} + \beta _2 *{\rm Woody} \cr & + \beta _3 *{\rm EVI}_{{\rm South}} + \beta _4 *{\rm EVI}_{{\rm South}^2}} $$

We used a negative binomial regression model to estimate counts of zebra GPS fixes in cells (y i ) as a function of waterholes (Water), the mean cell fraction covered by woody vegetation (Woody), and a second-order polynomial of mean enhanced vegetation index values; a is the regression intercept and β is the regression coefficient of the associated independent variable.

(2) $$\eqalign{{\rm Prob}.{\rm Stopover}_{{\rm South}} = \displaystyle{\matrix{e[a + \beta _1 *{\rm Water} + \beta _2 *{\rm Woody} \cr +\, \beta _3 *{\rm EVI}_{{\rm South}} + \beta _4 *{\rm EVI}_{{\rm South}^2} ] } \over \matrix{1 + e[a + \beta _1 *{\rm Water} + \beta _2 *{\rm Woody} \cr +\, \beta _3 *{\rm EVI}_{{\rm South}} + \beta _4 *{\rm EVI}_{{\rm South}^2} ] }} $$

We used a logistic regression model to investigate the probability of a cell being used as a stopover location as a function of waterholes (Water), the mean cell fraction covered by woody vegetation (Woody), and a second-order polynomial of mean enhanced vegetation index values; a is the regression intercept and β is the regression coefficient of the associated independent variable.

Northward migration

(3) $${\rm Prob}(Y = y_i ) = \displaystyle{{e^{ - u_i} [u_i ]y_i} \over {y_i !}},y_i = 0,1,2,{\rm \ldots} $$

where

$$\eqalign{{\rm In}(u_i ) =\ & a + \beta _1 *{\rm Water} + \beta _2 *{\rm Woody} \cr & + \beta _3 *{\rm EVI}_{{\rm North}} + \beta _4 *{\rm EVI}_{{\rm North}^2}} $$

Variable descriptions are as in (1), with the counts of GPS fixes and the enhanced vegetation index measurements made during the northward migration phase.

(4) $$\eqalign{{\rm Prob}.{\rm Stopover}_{{\rm North}} = \displaystyle{\matrix{e[a + \beta _1 *{\rm Water} + \beta _2 *{\rm Woody} \cr + \beta _3 *{\rm EVI}_{{\rm North}} + \beta _4 *{\rm EVI}_{{\rm North}^2} ] } \over \matrix{1 + e[a + \beta _1 *{\rm Water} + \beta _2 *{\rm Woody} \cr +\, \beta _3 *{\rm EVI}_{{\rm North}} + \beta _4 *{\rm EVI}_{{\rm North}^2} ] }} $$

Variable descriptions are as in (2), with stopover and enhanced vegetation index measured during the northward migration phase.

Note that we used negative binomial rather than Poisson count data regression because the variance of zebra GPS counts in the cells was substantially greater than the mean of the counts, indicating overdispersed count data (Cameron & Trivedi, Reference Cameron and Trivedi1998). Residuals from each of the preliminary global regression models also showed significant first-order spatial autocorrelation in the residuals, and therefore in addition to the terms noted in equations (14) we included a first-order neighbourhood term that contained the sum of the count values of the nine cells in the immediate neighbourhood of the target cell (Dormann et al., Reference Dormann, McPherson, Araujo, Bivand, Bolliger and Carl2007). We hypothesized that after correcting for spatial autocorrelation, cells with more waterholes, an optimal enhanced vegetation index value and lower levels of woody vegetation would be used more frequently and would be more likely to be used as stopover points than those with lower levels of each. We tested these hypotheses using a multi-model inference approach, evaluating all 64 possible combinations of independent variables described in equations (14). Model-averaged coefficients were calculated based on each model's Akaike weight (Burnham & Anderson, Reference Burnham and Anderson1998), with means, standard errors, confidence limits, and tests of significance calculated using MuMin v. 1.9.13 (Barton, Reference Barton2013), and all other analyses performed in R v.2.15.1 (R Development Core Team, 2008).

Results

Adult zebras spent several weeks in the vicinity of the Chobe River floodplains, regularly crossing the river to and from Namibia and Botswana. On 26–27 November 2012 all individuals remaining in Namibia crossed the river south into Botswana, with three females moving off the floodplains and adjoining woodlands to an area c. 30 km south in Botswana (Fig. 1). On 3 December 2012 the first significant directional movement due south was observed in one individual, after which seven of the eight zebras moved south to the Nxai Pan, arriving during 14–27 December (Table 1). The remaining female zebra spent several weeks south of Seloko Plains, eventually arriving at Nxai Pan on 9 January 2013. By December 2013 the surviving GPS-collared zebras (n = 6) were again making their way south from the Chobe floodplains towards Nxai Pan.

Table 1 Summary statistics for migration movements of eight female zebras Equus quagga collared on the Chobe floodplains in Botswana and Namibia (Fig. 1).

* Died on 13 May 2013

The straight-line distance between the outer extent of the Chobe dry season range and the pan at Nxai Pan National Park is almost exactly 250 km. Most of this ground was covered during a 2-week period in December 2012 (median 14 days, range 7–37) during which movement was highly linear and unidirectional, with some individuals travelling > 50 km per day. Total distances travelled by individuals on their southward journey to Nxai Pan ranged from 256 km for an individual with highly directional movement to 424 km for an individual that deviated to Seloko Plains before eventually moving to Nxai Pan (mean 300 ± SD 54 km; Table 1).

The timing of the movement away from the Chobe floodplains appears to be related to the amount and timing of rainfall at the destination in Nxai Pan rather than at the Chobe origin. A major rainfall event (37 mm, nearly 10% of the region's annual precipitation) occurred in the Nxai Pan area on 4 December 2012, several days before the zebras began their rapid, unidirectional movement south (Fig. 2). There was no rainfall event of comparable magnitude in the Chobe region before the long-distance migration southward.

Fig. 2 Timing of daily rainfall and zebra migration during September 2012–May 2013. The thick black line represents the mean net displacement away from dry season home ranges for nine zebras. Vertical lines represent daily rainfall over dry season ranges (black) and the Nxai Pan wet season range (grey).

Zebras spent a mean of 73 days in the Nxai Pan area (range 53–96), mostly on the pan itself. We counted 1,534 zebras in the area during our aerial survey on 21 February 2013, although correspondence with safari guides in the park indicated that in the 2 weeks prior to our survey there had been more zebras in the vicinity of the pan.

Return trajectories of zebras moving from the Nxai Pan north to the Chobe River were less direct and took significantly longer. The mean return time was 85 days (vs 15 days for the southward journey; t = 4.33, P = 0.004) and mean distance travelled was 653 ± SE 153 km (vs 302 km; t = 5.79, P = 0.0006). The mean distance of actual movement trajectories of the seven zebras that survived to complete the entire return journey from the Chobe floodplains to Nxai Pan was 955 km (range 735–1,170 km).

Environmental characteristics of Nxai Pan overlapped with the two possible alternative wet season destinations for Chobe River zebras. Woody vegetation cover, a measure of broad habitat type, was similar at all three potential destinations (Nxai Pan, 18.0 ± SE 3.7%; Seloko Plains, 17.1 ± SE 2.6%; Savuti Marsh, 20.0 ± SE 5.3%), as was total rainfall during September 2012–March 2013 (Nxai Pan, 444 mm; Seloko Plains, 477 mm; Savuti Marsh, 470 mm). The area of Nxai Pan used by zebras had a consistently lower mean vegetation biomass than Savuti Marsh. However, biomass levels at Nxai Pan and Seloko Plains were similar from mid December to early February, which encompasses most of the wet season period that zebras were present at Nxai Pan (Fig. 3). Biomass levels at Nxai Pan dropped below those of Seloko Plains in mid February; at the same time the zebras began to move back north towards the Chobe River floodplains (Table 1).

Fig. 3 Temporal pattern of mean enhanced vegetation index (September 2012–May 2013) at the Nxai Pan wet season migration destination and two alternative destinations (Seloko Plains and Savuti Marsh). The ellipse indicates the near-identical values at Nxai Pan and Seloko Plains during the core wet season time when zebras were at Nxai Pan.

As suggested by initial analyses, averaged coefficients across all model combinations showed that spatial autocorrelation in zebra location data was strong, with coefficients on the first-order neighbourhood variable significantly greater than zero in all four analyses (Table 2). Controlling for this autocorrelation the probability of a 5 × 5 km cell being used as a stopover point (i.e. where zebras stopped for at least 48 hours during the migratory phases) appeared unrelated to waterhole abundance, enhanced vegetation index or woody vegetation cover (Table 2). The number of GPS relocations in cells along the southward migratory pathway also appeared unrelated to waterhole abundance, enhanced vegetation index or woody vegetation cover. However, cells with higher numbers of GPS fixes during the northward migratory phase contained higher densities (and presence) of waterholes than those used less frequently (Table 2). The top 10 models for each of our four analyses accounted for 68–100% of the cumulative Akaike model weight, indicating that a relatively small fraction of models (15%) had high explanatory power (Supplementary Table S1).

Table 2 Coefficients for regression models relating environmental characteristics to the count of zebra locations in 5 × 5 km cells along the migratory pathway (negative binomial regression model) and to the probability of a cell being a stopover point on the migration route (logistic regression model) during both northward and southward movement phases.

1 Coefficients significantly different from 0 (P < 0.05) are shown in bold.

2 Coefficients on Waterholes present and Waterholes abundant are relative to reference level of factor variable (Waterholes absent).

Discussion

Prior to the discovery of this migration it was unclear where Chobe River zebras moved to in the wet season, with local residents and wildlife professionals only aware that animals left at the beginning of the wet season and returned the following dry season (Chase, Reference Chase2011). Similarly, the appearance of zebras on the Nxai Pan in the wet season is an annual occurrence but there was no evidence that animals originated as far away as the Chobe River. The timing of the migration from the Chobe River appeared to coincide with heavy, episodic rainfall events early in the wet season at Nxai Pan. Other migratory ungulate species (Sinclair & Norton-Griffiths, Reference Sinclair and Norton-Griffiths1979), as well as zebras in the recently re-established migration from the Okavango Delta to the Makgadikgadi Pans (Bartlam-Brooks et al., Reference Bartlam-Brooks, Bonyongo and Harris2011), also appear to time migration to rainfall events, and zebras are able to orient long-range movements to environmental conditions outside their perceptory range (Brooks & Harris, Reference Brooks and Harris2008). However, as we have only observed two southward migration onsets we are unable to assess statistically the nature of the relationship between timing of rainfall and migration without a substantially longer time series.

Our aerial census showed that there were at least 1,500 zebras on Nxai Pan in the wet season of 2012–2013, with anecdotal information indicating this number may have been substantially higher a few weeks before our census. This figure is similar to dry season estimates of zebra numbers (1,558–3,593) on the Chobe floodplains during 2007–2012, which were also estimated from strip-transect surveys (MC, unpubl. data). Combined with the facts that all eight of the female zebras we collared in 2012 moved from the Chobe River to Nxai Pan, the six surviving zebras repeated this southward movement to Nxai Pan in December 2013, and zebra numbers on the Chobe floodplains during the wet season (255) are an order of magnitude lower than during the dry season (Chase, Reference Chase2013), this suggests that most of the zebras on the Chobe floodplains may be moving to Nxai Pan in the wet season. Further censuses and satellite tracking of more individuals over time could confirm the annual periodicity of the migration, as well as the fraction of Chobe animals that move to Nxai Pan vs those that remain or migrate elsewhere. Additional studies are also needed to determine whether other ungulates such as the wildebeest, long-distance migrants that move in tandem with zebras in other ecosystems (e.g. Serengeti, Liuwa Plains—Mussuma Transfrontier Conservation area; Harris et al., Reference Harris, Thirgood, Hopcraft, Cromsigt and Berger2009), perform the same Chobe–Nxai Pan migration.

The reasons for such a long-distance migration, with closer and apparently suitable alternative destinations by-passed, remain unknown. Our analyses of potential environmental drivers of the migration are preliminary and our hypothesis that the migration is a fixed behaviour rather than a variable response to environmental conditions needs further testing. Seloko Plains and Savuti Marsh have similar levels of woody vegetation, soil fertility and rainfall to Nxai Pan, and wet season vegetation biomass was similar for Seloko Plains and Nxai Pan. Vegetation biomass was higher at Savuti Marsh than Nxai Pan, which, if anything, suggests greater grazing availability and therefore higher suitability for zebras migrating from the Chobe floodplains. We were unable to quantify access to drinking water at any of the three locations but at the time of migration Savuti Marsh held water for the first time in 28 years and was therefore much wetter than either Seloko Plains or Nxai Pan. Access to water would therefore not have been a problem at that site, although it is possible that having been dry for so long the marsh is still an unfavoured destination that will need to be rediscovered by migratory ungulates. Also, we could not quantify predation risk, which is an important driver of migration in ungulates (Fryxell & Sinclair, Reference Fryxell and Sinclair1988), although lions Panthera leo, the main predators of zebras in this region, are present in sufficient numbers to draw tourists in both Savuti Marsh and Nxai Pan (Power & Compion, Reference Power and Compion2009; Johnson, Reference Johnson2010).

With no obvious environmental differences between Nxai Pan and the other two possible destinations for zebras from Chobe River, this migration may represent an ancient, conserved phenomenon that has a genetic basis and/or is an outcome of learned behaviour and cultural transmission between individuals and social groups (Alerstam, Reference Alerstam2006; Bolger et al., Reference Bolger, Newmark, Morrison and Doak2008). Research has shown that the Okavango–Makgadikgadi zebra migration in Botswana, itself a long-distance event involving c. 15,000 zebras, reoccurred during 2008–2009 for the first time following the removal of a veterinary fence that had blocked the route during 1968–2004. This completely separate migration suggests a conserved memory of an ancient route that may be genetically encoded (Bartlam-Brooks et al., Reference Bartlam-Brooks, Bonyongo and Harris2011), and raises the question of why zebras from Chobe River move to Nxai Pan rather than the nearby Makgadikgadi Pans (and vice versa). This fixed-route phenomenon has also been observed in pronghorn Antilocapra americana in the USA, where individuals in the Grand Teton National Park area have used the same seasonal migration pathway, as narrow as tens of metres in some places, since the Holocene era (Berger et al., Reference Berger, Cain and Berger2006).

Our analyses imply that zebra movement and stopovers along the migratory pathway between dry and wet season endpoints are unrelated to either spatial or temporal variation in vegetation, and cells with waterholes are not used preferentially as stopover points, although cells with waterholes had higher numbers of zebra GPS counts than those without for the northward (but not southward) migration. As with the actual and alternative wet season destinations there does not appear to be a strong response to fluctuating environmental conditions along the migratory pathway.

Comparing our data with those compiled in a review (Harris et al., Reference Harris, Thirgood, Hopcraft, Cromsigt and Berger2009) shows that the newly discovered zebra migration ranks as the longest large-mammal migration in Africa. The round-trip length of this migration is longer than long-distance movements of zebras and wildebeest in the Serengeti ecosystem (A.R.E. Sinclair, pers. comm.; Harris et al., Reference Harris, Thirgood, Hopcraft, Cromsigt and Berger2009). It is also longer than migrations of topi Damaliscus lunatus, Mongalla gazelle Gazella thomsonii albonotata, and white-eared kob Kobus kob in southern Sudan (Harris et al., Reference Harris, Thirgood, Hopcraft, Cromsigt and Berger2009) and of wildebeest in Angola and Zambia (Harris et al., Reference Harris, Thirgood, Hopcraft, Cromsigt and Berger2009), although data for these species and locations are incomplete. Although migrations of greater distances occurred in the past (Harris et al., Reference Harris, Thirgood, Hopcraft, Cromsigt and Berger2009) there are no extant large-mammal migrations in Africa that are longer than the one we document here.

This discovery is important in light of the effects of anthropogenic activities on large-mammal migrations in and around the study area during the last half-century. The erection of veterinary control fences during 1950–1980 disrupted many large-scale migrations, severing connectivity between the dry season grassland habitats of the Makgadikgadi Pans, Savuti Marsh and the Kalahari Desert and wet season strongholds in the Okavango Delta and the Chobe/Linyanti Rivers (Williamson & Williamson, Reference Williamson and Williamson1984; Mbaiwa & Mbaiwa, Reference Mbaiwa and Mbaiwa2006; Ferguson & Hanks, Reference Ferguson and Hanks2010). As a consequence, populations of migratory ungulate species such as the zebra, wildebeest and buffalo Syncerus caffer declined and, decades later, continue to show little sign of recovery (Perkins, Reference Perkins, Ferguson and Hanks2010). Nor are these disruptions consigned to the distant past; in 2005 a fence erected along the Boteti River excluded most of the river from the Makgadikgadi Pans National Park, rendering adjacent dry season habitat inaccessible and resulting in the deaths of > 300 migrating zebras that were confronted with an unknown barrier (Gibson, Reference Gibson, Ferguson and Hanks2010).

As with the Serengeti migration, the Namibia–Botswana zebra migration involves the crossing of an international border, and therefore cross-country cooperation will be required to ensure its persistence. Given the documented disappearance of numerous wildlife migrations, the speed at which migratory populations can be lost as a result of anthropogenic perturbations (Bolger et al., Reference Bolger, Newmark, Morrison and Doak2008; Wilcove, Reference Wilcove2009) and the growing recognition of wildlife-based economies as engines for rural development in sub-Saharan Africa (Lindsey et al., Reference Lindsey, Roulet and Romanach2007; Jones et al., Reference Jones, Davis, Diez, Diggle, Roe, Elliott, Sandbrook and Walpole2012) we suggest that ensuring the persistence of one of Africa's longest wildlife migrations should be a high conservation priority. The proposed Kavango–Zambezi Transfrontier Conservation Area (KAZA TFCA, 2013), spanning parts of Namibia, Botswana, Angola, Zambia and Zimbabwe, covers the entire migration route and therefore would provide an ideal platform from which to pursue targeted international conservation activities. These should include continued tracking of individuals to assess year-to-year variation in movements, an assessment of potential threats and impediments to migration (e.g. the construction of new fences or roads, the expansion of agriculture or human settlements), coordinated monitoring of population vital rates and the development of innovative mechanisms (such as payments for environmental services) that provide incentives for the conservation of the migration corridor.

In a human-dominated world in which wildlife migrations are disappearing, the discovery of a previously undocumented migration of zebras across hundreds of kilometres and over international borders provides hope that the conservation of large-scale ecological phenomena, among other biodiversity conservation goals, continues to be a worthwhile pursuit in the Anthropocene era (Caro et al., Reference Caro, Darwin, Forrester, Ledoux-Bloom and Wells2012).

Acknowledgements

We thank Hans Swartbooi, Lise Hanssen, Carl-Heinz Moeller, Michelle Kastern, Werner Kilian and Johannes Kapner for help in the field, Alicia Ellis, Colby Loucks, Sylvia Thompson, David Wilcove, Tony Sinclair, Dave Ward, Raymond Peters and two anonymous reviewers for their constructive comments, and the Namibia Ministry of Environment and Tourism, the Botswana Department of Wildlife and National Parks, the Eppley Foundation, the Paul G. Allen Family Foundation, and the San Diego Zoo Global Wildlife Conservancy for funding and logistical support.

Biographical sketches

Robin Naidoo is a conservation scientist who has worked in Africa on various issues for over a decade. Mike Chase is an expert on elephant ecology and conservation. Piet Beytell is a conservation biologist with the Namibian government. Pierre Du Preez is Namibia's chief conservation biologist and runs the country's rhino recovery programme. Kelly Landen is based in Botswana and helps run Elephants Without Borders. Greg Stuart-Hill has decades of experience in community-based conservation and rangeland management in southern Africa. Russell Taylor is a regional expert in the conservation of large wildlife and natural resources management.

Footnotes

* Also at: Institute for Conservation Research, San Diego Zoo Global, Escondido, USA
This paper contains supplementary material that can be found online at http://journals.cambridge.org

References

Alerstam, T. (2006) Conflicting evidence about long-distance animal navigation. Science, 313, 791794.CrossRefGoogle ScholarPubMed
Bartlam-Brooks, H.L.A., Bonyongo, M.C. & 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, 210216.Google Scholar
Barton, K. (2013) MuMIn: Multi-model inference. R v. 1.9.13. Http://CRAN.R-project.org/package=MuMIn [accessed 14 April 2014].Google Scholar
Berger, J., Cain, S.L. & Berger, K.M. (2006) Connecting the dots: an invariant migration corridor links the Holocene to the present. Biology Letters, 2, 528531.Google Scholar
Bolger, D.T., Newmark, W.D., Morrison, T.A. & Doak, D.F. (2008) The need for integrative approaches to understand and conserve migratory ungulates. Ecology Letters, 11, 6377.Google Scholar
Briers, R.A. (2002) Incorporating connectivity into reserve selection procedures. Biological Conservation, 103, 7783.CrossRefGoogle Scholar
Brooks, C.J. & Harris, S. (2008) Directed movement and orientation across a large natural landscape by zebras, Equus burchelli antiquorum . Animal Behaviour, 76, 277285.Google Scholar
Brower, L.P. & Malcolm, S.B. (1991) Animal migrations: endangered phenomena. American Zoologist, 31, 265276.Google Scholar
Bunnefeld, N., Börger, L., van Moorter, B., Rolandsen, C.M., Dettki, H., Solberg, E.J. & Ericsson, G. (2011) A model-driven approach to quantify migration patterns: individual, regional and yearly differences. Journal of Animal Ecology, 80, 466476.Google Scholar
Burnham, K.P. & Anderson, D. (1998) Model Selection and Multimodel Inference: a Practical Information-theoretic Approach. Springer, New York, USA.Google Scholar
Cameron, A.C. & Trivedi, P.K. (1998) Regression Analysis of Count Data. Cambridge University Press, New York, USA.Google Scholar
Caro, T., Darwin, J., Forrester, T., Ledoux-Bloom, C. & Wells, C. (2012) Conservation in the anthropocene. Conservation Biology, 26, 185188.Google Scholar
Chase, M.J. (2011) Status of Wildlife Populations and Land Degradation in Botswana's Forest Reserves and Chobe District. Elephants Without Borders, Kasane, Botswana.Google Scholar
Chase, M. (2013) Status of Wildlife Populations and Land Degradation in Botswana's Forest Reserves and Chobe District. Elephants Without Borders, Kasana, Botswana.Google Scholar
Crooks, K.R. & Sanjayan, M. (eds) (2006) Connectivity Conservation. Cambridge University Press, Cambridge, UK.Google Scholar
Department of Wildlife and National Parks (1995) Makgadikgadi/Nxai Pan Management Plan. Ministry of Commerce and Industry, Gabarone, Botswana.Google Scholar
Dingle, H. & Drake, V.A. (2007) What is migration? BioScience, 57, 113121.Google Scholar
Dobson, A.P., Borner, M., Sinclair, A.R.E., Hudson, P.J., Anderson, T.M., Bigurube, G. et al. (2010) Road will ruin Serengeti. Nature, 467, 272273.Google Scholar
Dormann, C.F., McPherson, J.M., Araujo, M.B., Bivand, R., Bolliger, J., Carl, G. et al. (2007) Methods to account for spatial autocorrelation in the analysis of species distributional data: a review. Ecography, 30, 609628.Google Scholar
Epps, C.W., Mutayoba, B.M., Gwin, L. & Brashares, J.S. (2011) An empirical evaluation of the African elephant as a focal species for connectivity planning in East Africa. Diversity and Distributions, 17, 603612.Google Scholar
Ferguson, K. & Hanks, J. (eds) (2010) Fencing Impacts: A Review of the Environmental, Social and Economic Impacts of Game and Veterinary Fencing in Africa with Particular Reference to the Great Limpopo and Kavango–Zambezi Transfrontier Conservation Areas. Mammal Research Institute, Pretoria, South Africa.Google Scholar
Fryxell, J.M. & Sinclair, A.R.E. (1988) Causes and consequences of migration by large herbivores. Trends in Ecology & Evolution, 3, 237241.Google Scholar
Fynn, R.W.S. & Bonyongo, M.C. (2011) Functional conservation areas and the future of Africa's wildlife. African Journal of Ecology, 49, 175188.Google Scholar
Gibson, D. (2010) Environmental assessment of Botswana's fences. In Fencing Impacts: A Review of the Environmental, Social and Economic Impacts of Game and Veterinary Fencing in Africa with Particular Reference to the Great Limpopo and Kavango-Zambezi Transfrontier Conservation Areas (eds Ferguson, K. & Hanks, J.), pp. 121127. Mammal Research Institute, Pretoria, South Africa.Google Scholar
Hansen, M.C., DeFries, R.S., Townshend, J.R.G., Carroll, M., Dimiceli, C. & Sohlberg, R.A. (2003) Global percent tree cover at a spatial resolution of 500 meters: first results of the MODIS vegetation continuous fields algorithm. Earth Interactions, 7, 115.2.0.CO;2>CrossRefGoogle Scholar
Harris, G., Thirgood, S., Hopcraft, G., Cromsigt, J.P.G.M. & Berger, J. (2009) Global decline in aggregated migrations of large terrestrial mammals. Endangered Species Research, 7, 5576.Google Scholar
Johnson, S. (2010) Chapter 3: Tourism and heritage development. In Makgadikgadi Framework Management Plan. Centre for Applied Research and Department of Environmental Affairs, Government of Botswana, Gabarone, Botswana.Google Scholar
Jones, B.T.B., Davis, A., Diez, L. & Diggle, R.W. (2012) Community-based natural resource management (CBNRM) and reducing poverty in Namibia. In Biodiversity Conservation and Poverty Alleviation: Exploring the Evidence for a Link (eds Roe, D., Elliott, J., Sandbrook, C. & Walpole, M.), pp. 191205. John Wiley & Sons, Ltd., London, UK.Google Scholar
KAZA TFCA (2013) Kavango–Zambezi Transfrontier Conservation Area. Http://www.kavangozambezi.org [accessed 5 April 2014].Google Scholar
Lindsay, K., Muir, A. & Robson, K. (1998) Pilot Botswana Range Inventory and Management Project. Ministry of Agriculture, Gabarone, Botswana.Google Scholar
Lindsey, P.A., Roulet, P.A. & Romanach, S.S. (2007) Economic and conservation significance of the trophy hunting industry in sub-Saharan Africa. Biological Conservation, 134, 455469.Google Scholar
Mbaiwa, J.E. & Mbaiwa, O.I. (2006) The effects of veterinary fences on wildlife populations in the Okavango Delta, Botswana. International Journal of Wilderness, 12, 1741.Google Scholar
Morales, J.M., Haydon, D.T., Frair, J., Holsiner, K.E. & Fryxell, J.M. (2004) Extracting more out of relocation data: building movement models as mixtures of random walks. Ecology, 85, 24362445.Google Scholar
Norton-Griffiths, M. (1978) Counting Animals. African Wildlife Leadership Foundation, Nairobi, Kenya.Google Scholar
Perkins, J. (2010) Fences and landscape scale degradation. In Fencing Impacts: A Review of the Environmental, Social and Economic Impacts of Game and Veterinary Fencing in Africa with Particular Reference to the Great Limpopo and Kavango-Zambezi Transfrontier Conservation Areas (eds Ferguson, K. & Hanks, J.), pp. 108120. Mammal Research Institute, Pretoria, South Africa.Google Scholar
Pettorelli, N., Vik, J.O., Mysterud, A., Gaillard, J.-M., Tucker, C.J. & Stenseth, N.C. (2005) Using the satellite-derived NDVI to assess ecological responses to environmental change. Trends in Ecology & Evolution, 20, 503510.Google Scholar
Power, R.J. & Compion, R.X.S. (2009) Lion predation on elephants in the Savuti, Chobe National Park, Botswana. African Zoology, 44, 3644.CrossRefGoogle Scholar
R Development Core Team (2008) R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria.Google Scholar
Sawyer, H. & Kauffman, M.J. (2011) Stopover ecology of a migratory ungulate. Journal of Animal Ecology, 80, 10781087.Google Scholar
Sianga, K., Bonyongo, M.C. & Fynn, R. (2013) Seasonal Movements of Plains Zebra Equus burchelli in Savuti-Mababe-Linyanti (SLE) Ecosystem in Northern Botswana. University of Botswana, Gaborone, Botswana.Google Scholar
Sinclair, A.R.E. & Norton-Griffiths, M. (eds) (1979) Serengeti: Dynamics of an Ecosystem. University of Chicago Press, Chicago, USA.Google Scholar
Vandewalle, M.E. (2000) Movement of Zebra and Wildebeest in Relation to Environmental Factors in Northern Botswana. University of the Witwatersrand, Johannesburg, South Africa.Google Scholar
Wilcove, D.S. (2009) No Way Home: The Decline of the World's Great Animal Migrations. Island Press, Washington, DC, USA.Google Scholar
Wilcove, D.S. & Wikelski, M. (2008) Going, going, gone: is animal migration disappearing? PLoS Biology, 6(7), e188.Google Scholar
Williamson, D. & Williamson, J. (1984) Botswana's fences and the depletion of Kalahari wildlife. Oryx, 18, 218222.Google Scholar
Figure 0

Fig. 1 Movement trajectories of eight female zebra Equus quagga collared on the Chobe floodplains in Botswana and Namibia from the late dry season (September–October) 2012 until June 2013. Grey polygons indicate actual (Nxai Pan National Park) and potential alternative (Savuti Marsh and Seloko Plains) migration destinations. The rectangle on the inset indicates the location of the main map in southern Africa.

Figure 1

Table 1 Summary statistics for migration movements of eight female zebras Equus quagga collared on the Chobe floodplains in Botswana and Namibia (Fig. 1).

Figure 2

Fig. 2 Timing of daily rainfall and zebra migration during September 2012–May 2013. The thick black line represents the mean net displacement away from dry season home ranges for nine zebras. Vertical lines represent daily rainfall over dry season ranges (black) and the Nxai Pan wet season range (grey).

Figure 3

Fig. 3 Temporal pattern of mean enhanced vegetation index (September 2012–May 2013) at the Nxai Pan wet season migration destination and two alternative destinations (Seloko Plains and Savuti Marsh). The ellipse indicates the near-identical values at Nxai Pan and Seloko Plains during the core wet season time when zebras were at Nxai Pan.

Figure 4

Table 2 Coefficients for regression models relating environmental characteristics to the count of zebra locations in 5 × 5 km cells along the migratory pathway (negative binomial regression model) and to the probability of a cell being a stopover point on the migration route (logistic regression model) during both northward and southward movement phases.

Supplementary material: PDF

Naidoo Supplementary Material

Table S1

Download Naidoo Supplementary Material(PDF)
PDF 329 KB