PT - JOURNAL ARTICLE AU - Rogan, Peter K. AU - Mucaki, Eliseos J. AU - Lu, Ruipeng AU - Shirley, Ben C. AU - Waller, Edward AU - Knoll, Joan H.M. TI - Meeting radiation dosimetry capacity requirements of population-scale exposures by geostatistical sampling AID - 10.1101/2020.04.08.20058446 DP - 2020 Jan 01 TA - medRxiv PG - 2020.04.08.20058446 4099 - http://medrxiv.org/content/early/2020/04/14/2020.04.08.20058446.short 4100 - http://medrxiv.org/content/early/2020/04/14/2020.04.08.20058446.full AB - Background Accurate radiation dose estimates are critical for determining eligibility for therapies by timely triaging of exposed individuals after large-scale radiation events. However, the universal assessment of a large population subjected to a nuclear spill incident or detonation is not feasible. Even with high-throughput dosimetry analysis, test volumes far exceed the capacities of first responders to measure radiation exposures directly, or to acquire and process samples for follow-on biodosimetry testing.Aim To significantly reduce data acquisition and processing requirements for triaging of treatment-eligible exposures in population-scale radiation incidents.Methods Physical radiation plumes modelled nuclear detonation scenarios of simulated exposures at 22 US locations. Models assumed only location of the epicenter and historical, prevailing wind directions/speeds. The spatial boundaries of graduated radiation exposures were determined by targeted, multistep geostatistical analysis of small population samples. Initially, locations proximate to these sites were randomly sampled (generally 0.1% of population). Empirical Bayesian kriging established radiation dose contour levels circumscribing these sites. Densification of each plume identified critical locations for additional sampling. After repeated kriging and densification, overlapping grids between each pair of contours of successive plumes were compared based on their diagonal Bray-Curtis distances and root-mean-square deviations, which provided criteria (<10% difference) to discontinue sampling.Results/Conclusions We modeled 30 scenarios, including 22 urban/high-density and 2 rural/low-density scenarios under various weather conditions. Multiple (3-10) rounds of sampling and kriging were required for the dosimetry maps to converge, requiring between 58 and 347 samples for different scenarios. On average, 70±10% of locations where populations are expected to receive an exposure ≥2Gy were identified. Under sub-optimal sampling conditions, the number of iterations and samples were increased and accuracy was reduced. Geostatistical mapping limits the number of required dose assessments, the time required, and radiation exposure to first responders. Geostatistical analysis will expedite triaging of acute radiation exposure in population-scale nuclear events.Competing Interest StatementPeter K. Rogan and Joan H. Knoll are founders and Ben Shirley is an employee of CytoGnomix, which partially supported this research. Cytognomix holds relevant patents and relevant patent applications.Funding StatementWe are grateful to the SOSCIP Smart Computing for Innovation Consortium (J.H.M.K., P.K.R.), Natural Sciences and Engineering Research Council of Canada (Engage Program; E.W., P.K.R), Ontario Centres of Excellence (Talent Edge Postdoctoral Fellowship Program; J.H.M.K., P.K.R.), CytoGnomix (P.K.R.) for support of this project. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.Author DeclarationsAll relevant ethical guidelines have been followed; any necessary IRB and/or ethics committee approvals have been obtained and details of the IRB/oversight body are included in the manuscript.YesAll necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived.YesI understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).YesI have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.YesAll data and materials are fully available without restriction. Data and Software from this study are available in the Zenodo repository below. This includes modified U.S. state and sub-division boundary files [in KML format], as well as the geographic coordinates and dose values (World Geodetic System [WGS] 1984) which define each HPAC and geostatistical-derived plume for all scenarios described. The programs include custom scripts used to preprocess the data for use with ArcMap GIS software. https://doi.org/10.5281/zenodo.3572574 http://dx.doi.org/10.17504/protocols.io.ba4nigve