key: cord-0908947-w7kzkmf3 authors: Riyapan, Pakwan; Shuaib, Sherif Eneye; Intarasit, Arthit title: A Mathematical Model of COVID-19 Pandemic: A Case Study of Bangkok, Thailand date: 2021-03-30 journal: Comput Math Methods Med DOI: 10.1155/2021/6664483 sha: 196fa068122c405fa730b66fb3685f19f1a849c6 doc_id: 908947 cord_uid: w7kzkmf3 In this study, we propose a new mathematical model and analyze it to understand the transmission dynamics of the COVID-19 pandemic in Bangkok, Thailand. It is divided into seven compartmental classes, namely, susceptible (S), exposed (E), symptomatically infected (I(s)), asymptomatically infected (I(a)), quarantined (Q), recovered (R), and death (D), respectively. The next-generation matrix approach was used to compute the basic reproduction number denoted as R(cvd19) of the proposed model. The results show that the disease-free equilibrium is globally asymptotically stable if R(cvd19) < 1. On the other hand, the global asymptotic stability of the endemic equilibrium occurs if R(cvd19) > 1. The mathematical analysis of the model is supported using numerical simulations. Moreover, the model's analysis and numerical results prove that the consistent use of face masks would go on a long way in reducing the COVID-19 pandemic. The world continues to battle with the coronavirus disease 2019 (COVID-19) caused by the novel coronavirus, SARS-CoV-2, which is regarded as a highly virulent virus that targets the human respiratory system. The pandemic started in late December 2019 with patients admitted to hospitals with an initial diagnosis of pneumonia. The admitted patients' sickness was linked to the seafood and wet animal market in Wuhan, Hubei Province, China [1] . On January 2, 2020, a total number of 41 admitted hospital patients were confirmed to be infected with COVID-19 [2] . On January 22, 2020, 571 COVID-19 cases were reported in 25 different provinces in China [1, 2] . On January 30, 2020, China had about 7734 confirmed COVID-19 cases, while 90 cases were reported in about 13 countries [1, 3] , including Canada, India, Germany, France, the United States, India, and the United Arab Emirates. As of October 31, 2020, a total of 4,667,780 COVID-19 cases (Africa: 1,776,595 cases, Asia: 13,461,293 cases, America: 20,546,580 cases, Europe: 9,840,736 cases, Oceania: 41,880 cases, and others: 696 cases) have been reported worldwide including 1,189,499 deaths (Africa: 42,688 deaths, Asia: 239,675 deaths, America: 640,513 deaths, Europe: 265,565 deaths, Oceania: 1,051 deaths, and others: 7 deaths) [4] . COVID-19 is transmitted from human to human via direct contact with contaminated surfaces and through respiratory droplets' inhalation from infected individuals [5] . Presently, there is no vaccine or antiviral treatment approved for the prevention or management of COVID-19 [6] . To effectively reduce the spread of COVID-19, governments have been implementing various control measures such as imposing strict, mandatory lockdowns and encouraging (and in some cases strictly enforcing) other measures such as individuals maintaining a minimum distance between themselves (social distancing), avoiding crowded events, imposing a maximum number of individuals in any gathering (religious and social), and the use of face masks in public [7] . To further help mitigate the spread of COVID-19, contact tracing of suspected infected cases has been stepped up in several countries and detected cases (asymptomatic and symptomatic) are quickly placed in isolation for prompt treatment [7] . In Thailand, the virus was first confirmed to exist on January 13, 2020 [8] , while the first reported local transmission was confirmed on January 31, 2020 [9] . The number of cases remained low throughout February, but it surged in mid-March. The rise in the number of cases was traced to several transmission clusters, the largest of which occurred with a Muay Thai fight at the military-run Lumpinee Boxing Stadium on March 6, 2020 [10] . Confirmed cases rose to over a hundred per day over the following week, and public venues and businesses were ordered to close in Bangkok and several other provinces [11] . Bangkok businesses' abrupt closure prompted tens of thousands of workers to travel to their hometowns [12] . The Prime Minister of Thailand, Prayut Chan-o-cha, declared a state of emergency, effective on March 26, 2020 [13] , and a curfew went into effect on April 3, 2020 [14] . All commercial international flights were suspended from April 4, and lockdown measures were implemented in varying degrees throughout the country. The rate of new cases gradually dropped throughout April, and by mid-May, locally transmitted infection rates had fallen to near-zero, and easing of restrictions was gradually implemented [15, 16] . Even though there had been no new domestic cases since mid-May, on August 21, Thailand extended its emergency degree until September 30 to prevent incoming aliens from overseas in many routes. In September, a prison inmate who had not been abroad was Thailand's first locally transmitted case in 100 days. Later in the month, Akbar Ismatullaev, a footballer, was infected with the virus after completing the 14-day state quarantine since he arrived nearly a month earlier. In October, foreign tourists entered Thailand for the first time in seven months [17] . A French tourist on Ko Samui in Surat Thani contracted the disease after passing the 14-day state quarantine. She developed a fever 17 days after arriving in the country. In October, foreign tourists entered Thailand for the first time in seven months under the Special Tourist Visa program [18] . Mathematical modeling is a valuable tool to control disease spread effectively. Several useful mathematical models have been formulated in the last few decades to study infectious diseases and develop helpful strategies for the efficient elimination of infection [19] [20] [21] [22] [23] [24] [25] . The compartmental models and real cases are more effective in providing valuable information about a particular disease outbreak. Several mathematical models have already been formulated in various countries to analyze the complex transmission pattern of the COVID-19 pandemic, using ordinary differential equations [26] [27] [28] , delay differential equations [29] , stochastic differential equations [30] , and fractional order Caputo derivative [31] [32] [33] [34] [35] [36] . Thailand has begun opening its borders to foreigners and with the growing cases in the USA and Europe, there is a concern that COVID-19 cases in Thailand may rise soon. Hence, this study is aimed at (i) formulating a mathematical model to understand the dynamics of the COVID-19 pandemic in Bangkok, Thailand, (ii) examining the impact of the control measures currently employed in Thailand, and (iii) determining if these measures will be effective in preventing COVID-19 cases in Thailand. The remainder of this article is structured as follows. The proposed model is presented in Mathematical Model Formulation and Description. The mathematical analysis of the model is presented in Analysis of the Model. The results obtained from numerical simulations of the model are provided in Numerical Simulation. Finally, the conclusion drawn from this study is given in Conclusion. The mathematical model of COVID-19 transmission formulated in this study was motivated by the study of [26] . The model proposed in [26] was constructed from the SEIR model and is comprised of six compartments with the infected compartment divided into three categories. However, in the present study, the model will be divided into seven compartments. The total human population to be considered is denoted as NðtÞ, and at any time, it comprises of the susceptible ðSÞ, exposed ðEÞ, symptomatically infected ð I s Þ, asymptomatically infected ðI a Þ, quarantined ðQÞ, recovered ðRÞ, and death ðDÞ compartments, respectively. The susceptible compartment comprises individuals living in the country or who have recently returned before the border's closure. The individuals exposed to COVID-19 and show signs of symptoms are moved to the symptomatic infectious compartment. In contrast, individuals who show no sign of symptoms are moved to the asymptomatic infectious compartment. There is a reduction in the risk of infection for the individuals in the susceptible compartment since they practice preventive measures such as social distancing, wearing a mask, and refraining from mass gatherings or meetings. Individuals move to the recovery compartment through recovery from both the quarantined and infected compartments, respectively. The parameters β s and β a , respectively, represent the effective contact rate (contacts capable of leading to COVID-19 transmission) for individuals in the symptomatically infectious and asymptomatically infectious compartments. The proportion of individuals who wear face masks correctly within a community is denoted as 0 < ψ ≤ 1 while 0 < ξ ≤ 1 represent the expected decrease in the risk of infection due to the face mask's use. The progression rate of exposed individuals is denoted as φ. A proportion 0 < θ ≤ 1 of exposed individuals showed no clinical symptoms of COVID-19 (and move to the compartment I a ) at the end of the incubation period. The remaining proportion 1 − θ shows clinical symptoms and moves to the I s compartment. The parametersλ s , λ a , and λ q represent the recovery rate for individuals in I s ,I a , and Q compartments, respectively. Similarly, α s and α a are the isolation rate of individuals. Finally, the parameters δ s and δ q represent the COVID-19-induced mortality rate for individuals in the asymptomatic infectious and quarantined compartments, respectively. In this study, the effect of social distancing and community lockdown 2 Computational and Mathematical Methods in Medicine measures will be measured based on the overall reduction in the community contact rate parameter's baseline values. The flowchart of the formulated model using all the above assumptions is given in Figure 1 . Additionally, all variables with their meaning and the parameters with their description are provided in Tables 1 and 2, respectively. From the above assumptions and variables, we formulate the model with nonnegative initial conditions: Since the total population is From (2), NðtÞ would approach a carrying capacity ρ/μ. Model (1) describes the human population, and thus, the model variables can be shown to be nonnegative for all time t ≥ 0 and that all solutions of the model (1) will remain pos-itive for all time t ≥ 0. Therefore, model (1) is mathematically well-posed, and its dynamics can be considered in the region below: In this section, the dynamical properties of model (1) are qualitatively studied. By using the technique of variable separation, dS/dt can be reduced to Then, the above equation is integrated to yield the solution below Since the initial value S 0 and the exponential functions in Equation (6) Proof. The addition of all equations in Equation (1) yields which can be rewritten as From Equation (8), NðtÞ approaches ρ/μ as t → ∞: Hence, the positive solutions of model (1) are bounded. The equilibrium points of model (1) are obtained by zeroing the right-hand side of all equations in model (1), resulting in The simplification of Equation (9) can yield many solutions. However, in this study, we consider two solutions: the disease-free equilibrium point (DFEP) and the endemic equilibrium point (EEP). The DFEP in this study is denoted as Furthermore, the EEP is denoted by The effective contact rate β a The effectiveness of social distancing ψ The proportion of individuals who use a face mask ξ The efficacy of face masks 1 − θ The fraction of exposed individuals who show clinical symptoms after the incubation period φ The rate of progression from the exposed compartment to the infectious compartment (1). To compute the model's basic reproduction number (BRN), the next-generation matrix approach is employed. In this study, we denote the BRN of model (1) as R cvd19 , which is defined as the number of secondary cases of COVID-19 infection arising from one individual infected with the COVID-19 disease. By using the notation in the study of [37] , the vectors F (denotes new infection) and V (transfer of individuals between compartments) are given as follows: From vectors F and V , the Jacobian F and V are computed below: The Jacobian matrices F and V evaluated at χ * eep yield: Hence, computing the Jacobian matrices F and V, the BRN of the model is where S * dfep is denoted in Equation (10). From Equation (15), the BRN of model (1) can be further simplified to Equation (16): Equation (16) comprises two reproduction numbers. The first BRN is ðβ s φS * dfep ð1 − ψξÞð1 − θÞÞ/ðφ + μðα s + δ s + λ s + μÞÞ, and it defines the number of new COVID-19 cases generated from symptomatically infectious human in compartment I s : The second BRN is ðθβ a φS * dfep ð1 − ψξÞÞ/ðφ + μ ðα a + λ a + μÞÞ, and it defines the number of new COVID-19 cases generated from asymptomatically infectious humans in compartment I a . Hence, mathematically, the BRN can be simplified to be To prove the global stability of the disease-free equilibrium, the Lyapunov function below is constructed: where From Equation (18) , the derivative is given below as The expansion of Equation (20) yields Then, Equation (21) becomes Equation (22) can be simplified to Hence, ∂L/∂t ≤ 0 if R cvd19 ≤ 1, and ∂L/∂t = 0 if E = 0: By LaSalle's Invariance Principle, we can conclude that the DFEP of model (1) is globally asymptotically stable in Ω Cvd19 whenever R cvd19 ≤ 1. The global asymptotic stability of χ * eep is discussed using the Lyapunov asymptotic stability theorem. From model (1), we will construct a Lyapunov function by following the study of Xu et al. [38] . (1) is given by the expression below: Computational and Mathematical Methods in Medicine From (2), all the solutions of (11) satisfy the equality Also, N = e −ut+C + ρ/μ, where C is the value that satisfies the condition N 0 = ρ/μ. Thus, L = ln ðN − N * + 1Þ ≥ 0. Therefore, L is a Lyapunov function for model (1) , and the endemic equilibrium χ * eep is globally asymptotically stable by the Lyapunov asymptotic stability theorem. To support the mathematical analysis of model (1), the numerical simulations are carried out using the deSolve package [39] with the fourth-order Runge-Kutta method in RStudio programming software version 1.1.442. The simulations are divided into three parts. Part 1 is to illustrate the numerical interpretation of the disease-free and endemic equilibrium points. Part 2 is to explore the varying effects of face masks. Part 3 is to find out model (1) fitting with real data. The parameter values used for the numerical simulations in Part 1 are provided in Table 3 . The results obtained from the numerical simulations are presented in Figures 2 and 3 , respectively. Figure 2 depicts a rise in the number of the susceptible population in the absence of the COVID-19 pandemic. We also note that when there is a rise in the susceptible population, the exposed population also increases. However, the increase in the exposed population is lesser compared to the susceptible population. Furthermore, no changes occurred in the quarantined and recovered population. This can be because in the absence of the COVID-19 pandemic (i.e., I s = 0 and I a = 0), there will be no infected population to spread the disease to the other compartments. Figure 3 shows a slight rise in the number of individuals in the exposed compartment and a gradual reduction in the number of the infected population when necessary interventions are used as a preventive measure to reduce the spread of the COVID-19 pandemic. Additionally, there is a rise in the number of quarantined and recovered individuals from the COVID-19 pandemic. The figure also finds that if the interventions are strictly followed, it can reduce the pandemic's spread. The Varying Effects in the Use of Face Mask. The numerical simulation in Part 2 examines how varying ψ value affects the basic reproduction number computed in this study. As earlier defined in Table 2 , ψ denotes the proportion of individuals who use face masks. The parameter values used in the numerical simulations for Part 2 are provided in Table 3 . However, different ψ values (ψ = 0:1, 0:5, and 0:7) are considered. By computing the basic reproduction numbers and using the parameter values in Table 3 , the obtained results are presented in Table 4 . From the results obtained in Table 4 , it is evident that if a larger number of people consistently use face masks in a community, then the COVID-19 pandemic can be reduced. In this section, the real data provided in Table 5 is fitted with model (1). The data were grouped monthly and obtained from the Department of Disease Control, Thailand [47] . The data span from January 2020 to December 2020. The monthly recorded data in Table 5 [47] were first interpolated into daily data to easily fit the real data with the simulated data. Afterward, model (1) was fitted to the daily interpolated data using a step size of 0.01. The result of the model fit is presented in Figure 4 . Table 3 when R cvd19 = 0:28738 < 1. Table 3 when R cvd19 = 1:40995 > 1. From Figure 4 , it can be seen that the real data have two peaks which occurred on days 360 and 1,073. The model proposed in this study was able to generate the peaks obtained in the real data. In this study, the nonlinear mathematical model was proposed and analyzed to understand the dynamics of the COVID-19 pandemic in Thailand. The equilibrium point relating to the formulated model was computed. Using the next generation matrix approach, the basic reproduction number denoted as relating to the model was also computed. Moreover, this study also showed that if the BRN is denoted as R cvd19 < 1, then the pandemic will die out. However, if R cvd19 > 1, then the pandemic will remain in the population. Additionally, the global asymptotical stability of the disease-free and endemic equilibrium points has been proved. Numerical simulations were carried out to support the model analysis. The real data were also fitted to the model for predicting the infected population cases in real life. The varying effects of the use of face masks were also explored in this study, and it was found that the continuous and appropriate use of face masks can prevent the spread of the COVID-19 pandemic. Presently, the research on a vaccine to prevent the COVID-19 pandemic has yielded excellent results, with Pfizer announcing that their vaccine has a 95% efficacy. However, it will take a while before the vaccines are made readily available in all countries worldwide. Therefore, the use of face masks should be made compulsory till the vaccines are available for everyone. We propose that future researchers implement the model proposed in this study to the second wave of infected cases in Thailand to explore the efficiency of the current measures used to prevent COVID-19. The data used to support the findings of this study are included within the article. The epidemiology and pathogenesis of coronavirus disease (COVID-19) outbreak Drug treatment options for the 2019-new coronavirus (2019-nCoV) The novel Chinese coronavirus (2019-nCoV) infections: challenges for fighting the storm The European center for disease prevention and control Presumed asymptomatic carrier transmission of COVID-19 An updated estimation of the risk of transmission of the novel coronavirus (2019-nCov) Coronavirus: The World in Lockdown in Maps and Charts COVID-19 Dashboard USA, COVID-19 Dashboard Thailand Confirms First Case of Wuhan Virus outside China, South China Morning Post First Thai Case with No Records of Travelling to China, Thairath Online COVID-19 cases in Thailand surge; spreads from four major Bangkok-area entertainment districts -boxing stadiums to 9 provinces Bangkok and five other provinces close malls and 25 other locations (with clip) Nearly 100,000 flock home on tour buses to flee COVID! Transport Minister orders 7-step procedure PM declares state of emergency effective March 26 in response to COVID-19, maintains there is no curfew yetThe Standard (in Thai) Curfew starts today Thai Traffic Back to Gridlock as Coronavirus Measures Ease, Reuters 39 Chinese Are Thailand's 1st Foreign Tourists in 7 Months Mathematical analysis of a tuberculosis model with imperfect vaccine A contribution to the mathematical theory of epidemics A mathematical model for the population dynamics of malaria with a temperature dependent control Backward bifurcation of an epidemiological model with saturated incidence, isolation and treatment functions Analysis of a co-infection model for HPV-TB A mathematical study of a model for HPV with two high-risk strains SIR epidemic model with Mittag-Leffler fractional derivative Mathematical Modeling and Analysis of COVID-19 Pandemic in Nigeria, medRxiv Modeling COVID-19: Forecasting and Analyzing the Dynamics of the Outbreak in Hubei and Turkey, medRxiv Mathematical modeling for novel coronavirus (COVID-19) and control Modeling the epidemic dynamics and control of COVID-19 outbreak in China A fractional-order SEIHDR model for COVID-19 with inter-city networked coupling effects Modelling the spread of COVID-19 with new fractal-fractional operators: can the lockdown save mankind before vaccination? Modeling the dynamics of novel coronavirus (2019-nCov) with fractional derivative Modeling and simulation of the novel coronavirus in Caputo derivative The dynamics of COVID-19 with quarantined and isolation Dynamics of fractional order COVID-19 model with a case study of Saudi Arabia Modeling and analysis of the dynamics of novel coronavirus (COVID-19) with Caputo fractional derivative On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations Global stability of a variation epidemic spreading model on complex networks Solving differential equations in R: package Desolve Mathematical Modeling and Analysis of COVID-19 Pandemic in Nigeria, medRxiv Impact of Non-Pharmaceutical Interventions (NPIs) to Reduce COVID-19 Mortality and Healthcare Demand, Imperial College COVID-19, Response Team Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV-2) Estimates of the severity of coronavirus disease 2019: a model-based analysis Public health responses to covid-19 outbreaks on cruise ships-worldwide To mask or not to mask: modeling the potential for face mask use by the general public to curtail the COVID-19 pandemic Impact of Non-Pharmaceutical Interventions (Npis) to Reduce COVID19 Mortality and Healthcare Demand The authors declare no conflict of interest.