Introduction
Over the last 40 years unprecedented land cover and land-use changes have occurred in sub-Saharan Africa, driven mainly by the conversion of woodlands into croplands, shrublands, and pasture (Lambin et al. Reference Lambin, Geist and Lepers2003, Brink and Eva Reference Brink and Eva2009). Rates of land conversion have been particularly high in West Africa’s Sudano-Sahelian savannas (Paré et al. Reference Paré, Söderberg, Sandewall and Ouadba2008, Brink and Eva Reference Brink and Eva2009, Ouedraogo et al. Reference Ouedraogo, Savadogo, Tigabu, Dayamba and Odén2010), where some of the highest human population densities and growth rates on the continent have been reported (UN World Population Prospects 2010). The resulting habitat impoverishment, caused by loss of nest sites and prey populations that depend on vegetation cover, has been suggested to have contributed significantly to the decline of raptor populations in this region (Thiollay Reference Thiollay2006, Reference Thiollay2007). In light of this immense anthropogenic pressure, the situation in western Africa offers valuable insight into the ecology, ranging and behaviour of raptors confronted with such conditions, which may be representative for future conditions in many parts of the continent.
Multiple studies have examined nest site selection and diet of Afrotropical raptors in habitats modified by human activities (Machange et al. Reference Machange, Jenkins and Navarro2005, Bamford et al. Reference Bamford, Monadjem and Hardy2009, Ogada and Kibuthu Reference Ogada and Kibuthu2009, Virani and Harper Reference Virani and Harper2009). No studies, to our knowledge, have examined foraging habitat selection at patch level scales by raptors in human-transformed landscapes. However, a thorough understanding thereof is a prerequisite to effective conservation, since maintenance of populations inside protected reserves alone may not guarantee population persistence (Chape et al. Reference Chape, Harrison, Spalding and Lysenko2005, Clerici et al. Reference Clerici, Bodini, Eva, Gregoire, Dulieu and Paolini2006). Raptors concentrate foraging in habitat patches where food availability and energy returns are highest (Nishimura Reference Nishimura1991, Reference Nishimura1994, Ward et al. Reference Ward, Gutiérrez and Noon1998), but morphological, ecological and behavioural adaptations underlie differential selection at the patch level (Thiollay and Clobert Reference Thiollay and Clobert1990). Sedentary raptors occupy year-round territories, probably heightening their sensitivity to anthropogenic land-use compared to more mobile species (Buij et al. Reference Buij, Croes and Komdeur2013a). This may lead to their disappearance or decline from intensely exploited, overgrazed areas with little vegetation cover and impoverished prey populations (Herremans and Herremans-Tonnoeyr Reference Herremans and Herremans-Tonnoeyr2000). Cultivated fields may provide abundant prey (e.g. rodents; Buij et al. Reference Buij, Van der Goes, Iongh, Gagare, Haccou, Komdeur and De Snoo2012), but can also lead to shifts in prey abundance and diet, with potential consequences for long-term productivity of sedentary raptors (Ogada and Kibuthu Reference Ogada and Kibuthu2009).
The Dark-chanting Goshawk Melierax metabates is a medium-sized, generalist diurnal raptor that inhabits West Africa's Sudano-Sahelian savannas, where it mainly still-hunts predominantly terrestrial vertebrate prey, most often from a prominent perch such as a tree (Steyn Reference Steyn1982). Like other sedentary raptors (Thiollay Reference Thiollay1978, Thiollay and Clobert Reference Thiollay and Clobert1990), the Dark-chanting Goshawk - with its high wing load, long tail and long legs - is morphologically adapted to capturing agile, often terrestrial prey in thick cover. Since human populations have appropriated an increasing proportion of scarce natural resources in agro-ecosystems, it seems likely that many areas have become unprofitable for Dark-chanting Goshawks due to scarcity of food. To examine the impact of habitat transformation on foraging in this Afrotropical raptor, we investigated differences in microhabitat selection, foraging effort and energy returns between male Dark-chanting Goshawks inhabiting natural and transformed savannas in northern Cameroon. We radio-tracked male Dark-chanting Goshawks in natural and transformed habitats, to determine foraging effort (time spent on perches and distance covered between perches) and energy returns. We examined microhabitat selection based on the premise that patch utilisation patterns reflect patch quality as in other sit-and-wait raptors (Nishimura Reference Nishimura1991), with prey-rich patches being visited more frequently and for longer periods than other patches. We predicted that (1) apart from perch availability, the amount of herbaceous layer and shrub cover were the most important factors determining probability of foraging-patch use by foraging Dark-chanting Goshawks, given the importance of vegetation cover for small vertebrate prey populations in arid savannas (e.g. Castellano and Valone Reference Castellano and Valone2006); (2) foraging patches are comparatively rare (i.e. further apart and of poorer quality) in transformed compared to natural habitat, which would necessitate larger ranges, longer inter-perch flights and greater foraging effort in transformed habitat; and (3) differences in diet composition may result in lower energy intake in transformed compared to natural habitats, which together with lower foraging effort in the latter, will give rise to significantly greater foraging returns per effort in natural compared to transformed habitats.
Methods
Study area
The study area (10°40’N–11°24’N and 14°15’E–14°50’E) was located inside the southern half of Waza National Park and the cultivated zone south of the park in the far north of Cameroon (Fig. 1). The degree of human influence on the habitat in the protected area (e.g. through livestock grazing, woodcutting; Scholte Reference Scholte2003) was negligible compared to the unprotected area. We therefore refer to Waza N.P. as “natural habitat” vs. “transformed habitat” for the unprotected area. The area was characterised by a Sudano-Sahelian climate, with a short rainy season (annual rainfall c.500 mm, mainly July-September), and woodland dominated by marula Sclerocarya birrea with African birch Anogeissus leiocarpus and desert date Balanites aegyptiaca, locally interspersed with dense clusters of red acacia Acacia seyal. The transformed habitat consisted of a mosaic of cultivation (10–80% of land surface; notably dry-season sorghum but also maize, tomatoes, groundnut, onions), human settlements (< 5%), and fragments of heavily exploited woodlands (5–35%), under severe pressure by livestock grazing and slash-and-burn activities for cultivation and wood harvesting (see also Buij et al. Reference Buij, Folkertsma, Kortekaas, Iongh and Komdeur2013b).
Location of territory centers of the 12 radio-tagged Dark-chanting Goshawk males in northern Cameroon.
Radio-tracking
We captured 12 adult (> 2 years old) male Goshawks in the natural (n = 6) and transformed habitats (n = 6) between 8 November 2009 and 9 February 2010. We repeatedly drove a predefined road network of 290 km in transformed habitat and 80 km in natural habitat and considered each detected Dark-chanting Goshawk for capture, until 12 males had been caught (see Fig. 1 for locations of tagged males). We selected adult, mated males during the middle of the dry season (the non-breeding season), to minimise variability due to sex, age, and differences in reproductive status, and because wet-season conditions prohibited the proposed monitoring. Mated males were distinguished from unmated males by the presence of adult females in their direct vicinity (< 300 m); the presence of a female could usually be determined within one hour of observation. Dark-chanting Goshawks were trapped with a dho-gaza net or bal-chatri trap (Bloom et al. Reference Bloom, Clark, Kidd, Bird and Bildstein2007), ringed, weighed to the nearest 0.1 g and equipped with a 10-g TW-3 single celled transmitter (Biotrack Ltd, Dorset, UK) mounted to the central tail feathers, and conforming to general safety precautions (Walls and Kenward Reference Walls, Kenward, Bird and Bildstein2007).
All twelve tagged Dark-chanting Goshawks were tracked for three consecutive days from sunrise to sunset (c.11.5-12 h) between 5 December 2009 and 1 March 2010, from a vehicle by three observers using a 3-element Yagi antenna and 4 MHz Sika receiver, binoculars and a telescope. Observation conditions were ideal during each 3-day tracking session, with clear skies and no sandstorms. Observers remained inside the vehicle and the distance to the bird was kept at a minimum of 70 m or further if visual contact could be maintained. Monitoring was alternated between natural and transformed habitats to minimise the potentially confounding effects of changes in vegetation, climate, and prey parameters as the dry season progressed. Dark-chanting Goshawk tracking sessions were followed by vegetation surveys in the bird’s range, after which the next Dark-chanting Goshawk was tracked. We compiled behaviour protocols and noted the time spent (in seconds) on each of the following five behavior categories: (1) perching: bird sitting in tree, other structures or on ground, while not engaged in feeding, (2) inter-perch flight: active flight between different perches, (3) circling: circling flight with little or no wing movement, (4) feeding: consuming prey, equivalent to meal duration for individual prey items, and (5) foraging on foot: bird walking on the ground. We also noted strikes on prey, i.e. rapid movements directed towards terrestrial, or rarely, avian prey. We distinguished three categories: 1) successful strikes, when prey was captured, 2) unsuccessful strikes, when a strike at prey did not result in capture, and 3) strikes with an unknown outcome (20.5% of 337 strikes; on average 5.7 and 5.8 strikes per bird in respectively transformed and natural habitats), i.e. strikes made at terrestrial prey but herbaceous vegetation obstructed a clear view of the bird on the ground, making it impossible to determine success and prey intake. Capture attempts classified to the latter category invariably involved small terrestrial prey (lizards, insects), which if captured was consumed on the ground, whereas large prey was always consumed above ground level. An additional 18 strikes (1.5 ± 0.63 SE strikes per bird) were classified into a fourth category: movements toward the ground that were rapidly aborted, i.e. insufficient to be classified as a strike at prey, or possibly belonged to categories 2 or 3 but could not be classified with certainty; these were excluded from calculations of strike frequency.
Incremental analysis during a pilot study based on movements of a radio-tracked male Dark-chanting Goshawk indicated that 95% of the outer minimum convex polygon (MCP) range (Harris et al. Reference Harris, Cresswell, Forde, Trewhella, Woollard and Wray1990) was attained based on all locations after three days compared to the asymptotic MCP-estimate reached at day five. Therefore, we assumed a three-day survey would provide a reasonable representation of short-term range size. We recorded the GPS coordinates of each perch after every displacement, and the tree height and the altitude above the ground of the perch were estimated to the nearest 0.5 m after the Dark-chanting Goshawk had left. Relocation and identification of the perch site were fairly straightforward due to the openness of the habitat, although birds occasionally remained out of sight (e.g. when circling). Linear flight is the most energy consuming activity of raptors (Mosher Reference Mosher1976); thus, to assess foraging effort, the length of linear displacements was measured using the home-range extension tool in ArcView GIS 3.2 software (Esri 1999).
Habitat measurements
The vegetation structure at six patches used by foraging Dark-chanting Goshawks was compared with six nearby unused patches within the MCP range of the bird (i.e. third-order selection; Johnson Reference Johnson1980), a method successfully used to provide insight into the importance of habitat elements for foraging raptors (Kenward Reference Kenward1982, Beier and Drennan Reference Beier and Drennan1997). The foraging patch or “hunting plot” was defined as a circle, with a perch at the centre that met with the following criteria: at least one capture attempt was executed from the perch during the three-day observation period, and the cumulative time spent on the perch was longer than the time spent on other perches used for foraging. We used the average strike distance from the perch at which Dark-chanting Goshawks attempted to capture prey as the radius of the hunting plot. This average strike distance (45 m ± 6.0 SE; range: 7–108 m) was based on 24 capture attempts by seven Dark-chanting Goshawks equally divided between habitats, obtained during a pilot phase prior to the study. Strike distances were comparable in length in natural (47.7 m ± 9.2 SE, n = 12) and transformed habitats (42.5 m ± 8.0 SE, n = 12; U = 67.0, P = 0.77). For each hunting plot, we allocated a paired control plot in a fixed direction (north) of the hunting plot. To be sure that we measured an unused area, the control plot was located > 90 m (i.e. twice the average strike distance) from all other used perches and > 216 m (twice the maximum observed strike distance) from the centre of each of the five other selected hunting plots. The centre of a control plot had to meet the criteria of a potentially suitable foraging perch, i.e. a woody plant of at least 3 m, representing the minimum perch height from which a successful strike at prey was observed.
For each of the 72 paired hunting and control plots, habitat characteristics (see Table 1) were sampled on the day following the three-day monitoring period. These features were hypothesised to determine prey abundance (tree density, shrub cover, herbaceous layer cover, agriculture cover) and prey accessibility (perch availability: tree height, tree density; prey accessibility: herbaceous layer cover and height) for perch-hunting Dark-chanting Goshawks, based on published information (Thiollay Reference Thiollay1978, Thiollay and Clobert Reference Thiollay and Clobert1990). Features of the rooted herbaceous layer were measured on 4 x 4 m squares located at 4-m intervals along four 45-m transects positioned towards the north, east, south and west and radiating out from the centre of the plot, thus covering a proportionally greater amount of habitat near the perch than further from it (following Beier and Drennan Reference Beier and Drennan1997). The woody plant cover was subdivided into three categories representing the woody plant strata (shrubs < 3 m, trees: > 3 m, > 6 m). The canopy dimensions of all shrubs in the plots were measured with a measuring pole, and the percentage cover calculated relative to the surface area of the plot. The plot was divided into four equal sections and the total surface area under cultivation estimated for each section; these areas were summed to calculate the percentage agriculture cover within the plot. All trees > 3 and tall trees > 6 m inside the plot were counted to estimate their density and the height of all tall trees was measured to the nearest 0.5 m.
Vegetation parameters measured on hunting plots and paired control plots.
1 We used a measuring pole with a 20-cm radius circular disc dropped onto the vegetation sensu Holmes (Reference Holmes1974).
Habitat characteristics were also examined at territory level to allow comparisons of vegetation composition between territories in natural and transformed habitats. To achieve this, we positioned a 1-km transect with a fixed band width of 200 m along the territory’s longest axis using ArcView GIS 3.2, adding a second transect perpendicular to the first if the longest axis did not attain 1 km. We determined tree density along the 1-km transect and shrub and agriculture cover along a 200-m subsection, and herbaceous layer cover and height on sixty 4x4-m squares as described for the plots.
Gross energy intake
We found that Dark-chanting Goshawks consumed the entire prey except for one occasion (1.6%; n = 63 prey items) where a snake was only partially eaten. We therefore assumed that gross energy intake corresponded with the energy value of the entire prey, henceforth called “meal energy value”. To enable estimation of prey body mass, we attempted to identify prey species and genus or family when identification at genus level was not possible. We estimated prey dimensions through comparison with the known length of the Dark-chanting Goshawk tarsus during feeding with the use of a telescope (Marti et al. Reference Marti, Bechard, Jaksic, Bird and Bildstein2007). Estimated body dimensions were snout-vent length (SVL) for lizards and snakes and body width and length, excluding tail, for small mammals. Insects were classified into three size categories: small, < 2 cm; medium, 2-4 cm; large: > 4 cm. To convert body size of consumed lizards (Scincidae, Agamidae spp.) into body mass, we developed SVL-body mass regression curves using specimens collected in the study area. For mice, body mass estimates were based on known body dimensions and mass of specimens of different size and family collected in the study area (R. Buij unpubl. data). Mean body mass of grasshoppers and birds was determined from specimens collected in the field. Body mass estimates were related to meal duration (in s) via linear regressions, which can be used as a proxy measure for prey mass (Masman et al. Reference Masman, Gordijn, Daan and Dijkstra1986). Feeding time was strongly related to body mass estimates for small mammals (F 1,8 = 38.5, r2 = 0.82, P < 0.001) and lizards (F 1,29 = 67.0, r2 = 0.71, P < 0.001), supporting the validity of our assessment.
Prey energy content, categorised to taxon and size, was determined for lizards, grasshoppers, birds, and small mammals using bomb calorimetry (kJ/g; Gessaman Reference Gessaman, Giron Pendleton, Millsap, Cline and Bird1987). Representative prey specimens collected in the study area were oven dried to constant mass (± 0.0001 g) before homogenisation. From each prey specimen, pellets (n = 1-4) were prepared (c.1 g) using a pellet press. The pellets were dried before weighing (± 0.0001 g), and combusted in an IKA C5000 adiabatic calorimeter, which was calibrated using a benzoic acid standard (Gessaman Reference Gessaman, Giron Pendleton, Millsap, Cline and Bird1987). Energy content of prey specimens is presented in Appendix S1 in the online supplementary material.
The energy value of snakes of different length was calculated based on body mass-length relationships of whip and grass snakes (Colubridae; Kaufman and Whitfield Gibbons Reference Kaufman and Whitfield Gibbons1975), which were the most commonly captured snakes, and published energy content (Secor and Nagy Reference Secor and Nagy1994). Dry body mass and caloric values for termites and beetles were similarly adopted from literature (Nagy et al. Reference Nagy, Huey and Bennett1984, Williams et al. Reference Williams, Anderson and Richardson1997, Allen Reference Allen1989, Oba et al. Reference Oba, Shintani, Nakamura, Ojika and Inouye2008). To correct for small prey intake of strikes with unknown outcome, we determined the success rate of capture attempts involving small terrestrial prey with known outcome (23%; n = 56), and the average energy value of the prey items involved (22 kJ). Henceforth, a meal energy value of 5.06 kJ (0.23 × 22 kJ) was ascribed to each capture attempt at terrestrial prey with uncertain outcome. Relaxing this assumption to a 50% deviation (-/+) from estimated consumption did not significantly affect the outcome of comparative analyses with regard to gross energy intake. To examine foraging returns per effort, gross energy intake (kJ) was expressed per strike, successful strike, distance covered (km), and the total time a Dark-chanting Goshawk was observed during the three-day track session (h).
Statistical analyses
Vegetation characteristics, average gross energy intake per unit effort (km traveled, h observed), and foraging data were averaged per individual range and compared with independent samples t-tests or Mann-Whitney U-tests. Standard univariate linear and quadratic regression models were used to investigate the relationships between range size and tree density, tree height, shrub cover, and herbaceous cover and height at range-level, considering the entire dataset (n = 12 Dark-chanting Goshawks). To test which habitat features at range-level were related to foraging success, we used linear and quadratic regression models with the mean number of successful strikes per hour as dependent variable and tree density and height, herbaceous layer cover and height within ranges as independent variables.
Dark-chanting Goshawk foraging data were analysed using Generalized Linear Mixed Models (GLMMs) implemented in R (version 2.6.1., function lmer in R package lme4). To test our assumption that hunting plots were indeed sites where the expected returns were greatest, we examined whether time spent on perches was positively related to hunting attempts made from those perches, taking into account inter-individual variation. Hereto, we used GLMMs with a negative binomial distribution and log link function, with time spent on perches as dependent variable, bird ID as random variable, and hunting attempt (yes = 1, no = 0) as categorical explanatory variable.
To examine which habitat variables best explained the probability of hunting patch use by Dark-chanting Goshawks, we constructed mixed effects logistic regression models with probability of use (0 = paired control plot [reference], 1 = hunting plot) as the dependent variable. Based on published information, we explored the following predictor variables: 1) percentage cover of shrubs, which provide refuges for small terrestrial vertebrates (lizards, rodents; Díaz Reference Díaz1992, Blaum et al. Reference Blaum, Rossmanith, Schwager and Jeltsch2007); 2) tree density, or tall tree density, based on positive relationships between arboreal lizard abundance and tree density and size (Cooper and Whiting Reference Cooper and Whiting2000, Meik et al. Reference Meik, Jeo, Mendelsson and Jenks2002); 3) herbaceous layer cover and height (averaged across plots), because ground cover is positively associated with abundance of terrestrial lizards and small mammals in arid savannas (Castellano and Valone Reference Castellano and Valone2006, Buij et al. Reference Buij, Kortekaas, Van Krimpen, Van Wijk, Van Der Zanden, Iongh, Heitkönig, Snoo and Komdeur2013c); 4) the percentage cover of cultivated fields, with often abundant rodent prey (Buij et al. Reference Buij, Van der Goes, Iongh, Gagare, Haccou, Komdeur and De Snoo2012); and 5) maximum tree height, because tree height is an important determinant of prey accessibility for perch-hunting raptors in savannas (Thiollay and Clobert Reference Thiollay and Clobert1990). We fitted models using all possible combinations of main effects and included interaction terms between tree height and herbaceous layer cover (respecting the marginality principle for modeling; Nelder Reference Nelder1977), because we expected that Dark-chanting Goshawks adapted perch height in response to differences in herbaceous layer cover (cf. Thiollay and Clobert Reference Thiollay and Clobert1990), affecting their patch use. Bird ID was included as random variable in the above-mentioned models to account for inter-individual variation. Before modelling, we verified collinearity between predictor variables (r > 0.7) using Spearman rank correlation tests and retained only those variables that had the lowest corrected Akaike information criterion (AICc) value from univariate models of the two predictors (Burnham and Anderson Reference Burnham and Anderson2002); this meant that tall tree density and herbaceous layer height were excluded from models. We calculated AICc to obtain the model that best fit the data and considered the following levels of empirical support for models: ∆ AICc < 2 substantial support, ∆ AICc > 4 considerably less support (Burnham and Anderson Reference Burnham and Anderson2002). We used the dredge function in R ( part of R package MuMIn, version 1.9.5; Barton Reference Barton2009) to estimate the relative importance of each predictor based on the Akaike weights (w i).
Differences in time budgets between habitats were examined using GLMMs with negative binomial distribution and log link function with the time spent on perching, inter-perch flight, circling, and foraging on foot as dependent variables. Time spent was defined here as the time period spent on a particular behaviour until the bird switched to another behaviour. We constructed a second set of models with the relative contribution of each behavioural category to total observation time, i.e. the percentage time spent per day on each behavioral category, as dependent variables. Habitat (natural/transformed) was included as categorical predictor variable in these models and bird ID as random variable. To test for differences in the length of individual linear displacements between habitats, GLMMs with negative binomial distribution and log link were used with bird ID as random variable, the length of displacements as dependent variable, and habitat type as explanatory variable. For differences between habitats in minimum inter-perch distance, i.e. the distance of each structure (tree, electricity pole) used as a perch to the nearest neighbouring structure used as a perch at any time during the three-day track period, GLMMs with normal distribution and identity link function were used, with bird ID as a random variable, the log transformed inter-perch distances as dependent variable, and habitat type as the categorical explanatory variable.
To test for associations between prey choice and habitat type, we used GLMMs with a multinomial distribution and log link function to examine associations between habitat (natural/transformed) and three prey categories (reptiles, rodents, and infrequent – and therefore lumped – “other prey” such as birds and insects). Prey type was included as dependent variable, habitat as categorical fixed effect, and bird ID was added as a categorical random variable to account for inter-individual differences in prey choice. To test for differences in gross energy return (rounded to kJ) per successful strike between habitats, we used GLMMs with negative binomial distribution and log link function with meal energy value as dependent variable, bird ID as a random variable, and habitat type as categorical explanatory variable. We used the same approach to test for differences in gross energy return per strike. All means are given ± SE.
Results
Habitat features of Dark-chanting Goshawk ranges
Mean range size was greater in transformed compared to natural habitat (Table 2). Comparisons of vegetation structural characteristics between habitat types at the territory level indicated on average higher overall tree densities and higher tall-tree densities in natural than transformed habitat (Table 3). Territories in natural habitats also had a higher herbaceous layer cover, but a lower shrub cover than in the transformed habitat. Maximum tree height tended to be greater in natural habitats. The proportion of herbaceous layer cover was comparable between habitats. MCP-range size was not significantly related to any measure of vegetation structure at the territory-size scale, also when relationships were considered separately for each habitat type. However, range size tended to decline with increasing tree density inside ranges (F 1,10 = 3.56, r2 = 0.26, P < 0.10).
Foraging statistics (mean ± SE) for 12 male Dark-chanting Goshawks in natural and transformed habitats between December 2009 and March 2010. Asterisks indicate significant differences between habitats with independent samples t-test; *P < 0.05.
1 Includes successful strikes (n natural = 20, n transformed = 43), unsuccessful strikes (n natural = 87, n transformed = 118), and strikes with unknown outcome (n natural = 35, n transformed = 34).
Vegetation features (mean ± SE) for territories of 12 male Dark-chanting Goshawks in natural and transformed habitats in northern Cameroon, December 2009-February 2010. Vegetation parameters were recorded along a 1-km transect with a fixed band width of 200 m centered on the territory. Asterisks indicate significant differences between habitats with independent samples t-test or Mann-Whitney U-tests; ***P < 0.01, **P < 0.05, *P < 0.10.
Comparison of hunting plots and paired control plots
Dark-chanting Goshawks spent more time at perches from which capture attempts were observed (41.7 ± 3.93 min) than at perches from which no attempts were observed (12.5 ± 1.03 min; F 1,869 = 162.4, P < 0.001), confirming our assumption that they spent more time in patches that had larger potential returns than in patches with lower potential returns. Table 4 summarizes the results of generalized mixed models of the probability of use of foraging patches by Goshawks, based on the habitat predictors. All top models with comparable support (∆AICc < 2) in both transformed and natural habitats included a maximum tree height effect on foraging patch use. Inspection of the parameter estimates for the best-supported models suggested that the probability of use of foraging patches by Dark-chanting Goshawks increased with maximum tree height in natural ( β = 0.241 ± 0.068) and transformed habitats ( β = 0.235 ± 0.078).
Mixed effects logistic regression models comparing probability of foraging patch-use with habitat features within territories of six male Dark-chanting Goshawks in natural and transformed habitats (n = 36 hunting plots vs. n = 36 paired control plots in each habitat type). All possible models with main effects were fitted: herbaceous layer cover (herb), shrub cover (shrub), agriculture cover (agri), tree density (treedens), and maximum tree height (treeheight), and the interaction of tree height with herbaceous layer cover. Bird ID (1-12) was added as a random effect. Only parameter estimates for models with ∆AICc < 4 from the best model are shown, in order of increasing ∆AICc, with models ∆ AICc < 2 having substantial support (Burnham and Anderson Reference Burnham and Anderson2002). Superscripts with estimates of intercept and effect size indicate significance levels: P < 0.10 (+); P < 0.05 (*); P < 0.01(**); and P < 0.001(***). The Akaike weight w i, i.e. the weight of evidence in support of the model, was calculated relative to all models.
Among models of high to moderate importance in transformed habitat (∆AICc < 4), those that included a maximum tree height effect on patch use comprised 97% of the total model weight. Herbaceous layer cover was another important predictor, included in 15 models within ∆AICc < 4 that together comprised 67% of total model weight. Models with shrub cover and agriculture cover each contributed 51% to the total model weight for ∆AICc < 4 models, those with an interaction between herbaceous layer cover and tree height 25%, and those with tree density 22%. The top models (∆AICc < 2; Table 4) provided support for positive relationships between foraging patch use and agriculture cover, shrub cover, and herbaceous layer cover,: this suggests that Dark-chanting Goshawks concentrated their foraging effort in patches where agriculture and ground vegetation cover were comparatively high. Three of these best-supported models also included an interaction effect between herbaceous layer cover and tree height, but there was little evidence that this affected patch use.
In natural habitat, models with a tree height effect comprised 82% of the total model weight for models of ∆AICc < 4. Models with shrub cover and tree density contributed 20% and 19% respectively to total model weight, and those with herbaceous layer cover 13%. These results suggest that habitat covariates differed in their importance between habitat types: in particular, herbaceous layer cover was a relatively more important predictor of foraging patch use in transformed than in natural habitats. The top models (∆AICc < 2; Table 4) suggested that foraging patch use differed with shrub cover inside patches, similar to transformed habitat, and with tree density; but evidence that these covariates affected patch use was weak.
Movements and foraging effort
Despite larger ranges in transformed habitat, the length of single linear displacements did not differ between natural (277 ± 9.06 m) and transformed habitats (287 ± 11.3 m; F 1,1337 = 0.17, P = 0.68). Relative and absolute time spent on perching, inter-perch flight, circling, and foraging on foot did not differ between habitats. Similarly, foraging efforts, i.e. the distance covered per hour and the number of strikes at prey per hour, did not differ between habitats (Table 2). There was also no difference in the number of perches used between territories in natural and transformed habitat (Table 2). But the minimum inter-perch distance, i.e. the distance of each perch to the nearest neighboring perch, was greater in transformed (89.9 ± 3.17 m) than natural habitat (63.5 ± 2.08 m; F 1,805 = 17.5, P < 0.001). The 1.4 times greater inter-perch distance in transformed habitat compared to natural habitat and the comparable strike distances in both habitats, suggests that Dark-chanting Goshawks in the natural habitat more often searched overlapping areas, i.e. they exploited their range more intensively than those in the transformed habitat.
Diet and gross energy intake
Body mass of Dark-chanting Goshawks did not differ between habitats and the same was true for a proxy measure of wing load (body mass/wing length; Table 2), suggesting that body condition and energy requirements of flight did not differ between individuals in natural and transformed habitats. We found no association between habitat type and the contribution of reptiles (70% of identified prey items recorded; n = 63), rodents (19%), and other prey (11%; birds and insects; F 2,59 = 0.46, P = 0.63) to the diet (Fig. 2). Despite the absence of dietary differences in terms of broad prey categories, the proportion of Agamidae spp. among reptile prey items identified to family level was higher in natural (33.3%; n = 12) than transformed habitat (9.4%; n = 32); whereas on average smaller lizards (Scincidae, Lacertidae, Gekkonidae spp.) were more commonly caught in transformed (71.9%) than natural habitat (41.7%). No such difference was noted for rodents (Muridae spp.). The gross energy return measured by effort or time did not differ between natural and transformed habitats (Table 2). Gross energy return per strike (successful and unsuccessful strikes combined) was also comparable between natural (32.8 ± 7.67 kJ/strike) and transformed habitats (26.1 ± 4.34 kJ/strike; F 1,266 = 1.19, P = 0.28). Dark-chanting Goshawks in natural and transformed habitat did differ by a higher frequency of successful strikes in the latter (Table 2). On the other hand, gross energy return per successful strike was higher in natural (153.2 ± 25.9 kJ/strike) than transformed habitat (79.2 ± 10.9 kJ/strike; F 1,61 = 7.06, P < 0.05), as a result of a higher meal energy value (especially of reptiles) in natural compared to transformed habitats (Fig. 3). The frequency of successful strikes declined with herbaceous layer height and tree density within ranges (Table 5), but was unrelated to herbaceous layer cover.
Percentage contribution of the broad prey categories to the diets of six Goshawks in natural (n = 20 prey items) and six Goshawks in transformed habitats (n = 43 items).
Box plots showing the distribution of energy values of reptile (n = 44) and rodent prey items (n = 12), for natural and transformed habitat. Medians are denoted by solid black lines while the top and bottom box edges denote the first and third quartile, whiskers denote the largest and smallest data, and outliers as separate points. Prey energy value differed significantly between habitats for reptiles (U12,32, z = -2.48, P < 0.05), but not for rodents (t10 = -1.07, P = 0.31).
Relationship between mean strike success rate (number successful strikes/hour) by Dark-chanting Goshawks (n = 12) and tree density and herbaceous layer height within ranges.
Discussion
Contrary to expectations, we did not detect a significant shift in the contribution of broad prey categories to the diet of Dark-chanting Goshawks in response to anthropogenic modification of their savanna habitat. Dark-chanting Goshawks in the human-transformed habitat neither adjusted their activity patterns, nor did they invest more energy in foraging. Rather, in transformed habitat higher prey capture rates compensated for the less profitable prey items (prey items with a lower energy value). Dark-chanting Goshawks in natural habitat were limited in their strike success rate by a relatively thick herbaceous layer and high tree density. Thus, anthropogenic land-use appeared to improve prey accessibility for Dark-chanting Goshawks, through tree clearance and a reduction of grass height resulting from livestock grazing, but reduced the quality of individual prey items. Although prey size has been found to vary with land-use in other raptors (Rodríguez et al. Reference Rodríguez, Johst and Bustamante2006, Cardador et al. Reference Cardador, Planas, Varea and Mañosa2012), to our knowledge this is the first study that demonstrates a deterioration of prey quality with land transformation in an African raptor. These results seem indicative of constraints on body mass of small vertebrates in transformed habitat, assuming that Dark-chanting Goshawks will attempt to maximise their foraging returns by focusing on the most profitable prey available to them (Stephens and Krebs Reference Stephens and Krebs1986). Such differences in body mass may have been imposed by various aspects of land use; for example, terrestrial vertebrates such as lizards were found to be more active where grazing pressure is high (Wasiolka et al. Reference Wasiolka, Blaum, Jeltsch and Henschel2009), which may negatively impact their condition and body mass (Amo et al. Reference Amo, Lopez and Martin2007). In addition, prey species composition or accessibility may change following human exploitation leading to structural vegetation changes of savannas (Meik et al. Reference Meik, Jeo, Mendelsson and Jenks2002), and larger, more profitable prey species might be more commonly encountered in natural habitat. Indeed, we found indications that Dark-chanting Goshawks captured proportionally more Agamid lizards, which are relatively heavy-bodied, in natural habitat than in transformed habitat, implying that high-quality, larger reptile species may be more readily available in natural habitats.
Forage-patch use by Dark-chanting Goshawks in natural and transformed habitats was strongly and positively influenced by the presence of tall trees, which provide cover and perches from which to hunt. In addition, we found evidence that tree density, herbaceous layer, shrub, and agriculture cover influenced patch use, although the importance of these covariates differed between habitats. Foraging patch selection may shift during the wet season, when vegetation cover and prey abundance and activity increase (Thiollay Reference Thiollay1978), but dry conditions prevail during 8 months of the year, suggesting that our results are representative for most of the year when resources are limited. Although we did not measure prey densities, patch-level selection for a well-developed ground vegetation layer in transformed habitats appears to have been primarily related to optimisation of prey abundance, or quality, within foraging patches; for example, shrubs and a well-developed grass layer provide important shelter, foraging, and thermal regulation sites for lizards and rodents (Díaz Reference Díaz1992, Vasquez et al. Reference Vasquez, Ebensperger and Bozinovic2002, Castellano and Valone Reference Castellano and Valone2006). In heavily grazed areas, patches of shrubs and grasses may become especially important refuges for terrestrial prey, which may explain their importance to foraging Dark-chanting Goshawks. Similarly, Dark-chanting Goshawks may have selected sorghum and millet fields for the comparatively high abundance of rodents associated with them (Buij et al. Reference Buij, Van der Goes, Iongh, Gagare, Haccou, Komdeur and De Snoo2012).
Our results suggest that the persistence of Dark-chanting Goshawk populations in agro-ecosystems is based on a trade-off between prey quality and prey capture rates, governed by anthropogenic activities such as livestock grazing and tree clearance. Goshawks may only be able to persist where profitable, exploitable foraging patches are sufficiently common, or where it is possible to increase territory size to meet nutritional demands (Newton Reference Newton1979, Kenward Reference Kenward1982). We found that Dark-chanting Goshawks required larger ranges to meet their food requirements in transformed habitats, where foraging patches were more widely spaced compared to natural habitats, despite the positive effect that some habitat modifications may have had for foraging. Our results further suggest that a moderate degree of human land-use, with livestock grazing reducing sward length, an increase in cultivation of food crops, openness, and shrubs, and the often selective preservation of tall trees, might be beneficial by maintaining or positively affecting prey accessibility and availability. Dark-chanting Goshawks are likely to decline, however, in areas where extensive vegetation clearance reduces the capacity of the habitat to support prey populations, or their exploitability. Such vegetation denudation is a continuing process in West Africa, which is likely to have contributed to the sharp decline of Dark-chanting Goshawks during the past four decades, as well as other sedentary raptors with comparable morphological traits and foraging ecology (Thiollay and Clobert Reference Thiollay and Clobert1990, Thiollay Reference Thiollay2006). Additional pressure from increased human disturbance at nest sites and direct or indirect effects of pesticide use in croplands may further impact on raptor populations (Buij et al. Reference Buij, Croes and Komdeur2013a). Future examinations are required to confirm this and to determine whether the observed differences in prey quality and range size negatively affect nestling condition, survival, and long-term reproductive output (Rodríguez et al. Reference Rodríguez, Johst and Bustamante2006, Buij et al. Reference Buij, Folkertsma, Kortekaas, Iongh and Komdeur2013b).
Management implications emanating from this study would be to conserve tall trees, pasture and shrubs in cultivated landscapes, thus maintaining adequate food and shelter for terrestrial prey animals exploited by Dark-chanting Goshawks. Unfortunately, land saturation and a general lack of financial means to enable sustainable cropping systems lead to heavily exploited agro-ecosystems increasingly deprived of natural vegetation cover (Lambin et al. Reference Lambin, Geist and Lepers2003, Gray Reference Gray2005). Appropriate agricultural intensification-related policy initiatives (Buij et al. Reference Buij, Croes and Komdeur2013a), encouraged and enforced by state and traditional authorities, including silvicultural practices with designated pastures protected from heavy exploitation, might represent the only alternative to the progressive impoverishment of West Africa’s savannas. Such action may help curb the gradual disappearance of sedentary Afrotropical raptors in the agro-ecosystems that increasingly dominate their distribution range.
Supplementary Materials
The supplementary materials for this article can be found at journals.cambridge.org/bci
Acknowledgements
This study was financially and logistically supported by the Institute of Environmental Sciences of the University of Leiden, and the Centre for Ecological and Evolutionary Studies, University of Groningen, the Netherlands. M. Bechard, M. Murphy, R. Limiñana, H. H. de Iongh, G. de Snoo, and two anonymous reviewers provided useful comments on a draft. G. Malan, S. Thomsett, A. Kemp, R. G. Bijlsma, M. Hammers, W. S. Clarke, and M. Virani provided valuable suggestions. H. Hamadou is thanked for his assistance with logistics.
