
Citation: Mario Santana-Cibrian, Manuel A. Acuña-Zegarra, Jorge X. Velasco-Hernandez. Lifting mobility restrictions and the effect of superspreading events on the short-term dynamics of COVID-19[J]. Mathematical Biosciences and Engineering, 2020, 17(5): 6240-6258. doi: 10.3934/mbe.2020330
[1] | Fen-fen Zhang, Zhen Jin . Effect of travel restrictions, contact tracing and vaccination on control of emerging infectious diseases: transmission of COVID-19 as a case study. Mathematical Biosciences and Engineering, 2022, 19(3): 3177-3201. doi: 10.3934/mbe.2022147 |
[2] | Boqiang Chen, Zhizhou Zhu, Qiong Li, Daihai He . Resurgence of different influenza types in China and the US in 2021. Mathematical Biosciences and Engineering, 2023, 20(4): 6327-6333. doi: 10.3934/mbe.2023273 |
[3] | Wenhan Guo, Yixin Xie, Alan E Lopez-Hernandez, Shengjie Sun, Lin Li . Electrostatic features for nucleocapsid proteins of SARS-CoV and SARS-CoV-2. Mathematical Biosciences and Engineering, 2021, 18(3): 2372-2383. doi: 10.3934/mbe.2021120 |
[4] | Sarafa A. Iyaniwura, Rabiu Musa, Jude D. Kong . A generalized distributed delay model of COVID-19: An endemic model with immunity waning. Mathematical Biosciences and Engineering, 2023, 20(3): 5379-5412. doi: 10.3934/mbe.2023249 |
[5] | Chenxi Dai, Jing Yang, Kaifa Wang . Evaluation of prevention and control interventions and its impact on the epidemic of coronavirus disease 2019 in Chongqing and Guizhou Provinces. Mathematical Biosciences and Engineering, 2020, 17(4): 2781-2791. doi: 10.3934/mbe.2020152 |
[6] | Abdelkarim Lamghari, Dramane Sam Idris Kanté, Aissam Jebrane, Abdelilah Hakim . Modeling the impact of distancing measures on infectious disease spread: a case study of COVID-19 in the Moroccan population. Mathematical Biosciences and Engineering, 2024, 21(3): 4370-4396. doi: 10.3934/mbe.2024193 |
[7] | A. M. Elaiw, Raghad S. Alsulami, A. D. Hobiny . Global dynamics of IAV/SARS-CoV-2 coinfection model with eclipse phase and antibody immunity. Mathematical Biosciences and Engineering, 2023, 20(2): 3873-3917. doi: 10.3934/mbe.2023182 |
[8] | Anichur Rahman, Muaz Rahman, Dipanjali Kundu, Md Razaul Karim, Shahab S. Band, Mehdi Sookhak . Study on IoT for SARS-CoV-2 with healthcare: present and future perspective. Mathematical Biosciences and Engineering, 2021, 18(6): 9697-9726. doi: 10.3934/mbe.2021475 |
[9] | Rahat Zarin, Usa Wannasingha Humphries, Amir Khan, Aeshah A. Raezah . Computational modeling of fractional COVID-19 model by Haar wavelet collocation Methods with real data. Mathematical Biosciences and Engineering, 2023, 20(6): 11281-11312. doi: 10.3934/mbe.2023500 |
[10] | Fang Wang, Lianying Cao, Xiaoji Song . Mathematical modeling of mutated COVID-19 transmission with quarantine, isolation and vaccination. Mathematical Biosciences and Engineering, 2022, 19(8): 8035-8056. doi: 10.3934/mbe.2022376 |
Latin America has recently become the world's epicenter for SARS-CoV-2 activity. Lack of infrastructure, weak economic performance and social inequality can, in part, explain this situation. In many countries such as Mexico, a large percentage of the economically active population works in informal jobs that largely depend on working in the street or plazas of towns and cities. This explains the need to reopen the economy in these countries but, on the other hand, ineffectively justifies the lifting of mobility restriction policies in the midst of the epidemic. However, it is not only the economically deprived countries of the American continent that are suffering the blunt of the COVID-19 pandemic. The United States has been unable to bring the spread of the disease under control. As has been the case in many Latin American countries, the USA, in general, reopened its economy too fast and without a clear strategy in place. This has brought a worrisome rebound of the epidemic that threatens the very way of life of this rich and powerful country. It is within this context that we present our work. There is the need for strategic modeling, geared to local, specific situations, that can provide projections to understand and evaluate the consequences of lifting mitigation measures that were put into place earlier this year, to fight the so-called first epidemic wave at the local (county or city) levels. Mathematical modeling has been a valuable tool to understand the dynamics of COVID-19 around the world. There is a large number of models published both in journals and public repositories. Several models share similarities that arise from the fact that these are based on the Kermack-McKendrick equations and/or have similar objectives (e.g., analysis of the lifting of mobility restrictions). We briefly mention here those that have direct relation to our approach, namely, the consideration of a decreasing effective contact rate during lockdown, and an explicit modelling of population classes induced by the nonpharmaceutical interventions (NPIs). Dehning et al. [1] propose a model with several change points in the contact rate, as we do in this work, and show that, in the case of Germany, this changes can be correlated to dates of publicly announced interventions. Sardar et al. [2] present other extension including lockdown and hospitalized classes and suggest that a high percentage of symptomatic individuals increases lockdown efficacy. Bertozzi et al. [3] presents three simple models to describe the early evolution of the COVID-19 pandemic: exponential growth, a self-exciting branching process, and the SIR compartmental model. The authors show scenarios where relaxing social distancing and non-pharmaceutical interventions results in the resurgence of the epidemic. However, neither of the two previous references include explicitly the population splitting into two different but inter-related systems as we do here. With a different approach, Badr et al. [4] look at the relationship between changes in mobility and the rate of new infections. It is shown that, in the absence of pharmaceutical measures, social-distancing is a good strategy against COVID-19 spread. Nonetheless, these authors do not explicitly explore the enhanced superspreading of the disease that the higher mobility induces nor its effects on the shape of the epidemic curve as we do here. The importance of NPIs in the mitigation of the epidemic has been explored in [5] where it is shown that increases in the compliance rate of facemask use, significantly decrease the effective contact rate.
In this paper, we explore the effect of atypical superspreading events and the effect of lifting mitigation measures on the shape of the epidemic curve. We particularly concentrate on the impact of what we call atypical mobility events that are and were associated with population events where some level of disease superspreading occurred. These short-duration, superspreading events have a non-negligible effect upon the shape and time evolution of the epidemic curve. In many places of the world, such as in the USA and, at present, in several regions of Mexico, we are seeing that the epidemic has entered a plateau, and when the plateau ends it gives way to a new epidemic growth phase. We apply our model to the case of Mexico City. We analyze the interrelation between the date of maximum incidence and heightened incidence events taking place before the projected maximum. The manuscript is organized as follows. In Section 2, we formulate the mathematical model and give some general results arising from it. In Section 3, we describe the Mexico City data that is used to exemplify our model results and present a brief analysis of the impact of mitigation measures. In Section 4, we use our proposed model to generate scenarios that describe the impact of superspreading events associated with holidays on the shape of the epidemic curve and explore the probable consequences of lifting mitigation measures. Finally, the discussion and conclusions are presented in Section 5.
Our model is an extension of the mathematical model first published in [6], and further developed in [7]. We extend it as follows: (ⅰ) we introduce a compartment for the reported infectious cases, (ⅱ) we introduce two time-dependent contact rates to approximate their reduction after the initial mitigation measures were implemented, and (ⅲ) we introduce a time-dependent confinement-compliance rate. These extensions allow us to provide a better account of the social-distancing effects and the role and impact of several short periods of increased transmission due to the lack of compliance of mobility restriction guidelines during the confinement.
Our mathematical model explicitly takes into account two time intervals, before and after the date when mitigation measures were implemented, T, as illustrated in Figure 1. The model is of SEIR-type with three infectious classes: Ia asymptomatically, Is symptomatically, and Ir reported (confirmed) infectious individuals. The dynamics for the first time interval are described by the following system of ordinary differential equations:
S′=−b(Is+Ia)SN∗,E′=b(Is+Ia)SN∗−γE,Ia′=ργE−ηaIa,Is′=(1−ρ)γE−(ηs+δs)Is,Ir′=δsIs−δrIr,R′=ηaIa+ηsIs+(1−μ)δrIr,D′=μδrIr, | (2.1) |
where N∗=S+E+Ia+Is+R. Note that the Ir compartment does not participate in transmission since, it is assumed, once confirmed, the cases are effectively isolated. b represents the effective contact rate, 1/γ the incubation period, ρ the proportion of asymptomatically infected individuals, 1/ηa and 1/ηs the periods from symptoms onset to recovery for symptomatically and asymptomatically infected individuals, δs the rate at which a symptomatically infected individuals becomes a reported (confirmed) case, 1/δr the time from confirmation to recovery or death and μ is the proportion of those reported cases that die. Currently, the role of asymptomatic individuals in the transmission dynamics is not yet clear. The contact rate is related to the probability of infection and to the number of contacts between infected and susceptible individuals. There is evidence that the viral load is similar for both asymptomatically and symptomatically infected individuals [8,9] which makes it reasonable to assume that the probabilities of infection by a symptomatic or by an asymptomatic individual are similar although not necessarily equal. However, since we assume that symptomatic and asymptomatic individuals have the same mobility, it is until the symptomatic cases are reported (Ir) that they start confinement and no longer participate in the transmission process. We assume that the proportion of symptomatic cases that is not reported (i.e., escape detection) has mild symptoms and behaves as an asymptomatic individual. This justifies our hypothesis that Ia and Is have the same effective contact rate b. This simplification helps to reduce the number of parameters in the model.
The basic reproductive number of system (2.1) is given by:
R0=bρηa+b(1−ρ)ηs, |
where the first term represents the number of new cases produced by asymptomatically infected individuals and the second term represents the number produced by symptomatically infected individuals. After mitigation measures are implemented at time T, the population is split into two groups: one constituted by those individuals who comply with the measures that is referred to as the confined group, and another constituted by those individuals who do not, that is called the unconfined group, either because they disobey the social-distancing guidelines or because they belong to an strategic sector of the economy. The dynamics in the second time interval are given by:
S′i=−βi(t)(Isi+Iai)SiN∗+(−1)iω(t)S1E′i=βi(t)(Isi+Iai)SiN∗−γEi+(−1)iω(t)E1Iai′=ργEi−ηaIai+(−1)iω(t)Ia1Isi′=(1−ρ)γEi−(ηs+δs)Isi+(−1)iω(t)Is1Ir′=δs(Is1+Is2)−δrIrR′=ηa(Ia1+Ia2)+ηs(Is1+Is2)+(1−μ)δrIrD′=μδrIr | (2.2) |
where N∗=∑2i=1(Si+Ei+Iai+Isi)+R. The index i=1 gives the population that complies with the control measures, while i=2 indicates those who do not. Note that the compartments Ir, R, and D are common to both groups. Following [6], the effective contact rates are different for each group and are given by:
βi(t)={b−(1−q1i)θ1b(t−T),T≤t<T+θ1,q1ib−(q1i−q2i)θ2b(t−T−θ1),T+θ1≤t<T+θ1+θ2,q2ib,T+θ1+θ2≤t. | (2.3) |
Equation (2.3) assumes decreasing contact rates where the parameters q2i≤q1i, i=1,2, represent second and first reductions of the baseline effective contact rate b. T is the time at which the effective contact rates start to act and θj, j=1,2, are the duration of the contact reduction interval. See Figure 2.
Finally, we define the compliance-failure rate ω(t) as a step function of the form
ω(t)={κω0,T∗≤t≤T∗+θω,ω0,otherwise, | (2.4) |
where κ>1 is the increase in the compliance-failure rate relative to the baseline value ω0 (κ=1), T∗ is the time at which ω, the perturbation, starts and θω is the duration of the perturbation. This model admits several periods during which ω(t) may act, but here we illustrate our results with only one period.
The inclusion of short-term superspreading events into Eqs 2.3 and 2.4 renders scenarios in which plateau-like behavior appears. We show three scenarios: (Ⅰ) long plateau with slight decrease and rebound; (Ⅱ) shorter plateau with a later decreasing trend in the epidemic curve; (Ⅲ) long plateau with a later decreasing epidemic curve. Each scenario is subdivided into three sub-cases defined as follows:
● Scenario X.a: Increase κ=3 of the baseline value ω0.
● Scenario X.b: Increase κ=5 of the baseline value ω0.
● Scenario X.c: Increase κ=7 of the baseline value ω0.
where X=I,II,III. Fixed parameter values are given in Table 1. The length of the interval during which the compliance-failure rate increases is 11 days (Eq 2.4).
Parameter | Value | Parameter | Value | Parameter | Value |
b | 0.887677 | γ−1 | 5.1 | T | 30 |
ρ | 0.698725 | η−1a | 6 | q11 | 0.6 |
μ | 0.11 | η−1s | 7 | θ1 | 18 |
ω | 0.005 | δ−1s | 4 | q12 | 0.7 |
q | 0.7 | δ−1r | 12.070691 | T∗ | 67 |
To qualitatively illustrate the model projections, we set q1i as in Table 1, and in (2.3), we set, for both groups, θ1 to be approximately two weeks and θ2 to be about 3 months with q21=0.3 (decrease of 70% for the confined group), q22=0.4 (decrease of 60% for the unconfined group).
Figure 3 illustrates the case where the perturbed epidemic curve shows a slight decrease immediately after the perturbation relative to the baseline curve. Here, it is observed that as κ increases, the incidence increases again. Ic is the worst-case scenario, where the growth rate keeps increasing. Scenario Ia illustrates the appearance of plateau-like behavior followed first by a slight decline of incidence and then growth again. The length of the plateau sets in for several weeks.
Setting q1i as in Table 1, in (2.3), we set, as before, θ1 to be about two weeks and θ2 to be about three months with q21=0.3 (decrease of 70% for the confined group), q22=0.3 (decrease of 70% for the unconfined group).
Figure 4 illustrates the case where the epidemic curve decreases after the end of the interval of maximum incidence. IIc is the worst-case scenario, where, after the perturbation, the epidemic peak is pushed higher and later in time. Scenario IIa, on the other hand, is comparatively benign. The peak is still reached as projected on the baseline case, and then the epidemic curve decreases at a slower rate than the baseline. Finally, scenario IIb shows the appearance of plateau-like behavior. In this case, when the peak is reached, the incidence curve does not show a significant decay but, rather, enters a sustained phase of maximum incidence that lasts several weeks.
Set q1i as in Table 1, and, in (2.3), we set, for both groups, θ1 to be about two weeks but now θ2 to be about four months with q21=0.2 (reduction of 20% in the confined group), q22=0.3 (reduction of 30% in the unconfined group).
Figure 5 illustrates the case where the epidemic curve decreases after the peak. Scenario IIIa shows the appearance of a plateau that lasts some months after the date of the baseline peak. On the other hand, scenarios IIIb and IIIc show a runaway epidemic with a much higher peak than that of the baseline curve.
Remark. Depending upon the scenario, parameter values θj and qij can be modified to explore it, including the effect of varying effective contact rates (Eq 2.3).
In this section, we illustrate our model results with the case study of Mexico City. We start the analysis with some context. On March 23, 2020, social-distancing measures were implemented in Mexico to slow the spread of the COVID-19 pandemic, mainly focused on closing schools and some non-essential activities. One week later, on March 30, 2020, a Sanitary Emergency was declared to last until April 30, 2020, a date that was later extended to June 01, 2020. These measures aimed to “flatten the curve”, meaning to lower the incidence to ensure that critical cases would remain under manageable levels. On April 16, 2020, the federal government announced that, for Mexico City, the day of maximum incidence would occur on May 08–10, 2020 [10,11]; consequently, the date for lifting mitigation measures was announced to start on June 01, 2020. However, the peak did not happen on the projected dates. Independent projections by Mexican scientists [6,12] and later, the factual evidence from data, set the dates of maximum incidence to be in an interval of about ten days around May 30, 2020, which puts the lifting of mobility restrictions indicated by the Federal government precisely in the likely days of higher transmission. Confronted with this scenario, the local Mexico City government quickly put into place a set of specific actions that have so far successfully kept the epidemic in the city within still manageable levels. Data shows that, by the end of May, the epidemic in Mexico City seemed to stabilize on a plateau that continued beyond June up to at least mid July. This plateau may have appeared as a consequence of two important holidays (children’s day and mother’s day) occurring approximately 10–14 days before the expected date for the peak of maximum incidence. We postulate that these events are an important factor that explains the observed quasi-stationary epidemic trend that the data shows at present (mid-July 2020). In this section, we attempt an explanation of the probable causes of such behavior. Furthermore, we explore the possible consequences of lifting the mitigation measures under diverse conditions.
COVID-19 epidemic data was provided by the Secretaría de Ciencia, Tecnología e Innovación of the Government of Mexico City through the COVID-19-CDMX database [13]. This data set contains details on all the confirmed and suspected COVID-19 individuals such as sex, age, residence place, date of symptoms onset, etc. We use records from February 22, 2020, to July 19, 2020.
The Richards model is an extension of a simple logistic growth model that is a standard tool commonly used to predict cumulative COVID-19 cases in China (see, for example, [14]). In this model, the curve of cumulative cases, C(t), is described by the solution of
C′(t)=rC(t)[1−(C(t)K)a], | (3.1) |
where r is the epidemic growth rate, K is the final epidemic size, and a is a parameter that accounts for the asymmetry of the epidemic curve. We estimate the Richards curve for two different periods: the first goes from February 22, 2020 to March 22, 2020, and the second from March 23, 2020 to April 30, 2020. The parameters a, r, and K for each period are estimated using Bayesian inference [15,16]. Technical details can be found in Appendix A. We limit the fitting of the Richards curve to April 30, since on that date an important event of increased transmission started that violated the hypothesis that population conditions for transmission were essentially unchanged since March 23, 2020. In terms of our model, unchanged conditions mean ω(t)=ω0 for all t.
Before the start of social distancing on March 23, 2020, the growth rate was approximately 0.76, with a 95% interval (0.28, 1.91), which indicates a very fast growing epidemic with a doubling time of about a day. After March 23, 2020, the median growth rate was 0.12, with a 95% interval of (0.11, 0.21), rendering a median doubling time of about 6 days. This gives us a median reduction of 84% of the epidemic growth rate in the early days of the implementation of social-distancing measures. The uncertainty around this value is significant, but our estimate gives evidence that the confinement strategy was indeed effective in reducing transmission during the first weeks of mitigation. This justifies the idea of including decreasing effective contact rates into our modeling framework. Unfortunately, however, those same social-distancing measures also pushed the interval of the likely occurrence of the epidemic peak towards the end of May, thus overlapping it with the date set to reopen the economy.
Remark: The interval of dates for maximum incidence is in agreement in three models. Two of them, published in early May [6,12], projected the peak within an interval from late May to early June; the third model is the projection of the Richards curve (see Figure 6).
To evaluate the impact of mobility events, we set up Eqs 2.1–2.2 under the scenario that 70% of the general population is confined (obeying mitigation restrictions) during the confinement period (second time interval of our model) and 30% is not [17]. The unconfined population includes workers with essential economic activities in government and industry and individuals who work in the informal economy. The time where confinement started, T, is March 23, 2020. Roughly following the trend of the instantaneous reproduction number Rt [18,19], the duration of the first reduction in the effective contact rate is θ1=18 days (Figure 7). The magnitudes of the first reduction in both effective contact rates are q11=0.6 (40% reduction in the confined group) and q12=0.7 (30% reduction in the unconfined group), respectively. The duration θ2 and magnitudes of the second reductions (q2i) depend on the scenario under consideration as described in the next section.
Once under confinement, individuals may abandon the isolation with a compliance-failure rate ω(t). This is the parameter we use as a proxy for population mobility that triggers superspreading events. ω(t) is a time-dependent rate: we assume that increased mobility lasts only for a period of θω days, with a background compliance-failure rate ω0 that in atypical events is increased by a factor κ (Eq 2.4). We take ω0=0.005/day as a baseline value [6]; that is, given the exponential nature of this rate in the baseline case, 50% of the confined populations will stay there at least for 4 months, while the other 50% may abandon confinement faster.
Non-pharmaceutical interventions (NPIs) have been used to reduce contact between individuals in many countries. However, short-term increases in population mobility have occurred within the confinement period. This increased mobility weakens the strength of the NPIs and, therefore, has an impact on disease transmission. Here we explore the consequences of increasing mobility to give a plausible explanation for the appearance of the plateau in the dynamics observed in several countries using, as an example, the case of Mexico City.
In Mexico, there were two important holidays (in terms of population mobility) within the period of confinement: April 30, 2020 (children’s day) and May 10, 2020 (mother’s day). Population mobility increased these days as evidenced by the instantaneous reproduction number (Figure 7). Note that Rt shows a slight increase just around April 30, 2020, and May 01, 2020. We center our attention on the effect of increased mobility within the period from April 29, 2020, to May 10, 2020.
Mexico City epidemic has shown a plateau-like behavior since middle-May (Figure 8). As Figure 7 shows, an increase in mobility occurred in the period from April 29, 2020, to May 10, 2020. Therefore, we set the start of the perturbation, T∗, as April 29, 2020, and its duration θw=11 days. The second reductions in both effective contact rates are set to q21=0.3 and q22=0.4, while θ2=80 days. Other parameter values are fixed and given in Table 1.
Figure 8 shows simulations for κ=3. We compare our results with the reported confirmed cases and deaths per day. Our projections go until July 01, 2020. Blue and yellow bars represent confirmed and suspect+confirmed cases, respectively. Observe that the mortality is not well described by our model (Figure 8). This observed decrease in mortality is intriguing. It could be due to enhanced treatment of grave cases, shifting of the morbidity towards age classes with a lower risk of death or, perhaps, incomplete mortality records and reporting time delays. Since we are comparing with the absolute number of deaths, testing has little impact in this case.
Currently, an important question relates to the consequences of lifting mitigation measures too early. The concern is that an increase in the number of cases may occur, as it has in the USA. In this section, we explore the possibility of this scenario for Mexico City.
Mexico's federal government developed an epidemiological panel to oversee the reactivation of economic activities in the different states of the country. The panel is based on an index that essentially combines incidence and mortality data, and it has been used since June 01, 2020. The aim was to gradually increase the population mobility to keep the hospital demand and capacity under control. It has four colors: red, orange, yellow, and green.
By the third week of July, the epidemiological panel for Mexico City has had two changes: it started under a red flag on June 01, 2020, the initial date for lifting of mobility restrictions, and it changed to an orange flag on June 29, 2020. A red flag allows only essential economic activities and brief trips to pharmacies, supermarkets, doctor, etc. An orange flag allows added non-essential activities with businesses staffed only to a 30% capacity, opening of public spaces with reduced capacity, etc [21]. On both dates, an unknown number of individuals abandoned confinement and returned to their daily pre-epidemic activities.
We model this process by updating the number of individuals at a fixed time T∗ in Eq (2.2), as x1(ti)=(1−q)x1(ti−1) and x2(ti)=qx1(ti−1)+x2(ti−1) where ti=T∗, x1 and x2 represent the susceptible, exposed, asymptomatically infected and symptomatically infected for the confined and unconfined groups, respectively.
Now we explore the possible consequences of lifting mitigation measures on June 29, 2020. We have chosen this date based on an actual change of the Mexico City epidemiological panel (a change from red to orange). We explore the short and medium-term effects of the lifting of the mitigation measures on the epidemic curve. For each case, we project three scenarios: 25%, 50%, or 75% of the confined population abandoning confinement. We present three curves that have a good fit to Mexico City data. These curves are obtained using the parameters shown in Tables 1 and 2.
Parameter | Blue line | Red line | Black line |
b | 0.89 | 0.91 | 0.91 |
ρ | 0.70 | 0.64 | 0.66 |
δ−1r | 12.07 | 13.84 | 12.26 |
Both effective contact rates (for the confined and unconfined groups) decrease until June 29, 2020, and then both remain constant. For this scenario, reductions in effective contact rates are described by q11=0.6,q21=0.3 (for the confined group), and q12=0.7,q22=0.4 (for the unconfined group).
Figure 9 shows the behavior of reported cases Ir. Solid curves are the baseline, while dotted curves represent the perturbed system (after lifting restrictions on June 29, 2020). Note that, as the percentage of the confined population is decreased, the incidence increases to higher levels. Also, for the scenario of 25% of the population leaving confinement (i.e., increasing its effective contact rate), we obtain a good data fit.
We explore the scenario where the effective contact rate continues to decline even after partially lifting mitigation measures. Here, we consider that the effective contact rates in both groups decrease until August 01, 2020, and then both remain constant. For this scenario, we set q11=0.6,q21=0.2 (confined group), and q12=0.7,q22=0.3 (unconfined group). Other parameter values are similar to those employed in the above subsection.
Figure 10 underlines the importance of a continuing decrease of the transmission rate after lifting confinement. Although the contingency measures are lifted on June 29, even in the worst-case scenario (Figure 10D), there can be a reduction of approximately 50% of newly reported cases (see Figure 9D) after lifting the mitigation measures.
Remark: Figures 9 and 10 show that, as the percentage of the population under confinement decreases, the number of new reported cases correspondingly increases. Similar behavior is observed for the dynamics of daily deaths.
The analysis of interventions, and the important characteristics that these determine both in the epidemics evolution has been discussed by several authors (e.g., [22,23,24,25]). For the case of Mexico City, the precise epidemiological situation of COVID-19 is difficult to accurately assess because, among other things, testing is limited, the positivity rate is high and there is substantial under-reporting of cases and mortality. Our case study uses the confirmed cases reported by the Mexico City government on its official web page. Results show that the growth rate of the epidemic was reduced upon the implementation of mitigation measures on March 23, 2020. However, for Mexico City, the federal government originally forecast the peak of maximum incidence to occur by May 8–10, 2020 and the date for lifting mitigation measures was set to start on June 1, 2020. The peak, however, occurred later in the month, as the data and independent models have shown [6,7,12]. The epidemic in Mexico City entered a period of constant incidence from around May 20, 2020 that, as mentioned in the previous sections, still continues to date. The end result was that the lifting of mobility restrictions coincided with this interval of maximum incidence. Events of high mobility and increased transmission occurring during confinement at times previous to the expected epidemic peak (e.g., children’s day and mother’s day in Mexico City), may have affected the epidemic curves, pushing them into the plateau behavior observed in the data, later reinforced by the lifting of mobility restrictions, as described in the different scenarios presented in this work. Our results are qualitative and strategic in nature but, we believe, they illuminate the importance of counting with reasonably accurate estimates for the occurrence of maximum incidence and the need for constant surveillance and evaluation of atypical high transmission events during confinement. A characteristic of our model that we think worth underlining is that the effective contact rates are continuous, time-varying functions that constitute a realistic approximation, since we are looking at population averages. For example, adopting safe behaviors, such as wearing face masks, involves a learning process; it does not suddenly occur. Rather, it takes time for the face masks to be adopted by a significant proportion of the general population. Our results indicate that, after lifting mobility restrictions, a decrease of the effective contact rate should continue in order to force the epidemic curve to make a downward turn. We interpret this continuing decrease in the effective contact rate as related to the use of face masks, social distancing, washing hands, etc. With respect to extensions on our work we are generalizing our approach using the principles layed off in [26].
Mathematical models are essential in the fight against COVID-19. They are tools for evaluating mitigation measures, estimating mortality and incidence, and projecting scenarios to help public health decision-makers in their very difficult and important task of controlling the epidemic as reviewed in [27]. In this paper, we have used mathematical models to evaluate and generate scenarios. Although precise forecasting is not our aim, we consider that these results can be helpful for decision-makers.
All authors acknowledge support from DGAPA-PAPIIT-UNAM grant IV100220 (convocatoria especial COVID-19). JXVH acknowledges support from grant PAPIIT-DGAP IN115720. We thank Dr. Héctor Benitez and Dr. William Lee Ardavín for their support and Ruth Corona-Moreno for her technical help. We also thank the Secretaría de Ciencia, Tecnología e Innovación of the Government of Mexico City for facilitating the data on the COVID-19 epidemic.
The authors declare no conflicts of interest.
Let Yj, for j=1,2,...,n, be the number of observed new daily cases at time tj, with tj given in days since the first reported case started symptoms. We assume that Yj follows a Negative Binomial distribution with mean value X(tj|a,r,K)=C(tj|a,r,K)−C(tj−1|a,r,K) and dispersion parameter α. C(t|a,r,K) is the solution of the Richards model (3.1). Assuming that, given the parameters, the observations Y1,Y2,…,Yn are conditionally independent, then
E[Yj|a,r,K,α]=X(tj|a,r,K) | (5.1) |
Var[Yj|a,r,K,α]=X(tj|a,r,K)+αX(tj|a,r,K)2. | (5.2) |
The Negative Binomial distribution allows us to control the variability of the data by considering over-dispersion, which is common for epidemiological data. If α=0, then we return to the Poisson model, which is often used in this context.
Let θ=(a,r,K,α) be the vector of parameters that we wish to estimate. The inclusion of the parameter α is related to the variability of the data, not to the Richards model, is necessary since, in practice, the variability is unknown. Then, the likelihood function, which represents how plausible the data is under the Negative Binomial assumption and the Richards model if we knew the parameters, is given by
π(y1,…,yn|θ)=n∏j=1Γ(yj+τ)Γ(yj)Γ(τ)(ττ+C(tj|a,r,K))τ(C(tj|a,r,K)τ+C(tj|a,r,K))yj. | (5.3) |
Consider that parameters a, r, K and α as random variables. Assuming prior independence, the joint prior distribution for vector θ is
π(θ)=π(a)π(r)π(K)π(α), |
where π(a) is the probability density function (pdf) of a Uniform(0, 2) distribution, π(r) is the pdf of a Uniform(0, 2), π(K) is the pdf of a Uniform(Kmin, Kmax), and π(α) is the pdf of a Gamma(shape = 2, scale = 0.1) distribution. To select the prior for parameter r, we note that previous estimations of r are close to 0.3 [14], and a Uniform(0, 2) represents a weekly informative prior, as it allows for a wide range of values of r. Also, there is no available prior information regarding the final size of the outbreak K. This is a critical parameter in the model and, in order to avoid bias, we assume a uniform prior over Kmin and Kmax. To set these last two values, we assume that the minimum number of confirmed cases is the maximum number of observed cumulative cases Y(tn) (for the studied period) times two; i.e., Kmin=yn×2. To set the upper bound for K, we consider a fraction of the total population Kmax=N×0.05, where N is the population size of Mexico. This fraction was determined based on the observations of other cities such as the New York where the total population size is similar to Mexico City.
Then the posterior distribution of the parameters of interest is
π(θ|y1,…,yn)∝π(y1,…,yn|θ)π(θ), |
and it does not have a simple form, since the likelihood function depends on the solution of the Richards model, which is not linear in the parameters. We analyze the posterior distribution using an MCMC algorithm called t-walk that does not require tuning [28]. This algorithm generates samples from the posterior distribution that can be used to estimate marginal posterior densities, mean, variance, quantiles, etc. We refer the reader to [15] for more details on MCMC methods and to [16] for an introduction to Bayesian inference with differential equations.
We estimate the vector θ using data from Mexico City for two different periods, from February 22, 2020 to March 22, 2020, and from March 23, 2020 to April 30, 2020. Table 3 shows the median posterior estimates and 95% probability intervals for the parameters in each period.
February 22 – March 22 | March 23 – April 30 | ||||||
Lower | Median | Upper | Lower | Median | Upper | ||
a | 0.02 | 0.05 | 0.38 | 0.12 | 0.27 | 0.45 | |
K | 1,186 | 26,156 | 167,504 | 52,201 | 87,608 | 174,583 | |
r | 0.28 | 0.76 | 1.91 | 0.11 | 0.12 | 0.21 |
[1] |
J. Dehning, J. Zierenberg, F. P. Spitzner, M. Wibral, J. P. Neto, M. Wilczekand, et al., Inferring change points in the spread of COVID-19 reveals the effectiveness of interventions, Science, 369 (2020), eabb9789. doi: 10.1126/science.abb9789
![]() |
[2] |
T. Sardar, S. S. Nadim, S. Rana, J. Chattopadhyay, Assessment of Lockdown Effect in Some States and Overall India: A Predictive Mathematical Study on COVID-19 Outbreak, Chaos Solitons Fractals, 139 (2020), 110078. doi: 10.1016/j.chaos.2020.110078
![]() |
[3] |
A. L. Bertozzi, E. Franco, G. Mohler, M. B. Short, D. Sledge, The challenges of modeling and forecasting the spread of COVID-19, Proc. Natl. Acad. U S A, 117 (2020), 16732-16738. doi: 10.1073/pnas.1914072117
![]() |
[4] | H. S. Badr, H. Du, M. Marshall, E. Dong, M. M. Squire, L. M. Gardner, Association between mobility patterns and COVID-19 transmission in the USA: a mathematical modelling study, Lancet Infect. Dis., DOI: 10.1016/S1473-3099(20)30553-3. |
[5] |
C. N. Ngonghala, E. Iboi, S. Eikenberry, M. Scotch, C. R. MacIntyre, M. H. Bonds, et al., Mathematical assessment of the impact of non-pharmaceutical interventions on curtailing the 2019 novel Coronavirus, Math. Biosci., 325 (2020), 108364. doi: 10.1016/j.mbs.2020.108364
![]() |
[6] |
M. A. Acuña-Zegarra, M. Santana-Cibrian, J. X. Velasco-Hernandez, Modeling behavioral change and COVID-19 containment in Mexico: A trade-off between lockdown and compliance, Math. Biosci., 325 (2020), 108370. doi: 10.1016/j.mbs.2020.108370
![]() |
[7] | M. Santana-Cibrian, M. A. Acuna-Zegarra, J. X. Velasco-Hernandez, Flattening the curve and the effect of atypical events on mitigation measures in Mexico: a modeling perspective, preprint, May 2020, medRxiv: 10.1101/2020.05.21.20109678. |
[8] | S. Lee, T. Kim, E. Lee, C. Lee, H. Kim, H. Rhee, et al., Clinical course and molecular viral shedding among asymptomatic and symptomatic patients with SARS-CoV-2 infection in a community treatment center in the Republic of Korea, JAMA Intern. Med., DOI: 10.1001/jamainternmed.2020.3862. |
[9] |
L. Zou, F. Ruan, M. Huang, L. Liang, H. Huang, Z. Hong, et al., SARS-CoV-2 viral load in upper respiratory specimens of infected patients, N. Engl. J. Med., 382 (2020), 1177-1179. doi: 10.1056/NEJMc2001737
![]() |
[10] | N. Lozano, Ssa estima pico máximo de casos por covid-19 entre 8 y 10 de mayo, 2020. Available from: https://www.politico.mx/minuta-politica/minuta-politica-gobierno-federal/ssa-estima-pico-m%C3%A1ximo-de-casos-por-covid-19-entre-8-y-10-de-mayo/, Político.mx. |
[11] | E. Castillo, En valle de méxico se prevé fin de epidemia por coronavirus el 25 de junio: López-gatell, 2020. Available from: https://www.milenio.com/politica/coronavirus-valle-mexico-preven-epidemia-25-junio, Milenio. |
[12] | R. H. Mena, J. X. Velasco-Hernandez, N. B. Mantilla-Beniers, G. A. Carranco-Sapiéns, L. Benet, D. Boyer, et al., Using posterior predictive distributions to analyse epidemic models: COVID-19 in Mexico City, Phys. Biol., (2020). DOI: 10.1088/1478-3975/abb115. |
[13] | Gobierno de la Ciudad de México, Datos abiertos, 2020. Available from: https://datos.cdmx.gob.mx/explore/dataset/base-covid-sinave/table/. |
[14] |
K. Roosa, Y. Lee, R. Luo, A. Kirpich, R. Rothenberg, J. Hyman, et al., Real-time forecasts of the covid-19 epidemic in China from February 5th to February 24th, 2020, Infect. Disease Model., 5 (2020), 256-263. doi: 10.1016/j.idm.2020.02.002
![]() |
[15] | C. P. Robert, G. Casella, Monte Carlo Statistical Methods, 2nd edition, Springer Texts in Statistics, Springer New York, 2004. |
[16] | J. P. Kaipio, E. Somersalo, Statistical and Computational Inverse Problems, vol. 160 of Applied Mathematical Sciences, Springer-Verlag, New York, 2005. |
[17] | Gobierno de la Ciudad de México, Información laboral, 2020. Available from: http://www.stps.gob.mx/gobmx/estadisticas/pdf/perfiles/perfil%20distrito%20federal.pdf, |
[18] |
R. N. Thompson, J. E. Stockwin, R. D. van Gaalen, J. A. Polonsky, Z. N. Kamvar, P. A. Demarsh, et al., Improved inference of time-varying reproduction numbers during infectious disease outbreaks, Epidemics, 29 (2019), 100356. doi: 10.1016/j.epidem.2019.100356
![]() |
[19] |
A. Cori, N. M. Ferguson, C. Fraser, S. Cauchemez, A New Framework and Software to Estimate Time-Varying Reproduction Numbers During Epidemics, Am. J. Epidemiol., 178 (2013), 1505-1512. doi: 10.1093/aje/kwt133
![]() |
[20] |
H. Nishiura, N. M. Linton, A. R. Akhmetzhanov, Serial interval of novel coronavirus (COVID-19) infections, Int. J. Infect. Dis., 93 (2020), 284-286. doi: 10.1016/j.ijid.2020.02.060
![]() |
[21] | Gobierno de la Ciudad de México, Sem áforo epidemiológico COVID-19, 2020. Available from: https://semaforo.covid19.cdmx.gob.mx/tablero/, 2020. |
[22] |
X. Rong, L. Yang, H. Chu, M. Fan, Effect of delay in diagnosis on transmission of COVID-19, Math. Biosci. Eng., 17 (2020), 2725-2740. doi: 10.3934/mbe.2020149
![]() |
[23] |
C. Yang, J. Wang, A mathematical model for the novel coronavirus epidemic in Wuhan, China, Math. Biosci. Eng., 17 (2020), 2708-2724. doi: 10.3934/mbe.2020148
![]() |
[24] |
L. Wang, J. Wang, H. Zhao, Y. Shi, K. Wang, P. Wu, et al., Modelling and assessing the effects of medical resources on transmission of novel coronavirus (COVID-19) in Wuhan, China, Math. Biosci. Eng., 17 (2020), 2936-2949. doi: 10.3934/mbe.2020165
![]() |
[25] |
G. Giordano, F. Blanchini, R. Bruno, P. Colaneri, A. Di Filippo, A. Di Matteo, et al., Modelling the covid-19 epidemic and implementation of population-wide interventions in Italy, Nat. Med., 26 (2020), 855-860. doi: 10.1038/s41591-020-0883-7
![]() |
[26] |
H. Zhao, Z. Feng, Staggered release policies for covid-19 control: Costs and benefits of relaxing restrictions by age and risk, Math. Biosci., 326 (2020), 108405. doi: 10.1016/j.mbs.2020.108405
![]() |
[27] |
E. Estrada, COVID-19 and SARS-CoV-2. Modeling the present, looking at the future, Phys. Rep., 869 (2020), 1-51. doi: 10.1016/j.physrep.2020.07.005
![]() |
[28] |
J. A. Christen, C. Fox, A general purpose sampling algorithm for continuous distributions (the t-walk), Bayesian Anal., 5 (2010), 263-281. doi: 10.1214/10-BA603
![]() |
1. | Manuel Adrian Acuña-Zegarra, Saúl Díaz-Infante, David Baca-Carrasco, Daniel Olmos Liceaga, COVID-19 optimal vaccination policies: A modeling study on efficacy, natural and vaccine-induced immunity responses, 2021, 00255564, 108614, 10.1016/j.mbs.2021.108614 | |
2. | Carlos G. Aguilar-Madera, Gilberto Espinosa-Paredes, E.C. Herrera-Hernández, Jorge A. Briones Carrillo, J. Valente Flores-Cano, Víctor Matías-Pérez, The spreading of Covid-19 in Mexico: A diffusional approach, 2021, 27, 22113797, 104555, 10.1016/j.rinp.2021.104555 | |
3. | Miaolei Li, Jian Zu, Zongfang Li, Mingwang Shen, Yan Li, Fanpu Ji, How to Reduce the Transmission Risk of COVID-19 More Effectively in New York City: An Age-Structured Model Study, 2021, 8, 2296-858X, 10.3389/fmed.2021.641205 | |
4. | Fernando Saldaña, Jorge X. Velasco-Hernández, Modeling the COVID-19 pandemic: a primer and overview of mathematical epidemiology, 2022, 79, 2254-3902, 225, 10.1007/s40324-021-00260-3 | |
5. | Mario Santana-Cibrian, Manuel A. Acuña-Zegarra, Marco Tulio Angulo, Andreu Comas-García, Esteban A. Hernández-Vargas, Jorge X. Velasco-Hernandez, 2022, Chapter 16, 978-3-030-72833-5, 543, 10.1007/978-3-030-72834-2_16 | |
6. | Jinxing Guan, Yang Zhao, Yongyue Wei, Sipeng Shen, Dongfang You, Ruyang Zhang, Theis Lange, Feng Chen, Transmission dynamics model and the coronavirus disease 2019 epidemic: applications and challenges, 2022, 2, 2749-9642, 89, 10.1515/mr-2021-0022 | |
7. | Philip J. Gerrish, Fernando Saldaña, Alexandre Colato, Erika E. Rodriguez Torres, T.Y. Okosun, Jorge X. Velasco-Hernández, 2023, 9780323950640, 93, 10.1016/B978-0-323-95064-0.00007-5 | |
8. | Mario Santana-Cibrian, Manuel A. Acuña-Zegarra, Mayra R. Tocto-Erazo, Ruth Corona-Moreno, 2023, 9780323950640, 43, 10.1016/B978-0-323-95064-0.00005-1 | |
9. | M. Acuña-Zegarra, Mario Santana-Cibrian, Carlos Hernandez-Vela, Ramsés Mena, Jorge Velasco-Hernández, 2024, 793, 978-1-4704-7606-9, 229, 10.1090/conm/793/15878 | |
10. | Hugo Flores-Arguedas, José Ariel Camacho-Gutiérrez, Fernando Saldaña, 2023, Chapter 43-1, 978-3-030-93954-0, 1, 10.1007/978-3-030-93954-0_43-1 | |
11. | Wenyu Du, Zhenghui Feng, Yi Zhao, The impact of lifting COVID-19 restrictions on influenza transmission across countries, 2025, 2025, 2731-4235, 10.1186/s13662-025-03872-5 |
Parameter | Value | Parameter | Value | Parameter | Value |
b | 0.887677 | γ−1 | 5.1 | T | 30 |
ρ | 0.698725 | η−1a | 6 | q11 | 0.6 |
μ | 0.11 | η−1s | 7 | θ1 | 18 |
ω | 0.005 | δ−1s | 4 | q12 | 0.7 |
q | 0.7 | δ−1r | 12.070691 | T∗ | 67 |
Parameter | Blue line | Red line | Black line |
b | 0.89 | 0.91 | 0.91 |
ρ | 0.70 | 0.64 | 0.66 |
δ−1r | 12.07 | 13.84 | 12.26 |
February 22 – March 22 | March 23 – April 30 | ||||||
Lower | Median | Upper | Lower | Median | Upper | ||
a | 0.02 | 0.05 | 0.38 | 0.12 | 0.27 | 0.45 | |
K | 1,186 | 26,156 | 167,504 | 52,201 | 87,608 | 174,583 | |
r | 0.28 | 0.76 | 1.91 | 0.11 | 0.12 | 0.21 |
Parameter | Value | Parameter | Value | Parameter | Value |
b | 0.887677 | γ−1 | 5.1 | T | 30 |
ρ | 0.698725 | η−1a | 6 | q11 | 0.6 |
μ | 0.11 | η−1s | 7 | θ1 | 18 |
ω | 0.005 | δ−1s | 4 | q12 | 0.7 |
q | 0.7 | δ−1r | 12.070691 | T∗ | 67 |
Parameter | Blue line | Red line | Black line |
b | 0.89 | 0.91 | 0.91 |
ρ | 0.70 | 0.64 | 0.66 |
δ−1r | 12.07 | 13.84 | 12.26 |
February 22 – March 22 | March 23 – April 30 | ||||||
Lower | Median | Upper | Lower | Median | Upper | ||
a | 0.02 | 0.05 | 0.38 | 0.12 | 0.27 | 0.45 | |
K | 1,186 | 26,156 | 167,504 | 52,201 | 87,608 | 174,583 | |
r | 0.28 | 0.76 | 1.91 | 0.11 | 0.12 | 0.21 |