Introduction
Managing conflicts involving wildlife is a serious challenge globally, especially when the conflicts pose a direct threat to human life (Woodroffe et al., Reference Woodroffe, Thirgood and Rabinowitz2005). Over large regions of the tropics and subtropics crocodilians are responsible for more attacks on people per year than any other large carnivore, with incidents reported for 65 countries (CrocBITE, 2018). They also attack domestic animals such as livestock and dogs, damage fishing nets, burrow through earthen dam walls, and cause general disruption to movement and livelihoods by restricting the use of waterways that are required for travel, fishing and household chores (Aust et al., Reference Aust, Boyle, Fergusson and Coulson2009; Lamarque et al., Reference Lamarque, Anderson, Fergusson, Lagrange, Osei-Owusu and Bakker2009).
Successfully mitigating harmful human–crocodilian interactions has important implications not just for human safety but also for conservation. Negative impacts on people have resulted in retaliatory killings and exacerbated issues such as overharvesting and habitat destruction (Fukuda et al., Reference Fukuda, Manolis, Saalfeld and Zuur2015), which can have wider consequences for ecosystems. As the apex predator of the aquatic environments they inhabit, and as carnivores with varied diets, crocodilians play a key role in local food webs, affecting abundance and behaviour of prey species in ways that can influence the function and community structure of ecosystems (van der Ploeg et al., Reference van der Ploeg, Cauillan-Cureg, van Weerd and Persoon2011). They also transport nutrients from land to water by preying on land animals, or across ecosystems, as in cases where they travel to feed in estuarine areas before returning to freshwater areas. Crocodilians are also ecosystem engineers, modifying their environment by digging dens, holes and tunnels, and creating nest sites. Alligator holes, for instance, retain water during dry seasons and serve as refuges for other aquatic species (Britton, 2009).
Understanding the drivers of crocodilian attacks can help minimize future incidents and inform conflict management. Previous analysis of crocodilian attack patterns has focused on identifying high-risk human demographics and behaviours; there has been limited analysis of potential abiotic environmental predictors of attack patterns, such as temperature and rainfall (Fukuda et al., Reference Fukuda, Manolis, Saalfeld and Zuur2015; Brien et al., Reference Brien, Gienger, Browne, Read, Joyce and Sullivan2017; Das & Jana, Reference Das and Jana2017). There are, however, strong theoretical grounds for the importance of abiotic conditions because temperature constrains crocodilian hunting behaviour via effects on ectotherm physiology, especially metabolic function (Hutton, Reference Hutton1987; Lance, Reference Lance2003; Seebacher & Franklin, Reference Seebacher and Franklin2005). Similarly, people's propensity to expose themselves to crocodilian-inhabited waterways is probably correlated with environmental conditions. Accordingly, despite crocodilian attacks occurring year-round in many tropical regions (e.g. India and Bangladesh), documented attack patterns are highly seasonal across Australia (Brien et al., Reference Brien, Gienger, Browne, Read, Joyce and Sullivan2017) and are correlated with mean monthly rainfall and temperature in South Africa, Eswatini and the southern USA (Woodward et al., Reference Woodward, Leone, Dutton, Hord and Waller2014; Pooley, Reference Pooley2015). To date, studies have only explored the contribution of abiotic factors at a coarse temporal and spatial resolution. In particular, they have exclusively used average monthly or seasonal data, which do not capture the impact of climatic changes across days. Furthermore, these data have been used to represent regional environmental conditions, but have typically been sourced from single weather stations, limiting the ability to capture variation across space.
Here, we use environmental niche modelling methods to examine the relationships between abiotic and biotic variables and the occurrence of crocodilian attacks on people at a greater spatial and temporal resolution than in previous studies. Environmental niche modelling describes methods for inferring the occurrence of biological phenomena in environmentally heterogeneous areas and has previously been utilized to analyse human conflict with jaguars in Brazil and Mexico (Zarco-González et al., Reference Zarco-González, Monroy-Vilchis and Alaníz2013; Carvalho et al., Reference Carvalho, Zarco-González, Monroy-Vilchis and Morato2015), wolves in Iran (Behdarvand et al., Reference Behdarvand, Kaboli, Ahmadi, Nourani, Salman Mahini and Asadi Aghbolaghi2014), leopards in Pakistan (Kabir et al., Reference Kabir, Ghoddousi, Awan and Awan2014) and monkeys in India (Beisner et al., Reference Beisner, Heagerty, Seil, Balasubramaniam, Atwill and Gupta2015). Crocodilian attacks can theoretically be thought of as an entity that can only occupy a specific environmental niche defined by particular biophysical limits; thus, it is possible to use environmental niche modelling to illustrate the likelihood of attack occurrence over time in a heterogeneous environment.
We focus on two species of crocodilian that come into regular conflict with humans: the East African Nile crocodile Crocodylus niloticus in South Africa and Eswatini, which is estimated to be responsible for more human fatalities per year than all other crocodilians combined (Dunham et al., Reference Dunham, Ghiurghi, Cumbi and Urbano2010; Wallace et al., Reference Wallace, Leslie and Coulson2011), and the American alligator in Florida Alligator mississippiensis, which is one of few indigenous predators that poses a substantial threat to human life in the southern USA (Woodward et al., Reference Woodward, Leone, Dutton, Hord and Waller2014). We discuss the potential of our findings to contribute to the management of crocodilians, especially in subtropical and temperate areas where climatic conditions are similar to those in our study regions. We also highlight the benefits for human safety and conservation arising from this study, as well as the potential application of environmental niche modelling to the analysis of conflict with other species, both ectothermic and endothermic.
Methods
Attack and biophysical data
Information on attacks by Nile crocodiles in South Africa and Eswatini (1951–2016) was obtained from archival searches including the personal archive of Tony (A.C.) Pooley (conservation warden responsible for monitoring crocodile attacks in Natal Province, 1966–1984, who was active in investigations until 2003), attack records kept at the St Lucia Crocodile Centre, Mtubatuba, South Africa, provincial and national newspapers (print and online) in KwaZulu-Natal and South Africa, and the Times of Swaziland archive in Mbabane, Eswatini. Online searches used the search term ‘crocodile’ paired with ‘attack’, ‘bite’ or ‘victim’, in English and Afrikaans. Information on American alligator attacks in Florida (1971–2014) was provided by the Fish & Wildlife Conservation Commission (FWC, 2016). Incidents involving a bite were reported to the Commission by telephone, by victims or their families, medical caregivers, or law enforcement agencies. In most instances post 1970, the Commission's law enforcement officers conducted investigations and interviewed victims and witnesses. Attacks by crocodilians were excluded if they were not witnessed or lacked forensic support, if spatial or temporal data were lacking, or if they were classified as provoked.
The attack records had different spatial resolutions: geographical coordinates were available for the Nile crocodile attacks (Fig. 1), whereas for the American alligator attacks the Florida county in which the attack occurred was recorded (Fig. 2). ArcGIS 10.4 (Esri, Redlands, USA) was used to ensure Nile crocodile attack coordinates fell within the IUCN predicted distribution of Nile crocodiles in South Africa and Eswatini (IUCN, 1996), and within 1 km of a waterbody (Defence Mapping Agency, 1992).
We compiled a dataset of biotic and abiotic environmental predictors that could be spatially and temporally linked to attacks. The climatic variables were mean daily temperature and rainfall, and mean monthly rainfall (20-year means). Historical data were not available for crocodile and alligator population densities across the study regions, but human population density was included as a biotic predictor because the majority of attacks in our study areas were on local residents (Woodward et al., Reference Woodward, Leone, Dutton, Hord and Waller2014; Pooley, Reference Pooley2015). Year and month were included, to account for temporal trends.
Daily rainfall and minimum and maximum temperature data were obtained from Young et al. (Reference Young, Knapp, Inamdar, Rossow and Hankins2017). We gathered data from all weather stations within a 300 km perimeter of the distribution (IUCN, 1996) of Nile crocodiles in South Africa and Eswatini (n = 514), and all stations within Florida (n = 1,131). To account for collinearity, mean temperature was calculated from minimum and maximum values. For the Nile crocodile dataset, missing temperature data were linearly interpolated when the gap between records did not exceed 3 days, increasing coverage by 9.99%. Average monthly rainfall was calculated as means of all days within each month. Climatic variables were linked to each Nile crocodile attack record from the closest NOAA station. Temperature measurements were adjusted based on the difference between the attack elevation and station elevation (Fick & Hijmans, Reference Fick and Hijmans2017), according to the temperature lapse rate defined by the International Civil Aviation Organization (ICAO, 2017). Mean climatic values were calculated for each county in Florida from all encompassed stations and were linked to attack records according to date and county.
Human population density in South Africa was estimated for the magisterial district encompassing each attack from the closest of seven historical censuses (1951, 1960, 1970, 1980, 1991, 1996, 2001; digitized and geo-referenced by Giraut & Vacchiani-Marcuzzo, Reference Giraut and Vacchiani-Marcuzzo2009). The annual human population density of Eswatini was obtained from UN (2017). Florida county human population densities were estimated using contemporaneous data from the year of the attack provided by the Office of Economic and Demographic Research (OEDR, 2017). In addition, the total wetland area in each county, taken to be a measure of potential alligator habitat (and used when selecting background points, see below), was estimated from the National Wetlands Inventory (U.S. Fish & Wildlife Service, 2017).
The final crocodile dataset comprised 188 attacks recorded in South Africa and Eswatini between September 1951 and December 2016. The final alligator dataset comprised 335 attacks recorded in Florida between September 1971 and December 2014. Both datasets included mean daily temperature and rainfall, mean monthly rainfall, human population density, month and year as variables for modelling.
Background points
Attack records constitute presence only data, and there is no reliable measure of human–crocodilian exposure to quantify when attacks did not occur. However, when absence data are not available, models can be used with background points that are randomly sampled to characterize the null distribution of each variable in the model (Elith et al., Reference Elith, Leathwick and Hastie2008; Barbet-Massin et al., Reference Barbet-Massin, Jiguet, Albert and Thuiller2012).
Background point sampling influences the results of environmental niche models (Barbet-Massin et al., Reference Barbet-Massin, Jiguet, Albert and Thuiller2012; Phillips et al., Reference Phillips, Dudík, Elith, Graham, Lehmann, Leathwick and Ferrier2009). Therefore, we used two background sampling methods for each species to compare and assess the influence of background point sampling on the results (background sampling methods A and B). Each method was used to create datasets of 10,000 background points for Nile crocodiles and American alligators that could be sampled for model fitting. For method A we stratified background point sampling according to where attacks have been reported: 10,000 Nile crocodile attack coordinates were randomly sampled and randomly assigned dates between 1951 and 2016, and 10,000 American alligator background points were randomly sampled from Florida counties proportional to the number of reported attacks and randomly assigned dates between 1971 and 2014. For method B we sampled background points across the species' predicted ranges: 10,000 Nile crocodile background points were randomly selected as coordinates that fall on water bodies within the IUCN distribution of Nile crocodiles in South Africa and Eswatini, and 10,000 American alligator background points were sampled from all counties in Florida (Woodward et al., Reference Woodward, Leone, Dutton, Hord and Waller2014). To make this method comparable between the species, the probability of a Florida county being sampled was proportional to its total wetland area (i.e. counties with more wetland area had a proportionally higher probability of selection). Covariates were linked to background points following the same process as for attack records.
Model fitting
Spearman's Rank correlation was calculated between all variables prior to model fitting (Supplementary Table 1). We first assessed attack seasonality for both species by calculating Spearman's Rank correlations between the total number of attacks per month, and mean monthly temperature and monthly rainfall, calculated using the background datasets. We analysed attack patterns at a finer spatial and temporal resolution using boosted regression trees, which were fit to express nonlinearities in the data (Elith et al., Reference Elith, Graham, Anderson, Dudík, Ferrier and Guisan2006).
We sampled background points from the background datasets following Barbet-Massin et al. (Reference Barbet-Massin, Jiguet, Albert and Thuiller2012). Boosted regression trees had an equal ratio of background points to attacks, and were fit with 10-fold cross validation, as recommended by Miller (Reference Miller2015), and model fit was evaluated as the mean area under curve of the receiver operating characteristics (ROC) plot. Boosted regression tree model parameters that maximized ROC were selected. Each tree had a learning rate of 0.005, a tree complexity of 5 and a bag fraction of 0.5. The number of boosted trees varied for each model but were selected to minimize holdout deviance. Boosted regression tree model contribution was assessed based on the number of times the variable was selected for splitting, weighted by the squared improvement to the model as a result of each split, and averaged over all trees (Friedman & Meulman, Reference Friedman and Meulman2003). Background point sampling affects ROC and variable importance and therefore we used a bootstrapping technique in which each model was fit 100 times with independently sampled background points, and variable importance and ROC were measured across each iteration. All analyses were conducted in R 3.4.3 (R Core Team, 2018) and boosted regression trees were fit using the gbm 2.1.3 package (Ridgeway, Reference Ridgeway2010).
Results
Both species have a seasonal pattern of attack frequency (Fig. 3). There is a significant correlation between number of attacks and mean monthly temperature over the study period for alligators (ω(10) = 0.95, P < 0.001) in Florida and crocodiles (ρ(10) = 0.98, P < 0.001) in South Africa and Eswatini, and a significant correlation between number of attacks and mean monthly rainfall for alligators (ρ(10) = 0.82, P < 0.001) and crocodiles (ρ(10) = 0.95, P < 0.001).
The crocodile boosted regression trees had a mean cross-validation ROC value of 0.83 ± SD 0.02 using background dataset A, and 0.85 ± SD 0.01 using background dataset B. The alligator boosted regression trees had a mean cross-validation ROC value of 0.69 ± SD 0.02 using background dataset A, and 0.77 ± SD 0.02 using background dataset B.
With background point sampling stratified by attack locations (background dataset A), mean daily temperature was the most important predictor of attack occurrence based on model contribution for both Nile crocodiles (41.16 ± SD 3.73) and American alligators (42.97 ± SD 4.16), and this was consistent for each model fitting iteration (Fig. 4). The distributions of model contribution overlap for daily rainfall, mean monthly rainfall, human population density, year and month for both species, highlighting the variability introduced by background sampling (Fig. 4). However, on average, the second and third most important predictors of Nile crocodile attack occurrence were human population density (17.86 ± SD 2.58), and mean monthly rainfall (14.80 ± 2.71), and of American alligator attack occurrence, year (14.87 ± SD 1.98) and human population density (11.90 ± SD 1.95).
With background points sampled across the species' predicted ranges (background dataset B), mean daily temperature and human population density were the strongest predictors of attack occurrence. For American alligators, mean boosted regression tree contribution was 32.34 ± SD 3.42 and 34.85 ± SD 3.44 for daily average temperature and human population density, respectively. For Nile crocodiles, mean boosted regression tree contribution was 49.42 ± SD 3.83 and 16.93 ± SD 2.67 for mean daily temperature and human population density, respectively. The relative increase in the strength of human population density as a predictor when background points are sampled across the species' predicted ranges (background dataset B) compared with when sampling is stratified by attack location (background dataset A) is probably the result of an increase in attack likelihood variance across space. This would increase the relative importance of the covariate that is the strongest spatial predictor. In our case, this is human population density, which varies greatly across state county and magisterial district. Human population density has a non-linear relationship with attack occurrence across the predicted ranges of both species. This is probably a result of the patchy distribution of crocodilians and people, as attacks are heavily dependent on the relative local abundance of both and can only occur when there is exposure between the two. The lack of historical fine-scale human and crocodilian population density data, or other proxies for encounter rates between people and crocodilians, probably limits the accuracy of our models, and it is possible that the inclusion of further covariates would affect the relative importance of predictors.
Boosted regression trees indicate that the probability of attack occurrence increases for Nile crocodiles and American alligators as temperature increases. Boosted regression tree partial dependence plots highlight the relationship between mean daily temperature and attack likelihood after accounting for the average effects of all other variables in the model. These indicate that attack likelihood increases markedly when mean daily temperature exceeds 18 °C for both species, and that it plateaus above c. 28 °C (Fig. 5). This does not mean attacks could not occur at temperatures below 18 °C. For example, the minimum daily average temperature linked to Nile crocodile and American alligator attacks in our dataset were 10.6 °C and 12.5 °C, respectively.
Discussion
Using environmental niche modelling, we explored the influence of abiotic and biotic environmental variables on crocodilian attack occurrence at a finer spatial and temporal resolution than in previous studies. Three key findings emerged for both Nile crocodiles and American alligators. Firstly, temperature was the most important abiotic temporal predictor of attack occurrence. Secondly, attack likelihood increased markedly when mean daily temperature exceeded 18 °C. Thirdly, the probability of attacks was highest above 28 °C, where it plateaued. We discuss three non-mutually exclusive hypotheses that may explain these findings.
From the perspective of crocodilian behaviour, attack likelihood may increase with ambient temperature for physiological reasons. As ectotherms, the ability of crocodilians to digest food and to hunt is dependent on ambient temperatures (Hutton, Reference Hutton1987; Emshwiller & Gleeson, Reference Emshwiller and Gleeson1997; Seebacher & Franklin, Reference Seebacher and Franklin2005). For instance, in alligators, digestive efficiency increases significantly as body temperature rises from 25 to 31 °C (Coulson & Coulson, Reference Coulson and Coulson1986), and they stop eating when ambient temperatures drop below 16 °C (Lance, Reference Lance2003). Hossain et al. (Reference Hossain, Jaman, Ahmed, Rahman and Uddin2013) found comparable patterns for captive saltwater crocodiles. In our study, a remarkably similar relationship was evident, including comparable minimum and maximum temperature thresholds, which offers strong support for the hypothesis that ambient temperature influences attack patterns via its effect on crocodilian physiology (Seebacher & Franklin, Reference Seebacher and Franklin2005; Mazzotti et al., Reference Mazzotti, Cherkiss, Parry, Beauchamp, Rochford and Smith2016).
A related hypothesis proposes that increased attacks at higher temperatures reflect heightened crocodilian aggression during the breeding season, mediated either by the effect of testosterone on male behaviour or female nest-guarding (Pooley et al., Reference Pooley, Hines, Shield and Ross1992). However, studies on alligators indicated that testosterone was not correlated with attack frequency (Woodward et al., Reference Woodward, Leone, Dutton, Hord and Waller2014), and the limited research on hormones and breeding undertaken on captive Nile crocodiles suggests that heightened hormone levels and breeding occur during the late winter months, when attack occurrence is low (Kofron, Reference Kofron1990; J. Myburgh, pers. comm., 2018). Similarly, in the USA, limited data suggest elevated testosterone levels in the early breeding season (February–March) that decline sharply over summer (Hamlin et al., Reference Hamlin, Lowers and Gillette2011). Moreover, for Nile crocodiles, nest-guarding in most regions occurs in the hottest months, meaning many of the larger adult females are not hunting and remain on or near their nests, reducing the likelihood of attacks by these animals (Pooley et al., Reference Pooley, Hines, Shield and Ross1992; Kofron, Reference Kofron1990; Combrink et al., Reference Combrink, Jonathan, Warner and Downs2017).
From a human behavioural perspective, attack likelihood may increase with ambient temperature as a result of people spending more time in and around water in hot weather, increasing encounters with crocodilians. In South Africa and Eswatini, nearly half of all attacks occurred on weekends and holidays, suggesting that human activity patterns are important. In further support of this, data collected from our study area in southern Africa suggest that Nile crocodile attacks occur most often during activities such as crossing rivers, doing domestic chores, fishing and swimming (Pooley et al., Reference Pooley, Botha, Combrink and Powell2019), and similar observations have been made in Zambia (Wallace et al., Reference Wallace, Leslie and Coulson2011). This hypothesis is clearly applicable to temperature-dependent activities, such as swimming, but is less so for river crossing and domestic chores, which occur year-round (although they may differ in duration and method at different temperatures). Unfortunately, because of the lack of direct data on human activity levels for the full duration of the study periods, we are currently unable to quantify how changes in human behaviour contribute to attack patterns. However, this highlights an important subject for future research.
Variable importance and predictive accuracy (measured as ROC) of the models were influenced by how background points were sampled. Predictive accuracy and the relative importance of human population density as a predictor were increased for the American alligator attack model when background points were sampled from across all counties in Florida (stratified by total wetland area), as opposed to stratified by counties in which attacks have occurred. This is probably a result of spatial variation in attack likelihood between counties, which is captured by human population density because of the close association between human population density and county. This is supported by the non-linear relationship observed between attack likelihood and human population density across Florida. Spatial variation in attack likelihood is probably a result of the patchy distribution of American alligators across their predicted range and resulting variation in human–crocodilian exposure rates.
Our findings have implications for crocodilian conservation. In particular, improving our understanding of the climatic conditions under which negative interactions between people and crocodilians are more common could inform future strategies for conflict management, including educating local communities in high-risk areas about the relationship between temperature and attack occurrence. Our findings could also discourage antagonistic approaches that threaten crocodilian welfare, such as removal, extermination and nest destruction (Woodward et al., Reference Woodward, Leone, Dutton, Hord and Waller2014). Behavioural modification strategies are likely to be more practical for irregular activities that can be scheduled during the winter, such as dam repair (see CrocBITE, 2018, for reports of attacks on divers). However, the use of crocodilian-inhabited waterways for domestic chores and livelihoods could also be modified by organizing essential activities during cooler periods. As with all conflict mitigation strategies, this comes with the caveat that attacks may still be possible at colder temperatures because of individual variation in crocodilian behaviour and physiology or the risk of accidental encounters (e.g. stumbling upon a nest or treading on a crocodilian). This is especially true when temperatures are relatively low at dawn, dusk or during the night, as crocodilian activity often increases during these periods and humans or other prey are more vulnerable to attack because of poor light conditions (CrocBITE, 2018). More generally, in regions with less seasonal temperature variation, the behavioural modification strategies described above may not be effective. For example, in India and Bangladesh, mean daily temperature rarely drops below 23 °C, which is within the physiological range at which crocodilians hunt (Das & Jana, Reference Das and Jana2017).
Seasonal limitations on crocodilian attacks linked to temperature could be affected markedly by predicted trends in global warming. For instance, average temperatures in Africa are predicted to increase by as much as 5–8 °C by 2050 (Ziervogel et al., Reference Ziervogel, New, Archer van Garderen, Midgley, Taylor and Hamann2014; Munzhedzi, Reference Munzhedzi2018). This is likely to facilitate more regular hunting by crocodiles, potentially to the point that it becomes commonplace throughout the year, as it is in equatorial Africa, India and Bangladesh. Ongoing monitoring of the influence of temperature changes on crocodilian behaviour and their effect on attack patterns will be important.
Our findings also have implications beyond crocodilian conflict management. The finding that crocodilian attack patterns appear to be influenced by inherent physiological constraints that are common to all ectotherms suggests that similar methods could be applied to understand negative interactions with other species. For example, snakes commonly bite humans, their pets and livestock, and the World Health Organization (2018) estimates that c. 5 million snakebites occur annually, resulting in 81,000–138,000 deaths. This issue has been neglected, resulting in the slow development of mitigation strategies (Mohapatra et al., Reference Mohapatra, Warrell, Suraweera, Bhatia, Dhingra and Jotkar2011). However, as with crocodilians, there is both theoretical and observational evidence that behaviour in snakes is influenced by climatic factors, including temperature constraints on metabolism (Saint Girons, Reference Saint Girons1980; Wang et al., Reference Wang, Zaar, Arvedsen, Vedel-Smith and Overgaard2002). Seasonal snakebite patterns peak in the summer months in Bangladesh (Rahman et al., Reference Rahman, Faiz, Selim, Rahman, Basher and Jones2010) and Nepal (Longkumer et al., Reference Longkumer, Armstrong and Finny2017). Our methods could also be applied to local issues that affect small communities, such as attacks by Komodo dragons on people, pets and livestock in Indonesia. These attacks were reported to have affected 33.3% of people in a study of affected areas (Ardiantiono et al., Reference Ardiantiono, Imansyah and Andayani2015).
In addition, the current approach of linking attacks by endothermic predators to spatial variables such as topography and urban land use (Zarco-González et al., Reference Zarco-González, Monroy-Vilchis and Alaníz2013; Behdarvand et al., Reference Behdarvand, Kaboli, Ahmadi, Nourani, Salman Mahini and Asadi Aghbolaghi2014; Kabir et al., Reference Kabir, Ghoddousi, Awan and Awan2014; Carvalho et al., Reference Carvalho, Zarco-González, Monroy-Vilchis and Morato2015) could be expanded to include abiotic predictors. The heat dissipation hypothesis suggests that activity in large species, including dangerous predators such as the tiger Panthera tigris tigris, leopard Panthera pardus saxicolor, lion Panthera leo and hunting dog Lycaon pictus, will be limited during hotter periods because of the difficulty of thermoregulating (Acharya et al., Reference Acharya, Paudel, Neupane and Köhl2016; Smith & Kok, Reference Smith and Kok2006; Creel et al., Reference Creel, Creel, Creel and Creel2016; Ghoddousi et al., Reference Ghoddousi, Soofi, Kh Hamidi, Lumetsberger, Egli and Khorozyan2016; Robertson et al., Reference Robertson, Roodbol, Bowles, Dures and Rowcliffe2019). A recent study exploring injuries and fatalities caused by large mammals in Nepal found that attacks on humans by elephants Elephas maximus, rhinoceroses Rhinoceros unicornis and tigers were all significantly higher in the winter, potentially indicating increased aggression among both predators and large herbivores at colder temperatures when vigorous activity becomes less problematic (Acharya et al., Reference Acharya, Paudel, Neupane and Köhl2016). As with previous work on crocodilians, Acharya et al. used monthly climatic data that only offered a coarse temporal resolution. Our methods, therefore, could provide a template for the study of human–endotherm conflict.
In conclusion, our findings offer strong support for the hypothesis that seasonality of crocodilian attacks is determined predominantly by temperature in regions where fluctuations are large enough to have significant impacts on crocodilian behaviour. We provide the first evidence of this at a high spatial and temporal resolution. Our analysis demonstrates that attacks increase predictably across a given temperature range and appear to be constrained by a biologically-relevant threshold. This has the potential to inform conflict management and conservation. The strong theoretical grounding for the temperature–physiology relationship suggests that our approach could be of value in casting light on the dynamics of conflict between people and other species, including both ectotherms and endotherms.
Acknowledgements
We thank Allan Woodward and the Florida Fish & Wildlife Conservation Commission for alligator attack data. Prior to this study, SP worked on early ideas on the seasonality of crocodile attacks related to biophysical variables with Joshua Potter and Aiden Keane.
Author contributions
Conception: SP, GP; data collection: SP; environmental and spatial data collection, data preparation: GP, JJW; analyses: GP, ST; writing: TMMV, GP; revisions: all authors.
Conflicts of interest
None.
Ethical standards
This research abided by the Oryx guidelines on ethical standards, and did not involve human or animal test subjects or the collection of specimens.