Hostname: page-component-586b7cd67f-t7fkt Total loading time: 0 Render date: 2024-11-28T07:16:56.579Z Has data issue: false hasContentIssue false

Response of large galliforms and tinamous (Cracidae, Phasianidae, Tinamidae) to habitat loss and fragmentation in northern Guatemala

Published online by Cambridge University Press:  04 October 2012

Daniel H. Thornton*
Affiliation:
Department of Wildlife Ecology and Conservation, University of Florida, Gainesville, Florida, USA.
Lyn C. Branch
Affiliation:
Department of Wildlife Ecology and Conservation, University of Florida, Gainesville, Florida, USA.
Melvin E. Sunquist
Affiliation:
Department of Wildlife Ecology and Conservation, University of Florida, Gainesville, Florida, USA.
*
(Corresponding author) E-mail [email protected]
Rights & Permissions [Opens in a new window]

Abstract

The potential conservation value of fragmented or countryside landscapes in the tropics is being increasingly recognized. However, the degree to which fragmented landscapes can support species and the key patch and landscape features that promote population persistence remain poorly understood for elusive species such as ground-dwelling birds. We examined the presence/absence of seven species of galliforms and tinamous in 50 forest patches of 2.9–445 ha in northern Guatemala using camera traps and audiovisual surveying. After accounting for differences in detectability among species we found great variation in patterns of vulnerability of these species to habitat loss and fragmentation, with the three largest species being the most vulnerable. Distribution patterns of species among patches was influenced more strongly by measures of landscape context, such as the amount and configuration of habitat in the surrounding landscape, than within-patch variation in vegetation structure or disturbance. Our results indicate that large-bodied game birds may be particularly sensitive to habitat loss and fragmentation and emphasize that management efforts for these species need to go beyond consideration of local, within-patch factors to consider the impact of processes in the surrounding landscape. Our findings also demonstrate the utility of camera traps as a methodology for surveying large terrestrial bird species in fragmented landscapes.

Type
Papers
Copyright
Copyright © Fauna & Flora International 2012

Introduction

Clearing of tropical forest for food production and urban expansion is proceeding at a tremendous rate (Asner et al., Reference Asner, Rudel, Aide, DeFries and Emerson2009), resulting in landscapes that contain small patches of forest surrounded by a matrix of agriculture, cattle pasture and urban lands. Variation exists in the ability of species to use forest patches embedded in human-dominated environments, with some species being much more likely than others to persist within fragmented landscapes (Laurance, Reference Laurance1991, Reference Laurance2002; Gascon et al., Reference Gascon, Lovejoy, Bierregaard, Malcolm, Stouffer and Vasconcelos1999). The identification of which species are most at risk and the characteristics of patches and landscapes that influence the likelihood of persistence within fragmented environments are key conservation issues that need to be resolved for planning effective conservation strategies.

Previous research has shown that large terrestrial and/or game species such as galliforms are vulnerable to hunting pressure and habitat disturbance (Thiollay, Reference Thiollay1999, Reference Thiollay2005) and may disappear or become rare in single isolated forest remnants (Robinson, Reference Robinson1999). However, beyond interview-based occupancy research (Urquiza-Haas et al., Reference Urquiza-Haas, Peres and Dolman2009) little work has focused on understanding the factors influencing distribution of these species across large, highly fragmented landscapes. Given the site-specific nature of the response of birds to fragmentation (Sigel et al., Reference Sigel, Robinson and Sherry2010) more data from a variety of systems are needed to understand better how these species cope with environmental change. Because large birds such as cracids and tinamous often constitute a substantial portion of avian biomass in undisturbed neotropical forests (Terborgh et al., Reference Terborgh, Robinson, Parker, Munn and Pierpont1990; Robinson et al., Reference Robinson, Brawn and Robinson2000) and play important roles as seed predators and seed dispersers (Brooks & Strahl, Reference Brooks and Strahl2000; Érard et al., Reference Érard, Théry and Sabatier2007; Bertsch & Barreto, Reference Bertsch and Barreto2008b) their ability to persist within fragmented landscapes may have ramifications for the long-term trajectory of bird communities and forest composition within forest remnants.

In addition to traditional survey methods camera traps provide another methodology for documenting large-scale distribution patterns of terrestrial birds (O’Brien & Kinnaird, Reference O'Brien and Kinnaird2008). Although camera traps have been used primarily to document the presence or abundance of elusive terrestrial mammals (e.g. Maffei et al., Reference Maffei, Cuellar and Noss2004; Dillon & Kelly, Reference Dillon and Kelly2008; Thornton et al., Reference Thornton, Branch and Sunquist2011a) they can also provide a wealth of information on the occupancy of large ground-dwelling birds, even when those species are not the focal subject of a particular study (O’Brien & Kinnaird, Reference O'Brien and Kinnaird2008). To date, however, the use of camera traps in this regard has been limited (O'Brien & Kinnaird, Reference O'Brien and Kinnaird2008).

We studied the response of three cracids (great curassow Crax rubra, crested guan Penelope purpurascens, plain chachalaca Ortalis vetula), ocellated turkey Meleagris ocellata and three tinamous (great tinamou Tinamus major, slaty-breasted tinamou Crypturellus boucardi, little tinamou Crypturellus soui) to forest loss and fragmentation in northern Guatemala by examining occupancy patterns of these species in forest patches. We assessed the vulnerability of these species to habitat loss and fragmentation and investigated the influence of patch- and landscape-level factors on distribution patterns of species among forest patches.

Study area

We conducted this study on a 250,000-ha area in the Petén region of northern Guatemala (Fig. 1). The northernmost part of the study area is located within the buffer zone of the Maya Biosphere Reserve (a UNESCO world heritage site and the central section of the largest contiguous tropical forest in Central America) and the remainder on private lands further to the south. This area was a formerly contiguous forest but now consists of a diverse mosaic of land-use types (V.H. Ramos, pers. comm.). Forest cover in this area consists primarily of subtropical humid rainforest interspersed with scattered small patches of seasonally inundated bajo forest. Mean temperature in this region is 21–24°C and mean annual precipitation 1,350 mm, with a marked dry season from late December to May when mean monthly rainfall is only 60 mm.

Fig. 1 Location of the 50 sampled patches within the study area (Thornton et al., Reference Thornton, Branch and Sunquist2011b). Inset shows approximate location of the main map in northern Guatemala.

Methods

Site selection

Our study sites were 50 primary rainforest patches (hereafter focal patches) of 2.9–445.5 ha (mean 68.5±SE 98.7 ha). We did not sample rainforest patches that were severely degraded by logging or fire (loss of >25% canopy cover) but almost all patches included in the study were lightly degraded. Because of the size of our overall study area and difficulty of travel in certain areas we did not select patches randomly across the entire study landscape. Instead, we divided the area into five zones that contained multiple patches and surveyed the zones in a random order. Within each zone we randomly selected patches for study until a maximum of 15 patches had been surveyed.

Occupancy surveys

We used camera trapping and audiovisual surveys to detect bird species within each patch. Sampling was from January to mid September. Mid–late wet season sampling (mid September to December) was avoided because of problems with camera-trap performance. Camera traps were set to survey for mid- to large-sized mammals (Thornton et al., Reference Thornton, Branch and Sunquist2011a) but also recorded the presence of all seven bird species considered here. We placed cameras along roads and small and large game trails, and at water holes and in other areas containing substantial signs of use by animals. Although the proportion of cameras placed at each of these sites varied among patches, each patch had at least one camera placed at each type of site. We placed more cameras in larger patches: seven cameras in patches <10 ha, 10 in patches ≥10 and <20 ha, 14 in patches ≥20 and <40 ha, 17 in patches ≥40 and <80 ha, 20 in patches ≥80 and <160 ha, 25 in patches ≥160 and <320 ha, and 28 in patches ≥320 ha. We placed camera traps at least 10 m from the edge of patches, with the sensor c. 10–20 cm off the ground so that smaller species could not avoid detection by walking under the sensor. We used the following minimum spatial separation between cameras so that we would cover a greater area in larger patches: 50 m for patches <20 ha, 100 m for patches <80 ha, and 150 m for all other patches. Cameras were located in an approximate grid formation, with the caveat that we sought out the best locations possible for placement of each camera. We could typically survey 1–3 patches simultaneously (depending on patch size), using 30–40 cameras during each sample period.

We deployed camera traps for a 16-day period in each patch. A photograph of a particular species at any of the cameras within a patch was considered an indication of presence. We recorded presence/absence for each species within each patch after every 4-day interval. By dividing the 16-day period into 4-day sessions we created a series of repeat detection/non-detection data (i.e. a detection history) for use in modelling detection probabilities for each species (MacKenzie et al., Reference MacKenzie, Nichols, Royle, Pollock, Bailey and Hines2006).

We also surveyed patches using audiovisual surveys in the early morning (between sunrise and 3 hours after sunrise), recording bird sightings and vocalizations. We acknowledge that the timing of surveys may not have been ideal for our study species, which are most vocal pre-dawn (Terborgh et al., Reference Terborgh, Robinson, Parker, Munn and Pierpont1990). However, as the overall aim of the audiovisual surveys was to document a wide range of primarily diurnal/crepuscular birds and mammals (Thornton et al., Reference Thornton, Branch and Sunquist2011a) we chose this timing as most appropriate. We sampled a subset (n=17) of patches in pre-dawn hours but, because of the low numbers of birds or mammals detected, these surveys were not continued and are not included in the analysis.

Audiovisual surveys were repeated five times within a 2-week period for each patch, resulting in a series of detection/non-detection data. Surveys were not performed in the rain because of low detectability in such conditions. Because we were recording the presence (not density) of species we did not cut transects but rather sought to walk through as much of the patch as possible to increase detection probabilities. We surveyed patches by walking along small roads, footpaths and game trails, and by walking through sections without any obvious trails. We walked at c. 1 km per hour and recorded direct observations of animals and vocalizations as indications of presence within the patch. For patches <10 ha we were able to walk through most or all of the patch in each session, with a new randomly chosen starting point for each survey. For patches too large to survey completely in one session we divided the site into 2–4 sections and randomly chose a section, a starting point and a compass direction for each session. We attempted to follow compass directions as closely as possible given the terrain, availability of trails and degree of understorey cover. We repeated this process until we had five surveys for the patch. The total distance walked over the five visits within each patch varied positively with patch size: 4 km walked within sites <10 ha, 5 km within sites ≥10 and <40 ha, 7.5 km within sites ≥40 and <80 ha, 9 km within sites ≥80 and <160 ha, and 10 km within sites ≥160 ha. We did not sample patches concurrently with cameras and audiovisual surveys because of concerns that our presence within patches would bias detections on cameras. For most patches we performed audiovisual surveys either 2 weeks before or 2 weeks after camera trapping. However, for a small subset of patches >1 month elapsed between camera trapping and audiovisual surveys because of weather, access and other logistical problems.

Habitat measurements: patch variables

For each focal patch we calculated focal patch size and several indices of within-patch structure and disturbance (Table 1). We estimated the basal area of trees using the point-centred quarter method (PCQM; Cottom et al., Reference Cottom, Curtis and Hale1953) along randomly placed transects within each focal patch. In each quadrat of the PCQM sample we recorded distance to the nearest stem in two size classes: 0–10 cm diameter at breast height (DBH) and ≥10 cm DBH. We also assessed evidence for past fire along the edge and interior of the patch and created a binary variable with 0=no fire effects and 1=fire along edge and/or <25% of interior. We assessed cattle disturbance in the patches in an informal manner by noting observations of cattle, cattle paths, dung and tracks, and created a binary variable with 0=cattle not recently present within the patch and 1=cattle recently present.

Table 1 Patch and landscape variables for the 50 focal patches in the study area in northern Guatemala (Fig. 1).

* Proportion of patches with evidence of fire or cattle disturbance

Habitat measurements: landscape context

We created a vegetation map by performing an unsupervised classification of 2003 Landsat Enhanced Thematic Mapper Plus (ETM+) images with ERDAS Imagine v. 9.0 (Leica Geosystems, Atlanta, USA) that identified four classes: water, pasture/agriculture, regenerating forest ≤15 years old, and primary forest/regenerating forest >15 years old. Primary forest and older regenerating forest were combined into one class because of our inability to separate these classes with Landsat imagery. The overall classification accuracy was 80.3%.

We measured landscape variables within 2-km buffers from the edge of focal patches and used FRAGSTATS (McGarigal et al., Reference McGarigal, Cushman, Neel and Ene2002) to calculate measures of habitat amount and configuration (Table 1). Proportion of forest in the landscape represented habitat amount, and we chose density of patches as our measure of habitat configuration because it was easily interpretable as an index of fragmentation and was not correlated with proportion of forest. We also calculated proportion of regenerating forest as a measure of low quality habitat around patches (Table 1).

We calculated two other metrics reflecting the landscape context of focal patches (Table 1). We calculated distance from the focal patch to the nearest continuous forest site >10,000 ha (i.e. the Maya Biosphere Reserve boundary) as a measure of the potential effect of a nearby source, and distance from the focal patch to the nearest human community as a measure of hunting pressure and/or human disturbance. Finally, we determined the season of sampling (dry vs early wet) for each patch. Correlation coefficients between all continuous potential predictor variables included in the analysis (patch, landscape, and additional variables) were <0.5.

Data analysis

We used detection/non-detection data from the 50 focal forest patches to model patch occupancy and detection probabilities with PRESENCE v. 2.0 (Hines, Reference Hines2006). Based on very limited data from our study species and closely related species we cannot rule out the possibility that some of the larger species may have had home ranges not completely enclosed by the sampling units, or were capable of moving between forest patches (Poulsen, Reference Poulsen1994; Gonzalez et al., Reference Gonzalez, Quigley and Taylor1998; Bertsch & Barreto, Reference Bertsch and Barreto2008a). The occupancy estimator is therefore best interpreted as probability of use of a patch rather than probability of occupancy, and the detection parameter is best interpreted as a combination of two components: the probability that the species was within the patch when sampling occurred and the probability the species was detected, given it was within the patch (MacKenzie et al., Reference MacKenzie, Nichols, Royle, Pollock, Bailey and Hines2006). For ease of presentation we use traditional occupancy terminology.

For six of the seven species we modelled the influence of covariates on detection and occupancy probability. We excluded ocellated turkeys from an analysis of covariate influence because of sparse data (detected in only two patches). However, we did develop a model for ocellated turkeys based on constant probabilities of detection and occupancy to obtain an estimate of proportion of patches occupied for this species. For the remaining species we first determined detection probabilities based on the camera-trapping and audiovisual survey data. The dataset that gave the higher overall detection probability was used in subsequent analyses. We did not combine data from cameras and audiovisual surveys for estimating occupancy and detection (multimethod occupancy modelling; Nichols et al., Reference Nichols, Bailey, O'Connell, Talancy, Grant and Gilbert2008) because camera trapping and audiovisual surveys were not conducted at the same time. We then modelled detection as a function of sampling effort (number of camera-trap nights or number of km walked) and season of survey, holding occupancy constant. The best-fit detection model for each species was determined using the corrected form of Akaike's Information Criterion adjusted for small sample sizes (AICc), which is appropriate when the sample size to parameter ratio is <40 (Burnham & Anderson, Reference Burnham and Anderson2002). This best-fit detection model was then used in all subsequent analyses of occupancy probabilities.

We tested a subset of 29 a priori occupancy models (Appendix) developed from the 11 patch and landscape variables. Candidate models included three models at the patch level: focal patch size, focal patch structure (basal area of small trees/sapling and basal area of large trees), and focal patch disturbance (fire and cattle disturbance). Candidate models also included four models incorporating the effects of landscape context: habitat amount (percentage forest remaining in the buffer around focal patches), habitat configuration (density of forest patches in the buffer), matrix structure (percentage regenerating forest in the buffer), and distance effects (distance to nearest human community and distance to reserve boundary). We tested these patch and landscape models singly and also tested all possible combinations of two models considered simultaneously. We identified best-fit models and best-fit model sets (all models that fell within ∆4 AICc values of the best-fit model) for each species. We did not test more complex models because of sample size limitations. Overall patch occupancy rates were calculated for each species and used as an index of vulnerability to habitat loss and fragmentation. Parameter estimates of covariate influence were calculated based on the best-fit models. For occupancy and detection covariates that appeared in more than one model in the best-fit model sets we calculated model-averaged parameter estimates and unconditional standard errors (Burnham & Anderson, Reference Burnham and Anderson2002). We assessed overdispersion in best-fit models for each species with a goodness-of-fit test procedure implemented in PRESENCE (MacKenzie & Bailey, Reference MacKenzie and Bailey2004). We found no indications of overdispersion for any of the species.

Autocorrelation is problematic for classic statistical tests, such as regression, that rely on independently distributed errors and may lead to erroneous conclusions regarding the significance of covariates in studies of species–environment relationships (Lichstein et al., Reference Lichstein, Simons, Shriner and Franzreb2002). We therefore evaluated model residuals of our best-fit models in each analysis for evidence of spatial autocorrelation using Moran's I correlograms (Fortin & Dale, Reference Fortin and Dale2005). We used the following lag distances for calculation of Moran's I: 0–4, >4–8, >8–12, >12–16 and >16–20 km. A Monte Carlo randomization test with 999 permutations was used to determine the probability of obtaining an I value as large as the observed value in each lag category. Significance was tested using Bonferonni corrections, where the ith lag was tested at 0.05/I (Fortin & Dale, Reference Fortin and Dale2005).

Results

We obtained 521 photographs of the seven species of terrestrial birds in 12,963 camera-trap nights. We also obtained 850 photographs of 25 other bird species. More than half of these photographs were of two terrestrial ground doves (Geotrygon montana and Leptotila verreauxi). Detection probabilities varied among species (0.22–0.59; Table 2). Sampling effort influenced probability of detection for slaty-breasted tinamous, great tinamous and plain chachalacas. Detection probabilities were higher from camera trapping than for audiovisual surveys for great curassow, great tinamou, slaty-breasted tinamou and ocellated turkey, and vice versa for crested guan, little tinamou and plain chachalaca.

Table 2 Estimates of detectability and proportion of patches occupied for the seven large terrestrial birds detected in our study.

1 Estimated proportion of patches occupied not corrected for detectability (number of fragments where species was detected at least once/total number of fragments)

2 Estimated proportion of patches occupied corrected for detectability (sum of individual occupancy probabilities from each patch/total number of patches, based on parameter estimates from best-fit models)

Species displayed large differences in the proportion of patches occupied (Table 2). Great curassow, crested guan, ocellated turkey and great tinamou had low overall patch occupancy rates. By contrast, the slaty-breasted tinamou, little tinamou and plain chachalaca had relatively high rates of patch occupancy (Table 2).

Species responded strongly to a variety of covariates. Best-fit models had very high AIC weights (>0.5) for five of the six species, indicating low model selection uncertainty for the set of models tested (Table 3). Plain chachalaca was the only species with a high degree of model selection uncertainty and also the only species where the null model was within Δ4 AICc values of the best-fit model. Based on parameter estimates from best-fit model sets greater amounts of habitat in the landscape surrounding patches increased occupancy probabilities for two species (great curassow and little tinamou; Table 4). A greater density of forest patches in the landscape surrounding patches decreased occupancy probabilities of great tinamou and crested guan. Although density of forest patches increased occupancy probabilities of plain chachalaca, unconditional standard errors were larger than model-averaged parameter estimates, suggesting that this variable had a weak effect on occupancy. Similarly, although proportion of regenerating forest in the landscape appeared in the best-fit model for slaty-breasted tinamous, the standard error was high relative to the parameter estimate. Larger focal patch size was associated with increasing rates of occupancy for three species (great tinamou, crested guan and slaty-breasted tinamou; Table 4). Increasing distance to the nearest human community decreased occupancy of little tinamous and increased occupancy of great curassows. Increasing distance to source increased occupancy of both little tinamous and great curassows. We did not find evidence of significant autocorrelation in the residuals of best-fit models for any of the species.

Table 3 Patch occupancy models for galliforms and tinamous. All models that fall within Δ4 AICc units of the best-fit model are shown (except for chachalacas where only models within Δ2 AICc values of best-fit models are shown because of the large number of models within Δ4 AICc values).

1 P, Probability of detection; ψ, probability of occupancy. Explanatory variables are: effort, sampling effort in patch; FPS, size of the sampled patch; PF, proportion of forest in the landscape (i.e. habitat amount); DFP, density of forest patches in the landscape (i.e. habitat configuration); PREG, proportion of regenerating forest in the landscape; DISTC, distance to nearest human community; DISTS, distance to boundary of the Maya Biosphere Reserve

2 Number of parameters in the model

3 Maximized log-likelihood

4 Model weights

Table 4 Parameter estimates (±SE) for all explanatory variables that appeared in best-fit model sets. For variables that appeared in more than one model within the best-fit model set for a single species, the values shown are model-averaged parameter estimates and unconditional standard errors.

Discussion

The four largest species in our study (great curassow, crested guan, ocellated turkey and great tinamou) were the least common in forest patches and therefore the most vulnerable to forest loss and fragmentation. By contrast the three smallest species (little tinamou, slaty-breasted tinamou and plain chachalaca) were the least vulnerable and occupied over half of the surveyed patches. Although patterns of relative use of patches may be misleading, as certain species may naturally be more rare and therefore occupy fewer patches (i.e. the problem of passive sampling effects), a subsidiary analysis correcting for these effects by comparing occupancy probabilities of forest patches with continuous forest sites confirmed our results (D. Thornton, unpubl. data). Patch occupancy patterns of almost all species in our study were influenced strongly by measures of landscape context, with relatively little influence of habitat structure or disturbance within patches. Focal patch size also exerted the expected strong positive influence on occupancy patterns.

Patterns in relative vulnerability of the studied species to habitat loss and fragmentation largely agree with interview-based occupancy surveys in the nearby Yucatán (Urquiza-Haas et al., Reference Urquiza-Haas, Peres and Dolman2009). Larger species in the Yucatán study (great curassows, crested guan, great tinamous and ocellated turkeys) had lower occupancy rates than the smaller tinamous or plain chachalaca. Although little additional work detailing patch occupancy rates of our study species or their close relatives has been published, large galliforms and tinamous are often rare or absent in isolated forest remnants (e.g. Robinson, Reference Robinson1999; Pearson et al., Reference Pearson, Anderson, Mitchell, Rosenberg, Navarrete and Coopmans2009; Manica et al., Reference Manica, Telles and Dias2010). In addition to large body size, which is often implicated as being closely associated with vulnerability to fragmentation in tropical birds (Kattan et al., Reference Kattan, Alvarezlopez and Giraldo1994; Lees & Peres, Reference Lees and Peres2010), these species share other traits that are often associated with vulnerability in neotropical avian species, such as a frugivore–granivore diet and forest dependency (Vetter et al., Reference Vetter, Hansbauer, Vegvari and Storch2011). These species are also often heavily hunted, which can increase vulnerability in fragmented landscapes (Peres, Reference Peres2001).

The strong negative effect of habitat configuration (density of forest patches) on crested guans and great tinamous suggests that more fragmented landscapes are less suitable for these species. The response of these two species contrasts with many studies of temperate bird fauna, which have almost universally shown that habitat configuration has relatively weak effects on species distribution and abundance patterns compared with the influence of habitat area (Fahrig, Reference Fahrig2003; Smith et al., Reference Smith, Koper, Francis and Fahrig2009), but is consistent with other recent work on tropical birds (Cerezo et al., Reference Cerezo, Perelman and Robbins2010). A negative response to increasingly fragmented habitat, holding habitat amount constant, could be driven by sensitivity to edges as more fragmented habitats tend to have a higher density of edges. Tropical birds respond strongly to edges (Laurance et al., Reference Laurance, Lovejoy, Vasconcelos, Bruna, Didham and Stouffer2002; Lindell et al., Reference Lindell, Riffell, Kaiser, Battin, Smith and Sisk2007) and may be more likely to respond in a negative manner to edges than temperate species because of factors such as a shorter history of exposure to patchiness, narrower resource requirements, or other factors (Lindell et al., Reference Lindell, Riffell, Kaiser, Battin, Smith and Sisk2007). The relative influence of habitat amount vs configuration on tropical bird species should be further examined on multiple large and independent landscapes as our study evaluated these factors only in the neighbourhood around patches (a patch–landscape study sensu McGarigal & Cushman, Reference McGarigal and Cushman2002).

The positive effect of the amount of habitat in the landscape and also of focal patch size on occupancy patterns is consistent with many other studies of bird response to fragmentation (e.g. Winter et al., Reference Winter, Johnson, Shaffer, Donovan and Svedarsky2006; Evans et al., Reference Evans, Newson and Gaston2009) and may be driven by positive responses to increasing habitat area or avoidance of edges (Reino et al., Reference Reino, Beja, Osborne, Morgado, Fabiao and Rotenberry2009). The differential effect of distance to human community on species in our study adds to recent work showing that human disturbance can exert a strong impact on occupancy patterns of avian species in fragmented neotropical habitats, independent of landscape structure (Urquiza-Haas et al., Reference Urquiza-Haas, Peres and Dolman2009; Aubad et al., Reference Aubad, Aragon and Rodriguez2010). Patches closer to human habitations are often used more frequently for hunting (Parry et al., Reference Parry, Barlow and Peres2009) and experience higher levels of other forms of human disturbance (e.g. collection of firewood; Aubad et al., Reference Aubad, Aragon and Rodriguez2010).

The general lack of response by species in our study to within-patch vegetation structure or disturbance is counter to many studies of tropical bird species that find an effect of patch quality on distribution patterns among patches (e.g. Lees & Peres, Reference Lees and Peres2010). This finding may indicate a greater sensitivity of these large-bodied species to coarse-grained landscape factors rather than fine-grained within-patch factors, or that the within-patch factors we measured do not adequately represent patch quality for our species. Our exclusion of heavily disturbed patches from analysis, as well as our binary (instead of continuous) measures of cattle and fire disturbance, may have also contributed to a lack of effect of within-patch variables on occupancy patterns.

Conclusions and management recommendations

In our study large terrestrial species were highly vulnerable to habitat loss and fragmentation in northern Guatemala. Because we were able to incorporate detectability in our analysis, our findings do not reflect detection differences among species but instead depict real patterns in vulnerability. Although none of the species that showed a marked negative response are currently considered threatened, all three species are declining across their range, with the great curassow categorized as Vulnerable and the ocellated turkey as Near Threatened on the IUCN Red List (IUCN, 2011). Significantly, deforestation in our study area began only 30–40 years ago yet this was ample time to see substantial declines in occupancy of the study species. Synergistic effects of reduced habitat area and increased hunting pressure could quickly drive the decline of these and similar large-bodied species in fragmented landscapes (Peres, Reference Peres2001). Research to separate effects of hunting pressure vs landscape conversion on large bird species may be useful to identify the most effective management options for such species in human-dominated landscapes.

The strong response of our study species to landscape context emphasizes the importance of management efforts that go beyond consideration of local factors such as patch size or structure (Villard et al., Reference Villard, Trzcinski and Merriam1999). Consideration of the surrounding landscape may be particularly important for highly vagile/large-bodied species such as those in our study. The strong response to habitat configuration suggests that, in addition to preserving sufficient amounts of habitat, increased attention to management of landscape pattern may be necessary to preserve tropical bird species. This could take the form of conservation planning in conjunction with outreach efforts to landowners (e.g. Thomson et al., Reference Thomson, Moilanen, Vesk, Bennett and MacNally2009) to maintain forest cohesion across the landscape.

To date camera trapping has been an underutilized methodology for surveying bird species (O'Brien & Kinnaird, Reference O'Brien and Kinnaird2008). Our findings demonstrate that this methodology provides a complement to audiovisual surveys and results in relatively high detection rates for large terrestrial species. Our success at photo-capturing terrestrial birds may be related to the fact that we placed cameras at relatively low heights (10–20 cm) off the ground so that smaller animals would not escape detection. Significantly, this methodology did not affect our ability to detect larger mammals, although it would have affected our ability to individually identify larger mammals had this been a goal of our study. We support the recommendation of O'Brien & Kinnaird (Reference O'Brien and Kinnaird2008) that camera trapping provides an effective means to estimate occupancy parameters for elusive terrestrial bird species.

Acknowledgements

We thank the Wildlife Conservation Society–Guatemala, the Consejo Nacional para Áreas Protegidas, and the administration of Tikal National Park for supporting this work, and Nery Jurado and Demetrio Córdova for assistance with data collection in the field. Funding for this work was provided by a National Science Foundation Doctoral Dissertation Enhancement Grant, Wildlife Conservation Society–Guatemala, American Society of Mammalogists, University of Florida–Integrative Graduate Education and Research Traineeship Working Forest in the Tropics, and Idea-Wild.

Appendix The 29 occupancy models tested for each species. Note that these models were fit after a best-fit detection model was identified for each species, based on the sampling effort of covariates and season of survey.

* FPS, size of the sampled patch; PF, proportion of forest in the landscape (i.e. habitat amount); DFP, density of forest patches in the landscape (i.e. habitat configuration); PREG, proportion of regenerating forest in the landscape; BALRG, basal area of large trees; BASM, basal area of small trees/saplings; CAT, presence of cattle within patch (binary); FIRE, indications of past fire within patch (binary); DISTC, distance to nearest human community; DISTS, distance to boundary of the Maya Biosphere Reserve

Biographical sketches

Daniel Thornton is interested in understanding the multiscaled response of species to landscape change and the application of distribution and occupancy modelling to the conservation of large vertebrates. Lyn Branch has broad interests in the ecology and conservation of wildlife in Latin America, with a focus on impacts of landscape change on wildlife populations and synergisms between these processes and other anthropogenic factors that affect wildlife. Mel Sunquist has an interest in mammalian ecology and conservation, with a focus on how best to conserve carnivores in increasingly human-dominated landscapes.

Footnotes

*

Current address: Department of Biology, Trent University, 1600 West Bank Drive, Peterborough, Ontario K9J7B8, Canada

* FPS, size of the sampled patch; PF, proportion of forest in the landscape (i.e. habitat amount); DFP, density of forest patches in the landscape (i.e. habitat configuration); PREG, proportion of regenerating forest in the landscape; BALRG, basal area of large trees; BASM, basal area of small trees/saplings; CAT, presence of cattle within patch (binary); FIRE, indications of past fire within patch (binary); DISTC, distance to nearest human community; DISTS, distance to boundary of the Maya Biosphere Reserve

References

Asner, G.P., Rudel, T.K., Aide, T.M., DeFries, R. & Emerson, R. (2009) A contemporary assessment of change in humid tropical forests. Conservation Biology, 23, 13861395.CrossRefGoogle ScholarPubMed
Aubad, J., Aragon, P. & Rodriguez, M.A. (2010) Human access and landscape structure effects on Andean forest bird richness. Acta Oecologia, 36, 396402.CrossRefGoogle Scholar
Bertsch, C. & Barreto, G.R. (2008a) Abundance and home range of the yellow-knobbed curassow (Crax daubentoni) in the Venezuelan central llanos. Ornitologia Neotropical, 19, 287293.Google Scholar
Bertsch, C. & Barreto, G.R. (2008b) Diet of the yellow-knobbed curassow in the central Venezuelan llanos. The Wilson Journal of Ornithology, 120, 767777.CrossRefGoogle Scholar
Brooks, D.M. & Strahl, S.D. (2000) Curassows, Guans, and Chachalacas: Status Survey and Conservation Action Plan for Cracids 2000–2004. IUCN/Species Survival Commission Cracid Specialist Group, Gland, Switzerland.Google Scholar
Burnham, K. & Anderson, D. (2002) Model Selection and Multimodel Inference. Springer-Verlag, New York, USA.Google Scholar
Cerezo, A., Perelman, S. & Robbins, C.S. (2010) Landscape-level impact of tropical forest loss and fragmentation on bird occurrence in eastern Guatemala. Ecological Modelling, 221, 512526.CrossRefGoogle Scholar
Cottom, G., Curtis, J.T. & Hale, B.W. (1953) Some sampling characteristics of a population of randomly dispersed individuals. Ecology, 34, 741757.CrossRefGoogle Scholar
Dillon, A. & Kelly, M.J. (2008) Ocelot home range, overlap and density: comparing radio telemetry with camera trapping. Journal of Zoology, 275, 391398.CrossRefGoogle Scholar
Érard, C., Théry, M. & Sabatier, D. (2007) Fruit characters in the diet of syntopic large frugivorous forest bird species in French Guiana. Revue D'Ecologie–La Terre Et La Vie, 62, 323350.Google Scholar
Evans, K.L., Newson, S.E. & Gaston, K.J. (2009) Habitat influences on urban avian assemblages. Ibis, 151, 1939.CrossRefGoogle Scholar
Fahrig, L. (2003) Effects of habitat fragmentation on biodiversity. Annual Review of Ecology, Evolution, and Systematics, 34, 487515.CrossRefGoogle Scholar
Fortin, M.-J. & Dale, M. (2005) Spatial Analysis: A Guide for Ecologists. Cambridge University Press, Cambridge, UK.CrossRefGoogle Scholar
Gascon, C., Lovejoy, T.E., Bierregaard, R.O. Jr, Malcolm, J.R., Stouffer, P.C., Vasconcelos, H.L. et al. (1999) Matrix habitat and species richness in tropical forest remnants. Biological Conservation, 91, 223229.CrossRefGoogle Scholar
Gonzalez, M.J., Quigley, H. & Taylor, A.C. (1998) Habitat use and reproductive ecology of the ocellated turkey in Tikal National Park, Guatemala. The Wilson Bulletin, 110, 505510.Google Scholar
Hines, J.E. (2006) presence 2 – Software to Estimate Patch Occupancy and Related Parameters. US Geological Service–Patuxent Wildlife Research Center. Http://www.mbr-pwrc.usgs.gov/software/presence.html [accessed 29 March 2010].Google Scholar
IUCN (2011) IUCN Red List of Threatened Species v. 2011.2. IUCN, Gland, Switzerland. Http://www.iucnredlist.org [accessed 12 March 2012].Google Scholar
Kattan, G.H., Alvarezlopez, H. & Giraldo, M. (1994) Forest fragmentation and bird extinctions: San Antonio eighty years later. Conservation Biology, 8, 138146.CrossRefGoogle Scholar
Laurance, W.F. (1991) Ecological correlates of extinction proneness in Australian tropical rain forest mammals. Conservation Biology, 5, 7989.CrossRefGoogle Scholar
Laurance, W.F. (2002) Hyperdynamism in fragmented habitats. Journal of Vegetation Science, 13, 595602.CrossRefGoogle Scholar
Laurance, W.F., Lovejoy, T.E., Vasconcelos, H.L., Bruna, E.M., Didham, R.K., Stouffer, P.C. et al. (2002) Ecosystem decay of Amazonian forest fragments: a 22-year investigation. Conservation Biology, 16, 605618.CrossRefGoogle Scholar
Lees, A.C. & Peres, C.A. (2010) Habitat and life history determinants of antbird occurrence in variable-sized Amazonian forest fragments. Biotropica, 42, 614621.CrossRefGoogle Scholar
Lichstein, J.W., Simons, T.R., Shriner, S.A. & Franzreb, K.E. (2002) Spatial autocorrelation and autoregressive models in ecology. Ecological Monographs, 72, 445463.CrossRefGoogle Scholar
Lindell, C.A., Riffell, S.K., Kaiser, S.A., Battin, A.L., Smith, M.L. & Sisk, T.D. (2007) Edge responses of tropical and temperate birds. The Wilson Journal of Ornithology, 119, 205220.CrossRefGoogle Scholar
MacKenzie, D.I. & Bailey, L.L. (2004) Assessing the fit of site-occupancy models. Journal of Agricultural, Biological, and Environmental Statistics, 9, 300318.CrossRefGoogle Scholar
MacKenzie, D.I., Nichols, J.D., Royle, A.J., Pollock, K.H., Bailey, L.L. & Hines, J.E. (2006) Occupancy Estimation and Modeling: Inferring Patterns and Dynamics of Species Occurrence. Elsevier, Boston, USA.Google Scholar
Maffei, L., Cuellar, E. & Noss, A. (2004) One thousand jaguars (Panthera onca) in Bolivia's Chaco? Camera trapping in the Kaa-Iya National Park. Journal of Zoology, 262, 295304.CrossRefGoogle Scholar
Manica, L.T., Telles, M. & Dias, M. (2010) Bird richness and composition in a cerrado fragment in the state of São Paulo. Brazilian Journal of Biology, 70, 243254.CrossRefGoogle Scholar
McGarigal, K. & Cushman, S.A. (2002) Comparative evaluation of experimental approaches to the study of habitat fragmentation effects. Ecological Applications, 12, 335345.CrossRefGoogle Scholar
McGarigal, K., Cushman, S.A., Neel, M.C. & Ene, E. (2002) FRAGSTATS: Spatial Pattern Analysis Program for Categorical Maps. Http://www.umass.edu/landeco/research/fragstats/fragstats.html [accessed 13 March 2012].Google Scholar
Nichols, J.D., Bailey, L.L., O'Connell, A.F., Talancy, N.W., Grant, E.H.C., Gilbert, A.T. et al. (2008) Multi-scale occupancy estimation and modelling using multiple detection methods. Journal of Applied Ecology, 45, 13211329.CrossRefGoogle Scholar
O'Brien, T.G. & Kinnaird, M.F. (2008) A picture is worth a thousand words: the application of camera trapping to the study of birds. Bird Conservation International, 18, S144S162.CrossRefGoogle Scholar
Parry, L., Barlow, J. & Peres, C.A. (2009) Allocation of hunting effort by Amazonian smallholders: implications for conserving wildlife in mixed-use landscapes. Biological Conservation, 142, 17771786.CrossRefGoogle Scholar
Pearson, D.L., Anderson, C.D., Mitchell, B.R., Rosenberg, M.S., Navarrete, R. & Coopmans, P. (2009) Testing hypotheses of bird extinctions at Rio Palenque, Ecuador, with informal species lists. Conservation Biology, 24, 500510.CrossRefGoogle ScholarPubMed
Peres, C.A. (2001) Synergistic effects of subsistence hunting and habitat fragmentation on Amazonian forest vertebrates. Conservation Biology, 15, 14901505.CrossRefGoogle Scholar
Poulsen, B.O. (1994) Movements of single birds and mixed-species flocks between isolated fragments of cloud forest in Ecuador. Studies on Neotropical Fauna and Environment, 29, 149160.CrossRefGoogle Scholar
Reino, L., Beja, P., Osborne, P.E., Morgado, R., Fabiao, A. & Rotenberry, J.T. (2009) Distance to edges, edge contrast and landscape fragmentation: interactions affecting farmland birds around forest plantations. Biological Conservation, 142, 824838.CrossRefGoogle Scholar
Robinson, J.G., Brawn, J.D. & Robinson, S.K. (2000) Forest bird community structure in central Panama: influence of spatial scale and biogeography. Ecological Monographs, 70, 209235.CrossRefGoogle Scholar
Robinson, W.D. (1999) Long-term changes in the avifauna of Barro Colorado Island, Panama, a tropical forest isolate. Conservation Biology, 13, 8597.CrossRefGoogle Scholar
Sigel, B.J., Robinson, W.D. & Sherry, T.W. (2010) Comparing bird community responses to forest fragmentation in two lowland Central American reserves. Biological Conservation, 143, 340350.CrossRefGoogle Scholar
Smith, A.C., Koper, N., Francis, C.M. & Fahrig, L. (2009) Confronting collinearity: comparing methods for disentangling the effects of habitat loss and fragmentation. Landscape Ecology, 24, 12711285.CrossRefGoogle Scholar
Terborgh, J., Robinson, S.K., Parker, T.A. III, Munn, C.A. & Pierpont, N. (1990) Structure and organization of an Amazonian forest bird community. Ecological Monographs, 60, 213238.CrossRefGoogle Scholar
Thiollay, J.-M. (1999) Responses of an avian community to rain forest degradation. Biodiversity and Conservation, 8, 513534.CrossRefGoogle Scholar
Thiollay, J.-M. (2005) Effects of hunting on Guianan forest game birds. Biodiversity and Conservation, 14, 11211135.CrossRefGoogle Scholar
Thomson, J.R., Moilanen, A.J., Vesk, P.A., Bennett, A.F. & MacNally, R. (2009) Where and when to revegetate: a quantitative method for scheduling landscape reconstruction. Ecological Applications, 19, 817828.CrossRefGoogle ScholarPubMed
Thornton, D.H., Branch, L.C. & Sunquist, M.E. (2011a) Passive sampling effects and landscape location alter associations between species traits and response to fragmentation. Ecological Applications, 21, 817829.CrossRefGoogle ScholarPubMed
Thornton, D.H., Branch, L.C. & Sunquist, M.E. (2011b) The relative influence of habitat loss and fragmentation: do tropical mammals meet the temperate paradigm? Ecological Applications, 21, 23242333.CrossRefGoogle ScholarPubMed
Urquiza-Haas, T., Peres, C.A. & Dolman, P.M. (2009) Regional scale effects of human density and forest disturbance on large-bodied vertebrates throughout the Yucatán Peninsula, Mexico. Biological Conservation, 142, 134148.CrossRefGoogle Scholar
Vetter, D., Hansbauer, M., Vegvari, Z. & Storch, I. (2011) Predictors of forest fragmentation sensitivity in Neotropical vertebrates: a quantitative review. Ecography, 34, 18.CrossRefGoogle Scholar
Villard, M.A., Trzcinski, M.K. & Merriam, G. (1999) Fragmentation effects on forest birds: relative influence of woodland cover and configuration on landscape occupancy. Conservation Biology, 13, 774783.CrossRefGoogle Scholar
Winter, M., Johnson, D.H., Shaffer, J.A., Donovan, T.M. & Svedarsky, W.D. (2006) Patch size and landscape effects on density and nesting success of grassland birds. Journal of Wildlife Management, 70, 158172.CrossRefGoogle Scholar
Figure 0

Fig. 1 Location of the 50 sampled patches within the study area (Thornton et al., 2011b). Inset shows approximate location of the main map in northern Guatemala.

Figure 1

Table 1 Patch and landscape variables for the 50 focal patches in the study area in northern Guatemala (Fig. 1).

Figure 2

Table 2 Estimates of detectability and proportion of patches occupied for the seven large terrestrial birds detected in our study.

Figure 3

Table 3 Patch occupancy models for galliforms and tinamous. All models that fall within Δ4 AICc units of the best-fit model are shown (except for chachalacas where only models within Δ2 AICc values of best-fit models are shown because of the large number of models within Δ4 AICc values).

Figure 4

Table 4 Parameter estimates (±SE) for all explanatory variables that appeared in best-fit model sets. For variables that appeared in more than one model within the best-fit model set for a single species, the values shown are model-averaged parameter estimates and unconditional standard errors.

Figure 5

Appendix The 29 occupancy models tested for each species. Note that these models were fit after a best-fit detection model was identified for each species, based on the sampling effort of covariates and season of survey.