- Research article
- Open Access
Assessing changing weather and the El Niño Southern Oscillation impacts on cattle rabies outbreaks and mortality in Costa Rica (1985–2016)
BMC Veterinary Researchvolume 14, Article number: 285 (2018)
Rabies is a major zoonotic disease affecting humans, domestic and wildlife mammals. Cattle are the most important domestic animals impacted by rabies virus in the New World, leading to thousands of cattle deaths per year and eliciting large economic losses. In the New World, virus transmission in cattle is primarily associated with Desmodus rotundus, the common vampire bat. This study analyses the association of weather fluctuations and the El Niño Southern Oscillation (ENSO), with the occurrence and magnitude, in terms of associated mortality, of cattle rabies outbreaks. Data from the 100 cattle rabies outbreaks recorded between 1985 and 2016 in Costa Rica were analyzed. Periodograms for time series of rabies outbreaks and the El Niño 4 index were estimated. Seasonality was studied using a seasonal boxplot. The association between epidemiological and climatic time series was studied via cross wavelet coherence analysis. Retrospective space-time scan cluster analyses were also performed. Finally, seasonal autoregressive time series models were fitted to study linear associations between monthly number of outbreaks, monthly mortality rates and the El Niño 4 index, temperature, and rainfall.
Large rabies mortality occurred towards the Atlantic basin of the country. Outbreak occurrence and size were not directly associated with ENSO, but were sensitive to weather variables impacted by ENSO. Both, ENSO phases and rabies outbreaks, showed a similar 5 year period in their oscillations. Cattle rabies mortality and outbreak occurrence increased with temperature, whereas outbreak occurrence decreased with rainfall. These results suggest that special weather conditions might favor the occurrence of cattle rabies outbreaks.
Further efforts are necessary to articulate the mechanisms underpinning the association between weather changes and cattle rabies outbreaks. One hypothesis is that exacerbation of cattle rabies outbreaks might be mediated by impacts of weather conditions on common vampire bat movement and access to food resources on its natural habitats. Further eco-epidemiological field studies could help to understand rabies virus transmission ecology, and to propose sound interventions to control this major veterinary public health problem.
Rabies is a major zoonotic disease worldwide, with a high human and domestic animal death toll. Rabies virus is a negative-sense single stranded ribonucleic acid (RNA) virus belonging to the Lyssavirus genus in the Rhabdoviridae family . Recently, it was estimated that globally canine rabies causes approximately 59,000 human deaths . Furthermore, over 99% of all virus transmission to humans comes from dogs, the rest mainly coming from bats .
The problem of humans being at risk of rabies virus transmission has existed in the New World for centuries [4, 5]. Nowadays, due to an elimination programme for dog transmitted human rabies initiated in 1983 [6, 7], the main mode of zoonotic rabies transmission in the New World is via bites from infected common vampire bats to other mammal species, including humans and livestock . This was also the case in Costa Rica, where dog rabies has not been a problem for decades .
The common vampire bat, Desmodus (D.) rotundus, is the major reservoir and vector of rabies virus in the New World [9,10,11,12,13,14]. D. rotundus bloodfeeds on mammals and transmits rabies virus, with livestock animals being greatly affected by this lethal disease . In fact, common vampire bat bites and rabies transmission are limiting factors for livestock, mainly cattle production, as shown in an economic evaluation for Mexico [6, 11]. In 1968 it was estimated that over 500,000 cattle died from bat-transmitted rabies in Latin America . Following the establishment of bat control methods and cattle rabies vaccination campaigns, rabies virus associated cattle mortality has plummeted [6, 16]. By 1983, rabies was responsible for 9904 cattle deaths, for 1831 in 1993 and for 1580 in 2006 in Latin America [6, 16]. Nevertheless, in Latin America, the decline has not been always monotonic, with surges in rabies virus associated cattle mortality, like in 2000, when 6088 cattle heads died by rabies, and in 2002, when cattle deaths totaled 3327 animals [13, 16]. The economic loss associated with cattle deaths by rabies has also been significant all over Latin America, adding to an estimate of 100 million US dollars in 1966  and between 44 and 50 million US dollars by 1986 .
The epidemiology of numerous cattle rabies outbreaks has been studied across Latin America [4, 7, 9], in places so far apart like northeastern Mexico  and southern Brazil . In Costa Rica, between 1985 and 2014, over 70 cattle rabies outbreaks were reported, including 723 cattle deaths, placing a considerable burden on the local livestock industry . All of these outbreaks have been linked by the unequivocal presence of common vampire bats, and this might be enhanced by a couple of factors. The first factor is that common vampire bat dispersal is directly linked with a preference for bovine meals, probably because fenced-in cattle are a predictable resource when compared with free-ranging natural hosts . This situation generates perfect conditions for transmission in domestic animals not observed in wildlife species or in free range grazing livestock herds . The second factor is that common vampire bats are organisms sensitive to weather changes, despite having an endothermic nature. For example, it is known that temperature influences the geographic distribution of the common vampire bat . Common vampire bats tend to be limited to low elevations in the tropics, mainly because of food availability, but also because of limitations on the maximum meal weight they can carry in flight, and the low energetic cost for keeping common vampire bat body temperature constant . Common vampire bats move between caves depending on the presence of roosts and climatic conditions . For example, high humidity seems essential for common vampire bat roosting, in dry areas common vampire bats roost in water wells and leave the wells as they become dry, concentrating in wells with water during the dry season, and dispersing to other wells as they get filled with water during the rainy season . Moreover, foraging time  and roost size  of common vampire bats seem to follow seasonal patterns. For example, common vampire bats spend more time attached to their prey during the dry season than they do during the wet season, likely ingesting more blood during the dry season . Although controversial [27, 28], the birth of common vampire bats seems to be seasonal, despite the fixed 205 day gestation period recorded for vampires. However, when vampire births occur during the rainy season, rabies transmission could be enhanced via increased bat movement . Nevertheless, little attention has been given to the potential impacts that anomalous weather conditions could have on the likelihood of cattle rabies outbreaks, despite abundant observations have shown that common vampire bat biology, and rabies transmission, could be sensitive to environmental changes.
ENSO can lead to extreme climatic conditions, such as severe droughts and floods, and more generally can alter weather patterns in a manner that creates ideal ecological conditions for disease transmission [29,30,31]. In Central America, several studies have found a differential impact of ENSO on the incidence of malaria [32, 33] and leishmaniasis , making this region ideal to test hypotheses about differential ENSO impacts on disease transmission. To the best of our knowledge, only one study has analysed rabies cases in relation to ENSO, where 416 human rabies cases mainly associated to dogs were recorded in Venezuela during 2002–2004 . This study showed an increased human rabies incidence during the cold ENSO phase (La Niña in 2004) when compared to the hot ENSO phase (El Niño in 2002) , a suggestive observation despite data limitations for statistical inferences. Thus, provided that global climate change, here defined as the increase in frequency and intensity of extreme weather events driven by phenomena like ENSO, is known to affect livestock infectious diseases , such as Rift Valley Fever, and many other infectious diseases [37, 38] and the existence of detailed records about cattle rabies outbreaks in Costa Rica  creates an ideal scenario to test the hypothesis that ENSO phases might have differential impacts on the occurrence and size of cattle rabies outbreaks. Here, it is studied if the frequency of rabies outbreaks and cattle mortality per outbreak are associated with ENSO driven changes in weather patterns in Costa Rica between 1985 and 2016.
Costa Rica is located in Central America, bordering Nicaragua to the North, Panama to the Southeast, the Pacific Ocean to the West, and the Caribbean Sea (Atlantic Ocean) to the East. This small country of 51,100 km2 has a tropical climate. According to the last agricultural census, in 2014, there were around 37,000 cattle farms with 1,300,000 cattle heads in Costa Rica . The number of cattle heads has decreased 37.5% since the last census in 1984. Several factors have driven this herd size decline, including a cattle crisis in the 1980s, mainly driven by a price reduction for meat . In addition, the creation of national parks and biological reserves since the 1970s, in accordance with Costa Rica’s land use development strategy, promoted land use change from cattle ranching to biodiversity conservation areas [41, 42].
Data on cattle rabies cases from 1985 to 2016 from Servicio Nacional de Salud Animal (SENASA), the Costa Rican Veterinary Authority were used. Their database includes the data of all registered outbreaks, the number and species of dead animals attributed to rabies, and the georeferenced location of the outbreak. Meanwhile, the geolocation of outbreaks was variable, as in some instances it was based on the centroid of the county where the outbreak happened, while in other instances it was based on the coordinates of the cattle ranch itself, measured with a Global Positioning System (GPS) .
Detailed information on rabies diagnostic testing was available from the laboratory database . During the study period diagnostic tools became more robust by the inclusion of more sensitive and specific methods. Starting in 1983 with the detection of Negri bodies in brain samples processed with Sellers staining, the diagnostic method changed to mice inoculation and immune fluorescence testing in 1993. The latter method was enhanced by the addition of several deoxyribonucleic acid (DNA) based molecular techniques in 2013, described in more detail elsewhere .
In this study, an outbreak was defined as the occurrence of at least one reported confirmed rabies virus cattle death at one farm. For the analysis, three time series were constructed using information from the 100 rabies outbreaks that were reported in cattle herds of Costa Rica between January 1985 and December 2016.
The first time series was a monthly record of cattle deaths due to rabies virus infection. The second time series recorded the monthly number of outbreaks, while the third series had the monthly mortality by rabies virus infection. In this last time series, mortality was defined as the percent of deaths, due to rabies virus infection, in the farm where the outbreak was detected, i.e., for each farm was estimated as follows: (number of recorded deaths due to rabies * 100) / number of cattle heads present in the farm. When two or more outbreaks occurred in the same month, mortality was estimated by adding deaths and total cattle heads across all affected farms.
For 33 of the outbreaks, all before 2000, the number of cattle heads in farms where deaths by rabies were recorded needed to be imputed. This was done because spatially explicit data about cattle heads in Costa Rica are only available after 2013, the time when the agricultural holdings database (Sistema de Reconocimentos Ambientales, SIREA) was established [8, 43], and also to consider that cattle density has been reduced in the country. Therefore, under the valid assumption that neighbouring farms have similar herd sizes , a fully cross-validated regression tree was employed to estimate the missing herd sizes as function of farm coordinates and altitude. Briefly, the regression tree is a quantitative method where a set of rules are derived by looking at patterns of association between the independent variables and a response variable, having the ability to capture non-linear relationships between variables when compared with linear regression . Full cross-validation means that parameters are estimated by fitting a set of models where, one at a time, each observation is left out of the model fitting and then that observation is predicted with parameters estimated from all other observations so that parameters are chosen by minimizing the difference between all predicted values and observations .
For the analysis, the definition of cold and hot ENSO phases from the National Oceanic and Atmospheric Administration (NOAA) was used . This definition is based on the Oceanic Niño Index, the 3 month running average of temperature anomalies in the Niño 3.4 region (5oN-5oS, 120o-170oW), and where anomalies are based on 30 year base periods, which are updated every 7 years. The analyses employed monthly anomalies from the El Niño 4 index  which are based on measurements from the following area: 5oN-5oS, 160°E-150oW. This selection was based on previous research [32, 34, 49, 50] that has shown this index has the highest association with weather anomalies in Central America. For temperature and rainfall, gridded data, available at the Royal Netherlands Meteorological Institute (KNMI) climate explorer , were employed. Data were downloaded for the land surface contained in the area enclosed by the following coordinates: 11.00oN-8.55oN, 86.00oW-83.00oW, which roughly corresponds to Costa Rica and contains all the recorded outbreaks between 1985 and 2016. More specifically, data from the following two databases at KNMI were used for: (i) temperature: NOAA data from the Global Historical Climatology Network version 2 and the Climate Anomaly Monitoring System (GHCN_CAMS 2 m model), with a spatial resolution of 0.5° ; (ii) rainfall: data from the NOAA Climate Data Record (CDR) of Satellite-Gauge Precipitation from the Global Precipitation Climatology Project (GPCP), V2.3, with a spatial resolution of 2.5° .
Cattle rabies outbreak and mortality cycles were studied with frequency domain time series methods, which are tools to study the cyclic behaviour of time series . Periodograms, which show the distribution of power (i.e., variance) among different frequencies, i.e., the inverse of the period for a given cycle, were estimated by taking the Fourier transform of the time series. In a periodogram, a peak indicates a dominant frequency in the cyclic behaviour of a time series . Periodograms were estimated for the annual time series rabies outbreaks and the El Niño 4 index, assuming that these two time series were stationary, i.e., with a constant mean and variance . The seasonality of the epidemiological and meteorological data was assessed using a seasonal boxplot, i.e., a graph where boxplots for the 12 months in the year are sequentially drawn . Meanwhile, the association between epidemiological and climatic time series was studied using cross wavelet coherence analysis. This is a time frequency time series analysis technique that can assess the non-stationary, i.e., changing through time, association between two time series. In this analysis, the association between cycles of different periods between two time series is studied through time, showing periods of time when the cycles are significantly coupled or not [48,49,50, 54,55,56].
Retrospective space-time cluster analyses were performed with the scan statistic . Discrete poisson models  were implemented for the study period 1985–2016. For the analyses elliptic scanning windows were used. The elliptic window is the preferred choice when looking for clusters over a space that is asymmetric, i.e., like the map of Costa Rica where land occurs along a southeast to northwest axis, and the elliptic window has the advantage of converging to a circle when there is a cluster around a focal point, or an ellipse when clusters follow a pattern similar to a line . The scan statistic was employed to search for high rates. Following the suggestion by Kulldorf et al. , the scan statistic was estimated with a scanning window whose maximum radius for the long axis of the ellipse generated an area covering half of the points where cattle rabies outbreaks were located. The number of reported rabies deaths and the population at risk for each outbreak were used in order to detect whether outbreaks were purely at random, or if clusters of high mortality existed. Four models were estimated: two models were set up with an annual aggregation of cases, and two with a monthly aggregation of cases. For the models of annual and monthly data aggregation, one model considered observations including the imputed values for the estimated populations with the regression trees, and one included only information for outbreaks where the population at risk had also been registered.
Finally, seasonal autoregressive (SAR) time series models were developed to study the linear association between the monthly number of outbreaks, the monthly mortality rates and the El Niño 4 index, temperature, and rainfall. The protocol to develop SAR models started by fitting “null models” that accounted for periodicities of the focal time series (outbreaks or mortality) by appropriately considering the correlation structure of the time series as inferred from the inspection of the auto-correlation function (ACF) and partial auto-correlation function (PACF) of the time series. The output of ACF is a plot that depicts the correlation of a time series with itself at different time lags, while the PACF output depicts a similar correlation but only considering consecutive time lags . The resulting “null models” were then used to pre-whiten, also known as filtering, the time series of the climatic covariates. This process removes any common auto-correlation structure from the climatic time series; thus, when a cross-correlation function (CCF) is estimated from the residuals of the null model and the pre-whiten time series, no spurious association — emerging from the time series having a similar auto-correlation structure — is found. The CCF is a graphical representation of the correlation between two time series at different lags . Following the inspection of the cross-correlation function, we built “full models” that considered all covariates at the significant lags (p < 0.05) with the highest correlations. These full models were then simplified by a process of backward elimination , where the model with the lowest Akaike Information Criterion (AIC) from a set of models with the same number of parameter was further simplified until a minimum AIC was reached . The AIC is a trade-off function that considers the number of parameters and the likelihood of a model, allowing the selection of models based on AIC minimization . To ease model interpretation, the mean of climatic time series were removed when fitting the SAR time series models , something that equals intercepts to the average value of the focal time series studied .
All geographic information systems procedures, mapping, regression trees, cross-wavelet analyses and the time series models were done with the statistical software R version 3.4.0. The scan cluster analyses were performed with SaTScan .
During the study period (1985–2016), there were 9 El Niño events, the two biggest events, i.e., those accounting for the most extreme ENSO fluctuations, occurring in 1998 and 2015 (Fig. 1a). During this time period, corresponding to the peak power value of approximately 0.2 cycles per year, ENSO had an approximate oscillation period of 5 years, as suggested by a periodogram of the monthly anomalies from the El Niño 4 index (Fig. 1b). During the same time period, rabies outbreaks in cattle occurred all over Costa Rica (Fig. 2a). As shown in Fig. 2a, most of the outbreaks had between 1 and 5 rabies deaths, but there were two large outbreaks of 139 and 194 deaths in 1985 and 2003, respectively (see Additional file 1: Figure S1A), which occurred in the Atlantic basin of the country. Figure 2a shows that outbreaks were very common between 1986 and 1991, 2001 and 2005 as well as between 2011 and 2015. Figure 2a also shows that most outbreaks occurred in the lowlands of Costa Rica. Additional file 2: Figure S2 is a video showing the month, year and size of each individual outbreak and ENSO phase. Further inspection of cattle rabies death records suggests that most deaths have occurred in April and July (Additional file 1: Figure S1B), which have been associated with the large outbreaks of 1985 and 2003 (Additional file 1: Figure S1A and C). Nevertheless, there is no seasonality in the number of rabies deaths recorded, based on the median values (zero) of the boxplots for cattle mortality in the whole territory of Costa Rica (Additional file 1: Figure S1D). Indeed, it seems that rabies deaths randomly happen throughout the year.
Figure 2b shows results from the annual scan spatio-temporal cluster analysis about rabies mortality in cattle. Two datasets were considered, a dataset for which denominators, i.e., cattle herd size were all known, and a dataset for which some denominators were imputed. In the two analyses, clusters of cattle rabies deaths were identified in northeastern Costa Rica, in the Atlantic basin. There was an important spatial overlap between the only cluster identified with the dataset without imputed denominators (which only spanned 2003) and Cluster 1 from the analysis with imputed denominators (which spanned from 1985 to 1999). Cluster 2 coincided temporally (from 2003 to 2015) and to a lesser extent spatially with the cluster from the analysis without imputed data. Similar results were observed when the cluster analysis was performed for monthly data of cattle rabies deaths (Additional file 3: Figure S3), where the 2003 outbreak of 194 dead cows was identified as an important cluster independently of considering the imputed denominators, and clusters were also mainly located in the Atlantic basin of Costa Rica. The calculated relative risks, where relative risk is an adimensional proportion defined as the ratio between the probability of cattle death by rabies happening inside the spatial cluster when compared with the probability of such event outside the cluster area, reflect the increased risks for cattle deaths by rabies virus infections in these identified areas when compared with the surrounding landscape.
The incidence of outbreaks, independently of the number of rabies deaths, was variable across the study period, with one to nine outbreaks per year occurring during the studied years, and with cyclical fadeouts every few years up to 2002 (Fig. 3a).
There is a periodicity of approximately 5 years in the number of annual outbreaks as suggested by the annual outbreak number time series periodogram (Fig. 3b), the same oscillation period observed for ENSO. Outbreaks occurred more frequently in the months of August, September and May (Fig. 3c), and this seems to be an incipient seasonal pattern, yet at any given month no outbreaks regularly happened (Fig. 3d).
The lack of seasonality in cattle rabies deaths, and outbreaks is in sharp contrast with the seasonal weather of Costa Rica. Temperature (Fig. 4a) and rainfall (Fig. 4b) show a clear unimodal pattern, with two seasons. Temperature reaches a maximum in April and is lowest in December (Fig. 4a). For rainfall, there is a dry season between December and April, followed by a wet season the rest of the year (Fig. 4b). Nevertheless, mortality associated with cattle rabies outbreaks (Fig. 4c), although proportionally higher in May, August and September in specific years, is not regularly higher in any given month of the year.
The signature of ENSO in temperature (Fig. 5a) suggests that the hot ENSO phase is associated with higher than average temperatures in Costa Rica. Meanwhile, the cold phase of ENSO is associated with increased rainfall in Costa Rica (Fig. 5b). After 1998, outbreaks of cattle rabies seemed to be more prone to occur during the hot and cold phases of ENSO (Fig. 5c), and some cattle rabies outbreaks with high mortality rates occurred during the cold phase of ENSO (Fig. 5d).
The cross wavelet coherence analysis suggests that interannual cycles of temperature and ENSO are highly coherent at periods of 2–5 years over the studied period (Fig. 6a). For rainfall (Fig. 6b), cycles of 3–6 years of ENSO and rainfall were highly coherent around the 1998 ENSO event, and have become increasingly coherent for the same periods (2–4 years) since 2006. Cattle rabies outbreaks have been coherent at larger time scales, 8–10 years, with ENSO (Fig. 6c), and cycles have become increasingly associated at scales of 2–4 years after 2010. Monthly cattle rabies mortality shows no major significant association with ENSO during the studied period (Fig. 6d). The cross wavelet coherence analysis suggests that number of outbreaks (Fig. 6e) and cattle rabies mortality (Fig. 6f) between 1998 and 2005, for interannual cycles of periods of 4–6 years, were highly associated with temperature. Moreover, for outbreak numbers, the association was also significant after 2010 for cycles of periods of 2–4 years (Fig. 6e). Rainfall was associated with number of outbreaks (Fig. 6g) for cycles of periods of 4–6 years from 2006 to 2008. Meanwhile, there was no major significant association between cycles in rainfall and mortality (Fig. 6h) during the study period.
The number of outbreaks time series was not strongly autocorrelated, i.e., there was no clear pattern of significantly (p < 0.05) decreasing auto-correlation with time lag (Additional file 4: Figure S4A), and showed periodicities every 2 months (Additional file 4: Figure S4B), which was considered when fitting a null SAR time series model (Additional file 5: Table S1). Moreover this time series had no significant linear association with ENSO 4 over the studied period and up to 22 months of lag (Additional file 4: Figure S4C). A similar pattern of low auto-correlation was observed for the outbreak mortality time series (Additional file 4: Figure S4D), which showed periodicities of 2 months (Additional file 4: Figure S4E), also considered when fitting a null SAR time series model (Additional file 6: Table S2). The mortality time series also had no significant association with ENSO 4 up to 22 months of lag (Additional file 4: Figure S4F). When considering covariates to develop a full model for the number of outbreaks time series (Additional file 5: Table S1), lags 4 and 17 of temperature were considered (Fig. 7a) and lags 6, 11 and 16 of rainfall (Fig. 7b). Meanwhile, the full model for the mortality rate only considered lag 17 of temperature as covariate (Fig. 7c), since no significant (p < 0.05) association was found with rainfall (Fig. 7d).
The best time series model for the number of outbreaks (Table 1) had temperature with 17 months of lag and rainfall with 16 months of lag as covariate, having a positive coefficient for temperature and a negative one for rainfall. For cattle rabies mortality, the best model considered temperature with 17 months of lag as covariate (Table 2) with a positive coefficient. These results indicate that outbreaks and mortality are more likely to increase following hot-spells. Outbreaks are also less prone if the hot spell is followed by rain in an amount that offsets the positive impact of hot temperatures from the previous time lag (Table 1).
The similar period, around 5 years, of interannual oscillations in ENSO and cattle rabies outbreaks is suggestive of an association between changes in the occurrence of rabies outbreaks and ENSO phases. This type of association has been widely found for several diseases, e.g., malaria , dengue , cutaneous leishmaniasis [34, 49] in Central America, but also elsewhere, e.g. Asia , Africa [29, 63, 64] and Oceania . Although this formal analysis did not find signatures of ENSO on outbreak occurrence and cattle mortality dynamics, these two time series were associated with temperature and rainfall (outbreak occurrence only). As further shown by the analysis, the dynamics of temperature and rainfall in the studied area had interannual cycles associated with ENSO, thus implying the observed dynamics might be sensitive to ENSO. In the studied data, no seasonality in the number of monthly rabies deaths and recorded outbreaks was observed.
Clusters identified with the scan statistic were mainly located in the Atlantic basin of Costa Rica, something suggested in a preceding description of most of the studied dataset . The largest outbreaks and high mortality clusters occurred in the Atlantic basin of Costa Rica, a fact that might be related with common vampire bat ease of movement or dispersal, in an area slightly more humid  than the Pacific basin of Costa Rica . Similarly, cattle rabies geolocation patterns might reflect common vampire bat abundance patterns, which are known to be associated with specific environments, defined by the niche of bats [13, 67, 68], and where vampire abundance could increase rabies virus transmission [14, 69]. An additional reason could be that surveillance and prevention might be less effective in the Atlantic basin of Costa Rica. For example, underreporting could delay appropriate action and promote a reduced coverage of cattle rabies vaccination, which in turn might lead to larger outbreaks, given an increased susceptibility in host populations without rabies vaccination, as predicted by mathematical models dealing with risk perception and infectious disease transmission prevention . Actually, the large 1985 outbreak occurred when no control programme was in operation. In 2003, high cattle mortality was likely fuelled by late reports to the veterinary authority, as well as, a difficult access to the study area, a situation that was exacerbated by heavy rainfall which delayed, by at least a month, cattle vaccination. We suspect common vampire bats were exceptionally abundant, and likely linked with a high rabies virus prevalence. The fact that the only rabies virus isolation from a common vampire bat in Costa Rica was made where the 2003 outbreak occurred might support this further, but pure coincidence can not be ruled out as frequency of bat testing for rabies virus is not available.
This analysis indicated that temperature might have a positive impact on the transmission of rabies virus to cattle, as suggested by the positive association of temperature with outbreak occurrence and cattle rabies mortality. The lag of this association roughly corresponds to two common vampire bat generations (16 to 17 months), and is in accordance with observations about increases in survival and recruitment of common vampire bats following hot temperatures . Specifically, common vampire bat pre-weaning time is around 3 months, a period when vampire pups mainly feed on milk . Suckling vampire pups are more likely to survive at higher temperatures given the increased blood foraging in hotter/dryer conditions by adult common vampire bats . Here, it is important to note that food deprivation for two to 3 days (e.g. through adverse weather conditions such as strong rains) might become fatal for young bats , which also have a high failure rate to obtain blood meals on their own . Therefore, any change that increases the survival of young vampire pups might be critical to increase rabies virus transmission by increasing the size of common vampire bat populations. Pre-weaning vampire pups also start to receive regurgitated blood during their first 3 months of age, a key event for rabies virus transmission within common vampire bats [71, 74]. Indeed, food deprivation in common vampire bats can be compensated through reciprocal blood sharing by roost mates . Thus, a likely increase in vampire pup survival during the pre-weaning period, followed by the six to 7 months common vampire bats need to reach adulthood, when bats are between 9 and 10 months old, and also become sexually mature, with a diet exclusively based on fresh blood, could lead to a larger cohort of reproducing common vampire bats [14, 75]. This larger cohort then takes seven additional months (thus adding up to around 16 to 17 months since birth) for the gestation of a new D. rotondus vampire cohort  which will require additional blood as they try to raise their pups. Most likely this will increase blood foraging pressure, thus, rendering plausible an increase in rabies virus transmission 16–17 months after hotter temperatures might have promoted breeding and increased the survival of pre-weaning vampire pups.
Meanwhile, a mechanism where the delays are due to changes in the rabies virus incubation period in cattle or common vampire bats is unlikely, provided experimental studies have shown incubation periods in the order of one to 2 months in cattle , and two to 4 weeks in bats . Nevertheless, rabies virus incubation period is known to be highly variable depending upon inoculation site  and viral loads . The variability in rabies virus incubation period could be one of the reasons behind the lack of a seasonality in cattle rabies outbreaks, although the occurrence of large outbreaks where the onset of clinical signs is more or less synchronous suggests that, at least in outbreak foci, cattle have a homogenous incubation period . The occurrence time for the two largest cattle rabies outbreaks, which happened during April and July, and the frequent occurrence of outbreaks during the months of August, September and May, further suggest that common vampire bat ecology plays a major role explaining the lagged effects of temperature on cattle rabies outbreaks. For example, Turner  found that most bats were pregnant or lactating in the months of April–May and July–October, an observation coinciding with records from Costa Rica , suggesting increased transmission might be related to maternal care of pre-weaning common vampire bats, a possibility supported by the increased rabies virus exposure of juvenile and subadult common vampire bats, which have shown higher rabies virus seroprevalence than adults .
In contrast to the positive impacts of temperature on outbreak occurrence and size, this analysis showed that rainfall had a negative impact on cattle rabies outbreak frequency in the SAR time series models. This negative association might emerge from rainfall impacts on bat foraging. For example, a study found that rainfall increased flight metabolism, while reducing or ceasing foraging activities . Reductions in vampire foraging time could then diminish the exposure of cattle to rabies virus infected common vampire bats. Similarly, the location of cattle by common vampire bats could become more difficult following rainy periods, when a reduction in cattle surface body temperature could interfere with common vampire bat thermoreceptors used to locate their prey . Furthermore, abundant rain after a hot spell might diminish vampire fecundity, via a reduction in common vampire bat movement and mating between roosts . For example, a study in Costa Rica showed that common vampire bats stayed away from their roost more time during the dry than the wet season .
Limitations of this study include the following: (i) The mortality analysis considered estimated denominators in a fraction of the data. Nevertheless, as shown by the scan cluster analysis, the results were robust when those “estimated” denominators were ignored. (ii) Data on cattle mortality were less robust compared to outbreak data which may hide the true number of rabies cases. (iii) The impact of rabies vaccination on cattle mortality patterns was not assessed due to lacking systematic information about vaccination rates in affected herds. (iv) It was assumed that rabies virus was only transmitted by common vampire bats [12, 20]. However, it is possible that other wildlife species might have been involved in rabies virus transmission, as have been documented for raccoons, skunks and foxes [5, 76, 83,84,85]. In Costa Rica, a cattle rabies case was associated with the rabies virus strain of Tadarida brasiliensis, an insectivorous bat species . In that sense, sequencing rabies viruses isolated in cattle rabies outbreaks might help to identify other species potentially involved in rabies virus transmission. Phylogenetic analyses could also help to better understand geographical and temporal structures of outbreaks, and whether common vampire bats and other wildlife host species are involved in local enzootic cycles or epizootic waves [73, 83]. Further efforts are also necessary to understand rabies virus spillovers across species and the role that bat immunity has on the persistence of rabies virus among common vampire bats under changing environments, to fully untangle the impacts that changing rainfall and temperature patterns have on common vampire bat population structure and dynamics, as well as, rabies virus transmission.
Finally, the results of this study clearly indicate that the occurrence, and mortality of cattle rabies outbreaks are associated with weather fluctuations in Costa Rica. Further eco-epidemiological field studies are necessary to articulate the mechanisms behind the associations observed in this study, especially regarding the hypothesis that outbreaks, and their magnitude, might be driven by changes in the population dynamics of vampire bats triggered by ENSO phases. The resulting information from those studies will be useful not only to understand rabies virus transmission ecology, but also to propose appropriate intervention measures to control this major veterinary public health problem.
Further efforts are necessary to articulate the mechanisms underpinning the association between weather changes and cattle rabies outbreaks. One hypothesis is that exacerbation of cattle rabies outbreaks might be mediated by impacts of weather conditions on common vampire bat movement and access to food resources on its natural habitats. Further eco-epidemiological field studies could help to understand rabies virus transmission ecology, and to propose sound interventions to control this major veterinary public health problem.
Akaike Information Criterion
Climate Data Record
- D. :
El Niño Southern Oscillation
- GHCN_CAMS 2 m model:
Global Historical Climatology Network version 2 and the Climate Anomaly Monitoring System
Global Precipitation Climatology Project
Global Positioning System
Royal Netherlands Meteorological Institute
Partial auto-correlation function
Software for spatial, temporal, or space-time scan statistics
Servicio Nacional de Salud Animal
Sistema de Reconocimentos Ambientales
Schnell MJ, McGettigan JP, Wirblich C, Papaneri A. The cell biology of rabies virus: using stealth to reach the brain. Nat Rev Micro. 2010;8(1):51–61.
Hampson K, Coudeville L, Lembo T, Sambo M, Kieffer A, Attlan M, Barrat J, Blanton JD, Briggs DJ, Cleaveland S, et al. Estimating the global burden of endemic canine rabies. PLoS Negl Trop Dis. 2015;9(4):e0003709.
Jackson AC. Human rabies: a 2016 update. Curr Infect Dis Rep. 2016;18(11):38.
Schneider MC, Romijn PC, Uieda W, Tamayo H, DFd S, Belotto A, JBd S, Leanes LF. Rabies transmitted by vampire bats to humans: an emerging zoonotic disease in Latin America? Rev Panam Salud Publica. 2009;25(3):260.
Krebs JW, Wilson ML, Childs JE. Rabies- epidemiology, prevention, and future research. J Mammal. 1995;76(3):681–94.
Belotto A, Leanes L, Schneider M, Tamayo H, Correa E. Overview of rabies in the Americas. Virus Res. 2005;111(1):5–12.
Vigilato MAN, Cosivi O, Knöbl T, Clavijo A, Silva HMT. Rabies update for Latin America and the Caribbean. Emerg Infect Dis. 2013;19(4):678–9.
Hutter SE, Brugger K, Sancho Vargas VH, González R, Aguilar O, León B, Tichy A, Firth CL, Rubel F. Rabies in Costa Rica: documentation of the surveillance program and the endemic situation from 1985 to 2014. Vector-Borne Zoonotic Dis. 2016;16(5):334–41.
Ellison JA, Gilbert AT, Recuenco S, Moran D, Alvarez DA, Kuzmina N, Garcia DL, Peruski LF, Mendonça MT, Lindblade KA, et al. Bat rabies in Guatemala. PLoS Negl Trop Dis. 2014;8(7):e3070.
Pawan JL. The transmission of paralytic rabies in Trinidad by the vampire bat (Desmodus rotundus Murinus Wagner, 1840). Ann Trop Med Parasitol. 1936;30(1):101–30.
Anderson A, Shwiff S, Gebhardt K, Ramírez AJ, Shwiff S, Kohler D, Lecuona L. Economic evaluation of vampire bat (Desmodus rotundus) rabies prevention in Mexico. Transbound Emerg Dis. 2014;61(2):140–6.
Johnson N, Aréchiga-Ceballos N, Aguilar-Setien A. Vampire bat rabies: ecology, epidemiology and control. Viruses. 2014;6(5):1911–28.
Lee DN, Papeş M, Van Den Bussche RA. Present and potential future distribution of common vampire bats in the Americas and the associated risk to cattle. PLoS One. 2012;7(8):e42466.
Streicker DG, Recuenco S, Valderrama W, Gomez Benavides J, Vargas I, Pacheco V, Condori Condori RE, Montgomery J, Rupprecht CE, Rohani P, et al. Ecological and anthropogenic drivers of rabies exposure in vampire bats: implications for transmission and control. Proc R Soc B Biol Sci. 2012;279(1742):3384–92.
Mayen F. Haematophagous bats in Brazil, their role in rabies transmission, impact on public health, livestock industry and alternatives to an indiscriminate reduction of bat population. J Vet Med B Infect Dis Vet Public Health. 2003;50(10):469–72.
Arellano-Sota C. Vampire bat-transmitted rabies in cattle. Rev Infect Dis. 1988;10:S707.
Steele J. International aspects of veterinary medicine and its relation to health nutrition and human welfare. Mil Med. 1966;131(9 P 1):765–78.
Acha PN, Szyfres B. Zoonosis y enfermedades transmisibles comunes al hombre ya los animales. Pub Científica. 2003;284–92.
Martínez-Burnes J, López A, Medellín J, Haines D, Loza E, Martínez M. An outbreak of vampire bat-transmitted rabies in cattle in northeastern Mexico. Can Vet J. 1997;38(3):175–7.
Kobayashi Y, Ogawa A, Sato G, Sato T, Itou T, Samara SI, Carvalho AAB, Nociti DP, Ito FH, Sakai T. Geographical distribution of vampire bat-related cattle rabies in Brazil. J Vet Med Sci. 2006;68(10):1097–100.
Voigt CC, Kelm DH. Host preference of the common vampire bat (Desmodus rotundus; Chiroptera) assessed by stable isotopes. J Mammal. 2006;87(1):1–6.
Wallace R. Big farms make big flu: dispatches on influenza, agribusiness, and the nature of science: NYU Press; 2016.
McNab BK. Energetics and the distribution of vampires. J Mammal. 1973;54(1):131–44.
Trajano E. Movements of cave bats in Southeastern Brazil, with emphasis on the population ecology of the common vampire bat, Desmodus rotundus (Chiroptera). Biotropica. 1996;28(1):121.
Lord RD. Seasonal reproduction of vampire bats and its relation to seasonality of bovine rabies. J Wildl Dis. 1992;28(2):292–4.
Young AM. Foraging of vampire bats (Desmodus rotundus) in Atlantic wet lowland Costa Rica. Revista de Biologia Tropical. 1971;18:73–88.
Wainwright M. The Nat Hist of Costa Rican mammals. Florida, US: Zona Tropical; 2002.
Turner DC: The vampire bat; a field study in behavior and Ecology 1975.
Anyamba A, Linthicum KJ, Small JL, Collins KM, Tucker CJ, Pak EW, Britch SC, Eastman JR, Pinzon JE, Russell KL. Climate teleconnections and recent patterns of human and animal disease outbreaks. PLoS Negl Trop Dis. 2012;6(1):e1465.
Patz JA, Campbell-Lendrum D, Holloway T, Foley JA. Impact of regional climate change on human health. Nature. 2005;438(7066):310–7.
Chretien J-P, Anyamba A, Small J, Britch S, Sanchez JL, Halbach AC, Tucker C, Linthicum KJ: Global climate anomalies and potential infectious disease risks: 2014-2015. PLoS Currents 2015. doi: https://doi.org/10.1371/currents.outbreaks.95fbc4a8fb4695e049baabfc2fc8289f.
Hurtado LA, Cáceres L, Chaves LF, Calzada JE. When climate change couples social neglect: malaria dynamics in Panamá. Emerging Microbes & Infections. 2014;3(4):e27.
Hurtado LA, Calzada JE, Rigg CA, Castillo M, Chaves LF. Climatic fluctuations and malaria transmission dynamics, prior to elimination, in Guna Yala República de Panamá. Malaria Journal. 2018;17(1):85.
Yamada K, Valderrama A, Gottdenker N, Cerezo L, Minakawa N, Saldaña A, Calzada JE, Chaves LF. Macroecological patterns of American cutaneous leishmaniasis transmission across the health areas of Panamá (1980–2012). Parasite Epidemiology and Control. 2016;1(2):42–55.
Rifakis PM, Benitez JA, Rodriguez-Morales AJ, Dickson SM, De-La-Paz-Pineda J. Ecoepidemiological and social factors related to rabies incidence in Venezuela during 2002-2004. International Journal of Biomedical Science : IJBS. 2006;2(1):1–6.
Bett B, Kiunga P, Gachohi J, Sindato C, Mbotha D, Robinson T, Lindahl J, Grace D. Effects of climate change on the occurrence and distribution of livestock diseases. Preventive veterinary medicine. 2017;137:119–29.
Patz JA, Olson SH. Climate change and health: global to local influences on disease risk. Ann Trop Med Parasitol. 2006;100:535–49.
Liang L, Gong P. Climate change and human infectious diseases: a synthesis of research findings from global and spatio-temporal perspectives. Environ Int. 2017;103:99–108.
VI Censo Nacional Agropecuario. [http://www.inec.go.cr]. Accessed 21 Dec 2017.
León Sáenz J. Historia económica de Costa Rica en el Siglo XX. Vol.2. La economía rural. San Jose: Universidad de Costa Rica; 2012.
Rosero-Bixby L, Palloni A. Population and deforestation in Costa Rica. Popul Environ. 1998;20(2):149–85.
Arturo Sánchez-Azofeifa G, Daily GC, Pfaff ASP, Busch C. Integrity and isolation of Costa Rica's national parks and biological reserves: examining the dynamics of land-cover change. Biol Conserv. 2003;109(1):123–35.
SIREA (Sistema Integrado de Registro de Establecimientos Agropecuarios). [Information provided by Bernal Leon].
Olden Julian D, Lawler Joshua J, Poff NL. Machine learning methods without tears: a primer for ecologists. Q Rev Biol. 2008;83(2):171–93.
Kuhn M, Johnson K. Applied predictive modeling. New York: Springer; 2013.
Cold & Warm Episodes by Season. [www.cpc.ncep.noaa.gov/products/analysis_monitoring/ensostuff/ensoyears.shtml]. Accessed 21 Dec 2017.
NOAA Climate Prediction Center. [http://www.cpc.ncep.noaa.gov/data/indices/ersst3b.nino.mth.81-10.ascii]. Accessed 21 Dec 2017.
Maraun D, Kurths J. Cross wavelet analysis: significance testing and pitfalls. Nonlinear Proc Geophys. 2004;11:505–14.
Chaves LF, Calzada JE, Valderrama A, Saldaña A. Cutaneous leishmaniasis and sand fly fluctuations are associated with El Niño in Panamá. PLoS Negl Trop Dis. 2014;8(10):e3210.
Chaves LF, Pascual M. Climate cycles and forecasts of cutaneous leishmaniasis, a nonstationary vector-borne disease. PLoS Med. 2006;3(8):e295.
KMNI Climate Explorer. [http://climexp.knmi.nl/start/cgi]. Accessed 21 Dec 2017.
GHCN_CAMS Gridded 2m Temperature (Land). [https://www.esrl.noaa.gov/psd/data/gridded/data.ghcncams.html]. Accessed 21 Dec 2017.
GPCP Version 2.3 Combined Precipitation Data Set. [https://www.esrl.noaa.gov/psd/data/gridded/data.gpcp.html]. Accessed 21 Dec 2017.
Shumway RH, Stoffer DS. Time series analysis and its applications. In: Springer Texts in Statistics. New York: Springer New York; 2011.
Torrence C, Compo G. A practical guide to wavelet analysis. Bull Am Meteor Soc. 1998;79:61–78.
Cazelles B, Chavez M, Berteaux D, Ménard F, Vik J, Jenouvrier S, Stenseth N. Wavelet analysis of ecological time series. Oecologia. 2008;156(2):287–304.
Kulldorff M, Rand K, Gherman G, Williams G, DeFrancesco D. SaTScan–software for the spatial and space–time scan statistics, version 2.1. Bethesda, Madison, National Cancer Institute; 1998.
Kulldorff M, Heffernan R, Hartman J, Assunção R, Mostashari F. A space–time permutation scan statistic for disease outbreak detection. PLoS Med. 2005;2(3):e59.
Kulldorff M, Huang L, Pickle L, Duczmal L. An elliptic spatial scan statistic. Stat Med. 2006;25(22):3929–43.
Faraway JJ: Linear models with R: CRC press; 2014.
Fuller D, Troyo A, Beier J. El Nino southern oscillation and vegetation dynamics as predictors of dengue fever cases in Costa Rica. Environ Res Lett. 2009;4(1):014011.
Chuang T-W, Chaves LF, Chen P-J. Effects of local and regional climatic fluctuations on dengue outbreaks in southern Taiwan. PLoS One. 2017;12(6):e0178698.
Chaves LF, Satake A, Hashizume M, Minakawa N. Indian ocean dipole and rainfall drive a Moran effect in East Africa malaria transmission. J Infect Dis. 2012;205(12):1885–91.
Chaves LF, Hashizume M, Satake A, Minakawa N. Regime shifts and heterogeneous trends in malaria time series from western Kenya highlands. Parasitology. 2012;139:14–25.
Chaves LF, Kaneko A, Taleo G, Pascual M, Wilson ML. Malaria transmission pattern resilience to climatic variability is mediated by insecticide-treated nets. Malar J. 2008;7(1):100.
Vargas G: Geografía de Costa Rica: EUNED; 2006.
Ramoni-Perazzi P, Muñoz-Romo M, Chaves LF, Kunz TH. Range prediction for the giant fruit-eating bat, Artibeus amplus (Phyllostomidae: Stenodermatinae) in South America. Stud Neotropical Fauna Environ. 2012;47(2):87–103.
Escobar LE, Peterson AT, Papeş M, Favi M, Yung V, Restif O, Qiao H, Medina-Vogel G. Ecological approaches in veterinary epidemiology: mapping the risk of bat-borne rabies using vegetation indices and night-time light satellite imagery. Vet Res. 2015;46(1):92.
Delpietro HA, Marchevsky N, Simonetti E. Relative population densities and predation of the common vampire bat (Desmodus rotundus) in natural and cattle-raising areas in north-East Argentina. Preventive Veterinary Medicine. 1992;14(1–2):13–20.
Predescu M, Levins R, Awerbuch-Friedlander T. Analysis of a nonlinear system for community intervention in mosquito control. Discrete and Continuous Dynamical Systems-Series B. 2006;6(3):605–22.
Schmidt U, Manske U. Die Jugendentwicklung der Vampirfledermaeuse. Z Saeugetierk. 1973;38:14733.
Freitas MB, Welker AF, Millan SF, Pinheiro EC. Metabolic responses induced by fasting in the common vampire bat Desmodus rotundus. J. Comp. Physiol. B. 2003;173(8):703–7.
Wilkinson GS. Reciprocal food sharing in the vampire bat. Nature. 1984;308(5955):181–4.
Reid F. A field guide to the mammals of Central America and Southeast Mexico. Oxford: Oxford University Press; 1997.
Greenhall AM, Schmidt U. Natural history of vampire bats. Florida: CRC Press; 1988.
Schmidt U. Die Tragzeit der Vampirfledermäuse (Desmodus rotundus). Z Saugetierk. 1974;3:9129–32.
Martell M, Di Batalla C, Baer G, Acuna J. Experimental bovine paralytic rabies—“derriengue”. Vet Rec. 1974;95(23):527–30.
Moreno JA, Baer GM. Experimental rabies in the vampire bat. The American Journal of Tropical Medicine and Hygiene. 1980;29(2):254–9.
Almeida MF, Martorelli LFA, Aires CC, Sallum PC, Durigon EL, Massad E. Experimental rabies infection in haematophagous bats Desmodus rotundus. Epidemiol Infect. 2005;133(3):523–7.
Hudson LC, Weinstock D, Jordan T, Bold-Fletcher NO. Clinical features of experimentally induced rabies in cattle and sheep. J Veterinary Med Ser B. 1996;43(1–10):85–95.
Voigt CC, Schneeberger K, Voigt-Heucke SL, Lewanzik D: Rain increases the energy cost of bat flight. Biol Lett 2011:rsbl20110313.
Kürten L, Schmidt U. Thermoperception in the common vampire bat (Desmodus rotundus). J. Comp. Physiol. A Neuroethol. Sens. Neural Behav. Physiol. 1982;146(2):223–8.
Smith DL, Lucey B, Waller LA, Childs JE, Real LA. Predicting the spatial dynamics of rabies epidemics on heterogeneous landscapes. Proc Natl Acad Sci U S A. 2002;99(6):3668–72.
Badilla X, Pérez-Herra V, Quirós L, Morice A, Jiménez E, Sáenz E, Salazar F, Fernández R, Orciari L, Yager P. Human rabies: a reemerging disease in Costa Rica? Emerg Infect Dis. 2003;9(6):721.
Escobar L, Peterson AT, Favi M, Yung V, Medina-Vogel G. Bat-borne rabies in Latin America. Rev Inst Med Trop Sao Paulo. 2015;57:63–72.
Data Sources & Accuracy for National Park Service Maps. [https://www.nps.gov/hfc/carto/data-sources.cfm]. Accessed 21 Dec 2017.
The authors wish to acknowledge the Ministry of Agriculture and Livestock of Costa Rica, as well as the SENASA officials which contributed to the rabies surveillance programme and made this study possible.
Availability of data and materials
Climatic data were taken from open access databases and referenced in the manuscript. Aggregated data on rabies are included in the manuscript. The detailed datasets related to rabies analysed during the current study are not publicly available due to confidentiality reasons but are available from the corresponding author on reasonable request.
Not applicable. Data on rabies were collected within the national surveillance programme and collected in compliance with national guidelines. The Costa Rican law permits the veterinary authorities to take any samples necessary for national surveillance and research purposes. In addition, by law (http://www.mag.go.cr/legislacion/2006/ley-8495.pdf) SENASA (represented by author BL) has the attribution to do the analysis and to inform the public on the outcome.
Consent for publication
Not applicable. Manuscript does not contain information on individuals.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Temporal patterns of cattle rabies outbreak mortality in Costa Rica (1985–2016). (A) Annual time series of cattle deaths from 1985 to 2016. Peak values occurred in 1985 and 2003, with 149 and 193 cattle deaths, respectively. (B) Bar chart of monthly cumulative cattle deaths from 1985-2016. The highest number of deaths occurred in April and July. (C) Monthly time series of cattle deaths from 1985 to 2016. It reflects both outbreaks during a year and cattle rabies deaths per month. Peak values occurred in 1985 and 2003. Note that in year 1995, although there is no peak value of cattle rabies death, there were many small outbreaks. (D) Boxplots of monthly cattle rabies deaths from 1985 to 2016. For all months from January to December, the median value of the boxplots is 0, which suggests that there is no seasonality in monthly cattle rabies deaths. (PDF 19 kb)
Figure S2. Video showing the geolocation and size of cattle rabies outbreaks (1985-2016). In the video the current outbreak is shown in red while past outbreaks are in blue. Circles are proportional to the number of cattle rabies deaths (for guidance see the inset legend at the bottom of the map) and for every outbreak, the month and year, the exact number of cattle rabies deaths, and the El Niño Southern Oscillation (ENSO) phase are indicated in the upper left corner of the map. (MP4 48445 kb)
Figure S3. Monthly scan cluster analysis. Cluster 1 (circle with dash line in green color) and cluster 2 (circle with dash line in red color) were calculated with data for which denominators were imputed; cluster 3 (circle with dash line in purple color) and cluster 4 (circle with dash line in pink color) were calculated with no imputed denominators. Relative risks (RR) reflect the magnitude of the difference of the deaths in the cluster to regular deaths, which is the expected death assuming they occurred as driven by a homogenous Poisson process. This figure was made using a public domain map from the US National Park Service (https://www.nps.gov/hfc/carto/data-sources.cfm) as background. (PDF 97 kb)
Figure S4. Correlation functions. (A) Auto-correlation function (ACF) of the monthly cattle rabies outbreak number time series; it shows the data were not strongly auto-correlated. (B) Partial auto-correlation function (PACF) of the monthly cattle rabies outbreaks time series data shows that there is a periodicity of 3 months. (C) Cross-correlation function (CCF) of the monthly cattle rabies outbreak number time series data with the ENSO4 index data. This CCF shows no significant linear association over the studied period, when considering up to 22 months of lag. (D) ACF of the monthly outbreak mortality time series. This time series shows low autocorrelation. (E) PACF of the monthly outbreak mortality time series. This time series shows a periodicity of 2 months. (F) CCF of the monthly outbreak mortality time series data with the ENSO4 index data. This CCF shows no significant linear association over the studied period, when considering up to 22 months of lag. In all panels blue dashed lines indicate the 95% confidence limits for correlations that can be expected by random. (PDF 17 kb)
Table S1. Selection of the best monthly cattle rabies outbreaks time series model. Columns indicate the type of model (models): full or the backward elimination round. The Akaike Information Criterion (AIC) is a model selection criterion which is minimized by the best model. The AIC for the best models of each selection round are bolded. o and x indicate, respectively, the presence or absence of a variable in a model. Temp is an abbreviation for temperature. Time lags are in months. (PDF 187 kb)
Table S2. Selection of the best monthly cattle mortality time series model. Columns indicate the type of model: Null, or full model. The Akaike Information Criterion (AIC) is a model selection criterion which is minimized by the best model. The best model has its AIC bolded. o and x indicate, respectively, the presence or absence of a variable in a model. Temp is an abbreviation for temperature. Time lags are in months. (PDF 174 kb)