Benefits of migration in a partially-migratory tropical ungulate
© Gaidet and Lecomte; licensee BioMed Central Ltd. 2013
Received: 14 December 2012
Accepted: 26 September 2013
Published: 30 September 2013
Partial migration, where one portion of a population conducts seasonal migrations while the other remains on a single range, is common in wild ungulate populations. However the relative costs and benefits associated with the distinct strategies adopted by coexisting migrant and resident individuals have rarely been investigated. Here we compare the body condition of migrants and residents in a partially migratory population of impalas (Aepyceros melampus) in Zimbabwe. The study was conducted during two consecutive years with highly contrasted population densities (16.4 and 8.6 indiv/km2) due to harvesting.
We first identify a population substructure with a north–south sub-division in two spatial units related to distinct soils and vegetation cover. Impalas in the north range had a consistently higher diet quality and body condition than those in the south range. At the beginning of the dry season about one third of the individuals migrated from the lower (i.e. south) to the higher (i.e. north) diet quality range. This partial migration pattern was consistent between the consecutive years, and most individuals showed constancy to their moving strategy (migrant or resident). In both years, these migrants had a significantly higher body condition at the end of the dry season than the south residents that remained year-round in the lower diet quality range. Diet quality and body condition of impalas were higher in the year of lower density; however we did not detect any evidence for density-dependence in migration propensity, at the individual or population levels, nor in the benefit associated with migration.
Our findings provide rare evidence for a significant relationship between body condition and seasonal migration strategy in wild ungulates in relation to a difference in the quality of resources acquired between distinct seasonal ranges. This study also constitutes rare evidence of partial migration in a tropical ungulate population.
KeywordsPartial migration Habitat–performance relationships Density dependence Impala Individual variability Spatial heterogeneity Diet quality Savannah
In spatially and seasonally heterogeneous environments, individuals have access to resources of different quality if limited dispersal or prospection abilities, social constraints, attachment to maternal environment or ecological barriers prevent them to be distributed proportionally to resource quality and availability [1–5]. Within populations of free-ranging herbivores, the individual differences in the capacity to access and acquire forage resources contribute to the inter-individual variations in phenotypic quality and performances including lifetime reproductive success [1, 5–8].
Seasonal migration is a large-scale habitat selection strategy frequently adopted by wild ungulates as a response to spatial variation in seasonal resource limitation [9–11]. In many ungulate populations, only a fraction of the population migrates between distinct seasonal ranges while the others remain sedentary year-round (i.e. partial migration ). This difference in strategy to track resources among co-existing migrant and resident individuals may contribute to shape individual variations in performance. Partial migration studies in ungulates have mostly focused so far on the characterisation of the biotic and abiotic differences between seasonal ranges used by residents and migrants to identify the determinants of the migration strategy (e.g. [13–15]). Incentive for migration has been most commonly related to access to higher forage quality [11, 16–20] and, in some cases, to a reduced predation exposure [21, 22] or human disturbance . However, comparative measures of the relative performances of migrant and resident individuals have rarely been conducted to evaluate the costs and benefits associated with distinct migration tactics [17, 18, 23].
In this study we investigated among-individual variations in nutrition and body condition in a partially migratory population of impalas (Aepyceros melampus) in a semi-arid savannah in Zimbabwe, Southern Africa. Nutrition and body condition are major determinants of survival and reproductive success in ungulates [24, 25] and are, therefore, commonly used proxies of individual performance. Impala is a common medium-sized African ungulate with a relatively small home-range (c. 0.5-5 km2[26–29]). This species is highly gregarious but has a typically unstable social group system . We here tested the existence of individual variations in body condition in relation to differences in resource acquisition between impalas that used different seasonal ranges and different migration strategies. We evaluated differences in the quality of food intake using two conventional faecal indicators (i.e. the nitrogen and the fibre concentration in the faeces) related to the crude protein content and the digestibility of the diet [31, 32].
Rainfall is highly seasonal and annually highly variable (mean: 770 mm/year, min.-max.:355-1,300 mm/year over 1981-2003, Zimbabwean Department of Meteorological Services). We considered two major seasons: a wet season from November to May corresponding to the period of major rainfall and vegetation growth, and a dry season from June to October corresponding to the period of vegetation senescence and restricted water surface availability. During the wet season drinking water was available from various water sources (pond, water points, river pools) spread over the entire ranch that filled up with local rainfall. The availability of drinking water at all these water points was monitored every month from May to October. At the end of the dry season (i.e., in October) surface water was restricted to only three and two water points in 2002 and 2003 respectively (the water point located on the river bed being dry in October 2003) (Figure 1). Predation was not measured directly, but the presence of potential predators (mainly leopard Panthera pardus, caracal Felis caracal and side-striped jackal Canis adustus) was assessed during the study period through spoor detection and direct observation on the network of dust roads (60 km over the entire ranch) continuously covered for management and field research activities (at least twice a week) .
A small resident population of about 100 impalas was originally present in the study area as estimated from line-transect counts . A total of 402 impalas were introduced in 1999 and 2000. All introduced individuals were marked with a numbered ear tag and were released after a few weeks of captivity in a boma (acclimatization enclosure) situated in the south of the ranch for adaptating to their new environment (Figure 1). The study was conducted from January 2002 to October 2003, by when all marked impalas were adults (≥ three years old). The population density of impalas estimated at the end of the dry season (October) from the proportion of marked animals observed at water points decreased from 16.4 ind./km2 in 2002 to 8.6 ind./km2 in 2003 (Figure 1) due to a combination of hunting, poaching and potential emigration .
Individual seasonal ranges
The location of marked impalas was monitored through direct observation every month (5 days per month on average) from January 2002 to October 2003. The entire study area was systemically prospected by bicycle along the dense road network to locate marked impalas. Location coordinates were recorded using a GPS and mapped through a GIS. Individuals were identified by their tag number using a telescope (20-60 × 80). The utilisation a telescope with a tripod for reading individual tags requires animals to stand still. The study population in this game ranch was hunted and impalas were shy. Hunting was conducted throughout the study area from a car at night (using a spot light) along the road network. As a result, impalas routinely flee when approached by the observer and the tag ID of marked impalas observed during the prospection was often not identified. The relatively low number of locations collected for each individual per season prevented us from using conventional home range metrics to characterise individual distribution patterns. Previous studies of the social structure of impalas have identified a clear sub-division of the population into stable clans but with unstable social groups characterised by frequent fission and fusion events [28, 30]. These studies revealed that the centre of activity of impalas from the same clan are strongly clustered, with an extensive overlap (c. 75% on average) between home ranges of clan members and a low overlap (c. 4%) between individuals from adjacent clans .
We here used a spatial cluster analysis method to characterise the discontinuous distribution of the study population and identify the subdivision units. We used a hierarchical clustering with the Ward’s agglomeration method using Euclidian distance between location coordinates of all marked impalas (‘stats’ R package, ‘hclust’ function). This sums-of-squares agglomerative clustering is based on the minimisation of within-cluster versus between-cluster variance. The distance between the nodes of the hierarchical tree was used to identify the relevant number of spatial clusters. The cluster analysis was conducted on all the locations collected during the study period (2002 and 2003) (Figure 2). Note that observations collected at water points during the body condition monitoring (see below) were not included in the spatial cluster analysis to avoid bias related to the attractive effect of water points during the dry season.
We analysed the consistency of assignment to a spatial cluster for each individual within and between seasons to evaluate the spatial fidelity to a seasonal range. A seasonal range could be attributed to only a portion of individuals from the marked population in each season and in each year since our monitoring system based on opportunistic observations did not systematically provide locations for all individuals. An impala was defined as resident when it was located in the same range in both wet and dry seasons and was defined as migrant when it used a distinct seasonal range. We categorised the direction of migration as the bearing from the wet season range to the dry season range.
We analysed the variation in the proportion of migrants between years and between males and females using a Generalized Linear Mixed Model (GLMM) with individual migration status as a binary-dependent variable (migrant vs. resident) and year and sex as fixed effects. The model was fitted with a binomial error distribution and a logit link function (lme4 R package, glmer procedure). A migration status could be attributed in the two consecutive years for about half the individuals that we monitored. To account for the effect of pseudo-replication associated with repeated measures from the same individual we included individual identity as a random effect in models. We used a likelihood ratio test (LRT) to test the significance of each variable in the model, computing the χ2 of the LRT between the model retaining and the model excluding the variable. We used the Pearson χ2-test to assess the goodness-of-fit of the selected models.
We evaluated the diet nutritional quality of impalas through the analysis of the nitrogen and acid detergent fibre concentrations in the faeces. Faecal nitrogen (FN) is a widely used indicator of diet quality in ungulates reflecting (at least broadly) crude protein content of the diet . Acid detergent fibre (ADF) that reflects mostly the ligno-cellulose fraction of the diet is inversely related to forage digestibility hence nutrient accessibility .
Faecal droppings of impalas were collected over a three week-period on three occasions each year: in the wet season (late March - early April), and at the beginning (July) and end (October) of the dry season. Samples were collected at various sites distributed over the range where marked impalas had been located. Coordinates of the sampling sites were recorded using a GPS. An average of 75 samples of fresh faeces was collected per sampling occasion. Each sample consisted of 10 to 25 pellets collected from an isolated dropping considered to represent a distinct individual. Morphological characteristics of faecal pellets permit droppings of impalas to be distinguished from those of other ungulate species present in the study area. Just after collection, faecal samples were dried at 60°C for 48 hours and mill-ground into a homogeneous powder.
Faecal samples (n = 447) were analysed separately for FN and ADF concentrations using the near infrared reflectance spectroscopy (NIRS), a cost-effective method increasingly used to measure nutritional content of free-ranging herbivores in particular ruminants . NIRS estimation of FN and ADF concentrations relies on a specific calibration equation between spectral and chemical measurement from a subset of reference samples. Near infrared absorbance spectra were collected for all samples on a monochromator (NIRS 6500, Foss NIRSystems, Laurel, MD, USA) every 2 nm in the 1102 to 2498 nm range. A representative subset of 20% of the samples (n = 89) were selected among the data set based on the Mahalanobis distances between samples. The selection procedure was performed using the SELECT procedure in WINISI software (WINISI III, Infrasoft International State College, PA, USA). These reference samples were analysed for FN concentration using the Dumas elementary analysis and for ADF concentration using the Dosi Fibre method . Concentrations were expressed in percentage of faeces dry mass. A multi-linear regression model was established for each constituent between the spectral absorbency data and the chemical values of the reference samples, according to a partial least squares regression procedure. The coefficients of determination of the calibration equation indicate a good predictive capacity of the model (R2 = 0.94 and 0.95 for FN and ADF respectively; standard errors of calibration (SEC) = 0.61 and 1.68 for FN and ADF respectively) . The model performance was evaluated with cross validation procedures on 6 groups. The standard errors of cross validation (SECV) indicate a high degree of precision in the estimates of the concentration of both constituents (SECV = 0.78 and 2.37 for FN and ADF respectively), with highly significant coefficients of determination of the cross validation (mean R2 = 0.90 and 0.90 for FN and ADF respectively). Finally the performance to deviation ratio (RPDcv = SD/SECV) was satisfactory according to the criteria of  (3.2 and 3.1 for FN and ADF respectively). The predictive models built on reference samples were used to estimate FN and ADF for all faecal samples collected.
An analysis of variance (Anova) was used to test for the effects of year, season and range (i.e., the north and south population spatial units identified in the previous spatial cluster analysis) on the concentration of FN and ADF. A potential difference in the effect of range between years or between seasons was also investigated (two-way interactions between range and year and between range and season). A backward stepwise selection procedure was used with successive removals of non-significant interactions and factors. Residuals were tested for normality using the Shapiro-Wilk test.
The body condition of impala was estimated through visual assessment at the end of the dry season when animals were at their lowest annual condition. Body condition was estimated by the same observer (NG) and at the same period on both years (second and third week of October) for comparative purposes. Observations were conducted from a blind or a tree-platform established near water points, offering prolonged and clear observation of animals in open area and at a close range (30-50 meters). A preliminary calibration session was conducted in September 2002 to test the visual scoring criteria on the impala population. Care was taken to classify animals only when observed from the side and standing on flat ground, since sunlight incidence and muscle solicitation while bending on sand ridges for drinking could affect body appearance, hence visual scoring . All the water points (three in 2002 and two in 2003) were successively monitored from 6 am to 6 pm for two to four consecutive days. Visual Body Condition (VBC) was assessed in two categories using criteria of progressive depletion of fat deposit on the animal’s hindquarters adapted from . Impalas were first identified from their ear tag number and then scored as “fair” (roundness or slightly angular appearance of rump and flank) or “poor” (side view of the tail, strongly angular to extremely prominent appearance of the lumbar vertebrae or the exterior line of the ischium).
The repeatability of VBC scoring was evaluated using the Cohen’s Kappa statistic for individuals that were independently re-scored during consecutive days in each year. GLMMs were also used to analyse individual variation in body condition, with VBC score as a binary-dependent variable (poor vs. fair). About one third of impalas were scored for body condition in both years, therefore the individual identity was included as a random factor in models to account for the repeated measures. Our analysis consisted of two steps. We first tested for the effect of the dry season range (north versus south), year and sex on body condition. We tested the interaction between range and year to evaluate the potential change in the between-range differences between the two consecutive years of high and low density. In a second step we investigated whether individual variation in body condition was related to their migration strategies. We similarly tested for the effect of the migration strategy, year, sex and the interaction between migration strategy and year on body condition. For all models a likelihood ratio test (LRT) was used to test the significance of each variable in the model. We used the Pearson χ2-test to assess the goodness-of-fit of the selected models.
Seasonal range and migration patterns
The prospection of the study area provided a total of 480 locations from 171 individual impalas (144 females and 27 males) during the first year and 618 locations from 165 individuals (128 females and 37 males) during the second year. The cluster analysis of the spatial distribution of all these locations revealed a clear partition into two spatial units, with a drastic drop in the cluster criterion value found after the two-cluster level (Figure 2a). The distribution pattern of the two spatial units showed a north–south subdivision in both wet and dry seasons (Figure 2b). The range of the north unit was dominated by a mopane woodland and the range of the south unit was dominated by a mixed woodland (Figure 2c). The presence of predators (leopard, caracal and side-striped jackal) was detected by spoor or directs observation on both spatial units. These two spatial units were located on both sides of the river bed that cross the study site, separated by an area of mixed shrubland where impalas (both marked and unmarked individuals) were rarely observed. The same north–south subdivision was observed in both years and was therefore retained to investigate individual distribution patterns. In both years, a permanent water point was available in both the north and the south spatial units (Figure 1).
Each individual was located between one to ten times per season. Individuals were seen on average on a similar number of occasions in the north range (mean = 1.94 locations per individual per season) and in the south range (mean = 2.19 locations). All the locations from an individual recorded within a given season were generally located in the same spatial unit. The consistency of individual assignment to the north or the south range within a given season was confirmed for almost all impalas (96%, n = 55 individuals) for which five to ten locations could be collected during a single season. Only two females were observed concurrently in the two distinct ranges within the same season. This high consistency in the location of individuals at the range scale justifies the attribution of a seasonal range for individuals that could be observed on only one occasion per season.
We also found a good agreement between the individual assignment to a spatial unit and the water point where these individuals were observed during the body condition monitoring. In both years all observations of marked impalas recorded at the water point located in the north of the study site (n = 93 individuals) corresponded to individuals that had been independently classified in the north-range unit for the dry season. Likewise all observations recorded at the water point located in the south (n = 56 individuals) correspond to individuals that had been classified in the south-range unit. In 2002 most observations (n = 228 individuals) were collected at the water point located on the river bed (Figure 1) where impalas from both the north and the south ranges were observed.
Our analysis of the variation in the proportion of migrants using a GLMM method revealed that most of the variance in migration status was related to differences between individuals, with a strong effect of repeated measures for individuals monitored in both years (variance estimation of the individual id. random effect = 74). Most of these individuals show a constant migration strategy over the two years: 45 out of 51 of impalas classified as resident in the first year remained resident in the second year, and 9 out of 10 impalas remained migrants (constant south to north migration). Because of lack of convergence of the GLMM (potentially due to the large number of levels within the random factor and the low number of repetitions, maximum 2 years, within each level), we re-ran the analysis using a fixed generalized linear model (GLM) on a partial data set where individuals that had been monitored in both years were randomly selected from one year only. Contrary to our prediction, we found no significant difference in the proportion of migrants between 2002 (31.6%) and 2003 (40.3%) (LRT: χ2 = 1.63, p = 0.20; difference of 0.54 ± 0.42 on a logit scale) (Figure 4). Further, no difference in the proportion of migrants was found between females (35.8%) and males (32.7%) (LRT: χ2 = 0.12, p = 0.73; difference of 0.21 ± 0.64 on a logit scale) (null model: AIC = 135.31; Pearson χ2 goodness-of-fit test = 103.0, ddl = 102, p = 0.45).
We classified the faeces sampling sites from their coordinates into the north and south ranges. We found a significant difference in nitrogen and fibre concentrations between the two ranges (Range: F = 4.45, p < 0.001 for FN; F = 32.32, p < 0.001 for ADF). The diet quality was higher in the north than in the south range in all seasons and in both years (interactions between range and year or between range and season were not significant) (Figure 5). All these relationships remained significant after removing outliers (mean ± 2 standard deviations) to satisfy the assumption of normality of residuals (Shapiro–Wilk, W = 0.99, p = 0.16 for FN; W = 0.99, p = 0.48 for ADF).
Variation in body condition
Body condition was assessed for a total of 153 individuals (126 females and 27 males) during the first dry season (2002) and for 63 individuals (58 females and 5 males) during the second dry season (2003). Repeatability of individual VBC scoring in a given dry season was relatively high: 79% of impalas (n = 126) that had been independently scored during consecutive days at the water points in the same year were consistently assigned to the same VBC class, and the Kappa statistics indicate a substantial agreement beyond chance between repeated scoring (Kappa test = 0.60, p < 0.001). The remaining individuals that had been assigned to different VBC classes in a given dry season were classified to their most frequent VBC class or (in case of paired number of VBC scores) were arbitrarily assigned to the class attributed at the first occasion. Removing these individuals (i.e., with a non steady VBC scores) from the dataset did not change the relationship between individual variation in body condition and explanatory variables tested in the following analyses, so we only reported results from the full dataset.
In our second analysis we found a significant relationship between migration strategy and body condition (LRT: χ2 = 13.47, p = 0.001) in addition to the between-year difference in body condition (LRT: χ2 = 18.29, p < 0.001) (Figure 6). Impalas that migrated from the south to the north range at the beginning of the dry season (i.e. S.-N. migrants) had a higher body condition than those that remained in the south range (i.e. south residents) (difference of 2.17 ± 0.68 on a logit scale, p < 0.001). No significant difference was found in body condition between the south to north migrants and the north residents (difference of 0.44 ± 0.62 on a logit scale, p = 0.48). The same pattern was observed in both years with no significant interaction between year and migration strategy (LRT: χ2 = 0.79, p = 0.67). The single individual that migrated from the north to the south range at the beginning of the dry season (i.e. N.-S. migrant) for which a body condition score had been attributed was classified as poor.
Our results provide rare evidence for a significant relationship between seasonal migration strategy and body condition in wild ungulates in relation to a difference in the quality of resources acquired between distinct seasonal ranges. The diet quality and body condition of impalas present in the north range was consistently higher than in impalas present in the south range. Individuals were consistently observed in one of the two ranges within a season but about one third of the individuals migrated from the lower (i.e. south) to the higher (i.e. north) diet quality range at the beginning of the dry season. This partial migration pattern was consistent between two consecutive years, and most individuals showed constancy to their moving strategy (migrant or resident). In both years migrants had a significantly higher body condition at the end of the dry season than individuals that remained year-round in the lower diet quality range (i.e. south resident).
The diet quality and the body condition of impalas were higher in the year of lower density suggesting a concomitant density dependent response in the nutritional and body condition of impalas. However, we did not detect any evidence for an influence of density on migration propensity and on the benefit associated with migration. First, we found no between-year difference in the proportion of migrants in the population and we observed a relatively stable migration strategy at the individual scale. Second, we did not detect a more contrasted spatial variation in diet quality or body condition between impalas that used distinct seasonal ranges or migration strategies in the year of higher density. This suggests that there is no density-dependence in the differences of body condition among individuals using these different space use tactics.
We acknowledge that the monitoring of individual space use through direct-visual identification of marked animals instead of telemetry (using radio-tracking or GPS technology) greatly limited the number of locations that we could collect and that this may be a limitation in our study. Though we cannot completely rule out the possibility for animals to have used a range elsewhere, the limited number of locations per individual was more likely due to the difficulty of identifying marked animals rather than a failure to detect marked animals during the monitoring. The entire study area was prospected to locate marked impalas but all marked impalas observed during the prospection could not be identified owing to long fleeing distance and the density of the vegetation cover. The impala population was hunted in this game ranch (by car at night throughout the study area) therefore routinely flee when approached by the observer during our monitoring.
Our data suggests that individuals were faithful between years to a seasonal range and to a migration strategy. Cautious interpretation, however, is needed because of the small number of locations collected for each individuals and the underlying assumption that they are faithful to a range during a given season. The consistency of relocation within the same seasonal range was high (96%) among individuals the most frequently observed (n = 5 to 10 locations per season), including both residents and migrants. Moreover, the attribution of a spatial tactic (from locations collected during the prospection over the study area) and a body condition score (from observation during the water point monitoring) was conducted independently through two separate protocols. Though errors in the attribution to a seasonal range or a migration strategy cannot be completely ruled out, this bias should have affected individuals independently to their body condition score, and should not have affected (apart from the introduced noise) the relationship between individual spatial tactic and body condition.
Spatio-temporal variations in diet quality and body condition
The subdivision of the study population in two spatial units that we found is consistent with the discontinuous  or the clumped distribution  of other impala populations previously described in patchy habitats. Our study area lies on sedimentary sandstone foundations that have led to the formation of soils of varied quality  and distinct vegetation covers (, Figure 2). These two ranges are separated by a river bed (Figure 2) that though is dry most of the year may act as a natural frontier between impalas’ ranges.
Both faecal indicators that we monitored indicated a typical decrease in the diet quality of impalas during the dry season  as the digestibility and the nutritive value of forage decreased with the senescence and the progressive depletion of forage resources. The significantly lower diet quality that we found during the year of higher population density, together with the concurrent between-year difference in body condition and recruitment previously reported , provides evidence for a resource-limitation control of the performance of our study population. Forage limitation has been similarly proposed as the main regulating factor of most of the African ungulate populations that have been studied so far (white-eared kob Kobus kob; greater kudu Tragelaphus strepsiceros; wildebeest Connochaetes taurinus).
A growing number of empirical studies have evidenced a relationship between the spatial heterogeneity in resource quality and the variation in performance among individuals or sub-population units (see  for a thorough review). In our study we found a positive co-variation in the diet quality and the body condition between the north and the south ranges, suggesting a spatial variation in habitat quality. The higher diet quality and body condition measured in impala present in the north range, in both years and in all seasons, is consistent with the spatial variation in soil quality and associated vegetation. The north range lies on alluvial terraces that form rich and well structured medium to deep profile lime soils ; it is mostly covered by mopane woodland, a tall, low-density but high-cover woodland with a species-rich herbaceous stratum dominated by annual grasses . The south range consists of brown semi-alluvial sandy soils of lower quality; it is mostly covered by a mixed woodland with a high under-storey cover and lower species richness dominated by perennial grasses.
Partial migration in impalas
Partial migration has been commonly reported in temperate or boreal ungulate populations (e.g. [13, 15–17, 20, 21, 23]) but it has only recently been described in a few tropical ungulates (African buffalo Syncerus caffer; wildebeest ). Our findings constitute the first description to our knowledge of a partial migration pattern in an impala population.
Impala is not generally considered as a migratory species  but some seasonal movements of impalas between neighbouring habitats have been described. The patterns of movements were either a small seasonal shift (0.5 km) in the home range centre (Sengwa Wildlife Research Area, Zimbabwe, ), a periodic abandonment of the usual home range (with movements of up to 10 km) during periods of severe food shortage (Serengeti NP, Tanzania, ), or more extensive movements (5-50 km) leading to distant wet and dry season ranges (Southern low-veld and Mana Pools NP, Zimbabwe [50, 51]). In our study about one third of the impala population shifted between two adjacent ranges and this pattern was consistent during two consecutive years. Migration in ungulates is usually viewed as long-distance movements. However, following Ball et al. (2001) , we considered that short-distance seasonal movements are properly termed ‘migration’ if there is a clear shift between spatially disjunct seasonal ranges, regardless of how close these might be.
Such a partial migration has important consequences for the degree of mixing within a population that has not been reported for impala so far. Though instability in group membership has long been recognised in impala , the social structure has been reported to be highly cohesive year-round at the clan level . A high level of association has been found between members of the same clan with only occasional transfer (2% per annum ) of females between adjacent clans. Conversely in our study about one third of individuals (both females and males) changed their seasonal range. Interestingly we observed no segregation between resident and migrant impalas, which were found in mixed herds while sharing the same seasonal range similar to that described in another gregarious ungulate species, the elk Cervus elaphus. Social constraints to individual exchanges between population spatial units have been reported in various gregarious ungulate species [1, 15]. In impala, the female social status as well as the mother-offspring bond after weaning is of little significance and agonistic interactions between females are scarce . In Southern Africa, most males are only territorial for a few weeks during the time of the annual rut (i.e. in April) and generally associate with females in mixed herd throughout most of the year . The seasonal social mixing that we observed may have been facilitated by the typical unstable group units of impalas with a fluid membership and by the absence of dominance hierarchy . The period of rut did not coincide with the timing of migration (migrants were observed for the first time in their dry season range in July). This suggests that migration was not associated with a mating dispersal in this population.
Several hypotheses have been proposed to explain partial migration in animals ([12, 15, 52]). In our study the direction of migration at the onset of the dry season was correlated with the difference in diet quality between the two ranges for almost all migrants. Access to higher quality forage is generally considered to be one of the major driving forces of seasonal migration [9–12]. At the beginning of the wet season most migrants however returned to their previous wet season range in the south associated with a relatively lower diet quality. The differences in diet quality between seasons were more contrasted than the differences between ranges. In tropical ecosystems food resource acquisition is likely to be much less constraining during the wet season when forage is of higher quality and availability. The nutritional cost associated with the relatively lower forage quality of the south range should be lower in the wet season and may be balanced by other benefits. The high fidelity to a summer or winter range documented in temperate ungulates (e.g. ) or to a wet season range in a tropical ungulate (e.g. ) is considered to provide benefits related to familiarity with the location of resources, predators or conspecifics . The impalas that we monitored had all been introduced at the same release site located in the south range (Figure 1) which therefore may be more familiar than the north range and may provide some predator-avoidance benefits. This is supported by the fact that migrants were first observed in the south range in November just before the calving period (occurring in December), a time when females may be particularly sensitive to the risk of predation . It also has to be noted that the migration movements performed by these impalas are relatively small (a few kilometres). This decreases the risks and associated costs of migration movements (increased energy expenditure or predation exposure along migration pathways) and probably contributes to the benefit of migration observed in this population.
Translocation has been shown to alter the movement patterns of wild ungulates only during the first years after introduction into a novel environment . In our study, monitoring started more than two years after the last impalas had been introduced in the study area. The fact that most individuals used the same seasonal range and displayed the same migration strategy for two consecutive years suggests that the translocated impalas had adopted a stable space-use pattern and that the shift in seasonal range that we observed were actual migrations rather than dispersive or exploratory movements.
Access to surface water has been proposed as a potential determinant of migration in tropical ungulates . Impala is generally considered to be highly water-dependent  and to be forced to make daily excursions for water during the dry season . In our study, despite the drastic reduction in the distribution of water availability between the wet and the dry seasons (see the Methods), we observed no marked changed between the wet and the dry seasons in the distribution of impalas over the north or the south ranges (Figure 3). Access to drinking water did not appear to have shaped migratory behaviour since permanent water sources were equally available in both the north and the south ranges, and impalas were always located at <3 km of a water source. Moreover our continuous monitoring of marked impalas at water points during consecutive days indicates that, at the end of the dry season, impalas were visiting water points only every two or three days rather than daily .
In our study, migration was associated with a benefit in terms of access to a range providing, on average, a higher diet quality and higher performance in terms of body condition. Interestingly, the increase in forage competition in the north range associated with the influx of migrants during the dry season did not offset the advantages of access to forage of higher quality as revealed by the higher diet quality measured in the north range. Migrants that had spent the wet season in the lower quality south range had a similar body condition at the end of the dry season to the individuals that remained year-round in the higher-quality range (i.e. north resident). This suggests that migrants did not experience a seasonal carry-over effect, i.e. an influence of the habitat quality in the preceding season . This result is consistent with the low potential of impalas to store fat reserves . Alternatively, there might have been some individual differences in resource acquisition during the wet season among impalas present in the south range, and only individuals of higher phenotypic quality may have decided to migrate. Investigating the relationship between the individual physiologic, genetic or social characteristics and their migration propensity would be a relevant direction for future researches.
Few studies have explored so far the relationship between migration strategies and fitness-related traits in ungulates [12, 55]. In two recent studies, migratory females had either a higher calving success than residents (in reindeers Rangifer tarandus) or a higher fecundity and calf weight (in elk Cervus elaphus). In contrast, a higher mortality rate (in elk ; and in mule deer Odocoileus hemionus) or lower pregnancy and recruitment rates (in elk ) were found in migratory compared to resident females. High mortality risk during the migratory route, associated with energy travelling expenditure, injuries or higher predation exposure, have been documented in a few other cases (in elk ; in reindeers ). Such a distinct trade-off between survival and reproduction between co-existing migrants and residents, if it is confirmed for other ungulates, would be crucial for shaping the population dynamics of partially migratory ungulate populations  and for predicting their sensitivity to migration disruption [23, 55]. The individual variation in space use tactics increases the plasticity of ungulate population to changes in environmental conditions. An important aspect of future research, with high relevance to conservation of migratory ungulates, would be to examine the individual flexibility to switch between migration tactics in response to change in habitat quality, predation or human disturbance.
This survey was realised within the framework of the Biodiversity Project, Zimbabwe, funded by the FFEM (Fonds Français pour l’Environnement Mondial). We are grateful to C. Muzeza - manager of the Chivaraidze Game Ranch and his collaborators for permission and assistance to work at the ranch. We would like to thank S. Le Doze and B. Butete for assistance in field activities, D. Bastianelli and L. Bonnal for chemical and SPIR analysis, T. Dodman for the English revision and M. Gély for assistance in map production. We thank Professor J. M. Fryxell and Dr. M. Hewison, as well as three other anonymous reviewers, for their insightful suggestions for improvement. This work benefited from comments of Professor E. J. Milner-Gulland, Dr. J. Clobert, Dr. H. Fritz and Dr. J.-M. Gaillard on an earlier draft of this manuscript.
- Conradt L, Clutton-Brock TH, Guinness FE: The relationship between habitat choice and lifetime reproductive success in female red deer. Oecologia. 1999, 120: 218-224. 10.1007/s004420050851.View ArticleGoogle Scholar
- Pettorelli N, Gaillard JM, Duncan P, Maillard D, Van Laere G, Delorme D: Age and density modify the effects of habitat quality on survival and movements of roe deer. Ecology. 2003, 84: 3307-3316. 10.1890/02-0602.View ArticleGoogle Scholar
- Nilsen EB, Linnell JDC, Andersen R: Individual access to preferred habitat affects fitness components in female roe deer. J Anim Ecol. 2004, 73: 44-50. 10.1111/j.1365-2656.2004.00790.x.View ArticleGoogle Scholar
- Jones OR, Pilkington JG, Crawley MJ: Distribution of a naturally fluctuating ungulate population among heterogeneous plant communities: ideal and free?. J Anim Ecol. 2006, 75: 1387-1392. 10.1111/j.1365-2656.2006.01163.x.View ArticlePubMedGoogle Scholar
- McLoughlin PD, Gaillard JM, Boyce MS, Bonenfant C, Duncan P, Delorme D, Van Moorter B, Messier F, Saïd S, Klein F: Lifetime reproductive success and composition of the home range in a large herbivore. Ecology. 2007, 88: 3192-3201. 10.1890/06-1974.1.View ArticlePubMedGoogle Scholar
- Coulson T, Albon SD, Guinness FE, Pemberton J, Clutton-Brock TH: Population substructure, local density, and calf winter survival in red deer. Ecology. 1997, 78: 852-863.View ArticleGoogle Scholar
- Pettorelli N, Gaillard JM, Yoccoz NG, Duncan P, Maillard D, Delorme D, Van Laere G, Toïgo C: The response of fawn survival to changes in habitat quality varies according to cohort quality and spatial scale. J Anim Ecol. 2005, 74: 972-981. 10.1111/j.1365-2656.2005.00988.x.View ArticleGoogle Scholar
- McLoughlin PD, Boyce MS, Coulson T, Clutton-Brock T: Lifetime reproductive success and density-dependent, multi-variable resource selection. Proc R Soc B. 2006, 273: 1449-1454. 10.1098/rspb.2006.3486.PubMed CentralView ArticlePubMedGoogle Scholar
- Festa-Bianchet M: Seasonal range selection in bighorn sheep: conflicts between forage quality, forage quantity, and predator avoidance. Oecologia. 1988, 75: 580-586. 10.1007/BF00776423.View ArticleGoogle Scholar
- Fryxell JM, Sinclair ARE: Causes and consequences of migration by large herbivores. Trends Ecol Evol. 1988, 3: 237-241. 10.1016/0169-5347(88)90166-8.View ArticlePubMedGoogle Scholar
- Albon SD, Langvatn R: Plant phenology and the benefits of migration in a temperate ungulate. Oikos. 1992, 65: 502-513. 10.2307/3545568.View ArticleGoogle Scholar
- Chapman BB, Brönmark C, Nilsson JÅ, Hansson LA: The ecology and evolution of partial migration. Oikos. 2011, 120: 1764-1775. 10.1111/j.1600-0706.2011.20131.x.View ArticleGoogle Scholar
- Ball JP, Nordengren C, Wallin K: Partial migration by large ungulates: characteristics of seasonal moose Alces alces ranges in northern Sweden. Wildlife Biol. 2001, 7: 39-47.Google Scholar
- Cagnacci F, Focardi S, Heurich M, Stache A, Hewison AJM, Morellet N, Kjellander P, Linnell JDC, Mysterud A, Neteler M, Delucchi L, Ossi F, Urbano F: Partial migration in roe deer: migratory and resident tactics are end points of a behavioural gradient determined by ecological factors. Oikos. 2011, 120: 1790-1802. 10.1111/j.1600-0706.2011.19441.x.View ArticleGoogle Scholar
- Mysterud A, Loe LE, Zimmermann B, Bischof R, Veiberg V, Meisingseth EL: Partial migration in expanding red deer populations at northern latitudes - a role for density dependence?. Oikos. 2011, 120: 1817-1825. 10.1111/j.1600-0706.2011.19439.x.View ArticleGoogle Scholar
- Nicholson MC, Bowyer RT, Kie JG: Habitat selection and survival of mule deer: tradeoffs associated with migration. J Mammal. 1997, 78: 483-504. 10.2307/1382900.View ArticleGoogle Scholar
- Hansen BB, Aanes R, Sæther BE: Partial seasonal migration in high-arctic Svalbard reindeer (Rangifer tarandus platyrhynchus). Can J Zool. 2010, 88: 1202-1209. 10.1139/Z10-086.View ArticleGoogle Scholar
- Hebblewhite M, Merrill EH: Demographic balancing of migrant and resident elk in a partially migratory population through forage–predation tradeoffs. Oikos. 2011, 120: 1860-1870. 10.1111/j.1600-0706.2011.19436.x.View ArticleGoogle Scholar
- Naidoo R, et al: Home on the range: factors explaining partial migration of african buffalo in a tropical environment. PLoS ONE. 2012, 7: e36527-10.1371/journal.pone.0036527.PubMed CentralView ArticlePubMedGoogle Scholar
- Bischof R, Loe LE, Meisingset EL, Zimmermann B, Van Moorter B, Mysterud A: A migratory northern ungulate in the pursuit of spring: jumping or surfing the green wave?. Am Naturalist. 2012, 180: 407-424. 10.1086/667590.View ArticleGoogle Scholar
- Hebblewhite M, Merrill EH: Trade-offs between predation risk and forage differ between migrant strategies in a migratory ungulate. Ecology. 2009, 90: 3445-3454. 10.1890/08-2090.1.View ArticlePubMedGoogle Scholar
- Robinson BG, Hebblewhite M, Merrill EH: Are migrant and resident elk (Cervus elaphus) exposed to similar forage and predation risk on their sympatric winter range?. Oecologia. 2010, 164: 265-275. 10.1007/s00442-010-1620-6.View ArticlePubMedGoogle Scholar
- Middleton AD, Kauffman MJ, McWhirter DE, Cook JG, Cook RC, Nelson AA, Jimenez MD, Klaver RW: Animal migration amid shifting patterns of phenology and predation: lessons from a Yellowstone elk herd. Ecology. 2013, 94: 1245-1256. 10.1890/11-2298.1.View ArticlePubMedGoogle Scholar
- Gaillard JM, Festa-Bianchet M, Yoccoz NG, Loison A, Toïgo C: Temporal variation in fitness components and population dynamics of large herbivores. Annu Rev Ecol Syst. 2000, 31: 367-393. 10.1146/annurev.ecolsys.31.1.367.View ArticleGoogle Scholar
- Parker KL, Barboza PS, Gillingham MP: Nutrition integrates environmental responses of ungulates. Funct Ecol. 2009, 23: 57-69. 10.1111/j.1365-2435.2009.01528.x.View ArticleGoogle Scholar
- Jarman MV: Attachment to home area in impala. E Afr Wildlife J. 1970, 8: 198-200. 10.1111/j.1365-2028.1970.tb00841.x.View ArticleGoogle Scholar
- Leuthold W: Observations on the social organisation of Impala (Aepyceros melampus). Z Tierpsychol. 1970, 27: 693-721.View ArticleGoogle Scholar
- Murray MG: Home range dispersal and the clan system of impala. Afr J Ecol. 1982, 20: 253-269. 10.1111/j.1365-2028.1982.tb00301.x.View ArticleGoogle Scholar
- du Toit JT: Home range - body mass relations: a field study on African browsing ruminant. Oecologia. 1990, 85: 301-303. 10.1007/BF00319416.View ArticleGoogle Scholar
- Murray MG: Structure of association in impala Aepyceros melampus. Behav Ecol Sociobiol. 1981, 9: 23-33. 10.1007/BF00299849.View ArticleGoogle Scholar
- van Soest PJ: Nutritional Ecology of the Ruminant. 1982, Corvallis Oregon: O and B BooksGoogle Scholar
- Leslie DM, Bowyer RT, Jenks J: Facts from feces: Nitrogen still measures up as a nutritional index for mammalian herbivores. J Wildlife Manage. 2008, 72: 1420-1433. 10.2193/2007-404.View ArticleGoogle Scholar
- Gaidet N, Gaillard JM: Density-dependent body condition and recruitment in a tropical ungulate. Can J Zool. 2008, 86: 24-32. 10.1139/Z07-111.View ArticleGoogle Scholar
- Poilecot P, Gaidet N: A quantitative study of the grass and woody layers of a mopane (Colophospermum mopane) savannah in the mid-Zambezi valley, Zimbabwe. Afr J Ecol. 2010, 49: 150-164.View ArticleGoogle Scholar
- De Wispelaere G: Natural Resource Management within Multispecies Systems in the Mid-Zambezi Valley, Zimbabwe. Report No. 01-32. 2001, Montpellier, France: Commission Européenne/CIRADGoogle Scholar
- Gaidet N, Fritz H, Nyahuma C: A participatory counting method to monitor populations of large mammals in non-protected areas: a case study of bicycle counts in the Zambezi Valley, Zimbabwe. Biodivers Conserv. 2003, 12: 1571-1585. 10.1023/A:1023646012700.View ArticleGoogle Scholar
- Dixon RM, Coates DB: Near infrared spectroscopy of faeces to evaluate the nutrition and physiology of herbivores. J Near Infrared Spectrosc. 2009, 17: 1-31. 10.1255/jnirs.822.View ArticleGoogle Scholar
- Williams PC: Implementation of near-infrared technology. Near-infrared Technology in the Agricultural and Food Industries. Edited by: Williams P, Norris K. 2001, Minnesota: American Association of Cereal Chemists, 145-169.Google Scholar
- Gasaway WC, Gasaway KT, Berry HH: Persistent low densities of plains ungulates in Etosha National Park, Namibia: testing the food-regulating hypothesis. Can J Zool. 1996, 74: 1556-1572. 10.1139/z96-170.View ArticleGoogle Scholar
- Riney T: A field technique for assessing physical condition in some ungulates. J Wildlife Manage. 1960, 24: 92-94. 10.2307/3797362.View ArticleGoogle Scholar
- Project B: Mankind and Animals in the Mid Zambezi Valley. 2001, Montpellier, France: Zimbabwe. CIRAD- EmvtGoogle Scholar
- Dunham KM: The diet of impala (Aepyceros melampus) in the Sengwa wildlife research area, Rhodesia. J Zool. 1980, 192: 41-57.View ArticleGoogle Scholar
- Fryxell JM: Food limitation and demography of a migratory antelope, the white-eared kob. Oecologia. 1987, 72: 83-91. 10.1007/BF00385049.View ArticleGoogle Scholar
- Owen-Smith N: Demography of a large herbivore, the greater kudu Tragelaphus strepsiceros in relation to rainfall. J Anim Ecol. 1990, 59: 893-913. 10.2307/5021.View ArticleGoogle Scholar
- Mduma SAR, Sinclair ARE, Hilborn R: Food regulates the Serengeti wildebeest: a 40-year record. J Anim Ecol. 1999, 68: 1101-1122. 10.1046/j.1365-2656.1999.00352.x.View ArticleGoogle Scholar
- Gaillard JM, Hebblewhite M, Loison A, Fuller M, Powell R, Basille M, Van Moorter B: Habitat–performance relationships: finding the right metric at a given spatial scale. Phil Trans R Soc B. 2010, 365: 2255-2265. 10.1098/rstb.2010.0085.PubMed CentralView ArticlePubMedGoogle Scholar
- Morrison TA, Bolger DT: Wet season range fidelity in a tropical migratory ungulate. J Anim Ecol. 2012, 81: 543-552. 10.1111/j.1365-2656.2011.01941.x.View ArticlePubMedGoogle Scholar
- Estes R: The Behaviour Guide to African Mammals. 1991, Berkeley: University of California PressGoogle Scholar
- Jarman MV:Impala Social Behavior. Territory, Hierarchy, Mating and the Use of Space. Advances in Ethology, Vol. 21. Edited by: Wickler W, Curio E. 1979, Berlin: Parey,Google Scholar
- Dasmann RF, Mossman AS: Population studies of impala in southern Rhodesia. J Mammal. 1962, 43: 375-395. 10.2307/1376947.View ArticleGoogle Scholar
- Jarman PJ: Seasonal distribution of large mammal populations in the unflooded Middle Zambezi Valley. J Appl Ecol. 1972, 9: 283-299. 10.2307/2402062.View ArticleGoogle Scholar
- Fryxell JM, Holt RD: Environmental change and the evolution of migration. Ecology. 2013, 94: 1274-1279. 10.1890/12-0668.1.View ArticlePubMedGoogle Scholar
- Switzer PV: Past reproductive success affects future habitat selection. Behav Ecol Sociobiol. 1997, 40: 307-312. 10.1007/s002650050346.View ArticleGoogle Scholar
- Fryxell J, Hazell M, Borger L, Dalziel BD, Haydon DT, Morales JM, McIntosh T, Rosatte RC: Multiple movement modes by large herbivores at multiple spatiotemporal scales. Proc Natl Acad Sci USA. 2008, 105: 19114-19119. 10.1073/pnas.0801737105.PubMed CentralView ArticlePubMedGoogle Scholar
- Bolger DT, Newmark WD, Morrison TA, Doak DF: The need for integrative approaches to understand and conserve migratory ungulates. Ecol Lett. 2008, 11: 63-77.PubMedGoogle Scholar
- Fairall N: Production parameters of the impala, Aepyceros melampus. S Afr J Anim Sci. 1983, 13: 176-179.Google Scholar