key: cord-0868266-655eicxb authors: Wang, Yulan; Li, Bernard; Gouripeddi, Ramkiran; Facelli, Julio C. title: Human Activity Pattern Implications for Modeling SARS-CoV-2 Transmission date: 2020-12-08 journal: Comput Methods Programs Biomed DOI: 10.1016/j.cmpb.2020.105896 sha: 05a36f152c3c15aa918b38041865fd900a4f75a4 doc_id: 868266 cord_uid: 655eicxb BACKGROUND AND OBJECTIVES: : SARS-CoV-2 emerged in December 2019 and rapidly spread into a global pandemic. Designing optimal community responses (social distancing, vaccination) is dependent on the stage of the disease progression, discovery of asymptomatic individuals, changes in virulence of the pathogen, and current levels of herd immunity. Community strategies may have severe and undesirable social and economic side effects. Modeling is the only available scientific approach to develop effective strategies that can minimize these unwanted side effects while retaining the effectiveness of the interventions. METHODS: : We extended the agent-based model, SpatioTemporal Human Activity Model (STHAM), for simulating SARS-CoV-2 transmission dynamics. RESULTS: : Here we present preliminary STHAM simulation results that reproduce the overall trends observed in the Wasatch Front (Utah, United States of America) for the general population. The results presented here clearly indicate that human activity patterns are important in predicting the rate of infection for different demographic groups in the population. CONCLUSIONS: : future work in pandemic simulations should use empirical human activity data for agent-based techniques. SARS-CoV-2 emerged in December 2019 and rapidly spread into a global pandemic (1) (2) (3) . Often considered strategies for reducing local SARS-CoV-2 spread are social distancing (4-6) and more aggressive lockdowns (7) . The impacts of social distancing and other public health measures including vaccination (8, 9) and prophylactic medications (when these become available), for designing optimal community responses will be dependent on the stage of disease progression within the community (10), discovery of asymptomatic individuals (11, 12) , changes in virulence of the pathogen (13) and current levels of herd immunity. In some countries, there have been efforts to let the disease spread to develop a level of herd immunity, as well as calls for younger and low-risk individuals to join back into the workforce. The efficacy as well as the undesired effects of these approaches are largely not known (13) . As observed in previous pandemics, social determinants of health (SDOH) are likely to influence the ability of individuals and families to follow SARS-CoV-2 public health guidelines and may alter health outcomes (14, 15) . Impacts of pandemics are amplified in lower socioeconomic, densely populated urban areas compared to distributed rural areas (14, 16) . SDOH like education levels, household size, age-structure, population density, workplaceschoolhouseholdneighborhood structures, and area-level poverty, along with behaviors such as cigarette smoking that impact health, can directly and indirectly affect both immediate, and long term outcomes in pandemics (17, 18) . The effects of selective movement of individuals across time and demographic groups in the evolution of disease are difficult to ascertain. Modeling is the only viable scientific approach to develop effective strategies that can minimize these undesirable side effects and show the effectiveness of the interventions (19, 20) . With the principal cause of SARS-CoV-2 spread being movement and interaction of individuals within a community (21) , modeling should be informed by human activity patterns as different humans will come in contact with infected individuals depending upon their daily activities. Numerous simulation models have been reported, and most of them have been implemented in EpiModel (22) . Many papers have presented simulations of SARS-CoV-2 spread using a variety of mathematical and computational models; examples of this work can be found in Refs. (23) (24) (25) (26) (27) (28) (29) (30) (31) (32) . Of relevance to the work reported here, several agent base models have been used to predict SARS-CoV-2 spread (33, 34) . Still, to the authors' knowledge, none of these models have incorporated empirical information about human activity patterns in their simulations. We have shown that human activity patterns, as recorded by the American Time Use Survey (ATUS) (35) , can be classified and used to simulate human activity in the context of estimating personal exposure to PM 2.5 (36, 37) . Based on this classification of human activities, we have developed the SpatioTemporal Human Activity Model, or STHAM (37), to generate and characterize travel trajectories and activity patterns (38) (39) (40) , and to integrate exposure-agent The STHAM simulation method is described in detail in Ref. (37) . Here for completeness, we briefly describe the key steps to produce a STHAM simulation. The STHAM initialization process is detailed in Figure 1 . At the core of the STHAM is the agent, which is represented as a single data record with demographic properties representing a hypothetical person. These properties (e.g., age, gender, census block location) are assigned to each agent randomly using the aggregate statistics at the census block level as constraints for the 2010 United States Census (42-44) so that the generated agents and the households to which they are assigned match the data provided in the census. Each household is then assigned a location in a round-robin fashion based on a list of physical postal addresses (45) that are matched to each census block. Next, each agent is assigned an employment status and, if relevant, an employment location, based on data from the Longitudinal Employee Household Dynamics Program (under the US Census Bureau) (LEHD) (46). The LEHD also provides an approximation for employment location at the census block level, which allows the model to simulate commuter patterns accurately. Finally, school participation is assigned using age, and gender stratified enrollment rates from the American Community Survey (ACS) (43). there is a significant intersection between workers and students, agents are classified as Worker and Student at the same time, and the sum of these three classes is larger than the total population. Based on the STHAM model, we calculated the daily activity patterns for all agents, including activity categories and activity locations. There are two kinds of daily life activities patterns in our model: one for weekdays and another is the weekends. To simulate one week, we run the model with five weekdays and two weekends. We divided the 24 hours in each day into 288 five-minute segments and aggregated agents within a location at a five-minute interval. For any colocation, if there are one or more infected agents, that location is considered as an infected location, and any susceptible agents (never infected before) in these locations have a non-zero probability of being infected. This follows the current understanding that there have been no cases of SARS-CoV-2 reinfections, but could be modified as knowledge of SARS-CoV-2 immunity evolves (47) . A susceptible agent in an infected location has its activity contact index randomly sampled by a uniform distribution from 0 -1. If the agent's contact index is lower than the transmission rate, this agent's status changes to infected. The simulation used a transmission rate defined by a Gaussian distribution and having a mean of 2.5/(288*14*30), and a standard error of 0.00001. We used a basic reproduction number (R 0 ) of 2.5 over the span of 14 days in the model, which is based on current best estimates for R 0 and the average length of disease as per the Centers for Disease Control and Prevention (48). The STHAM model provides an average daily inter-agent contact rate of 30 (36) . Therefore, the mean value of the transmission rate is set to 2.5/(288*14*30). The infection process repeats every five-minutes, and after oneday simulation, the list of infected agents is updated by adding the new infected agents and excluding the recovered agents. Infected individuals remain infectious for 14 days. After this period, they are neither infective nor susceptible to reinfection. We arbitrary assigned five Workers, three Students, and two OtherAgents to be infected at the beginning of the simulations, and ignored all inbound and outbound travel from the region. We simulated the epidemic for 60 days. The EpiModel (22) package utilizes a stochastic individual contact model (ICM), which simulates agent actions at an individual level for each time step. The ICM has a threecompartment model, where individuals in a population are assigned to compartments with labels -S, I, and R (Susceptible, Infected, and Recovered). The progression between compartments is dictated by user-specified percentages that are used to define random draws from a binomial distribution. Progression between compartments depends on actions taken by infected agents. In EpiModel, this is defined as the act_rate, or the number of both effective and ineffective contacts made by an infected agent. For the scenario reported here, the act_rate was set to 0.25 acts/day, which is consistent with the infecting probability used in the STHAM simulations. After an individual spends 14 days in the infected compartment, they are moved to the recovered compartment and treated as non-acting agents, which are neither infective nor susceptible to reinfection. A simulation of 10 weeks was run ten times, using randomly selected seeds, and the results for each simulation were aggregated and averaged to obtain the final results reported here. For starting parameters, there were 2,000,000 susceptible agents, ten infected agents, and 0 recovered agents assigned at the start of the simulation. The number of infected individuals in the Wasatch Front, which is the most densely populated region of the State of Utah, United States of America, with a population of approximately 2,000,000 individuals, was obtained from the Utah Department of Health (UDOH). The first day in which the cumulative number of infected cases reported to UDOH was over ten individuals was March 15 th 2020, with 18 cases. The raw values obtained from UDOH were scaled by a factor 1.8 to be consistent with the simulated results obtained with an initial population of 10 infected individuals at day zero, which was set to March 15 th 2020. While this study demonstrates the need of using empirical activity patterns to obtain higher accuracy in modeling epidemics, the study also has limitations that need to be recognized. The activities and mobility used here are those expected from people in the United States of America under normal times. These may be significantly different activity patterns than those expected during a pandemic, and better models based on data of human activity collected during an epidemic (e.g., COVID-19 Community Mobility Reports (50)) are necessary for developing more accurate predicted models. The demographic categories used here are quite broad, and future work should be the focus in developing more granular models for the demographics, these models should be grounded in empirical data, which can be transformed into probabilistic agent behaviors using the STHAM approach (36, 37) . For simplicity, here we used a SIR compartment model (Susceptible-Infected-Recovered), but we acknowledge that using a SEIR (Susceptible-Exposed-Infectious-Recovered) model could provide more accurate results (51-53). Our modeling software is able to accommodate any infection model including different rates of symptomatic/asymptomatic rates and scenarios that are potentially dependent on the age of the agent, load of exposure, immunization status, use of precautionary measures (e.g. agents wearing masks) among other factors. As next steps, we will also evaluate the effects of the environment (54, 55) and selective movement of individuals across time and demographic groups in the evolution of the disease for the State of Utah with different mitigation strategies and demonstrate that our model, taking into account the empirical evidence of human mobility, produces more detailed information than models homogeneously treating individual agents. The principal cause of the spread of SARS-COV-2 is the movement and interaction of individuals within a community. The STHAM model (36, 37) assigns different mobility patterns to an individual belonging to different demographic groups as defined by empirical activity sequences derived from the ATUS (35) . Here, we show that the STHAM modeling approach can reproduce the overall trends of the spread of SARS-COV-2 in the Utah Wasatch front, but also show that different demographic groups have significantly different infection rates. Future work should focus on using empirical data on human activity patterns for pandemic simulations using agent-based techniques and including more categories of agents to get higher specificity on the lead activates that contribute to the spread of the virus. Developing this kind of models will allow policymakers to make contention strategies than may be much more granular and leading to less undesirable effects. JCF and RG conceived the research; YW and BL performed the calculations and all authors analyzed the data, developed the manuscript and approved the submission. All authors agree to be accountable for the contents of the paper. This research was partially funded by a Covid-19 -Navigating the Uncharted WHO Declares COVID-19 a Pandemic COVID-19: Knowns, Unknowns, and Questions. mSphere Public health interventions and epidemic intensity during the 1918 influenza pandemic National Academies of Sciences E, Medicine. Rapid Expert Consultation on Social Distancing for the COVID-19 Pandemic Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus-Infected Pneumonia Preliminary Identification of Potential Vaccine Targets for the COVID-19 Coronavirus (SARS-CoV-2) Based on SARS-CoV Immunological Studies Immune responses in COVID-19 and potential vaccines: Lessons learned from SARS and MERS epidemic Fundamental principles of epidemic spread highlight the immediate need for large-scale serological surveys to assess the stage of the SARS-CoV-2 epidemic Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV-2) Covid-19: identifying and isolating asymptomatic people helped eliminate virus in Italian village Covid-19: experts question analysis suggesting half UK population has been infected Socio-economic disparities in mortality due to pandemic influenza in England Protecting vulnerable populations from pandemic influenza in the United States: a strategic imperative Does Density Aggravate the COVID-19 Pandemic The social determinants of health and pandemic H1N1 2009 influenza severity Social Determinants of Health, Disaster Vulnerability, Severe and Morbid Obesity in Adults: Triple Jeopardy? The Next Epidemic -Lessons from Ebola Innovation for Pandemics COVID-19 infection: Origin, transmission, and characteristics of human coronaviruses EpiModel: An R Package for Mathematical Modeling of Infectious Disease over Networks Stability analysis and numerical simulation of SEIR model for pandemic COVID-19 spread in Indonesia A fractional differential equation model for the COVID-19 transmission by using the Caputo-Fabrizio derivative A review on COVID-19: Origin, spread, symptoms, treatment, and prevention Evolutionary modelling of the COVID-19 pandemic in fifteen most affected countries A mathematical model of COVID-19 using fractional derivative: outbreak in India with dynamics of transmission and control. Advances in Difference Equations A mathematical model for COVID-19 transmission by using the Caputo fractional derivative Study of transmission dynamics of novel COVID-19 by using mathematical model Simulating the spread of COVID-19 via a spatially-resolved susceptible-exposed-infected-recovered-deceased (SEIRD) model with heterogeneous diffusion A numerical simulation of fractional order mathematical modeling of COVID-19 disease in case of Wuhan China COVID-19) Mathematical Modeling. Centers for Disease Control and Prevention An agent-based model to evaluate the COVID-19 transmission risks in facilities TRACE (Testing Responses through Agent-based Computational Epidemiology Generation and Classification of Activity Sequences for Spatiotemporal Modeling of Human Populations STHAM: an agent based model for simulating human exposure across high resolution spatiotemporal domains SACSIM: An applied activity-based model system with fine-level spatial and temporal resolution Comprehensive Econometric Microsimulator for Daily Activity-Travel Patterns Implementation Framework and Development Trajectory of FEATHERS Activity-Based Simulation Platform An agent-based model for estimating human activity patterns on the Wasatch Front. Air Quality: Science for Solutions; Salt Lake City2017. 42. 2015 TIGER/Line Shapefiles (machine-readable data files) Census demographic Profile summary File (machine-readable data files). United States Census Bureau / American FactFinder LEHD Origin-Destination Employment Statistics (LODES) Dataset Version 7.3 (machinereadable data files) COVID-19 and Postinfection Immunity: Limited Evidence, Many Remaining Questions Estimation of Unreported Novel Coronavirus (SARS-CoV-2) Infections from Reported Deaths: A Susceptible-Exposed-Infectious-Recovered-Dead Model Modelling the COVID-19 epidemic and implementation of population-wide interventions in Italy Spread and dynamics of the COVID-19 epidemic in Italy: Effects of emergency containment measures Modeling the Novel Coronavirus (SARS-CoV-2) Outbreak in Sicily, Italy Will Coronavirus Pandemic Diminish by Summer? Spread of SARS-CoV-2 Coronavirus likely to be constrained by climate Special thanks to the State of Utah Data and Analysis and Support Team that was composed of individuals from the Department of Health, Governor's Office of Management and Budgets, Department of Technology Services, Local Health Departments, and the Legislature for providing the relevant data. We would also like to thank Dr. Albert Lund for his inputs in extending his original work on STHAM for SARS-CoV-2 transmission.