Controlling disease outbreaks in wildlife using limited culling: modelling classical swine fever incursions in wild pigs in Australia
Veterinary Researchvolume 43, Article number: 3 (2012)
Disease modelling is one approach for providing new insights into wildlife disease epidemiology. This paper describes a spatio-temporal, stochastic, susceptible- exposed-infected-recovered process model that simulates the potential spread of classical swine fever through a documented, large and free living wild pig population following a simulated incursion. The study area (300 000 km2) was in northern Australia. Published data on wild pig ecology from Australia, and international Classical Swine Fever data was used to parameterise the model. Sensitivity analyses revealed that herd density (best estimate 1-3 pigs km-2), daily herd movement distances (best estimate approximately 1 km), probability of infection transmission between herds (best estimate 0.75) and disease related herd mortality (best estimate 42%) were highly influential on epidemic size but that extraordinary movements of pigs and the yearly home range size of a pig herd were not. CSF generally established (98% of simulations) following a single point introduction. CSF spread at approximately 9 km2 per day with low incidence rates (< 2 herds per day) in an epidemic wave along contiguous habitat for several years, before dying out (when the epidemic arrived at the end of a contiguous sub-population or at a low density wild pig area). The low incidence rate indicates that surveillance for wildlife disease epidemics caused by short lived infections will be most efficient when surveillance is based on detection and investigation of clinical events, although this may not always be practical. Epidemics could be contained and eradicated with culling (aerial shooting) or vaccination when these were adequately implemented. It was apparent that the spatial structure, ecology and behaviour of wild populations must be accounted for during disease management in wildlife. An important finding was that it may only be necessary to cull or vaccinate relatively small proportions of a population to successfully contain and eradicate some wildlife disease epidemics.
Wildlife infectious disease can have enormous ecological, biodiversity and societal impacts [1–4]. However, management responses required for mitigation are frequently limited by poor understanding of wildlife disease epidemiology.
Disease modelling is one approach for providing new insights into wildlife disease epidemiology and has allowed important conceptual advances in wildlife disease management . Mathematical modelling was an early method used (and is still widely applied) [6–9]. However, application of this method has often been simplistic, not incorporating many of the major ecological factors that affect disease epidemiology . Furthermore, one of the key concepts in mathematical models - the existence of a threshold level of host abundance required for invasion or persistence of infection - originated in human health and is poorly supported by evidence from wildlife disease studies .
With the improvement of information technology, process models (or simulation models) have been advocated by some as a method of more realistically representing the complexity of real world animal health problems [12, 13]. Process models can capture great complexity, thus enhancing our ability to model complex situations. These models have been widely applied in animal health generally, but relatively less commonly in wildlife disease epidemiology, with some exceptions [14–19].
To take advantage of the great complexity that process models can incorporate, a good understanding of the "process" (host-infection interaction) is required. Sus scrofa, commonly known as wild boar, feral pig, wild hog and wild pig (herein referred to as wild pig) is an important international wildlife species found on every continent except Antarctica . Considerable research has been conducted internationally on wild pig ecology [21–25], and this research can be harnessed to construct detailed process models to study disease epidemiology in this species. Wild pigs have been involved in the transmission or maintenance of many agriculturally important infectious trans-boundary diseases such as African swine fever , pseudorabies  and foot-and-mouth disease , as well as the spread of important endemic zoonoses such as Brucella suis. Classical swine fever (CSF) is another important trans-boundary agricultural disease of domestic and wild pigs . Outbreaks of CSF in Europe have cost many billions of dollars to eradicate , and cause ongoing costs in areas where it is endemic. CSF also has a wide geographic distribution, being found in Asia, Europe, parts of Africa and central and South America , but not Australia. Wild pigs are frequently important in the epidemiology of CSF , but the issue is complex as demonstrated by Boklund et al.  who investigated the potential role of wild boar in CSF epidemics in Denmark.
Some limited epidemiological modelling of CSF in wild pigs has been conducted. Hone and Yip  estimated model parameters with field data and used a mathematical modelling approach to study CSF in wild pig populations. They found that CSF will establish in a small population of wild pigs. Milne et al. , using a process modelling approach found that seasonality is important in dispersal of CSF during epidemics but made some significant logical errors in formulating their model (for example, that wild pigs will be attracted to water only every 4-8 days in the extremely hot Australian sub-tropics). Kramer-Schadt et al.  conducted a review and used a conceptual model to putatively identify the reasons that CSF can persist in some populations. They found virulence of CSF and the size and structure of a wild pig population to be important. They used spatial modelling to show that individual level variation in infection persistence and production of new susceptible individuals was important for disease persistence . Boklund et al.  found a complex epidemiology for CSF outbreaks where wild boar and domestic herds co-exist.
The objective of the research reported in this paper was to enhance knowledge of wildlife disease ecology and assess some control techniques for eradicating disease in wildlife. The paper first describes a spatio-temporal, stochastic, susceptible-exposed-infected-recovered process model that simulates the potential spread of classical swine fever through a well documented, large and free living wild pig population in Australia (which is free of CSF). Results are then used to explore disease ecology and control of CSF in free living wild pigs.
Materials and methods
This study focused on a large wild pig population in a remote area of north-west Australia. It simulated the introduction of a virulent CSF virus into the population to explore epidemic behaviour, disease ecology and various epidemic control options. Three simulation models were developed. Model 1 was a non-spatial within-herd model in which the unit of interest was individual wild pigs (this model is summarised in Additional file 1; results are presented in Table 1). This model was simply used to estimate herd-based epidemiological parameters (i.e. convert individual parameters such as individual infectious period to a herd based parameter). These herd based parameters were then used in a between-herd model (Model 2) which was the main focus of the paper. Model 2 simulated the spread of CSF across a population of wild pig herds in time and space. In Model 2, herds ranged from individual boars to a group of co-mingling wild pigs occupying a territory or home range. The logic of Model 2 was structured on the recommendations of Cowled and Garner  who stated that a number of factors should be accounted for during disease modelling in wild pigs, including distribution and habitat connectivity, density, movements, social and group structure and age structure. Model 3 was a non-spatial herd model designed to replicate model 2, except that model 3 assumed no spatial relationships between wild pig herds. Comparison between model 2 and 3 allowed consideration of the importance of spatial relationships during simulated epidemics. See Additional file 1 for a method summary of model 3.
Study area, biology and distribution of wild pigs
The between-herd model was structured on a population of wild pigs in the Kimberley region of north-west Australia. The Kimberley region is a large (approximately 300 000 km2), remote and sparsely populated pastoral (cattle) region. Pigs were introduced by European settlers during the late 19th century and subsequently became wild . Questionnaire surveys were conducted across the Kimberley region to estimate wild pig distributions and densities and have been previously reported [40, 41]. Wild pigs are currently found across approximately 26 000 km2 of the Kimberley region. The population chosen for disease introduction was located in the Fitzroy River area (see study area and introduction site in Figure 1). Other researchers have investigated the biology of wild pigs in the region  and the population structure was typical of other wild pig populations, for example with groups comprised mostly of solitary boars or herds containing adult females and juveniles [42, 43]. Average group sizes were generally small (mostly herds of 12 or less, but up to 30 pigs). In high density habitat, wild pigs may be present at approximately 3-8 pigs per km2 (the range in density depending on whether an edge effect is taken into account) .
Population at risk and habitat contiguity
Within the known wild pig distribution, permanent water sources (either linear water sources such as rivers or point sources such as dams) were identified and buffered (by 2 km) within a GIS (Mapinfo® v. 10.5). These polygons represented the core habitat of wild pigs in the study area. Thus, a 16 701 km2 area of "core" habitat was identified in the overall distribution of wild pigs (of 26 000 km2). This method was chosen to refine the distribution of wild pigs because the Kimberley region is a tropical ecosystem, is very hot (November mean maximum and minimum daily temperature are 41.0 and 25.6°C respectively ), and in accordance with field observations it is recognised that wild pigs in these conditions require at least daily access to water for survival . Saunders and Kay  demonstrated daily home range lengths of approximately 2 km, and it was therefore assumed that wild pig home ranges must be located within 2 km of permanent water for survival (that is, a large wild pig home range must contain some of this high quality water habitat). Permanent water was identified using a spatial dataset from Geoscience Australia , with additional data layers of artificial water sources (stock water) supplied from the Department of Agriculture and Food, Western Australia (unpublished data).
The total number of wild pigs within each polygon was estimated by classifying the polygon as having a high, medium or low density of wild pigs . Thus the 16 701 km2 area was divided into three polygons of 7563 km2 (low density), 469 km2 (medium density) and 8639 km2 (high density). Relative densities were then quantified using published estimates. Choquenot et al.  reviewed wild pig densities in various habitats in northern Australia. These ranged from 1-20 pigs km-2. Twigg et al.  estimated 3 pigs km-2 within a high density region of this study site. Given our study site was in a semi-arid region it was assumed that densities would be at the lower end of the range listed in Choquenot et al. . Thus, estimates used for population densities were 1, 2 and 3 pigs per km2 in low, medium and high density pig habitat, respectively. Our upper estimate is thus consistent with an edge effect modified estimate from prior work in our study area  and consistent with the lower estimate of Choquenot et al. .
The total population in each polygon was divided into groups. Group sizes and structure were estimated based on both published literature and unpublished data from the study area. Caley  found approximately 12% of a trapped population were males greater than 18 months of age. Thus 12% of the population was assumed to be solitary males. The remaining population was divided into social groups (female groups) with simulation using a B pert distribution of group sizes (minimum group size = 5, most likely = 7, maximum = 45) based on prior research on group sizes [24, 42] and allowing for a slightly greater range due to the limited size of the study in Twigg et al. . These groups were dispersed randomly in the core habitat (2 km buffered permanent water) whilst maintaining low, medium and high relative density classifications. This formed the final population data base for use in Model 2. The simulated distribution of wild pig herds (and solitary males) throughout the Kimberley region is shown in Figure 1.
Classical swine fever
There are several reports detailing mortality and morbidity rates associated with CSF outbreaks in wild pigs. There appears to be a wide spectrum of clinical outcomes seen, with some outbreaks in wild pigs leading to very high mortality and morbidity, [34, 48–50] with other caused by low or moderate virulence strains, especially in Europe [51–54].
South-east Asia has regions that are endemically infected with CSF but is also in relatively close proximity to our Australian study site. For geographical reasons south-east Asia may thus represent a potential source of an outbreak of CSF in wild pigs at our study site. The virulence of South-east Asian CSF strains are largely undocumented, but outbreaks in some islands with highly susceptible pig populations have lead to substantial mortality events and may therefore be due to moderate or highly virulent strains (Jenny-Ann Toribio, personal communication, July 2011).
Therefore, in this study, it was assumed that the virus was highly virulent, although a sensitivity scenario assuming a lower virulence was also conducted. In the highly virulent simulations, an individual case fatality rate of 90% associated with highly virulent CSF infection was assumed (within-herd model-see Additional file 1). This resulted in 42% of herds having all members die. Table 1 summarises the key epidemiological parameter estimates used in Model 2.
Model 2: description of between-herd model
The model was reminiscent of a previous disease model in domestic animal populations [55–57]. Applications were coded in MapBasic®, and implemented in Mapinfo®. These software environments together represent a sophisticated and customisable geographical information system (GIS).
Classical swine fever transmission
The model's treatment of virus transmission can be considered in two ways, transition of individual herds between disease states temporally, and between herd transmission.
When a susceptible herd becomes exposed to virus it may become infected whereby it will progress through a latent, infectious and recovered sequence, although the herd may cease to exist if all members are killed by the infection (in the model this is simulated according to a probability derived from Model 1). The time spent in each state is stochastically determined using probability distributions. The time step used in the model is one day. To represent disease transmission between herds it is necessary to consider the spatial distribution of the wild pig population (see sections below). Local herds have the chance of coming into contact where their daily home ranges overlap, and if one of these herds is infected then there is a probability of virus transmission between an infected and un-infected herd. In the absence of good field information, this probability was arbitrarily set at 0.75, with a thorough sensitivity analysis undertaken to determine how the model outcomes change in response to changes in this parameter. Figure 2 shows diagrammatically how disease transmission between a single infected herd and neighbouring susceptible herd occurs. In this way, CSF can spread through a wild pig population comprising contiguous herds. However, CSF (particularly virulent strains) frequently produces severe clinical disease  that could be expected to affect the activity of feral pigs . To allow for this, when wild pig herds are in the infected state, their mobility is assumed to decline to just 10% of their normal daily movement and home range (but sensitivity analyses occurred) (see Figure 2). This limits the probability of overlapping home ranges and thereby reduces the chance of transmitting virus to nearby uninfected herds.
Model logic (movements and home ranges)
Most ecological studies in Australia have demonstrated that wild pigs are relatively sedentary, within fixed home ranges, displaying little or no dispersal but moving small distances daily within their larger home range [24, 61]. Females move smaller distances than males. However, a small but potentially epidemiologically important proportion of pigs may disperse longer distances [61–63].
To capture movements, each pig herd (including solitary males) was assigned an annual home range. This was simply a circle (buffer) around each location. There was considerable overlap in annual home ranges, reflecting the overlap that occurs in the field . Pig herds were randomly moved a linear distance each day within each home range. Ninety-five percent of pig herds were constrained to movements within their own home range. It was assumed that since pig herds usually consist of females and offspring then published daily movements of females will represent the movements of groups, whilst the published daily movements of males will represent the movements of solitary pigs. Daily movement distances were estimated from the daily home ranges described by Caley . After their daily movement, a daily home range of approximately 1 km2 was structured around the herd's final location for the day [45, 47]. See Table 2 for a description of ecological parameters used in Model 2.
A small percentage (5%) of pig herds were allowed to move a normal daily distance but were unrestrained by their home range and were able to move to adjacent core water habitat. In effect, this allowed a small proportion of pigs to disperse or display extraordinary movements. Giles  showed that groups of wild pigs can move 20-30 km during short periods in response to natural events such as flooding. Saunders and Bryant  showed that wild pigs can move from a study area in response to persecution, and that long range dispersal can occur. Caley  showed that a small proportion of wild pigs may move 20-30 km over several years (although the majority stay within their home ranges).
Model logic (surveillance)
A surveillance module was included to allow assessment of the surveillance of wild pigs for disease after an outbreak was discovered. The aim of surveillance was to delineate the infected area of the pig population following detection of the disease incursion . Surveillance was simulated using realistic surveillance strategies and assumptions, and integrated with the control modules (see below).
A time to first detection of the outbreak was selected. At that point an index case was randomly chosen to be found, from all the infectious or recovered herds present at that time. A six week time to detection was arbitrarily chosen for simulations -- this is similar to other published estimates of potential time to first detection of FMD in Australian wild pigs  and detection of The Netherlands domestic outbreak . Surveillance was assumed to begin three days after detection of the index case to allow organisation of surveillance resources. Surveillance was then structured around a user defined grid (a 10 × 10 km grid structured across the wild pig distribution in the Kimberley region was used). Surveillance was assumed to be conducted by aerial shooting from helicopters. Aerial shooting from helicopters is a well researched, effective and humane wild pig control and surveillance tool used in Australia . The number of helicopters to be used (3), how many individual pigs that could be sampled by one helicopter team each day (70) and the area a helicopter can search each day (200 km2) were selected (based on author experience). Assuming 4-5 pigs are selected from each sampled herd (to give 95% confidence of detecting disease, where prevalence is assumed 50% ), the number of pig herds that can be sampled each day was determined. A sampling intensity was thus calculated, based on the number of herds that can be sampled each day and the average population of herds within a grid. Whether any given herd within a grid cell was actually sampled was probabilistically determined from the sampling intensity using Monte Carlo methods. A sampled herd was probabilistically categorised as infected based on defined test sensitivity (95%).
A circle two grid cells in radius surrounding an index cell (i.e. the index cell is the grid cell that contains the index case) was buffered and grid cells within this circle selected for surveillance. These cells were progressively sampled from closest to the index cell to furthest, each day depending on resources available. When all cells were sampled within the initial buffered region, and disease was present in at least one cell, the search area was expanded by another two grid cells and all grid cells were again sampled. This progressed until a final buffered area underwent surveillance and no infected cells were discovered. At this point an assumption was made that the epidemic was delineated.
Model logic (control)
Two control strategies were implemented in the model (aerial culling or vaccination) although only one of these methods could be selected during a single epidemic simulation. Control was assumed to begin after surveillance had finished delineating the infected area. For each control strategy, the infected area delineated during surveillance was buffered. The buffered area of land surrounding the infected area (herein control zone) was thus at least several pig home ranges wider than the "known" infected area. The control zone was constructed in this way to prevent migration of incubating or infected pig herds outside the infected area, and hence prevent spread of the epidemic to neighbouring susceptible populations.
To implement this, all grid cells within the control zone were ordered from the centre of the infected area outwards. If culling was the chosen control method, for each simulated day a portion of herds (priority from closest to the centre of the infected area to furthest) was culled. The proportion of the herds culled was defined by both the availability of control resources and the probability that a herd would be detected during aerial culling operations. The availability of control resources was measured by two parameters, the number of helicopters used for culling (4) and how many individual pigs could be culled each day (300) by each helicopter. The probability that a herd would be detected and culled during culling operations was user defined, with a default of 0.8 , meaning that 80% of herds would be randomly culled during baseline model simulations. The effect of culling a range of proportions of the population was assessed during experimentation and sensitivity analysis.
For the vaccination control option it was assumed an oral CSF vaccine that could be distributed aerially was available in Australia. In contrast to the culling option, pig herds were prioritised from furthest to closest to the zone centre. This approach was assumed to ensure that incubating pig herds had less opportunity to migrate infection beyond the immune buffer, before immunity developed. Similar to culling, there was a probability that each herd would be vaccinated, and the time to vaccinate herds in the control zone was determined by control resources available. A delay of 7-14 days until full herd vaccine immunity develops following vaccination was assumed [69, 70]. Once immune, it was assumed that herds could not transmit virus.
During each simulation in which control was instigated, epidemics were classified as successfully or unsuccessfully contained and eradicated. In each control simulation the total infected area of a controlled epidemic was calculated and compared with the identical but uncontrolled scenario. If the final epidemic area during a controlled scenario was less than the same epidemic without control, and if eradication occurred, it was assumed that control measures had contained and eradicated an epidemic.
Number of simulations
A vexed question for simulation modelling is how many simulations (or model runs) are required to produce a result of sufficient precision? Too many simulations are computationally inefficient. However, outputs from a stochastic simulation model have variability. If each simulation is considered one observation in a sample, it is important to have enough simulations (or a large enough sample size) to ensure that the estimate of the parameter of interest () approaches the true population value for the model (θ). If the sample size is large enough, the parameter estimate () generated from the model simulations converges with the true population value (θ)  for the model. Consistency (or convergence) can therefore be stated in relation to how the variance of the sample reduces to zero as the sample size approaches a theoretical infinity :
Thus in practical terms, when a sample size increases such that the variance is minimised, is close to the true θ. To estimate our sample size, we calculated the mean of the parameter-of-interest (after each simulation). We then determined the co-efficient of variation of this mean. At the point when the co-efficient of variation was less than 15% for 30 consecutive simulations we considered that convergence had occurred and that this number of simulations was adequate to estimate the parameter with precision. We repeated this process for every output parameter of the simulation model, and determined the maximum number of model simulations required across all output parameters. This number became our sample size (the number of simulations required).
Sensitivity analyses and detection of interaction
Best estimates (as assumed following the literature review and detailed above) for all input parameters were assessed during baseline runs. For all baseline runs, sensitivity and experimental analyses, infection was introduced into the same wild pig herd (see Figure 1) to ensure a valid comparison of outputs. The major ecological, epidemiological and population parameters were varied systematically, by multiplying the best estimates one at a time by 0.25, 0.5, 0.75, 1 (best estimate) 1.5 and 2. An exception was made for transmission probability, in which the 1.5- and 2-times factors were excluded and 1.33-times (probability = 0.99) included to ensure the probability remained less than one. Parameters selected for sensitivity analyses were:
CSF transmission probability (between herds with overlapping home ranges)
Herd mortality rate (proportion of herds with all individuals dying of CSF)
Home range size
Daily linear movement distances
Proportion of population that can move extra-ordinary distances
Density (pigs km-2)
Reduction in movement of a clinically affected herd (proportion)
Outputs measured are listed in Table 3. All output measures underwent pair-wise linear regression against the area of the infected land to determine whether parameter outputs were correlated and whether a single output measure could be used for comparison during sensitivity analyses. Subsequently, the total area infected was used as the output parameter to represent the scale of the epidemic.
Sensitivity analysis used an iterative process : first input parameters were screened to identify influential and non-influential parameters, and then influential parameters for interaction were identified. Scatter plots and regression (linear and polynomial) were used to identify influential parameters. The most influential parameters identified during one at a time sensitivity analyses underwent factorial experiments and were tested for interactions by ANOVA. A full factorial experiment was conducted in which the total area of the epidemic (km2) was the response variable and explanatory variables identified in one at a time analyses were varied at three different levels (half, baseline and double parameter estimates). The factorial experiments occurred during simulations in which disease was assumed detected at 6 weeks post-introduction and control using aerial culling was instigated after delineation of infected areas (see culling scenarios for more information). ANOVA included testing for effects and for pair-wise and three way interactions.
Scenarios analysed and model experimentation
Experimentation was conducted to examine the effects of disease ecology, surveillance and control options on the scale of an outbreak. Table 4 details the scenarios that were simulated.
Sample sizes required to achieve a low coefficient of variation for 30 consecutive simulations were calculated for all outputs. The maximum number of simulations required was 59 (for the output measure, total infected area). See Figure 3.
Simulation results using best parameter estimates (and a highly virulent virus)
Using the best estimates for parameters, CSF generally established (98% of simulations) following a single point introduction. CSF generally progressed in an epidemic wave down contiguous habitat for several years, before dying out (when the epidemic arrived at the end of a contiguous sub-population or at a low density pig area). The daily herd incidence rate was low, despite epidemics that lasted several years, across thousands of square km and cumulatively infecting thousands of herds.
Table 3 lists and defines outputs and the results of the simulations. Figure 4 is a typical epidemic curve for one simulation. Additional file 2 is a PowerPoint presentation that shows a week by week progression of a typical simulated epidemic (Model 2).
In general, output parameters were highly correlated with one another (r > 0.94), except for the number of immune animals and incidence rate (r < 0.1). Subsequently, the total area infected was arbitrarily chosen as the single output variable for comparison of parameter estimates during sensitivity analyses.
Four parameters appeared influential (Figure 5 and Table 5). As the density of herds, the daily linear distance that a herd could move and probability of disease transmission between herds with overlapping daily home ranges increased, epidemics were larger. However, as the probability that all members of a herd would die from CSF infection (a proxy for CSF virulence) increased, epidemics became smaller. Parameters that had little predictable effect on epidemic size were the proportion of the herds that moved extraordinary distances, the yearly home range size of a pig herd and the reduction in the movements of a herd when it became clinically affected.
Four influential factors were identified during one at a time sensitivity analyses. Since the density of herds and the daily linear distance a herd can move were highly correlated and changes in these variables produced near identical effects (see Figure 5), and to reduce the number of scenarios investigated during the factorial experiments the density of herds was not included during factorial experimentation. Thus, three influential variables (the daily linear distance a herd can move, probability of a herd dying due to CSF and probability of transmission of infection between overlapping herds) were varied at three levels resulting in 27 permutations for inclusion in analyses of the factorial experiment. The response variable (area of the epidemic) was not normally distributed and was log transformed, with the resulting transformed distribution being approximately normal (skewness = 0.298, kurtosis = 2.48).
In ANOVA all three main effects were significant and were significantly different from each other (P < 0.05) but no interactions were detected (P > 0.05).
Scenarios and experimental results
Based on the assumptions made, the surveillance program was able to consistently delineate the infected area within 1-2 weeks.
Culling was most successful at controlling the incursion when high proportions of herds were culled and control zones sizes were relatively large (see Table 6 and Figure 6), or when at least one variable was high (either a very wide culling zone or very high proportion of herds culled). Conversely, where culling proportions were low and control zone sizes small, containment and eradication attempts were less successful. If a realistic proportion of the herds could be culled - for example, > 60% of herds  - containment and eradication was achieved in all simulations in which the culling zone was sufficiently large (> 30 km width). Table 6 shows that there were critical combinations of culling proportion and control zone size above which disease could be contained and below which disease escaped in some simulations.
Similarly to culling, vaccination was most successful in scenarios in which a high proportion of herds were able to be vaccinated and vaccination zone widths were large (see Table 7 and Figure 6). Comparison of Table 6 with Table 7 shows similarities between the success of culling and vaccination. However, it was also evident that vaccination was less effective at containing and eradicating epidemics than culling, with a greater number of combinations of buffer sizes and widths resulting in some outbreaks unable to be contained. For example, assuming 99% cull at a buffer width of 10 km culling was successful but a 99% vaccination proportion with a 10 km vaccination zone was not. This is supported by Figure 6 in which the size of an epidemic was always smaller if culling was used rather than vaccination, even if containment and eradication was still achieved.
Low virulence strain
Outbreaks of a lower virulence CSF were longer in duration, spreading across a larger land area and infecting more herds (with a higher cumulative incidence). The lower virulence strain resulted in lower herd mortality but a similar incidence rate (see Table 4).
Under non-spatial modelling assumptions the epidemic progressed much more quickly than during spatial modelling (Table 3) - incidence rates of 61 versus approximately 2 herds/day, respectively. Outbreaks were of shorter duration (median 87 days versus several years, respectively), but with all herds being infected (100% versus 33% cumulative incidence, respectively). Control (assuming an 80% probability of herds being culled) was able to shorten the duration of epidemics (70 days until disease fadeout) but this required that a median of 3857 herds (73% of herds) be culled for eradication to occur (compared with 419 herds or 11% of herds required to be culled for disease eradication during spatial modelling).
This study highlights the critical importance of ecology, behaviour and spatial structuring of wildlife populations on the management of wildlife disease. It also provides some useful guidelines for containing and eradicating epidemics from widespread terrestrial vertebrate species such as the wild pig.
The spatial structuring evident in the wild pig population in our study area  had a large influence on simulated epidemics. For example, in our simulations, epidemics typically travelled in waves along the larger river systems along which wild pig meta-populations were located (see Additional file 2). This spatial population structuring resulted in a very low daily incidence rate (< 2 herds each day). In contrast, epidemics in non-spatial populations progressed much more quickly (~61 herds per day). This occurred because epidemics were not confined to that proportion of the population in the immediate vicinity of infected herds during non-spatial simulations, but instead could infect herds randomly throughout the entire population. That is, models that do not take realistic spatial structures into account may overestimate the rate at which a disease will spread and overestimate the size of an outbreak.
The spatial nature of epidemics is of critical importance when planning surveillance for wildlife disease. Despite typical epidemics infecting thousands of wild pig herds, across thousands of square kilometres and lasting several years, the spatial structure of the population meant the number of currently infected herds was just a handful each day and new infections were limited to those herds directly adjacent to the narrow epidemic front. This highlights potential difficulties in attempting to detect active infection in wild populations using agent identification techniques (such as culture or PCR) in which no carrier state exists and immunity to an infectious agent is long lived (e.g. FMD and CSF in wild pigs). This is in accordance with prior research . In addition, given the rapid turnover of individuals in many wild populations (especially wild pigs ) the decline in herd immunity may be rapid, leaving a relatively narrow window to conduct even serological investigations. Taken together this suggests that surveillance based on detection and investigation of clinical events is likely to be an efficient approach to finding short lived infections. Alternatively, in situations in which passive detection and reporting cannot be relied upon (e.g. little or limited opportunity to observe the population) then a more structured surveillance approach would benefit from a very good understanding of risk to allow targeted and therefore efficient sampling.
These simulation results also suggest that the spatial structuring and behaviour of wild animal populations should have an influence on the design of containment and eradication programs. Wild pigs are relatively sedentary with a high fidelity to a home range ; only short daily movements are observed [45, 47]. Populations in Australia are also predictably associated with riverine habitat [40, 75]. Our findings indicate that epidemics would only spread relatively slowly across the landscape (~9 km2 each day) and are containable with some relatively simple and well researched control methods such as aerial shooting . An adequate design in our modelling was a control zone width of approximately 30 km around the infected area, in which 60% of herds could be culled (although other combinations of control zone width and proportions able to be culled were also adequate). Indeed, following simulated introductions of disease into the largest and highest density sub-population of wild pigs in the Kimberley region (Fitzroy River populations), disease could be eradicated by culling just a median of 419 herds per outbreak (representing approximately 11% of herds in the sub-population of interest - the Fitzroy River population).
In comparison, when culling was instigated in the non-spatial model, very large numbers of wild pig herds (median 3857 herds) were required to be culled for disease fadeout to occur. This was because culling was not targeted to the immediate vicinity of the epidemic but was randomly implemented across the wild pig population. This is likely the reason that our findings differ from those of earlier mathematical modelling studies in wild pigs indicating that very high culling levels (e.g. 95%) may be required over short time periods or that relatively high culling rates may be required for a number of years (e.g. 49% per annum) to eradicate a disease such as FMD in wild pigs . However, our modelling results are in agreement with other simulation modelling studies informed by empirical field trials in other species. For example, in the Ethiopian wolf only small vaccination corridors were required to reduce disease transmission in spatially structured wildlife populations .
Disease control using vaccination was generally less effective than aerial culling in our simulations. Specifically, in every culling or vaccination combination (width of control zone and proportion of herds culled or vaccinated), culling lead to smaller epidemics than vaccination (see Figure 6) (although slight differences in program implementation should be noted: namely vaccination occurred from the outside in, culling from the inside out). Additionally, infection was routinely contained at lower combinations of control zone width and probability in the culling scenarios than the vaccination scenarios (Table 6 c.f. Table 7). The reasons for this could include the time taken for immunisation to become effective, resulting in pig herds continuing to become infected and continuing to transmit disease during vaccination programs compared to culling (in which a proportion of pig herds are immediately removed from the population). The effectiveness of culling provides countries like Australia, where wild pigs are a damaging, introduced invasive species and where control is mandated, a great advantage for controlling disease in wild pig populations. In other continents such as Europe where wild boar are a valued endemic species and effective aerial culling may not be as acceptable, there is greater reliance on less effective tools such as vaccination. This may be equally true during outbreaks in domestic pigs, in which ethical concerns force consideration of vaccination over culling .
Other authors have reviewed the factors contributing to persistence of CSF in wild boar populations . They found that attenuation of CSF viruses to moderate virulence, as well as the size and structure of wild boar populations may affect persistence, and recommended further spatially explicit modelling. Our modelling supports their hypothesis, with epidemics induced under assumptions of lower virulence resulting in longer and larger epidemics with higher cumulative incidences. However, in both low virulence and high virulence scenarios, without further introductions infection always died out after several years because as an epidemic front reached the end of a contiguous wild pig population there were insufficient susceptible hosts to maintain the infection. Infection could not spread back along the previous route of the epidemic because these herds had either died from CSF or were immune. These results concur with the majority of published studies investigating the nature of CSF outbreaks in European wild boar populations, in which CSF outbreaks generally fade out after several years [32, 50].
This research may assist in the elucidation or confirmation of some important factors that influence the size, scale and behaviour of epidemics. Our results suggest density, the daily linear distance a herd can move, the probability of herd death from CSF and the probability of transmission between herds can have a large influence on the size and scale of an epidemic. Addressing each of these in turn, it appears that density and daily movements are highly correlated. This indicates that at higher densities, or with greater daily movements, larger epidemics occur because infectious herds are more likely to come into contact with other susceptible herds. Arbitrarily increasing the probability of transmission also increased the size of epidemics for similar reasons. In contrast to this, as the probability that a herd will die following infection increased, epidemics declined in size. This was likely because herds were eliminated before they had a chance to transmit infection to nearby herds. However, having a small proportion of the population that are unrestrained by assumptions of home range fidelity or varying the size of a herd's annual home range made little difference to the overall scale of the epidemic. Several authors have found that Australian wild pigs have a high home range fidelity [61, 62, 78], but have also found that a small proportion of pigs can move much greater distances and are unrestrained to home ranges [61, 62]. Based on these modelling results it appears that these individuals would have little influence on disease spread, potentially because the probability that one of these individuals is incubating disease is low, and because clinical impacts of CSF reduce long range, aberrant movements anyway.
As with any simulation modelling, our results are very much dependant on the assumptions made. One of our main assumptions was that there are no artificial, human derived movements of infected wild pigs or contaminated fomites in the study area. Given that the Kimberley region is one of the least densely human populated regions in the world, and that there is no commercial pig production, we believe this assumption to be valid. However, were human mediated movements found to be important for transmitting infection, it is likely that control and surveillance programs would be made considerably more difficult. Although there is no field or published evidence of wild pigs chronically infected with CSF , Kramer-Schadt et al.  supposed that this may have been missed and assumed chronically infected pigs for the purposes of modelling. They used modelling to demonstrate that the existence of chronically infected wild boar may be a plausible mechanism to explain persistence of infection in a region. Other evidence  suggests that persistently infected piglets are critical in sustaining infection. Our modelling assumed chronically infected wild pigs, but only allowed for the longest recorded wild boar piglet infectious periods of 39 days,  which is shorter than has been recorded in domestic pigs .
In conclusion, our modelling has captured many of the important factors that are likely to influence epidemic behaviour in wild pig populations. Our results indicate that spatial structuring of wild pig populations is an extremely important feature. Density, daily movement distances, disease-induced herd mortality rates and transmission probabilities between adjacent herds are also important. It is also evident that control and surveillance programs should account for the spatial structuring of wild populations, and that it may only be necessary to cull or vaccinate relatively small proportions of a population to successfully contain and eradicate wildlife disease epidemics.
Jones KE, Patel NG, Levy MA, Storeygard A, Balk D, Gittleman JL, Daszak P: Global trends in emerging infectious diseases. Nature. 2008, 451: 990-993. 10.1038/nature06536.
Normile D: Driven to extinction. Science. 2008, 319: 1606-1609. 10.1126/science.319.5870.1606.
McCallum H, Jones M, Hawkins C, Hamede R, Lachish S, Sinn DL, Beeton N, Lazenby B: Transmission dynamics of Tasmanian devil facial tumor disease may lead to disease-induced extinction. Ecology. 2009, 90: 3379-3392. 10.1890/08-1763.1.
Li W, Shi Z, Yu M, Ren W, Smith C, Epstein JH, Wang H, Crameri G, Hu Z, Zhang H, Zhang J, McEachern J, Field H, Daszak P, Eaton BT, Zhang S, Wang LF: Bats are natural reservoirs of SARS-like coronaviruses. Science. 2005, 310: 676-679. 10.1126/science.1118391.
Anderson RM, May RM: Population biology of infectious diseases: Part I. Nature. 1979, 280: 361-367. 10.1038/280361a0.
Murray AG: A simple model to assess selection for treatment-resistant sea lice. Ecol Model. 2011, 222: 1854-1862. 10.1016/j.ecolmodel.2011.03.016.
Anderson RM, Jackson HC, May RM, Smith AM: Population dynamics of fox rabies in Europe. Nature. 1981, 289: 765-771. 10.1038/289765a0.
Pech RP, Hone J: A model of the dynamics and control of an outbreak of foot and mouth disease in feral pigs in Australia. J Appl Ecol. 1988, 25: 63-78. 10.2307/2403610.
Brauer F: Mathematical epidemiology is not an oxymoron. BMC Public Health. 2009, 9 (Suppl 1): S2-
Barlow ND: The ecology of wildlife disease control: simple models revisited. J Appl Ecol. 1996, 33: 303-314. 10.2307/2404752.
Lloyd-Smith JO, Cross PC, Briggs CJ, Daugherty M, Getz WM, Latto J, Sanchez MS, Smith AB, Swei A: Should we expect population thresholds for wildlife disease?. Trends Ecol Evol. 2005, 20: 511-519. 10.1016/j.tree.2005.07.004.
Hurd SH, Kaneene JB: The application of simulation models and systems analysis in epidemiology: a review. Prev Vet Med. 1993, 15: 81-99. 10.1016/0167-5877(93)90105-3.
Thulke HH: Application of recent approaches in modelling for Animal Health. Prev Vet Med. 2011, 99: 1-3. 10.1016/j.prevetmed.2011.01.007.
Ward MP, Laffan SW, Highfield LD: Modelling spread of foot-and-mouth disease in wild white-tailed deer and feral pig populations using a geographic-automata model and animal distributions. Prev Vet Med. 2009, 91: 55-63. 10.1016/j.prevetmed.2009.05.005.
Ward MP, Laffan SW, Highfield LD: The potential role of wild and feral animals as reservoirs of foot-and-mouth disease. Prev Vet Med. 2007, 80: 9-23. 10.1016/j.prevetmed.2007.01.009.
Highfield LD, Ward MP, Laffan SW, Norby B, Wagner GG: The impact of potential mitigation strategies on the predicted spread of foot and mouth disease in white-tailed deer in south Texas. Prev Vet Med. 2010, 94: 282-288. 10.1016/j.prevetmed.2010.01.015.
Thulke HH, Tischendorf L, Staubach C, Selhorst T, Jeltsch F, Müller T, Schlüter H, Wissel C: The spatio-temporal dynamics of a post-vaccination resurgence of rabies in foxes and emergency vaccination planning. Prev Vet Med. 2000, 47: 1-21. 10.1016/S0167-5877(00)00167-7.
Selhorst T, Thulke HH, Müller T: Cost-efficient vaccination of foxes (Vulpes vulpes) against rabies and the need for a new baiting strategy. Prev Vet Med. 2001, 51: 95-109. 10.1016/S0167-5877(01)00209-4.
Doran RJ, Laffan SW: Simulating the spatial dynamics of foot and mouth disease outbreaks in feral pigs and livestock in Queensland, Australia, using a susceptible-infected-recovered cellular automata model. Prev Vet Med. 2005, 70: 133-152. 10.1016/j.prevetmed.2005.03.002.
Sus scrofa. [http://www.iucnredlist.org]
Mayer JJ, Brisbin IL: Wild Pigs of the United States: Their History, Morphology, and Current Status. 1991, Athens: University of Georgia
Li CQ, Chang Q, Chen JQ, Zhang BW, Zhu LF, Zhou KY: Population structure and phylogeography of the wild boar Sus scrofa in Northeast Asia based on mitochondrial DNA control region variation analysis. Acta Zool Sinica. 2005, 51: 640-649.
Genov PW, Massei G, Kostova W: The utilization of wild boar (Sus scrofa L.) in Europe in theory and practice. Zeitschrift fur Jagdwissenschaft. 1994, 40: 263-267. 10.1007/BF02243448.
Choquenot D, McIlroy J, Korn T: Managing Vertebrate Pests: Feral Pigs. 1996, Bureau of Resource Sciences, Australian Government Publishing Service, Canberra
Gongora J, Fleming P, Spencer PBS, Mason R, Garkavenko O, Meyer JN, Droegemueller C, Lee JH, Moran C: Phylogenetic relationships of Australian and New Zealand feral pigs assessed by mitochondrial control region sequence and nuclear GPIP genotype. Mol Phylogenet Evol. 2004, 33: 339-348. 10.1016/j.ympev.2004.06.004.
Costard S, Wieland B, de Glanville W, Jori F, Rowlands R, Vosloo W, Roger F, Pfeiffer DU, Dixon LK: African swine fever: how can global spread be prevented?. Philos Trans R Soc London B Biol Sci. 2009, 364: 2683-2696. 10.1098/rstb.2009.0098.
Corn JL, Stallknecht DE, Mechlin NM, Luttrell MP, Fischer JR: Persistence of pseudorabies virus in feral swine populations. J Wildl Dis. 2004, 40: 307-310.
Khomenko S: Silence of the wild boar. Book Silence of the wild boar. 2011, 559-619.
Van Der Leek ML, Becker HN, Humphrey P, Adams CL, Belden RC, Frankenberger WB, Nicoletti PL: Prevalence of Brucella sp. antibodies in feral swine in Florida. J Wild Dis. 1993, 29: 410-415.
OIE: Classical Swine Fever: General Disease Information Sheets. Book Classical Swine Fever: General Disease Information Sheets. 2011, World Organisation for Animal Health
Meuwissen MPM, Horst SH, Huirne RBM, Dijkhuizen AA: A model to estimate the financial consequences of classical swine fever outbreaks: principles and outcomes. Prev Vet Med. 1999, 42: 249-270. 10.1016/S0167-5877(99)00079-3.
Laddomada A: Incidence and control of CSF in wild boar in Europe. Vet Microbiol. 2000, 73: 121-130. 10.1016/S0378-1135(00)00139-5.
Boklund A, Goldbach SG, Uttenthal A, Alban L: Simulating the spread of classical swine fever virus between a hypothetical wild-boar population and domestic pig herds in Denmark. Prev Vet Med. 2008, 85: 187-206. 10.1016/j.prevetmed.2008.01.012.
Hone J, Pech R, Yip P: Estimation of the dynamics and rate of transmission of classical swine fever (hog cholera) in wild pigs. Epidemiol Infect. 1992, 108: 377-386. 10.1017/S0950268800049840.
Milne G, Fermanis C, Johnston P: A mobility model for classical swine fever in feral pig populations. Vet Res. 2008, 39: 53-10.1051/vetres:2008029.
Kramer-Schadt S, Fernandez N, Thulke HH: Potential ecological and epidemiological factors affecting the persistence of classical swine fever in wild boar Sus scrofa populations. Mammal Rev. 2007, 37: 1-20. 10.1111/j.1365-2907.2007.00097.x.
Kramer-Schadt S, Fernandez N, Eisinger D, Grimm V, Thulke HH: Individual variations in infectiousness explain long-term disease persistence in wildlife populations. Oikos. 2009, 118: 199-208. 10.1111/j.1600-0706.2008.16582.x.
Cowled B, Garner G: A review of geospatial and ecological factors affecting disease spread in wild pigs: considerations for models of foot-and-mouth disease spread. Prev Vet Med. 2008, 87: 197-212. 10.1016/j.prevetmed.2008.03.012.
Pullar EM: The wild (feral) pigs of Australia: their origin, distribution and economic importance. Book The wild (feral) pigs of Australia: their origin, distribution and economic importance. Volume 18. 1953, National Museum, 7-23.
Cowled BD, Giannini F, Beckett SD, Woolnough A, Barry S, Randall L, Garner G: Feral pigs: predicting future distributions. Wildlife Res. 2009, 36: 242-251. 10.1071/WR08115.
Woolnough AP, West PB, Saunders GR: Institutional knowledge as a tool for pest animal management. Ecol Manage Restor. 2004, 5: 226-228.
Twigg LE, Lowe T, Martin G, Everett M: Feral pigs in north-western Australia: basic biology, bait consumption, and the efficacy of 1080 baits. Wildlife Res. 2005, 32: 281-296. 10.1071/WR04110.
Graves HB: Behavior and Ecology of Wild and Feral Swine (Sus Scrofa). J Anim Sci. 1984, 58: 482-492.
Climate Data Online. [http://www.bom.gov.au/climate/data/]
Saunders G, Kay B: Movements of feral pigs Sus-scrofa at Sunny Corner New South Wales Australia. Wildlife Res. 1991, 18: 49-62. 10.1071/WR9910049.
Geosience_Australia: GEODATA TOPO 250K Series 2 Topographic Data. Book GEODATA TOPO 250K Series 2 Topographic Data. 2006, Geoscience Australia
Caley P: The ecology and management of feral pigs in the "wet-dry" tropics of the Northern Territory. 1993, University of Canberra
Immediate notification: Classical Swine Fever, Israel. [http://web.oie.int/wahis/public.php?page=event_summary&reportid=7861]
Nettles VF, Corn JL, Erickson GA, Jessup DA: A survey of wild swine in the usa for evidence of hog cholera. J Wild Dis. 1989, 25: 61-65.
Artois M, Depner KR, Guberti V, Hars J, Rossi S, Rutili D: Classical swine fever (hog cholera) in wild boar in Europe. Rev Sci Tech. 2002, 21: 287-303.
Immediate notification: Classical Swine Fever Bulgaria. [http://web.oie.int/wahis/public.php?page=single_report&pop=1&reportid=8436]
Kaden V, Ziegler U, Lange E, Dedek J: Classical swine fever virus: clinical, virological, serological and hematological findings after infection of domestic pigs and wild boars with the field isolate "Spante" originating from wild boar. Berl Munch Tierarztl Wochenschr. 2000, 113: 412-416.
Kaden V, Lange E, Polster U, Klopfleisch R, Teifke JP: Studies on the virulence of two field isolates of the classical Swine Fever virus genotype 2.3 rostock in wild boars of different age groups. J Vet Med B Infect Dis Vet Public Health. 2004, 51: 202-208. 10.1111/j.1439-0450.2004.00759.x.
Immediate notification: Classical Swine Fever Russia. [http://web.oie.int/wahis/public.php?page=single_report&pop=1&reportid=8951]
Garner MG, Beckett SD: Modelling the spread of foot-and-mouth disease in Australia. Aust Vet J. 2005, 83: 758-766. 10.1111/j.1751-0813.2005.tb11589.x.
Dubé C, Stevenson MA, Garner MG, Sanson RL, Corso BA, Harvey N, Griffin J, Wilesmith JW, Estrada C: A comparison of predictions made by three simulation models of foot-and-mouth disease. N Z Vet J. 2007, 55: 280-288. 10.1080/00480169.2007.36782.
Beckett SD, Garner MG: Simulating disease spread within a geographic information system environment. Vet Ital. 2007, 43: 595-604.
Mapping & Analytical. [http://www.pbinsight.com/products/location-intelligence/applications/mapping-analytical]
Moennig V, Floegel-Niesmann G, Greiser-Wilke I: Clinical signs and epidemiology of classical swine fever: a review of new knowledge. Vet J. 2003, 165: 11-20. 10.1016/S1090-0233(02)00112-0.
Dahle J, Liess B: A review on classical swine fever infections in pigs epizootiology clinical disease and pathology. Comp Immunol Microbiol Infect Dis. 1992, 15: 203-211. 10.1016/0147-9571(92)90093-7.
Caley P: Movements, activity patterns and habitat use of feral pigs (Sus scrofa) in a tropical habitat. Wild Res. 1997, 24: 77-87. 10.1071/WR94075.
Giles JR: The ecology of the feral pig in western New South Wales. PhD thesis. 1980, University of Sydney
Spencer PBS, Lapidge SJ, Hampton JO, Pluske JR: The sociogenetic structure of a controlled feral pig population. Wild Res. 2005, 32: 297-304. 10.1071/WR04076.
Saunders G, Bryant H: The evaluation of a feral pig eradication program during a simulated exotic disease outbreak. Aust Wild Res. 1988, 15: 73-82. 10.1071/WR9880073.
Australia AH: Wild Animal Response Strategy (Version 3.2). Australian Veterinary Emergency Plan (AUSVETPLAN). Book Wild Animal Response Strategy (Version 3.2). Australian Veterinary Emergency Plan (AUSVETPLAN). 2005, Primary Industries Ministerial Council
Hone J, Pech R: Disease sureillance in wildlife with emphasis on detecting foot and mouth disease in feral pigs. J Environ Manage. 1990, 31: 173-184. 10.1016/S0301-4797(05)80007-6.
Elbers ARW, Stegeman A, Moser H, Ekker HM, Smak JA, Pluimers FH: The classical swine fever epidemic 1997-1998 in The Netherlands: descriptive epidemiology. Prev Vet Med. 1999, 42: 157-184. 10.1016/S0167-5877(99)00074-4.
Cannon RM, Roe RT: Livestock Disease Surveys: a Field Manual for Veterinarians. 1982, Canberra Australian Government Publishing Service
Faust A, Lange E, Kaden V: Efficacy of lyophilised C-strain vaccine after oral immunisation of domestic pigs and wild boar against classical swine fever: first results. Dtsch Tierarztl Wochenschr. 2007, 114: 412-417.
Kaden V, Lange B: Oral immunisation against classical swine fever (CSF): onset and duration of immunity. Vet Microbiol. 2001, 82: 301-310. 10.1016/S0378-1135(01)00400-X.
Wackerly DD, Mendenhall W, Schaeffer RL: Mathematical Ssatistics with applications. 2008, Belmont: Thomson Higher Education, 7
Frey HC, Patil SR: Identification and review of sensitivity analysis methods. Risk Anal. 2002, 22: 553-578. 10.1111/0272-4332.00039.
Rossi S, Artois M, Pontier D, Crucière C, Hars J, Barrat J, Pacholek X, Fromont E: Long-term monitoring of classical swine fever in wild boar (Sus scrofa sp.) using serological data. Vet Res. 2005, 36: 27-42. 10.1051/vetres:2004050.
Caley P: Population dynamics of feral pigs (Sus scrofa) in a tropical riverine habitat complex. Wild Res. 1993, 20: 625-636. 10.1071/WR9930625.
Dexter N: The influence of pasture distribution and temperature on habitat selection by feral pigs in semi-arid environment. Wild Res. 1998, 25: 547-559. 10.1071/WR97119.
Haydon DT, Randall DA, Matthews L, Knobel DL, Tallents LA, Gravenor MB, Williams SD, Pollinger JP, Cleaveland S, Woolhouse MEJ, Sillero-Zubiri C, Marino J, Macdonald DW, Laurenson MK: Low-coverage vaccination strategies for the conservation of endangered species. Nature. 2006, 443: 692-695. 10.1038/nature05177.
Mangen MJJ, Jalvingh AW, Nielen M, Mourits MCM, Klinkenberg D, Dijkhuizen AA: Spatial and stochastic simulation to compare two emergency-vaccination strategies with a marker vaccine in the 1997/1998 Dutch Classical Swine Fever epidemic. Prev Vet Med. 2001, 48: 177-200. 10.1016/S0167-5877(00)00195-1.
Mitchell J, Dorney W, Mayer R, McIlroy J: Migration of feral pigs (Sus scrofa) in rainforests of north Queensland: fact or fiction?. Wild Res. 2009, 36: 110-116. 10.1071/WR06066.
Kern B, Depner KR, Letz W, Rott M, Thalheim S, Nitschke B, Plagemann R, Liess B: Incidence of classical swine fever (CSF) in wild boar in a densely populated area indicating CSF virus persistence as a mechanism for virus perpetuation. Zentralbl Veterinarmed B. 1999, 46: 63-67.
Depner KR, Mueller A, Gruber A, Rodriguez A, Bickhardt K, Liess B: Classical swine fever in wild boar (Sus scrofa)-experimental infections and viral persistence. Dtsch Tierarztl Wochenschr. 1995, 102: 381-384.
Hone J: Note on seasonal changes in population density of feral pigs in three tropical habitats. Australian Wildlife Research. 1990, 17: 131-134. 10.1071/WR9900131.
Miller WM: A state-transition model of epidemic foot-and-mouth disease. In book. Potential Economic Impact of Foot-and-Mouth Disease in the United States. 1979, US Government Printing Office, Washington DC
Sargent RG: Verification and validation of simulation models. In Proceedings of the Winter Simulation Conference. Edited by: Medeiros DJ, Watson EF, Carson JS, Manivannan MS. 1998, Washington DC, USA, 121-130.
Garner M, Hamilton S: Principles of epidemiological modelling. Review Scientifique et Technique. In press
James AD, Rossiter PB: An epidemiological model of rinderpest I. Description of the model. Tropical Animal Health and Production. 1989, 21: 59-68. 10.1007/BF02297347.
Garner MG, Cowled B, East IJ, Moloney BJ, Kung NY: Evaluating the effectiveness of early vaccination in the control and eradication of equine influenza-- A modelling approach. Preventive Veterinary Medicine. 2011, 99: 15-27. 10.1016/j.prevetmed.2010.02.007.
Saunders G: The ecology and management of feral pigs in New South Wales. 1988, PhD thesis. Macquarie University
Depner KR, Muller A, Gruber A, Rodriguez A, Bickhardt K, Liess B: Classical swine fever in wild boar (Sus Scrofa) -Experimental Infections and Viral Persistence. Deutsche Tierarztliche Wochenschrift. 1995, 102: 381-384.
Leforban Y, Cariolet R: Characterization and pathogenicity for pigs of a hog cholera virus- strain isolated from wild boars. Annales De Recherches Veterinaires. 1992, 23: 93-100.
Artois M, Depner KR, Guberti V, Hars J, Rossi S, Rutili D: Classical swine fever (hog cholera) in wild boar in Europe. Revue Scientifique Et Technique De L Office International Des Epizooties. 2002, 21: 287-303.
Kaden V, Lange E, Steyer H: Does multiple oral vaccination of wild boar against classical swine fever (CSF) have a positive influence on the immunity?. Deutsche Tierarztliche Wochenschrift. 2004, 111: 63-67.
Matthaeus W, Korn G: Different types of neutralizing antibodies in swine serum following the inoculation of hog cholera virus. Zent Bak Teriol Parasitenk Infektionskrankh Hyg Abt I Orig. 1967, 204: 173-180.
Jezierski W: Longevity and mortality-rate in a population of wild boar. Acta Theriologica. 1977, 22: 337-348.
Vandeputte J, Too HL, Ng FK, Chen C, Chai KK, Liao GA: Adsorption of colostral antibodies against classical swine fever, persistence of maternal antibodies, and effect on response to vaccination in baby pigs. American Journal of Veterinary Research. 2001, 62: 1805-1811. 10.2460/ajvr.2001.62.1805.
Depner KR, Muller T, Lange E, Staubach C, Teuffert J: Transient classical swine fever virus infection in wild boar piglets partially protected by maternal antibodies. DTW (Deutsche Tieraerztliche Wochenschrift). 2000, 107: 66-68.
Kaden V, Lange E, Fischer U, Strebelow G: Oral immunisation of wild boar against classical swine fever: evaluation of the first field study in Germany. Veterinary Microbiology. 2000, 73: 239-252. 10.1016/S0378-1135(00)00148-6.
We gratefully acknowledge funding from the Cattle Council of Australia, Meat and Livestock Australia and the Australian Government Department of Agriculture, Fisheries and Forestry. This research was also supported under Australian Research Council's Linkage Projects funding scheme (LP100200110). We thank Sam Beckett and Ben Madin for fruitful discussions when designing our model.
The authors declare that they have no competing interests.
BC, GG and KN designed the model logic. BC and GG coded the model. BC, KN and MW conducted analyses. BC prepared the first draft of the manuscript with contributions from KN, MW and GG. MW and GG edited the manuscript. All authors read and approved the final manuscript.