Abstract
Wildlife species population density estimation is important from both ecological and management perspectives. The pellet-group counting method has been used to evaluate density of the white-tailed deer (
Introduction
Efficient and accurate population density estimates of wildlife species are an important aspect of management planning [1]. In areas where direct observation of animals is difficult, due to either low local abundance of the species, secretive behavior, low visibility, or human disturbance, the use of indirect methods such as the quantification of tracks and feces (pellet groups) can be more appropriate [2,3]. Pellet-group counting methods have been assessed to accurately evaluate the population density of both old and new world ungulates, lagomorphs, and elephants [e.g., 4–9]. In the field, pellet groups have been quantified in sampling units (e.g., circular plots, strip transects, ST) and for estimating decay rate (FSC, Fecal Standing Crop), or in re-sampling the same plots two or more times (FAR, Fecal Accumulation Rate) [10]. However, these methods do not take feces located outside the sampling units into account, and this can reduce the sampling size and affect the accuracy of density estimation, particularly in those sites with low abundance of animals. As an alternative, the line transect (LT) method has been applied to count pellet groups, and is a technique that has been shown to achieve an increase in precision of density estimates [11,12].
Besides its ecological roles as herbivore, seed disperser, and prey [13], the white-tailed deer (
Methods
The study was carried out in central Mexico, in the Mixteca region located in the southern part of Puebla state (18° 12′ N and 98° 46′ W; Fig. 1). The topography of the region is characterized by numerous hills and narrow valleys of elevations ranging from 600 to 2750 m. The mean annual precipitation and temperature are 817 mm and 21°C, respectively. The climate is warm sub-humid with rains in summer, and the predominant vegetation is seasonal-dry tropical forest. Estimates of white-tailed deer density (1 to 7 ind/km2) had been previously obtained in several locations in this region [17,18]. According to the findings of these studies, we selected a specific location (Jolalpan, Fig. 1) with the highest deer density in order to intensively sample pellet groups within a 1 km2 experimental plot. Habitat characteristics within this plot are relatively homogenous: the dominant vegetation is tropical dry forest with a well-developed understory that provides food and shelter to the deer, while the terrain is very irregular with pronounced slopes. As is normally the case throughout this region, the principal human economic activity is livestock production. Conversely, hunting activity is low in this location.

Geographic location of the study site (El Salado, Jolalpan) in the southern part of Puebla State, Mexico.
In May of 2006, we selected a random point location at which to establish a set of five parallel 1 km long transects, evenly spaced 250 m apart. These were used to apply the four quantification methods: FSC, FAR, ST, and LT. In the case of FSC and FAR, we marked 100 circular plots each of 9.3 m2 along each transect, spaced every 10 m (center to center) and counted all pellet groups within each plot. For ST, we defined a width of 1 m along the length of transect, considering this width to be adequate to census complete pellet groups in the dense vegetation present. For LT, we measured the perpendicular distance from the transect line to every fecal group observed along the length of transect. Particular attention was paid to meeting the assumptions of this method [19]. In May 2007, we repeated the process in the same transects and plots as recorded the previous year. Following this methodology, we obtained two snapshots for the FSC, ST, and LT methods. Since the FAR method employs an estimate of pellet group accumulation between two points in time, we identified and marked the pellet groups in the circular plots in January 2007 in order to quantify the new groups found after 60 and 120 days. The surveys were carried out by three to five people; to reduce observer bias, the same observer (ACS) counted the pellet groups for each treatment in each year. During this survey, we recorded the time necessary to count pellet groups in both the circular and transect methods.
To estimate population density (D, individuals km−2), we applied the equation proposed by Eberhardt and Van Etten [20] as: D = (NP × Dpg)/(T × dR), where NP = number of plots (circular or strip transects) per square kilometer, Dpg = mean pellet groups; T = either deposition time of fecal groups in the case of FAR, or mean decay time in FSC; and dR = defecation rate. In the case of FSC and FAR, Dpg was estimated considering each circular plot as a sampled unit. To estimate Dpg with the ST method, we used the equation: Dpg = n/2Lw, where n = number of fecal groups, w = transect width (1 m) and L = total transect length sampled. In this case, the standard error (Se) was estimated using each individual transect as a replicate. For the LT method, we used the equation: Dpg = nf(0)/2L, where f(0) = probability density function at 0 m. We used the DISTANCE 5.0 software to estimate Dpg, f(0) and Se [21]. We used likelihood ratio tests and Akaike's Information Criterion test as objective and quantitative methods for model selection (Uniform + cos) using non-truncated and truncated (5% and 10%) data.
For each method, we estimated the annual variance of density (D) according to Plumptre [22] as: Var(D) = [(NP × Dpg)/(T × dR)]2 × [(CV(Dpg))2 + (CV(dR))2 + (CV(T))2]. The estimates of CV(Dpg) for FSC and FAR were calculated following a negative binomial distribution to estimate Se with the equation: Se = √[(× + x2)/(k/n)], where k = parameter of the binomial negative, and n
To assess the effect of sampling effort (number and length of transects) on density precision, we re-sampled the database obtained in May 2006 with a code written in R [24]. We simulated (N = 100 replicates) sampling from one to 10 transects of 390 m and from one to five transects of 790 m in length. The simulation transects were placed randomly along the 5 km of field sampled transects. To determine the transect length (L) necessary in order to detect changes in population density using the LT method in particular, we followed Buckland et al. [19] as: L = L0 × CV(D) 2 / CVt(D)2, where L0 and CV(D) are the transect length and coefficient of variation of density of the pilot study, respectively, while CVt(D) is the selected coefficient of variation of population density (10%, 20%, and 30%). We estimated the required L, considering data from 2006 and 2007 as a pilot study.
Results
In May 2006, we counted 71, 135, and 326 fecal groups with the FSC, ST, and LT methods, respectively, while in May 2007, 48, 101, and 293 fecal groups were recorded. In May 2007, 21 and 48 fecal groups were quantified with the FAR method at 60 and 120 days of accumulation, respectively. Pellet decomposition rate was estimated as 123 ± 2.4 days.
As expected, the number of pellet-groups counted decay with the increase in perpendicular distance from the center of the transect (Fig. 2). According to LT, detectability probability must to be 1.0 at 0 m, but our observations indicate that complete observation of pellet groups was obtained between 0 to 1.0 m of perpendicular distance.

Histograms of the frequency of counted pellets-groups at different perpendicular distances in the LT method 2006 (a) and 2007 (b). 100% of the pellet-groups between 0 and 1 m were counted in the ST method.
According to data from January 2007, the mean sampling time per 500 m long transect to count pellet-groups in circular plots was 1.48 h (range 1.08 to 2.34 h), while in transects (LT) this was significantly different (Wilcoxon test, W = 48.5, P = 0.03) at 2.14 h (range 0.42 to 3.0 h). Thus, sampling in the transects was 1.2 times more time-consuming than in circular plots; however, the number of pellet groups counted was 6.6 times greater in the transects than in the circular plots (Fig. 3).
In 2006, the population density estimates were statistically similar between the FSC, ST, and LT methods (Appendix 1, F = 0.146, df = 2, 12, P = 0.87); and also in 2007, when the FAR estimation was included (F = 21.35, df = 4, 20, P = 0.11). Mean density estimates were similar between years in each method (paired

Relationships between sampling time and number of pellet-groups counted in transects (•) and circular plots (○) during January 2007. Sampling time is given considering a 500 m transect.

Relationship between simple size (number of pellet groups counted, n) and coefficient of variation (%) of the population density (CVD). Sampling pellet groups in FAR (▪), FSC (▫), ST (•) and LT (○) methods. The line represents the best-fitted linear model (coefficient of determination, R2).
The precision of the density estimates was dependent of the number of pellet groups counted (Fig. 4; r2 = 0.83; F = 28. 4; df = 6; P = 0.002). The methods were ranked in descending order of estimate precision as: LT, ST, FSC, and FAR.
The results of the simulations indicate that mean density did not vary as the number and length of transects increased in the FSC, ST, and LT methods (Fig. 5). In contrast, the precision of the estimates was found to improve with increased sampling effort. In the case of the LT method, estimates of transect length necessary to achieve 10%, 20%, and 30% of precision (CV%) were: 8.5, 2.1, and 0.9 km respectively, based on the data from 2006, while for 2007 these were 33.8, 8.5, and 3.8 km.
Discussion
Although fecal pellet counts have been widely used to index changes in deer abundance in forests, few studies have modeled the relationship between such indices and deer density [25]. For example, Forsyth

Simulation of the effects of sampling effort (transect length and number) on population density and precision (95% confidence interval) estimates for FSC (a), ST (b) and LT (c) methods over transects of 390 m (—♦—) and 790 m (—□—) in length.
To convert counts of pellet groups to population density, a defecation rate of 12.7 was initially used [20]. However, some other estimates have been: 13.7 [28], 19.6 (CV = 11.7) [29], 26.9 [30], and 34.0 [31] in the United States; 15.2 [32], 17.0 (CV = 9.6) [23], and 20.9 [33] in Mexico; and 14.0 in Venezuela [34]. The defecation rate varied depending on factors such as age, sex, and food quality. In particular, some data suggest that defecation rates are higher in wild deer in comparison with domesticated individuals [29,31]. In consequence, the defecation rate is a crucial parameter in the conversion of pellet groups to population density. In our study, we used the observed defecation rate in tame deer of the same subspecies close to the studied region; but, bias in population density exists in our estimates. We tried to reduce this bias by incorporating the variation of defecation rate. However, it is necessary to estimate the actual defecation rate in the studied site, which is data not easily obtained in wild animals.
Our results indicate that the four methods give similar mean population density estimates. However, in terms of precision, the results suggest that the LT is the most suitable method for fecal group quantification in terms of obtaining more precise density estimations of white-tailed deer in this tropical habitat. Marques
In contrast, the robustness of LT assumptions provides a better method for counting pellet groups in this type of habitat. However, very long transect lengths are necessary in order to reach a very precise density estimate (CV < 10%), a factor which could consume more sampling time in comparison to the ST and FSC methods. Although transect length is largely dependent on topographic, vegetation cover, and other conditions, Marques
Implications for conservation
The white-tailed deer is not considered an endangered species by the IUCN [38]; however, local extinctions as a result of over-hunting and habitat loss are common in many sites in the Neotropics [39]. Of the 38 recognized subspecies, at least 20 have a Neotropical distribution [39], but few (<6) have been studied [15,40], and population density data is not available for many sites [16]. Therefore, suggestions on methods and field procedures to estimate population densities of the white-tailed deer are a crucial aspect for the conservation of this species. This is particularly urgent in sites where subsistence, commercial, and game hunting of this deer are common. For example, in Mexico the white-tailed deer is one of the most important hunted species in Management Units for the Conservation of Wildlife (in Spanish: “Unidades de Manejo para la Conservación de la Vida Silvestre” or UMA). Important concerns and limitations of this model have been discussed [41], but one common aspect emphasized by these studies is the need for accurate density estimations. Another example, also in Mexico, is that even Protected Natural Areas are important conservation assets for this species and provide an opportunity to study deer in less perturbed sites [42], although few are of sufficient area to support a hypothetical minimum viable population [43]. In this case, accurate population density estimates are required for management purposes. A similar situation occurs in other Neotropical regions such as Costa Rica, Colombia, and Venezuela [16].

Two important variables to convert an index of abundance into population density are the rate of defecation and the rate of decomposition of fecal groups. An example of old (a) and fresh (b) fecal groups of white-tailed deer in captivity (c, d).
We strongly suggest replicating our study in other tropical areas in order to test these sampling methods. An important requirement is the comparison of these methods under different levels of white-tailed deer abundance. For example, when we revised methods of density estimation, quantification of fecal groups in circular plots was used successfully in tropical dry forests with higher population densities of > 11 white-tailed deer km−2 [37]; while strip transects were used in sites with medium densities of 6–9 individuals km−2 [44]. López-Tellez
Footnotes
Acknowledgments
This study was supported by the CONACYT and the Red de Biología y Conservación de Vertebrados, Instituto de Ecologia A. C. We are grateful for comments and suggestions by R. Guevara, E. Naranjo, G. Sánchez-Rojas and J. Pérez-Torres to an early draft of this work. Thanks also go to the people of Rancho El Salado, and to C. López-Tellez, L. A. Escobedo-Morales, A. González and J. López for their help during the field work. T. Perez helped us with logistic work. K. MacMillan revised the English transcription.
Appendix 1.
Attributes of woody vegetation structure and composition for plots in
Appendix 2.
Advantages (+) and disadvantages (−) of pellet-group counting methods evaluated in this study. Abbreviation: Fecal Accumulation Rate (FAR), Fecal Standing Crop (FSC), Strip Transect (ST) and Line Transect (LT).
