Road mortality is probably the best-known and visible impact of roads upon wildlife. Although several factors influence road-kill counts, carcass persistence time is considered the most important determinant underlying underestimates of road mortality. The present study aims to describe and model carcass persistence variability on the road for different taxonomic groups under different environmental conditions throughout the year; and also to assess the effect of sampling frequency on the relative variation in road-kill estimates registered within a survey.
Daily surveys of road-killed vertebrates were conducted over one year along four road sections with different traffic volumes. Survival analysis was then used to i) describe carcass persistence timings for overall and for specific animal groups; ii) assess optimal sampling designs according to research objectives; and iii) model the influence of road, animal and weather factors on carcass persistence probabilities. Most animal carcasses persisted on the road for the first day only, with some groups disappearing at very high rates. The advisable periodicity of road monitoring that minimizes bias in road mortality estimates is daily monitoring for bats (in the morning) and lizards (in the afternoon), daily monitoring for toads, small birds, small mammals, snakes, salamanders, and lagomorphs; 1 day-interval (alternate days) for large birds, birds of prey, hedgehogs, and freshwater turtles; and 2 day-interval for carnivores. Multiple factors influenced the persistence probabilities of vertebrate carcasses on the road. Overall, the persistence was much lower for small animals, on roads with lower traffic volumes, for carcasses located on road lanes, and during humid conditions and high temperatures during the wet season and dry seasons, respectively.
The guidance given here on monitoring frequencies is particularly relevant to provide conservation and transportation agencies with accurate numbers of road-kills, realistic mitigation measures, and detailed designs for road monitoring programs.
Citation: Santos SM, Carvalho F, Mira A (2011) How Long Do the Dead Survive on the Road? Carcass Persistence Probability and Implications for Road-Kill Monitoring Surveys. PLoS ONE 6(9): e25383. doi:10.1371/journal.pone.0025383
Editor: Brock Fenton, University of Western Ontario, Canada
Received: May 23, 2011; Accepted: September 2, 2011; Published: September 27, 2011
Copyright: © 2011 Santos et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was partially supported by a grant from the Fundação Luis Molina (SMS). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. No additional external funding received for this study.
Competing interests: The authors have declared that no competing interests exist.
Roads can exert severe impacts upon the long-term viability of animal populations , , either through direct killings that decrease the number of individuals (road mortality), or through habitat loss and fragmentation, and barrier effects increasing isolation of populations –. Road mortality is one of the best known and visible impacts of roads on animal populations, with millions of individuals from a wide range of taxonomic groups being killed every year , . The need for effective mitigation measures to minimize impacts of existing and future roads on wildlife populations – has thus lead to an increasing body of research relating the spatial patterns of road-kills with both ecological and road features , –. These studies rely primarily on estimates of road mortality, which are often based on a particular sampling scheme designed for a particular species, thus raising many questions regarding their accuracy and utility for comparative purposes and for guiding monitoring and mitigation plans targeting multiple species.
Several factors have been referred to affect the accuracy of road mortality estimates, including the rate at which the carcasses decompose, the time interval between the occurrence of mortality and road monitoring, the number of vehicles that pass over the carcass, the visibility of carcasses, the abundance and diversity of scavengers, the weather, and the accuracy and precision of the search method –. A common result of studies concerning persistence of carcasses in the field is that carcasses of small animals disappear sooner, usually within the first day of experiments, although some refer to non-road habitats , , . Some animal carcasses are removed by scavengers, such as corvids, domestic cats, polecats and foxes , ,  preventing detection in further surveys. We found only one study that discussed the effect of traffic on persistence. Heavy traffic decreases the access of scavengers to carcasses on the road , thus increasing carcass persistence. On the other hand, weather can influence counting numbers in several ways: poor visibility (e.g., fog) and heavy rain or wet conditions increase persistence rates of carcasses on the road by decreasing scavenger activity , , . Also, a count survey on foot is more effective than on a car vehicle, particularly for smaller animals , , although it is more costly in terms of time and manpower. Among the factors that influence the accuracy of estimates, the time of carcass persistence, defined as the time each animal carcass remained on the road, has been considered the most important factor inducing bias in road mortality estimates .
Although one can easily predict that carcass persistence should be longer for larger animals, and thus that road mortality estimates across different sized species may vary with the frequency of monitoring surveys, few studies have analyzed how carcass persistence probability throughout time may affect road mortality estimates across a large number of taxonomic groups and under different environmental conditions , , , . Moreover, most studies focusing on this subject have been based upon short duration experiences (one-two seasons) using captive-reared bird chicks , . The variability in the rate at which carcasses persist on the road in different situations, throughout the year, must be considered when planning road-kill monitoring surveys so that these data can be used as a reliable indicator of mortality , . In addition, an assessment of the effect of sampling frequency on the number of road-kills registered in a survey can be of high practical value. Establishing the frequency of road-kill surveys that produce higher accurate estimates at low cost is a major goal in most projects. The implementation of high-cost mitigation measures of new road projects often require monitoring schemes on the effectiveness of those measures, while environmental impact assessment studies must include the design of a monitoring program with a detailed description of the methods and periodicity of controls . The accuracy of road surveys, and the factors influencing the number of road carcasses effectively counted, is a critical point of discussion when striving to provide accurate estimates of wildlife collision rates to conservation and transportation agencies.
This study addressed this issue by using daily surveys of road-killed vertebrates and survival analysis to describe carcass persistence along four sections of roads with different characteristics. Specifically, we aimed to answer the following questions:
- How long may an animal carcass persist on a road after being hit by a vehicle?
- How does carcass persistence vary between different taxonomic groups?
- What is the periodicity of road mortality surveys that minimizes losses through undetected carcass removal?
- How does carcass persistence changes with road and animal characteristics? And with weather conditions?
In the present work “carcass disappearance” is defined as a carcass that is no longer available for detection from a moving vehicle, meaning that either it becomes unrecognizable for identification (at least to the studied taxonomic groups level), or is absent from the road. This allowed for the calculation of the number of days each road-killed animal remained on the road (persistence time). Overall, we expect that the analysis of different survival curves should differ greatly among taxonomic groups, allowing their classification based on the median expected persistence time, which in turn should be greatly explained by characteristics of road and animals, and by weather conditions. We further expect that our study may provide particularly useful insights regarding the numerical consequences of choosing different time intervals or accuracy targets for road monitoring studies.
Global persistence time of animal carcasses on the road
Initially, the persistence time of 4447 animal road-kills were analyzed. Small birds and salamanders were the most abundant groups, accounting for 45% and 18.9% of the casualties, respectively. Large birds and freshwater turtles were among the least represented, comprising just 1.1% and 0.5% of the data sample, respectively. The remaining 8 groups accounted each for somewhere between 9.3% (toads) and 2.0% (bats) of the sample (Table 1).
Table 1. Summary of results for persistence estimates for each taxonomic group and the “all taxa” data set (N: sample size; Median (95% CI): median persistence probabilities and corresponding 95% confidence intervals obtained with a Kaplan-Meier estimator; MPT: maximum persistence time recorded (in days); S(t = 1), S(t = 2), S(t = 7): estimate of persistence probability for 1-day, 2-day and 7-day intervals obtained with a Kaplan-Meier estimator).doi:10.1371/journal.pone.0025383.t001
Most animal carcasses persisted on the road for only one day (or less), after being killed by a vehicle. According to predictions, maximum persistence time varied among different taxonomic groups, ranging from 4 days (for lizards) to 158 days (among carnivores; Table 1).
The median of “all taxa” persistence time was one day and corresponded to a 0.5 probability of persistence. This probability decreased to 0.241 on the second day (i.e. 76% of road-killed animals did not persist on the road for longer than two days). These values indicate a low global persistence probability; with a substantial drop beyond the first day (Fig. 1; Table 1).
Figure 1. Overall survival curve for the “all taxa” data set (Kaplan-Meier estimate and corresponding 95% confidence intervals; N = 4447; the length of time axis is limited to a maximum of 15 days).doi:10.1371/journal.pone.0025383.g001
Differences in persistence time between taxonomic groups
Table 1 and Fig. 2 summarize the persistence time and simple survival functions for the taxonomic groups considered. Median persistence time varied between one (for seven of the 13 taxa) and nine days (for carnivores). All seven taxonomic groups with lower median time (toads, salamanders, lizards, snakes, small birds, bats, and small mammals) presented persistence probabilities lower than 0.50 after the first 24 h, meaning that, after that time, fewer than 50% of road-killed animals from these taxa were still on the road. Survival curves provided the classification of each taxonomic group according to their persistence time (median and maximum values). Very short persistence time (<1 day) was characteristic of lizards and bats (the lowest values observed); and short time (1 to 2 days) was registered for toads, salamanders, snakes, small birds, small mammals, and lagomorphs. Freshwater turtles, large birds, birds of prey, and hedgehogs had medium persistence time (3 to 6 days); and carnivores had the longest persistence time on the road (>7 days; see Fig. S1, S2, and S3 in Supporting Information).
Figure 2. Separate survival curves for the 13 taxonomic groups; Toads: toads (including frogs), Lbirds: large birds (more than 200 g, and excluding raptors), Sbirds: small birds (less than 200 g), Carniv: carnivores; Fturtles: freshwater turtles, Hedgeh: hedgehogs, Lizar: lizards, Lagom: lagomorphs, Smammals: small mammals, Bats: bats, Bprey: birds of prey (diurnal and nocturnal), Snake: snakes, Salam: salamanders (including newts); (Kaplan-Meier estimates; N = 4447; the length of time axis is limited to a maximum of 15 days).doi:10.1371/journal.pone.0025383.g002
Influence of monitoring periodicity on road-kill numbers
To discuss the implications of varying time frequencies in monitoring studies, three scenarios of sampling frequency were defined: 7-day, 2-day and 1-day intervals. In addition, two threshold persistence probabilities also were defined: 0.50 and 0.70. These values were chosen to represent two different possible monitoring goals (and distinct financial budgets), and were assumed to represent, respectively, that roughly 50% and 70% of road-kills are still recorded by a monitoring program. Thus, for a scenario of weekly road monitoring (7-day interval), carnivores and birds of prey were the only groups with nearly a 0.50 probability of persisting on the road. All other taxa exhibited lower probabilities, most of them less than 0.10. Within a scenario of road monitoring every 3 days (2 day-intervals between surveys), more than a 0.50 probability of persistence existed for freshwater turtles, large birds, birds of prey, hedgehogs and carnivores. For this monitoring frequency, all other groups had less than a 0.25 persistence probability. Considering the scenario of monitoring on alternate days (1 day-intervals), carnivores, freshwater turtles and hedgehogs had more than a 0.75 persistence probability; birds of prey, large birds and lagomorphs a probability between 0.50 and 0.75; salamanders, snakes, small mammals, small birds and toads a probability between 0.25 and 0.50; and bats and lizards less than a 0.25 probability of persistence (Table 1).
According to the previously defined monitoring goals and budgets, considering 50% to be an acceptable persistence probability (assuming a lower financial budget), the periodicity of road monitoring that minimizes bias on road mortality estimates is: daily monitoring for lizards (in the afternoon period) and bats (in the early morning), toads, small birds, small mammals, snakes, and salamanders; a one-day interval for lagomorphs; a two-day interval for freshwater turtles; a three-day interval for large birds and hedgehogs; a 4-day interval for birds of prey; and an 8-day interval for carnivores (Fig. 2; Fig. S1, S2, and S3). To achieve a goal of 70% persistence (assuming a greater financial budget), the corresponding periodicity of road monitoring would be: daily monitoring for lizards (in the afternoon) and bats (in the early morning); daily monitoring for toads, small birds, small mammals, snakes, salamanders, and lagomorphs; alternate days for large birds, birds of prey, hedgehogs, and freshwater turtles; and every three days for carnivores (Fig. 2; Fig. S1, S2, and S3). For a general monitoring program directed at capturing all vertebrate groups, daily monitoring would be required to attain at least a 50% persistence probability across all taxa (Fig. 1). Specific periods of the day are justified for monitoring lizards and bats, because they exhibited very short persistence time (<1 day). Accordingly, monitoring surveys should be conducted as close as possible to their activity (or mortality) periods: lizards are killed during daylight hours, while bats are killed overnight.
Influence of road, animal and weather on persistence time
The Cox model for the “all taxa” data set (stratified by taxonomic group) was highly significant (P<0.001), and explained 9.9% of the variance in carcass persistence time. According to this model, carcass persistence probability was lower for animals weighing less than 20 g when compared to higher body masses; in roads with very low traffic (less than 1000 vehicles/day); in paved lanes when compared to unpaved shoulders; and during periods with higher proportion of rainy days in the wet season (November–April), and with higher minimum temperatures in the dry season (May–October; Table 2).
Table 2. Multivariate Cox hazards models for the persistence time of “all taxa” data set, stratified by taxonomic group (β: coefficients, eβ: hazards ratio, eβ LCI 95%: 95% lower confidence interval for hazards ratio, eβ UCI 95%: 95% upper confidence interval for hazards ratio, z: Wald test, P value: significance of Wald test); performance parameters of the model: R2 = 0.099; Likelihood test = 464.9, df = 14, P-value<0.0001.doi:10.1371/journal.pone.0025383.t002
As an example of interpretation of the effect of body mass on persistence probabilities, based on estimates of the hazards ratio (or risk ratio, RR), the carcass of a <20 g animal was about half as likely to persist as a 100 to 1000 g carcass (RR respectively 0.89 and 0.45), and animals weighing more than 1000 g were almost 80% more likely to persist (RR = 0.21) (Table 2).
Overall, Cox models were significant for most taxonomic groups considered, explaining 6.4 to 66.8% of variance in carcass persistence time (bats and freshwater turtles, respectively), with an average value of 23.8% (n = 12 models; Tables 3, 4, 5, 6, and 7). It was not possible to build a significant model for lizards, since none of the available explanatory variables lowered the Akaike's Information Criterion (AIC; ) comparing with the null model. For bats, the best explanatory variable was not significant itself, nor did it contribute for a significant model. Despite this, the model for bats is presented (Tables 3 and 4).
Table 3. Summary of performance parameters of Cox proportional hazards models for the 12 taxonomic groups analyzed (R2: explained variance, LL test: Likelihood test, df: degrees of freedom, P-value: significance of Likelihood test).doi:10.1371/journal.pone.0025383.t003
Table 4. Multivariate Cox hazards models for Bats, Toads, Salamanders, and Small mammals (β: coefficients, eβ: hazards ratio, eβ LCI 95%: 95% lower confidence interval for hazards ratio, eβ UCI 95%: 95% upper confidence interval for hazards ratio, z: Wald test, P-value: significance of Wald test, for variables definition, see table 8).doi:10.1371/journal.pone.0025383.t004
Table 5. Multivariate Cox hazards models for Snakes, Small birds, and Lagomorphs (β: coefficients, eβ: hazards ratio, eβ LCI 95%: 95% lower confidence interval for hazards ratio, eβ UCI 95%: 95% upper confidence interval for hazards ratio, z: Wald test, P-value: significance of Wald test, for variables definition, see table 8).doi:10.1371/journal.pone.0025383.t005
Table 6. Multivariate Cox hazards models for Freshwater turtles, Large birds, and Hedgehogs (β: coefficients, eβ: hazards ratio, eβ LCI 95%: 95% lower confidence interval for hazards ratio, eβ UCI 95%: 95% upper confidence interval for hazards ratio, z: Wald test, P-value: significance of Wald test, for variables definition, see table 8).doi:10.1371/journal.pone.0025383.t006
Table 7. Multivariate Cox hazards models for Birds of prey and Carnivores (β: coefficients, eβ: hazards ratio, eβ LCI 95%: 95% lower confidence interval for hazards ratio, eβ UCI 95%: 95% upper confidence interval for hazards ratio, z: Wald test, P-value: significance of Wald test, for variables definition, see table 8).doi:10.1371/journal.pone.0025383.t007
Results of Cox models revealed that the weather related variables (“season”, “raindays”, “rainamount”, and “mintemp”; see Table 8 for variable details) were the most important, each of which being retained in six to eight models (out of the total of 12 models), covering taxa with short to long persistence time (Tables 4, 5, 6, and 7). Higher proportion of “raindays” decreased carcass persistence probability for toads, salamanders, small mammals, snakes, lagomorphs, and carnivores (Tables 4, 5, and 7). However, a larger “rainamount” only decreased the persistence probability for salamanders (with borderline significance; P = 0.053), while increasing persistence for small mammals, snakes, small birds, lagomorphs, hedgehogs, birds of prey, and carnivores (Tables 4, 5, 6, and 7).
Table 8. List of explanatory variables, their definition, and values.doi:10.1371/journal.pone.0025383.t008
The influence of temperature also varied between groups: higher “mintemp” decreased the persistence probability for snakes, small mammals (with borderline significance; P = 0.058), and freshwater turtles; while higher “mintemp” increased the persistence for birds of prey and carnivores (Tables 4, 5, 6, and 7). The persistence probability of large birds, birds of prey and carnivores decreased during the dry season (Tables 6 and 7). Some interactions were significant, particularly season with proportion of rainy days for small birds and hedgehogs, and season with minimum temperature for small birds and salamanders (interpreted as for the “all taxa” model; Tables 4, 5, and 6).
Two variables not pertaining to weather were also retained commonly in built multivariate Cox models: “bmass” and “rposition”. Road-killed animals with higher body masses showed higher persistence probabilities for salamanders, small mammals, and small and large birds (Tables 4, 5, and 6). The carcasses of small birds located at the center of the road were less likely to persist than those found on the lanes. Persistence probability also decreased on paved shoulders (when compared with lanes) for large birds, and increased for freshwater turtles and carnivores (Tables 6 and 7). On the other hand, the persistence probability on lanes relatively to unpaved shoulders was higher for freshwater turtles, and lower for large birds and carnivores (Tables 6 and 7). For bats, the position of carcasses on the road seemed to have some importance with persistence probability increasing in road center and decreasing in paved shoulders relative to lanes, though neither of the class coefficients was significant (Table 4).
“Traffic” and “bcondition” were selected only in three and two models, respectively (although some classes in these variables were not significant). Accordingly, persistence increased along road segments with 4 000 to 10 000 vehicles/day (when compared with the lowest class) for lagomorphs; and along segments with the greatest traffic (i.e., more than 10 000 vehicles/day when compared with the lowest class) among small birds and lagomorphs. Non-intact carcasses of carnivores and snakes were more likely to persist on the road, though this result is non-significant in the latter group (Tables 5 and 7).
Carcass persistence time on the road
Most animal carcasses on roads are quickly dismembered by passing vehicles, eaten or removed by scavengers and predators, or reduced to skeletons by ants and other decomposers , , . In the present study, most carcasses remained on the road for the first day only, with some groups disappearing at high rates over this first day. Our results are in accordance with several published short-term experiences, some of them applied to road persistence time, and some to other habitat types, most involving bird carcasses. One of the few studies that has tested persistence time in open fields across several species, identified a similar range of loss values (24 to 98% of carcasses disappeared within 24 h), which were explained through scavenger activity alone . The remaining literature on birds, includes a wider variety of results, from persistence time similar to our data for birds with variable sizes (10 to 62% loss in 24 h; , ), to lower loss estimates (5 to 40% loss; –). Other authors have estimated lower persistence time for both small birds , , ,  and medium-sized birds , , ,  than those reported here. There also are persistence estimates published for amphibians and snakes on roads, which are respectively similar to and lower than the values documented in the present study , , . As described, most available studies (road and non-road habitats) generated lower probabilities of carcass persistence over time relatively to our data. A possible explanation for this could be differences in experimental design, including: shorter observation time (e.g. a few weeks or just two seasons); frequent use of dead animals that are placed artificially (rather than wild animals actually killed by vehicles in that place); and small sample sizes. Globally, these differences might lead to higher disappearance rates. Inter-study differences might also be explained on the basis of different regional conditions, such as different climates, habitats, communities of scavengers and predators, or types of roads (see discussion below).
Classification of taxonomic groups according to persistence time
Persistence time was highly variable among the different taxonomic groups analyzed in our study, from very short (less than one day) to relatively long periods (more than seven days). Besides differences in body size, certain species traits may make them more or less likely to persist on the road. For example, animals that are covered by fur, spines or scales are more resistant to vehicles passing over them than amphibians , , , though some species of amphibian (e.g. Salamandra salamandra) may remain longer on the road due to their tough skin and unpalatability , . On the other hand, some species may be removed more frequently than others by scavengers and predators , or even persons. ANTWORTH and collaborators  suggest that long, linear snakes are more readily detected and recognized on the road as food than small birds. In addition, small carcasses may have a wider range of potential scavengers than larger species do, and may be more rapidly destroyed by invertebrates, like ants , . Another source of carcass removal can be related to “human clean-up” of roads. During the course of this study, there was only occasional removal of very large carcasses by road crews. However, a few situations during field work suggest that, occasionally, persons remove carcasses from the road: intact lagomorphs and partridges recently road-killed (for eating), and carnivores and birds of prey (for taxidermy and scientific studies; authors, pers. observ.).
Monitoring frequency in road-kill studies
A weekly schedule (or wider spaced interval) for road mortality monitoring studies should be avoided, even when larger species with higher persistence probabilities are the target. If the study is not too constrained by its budget, we suggest monitoring with 2-day intervals for carnivores; alternate days for large birds, birds of prey, hedgehogs, and freshwater turtles; and daily for all other groups. According to our results, studies focusing on the broad vertebrate community should be based on daily monitoring. Longer intervals between surveys will underestimate the number of road-kills for most species, especially smaller animals (nearly 80% of road casualties; ), creating a bias towards lower kill rates for such animals. For species of small size, which include many threatened taxa (like European bats), the required monitoring frequencies are much higher than those typically found in many published studies. This may have resulted in road-kills underestimation in many cases , , which in turn may have prevented effective mitigation. efforts regarding animal road mortality. Although there have been studies incorporating daily frequencies of road-kill sampling of general taxonomic groups , , , a large number of surveys have been conducted once weekly , , ,  or even bi-weekly , , . The present study shows that sampling intervals longer than one day can seriously underestimate road-kill numbers, especially for smaller animals with observed losses of almost 60% overall, and of more than 85% and 73% for bats and lizards, respectively. STEWART  and ERRITZOE et al.  both have suggested that, ideally, road-killed small birds should be monitored 2–3 times each day. If we could relate persistence time with classical information theory, such as the Shannon-Nyquist sampling theorem, we would also suggest sampling twice a day for those smaller taxa to reach minimal accurate estimates (i.e., assuming that persistence time is an estimate for the fundamental frequency of a “carcass signal”; , ). However, sampling two times a day for long periods implies very high costs (financial and manpower), and would only be advisable for specific situations.
When surveying only larger animals, like carnivores or owls, a sampling frequency of once every 15 days is commonly found in the literature , , although even more prolonged intervals have been adopted –. This monitoring frequency can lead to high rates of carcass loss; our own numbers suggest that losses near 60 and 75% could be observed for carnivores and birds of prey, respectively.
Factors influencing persistence time among vertebrates
Several factors exerted an overall effect on the persistence probability of wildlife carcasses on roads. As expected, carcass persistence probabilities were smaller for small-sized animals, because they are crushed and torn apart faster by the continuous impact of car wheels, and are more easily removed or consumed by scavengers and predators than larger carcasses .
Carcasses found along roads with lower versus higher traffic volumes, and located on road lanes versus unpaved shoulders were less likely to persist than their counterparts. Our results regarding road lanes agree with other authors' results for birds , . Although one might expect that carcasses located on roads with high traffic volume would have lower persistence probabilities, due to larger number of vehicles passing directly over them, our data revealed opposing results. This might be explained by the additional influence of scavenger and predator activity, removing road-killed animals from roads, particularly from road locations with better access . In fact, lower volumes of traffic allow easier access of avian and mammalian predators to eat or remove dead animals . Moreover, the location of carcasses on road lanes not only makes them more susceptible to being repeatedly damaged by vehicles, it also renders them more visible to avian predators that hunt regularly along roads (, ; author's pers. observ.).
There are several species that include carrion in their diet. The most common are corvids, birds of prey, and mammalian carnivores; but communities of invertebrate decomposers also are very relevant, due to their abundance and diversity; and hedgehogs and rats also are occasional consumers . The contribution of predators and scavengers to decreasing carcass persistence probability has been well documented , , ,  and their influence is also suggested by the present work. Several avian predators have been frequently observed eating or carrying road-killed prey in the study area; such scavengers include kites (Milvus milvus and M. migrans), buzzards (Buteo buteo), carrion crows (Covus corone), and magpies (Pica pica) (author's pers. observ.). Mammalian carnivores were also detected in the vicinity of the studied roads looking for food; these include domestic cats (Felis catus), foxes (Vulpes vulpes), and Egyptian mongooses (Herpestes ichneumon) (António Mira, unpublished data). In addition, several studies report the consumption of road carrion by most of these species , , , further supporting the importance of scavenging in observed road persistence probabilities. There are also some reported cases of carcass removal by regular persons or road crews, although in our study it is incidental and directed to certain species (see text above).
The present study also documented how carcasses experiencing more humid conditions during the wet season and higher temperatures during the dry season exhibit shorter persistence time. Rainy days during the cooler months of the wet season should promote faster carcass dismemberment by vehicles passing over them, comparing with drier and cold days, which facilitate the preservation of animal tissues . Although higher temperatures during the dry season could enhance the rapid desiccation of animal bodies, thereby increasing persistence probabilities , it also can promote higher activity rates of microorganisms and invertebrate decomposers (ants, coleoptera, maggots, etc.; , ), accelerating the disappearance of carcasses. Indeed, ants can transform a recently road-killed small bird into a skeleton in just one hour (F. Carvalho, pers. observ.). Our results are in accordance with those of several other authors, who have claimed that carcass persistence is lower in summer months than in spring or autumn, due to increased temperatures and the diversity of insect communities ,  or scavenger activity , . In addition, elevated temperatures during summer increase the formation of volatile and smelly chemicals that can attract scavengers and predators to the carcasses . On the other hand, predator and scavenging activity by vertebrates can increase during the dry season due to the greater energy needs of seasonal offspring and the later abundance of juveniles , .
Regarding the influence of explanatory variables on persistence probabilities among different taxonomic groups, certain seemingly-related variables occasionally exerted opposing effects (i.e., one variable increasing persistence while another variable decreased it). For example, carcass persistence probability among several taxa decreased with the proportion of rainy days, but increased with larger amounts of rain. Although ostensibly contradictory, one explanation for this is that, although humid days favor carcass softening and dismemberment, intensive precipitation can reduce predator activity, especially hunting flights by avian predators , . On the other hand, the effect of the amount of rain could influence some taxonomic groups differently, as certain species are not among preferred prey. For example, although most amphibians tend to be road-killed on rainy days, salamanders generally are avoided by most predators  and thus remain longer on the road than most toads and frogs under this weather conditions.
Although higher temperatures decreased persistence probability among most taxonomic groups, birds of prey and carnivores remained longer on the road under higher temperature conditions. The higher body masses and dimensions of these taxa could explain such differences. Larger carcasses are unlikely to be carried away by vertebrate scavengers ,  as few predators in the area can easily remove such weights. The most probable situation is that, after soft body parts have been eaten by scavengers and predators, the high temperatures promote hardening of the remaining body fur or feathered skin, which stay on the road in a desiccated form for prolonged periods of time (author's pers. observ.).
Carcass persistence probabilities were much reduced during rainy days in the wet season for small birds and hedgehogs, and at high temperatures during the dry season for small birds and salamanders. This is in accordance with the “all taxa” model and could reflect the specific climatic conditions that favor carcass persistence within each season (as previously discussed) or the seasonal differences in predator activity and abundance in the area.
Once again, differences in persistence probabilities between road shoulders and lanes should reflect different scavenging strategies by avian and mammalian predators. The high probabilities of carnivore carcass persistence on shoulders can be explained by the marked decrease in the amount of damage caused by vehicles passing over them in this part of the road, and the lower expected range of scavengers eating larger species , . Conversely, lagomorphs frequently are taken as prey by road scavengers , or even by persons, when recently killed (authors, pers. observ.). Feeding on a road-killed carcass poses its own risks for predators (and persons when stopping and getting out of the vehicle), this risk being less along roads with lower traffic intensity. Thus, higher levels of traffic intensity promoted higher persistence probabilities for lagomorphs.
Carnivores and snakes remained longer on the road as body remains (as opposed to intact animals). For snakes, it is commonly observed in the field that intact animals are more frequently carried away by predators, while remains get fixed to the pavement and skin scales persist for a long time (; authors, pers. observ.). For carnivores, this result might be explained by the partial consumption of carcasses by vertebrate scavengers that cannot remove higher weight animals from the road , .
To obtain realistic and complete pictures of wildlife road mortality, particular roads must be monitored in a systematic way and at regular intervals over at least a year . The present study highlights the numerical consequences of choosing different time intervals for road monitoring studies, focusing, for the first time, on a wide variety of European vertebrate groups in a large-sample study. Our results suggest that, although variable among taxonomic groups, persistence time of vertebrate carcasses on the road are globally short, thus requiring that road mortality estimates should be based on high quality data collected at fine-temporal scales (one day or even shorter intervals). The guidance given here regarding monitoring frequencies is particularly useful for studies aiming the recognition of road-kill patterns, and the identification of road mortality hotspots for impact assessment and validation of mitigation measures. In fact, recent results suggest that distinct frequencies of road monitoring can lead to differences in the spatial location of road mortality hotspots (A. Mira et al. in prep). The application of mitigation measures outside real mortality hotspots may be financially costly and environmentally ineffective in reducing wildlife mortality on roads.
The present study highlights, for the first time, the effects of multiple factors on the persistence of road-killed carcasses from a large and diverse sample of wild vertebrates. Our results also suggest that scavengers and predators may exert a strong influence on carcass persistence probabilities, and that their activity (and thus removal rate) might, in turn, depend upon the aforementioned factors (season, weather, location, and animal group). Moreover, differences in the diversity and abundance of predator communities between different geographic regions should lead to distinct carcass persistence probabilities. The carcass removal by scavengers and predators should be studied further in different regions and landscape contexts because, besides differences in population abundances, scavengers and predators with different sizes, periods of activity or food preferences must affect differently the probabilities of carcass persistence.
The results reported here may have practical implications in designing effective road mortality monitoring surveys. In particular, results showed that the accuracy of road monitoring counts can be greatly improved by adapting monitoring methods and frequencies according to local weather conditions, traffic volume, and size of target species. The yearly schedule of surveys may be adjusted to account for mean expected precipitation and temperatures, so as to account for varying carcass persistence time (e.g., more frequent surveys during rainy months). Also, the timetable of surveys may be attuned to the target taxa. However, different weather conditions through the world may affect different species differently, and monitoring frequencies should be adapted. Surveys directed towards smaller species should be of higher concern and must be conducted at higher frequencies, particularly in humid seasons as these conditions generally promote the shortest persistence time on roads. This can be critical under other climate regimes, such as monsoonal rains, where persistence time should be even shorter. Traffic volume is another factor to consider, as less traffic facilitates scavenger access to carcasses. Designing road mortality monitoring surveys should thus require adjustments in sampling efforts according to target groups and local environmental conditions. Besides increasing the accuracy of estimates, the use of standardized multiple species road-kill surveys will allow the comparison of results between studies, and the improvement of mitigation proposals.
We caution that the results from our work are not universal in their application, as the mean body size of some taxonomic groups discussed here may differ in other regions of the world. For instance, most of the European bat species are of small size (under 20 g; ) when compared to South American or Australian bats (e.g. several species of flying-foxes can weigh up to 1 000 g; ). Also, our work is limited to species under 10 kg weigh, which implies that carnivore results apply to medium-sized predators (a common scenario in most European countries).
Moreover, different weather conditions in distinct geographic regions are likely to influence carcass persistence probabilities. The overall trends reported here must verify in many regions. However, specific carcass removal rates may change, tending to be higher in warmer and/or wet regions. Thus, an adjustment of proposed survey intervals may be needed in regions where weather conditions are markedly different. Even so, we believe that our recommendations are valid in many circumstances, although being particularly useful in Mediterranean areas.
Most studies documenting the negative impact of roads make inferences from the number of road-kills documented during road surveys. We have quantified biases associated with different monitoring frequencies and taxonomic groups under different environmental conditions, highlighting how the numbers counted may not accurately reflect the actual number of road-kills.
Materials and Methods
This study was conducted in southern Portugal (38°32′24″ to 38°47′33″N, 08°13′33″ to −07°55′45″W), in an area between the cities of Montemor-o-Novo (west), Évora (east) and Arraiolos (north). The relief is smooth and undulating, with elevations ranging from ca. 150 m to 400 m above sea level. The study area is dominated by the typical Mediterranean forest known as “montado” (50.8%; cork and/or holm oak trees with agro-silvo-pastoral use; ) and agricultural areas (44.0%; mainly dry arable land, sparse “montado”, and olive orchards and vineyards). Weather is characteristic of Mediterranean climates with minimum and maximum mean temperatures of 5.8°C and 12.8°C during the winter (January), respectively, and 16.3°C and 30.2°C in the summer (July); annual rainfall averages 609.4 mm (Évora 1971–2000; ). Four road sections were selected within this area, with varying traffic volumes, and summing to 37 km in total length. Roads N4 and N114 (12 km and 9.5 km, respectively) are classified as national roads, while M529 and M370 are municipal roads (9 km and 6.5 km, respectively). All roads are two-lanes wide, without central barriers/dividers, except in two road crossings. The national roads have paved and unpaved shoulders; while the municipal roads have unpaved shoulders only. All animals used in the present study were already found dead (road-killed), and therefore an ethic approval is not required. All efforts were made to minimize suffering of animals found still alive after being hit by a vehicle, delivering them as soon as possible to wildlife recovering centers.
From December 2004 to February 2006, the four roads were surveyed daily by vehicle to detect road-killed vertebrates. Surveys began within 2 h of sunrise. The standard road sampling width corresponded to both lanes and shoulders (paved and unpaved). Surveys were conducted by one observer driving 20 km per hour and checking both sides of the road. Whenever a road-killed animal was detected, the species was identified in loco (or most accurate taxonomic position for carcass remains), and the geographic coordinate position was registered with a hand-held global position system (GPS) unit with 5 m-accuracy (Garmin eTrex Venture). The body condition of the carcass (full carcass or remains), and the position on the road (center, lanes, paved or unpaved shoulder) were registered. All the carcasses were left in the same position in which they were initially found, and during subsequent surveys their presence was rechecked to determine how long they lasted before disappearing. This allowed for the calculation of the persistence time.
Other explanatory variables were added later: classes of traffic volume for each road section, mean body mass and length of each species, and average meteorological conditions during the period of carcass persistence (proportion of days with rainfall, amount of rainfall, mean daily temperature, minimum daily temperature and maximum daily temperature) (see Table 8).
Traffic intensity values for different segments of the studied roads were obtained from national road reports  and from our own data (António Mira, unpublished data). The classes of traffic intensity were defined as in IUELL et al. : roads with 1000 or fewer vehicles/day were considered permeable to most species (M370); roads with 1001 to 4000 vehicles/day were deemed permeable to most species, but avoided by sensitive species (M529); roads with 4001 to 10 000 vehicles/day were considered strong barriers associated with high road mortality (N4 and a section of N114); and roads with more than 10 001 vehicles/day were assumed to be impermeable to most species (one section of N114).
Meteorological information was obtained from the meteorological weather station located in Mitra/Évora (Geophysics Center of Évora, University of Évora). Two variables concerning rain were defined: “raindays” quantifies the duration of periods of high humidity, and “rainamount” quantifies the amount of rain that falls over a certain time period. Accordingly, prolonged periods with little rain, a common event in the study area, tend to produce relatively small quantities of effective rainfall.
The original data were categorized into 13 data matrices, corresponding to the 13 most representative taxonomic groups of road-killed vertebrates: toads (anurans, including frogs; 4 to 125 g weigh); salamanders (urodeles, including newts; 25 to 30 g weigh); lizards (lacertids; 3 to 340 g weigh); snakes (colubrids; 15 to 240 g weigh); freshwater turtles (chelonids; 290 to 300 g weigh); small birds (passeriformes, coraciiformes and piciformes; with 8 to 200 g weigh); large birds (birds weighing between 200 g and 1200 g, excluding birds of prey); birds of prey (diurnal and nocturnal; accipitriforms and strigiforms; 175 to 1100 g weigh); bats (chiropters; 6 to 23 g weigh); small mammals (arvicolids, murids, soricids and talpids; 11 to 300 g weigh); lagomorphs (Oryctolagus cuniculus and Lepus granatensis; 1100 to 2300 g weigh); hedgehogs (Erinaceus europaeus; 850 g weigh); and mammalian carnivores (111 to 7300 g weigh). Global data (including the 13 taxonomic groups) were also analyzed (referred to in the text as “all taxa”).
Median carcass persistence was estimated and compared among different taxonomic groups using the Kaplan-Meier estimator . The survival curves and probabilities produced with this method allowed the comparison of overall persistence time among the 13 taxonomic groups.
The influence of explanatory variables (road and animal characteristics, and weather conditions; see Table 8) on the persistence probabilities of carcasses was assessed by means of Cox proportional hazards models . The interpretation of both the survival functions, described by the Kaplan-Meier estimator, and the Cox model correspond, in the present context, to the probability that a carcass remains on the road until a specified time, with low hazard ratios (and negative Cox model coefficients) indicating higher persistence probabilities.
Preliminary screening of variables was undertaken with exploratory plots and simple Cox survival models . Arcsin transformation was applied to “raindays”, and logarithmic transformation to “rainamount” for all data matrices. Logarithmic transformation was applied to “bmass” for the matrices of small birds, small mammals, and birds of prey . Pearson correlations among all pairs of continuous variables (or chi-square tests for categorical variables) were calculated to check for multicollinearity . For pairs of variables exhibiting correlation values greater than 0.7  only the strongest predictor in simple Cox models was used in further analysis.
Multiple Cox survival models were constructed for each data matrix (for the “all taxa” data set and each single taxonomic groups), and selection of the best model was based on AIC . Interaction terms among weather variables “season” and “traffic” were considered during model selection. The “all taxa” model included the taxonomic group as a stratified variable in order to account for different baseline hazards for each taxon .
The overall importance of variables in the Cox models was evaluated through the log-likelihood chi-square test, while coefficients for individual variables in the models were evaluated with the Wald test and 95% confidence intervals for the hazard ratios (not including the value of 1.0). Competing models (differences in AIC values less than 4) were evaluated with the chi-square test, and most complex models were retained only if a significant (P<0.05) decrease in AIC was observed. The number of explanatory variables in each model was limited to a number such that a ratio of 7 to 10 cases per variable could be achieved . The proportional hazards assumption was tested using a global chi-square test that compares transformed survival time with scaled Schoenfeld residuals, and plots the variation of coefficients for each variable in the model through time (the Kaplan-Meier transformation was chosen as it tends to spread residuals uniformly across the plot; ). Whenever the global chi-square test of proportional hazards assumption was significant (P<0.05), the plots of problematic variables were used to assess whether the degree of temporal deviation from estimated coefficients displayed no major deviation, meaning that the proportional hazards assumption was accepted . In order to evaluate each model's goodness-of-fit, we considered the overall likelihood ratio (LR) test, the proportion of variance explained (R2), and the plots of Martingale residuals against explanatory variables . Influential observations were assessed with plots of dfbetas residuals. Interpretation of strength of association between explanatory variables and carcass persistence probabilities was based upon respective hazards ratios, and only significant variable classes (Wald test) included in the models are highlighted in the Results section.
Kaplan-Meier estimates for individual survival functions of Lizards, Bats, Toads, Salamanders, Small birds, and Snakes, showing the persistence probability and 95% confidence intervals (the length of time axis is limited to a maximum of 15 days, when available, to allow comparison between groups).
Kaplan-Meier estimates for individual survival functions of Small mammals, Lagomorphs, Large birds, Birds of prey, Hedgehogs, and freshwater turtles, showing the persistence probability and 95% confidence intervals (the length of time axis is limited to a maximum of 15 days, when available, to allow comparison between groups).
Kaplan-Meier estimates for individual survival function of Carnivores showing the persistence probability and 95% confidence intervals (the length of time axis is limited to a maximum of 15 days, when available, to allow comparison between groups).
We would like to thank Sérgio Godinho and Luis Gomes for participating in our road-kill surveys, as well as the entire MOVE team for their suggestions and comments. We appreciate the comments and suggestions to improve the English and writing style made by Rui Lourenço, Ricardo Pita and Virginia Pimenta.
Conceived and designed the experiments: AM FC. Performed the experiments: AM FC. Analyzed the data: SMS. Wrote the paper: SMS AM FC.
- 1. Iuell B, Bekker H, Cuperus R, Dufek J, Fry G, et al. (2003) Wildlife and Traffic: A European Handbook for Identifying Conflicts and Designing Solutions. European Commission.
- 2. Benítez-López A, Alkemade R, Verweij P (2010) The impacts of roads and other infrastructure on mammal and bird populations: A meta-analysis. Biol Conserv 143: 1307–1316.
- 3. Trombulak SC, Frissell C (2000) Review of ecological effects of roads on terrestrial and aquatic communities. Conserv Biol 14: 18–30.
- 4. Seiler A (2001) Ecological Effects of Roads: A review. Introductory Research Essay No 9, Department of Conservation Biology, University of Agricultural Sciences, S-730 91 Riddarhyttan, Sweden.
- 5. Forman RTT, Sperling D, Bissonette JA, Clevenger AP, Cutshall CD, et al. (2003) Road Ecology: Science and Solutions. Washington, Covelo, London: Island Press.
- 6. Grilo C, Bissonette JA, Santos-Reis M (2009) Spatial-temporal patterns in Mediterranean carnivore road casualties: Consequences for mitigation. Biol Conserv 142: 301–313.
- 7. Clevenger AP, Waltho N (2005) Performance indices to identify attributes of highway crossing structures facilitating movement of large mammals. Biol Conserv 121: 453–464.
- 8. Ascensão F, Mira A (2007) Factors affecting culvert use by vertebrates along two stretches of road in southern Portugal. Ecol Res 22: 57–66.
- 9. Mata C, Hervás I, Herranz J, Suárez F, Malo JE (2008) Are motorway wildlife passages worth building? Vertebrate use of road-crossing structures on a Spanish motorway. J Environ Manage 88: 407–415.
- 10. Clevenger AP, Chruszcz B, Gunson KE (2003) Spatial patterns and factors influencing small vertebrate fauna road-kill aggregations. Biol Conserv 109: 15–26.
- 11. Gomes L, Grilo C, Silva C, Mira A (2009) Identification methods and deterministic factors of owl roadkill hotspot locations in Mediterranean landscapes. Ecol Res 24: 355–370.
- 12. Litvaitis JA, Tash JP (2008) An Approach Toward Understanding Wildlife-Vehicle Collisions. Environ Manage 42: 688–697.
- 13. Carvalho F, Mira A (2010) Comparing annual vertebrate road kills over two time periods, 9 years apart: a case study in Mediterranean farmland. Eur J Wildlife Res 57: 157–174.
- 14. Wobeser G, Wobeser A (1992) Carcass disappearance and estimation of mortality in a simulated die-off of small birds. J Wildlife Dis 28: 548–554.
- 15. Slater FM (2002) An assessment of wildlife road casualties - the potential discrepancy betwen numbers counted and numbers killed. Web Ecol 3: 33–42.
- 16. Prosser P, Nattrass C, Prosser C (2008) Rate of removal of bird carcasses in arable farmland by predators and scavengers. Ecotox Environ Safe 71: 601–608.
- 17. Enge KM, Wood K (2002) A pedestrian road survey of an upland snake community in Florida. Southeast Nat 1: 365–380.
- 18. Coelho IP, Kindel A, Coelho AVP (2008) Roadkills of vertebrate species on two highways through the Atlantic Forest Biosphere Reserve, southern Brazil. Eur J Wildlife Res 54: 689–699.
- 19. Balcomb R (1986) Songbird Carcasses Disappear Rapidly from Agricultural Fields. Auk 103: 817–820.
- 20. Erritzoe J, Mazgajski TD, Rejt L (2003) Bird casualties on European roads- a review. Acta Ornithol 38: 77–93.
- 21. Hels T, Buchwald E (2001) The effect of road kills on amphibian populations. Biol Conserv 99: 331–340.
- 22. Barthelmess EL, Brooks MS (2010) The influence of body-size and diet on road-kill trends in mammals. Biodiv Conserv 19: 1611–1629.
- 23. Antworth RL, Pike DA, Stenvens EE (2005) Hit and Run: Effects of Scavenging on Estimates of Roadkilled Vertebrates. Southeast Nat 4: 647–656.
- 24. Burnham KP, Anderson DR (2002) Model selection and multimodel inference: a practical information–theoretic approach. New York: Springer-Verlag.
- 25. Linz GM, Bergman DL, Bleier WJ (1997) Estimating Survival of Song Bird Carcasses In Crops and Woodlots. Prairie Nat 29: 7–13.
- 26. Kostecke RM, Linz GM, Bleier WJ (2001) Survival of avian carcasses and photographic evidence of predators and scavengers. J Field Ornithol 72: 439–447.
- 27. Tobin ME, Dolbeer RA (1990) Disappearance and recoverability of songbird carcasses in fruit orchards. J Field Ornithol 61: 237–242.
- 28. Taylor BD, Goldingay RL (2004) Wildlife road-kills on three major roads in north-eastern New South Wales. Wildlife Res 31: 83–91.
- 29. Ward MR, Stallknecht DE, Willis J, Conroy MJ, Davidson WR (2006) Wild bird mortality and West Nile virus surveillance: biases associated with detection, reporting, and carcass persistence. J Wildlife Dis 42: 92–106.
- 30. Ponce C, Alonso JC, Argandona G, Garcıa Fernandez A, Carrasco M (2010) Carcass removal by scavengers and search accuracy affect bird mortality estimates at power lines. Animal Conserv. doi:10.1111/j.1469-1795.2010.00387.x.
- 31. Stewart PA (1971) Persistence of remains of birds killed on motor highways. Wilson Bul 83: 203–204.
- 32. Byrd GV, Reynolds JH, Flint PL (2009) Persistence rates and detection probabilities of bird carcasses on beaches of Unalaska Island, Alaska, following the wreck of the M/V Selendang Ayu. Marine Ornithol 37: 197–204.
- 33. Langen TA, Machniak A, Crowe EK, Mangan C, Marker DF, et al. (2007) Methodologies for Surveying Herpetofauna Mortality on Rural Highways. J Wildlife Manage 71: 1361–1368.
- 34. Kiesecker JM, Chivers DP, Blaustein AR (1996) The use of chemical cues in predator recognition by western toad tadpoles. Anim Behav 52: 1237–1245.
- 35. Lesiński G (2008) Linear landscape elements and bat casualties on roads - an example. Ann Zool Fen 45: 277–280.
- 36. Russell Al, Butchkoski CM, Saidak L, McCracken GF (2009) Road-killed bats, highway design and the commuting ecology of bats. Endang Spec Res 8: 49–60.
- 37. Dodd CK, Barichivich WJ, Smith LL (2004) Effectiveness of a barrier wall and culverts in reducing wildlife mortality on a heavily traveled highway in Florida. Biol Conserv 118: 619–631.
- 38. Ramp D, Caldwell J, Edwards KA, Warton D, Croft DB (2005) Modelling of wildlife fatality hotspots along the snowy mountain highway in New South Wales, Australia. Biol Conserv 126: 474–490.
- 39. Lodé T (2000) Effect of a Motorway on Mortality and Isolation of Wildlife Populations. AMBIO 29: 163–166.
- 40. Hell P, Plavy R, Slamecka J, Gasparik J (2005) Losses of mammals (Mammalia) and birds (Aves) on roads in the Slovak part of the Danube Basin. Eur J Wildlife Res 51: 35–40.
- 41. Petronilho JMS, Dias CMM (2005) Impact of Two Forest Roads upon Wildlife after a Road Pavement Change in a Coastal Area in the Center of Portugal. Wildlife Biol Practice 1: 128–139.
- 42. Gryz J, Krauze D (2008) Mortality of vertebrates on a road crossing the Biebrza Valley (NE Poland). Eur J Wildlife Res 54: 709–714.
- 43. Nyquist H (1928) Certain topics in telegraph transmission theory. EE Trans 47: 617–644.
- 44. Shannon CE (1949) Communications in the presence of noise. Proc IRE 37: 10–21.
- 45. Fajardo I (2001) Monitoring non-natural mortality in the barn owl (Tyto alba), as an indicator of land use and social awareness in Spain. Biol Conserv 97: 143–149.
- 46. Barrientos R, Bolonio L (2008) The presence of rabbits adjacent to roads increases polecat road mortality. Biodiv Conserv 18: 405–41.
- 47. Smith-Patten BD, Patten MA (2008) Diversity, Seasonality, and Context of Mammalian Roadkills in the Southern Great Plains. Environ Manage 41: 844–852.
- 48. Lambertucci SA, Speziale KL, Rogers TE, Morales JM (2009) How do roads affect the habitat use of an assemblage of scavenging raptors? Biodiv Conserv 18: 2063–2074.
- 49. Meunier FD, Verheyden C, Jouventin P (2000) Use of roadsides by diurnal raptors in agricultural landscapes. Biol Conserv 92: 291–298.
- 50. Witmer GW, Pipas MJ, Campbell DL (1995) Effectiveness of Search Patterns for Recovery of Animal Carcasses in Relation to Pocket Gopher Infestation Control. Int Biodeter Biodegr 36: 177–187.
- 51. Kocárek P (2003) Decomposition and Coleoptera succession on exposed carrion of small mammal in Opava, the Czech Republic. Eur J Soil Biol 39: 31–45.
- 52. Preston CR (1981) Environmental influence on soaring in wintering red-tailed hawks. Wilson Bul 93: 350–356.
- 53. Sergio F (2003) From individual behaviour to population pattern: weather-dependent foraging and breeding performance in black kites. Anim Beahav 66: 1109–1117.
- 54. Sharanowski BJ, Walker EG, Anderson GS (2008) Insect succession and decomposition patterns on shaded and sunlit carrion in Saskatchewan in three different seasons. Forensic Sci Int 179: 219–240.
- 55. Cramp S (1980) The Birds of the Western Paleartic, vol II. Oxford: Oxford University Press.
- 56. Cavallini P (1996) Variation in the social system of the red fox. Ethol Ecol Evol 8: 323–342.
- 57. Brodie ED, Smatresk NJ (1990) The antipredator arsenal of fire salamanders: spraying of secretions from highly pressurized dorsal skin glands. Herpetologica 46: 1–7.
- 58. Dietz C, von Helversen O, Nill D (2009) Bats of Britain and Europe and Northwest Africa. London: A & C Black Publishers Ltd.
- 59. Hall LS, Richards G (2000) Flying-foxes, Fruit and Blossom Bats of Australia. Sydney: UNSW Press.
- 60. Pinto-Correia T (1993) Threatened landscape in Alentejo, Portugal: the ‘montado’ and other ‘agro-silvo-pastoral’ systems. Landsc Urban Plan 24: 43–48.
- 61. IM (2010) Instituto de Metereologia website, IP, Portugal. Available: http://www.meteo.pt. Assessed 2010 July 2.
- 62. EP (2005) Recenseamento do tráfego – Évora. Estradas de Portugal, E.P.E.
- 63. Cramp S, Simmons KEL, Perrins CM, editors. (1977–1994) The birds of the Western Palearctic, vol 1–9. Oxford: Oxford University Press.
- 64. Mathias ML (1999) Guia dos Mamíferos Terrestres de Portugal Continental, Açores e Madeira. Lisboa: Instituto da Conservação da Natureza.
- 65. Palomo LJ, Gisbert J, Blanco JC (2007) Atlas y libro rojo de los mamíferos terrestres de España. Madrid: Ministerio de Medio Ambiente.
- 66. Ferrand de Almeida N, Ferrand de Almeida P, Gonçalves H, Sequeira F, Teixeira J, et al. (2001) Guia FAPAS Anfíbios e répteis de Portugal. Fundo para a Protecção dos Animais Selvagens, Câmara Municipal do Porto.
- 67. Salvador A (2010) Enciclopédia virtual de los vertebrados españoles website. Museo Nacional de Ciencias naturales (CSIC). Available: http://www.vertebradosibericos.org. Assessed 2010 April 1.
- 68. Therneau TM, Grambsch PM (2000) Modeling Survival Data: extending the Cox model. New York: Springer-Verlag.
- 69. Cox DR (1972) Regression models and life tables (with discussion). J R Stat Soc B 74: 187–220.
- 70. Kleinbaum DG, Klein M (2005) Survival Analysis, a Self-Learning Text, 2nd edition. New York: Springer.
- 71. Zuur AF, Ieno EN, Smith GM (2007) Analysing Ecological Data. New York: Springer.
- 72. Tabachnick BG, Fidell LS (2001) Using multivariate statistics. 4th ed. Boston: Allyn & Bacon.
- 73. Vaida F, Xu R (2000) Proportional hazards model with random effects. Stat Med 19: 3309–3324.
- 74. Babyak MA (2004) What You See May Not Be What You Get: A Brief, Nontechnical Introduction to Overfitting in Regression-Type Models. Psychosom Med 66: 411–421.
- 75. R Development Core Team (2009) R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Available: http://http://www.R-project.org. Assessed 2010 April 1.
- 76. Therneau T (2009) survival: Survival analysis, including penalised likelihood. R package version 2.35-8. http://CRAN.R-project.org/package=survival.