key: cord-0798062-rq3k9bc5 authors: Shah, Nita H.; Sheoran, Nisha; Jayswal, Ekta; Shukla, Dhairya; Shukla, Nehal; Shukla, Jagdish; Shah, Yash title: Modelling COVID-19 transmission in the United States through interstate and foreign travels and evaluating impact of governmental public health interventions date: 2020-12-22 journal: J Math Anal Appl DOI: 10.1016/j.jmaa.2020.124896 sha: ce6776d89edf08f0ca5f53cc97be94c11fb58493 doc_id: 798062 cord_uid: rq3k9bc5 Background: The first case of COVID-19 was reported in Wuhan, China in December 2019. The disease has spread to 210 countries and has been labelled as a pandemic by the World Health Organization (WHO). Modelling, evaluating, and predicting the rate of disease transmission is crucial in understanding optimal methods for prevention and control. Our aim is to assess the impact of interstate and foreign travel and public health interventions implemented by the United States government in response to the COVID-19 pandemic. Methods: A disjoint mutually exclusive compartmental model was developed to study transmission dynamics of the novel coronavirus. A system of nonlinear differential equations was formulated and the basic reproduction number [Formula: see text] was computed. Stability of the model was evaluated at the equilibrium points. Optimal controls were applied in the form of travel restrictions and quarantine. Numerical simulations were conducted. Results: Analysis shows that the model is locally asymptomatically stable, at endemic and foreigners free equilibrium points. Without any mitigation measures, infectivity and subsequent hospitalization of the population increased. When interstate and foreign travel was restricted and the population placed under quarantine, the probability of exposure and subsequent infection decreased significantly; furthermore, the recovery rate increased substantially. Conclusion: Interstate and foreign travel restrictions, in addition to quarantine, are necessary in effectively controlling the pandemic. The United States has controlled COVID-19 spread by implementing quarantine and restricting foreign travel. The government can further strengthen restrictions and reduce spread within the nation more effectively by implementing restrictions on interstate travel. Background: The first case of COVID-19 was reported in Wuhan, China in December 2019. The disease has spread to 210 countries and has been labelled as a pandemic by the World Health Organization (WHO). Modelling, evaluating, and predicting the rate of disease transmission is crucial in understanding optimal methods for prevention and control. Our aim is to assess the impact of interstate and foreign travel and public health interventions implemented by the United States government in response to the COVID-19 pandemic. Methods: A disjoint mutually exclusive compartmental model was developed to study transmission dynamics of the novel coronavirus. A system of nonlinear differential equations was formulated and the basic reproduction number R 0 was computed. Stability of the model was evaluated at the equilibrium points. Optimal controls were applied in the form of travel restrictions and quarantine. Numerical simulations were conducted. Results: Analysis shows that the model is locally asymptomatically stable, at endemic and foreigners free equilibrium points. Without any mitigation measures, infectivity and subsequent hospitalization of the population increased. When interstate and foreign travel was restricted and the population placed under quarantine, the probability of exposure and subsequent infection decreased significantly; furthermore, the recovery rate increased substantially. Conclusion: Interstate and foreign travel restrictions, in addition to quarantine, are necessary in effectively controlling the pandemic. The United States has controlled COVID-19 spread by implementing quarantine and restricting foreign travel. The government can further strengthen restrictions and reduce spread within the nation more effectively by implementing restrictions on interstate travel. © 2020 Elsevier Inc. All rights reserved. 29th March Social distancing extended through 30th April. 3rd April All American wear non-medical, fabric or cloth masks to prevent asymptomatic spread of coronavirus. In December 2019, an unidentified pneumonia was found in Wuhan, Hubei province, China. The responsible virus was later identified as the novel coronavirus 2019, and the disease as coronavirus disease 2019 (COVID-19) [4] . COVID-19 is transmitted via direct contact with an infected person through respiratory droplets when a person coughs or sneezes, or by indirect contact with contaminated surfaces with respiratory droplets from infected person and then touching their eyes, nose or mouth [6] . Furthermore, the virus remains on surfaces from a few hours to several days and has an incubation period between 1-14 days. Consequently, the disease spread rapidly from Wuhan to all parts of the country and overseas. Introduction and spread of COVID-19 within the United States is a direct result of transmission through foreign and interstate travel. The first known case of COVID-19 in the U.S.A. was confirmed on 20th January 2020 in a 35-year-old individual who had travelled from Wuhan to Washington state [6] . Soon, cases started appearing and rising in many other states of the U.S.A. including New York, New Jersey, Illinois, Florida, Georgia, Texas, Pennsylvania due to interstate travel. The CDC alarmed that hospitals may get overwhelmed by a large number of people seeking care at the same time due to widespread transmission of disease which may lead to otherwise preventable deaths (2020) [2] . In response to the oncoming epidemic the US government implemented the following regulations (Table 1) . While the United States has implemented numerous public health interventions, it has not implemented a ban on interstate travel. According to the World Health Organization (WHO) [12] . New cases of COVID-19 have emerged in 210 countries with 1,733, 945 confirmed cases and 106, 518 confirmed deaths globally as of 10th April 2020. The United States has implemented quarantine measures, close contact tracing, early testing for individuals with symptoms, hospitalization if needed, and closing of teaching institutes and non-essential businesses. Studies have shown that in other countries, the complete lockdown of travel has decreased the spread of the disease in the surrounding states ([1]; [10] ). In order to prevent the transmission of COVID-19 Rate at which foreigner quarantine themselves 0.000015 β 7 Rate at which foreigner gets infected 0.000001 β 8 Rate at which quarantine humans gets infected by COVID-19 0.0025 Rate at which infected humans gets hospitalized 0.0037 β 10 Rate at which exposed humans gets hospitalized 0.000002 β 11 Rate at which hospitalized humans gets recovered 0.0000003 μ Death rate 0.000119 μ c Death rate due to COVID-19 0.0027 within the US, the mode of transmission must first be modelled and understood. Mathematical modelling is ideal for evaluating and predicting the rate of disease transmission. Data-driven mathematical modelling plays an important role in epidemic mitigation, in preparedness for future epidemic and in the evaluation of control effectiveness [13] . In this study, we adopted a disjoint mutually exclusive compartmental model to shed light on the transmission dynamics from foreign and interstate travel of the novel coronavirus and our aim is to assess the impact of public health interventions on infection by measuring basic reproduction number, contact rate, newly confirmed cases, total confirmed cases, total death. Our estimated parameters are largely in line with World Health Organization estimates and previous studies (2019). Mathematical modelling plays a vital role in determining dynamics of diseases. In this paper we consider a disjoint mutually exclusive compartmental model with compartments as follows: Exposed to COVID-19 E i.e. this compartment consists of individuals (Both foreign population and interstate population) who are exposed to COVID-19, next compartment is I S i.e. transmission of COVID-19 through interstate travel, COVID-19 transmission through foreigners F -this compartment includes U.S. population which is exposed to COVID-19, Quarantined class Q, COVID-19 Infected I -this class includes infected population as well infectious population, hospitalized H -this class includes hospitalization of both COVID-19 infectives and also those who are exposed to COVID-19 and last compartment includes recovered population from hospitalized population denoted by R. Notations and parametric values used in the formulation of dynamical system model are given in the following Table 2 . This model considers new recruitment in the exposed class at the rate B and all the compartments have mortality rate μ. Here β 1 is the US population exposed to COVID-19 via interstate travel, and β 2 is the rate at which the US population gets exposed to COVID-19 through contact with foreigners. Next, β 3 is the rate of COVID-19 exposure through foreigners engaged in interstate travel. The US population engaging in interstate travel and foreigners quarantine themselves at the rate β 4 and β 6 respectively. Similarly, after getting exposed to COVID-19, US population engaging in interstate travel and foreigner population gets the infection joining infectious class I with the rate β 5 and β 7 respectively. Quarantined humans also get the infection at the rate β 8 . Next, we assume infected population gets hospitalized at the rate β 9 joining H. We also assume population gets admitted to the hospital at the initial exposure of the disease at the rate β 10 . Hospitalized patients after undergoing treatment gets recover joining R at the rate β 11 . We take into consideration death due to COVID-19 μ c , when the individual is hospitalized. The Fig. 1 gives rise to the following set of non-linear ordinary differential equations Adding all the differential equations of model, we get, μ . Then, Feasible Region for the system is defined as with This system has following equilibrium points i. Foreigner free equilibrium point ii. Endemic equilibrium point Basic reproduction number is defined as the total number of secondary infections in a total susceptible population. Here, we calculate the reproduction number using Diekmann et al., when the disease in its endemic stage i.e. for this model it is defined as percentage of population infected by a single infection in a totally exposure situation [5] . We also compute the value of reproduction number R F when there are no foreigners present in the total population. The value of R E * = 81% and R F = 1.10. In this section we study the stability analysis of the model. Here we study Local stability of all the equilibrium points using Routh-Hurwitz criterion by Routh 1877 [9] . Theorem 1. The foreigner free equilibrium point is locally asymptotically stable if (β 4 + β 5 + μ) < max Bβ 1 β 10 +μ , β 1 (β 3 +β 6 +β 7 +μ) Proof. The Jacobian of system (1) at Foreigner free equilibrium is as follows The eigen values of the Jacobian J F are If it has imaginary roots i.e. ξ < 0. Then we have negative real part. Hence the theorem. But if ξ ≥ 0, then eigen values are negative if (β 4 + β 5 + μ) < max Bβ 1 β 10 +μ , β 1 (β 3 +β 6 +β 7 +μ) β 2 . Hence the Foreigner free equilibrium point is locally asymptotically stable. Proof. The Jacobian matrix of system (1) for endemic equilibrium is given by where, a 11 = β 1 I * S +β 2 F * +β 10 +μ, a 22 = −β 1 E * +β 4 +β 5 +μ, a 33 = −β 2 E * +β 3 +β 6 +β 7 +μ, a 44 = β 8 +μ, a 55 = β 9 + μ, a 66 = β 11 + μ C + μ. The characteristic polynomial for Jacobian J * is (a 44 a 55 a 66 + a 44 a 55 + a 44 a 66 + a 55 a 66 ) a 44 + a 55 ) + a 44 a 55 ) + a 11 a 22 ((a 33 a 44 a 55 (a 66 + μ)) + a 66 μ(a 33 (a 44 + a 55 ) + a 44 a 55 )) + a 33 a 44 a 55 a 66 μ(a 11 + a 22 ) + (a 66 + μ)E * a 44 a 55 (F * β 2 2 a 22 + I S β 2 1 a 33 ) b 2 = E * F * β 2 2 ((a 66 + μ)(a 22 (a 44 + a 55 ) + a 44 a 55 ) + a 66 μ(a 22 + a 44 + a 55 ) + a 22 a 44 a 55 ) + E * F * β 1 β 2 β 3 ((a 66 + μ)(a 44 + a 55 ) + a 66 μ + a 44 a 55 ) + E * I * S β 2 1 ((a 66 + μ)(a 33 (a 44 + a 55 ) + a 44 a 55 ) + a 66 μ(a 33 + a 44 + a 55 ) + a 33 a 44 a 55 ) + (a 66 + μ)(a 11 a 22 (a 44 + a 33 (a 44 + a 55 )) + a 33 a 44 a 55 (a 11 + a 22 )) + (a 44 + a 55 )(a 11 a 66 μ(a 22 + a 33 ) + a 22 a 33 a 66 μ) + a 44 a 55 a 66 μ(a 22 + a 33 ) + a 11 a 66 μ(a 44 a 55 + a 22 a 33 ) + a 11 a 22 a 33 a 44 a 55 ) b 3 = E * F * β 2 2 ((a 66 + μ)(a 22 + a 44 + a 55 ) + a 66 μ + a 22 (a 44 + a 55 )) + E * F * β 1 β 2 β 3 (a 44 + a 55 + a 66 + μ) + E * I * S β 2 1 ((a 66 + μ)(a 33 + a 44 + a 55 ) + a 33 (a 44 + a 55 ) + a 66 μ) + (a 66 + μ)((a 44 + a 55 )(a 11 (a 22 + a 33 ) + a 22 a 33 ) + a 44 a 55 (a 11 + a 22 + a 33 ) + a 11 a 22 a 33 ) + (a 44 + a 55 )(a 66 μ(a 11 + a 22 + a 33 ) + a 11 a 22 a 33 ) + (a 44 a 55 + a 66 μ)(a 11 (a 22 + a 33 ) + a 22 a 33 ) + a 44 a 55 a 66 μ b 4 = E * (F * β 2 2 a 22 + I * S β 2 1 a 33 ) + E * F * β 1 β 2 β 3 + a 11 a 22 a 33 ) + a 66 μ(a 22 + a 33 + a 44 + a 55 ) + a 55 (a 66 + μ)(a 11 + a 22 + a 33 + a 44 ) + a 44 (a 55 + a 66 + μ)(a 11 + a 22 + a 33 ) + (a 44 + a 55 + a 66 + μ)(E * (F * β 2 2 + I * S β 2 1 ) + a 11 (a 22 + a 33 ) + a 22 a 33 ) b 5 = E * (F * β 2 2 + I * S β 2 1 ) + (a 11 + μ)(a 22 + a 33 + a 44 + a 55 + a 66 ) + a 55 a 66 + (a 55 + a 66 )(a 22 + a 33 + a 44 ) + a 44 (a 22 + a 33 ) + a 22 a 33 b 6 = a 11 + a 22 + a 33 + a 44 + a 55 + a 66 + μ Here, all the eigen values are negative if a 11 > 0, a 22 > 0, a 33 > 0, a 44 > 0, a 55 > 0, a 66 > 0 i.e. E * < max β 4 +β 5 +μ β 1 , β 3 +β 6 +β 7 +μ β 2 . Then by Routh-Hurwitz criterion we say the endemic equilibrium point is locally asymptotically stable. The novel corona virus is spread through human contact with infected individuals. Therefore, one can put control on respective situation to prevent its spreading. Control description: u 1 : To prevent exposed foreign individuals in the interstate u 2 : Exposed interstate individuals should be quarantined u 3 : Exposed foreign individuals should be quarantined u 4 : Infected individuals should be quarantined The objective function is, where, Λ denotes set of all compartmental variables, Related Langrangian function is given by, The adjoint equation variables, λ i = (λ 1 , λ 2 , λ 3 , λ 4 , λ 5 , λ 6 , λ 7 ) for the system is calculated by taking partial derivatives of the Langrangian function with respect to each compartment variable. This calculation leads with resulting conditions as (Pontryagin, 1986) [8] , Based on analytical results, numerical simulation is given in next section. In this section we discuss the simulation performed for the system (1) From Fig. 2 , we observe 30% of interstate population is exposed to COVID-19 in 17.2 days. Whereas 24.55% of foreigner's population is exposed to COVID-19 in 21.8 days. 21% of foreigners come in contact with Interstate individuals in 7.1 days which increases the infectives of interstate to 22.78% in 9.2 days. Also 27.59% of interstate population gets hospitalized in 14.5 days. Scatter plotting is shown in Fig. 3 . Combined effect of group of three compartments is revealed in each plot. Fig. 3 (a) depicts that; more infected interstate and foreign individuals will be hospitalised at higher rate of level. Fig. 3(b) shows that, individuals who travelled more will be quarantined. From Fig. 3(c) , one can say that infected foreigner would be quarantined at higher rate. Infectedness in quarantined individuals increases which leads to the hospitalization of individuals as observed in Fig. 3(d). Fig. 3 (e) describes that how interstate infected individuals are quarantined. Fig. 4 shows the periodic nature of the interstate class exposed to the virus COVID-19. It indicates that interstate population is exposed again and again to the disease. It happens if the lockdown, social distancing is not followed as per government system. Which shows the importance of the government action taken against COVID-19 to protect the population. Fig. 5 shows the stability of the respective compartments at endemic equilibrium point. Since the government has decided to quarantine foreigners as soon as they arrive in their countries this makes the system stable as they are not exposed much to the COVID-19. 6a, 6b shows the trajectory at the endemic equilibrium point for the system (1). Here we observe the importance of quarantine as the system is stable when interstate and foreigners are quarantined. From Fig. 7 , we observe foreigner are moving towards interstate population. Fig. 8 and Fig. 9 illustrates the flow of interstate population and foreigners with the COVDI-19 infection. It shows that the interstate population gets the infection at a slower rate as compared to foreigners. Fig. 10 shows the flow of interstate and foreigners towards hospitalization. Foreigners gets hospitalized at faster rate than interstate population. Fig. 11 (a)-(g) show the oscillating behaviour of each compartment. As the epidemic nature of disease increases, this can oscillate the whole situation. In some intervals of data, exposed individuals increase ( Fig. 11(a) ) who are either interstate (Fig. 11(b) ) or foreigner (Fig. 11(c) ). If quarantined individuals do Fig. 4 . We plot phase diagram of interstate class when exposed to COVID-19 observing again and again exposure of interstate population to COVID-19. not follow quarantines rules which have been observed in Fig. 11(d) . This leads to a greater number of infected individuals (Fig. 11(e) ) hence they should be hospitalised (Fig. 11(f) ) which effects on recovery rate (Fig. 11(f) ). (1) is observed on each compartment. Here it can be seen that after control is applied, population of each compartment decreases. The above oscillating nature of the model is controlled by the Fig. 12 (a)-(g). All the four controls are effective to our system (1). In the presence of all the controls we observe decrease in the number of exposed individuals. Quarantining interstate and foreign individuals also reduce the infection when controls are applied. Fig. 14 a, we observe out of 26% of interstate population 17% is infected and from Fig. 14 b , among 19% of foreigners we have 18% infected population. The Fig. 13 , shows intensity of mortality due to COVID-19 among interstate and foreign travellers. The Fig. 14 clearly shows the infected population of foreigner is more than that of interstate population which shows the importance of complete ban on air arrivals. From the Fig. 15 it can be observed that 7% of population is exposed to COVID-19. Interstate population share the largest percentage. 14% of the population is quarantined including foreigners and interstate population. Similarly, the infection is 13%. The hospitalization is done at 15%. Of the total population recovery is 16%. Our model indicates that foreigners exhibit a larger infected population, hospitalization rate, and infection rate when compared to the interstate population. Moreover, as foreign individuals contact interstate individuals, the rate of infection within the interstate population increases significantly. To the best of our knowledge, our study is the first to create a disjoint mutually exclusive compartmental model. Our model suggests that both foreign and interstate travel lead to increased risk of infection within the United States population. Consequently, we validate the effectiveness of quarantine as a public health intervention model by the US government and encourage implementation of efforts to mitigate interstate travel. There are multiple reasons that foreigners have increased risk of obtaining COVID-19 when compared to interstate population. First, as people travel, they risk exposing themselves to a greater number of other individuals. The WHO indicates that transmission of COVID-19 occurs primarily through droplet transmission (2019). Most methods of international travel, including airways, railways, and waterways, crowd individuals in compact and enclosed spaces. Being in close contact with individuals with respiratory symptoms in an enclosed environment increases the risk of being exposed to infected mucosae [11] . Second, the guidelines and strategies for addressing the epidemic differ among countries. For example, while India has enforced total lockdown, the US government has not mandated enforced lockdown [3] . Consequently, when individuals from countries with different regulations arrive, they may be infected and increase the incidence of COVID-19. Finally, the vaccination standards differ among countries. In particular, BCG vaccine, believed to confer protective effects against COVID-19 is recommended in some countries, but not the US [7] . As a result, future research and modelling is necessary to determine the protective effects of the BCG vaccine, and its potential to reduce the incidence of COVID-19 within the United States. Given that 2019-nCoV is no longer contained within Wuhan, we recommend the United States government close their borders to both foreign and interstate travel. We recommend significant public health interventions at both international and interstate levels otherwise large cities with close inter-transport systems could become outbreak epicentres. Finally, we recommend preparedness plans and mitigation interventions be readied for quick deployment on both a state and federal level. Based on our model, compliance with these recommendations will effectively reduce the transmission of COVID-19 as a result of foreign and interstate travel. Mathematical models of isolation and quarantine Healthcare impact of COVID-19 epidemic in India: a stochastic mathematical model New SARS-like virus in China triggers alarm The construction of next-generation matrices for compartmental epidemic models First case of 2019 novel coronavirus in the united states Correlation between universal BCG vaccination policy and reduced morbidity and mortality for COVID-19: an epidemiological study Mathematical Theory of Optimal Processes, Routledge A Treatise on the Stability of a Given State of Motion: Particularly Steady Motion Estimation of the transmission risk of the 2019-nCOV and its implication for public health interventions Global transport networks and infectious disease spread Nowcasting and forecasting the potential domestic and international spread of the 2019-ncov outbreak originating in Wuhan, China: a modelling study The first three authors thank DST-FIST file # MSI-097 for technical support to the department. Third author (ENJ) is funded by UGC granted National Fellowship for Other Backward Classes (NFO-2018-19-OBC-GUJ-71790).