Estimating the Location and Spatial Extent of a Covert Anthrax Release

Rapidly identifying the features of a covert release of an agent such as anthrax could help to inform the planning of public health mitigation strategies. Previous studies have sought to estimate the time and size of a bioterror attack based on the symptomatic onset dates of early cases. We extend the scope of these methods by proposing a method for characterizing the time, strength, and also the location of an aerosolized pathogen release. A back-calculation method is developed allowing the characterization of the release based on the data on the first few observed cases of the subsequent outbreak, meteorological data, population densities, and data on population travel patterns. We evaluate this method on small simulated anthrax outbreaks (about 25–35 cases) and show that it could date and localize a release after a few cases have been observed, although misspecifications of the spore dispersion model, or the within-host dynamics model, on which the method relies can bias the estimates. Our method could also provide an estimate of the outbreak's geographical extent and, as a consequence, could help to identify populations at risk and, therefore, requiring prophylactic treatment. Our analysis demonstrates that while estimates based on the first ten or 15 observed cases were more accurate and less sensitive to model misspecifications than those based on five cases, overall mortality is minimized by targeting prophylactic treatment early on the basis of estimates made using data on the first five cases. The method we propose could provide early estimates of the time, strength, and location of an aerosolized anthrax release and the geographical extent of the subsequent outbreak. In addition, estimates of release features could be used to parameterize more detailed models allowing the simulation of control strategies and intervention logistics.Author Summary: Releasing highly pathogenic organisms into an urban population is a form of bioterrorism that could result in a large number of casualties. The first indication that a covert open-air release has occurred is quite likely to be individuals reporting for medical attention. If such an attack is suspected, then public health authorities would attempt to identify those individuals who have been infected in order to provide rapid treatment with the aim of reducing the possibility of disease and potential death. Aiming treatment at too small an area might miss individuals infected further down and/or up wind, whereas issues surrounding both treatment resources and serious side effects may rule out mass treatment campaigns of large sections of the population. Our work provides scientific robustness to firstly estimate where and when an aerosolized release has occurred and secondly identify the most critically affected geographic areas. In order to use this statistical tool during an outbreak, public health workers would only need to collect the time of symptomatic onset and the home and work locations of early cases; recent weather information would also be required. Although the accuracy of the estimates is likely to improve as more cases appear, treating individuals based on early estimates might prove more beneficial since time would be of the essence.