Skip Navigation


Behavioral Ecology Advance Access originally published online on May 2, 2007
Behavioral Ecology 2007 18(4):635-644; doi:10.1093/beheco/arm028
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Lay Summary
Right arrow All Versions of this Article:
18/4/635    most recent
arm028v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (5)
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Ryan, S.
Right arrow Articles by Getz, W.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Ryan, S.
Right arrow Articles by Getz, W.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© The Author 2007. Published by Oxford University Press on behalf of the International Society for Behavioral Ecology. All rights reserved. For permissions, please e-mail: journals.permissions@oxfordjournals.org

Ecological cues, gestation length, and birth timing in African buffalo (Syncerus caffer)

SJ Ryana,b, CU Knechtelc and WM Getza,d

a Department of Environmental Science, Policy and Management, University of California, Berkeley, CA 94720, USA b Museum of Vertebrate Zoology, 3101 VLSB, University of California, Berkeley, CA 94720, USA c Centre for Wildlife Management, University of Pretoria, Pretoria 0002, South Africa d Mammal Research Institute, Department of Zoology and Entomology, University of Pretoria, Pretoria 0002, South Africa

Address correspondence to S.J. Ryan. E-mail: sjryan{at}stanford.edu.

Received 29 August 2006; revised 16 February 2007; accepted 25 February 2007.


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 Methods
 RESULTS
 DISCUSSION
 REFERENCES
 
We examined annual variation in the timing of conception and parturition in the African buffalo (Syncerus caffer) and the synchrony of birth timing with resource cues, using 8 years of monthly birth, rainfall, and vegetation data, measured as Normalized Difference Vegetation Index (NDVI). Monthly births had the strongest significant correlations with NDVI and rainfall levels 12 and 13 months in the past, respectively. In addition, the synchrony of current year births corresponds most strongly to the synchrony of the previous year's NDVI distribution. Because the gestation period of buffalo has been estimated to be around 11 months, these findings suggest that improved protein levels, occurring approximately a month after the first green flush of the wet season, are either a trigger for conception or conception has evolved to be synchronous with correlated environmental cues that ensure females enter a period of peak body condition around the time of conception and/or parturition. With a gestation period of approximately 340 days, parturition occurs to take advantage of the period when forage has its highest protein content. A comparative analysis of gestation periods within the subfamily Bovinae indicates that African buffalo have a protracted gestation for their body size, which we suggest is an adaptation to their seasonal environment. We also found that interannual variation in the birth distribution suggests a degree of plasticity in the date of conception, and variation in the number of calves born each year suggest further synchrony at a timescale longer than a single year.

Key words: birth synchrony, breeding ecology, NDVI, phenology.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 Methods
 RESULTS
 DISCUSSION
 REFERENCES
 
Multiple hypotheses have been offered as explanations for birth synchrony in ungulates: 1) predator swamping through dilution or group protection of young (Estes 1976Go; Sekulic 1978Go; Estes RD and Estes RK 1979Go; Keech et al. 2000Go; Testa 2002Go); 2) optimization of social development as a cohort within a group (Brown 1985Go; Pfeifer 1985Go); or 3) timing births to coincide with periods of optimal environmental conditions (Rutberg 1987Go). These hypotheses are not necessarily mutually exclusive and may act in concert. Here we focus on the third hypothesis, which points to the role of resource availability as driving the timing of parturition in African buffalo (Syncerus caffer). Ungulates may time parturition to coincide with the availability of high-quality vegetation to support the energetic costs of lactation (Gaillard, Delorme, and Jullien 1993Go), to maximize the growth rate of neonates to assist predator avoidance (Clutton-Brock et al. 1982Go), and to replenish maternal reserves for future reproduction while conditions are optimal (Gaillard, Delorme, Jullien, and Tatin 1993Go; Keech et al. 2000Go). However, we hypothesize that buffalo are also conceiving while conditions are optimal, and due to a protracted gestation period of slightly more than 11 months, additionally benefit from good conditions at parturition.

African buffalo exhibit similar seasonal social ecology to many other ungulates, in which males form bachelor groups in the nonbreeding season to forage and improve their condition and then compete for mates in the breeding herd during the breeding season (Clutton-Brock et al. 1979Go, 1982Go; Halley and Mari 2004Go). The social ecology of buffalo seasonal breeding aggregation depends largely on resource availability: a threshold appears to exist at which the foraging cost for males of remaining in the group outweighs the reproductive opportunity (Turner 2003Go; Turner et al. 2005Go). In arid and semiarid systems, it is suggested that herbivore aggregation is also limited by a patchy spatial availability of resources (Fryxell 1995Go) whose quality varies with very local rainfall. This seasonal and spatial restriction by resources may act to limit male reproductive opportunity within the breeding herd due to reduced feeding availability in a mixed herd and competing with other males for mating opportunities. In addition, buffalo cows must attain sufficient body condition to conceive and contend with the competing potential energetic costs of gestation, lactation, or tending of the prior year's calf (Prins 1996Go). African buffalo are highly gregarious and appear to have a social system within the breeding herd in which females and their offspring move together within a looser female hierarchy (Mloszewski 1983Go; Prins 1996Go), using a set of vocalizations and postures to maintain herd cohesion and order (Mloszewski 1983Go; Estes 1991Go). Calves have been reported to suckle for up to 18 months in some study areas (Prins 1996Go).

In accord with the principles of optimal foraging theory, to maximize reproductive output, an organism should maximize food quality: for herbivores, this means maximizing protein intake (Kennish 1997Go). For grazers in a semiarid ecosystem, protein is a limiting nutrient (Bell 1971Go; Owen-Smith 1982Go; Owen-Smith and Novellie 1982Go). Although described as a "supreme bulk grazer," and therefore able to consume a high quantity of matter (Owen-Smith and Cumming 1993Go), buffalo still lose condition as grass biomass increases due to a decrease in the protein content of their diet (Fryxell 1995Go). Prins (1996)Go found that African buffalo in Manyara National Park, Tanzania, conceive at peak cow condition and suggested that the long gestation of buffalo additionally means that calving coincides with good nutritional condition. However, a direct connection between maximum protein availability and conception has not previously been made for buffalo.

African buffalo have birthing seasons reported to coincide with the wet season in multiple sites (Grimsdell 1973Go; Sinclair 1977Go; Mloszewski 1983Go; Prins 1996Go; Turner 2003Go) and a gestation period with a mean of around 340 days (Bertschinger 1996Go). In a study of buffalo diet in Kruger National Park, it was shown that remotely sensed vegetation greenness (Normalized Difference Vegetation Index, NDVI) correlated highly with the proceeding month's level of fecal nitrogen in buffalo (Ryan 2006Go), a widely used proxy for dietary protein intake in both domestic and wild ungulates (Holechek et al. 1982Go; Leslie and Starkey 1985Go; Wofford et al. 1985Go; Bartiaux-Thill and Oger 1986Go; Hobbs 1987Go; Leite and Stuth 1990Go, 1992, 1995Go; Nunez-hernandez et al. 1992Go; Lyons et al. 1995Go; Wehausen 1995Go; Grant et al. 1996Go; Kucera 1997Go). NDVI provides us with an important measure of habitat quality in this system because the patchiness of rainfall and the difference between sustained low rainfall and periodic or episodic events are hard to distinguish with precipitation data (Scanlon et al. 2005Go). As reviewed by Pettorelli et al. (2005)Go, NDVI is a well-established proxy for vegetation productivity and phenology. The body condition of adult female buffalo, both with and without a calf at foot, coincided with peak monthly remotely sensed NDVI values in Kruger National Park (Ryan 2006Go).

Therefore, we hypothesize that conception occurs around peak protein content in the vegetation, 1 month after the peak NDVI, when buffalo cows are in peak condition, whereas parturition occurs within a month of the peak protein content 1 year later. This hypothesis is supported by our data; however, insufficient information exist to fully distinguish the individual roles of rainfall and protein content in forage.

To further explore the observed phenomenon that African buffalo conceive and give birth in optimal conditions, we compare their gestation period with other members of the subfamily Bovinae to assess whether they demonstrate a protracted gestation that allows for a synchronized annual cycle of conception and parturition. This is consonant with a hypothesis that although this is a "relaxed" (Rutberg 1987Go) seasonal environment as compared with Northern climes, there has been sufficient selective pressure for African buffalo to evolve extended gestation beyond a physiological expectation for their bodyweight.

The occurrence of droughts in this system can alter the timing of wet season onset, which can result in lowered vegetation production and palatability of standing grass. The plasticity of the breeding season in mammals subject to a highly variable climate may be greater than those that have evolved with reliable resource cues. In a study by Sinclair et al. (2000)Go and studies by Rutberg (1984)Go, 1987Go), a data set was summarized from a period in 1977 (Sinclair 1977Go), in which 49 birth events were recorded. With our data that includes 786 births during a period spanning 1993–2001, we are able to evaluate the potential for a response to annual variations in ecological resource cues.


    Methods
 TOP
 ABSTRACT
 INTRODUCTION
 Methods
 RESULTS
 DISCUSSION
 REFERENCES
 
Study site
The Klaserie Private Nature Reserve (KPNR) is located in the Limpopo Province of South Africa, bordering Kruger National Park on its western boundary (24°3–22'S, 31°2–19'E; 303–535 m above sea level; 57 800 ha) (Figure 1). The main geological substrate is granitic gneiss, part of the extensive granalitic system underlying most of the country. The climate follows a subtropical savanna pattern: both annual temperature and rainfall follow a unimodal distribution (Figure 2). Mean annual rainfall for this study period is 574 mm. We define 2 main seasons based on rainfall and temperature: a hot wet "summer" season from October to March and a cool dry "winter" season from April to September. In this study, our years are from October to October, with the year (19XX) assigned as the calendar year for January. KPNR contains the range of 3 main herds of buffalo. The range and habitat selection of these buffalo are summarized in a prior study (Ryan et al. 2006Go).


Figure 1
View larger version (26K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 1 The study site: Klaserie Private Nature Reserve, South Africa.

 

Figure 2
View larger version (28K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 2 Rainfall and temperature (mean ± SE) of Klaserie Private Nature Reserve, 1991–2000 (adapted from Ryan et al. 2006Go).

 
Births
A data set containing 786 calving events by month for one herd were recorded by one of us (C.U.K.) between 1993 and 2001.

Precipitation
Mean monthly rainfall was taken from a maximum of 7 stations throughout KPNR. These data were transcribed from logbooks into Excel and monthly averages calculated to create a monthly distribution of rainfall through the study period (Figure 3).


Figure 3
View larger version (33K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 3 Monthly births (foreground, white), rainfall (millimeter) (background, black) (a), and NDVI (background, gray) (b) for the study period.

 
Remotely Sensed Vegetation Data
NDVI data from National Oceanographic and Atmospheric Association's (NOAA's) Advanced Very High Resolution Radiometer (AVHRR) sensor aboard weather satellites processed by the Global Inventory Monitoring and Modelling Studies group (GIMMS) at the National Aeronautical and Space Administration (NASA) at a scale of 8 km2 were downloaded from the Africa Data Dissemination Service (ADDS) (Pinzon 2002Go; Pinzon et al. 2004Go; Tucker et al. forthcoming). We use these data as a measure of the change in vegetation quality in the reserve. These data were imported into Agricultural Research Centre (ARC)/INFO 8.0 (Environmental Systems Research Institute 2000aGo) and ArcGIS 8.3 (Environmental Systems Research Institute 2000bGo), using ArcMap for visual assessment of available images. For viable images (no cloud cover over KPNR, no image corruption), 0–3 images per month were used to create monthly average grids of data clipped to the KPNR, using an Arc Macro Language (AML) script. We converted the raw 8-bit data (256) to the conventional scale (–1 to 1), where vegetated land usually takes on values in the range 0.1–0.7 (Pinzon 2002Go; Pinzon et al. 2004Go; Tucker et al. forthcoming). These monthly average values were averaged across the reserve to derive a monthly mean value in a vegetation quality distribution (Figure 3).

Birth timing
To find the peak of the birth distribution, monthly average proportions of births in the years 1993–2001 were compared using an ANOVA with a post hoc Tukey–Kramer's HSD test for means. A generalized least-squares model was also fit to the births versus month and year in SPLUS for Windows to compare annual and intra-annual variance in births. To examine breeding at a finer scale, we took a set of 28 individually identified females that were followed through 7 consecutive years and calculated both an individual fecundity and a group proportion of birthing within a year. We did not have information on calf survival for these individuals, as would be necessary to formally test the effects of lactation or tending on conception, but we informally describe calving sequence for these individuals.

Environmental cue timing
We cross-correlated the resource cue (precipitation, vegetation) against the births over a 15-month window. We used a conservative backward stepwise regression (P > 0.05 to exit) to test the distribution of births against each of the lagged resource cues (0–15 month lagged) to identify the significant lags. We chose a backward stepwise regression as appropriate because we had prior assumptions from the cross-correlation about which time lags we were modeling, and a backward stepwise regression is initiated with all parameters (in our case, time-lagged cues) in a full model. In a mixed stepwise regression (0.05 to enter and 0.05 to exit), we achieved the same results, suggesting that we did not lose any significant parameter contributions in the selection process. In addition, we tested the parameter "year" in all models, with no significant results, so we do not report it in the model fit for this parameter.

Evenness of births and environmental cues
To quantify the evenness of births and environmental cues each year, we applied the method of Bronikowski and Webb (1996)Go, which has been used in similar studies comparing ungulate behavioral responses to rainfall variability (Sinclair et al. 2000Go; Brashares and Arcese 2002Go). The evenness index Jq is given as

Formula
where t = 1994, ..., 2001 and i = 1, ..., 12 denotes the month of the year (respectively, from October of year t – 1 to September of year t); q = b, r, and g, respectively, represents births, rain, and NDVI (greenness); and pFormula(t) is the proportion of births or environmental cue in the ith month of the tth year. Evenness ranges from 0 (e.g., all births or rain occur in one month) to 1 (same levels each month).

For purposes of comparison, we followed the methods of Sinclair et al. (2000)Go and compared the evenness measure of births with that of each resource, both within the same year and the preceding year, using a ratio of Jbirths/Jresource. The ratio that most closely approximates 1 is the indication of similarity of birth and resource similarity (Table 1).


View this table:
[in this window]
[in a new window]

 
Table 1 Synchrony of births and resource cues

 
Evidence of protracted gestation
We collated data for gestation period in days (d), and bodyweight (w) in kilograms for 24 of 25 members of the subfamily Bovinae (Myers et al. 2006Go), as taken from the Cetartiodactyla supertree of Price et al. (2005)Go, using the software TREEVIEW and TREEedit (Page 1996Go). The phylogeny of the Bovinae is not well resolved, in particular, a reconstruction of the tribe bovini containing the African buffalo has not been resolved (Buntjer et al. 2002Go), although work is underway using mitochondrial DNA and Y-chromosomal sequences (Lenstra JA, personal communication), so a fully corrected analysis is not currently feasible. The genus Syncerus is monophyletic, and as this is an analysis at the subfamily level, we were comfortable, using a least-squares fit of log(d) versus log(w) at the genus level, to derive an allometric scaling equation (power law), d = bwa, where a is the regression coefficient or slope and b is the intercept (Martin et al. 2005Go). To predict the expected gestation period for buffalo, we used the other 23 members of the Bovinae for which we have data to create the power law model, using generic means, and then interpolated to examine where buffalo occur on the prediction interval. Because we had few data points in the regression (n = 9 genera), the prediction interval was quite wide, though we still retained the 2-tailed 95% prediction interval (but see [Garland and Ives 2000Go]) (Figure 4).


Figure 4
View larger version (10K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 4 Average monthly proportion of births (mean ± SE) from 1993 to 2001.

 
All statistics were run in SAS JMP 5.1 (SAS 1989–2003) or SPLUS 6.0 for Windows.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 Methods
 RESULTS
 DISCUSSION
 REFERENCES
 
Birth timing
We found that across this study period, January is the peak month of birthing (Figure 5). An ANOVA with post hoc means comparisons showed that there is a 3-month peak from January to March wherein the mean proportion of births occurring in January is significantly higher than all other months, and February and March are significantly different from all other months but not from each other (Fratio = 23.96, degrees of freedom [df] = 11, 95, P < 0.0001; Tukey–Kramer HSD, {alpha} = 0.05). The generalized least-squares model showed that the intra-annual variability of births was significant, whereas interannual variance was not (ANOVA: model df = 83, intercept df = 1, Fratio = 160.93, P < 0.0001; Fratio = 15.17, P < 0.0001; month: df= 11, Fratio = 23.97, P < 0.0001; year: df = 1, Fratio = 1.06, P = 0.31).


Figure 5
View larger version (12K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 5 The relationship of gestation period (days) to female bodyweight (kilogram) for the subfamily Bovinae excluding African buffalo. The power law of log gestation (days) against log bodyweight (kilogram) at the level of genus is shown (solid line) with 95% PIs (dashed lines), with raw data for the 9 Bovinae genera, depicting the African buffalo (Syncerus caffer) (pointer).

 
The group of 28 individually identified breeding females had an average individual fecundity (mean ± standard error [SE]) of 0.43 ± 0.02 calves/year. We found that the average proportion of females calving in a year is [mean (95% confidence interval (CI))] 0.43 (0.20–0.65). There was considerable interannual variability in the calving proportion and, in particular, 1998 appeared to be a particularly good calving year with 25 (89%) of the females calving (Figure 6). We bootstrapped the data from 28 females 10 000 times by rearranging each of the rows at random (i.e., each female's reproductive output for the 7-year period was held constant, but the years in which they had their calves was randomized). In none of these runs did 23 or more females reproduce at the same time in any of the 7 years, and in only one did 22 females reproduce at the same time (the mode of this distribution is 12). This implies that the probability of observing 25 or more by chance of the 28 individuals reproducing, as was the case in 1998, is P < 0.0001.


Figure 6
View larger version (13K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Figure 6 The proportion of 28 individually identified females calving in a given year (gray columns), with mean (solid line) and 95% CI (dashed lines)—note the high proportion of 89% in 1998.

 
The odds of having a calf in 2 consecutive years was 0.5 (prior year) and 0.58 (proceeding year), and of the females calving in 1998, 24% had calved the year previously and 20% calved in the following year. Seven individuals were observed to have calves in 3 consecutive years within this data set.

Environmental cue timing
Cross-correlations of the birth distribution from 1996 to 2001 against monthly rainfall and greenness (NDVI) lost significance (and meaning) at a 14-month lag. The highest correlation (Pearson's r = 0.69, P < 0.0001) is a 12-month lag in NDVI (1 year prior). For monthly rainfall, the highest correlation occurs at a 13-month lag (Pearson's r = 0.57, P < 0.0001). Positive correlations occur for both resource cues at a 0- to 3-month lag and between 10 and 14 for NDVI and 11 and 15 for rain, corresponding to the current and prior wet season months. These results support the idea that conception is tied to resources but that the extended gestation implies that the same high-quality conditions are realized at parturition a mean of 340 days later.

A series of backward stepwise linear regressions (SLR) yielded the following models for each resource. For monthly rainfall (r), the SLR selected lag 2, lag 11, and lag 13 as best explanatory variables, with most support for lag 13 (Ftotal = 23.43, P = 0.0001; rlag 2: F = 13.66, P = 0.0134, rlag 11: F = 4.71, P = 0.03, rlag 13: F = 31.20, P < 0.0001); y = 0.06(rlag 2) + 0.03(rlag 11) + 0.09(rlag 13) + 0.32.

For NDVI (v) the SLR selected lag 1, lag 3, and lag 12 as most explanatory variables (Ftotal = 51.81, P = 0.0001; vlag 1: F = 31.00, P < 0.0001; vlag 3: F = 24.36, P < 0.0001; vlag 12: F = 12.28, P = 0.0007); y = 53.41(vlag 1) – 35.86 (vlag 3) + 30.96 (vlag 12) – 12.27.

Evenness of births and environmental cues
The distributions of births, rainfall, and NDVI across the study period are shown in Figure 3. The average evenness, Jb, for births over the study period was 0.55 ± 0.05 (mean ± SE), for rainfall it was 0.71 ± 0.02, and for NDVI it was 0.72 ± 0.5. The mean ratios of Jb/Jq suggest that the birth evenness is most similar to the prior year's vegetation (NDVI) measure of evenness (0.79 ± 0.10) (Table 1).

Evidence of protracted gestation
The power law equation derived from the least-squares fit of log transformed gestation period (d) against female bodyweight (w) at the genus level is d = 5.25w0.057; (a = 0.057 ± 0.02 [mean ± SE], t = 2.94, P = 0.025; b = 5.25 ± 0.097 [mean ± SE], t = 54.06, P < 0.0001]. The estimate for African buffalo gestation with the individual estimate 95% prediction interval (PI) is 272 (95% PI: 237–310) days. Thus, the reported buffalo gestation period with a mean of 340–343 days (Sinclair 1977Go; Bertschinger 1996Go) is outside the 95% CI of its predicted range (Figure 5).


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 Methods
 RESULTS
 DISCUSSION
 REFERENCES
 
The proposed mechanisms for birthing synchrony in ungulates (i.e., timing to good conditions for lactation, replenishing maternal reserves, and providing high-quality food for offspring) may all contribute to the synchrony and timing of births in this system.

In savanna ecosystems, the seasonality or variation of resource availability is largely predicated on rainfall. A cursory examination of Figure 3 shows the close relationship between rainfall and vegetation greenness in this system. However, it is clear that vegetation productivity and quality are subject to more than just precipitation (Scanlon et al. 2005Go). Vegetation greenness, measured as NDVI, lags the rainfall by 1 month, and crude dietary protein lags NDVI by 1 month. It is therefore unsurprising that the synchrony of births is similar to that of resources. We found moderate values of birth synchrony for this period (Jb = 0.55 ± 0.05, mean ± SE), slightly more synchronous than the value found by Sinclair et al. (2000)Go of 0.695. In comparison, the rainfall in our system was more uneven, Jr = 0.71 ± 0.02 (mean ± SE), versus 0.935 in the western Serengeti (Sinclair et al. 2000Go). Although a formal evenness analysis was not conducted, a comparison between Manyara and the Serengeti by Prins (1996)Go pointed out that at Manyara, 50% of the births occur over 17 weeks, whereas in the Serengeti, 50% occur over 7 weeks. Manyara has 2 consecutive wet seasons (short rains and long rains), and the standing groundwater is higher than in the Serengeti in the dry season, but the birth season still peaked in the long rains. Even in captive buffalo, birthing appeared to have a peak in the austral midsummer (Skinner et al. 2006Go). The seasonality of buffalo birth has been noted in populations of Cape buffalo, but data for the subspecies of forest buffalo, Syncerus caffer nanus, are very scant. It is assumed that as the same species, that reproduction is similar, and new data suggests that gestation is 11 months also, with potential seasonality in birthing (Korte L, personal communication).

Following the method of Sinclair et al. (2000)Go, we found that the synchrony of births most closely approximated to that of the prior year's vegetation measure distribution, suggesting that there is a strong link between the pattern of resource distributions in 1 year and the proceeding birth distribution. In a similar finding, Visscher et al. (2004)Go suggested that because buffalo conceive into the late wet season in Kruger National Park, they may not be responding to prior rainfall but anticipating the green flush of grass in 11 months. But how could this timing occur? The idea that buffalo can use cues in the environment is not novel, whereas they are bulk grazers and they do exhibit selective grazing, preferring certain grasses over others (Macandza 2004; Turner et al. 2005Go). As a herd, buffalo may be capable of assessing the quality of grazing patches. In the Serengeti, crude protein was shown to decline with the maturation of grass, and Fryxell (1995)Go cited experiments wherein herbivores selected patches of intermediate biomass, suggesting that they are behaviorally capable of responding to the trade-off between abundance and quality. This supports the idea that the new, green leaf grass is the most nutritious and likely to be selected preferentially. However, it will still take time after the onset of the wet season for there to be sufficient biomass for buffalo cows to attain sufficient condition to conceive. Thus there will be a lag between rainfall and sufficient biomass and protein in grass to support conception.

In a comparison of a suite of African ungulate species, Sinclair et al. (2000)Go posited that birthing synchrony in buffalo evolved as predator satiation. This echoed Rutberg (1984Go, 1987)Go who asserted that buffalo calves are precocious "followers" rather than "hiders": part of large groups in open habitat, promoting defence against predators. However, the observed length of the buffalo birthing season, both in this study and others (Sinclair 1977Go; Prins 1996Go; Sinclair et al. 2000Go), is far longer than that of other follower species, corresponding to Rutberg's (1987)Go alternate hypothesis that the climate enables a longer birthing season, with peaks occurring at times of good nutritional condition. In addition, Sinclair (1977)Go noted that buffalo calves tend to hide with their mothers in the bush, away from the herd, after birth, refuting this argument, although this is unusual for our system. Thus, it is more likely that the synchrony of buffalo parturition is driven by resource availability, which confers the additional advantage of predator swamping. We suggest that buffalo have individual behavioral adaptations to predation (hiding young calves, keeping them in the center of the herd, or providing them the protection of the whole herd) for protection of neonates, which may simply be a function of habitat. The open savanna would necessitate "following," whereas bushveld allows for "hiding."

We found that the availability and seasonality of resources around conception correlate with the birth distribution for this population. Due to an uncharacteristically long gestation period for the Bovinae of over 11 months, the condition of the mother at both conception and during lactation will be maximized at vegetation quality peaks. Our examination of gestation times within the subfamily Bovinae suggests that African buffalo have a protracted gestation period to take advantage of these resource peaks (see Figure 4). An adult female buffalo carries a calf for several months during the dry season in this system, followed by a period of lactation, in which her maternal resources will be depleted and her protein intake requirement may double (based on cattle estimates [Prins 1996Go]). Depending at what point in the gestation the female enters the dry season, the following sequence of life-history events could be affected. This points to a much larger and in-depth suite of questions regarding the sequence of energetics for which we do not currently have the data. An explanation for the longer gestation in African buffalo is that this is part of a larger suite of life-history traits, characterizing "K-strategists"—those species which have a late age of sexual maturity, a long life span, and large brain size (Lee 1997Go). This would emphasize the role of weaning and lactation time for buffalo, which can be quite long (Bertschinger 1996Go). However, the high rate of juvenile mortality in buffalo appears in contrast to this idea.

Over the 8 years of data presented here, births occurred in all months except October. The NDVI distribution is relatively even and does not change significantly for 3 months around its peak during the year, so we were cautious in making conclusive statements as to its role in cueing conception by specific month. However, we find that in our correlations with all resource cues, both in current and previous years, the strongest correlation corresponds to a 1-year lag with NDVI (an 11-month lag with protein). This represents the time at which buffalo breed. Thus, we infer that they come into breeding condition when the forage condition is optimal. Peak rainfall lags the greenest month of the year by a month, meaning that we see a birth peak month followed by a rainfall peak month. This was also found in buffalo in Hluhluwe-Umfolozi Park, to the south of Klaserie (Turner 2003Go), and is in fact the pattern seen in the study of Sinclair et al. (2000)Go in the Serengeti. In their study, "green biomass" was predicted by rainfall, and births lagged peak biomass by 1 month. We suggest that rather than parturition nearly corresponding to biomass, the Serengeti buffalo conceived at the Serengeti green biomass peak, and parturition was observed 11 months later. With a longer period over which to examine this relationship, it is possible to distinguish between the prior year and conditions at parturition. However, Sinclair et al. (2000)Go also found that the birthing peak for buffalo was around 5 months out of phase with the peak protein content of their forage. The rainfall season of Eastern Africa is different from that of Southern Africa leading to a different pattern of protein availability than we find in our system, and we therefore do not find the same corresponding lag. We posit that buffalo in our more subtropical system are reaching conception peak, whereas the protein content of available grasses is highest, which coincides with the peak growing season.

In years in which rainfall is low or has a very late onset, grass production and availability will also be low and/or late. As buffalo have evolved such that breeding season occurs during the wet season, years in which conception cues are not strong are likely to lead to less synchronized or disrupted birth patterns in the following year. The year 2000 had the highest rainfall of the entire study period. The birth distribution in 2001 had an early onset, with births occurring as early as November (Figure 3), which is unique in the data set but anecdotal given that this could have happened by chance. The population will likely reflect some abiotic effects that affect demography, even at this timescale. The area experienced a severe drought in 1992, 2 years prior to this study, in which there was a 60% decline in the buffalo population in the neighboring Kruger National Park (Mills et al. 1995Go), and a 53% decline in the buffalo counts for KPNR (Ryan et al. 2006Go). The lack of resources due to low rainfall years may drive prime females either into nonreproductive condition or cause them to leave the population, perhaps migrating to better habitat. This would likely result in fewer offspring in following years in this data set, regardless of phenological resource pattern. Although the data are sparse, 1995 and 1996 fall below the estimated lower 95% of births per year in this data set (95% CI: 79–117) at 56 and 77 births, respectively, which we can only suggest might be indicative of a lagged effect of the 1992 drought leading to low-quality breeding-age females (4+) and fewer births.

To attempt to account for the numerical fluidity of this population, we instead looked at a subset of data. An examination of the breeding success of 28 individually identified females with 7 years of complete records (1994–2001) reveals that, in 1998, 89% reproduced, compared with a mean proportion of 43% for the period (Figure 6). We found that within the data set, calving did not appear to be constrained by a birth in the previous year, with females equally or slightly more likely to give birth in consecutive years; a quarter of the females calved in 3 consecutive years. The lack of information on individual calf survival for these females does not allow us to fully examine whether these females were essentially released from the constraints of lactating by calf loss and thus able to regain condition prior to the next breeding season. Grimsdell (1973)Go noted that the duration of the calving interval in his study in Western Uganda depended on nutrition and that cows in better condition tended to produce more offspring. Calves can be weaned as early as 4–5 months of age (Bertschinger 1996Go); however, in this system, it has been observed that calves suckle up to at least the seventh month of the following pregnancy. This is supported by evidence given in Estes (1999)Go that females stop lactating in the seventh month of pregnancy despite observations of suckling continuing up to 15–18 months or until the next calf is born (Prins 1996Go; Apps 2000Go). This suggests that if a female can maintain body condition to meet the demands of gestation and lactation, a yearly calf is possible. The high rates of pregnancy, from 47% to 83% in 5 different parks reported for adult female buffalo (Bertschinger 1996Go), combined with high calf mortality (30–70%, Bertschinger 1996Go) suggest that early postpartum condition of both the mother and, perhaps more importantly, the resources, will drive population growth. Within this study period, 1998 had the lowest rainfall on record, without apparent immediate or following year effect on either the onset of births or the number of births in the overall population; however, of the 28 individuals that calved in 1998, only 20% produced calves in the following year. This suggests that the following dry season may have been too harsh for many of the cows to complete gestation again. This indicates that although the seasonality of resources may influence the seasonality of birth events in the following year, the interaction between resource quality and fecundity is far more complex and may even operate at a timescale beyond the scope of this paper.

In addition, it is not just breeding females that rely on the resources and resource cues. In a study of buffalo in Hluhluwe-Umfolozi, Turner et al. (2005)Go suggest that males who join breeding herds face an energetic trade-off between obtaining sufficient food and engaging in reproductive activity. As the grass growing season tapers and the quality of food is reduced, it is likely that both sexes must face this trade-off. A habitat selection study (Ryan et al. 2006Go) showed that these buffalo range further in the dry season, most likely in search of forage that is still green. The combination of spending more time searching and more time ruminating would likely curtail reproductive activities after the wet season, although this may be a simplistic explanation.

Our findings of the close relationship between buffalo breeding ecology and seasonal availability of resources have further implications during periods of extreme conditions. A study of Alaskan moose with highly synchronized parturition suggested that they have adapted timing and synchrony of births to long-term patterns of climate and could not adjust to proximate resource changes (Bowyer et al. 1998Go). The authors concluded that this lack of proximate adjustment in the face of global warming could compromise the long-term viability of the population. It has been argued that greater birthing synchrony and a shorter breeding season will occur in temperate or subarctic climates (Rutberg 1987Go; Gaillard, Delorme, Jullien, and Tatin 1993Go; Gaillard et al. 1998Go; Rubenstein and Wikelski 2003Go), whereas tropical and subtropical habitats enable longer birthing periods (Rutberg 1987Go), except where predation provides a strong selective force (Estes 1976Go). Buffalo, with their extended birth season and a demonstrated capacity to adjust to the onset of seasonal forage availability, are less likely to suffer such dramatic impacts as their northern counterparts.

Despite the availability of more data than prior studies of buffalo birth seasonality, our study still spans a short period for this savanna ecosystem. The study period reflects the wetter second half of an interdecadal rain cycle (Ogutu and Owen-Smith 2003Go), yet has considerable annual variation in rainfall within the study period. This may create noise in the study system from which it is hard to untangle behavioral breeding responses. In addition, individual females do not breed every year, thereby confounding the analyses over a short period of time. In low rainfall years, cows may build condition on any and all available biomass and then conceive toward the end of the breeding season. This may in turn cause disruption in consecutive years due to females having dependent yearling calves during the next mating opportunity. Another possibility is that some cows may "skip" a bad year and be in high condition very early in the following mating season. Only with the availability of data sets, spanning 4 decades or more, can we hope to begin to capture the impacts of an 18-year interdecadal cycle on buffalo demography in the South African lowveld.


    ACKNOWLEDGEMENTS
 
We wish to thank Tiffany Kim for entering calving and reproduction data; thanks to M.I. Westphal, S.E. Reed, W.C. Turner, and A.L. Bidlack for edits and comments. Thanks also to Prof. J.S. Brashares and E.A. Lacey for reviewing drafts and making edits and comments. Thanks to Collin Rowles, Warden, and to the owners of KPNR for the opportunity to do this research. Thanks to Mike Peel at the ARC in Nelspruit, South Africa, for data and information regarding KPNR maps. S.J.R. research was supported by EPA-STAR fellowship FP-916382. NSF Grant DEB-0090323 to W.M.G. additionally supported S.J.R. research and partially supported C.U.K. and W.M.G. research. C.U.K. research was additionally supported by several institutions, primarily the German Academic Exchange Service (DAAD) and Daimler-Chrysler South Africa. All work conducted in the KPNR was of a noninvasive nature and in compliance with the laws of South Africa.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 Methods
 RESULTS
 DISCUSSION
 REFERENCES
 
Apps P. Smithers' mammals of Southern Africa, a field guide (2000) Cape Town (South Africa): Struik Publishers. 364.

Bartiaux-Thill N, Oger R. The indirect estimation of the digestibility in cattle of herbage from Belgian permanent pasture. Grass Forage Sci (1986) 41:269–272.[CrossRef]

Bell HV. A grazing ecosystem in the Serengeti. Sci Am (1971) 225:86–93.

Bertschinger HJ. Reproduction in the African buffalo: a review. In: Symposium on the African buffalo as a game ranch animal—Penzhorn BL, ed. (1996) Ondestepoort (South Africa): South African Veterinary Association. 62–74.

Bowyer RT, Van Ballenberghe V, Kie JG. Timing and synchrony of parturition in Alaskan moose: long-term versus proximal effects of climate. J Mammal (1998) 79:1332–1344.[CrossRef][Web of Science]

Brashares JS, Arcese P. Role of forage, habitat and predation in the behavioural plasticity of a small African antelope. J Anim Ecol (2002) 71:626–638.[CrossRef]

Bronikowski A, Webb C. Appendix: a critical examination of rainfall variability measures used in behavioral ecology studies. Behav Ecol Sociobiol (1996) 39:27–30.[CrossRef][Web of Science]

Brown RE. The primitive ungulates: orders Tubulidentata, Proboscidea and Hyracoidea. In: Social odours in mammals—Brown RE, McDonald DW, eds. (1985) Oxford: Clarendon Press. 235–244.

Buntjer JB, Otsen M, Nijman IJ, Kuiper MTR, Lenstra JA. Phylogeny of bovine species based on AFLP fingerprinting. Heredity (2002) 88:46–51.[CrossRef][Web of Science][Medline]

Clutton-Brock TH, Albon SD, Gibson RM, Guinness FE. Logical stag—adaptive aspects of fighting in red deer (Cervus-Elaphus L). Anim Behav (1979) 27:211–225.[CrossRef][Web of Science]

Clutton-Brock TH, Guinness FE, Albon SD. Red deer: behavior and ecology of two sexes (1982) Edinburgh: Edinburgh University Press.

Environmental Systems Research Institute Inc. ARC/INFO 8 (2000a) 8.02 ed. Redlands (CA): ESRI.

Environmental Systems Research Institute Inc. ArcGIS 8.3 (2000b) Redlands (CA): ESRI.

Estes RD. The significance of breeding synchrony in the wildebeest. East Afr Wildl J (1976) 14:135–152.

Estes RD. The behavior guide to Africa mammals: including hoofed mammals, carnivores, primates (1991) London: University of California Press.

Estes RD. The safari companion, a guide to watching African mammals (1999) White River Junction (VT): Chelsea Green Publishing Company.

Estes RD, Estes RK. The birth and survival of wildebeest calves. Z Tierpsychol (1979) 50:45–95.[Web of Science][Medline]

Fryxell JM. Aggregation and mitigation by grazing ungulates in relation to resources and predators. In: Serengeti II: dynamics, management and conservation of an ecosystem—Sinclair ARE, Arcese P, eds. (1995) Chicago (IL): Chicago University Press. 257–273.

Gaillard JM, Delorme D, Jullien JM. Effects of cohort, sex, and birth date on body development of roe deer (Capreolus-Capreolus) Fawns. Oecologia (1993) 94:57–61.[CrossRef][Web of Science]

Gaillard JM, Delorme D, Jullien JM, Tatin D. Timing and synchrony of births in roe deer. J Mammal (1993) 74:738–744.[CrossRef][Web of Science]

Gaillard JM, Festa-Bianchet M, Yoccoz NG. Population dynamics of large herbivores: variable recruitment with constant adult survival. Trends Ecol Evol (1998) 13:58–63.[CrossRef]

Garland T, Ives AR. Using the past to predict the present: confidence intervals for regression equations in phylogenetic comparative methods. Am Nat (2000) 155:346–364.[CrossRef][Web of Science][Medline]

Grant CC, Biggs HC, Meissner HH, Basson PA. The usefulness of faecal phosphorus and nitrogen in interpreting differences in live-mass gain and the response to P supplementation ingrazing cattle in arid regions. Onderstepoort J Vet Res (1996) 63:121–126.[Web of Science][Medline]

Grimsdell JJR. Reproduction in the African buffalo, Syncerus caffer, in Western Uganda. J Reprod Fertil Suppl (1973) 19:303–318.[Medline]

Halley DJ, Mari M. Dry season affiliation of African buffalo bulls at the Chobe riverfront, Botswana. S Afr J Wildl Res (2004) 34:105–111.

Hobbs NT. Fecal indices to dietary quality: a critique. J Wildl Manag (1987) 51:317–320.[CrossRef]

Holechek JL, Vavra M, Arthur D. Relationships between performance, intake, diet nutritive quality and fecal nutritive quality of cattle on mountain range. J Range Manag (1982) 35:741–744.

Keech MA, Bowyer RT, Ver Hoef JM, Boertje RD, Dale BW, Stephenson TR. Life-history consequences of maternal condition in Alaskan moose. J Wildl Manag (2000) 64:450–462.[CrossRef]

Kennish R. Seasonal patterns of food availability: influences on the reproductive output and body condition of the herbivorous crab Grapsus abolineatus. Oecologia (1997) 109:209–218.[CrossRef][Web of Science]

Kucera TE. Fecal indicators, diet, and population parameters in mule deer. J Wildl Manag (1997) 61:550–560.[CrossRef]

Lee PC. The meanings of weaning: growth, lactation and life history. Evol Anthropol (1997) 5:87–96.[CrossRef]

Leite ER, Stuth JW. Value of multiple fecal indices for predicting diet quality and intake of steers. J Range Manag (1990) 43:139–143.

Leite ER, Stuth JW. Fecal NIRS equations to assess diet quality of free-ranging goats. Small Ruminant Res (1995) 15:223–230.

Leslie DMJ, Starkey EE. Fecal indices to dietary quality of cervids in old-growth forests. J Wildl Manag (1985) 49:142–146.[CrossRef]

Lyons RK, Stuth JW. Fecal NIRS equations for predicting diet quality of free-ranging cattle. J Range Manag (1992) 45:238–244.

Lyons RK, Stuth JW, Angerer JP. Technical note—fecal Nirs equation field validation. J Range Manag (1995) 48:380–382.

Macandza VA, Owen-Smith N, Cross PC. Forage selection by African buffalo (Syncerus caffer) through the dry season in two landscapes of the Kruger National Park. S Afr J Wildl Res (2004) 34:113–121.

Martin RD, Genoud M, Hemelrijk CK. Problems of allometric scaling analysis: examples from mammalian reproductive biology. J Exp Biol (2005) 208:1731–1747.[Abstract/Free Full Text]

Mills MG, Biggs HC, Whyte IJ. The relationship between rainfall, lion predation and population trends in African herbivores. Wildl Res (1995) 22:75–88.[CrossRef]

Mloszewski MJ. The behavior and ecology of the African buffalo (1983) Cambridge: Cambridge University Press.

Myers P, Espinosa R, Parr CS, Jones T, Hammond GS, Dewey TA. The animal diversity [Internet]. (2006) Available at: http://animaldiversity.org. Accessed 19 Apr 2007.

Nunez-hernandez G, Holechek JL, Arthun D, Tembo A, Wallace JD, Galyean ML, Cardenas M, Valdez R. Evaluation of fecal indicators for assessing energy and nitrogen status of cattle and goats. J Range Manag (1992) 45:143–147.

Ogutu JO, Owen-Smith N. ENSO, rainfall and temperature influences on extreme population declines among African savanna ungulates. Ecol Lett (2003) 6:412–419.[CrossRef][Web of Science]

Owen-Smith N. Factors influencing the consumption of plant products by large herbivores. In: Ecology of tropical savannas—Huntley BJ, Walker BH, eds. (1982) Berlin (Germany): Springer. 359–404.

Owen-Smith N, Cumming DHM. Comparative foraging strategies of grazing ungulates in African savannah grasslands. In: The XVII International Grasslands Congress, New Zealand Grasslands Association (1993) Palmerston North, New Zealand. 691–697.

Owen-Smith N, Novellie P. What should a clever ungulate eat? Am Nat (1982) 119:151–178.[CrossRef][Web of Science]

Page RDM. TREEVIEW: an application to display phylogenetic trees on personal computers. Comp Appl Biosci (1996) 12:357–358.[Medline]

Pettorelli N, Vik JO, Mysterud A, Gaillard J-M, Tucker CJ, Stenseth NC. Using the satellite-derived NDVI to assess ecological responses to environmental change. Trends Ecol Evol (2005) 20:503–510.[CrossRef][Medline]

Pfeifer S. Flehmen and dominance among captive adult female scimitar-horned oryx (Oryx dammah). J Mammal (1985) 66:160–163.[CrossRef][Web of Science]

Pinzon J. Using HHT successfully to uncouple seasonal and interannual components in remotely sensed data (2002) Orlando (FL): SCI. 2002.

Pinzon J, Brown ME, Tucker CJ. Satellite time series correction of orbital drift artifacts using empirical node decomposition (2004) Hilbert-Huang Transform: Introduction and Applications: Chapter 10, Part II.

Price SA, Bininda-Emonds ORP, Gittleman AL. A complete phylogeny of the whales, dolphins and even-toed hoofed mammals (Cetartiodactyla). Biol Rev (2005) 80:445–473.[Medline]

Prins HHT. Ecology and behaviour of the African buffalo: social inequality and decision making (1996) New York: Chapman & Hall.

Rubenstein DR, Wikelski M. Seasonal changes in food quality: a proximate cue for reproductive timing in marine iguanas. Ecology (2003) 84:3013–3023.[CrossRef][Web of Science]

Rutberg AT. Birth synchrony in American bison (Bison bison)—response to predation or season. J Mammal (1984) 65:418–423.[CrossRef][Web of Science]

Rutberg AT. Adaptive hypotheses of birth synchrony in ruminants—an interspecific test. Am Nat (1987) 130:692–710.[CrossRef][Web of Science]

Ryan SJ. Spatial ecology of African buffalo (Syncerus caffer) and their resources in a savanna landscape [PhD thesis] (2006) Berkeley (CA): University of California at Berkeley.

Ryan SJ, Knechtel CU, Getz WM. Range and habitat selection of African buffalo in South Africa. J Wildl Manag (2006) 70:764–776.[CrossRef]

SAS. JMP version 5 (Academic) (1989 –2003) 5.1 ed. Cary (NC): SAS Institute Inc.

Scanlon TM, Caylor KK, Manfreda S, Levin SA, Rodriguez-Iturbe I. Dynamic response of grass cover to rainfall variability: implications for the function and persistence of savanna ecosystems. Adv Water Res (2005) 28:291–302.[CrossRef]

Sekulic R. Seasonality of reproduction in the sable antelope. East Afr Wildl J (1978) 16:177–182.

Sinclair ARE. The African buffalo: a study of resource limitation of populations (1977) Chicago (IL): University of Chicago Press.

Sinclair ARE, Mduma SAR, Arcese P. What determines phenology and synchrony of ungulate breeding in Serengeti? Ecology (2000) 81:2100–2111.[CrossRef][Web of Science]

Skinner JD, Dott HM, Matthee A, Hunt L. Captive breeding of the white rhinoceros, Ceratotherium simum, and the Cape buffalo, Syncerus caffer. Onderstepoort J Vet Res (2006) 73:237–239.[Web of Science][Medline]

Testa JW. Does predation on neonates inherently select for earlier births? J Mammal (2002) 83:699–706.[CrossRef][Web of Science]

Tucker CJ, Pinzon JE, Brown ME, Slayback D, Pak EW, Mahoney R, Vermonte E, El Saleous N. An extended AVHRR 8-km NDVI dataset compatible with MODIS and SPOT vegetation NDVI data. Int J Remote Sens (2005) 26:4485–4498.[CrossRef]

Turner W. Activity patterns of male buffalo [masters dissertation] (2003) Johannesburg (South Africa): University of Witwatersrand.

Turner WC, Jolles AE, Owen-Smith N. Alternating sexual segregation during the mating season by male African buffalo (Syncerus caffer). J Zoolog (2005) 267:291–299.

Visscher DR, van Aarde RJ, Whyte I. Environmental and maternal correlates of foetal sex ratios in the African buffalo (Syncerus caffer) and savanna elephant (Loxodonta africana). J Zoolog (2004) 264:111–116.

Wehausen JD. Fecal measures of diet quality in wild and domestic ruminants. J Wildl Manag (1995) 59:816–823.[CrossRef]

Wofford H, Holechek JL, Galyean ML, Wallace JD, Cardenas M. Evaluation of fecal indices to predict cattle diet quality. J Range Manag (1985) 38:450–454.


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?



This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Lay Summary
Right arrow All Versions of this Article:
18/4/635    most recent
arm028v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (5)
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Ryan, S.
Right arrow Articles by Getz, W.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Ryan, S.
Right arrow Articles by Getz, W.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?