key: cord-0855962-08wm1498 authors: Mohammadi, Z.; Cojocaru, M. G.; Thommes, E. W. title: Human behaviour, NPI and mobility reduction effects on COVID-19 transmission in different countries of the world date: 2022-05-01 journal: nan DOI: 10.1101/2022.04.29.22274485 sha: ab63397f0122559273f1699d8f3fe4d591676cfa doc_id: 855962 cord_uid: 08wm1498 Background: The outbreak of Coronavirus disease, which originated in Wuhan, China in 2019, has affected the lives of billions of people globally. Throughout 2020, the reproduction number of COVID-19 was widely used by decision-makers to explain their strategies to control the pandemic. Methods: In this work, we deduce and analyze both initial and effective reproduction numbers for 12 diverse world regions between February and December of 2020. We consider mobility reductions, mask wearing and compliance with masks, mask efficacy values alongside other non-pharmaceutical interventions (NPIs) in each region to get further insights in how each of the above factored into each region's SARS-COV-2 transmission dynamic. Results: We quantify in each region the following reductions in the observed effective reproduction numbers of the pandemic: i) reduction due to decrease in mobility (as captured in Google mobility reports); ii) reduction due to mask wearing and mask compliance; iii) reduction due to other NPI's, over and above the ones identified in i) and ii). Conclusion: In most cases mobility reduction coming from nationwide lockdown measures has helped stave off the initial wave in countries who took these types of measures. Beyond the first waves, mask mandates and compliance, together with social-distancing measures (which we refer to as other NPI's) have allowed some control of subsequent disease spread. The methodology we propose here is novel and can be applied to other respiratory diseases such as influenza or RSV. The first known case of disease caused by SARS-CoV-2 was identified in Wuhan, China, in December 2019. The disease spread worldwide in a few weeks, leading to a pandemic still ongoing as of the winter of 2022. Since December 2019, the basic and effective reproduction number of COVID-19 have been continuously discussed by scientists, political decision-makers and the media (regular and social). The basic reproduction number of an infectious disease, denoted by R 0 , represents the expected number of new cases generated by a single infectious individual in a fully susceptible population. In epidemiology, a pandemic will be under control and the transmission will die out when R 0 < 1. Estimating the basic reproduction number for COVID-19 was challenging since it was depending on the behavioural activity in local populations. The basic reproduction number for COVID-19 was reported to take values from 2.2 [1] to 6.33 [2] in different countries. In general, although whole populations started with being susceptible to virus, the progression of the pandemic steadily decreased the number of susceptibles in each region. Consequently, the average number of secondary cases per infectious case changed as the population became immunized (by recovering or dying). Thus as the pandemic progressed, R 0 gives way to the effective reproductive number, denoted by R ef f , which measures the average number of secondary cases per infectious case in a population at any specific time (where only a fraction is susceptible) [3] . In epidemiology, R ef f is a monitoring indicator of progress in controlling a pandemic. It can also be a way to monitor the effectiveness of interventions (both the effect of immunity and of additional non-pharmaceutical interventions) during a pandemic. Transmission of SARS-COV-2 depends on the rate of person-to-person contact and on the probability of transmission given one meaningful contact between an infected and a susceptible individual. While the probability of transmission per contact is reduced by NPI's, principally the wearing of masks [4, 5] , the effectiveness of mask-wearing in preventing transmission of SARS-CoV-2 was not recognized in some world regions during the first several months of the pandemic, despite the effectiveness of public mask use in controlling the spread of the 2003 SARS [6, 7] . Another critical non-pharmaceutical intervention for slowing epidemic growth is social distancing which includes shelter-in-place requirements, prohibition of indoor gatherings, imposition of travel restrictions, school closures and workplace closures. Previous studies investigated the quantitative relationships between the COVID-19 epidemic parameters (for instance the total death toll) and social distancing efforts [8, 9, 10, 11] . Mobile data provide a unique opportunity to investigate to some degree the effectiveness of social distancing measures in reducing the effective reproduction number of COVID-19 [12, 13, 14] . Mobile data can be interpreted as a proxy of person-to-person contact reduction as a consequence of social distancing measures, although it generally does not capture short-range behavioral changes, such as maintaining a 6 foot spacing between individuals in public settings. The publicly available data on human mobility that is provided by Google, Apple, Facebook, etc have been used in several articles to evaluate the effectiveness of non-pharmaceutical interventions (NPIs) on the spread of COVID-19 [15, 16, 17, 18, 19, 20] . Existing pre-pandemic work, such as [21, 22, 23] , represents the effect of the population-level contact patterns on the infection transmission dynamics. The infection transmission rate varies in different countries as a result of the different social and economic structures and different contact patterns. Mistry et al [22] used the derived contact matrices to model the spread of airborne infectious diseases. The transmission rate of COVID-19 fundamentally depends on interpersonal interaction rates. Most countries considered different strategies for reducing the contact rate in order to control the spread of the virus. Feehan and Mahmud calculated agestructured contact matrices to quantify how much interpersonal contact has changed during the pandemic in the United States [24] . They estimated about 82% decline in interpersonal contact between March 22nd and April 8-th, 2020 (known as wave 0) and an increase in daily average contact rates over the subsequent waves. Other studies observed the decline in contact rates in China [25] , United Kingdom [26] , Luxembourg [27] , Italy, Belgium, France, and the Netherlands [28] throughout the pandemic. Prem et al. created synthetic contact matrices to represent the effect of intervention measures to reduce social mixing on outcomes of the COVID-19 epidemic in Wuhan, China [29, 30] . The age-specific social contact characterization also supports the possibility of suspecting differences in transmission patterns of COVID-19 outbreak among different age-groups [31, 32, 33, 20] . In this work, projected contact matrices provided by Prem et al. [21] are used. This work presents a SEIRL (Susceptible, Exposed, Infectious, Recovered and isoLated) model that uses incidence data, Google mobility data (as a modifier of effective contact rates), mask efficacy and mask compliance data to provide insight into the strategies employed for controlling the pandemic in each country under consideration. It explores and quantifies the effectiveness of non-pharmaceutical policy interventions across 12 diverse countries/regions, using their effective reproduction number R ef f . To accomplish these goals we estimate the R 0 values of each region (using the classic next-generation matrix analysis [34]) and then a time-dependent effective reproduction number, denoted by R ef f −data , using known repositories of incidence data for each country [35] . Using the projected contact rates from [21] extrapolated to each country of interest with publicly available demographic data, we also infer the time-dependent effective contact rates, and hence the time-dependent effective reproduction number of each country as a function of multiple Google mobility indices. We denote this by R ef f −mobil . Comparing the two estimates of R ef f we investigate the effect of NPI's, over and above changes in contact rates due to mobility and mask mandates, in each country's epidemiology throughout the year 2020 [1] . It is concluded in most cases that mobility reduction coming from stringent lockdown measures helped stave off the initial wave in countries which took these types of measures (see for instance Figure 2 ). While mobility increased in the second part of 2020 (see Figure 7) , mask mandates together with all other NPI measures (for instance social-distancing) have allowed some "control" of subsequent waves, in the sense that countries were able to maintain their effective reproduction numbers around 1. The structure of the paper is as follows: Section 2 presents the setup and methodology employed in the analyses of each of the twelve regions. Section 3 presents the results and discussion on quantification of the reduction in the effective reproduction numbers R ef f −mobile and R ef f −mobilemask , as well as the estimated effective reproduction numbers based on incidence data R ef f −data . The paper closes with a discussion of the results and a few directions for future work. We chose 10 countries around the world, one U.S. state (Florida) and one province of Canada (Ontario), based on a diversity of characteristics: population density, median age, urbanization of population and gross domestic product (GDP) in 2020 (theses various characteristics are presented in Table 1 Table 1 : Countries under consideration. [1] We chose to concentrate on the year 2020 specifically because there were no preventative or antiviral treatments known against this virus at that time, and countries had to rely on some combination of NPI's to fight its spread. All rights reserved. No reuse allowed without permission. (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 1, 2022. ; https://doi.org/10.1101/2022.04. 29.22274485 doi: medRxiv preprint The data are collected from IndexMundi . In each of these regions we are interested to model the pandemic evolution during the year 2020 via a SEIRL model described below. Susceptible individuals (denoted by S and considered to be the entire population initially) exposed to the virus enter the exposed (E) compartment for an average of 1/σ days before they become contagious, at which point they move into the I (infected) compartment. In general, there is an a proportion of infected individuals who will not develop symptoms (so-called asymptomatic, denoted by I A ), while the remaining 1 − a percentage make up the infected symptomatic individuals, denoted by I S . Thus here : An proportion of I S (t) will self-isolate into the L (isolated) compartment. They do so with a delay of 1/κ days, accounting for a test result wait time and/or individuals who may disregard minor symptoms initially. After 1/γ days, individuals recover from (or succumb to) their infection and move into the R (recovered) compartment. A diagram of the process we model is as follows 1: Using the expressions in (1) and model diagram (1), we obtain the following equations: (2) Our model parameters have been taken from literature (as can be seen in Table 2 ), with the exception of the rate of isolation , which we assume to be equal to 95%. All rights reserved. No reuse allowed without permission. (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. Time series of effective reproduction numbers using a near-disease-free equilibrium estimate and incidence data The Jacobian analysis near the disease-free equilibrium (DFE) which consists of S(0) = N and I(0) = 0) for the system (2) and the next generation matrix method ([34]) for computing the initial R 0 are given in Appendix. We have employed a similar type of approach on our papers [41, 42] , where the compartmental models were slightly different. We obtain its closed form expression: We can also estimate R 0 as a function of the growth factor near the DFE in each region using (9) (see more details in Appendix): To estimate the exponential growth factor for each region as a time series, we rely on weekly case incidence data (which we denote by inc(t)) for each region. We assume that inc(t) is given by an exponential curve of the type: In this case, we can compute a time-series of the exponential growth factor: , with inc(t) = 0. Using this time series of ρ(t) in equation (10) we obtain a time series for the effective reproduction number: , where s(t) represents the remaining fraction of susceptibles at each t, (5) i.e. the difference between the entire population and the current cumulative number of infected individuals: inc(τ ). Figure 2 represents the weekly changes in the effective reproduction numbers of incidence data, R ef f −data that is explained in section (2.1), throughout 2020 for each geographical location under study. The vertical lines in each panel represent the dates in which local governments have introduced nationwide measures: in most countries lockdown measures took effect, while in Sweden and Indonesia partial lockdown measures were in place. In Sweden, nationwide lockdown was considered to be a violation of people's freedom of movement. All rights reserved. No reuse allowed without permission. (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 1, 2022. ; https://doi.org/10.1101/2022.04. 29.22274485 doi: medRxiv preprint The government strategy was based on individual responsibility. Measures included were: border closures, recommendations about social distancing, and traveling. Local reports show a 50% drop in public transport usage in April for the Swedish counties [43], [44] . In Stockholm, a 30% drop in the number of cars [45] , and 70% fewer pedestrians [46] was reported in April 2020. Later in November, the Swedish government imposed mandatory restrictions as well (e.g. all gatherings of more than eight people were banned).Italy had one of the earliest implementations of lockdown, as early as February 23, 2020. We showcase all dates for all locations in Table 3 . Hard Time-series of effective reproduction numbers accounting for mobility data and projected contact rates We first devise a mechanism to mix the daily average contact rates in each region's population with the behavioural activity in that population. In Prem et. al. [21] the authors compute projected daily average contact rates for 157 countries. Specifically, they estimate the average contacts for categories of activities during a typical day, such as: home, work and other locations. They present their results for a population stratified by age, and divided into 16 age subgroups (See Appendix for the details). We amalgamated the average contact rate in home, work, and other locations, then we computed the weighted average of a given projected contact matrix in a region with the corresponding proportions of 5-year age population groups in 2020 to determine one single average contact rate. We consider this last value as the average baseline (pre-lockdown) contact rate in that region (for instance, for Ontario it was computed to be: contact av ≈ 11.04). All regions' contacts are reported below in Table 4 . Google reports contain changes in movement over time, compared to baseline (pre-lockdown) activity in six categories: retail/recreation, groceries/pharmacies, parks, transit stations, workplaces, and domiciles (COVID-19 Community Mobility Reports) [47]. We have used the Google index data in other works, see [20, 53] . To find the mobility-influenced, time-dependent contact rates post-lockdown, we considered the average contacts rate All rights reserved. No reuse allowed without permission. (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 1, 2022. ; https://doi.org/10.1101/2022.04. 29.22274485 doi: medRxiv preprint (10), while the vertical pink dash line represents the start of nationwide social distancing orders in each region under consideration. All rights reserved. No reuse allowed without permission. (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 1, 2022. ; https://doi.org/10.1101/2022.04. 29.22274485 doi: medRxiv preprint for the home, work, and other location categories (comprising retail/recreation and groceries/pharmacies) from Prem et al. [21] for each region. Next, we used these category rates to modify the same categories of mobility data as follows: Let us now refine our view of the effective reproduction number R ef f = R 0 s(t) from the perspective of the changes in contact rates due to mobility reduction. From the last section we know that , however now we have where p is the probability of transmission per meaningful contact. To estimate p, we need to have a handle on the values of R 0 from incidence data as well. To estimate R 0 from incidence data we used equation (10). The growth factor ρ is computed from the initial phase of (close to) exponential growth in the neighbourhood of the DFE, which corresponds to a phase of linear growth in log(inc), with slope ρ. We identify the initial, fastest phase of nearly unchecked growth in any given region with the help of a piecewise linear fit to the log of the incidence. We utilize the R function dpseg(), which is a part of the dpseg package, https://cran.r-project. org/web/packages/dpseg/index.html. This function uses a dynamic programming algorithm to generate an optimal piecewise linear fit to a time series, which balances goodness of fit against an (adjustable) penalty for each additional segment. We then identify the earliest segment with the steepest positive slope (largest ρ) as corresponding to the initial near-unchecked exponential growth phase. We report the values we obtain in Table 2 . [2] [2] The last column in Table 2 highlights the number of weeks that dpseg is assigning for the steepest slope. The weeks are numbered from the first week with positive cases in each location, so near disease-free. All rights reserved. No reuse allowed without permission. (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. 3.1 Quantifying the effective reproduction number using incidence data To quantify the overall relative status of the pandemic at the various locations, we show in Figure 3 (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 1, 2022. ; https://doi.org/10.1101/2022.04. 29.22274485 doi: medRxiv preprint The color gradient signifies lowest values of the effective reproduction numbers in lightest color zones, and highest values in darkest zones. Figure 3 showcases a comparison of the R ef f −data between locations. We can immediately see the effect of the reduction in R ef f −data , across the board, in all countries between March and April, and the rest of the year. It is clear the drop achieved by initial lockdowns and mobility reduction allowed all countries to drop their R ef f −data to around 1, and interestingly most remained around this value for the rest of 2020. We include here ( Table 5 : Average R ef f −data at each location after the initial drop Qualitatively, we can understand why R ef f −data has generally fluctuated in the vicinity of 1: On the one hand, rapid growth of incidence causes alarm at both the policymaker and individual level, and typically leads to increased NPI directives as well as compliance. On the other hand, because these measures are hard to maintain and economically taxing, NPIs are generally relaxed not long after incidence starts to drop. In other words, it is the human behavioural element that we do not model here, but that we glimpse. We observe that the initial sharp reduction in contact rate due to mobility happened in the middle of March in most regions under study, when initial lockdown was in effect, except in Sweden. The contact rate gradually increased from early May when partial reopening was implemented by governments. The results of these measures are noticeable in the figure 4 and 7 upper right panels. The effect of mobility restrictions throughout 2020 for each geographical location under study is presented in Figure 5 below. We plot the theoretically estimate R ef f −mobile numbers together with the R ef f −data effective All rights reserved. No reuse allowed without permission. (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 1, 2022. ; https://doi.org/10.1101/2022.04. 29.22274485 doi: medRxiv preprint numbers at each location. Overall, under Google mobility reductions using formula 6, the mobility reduction effect is not enough, on its own, to explain the values R ef f −data . This is not surprising, as each location had a varying combination of NPI measures. Figures 5 (and 6 in an ensemble view) reveal the reduction in mobility in each region. It seems that the mobility decreased 51% and 54% in Nepal and Italy with respect to the baseline in the first and second weeks of April 2020, respectively. The mobility fell about 20% percent in Indonesia and Sweden. Large-scale social All rights reserved. No reuse allowed without permission. (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 1, 2022. ; https://doi.org/10.1101/2022.04. 29.22274485 doi: medRxiv preprint restrictions (sometimes called partial lockdown) were introduced by the Indonesian government in place of nationwide lockdown at the end of March 2020 while the Swedish authorities imposed some restrictions on gatherings in late November. Figure 6 : Reduction in mobility influenced by contact rate throughout 2020. The week that reduction happened is highlighted by "month/day" in each bar. The effects of mask mandates and mask compliance on further reduction of the effective reproduction numbers We will be looking at two main tools that populations can use to control the pandemic: social distancing and mask wearing. In our framework here, each one of these directly influences the transmission rate β. Denoting by mask ef f the efficacy of an average mask at preventing transmission and by compliance m the compliance with mask wearing let us imagine a meaningful contact between an infected and a susceptible individual: If an individual wears a mask and mask ef f = 0.3, then he/she has an increased protection against transmission. If the individual complies with mask wearing 50% of the time, then their protection due to mask wearing is, on average mask ef f · compliance m = 0.15 which implies that the per-contact transmission probability p will decrease with mask wearing: Recalling that the level of mask-wearing and social distancing both change over time, we estimate that Let us consider only mask-wearing for the time being. In this paper we use a value of mask ef f = 50%, as averaged based on estimates in [54] (see a more detailed discussion in Appendix. and the values of mask compliance from IHME (details of data sources in Data sources section above. ) . Using (7) in the estimate (6) All rights reserved. No reuse allowed without permission. (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 1, 2022. ; https://doi.org/ 10.1101/2022.04.29.22274485 doi: medRxiv preprint leads to: with data from Table 4 . Figure 7 shows the effective reproduction numbers as obtained from the incidence data (yellow curves) from Figure 2 , together with our theoretical estimates from Figure 5 (solid pink curves) using mobility data, to which now we add a new theoretically estimated R ef f −mobilemask number reflecting mask wearing data and formula 8 above. The red and black dashed lines in each panel show the lockdown and mask-wearing mandate dates, respectively, in each region as publicly available. We observe that the theoretically-estimated curve using both mobility reductions and mask wearing data (which we will now denote by R ef f −mobilemask ) is accounting for more of the reduction in transmission than the estimated R ef f −mobile curves of Figure 5 . Evidently, in some regions (Ontario and Argentina) we observe what it looks like an over-reduction in the values of R ef f −mobilemask as compared to R ef f −data , while in other countries, most notably Sweden, we notice essentially no contribution in further reduction from R ef f −mobile to R ef f −mobilemask . In Sweden's case, this is not surprising, as the mask wearing levels in the IHME data we used hover between 1% − 2% throughout 2020. Last but not least, this assumes a values of mask efficacy mask ef f = 50%. If this value is decreased, the R ef f −mobilemask curves will shift upwards, i.e. the reduction of R ef fmobile due to mask wearing will be smaller. As we highlighted in previous section, a host of other NPI's have been employed, at different strengths, across the regions. While it is clear from our investigations thus far that the mobility reduction (as reflected in Google mobility data), along with mask wearing and mask compliance helped tremendously in the de-escalation of the R ef f curves, we wish to further analyze the data to extract more information on the strength of other NPI's. Let us denote by compliance oN P I the compliance level in the local population with other NPI measures (by other we mean other than mobility and mask wearing). With this in mind, an average individual in a local population is protected over the course of their contacts proportional to what fraction of the time they also comply with other NPI measures: contact av = (1 − compliance oN P I (t)) · contact av which then means that This means that to quantify the effects of other NPI measures per region, we further look at the ratio: . We present these plots in Figure 8 below. Clearly when the ratio is less than 1, then the reduction in the R ef f that we quantified from observed data is stronger than the reduction of R ef f based on mobility and mask wearing and viceversa. Specifically we obtain the following: In the other case, under our assumptions, we get All rights reserved. No reuse allowed without permission. (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 1, 2022. ; https://doi.org/ 10.1101/2022.04.29.22274485 doi: medRxiv preprint Figure 7 : A comparison of the effective reproduction number as obtained from the incidence data (yellow curves), to our theoretical estimate from Section 3.1 (solid pink curves) using mobility data. Beyond the date of mask mandate enactment in each region, we show the theoretically-estimated incidence both with (solid pink) and without (dashed pink) the added effect of mask-wearing. The red and black dashed lines in each panel show the lockdown and mask-wearing mandate dates, respectively, in each region as publicly available. In cases where compliance oN P I < 0 we interpret it to mean that while the local population has been very compliant with mask wearing (assumed to be at 50% efficacy in all regions), they may not have been as observant towards other measures. All rights reserved. No reuse allowed without permission. (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 1, 2022. ; https://doi.org/10.1101/2022.04. 29.22274485 doi: medRxiv preprint One remark is worth to be made at this point: if one decreases the mask efficacy mask ef f to an average of 40%, respectively even lower to 30% (as in [54] ), then the R ef f −mobilemask curves of Figure 5 would scale equivalently upward by a constant, thus implying that the mask wearing effect is less effective and therefore that the ratio R ef f −data R ef f −mobil would be less than 1 (i.e., compliance oN P I (t) > 0) for most (respectively all) regions. Results in Sweden in particular stand out (see Figure 9 ). They indicate that other NPI measures were extremely effective in reducing the transmission rate of disease, in the absence of mandated lockdown periods and nearly 0 mask wearing. In the case of Sweden for instance (upper panel of Figure 9 ), compliance oN P I (t) All rights reserved. No reuse allowed without permission. (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 1, 2022. ; https://doi.org/10.1101/2022.04. 29.22274485 doi: medRxiv preprint is large starting in April 2020. Here compliance oN P I (t) > 0 at and over 50% in Sweden. In case of Indonesia, we estimate that compliance oN P I (t) at and over 30-40% in Indonesia. Both these ranges assume an average of mask ef f = 50% over the time period May-December 2020. Figure 9 : Effect of other NPI measures in Sweden and Indonesia. Mandatory masks were never introduced in Sweden until November 2020, but they are present earlier in Indonesia. Nationwide lockdown was never used in these countries in 2020. At the other end of the spectrum, Ontario (upper left corner panel in Figure 8 ) seems to display a ratio of R ef f −data R ef f −mobil ≈ 1 after the mask mandate date (vertical black dash line). Moreover, we have periods of time where the ration is larger than 1 somewhat, so then compliance oN P I might have been negative (in the case where mask ef f = 50%). This seems to indicate that, under our assumptions, the mobility reductions captured via Google indexes and the mask compliance levels have essentially captured the full picture of the pandemic evolution in this region. Similar trends can be seen in Argentina and Nepal (see Figure 8 and Figure 10 ). In the sections above, we have considered a fixed value of mask ef f = 50% across all regions, and we have commented on how a decrease of this value may affect the reductions of R ef f −mobile values. From formula 6 All rights reserved. No reuse allowed without permission. (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 1, 2022. ; https://doi.org/ 10.1101/2022.04.29.22274485 doi: medRxiv preprint we have for each region. Here, we perform a maximum-likelihood estimate of mask ef f for each region, such that the fit of R ef f −mobile to R ef f −data is optimized. In other words, if we considered NPIs to consist only of mask-wearing, what mask efficacy would best reproduce the observed time series of effective reproduction number in a given region? We present our results in a new plot with a small The graphs of the new R ef f −mobilemask time series are presented in Figure 11 (dark pink curves). Here we see again that Sweden stands out simply because mask wearing was not a policy the population had adopted All rights reserved. No reuse allowed without permission. (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 1, 2022. ; https://doi.org/ 10.1101/2022.04.29.22274485 doi: medRxiv preprint in 2020; even at a 100% efficacy, an adoption of 1-2% has negligible effect. Thus other NPI factors must have been in play. Among all the other countries, we see that high levels of mask wearing (for instance in Ontario, and similarly Romania and Argentina) correlates to realistic mask efficacy values (in [54] a realistic expected protection from the mostly cloth-type masks available widely in 2020 is between 30% to 50%) In regions such as Florida, Italy, Indonesia, South Africa, Ghana and Brazil the lower mask wearing levels would have needed to be compensated by higher mask efficacy levels. Since such levels of mask efficacy were not possible in 2020 for an average individual in any of these countries, then other NPI factors had to have been in play. This analysis furthers the importance of our analysis in the last subsection (Subsection 3. 3.2) where we manage to quantify the effect of other NPI factors (oN P I) in reductions of transmission in each of the 12 regions. This paper shows that estimated epidemiological parameters of an underlying SEIR(L) model can be used to compute a time series of effective reproduction numbers accounting for mobility and mask wearing (explicitly using available free data), which are then used to highlight differing pandemic trajectories in very different parts of the world. Our sample consisted of 12 regions (10 countries, 1 Canadian province and 1 US state) chosen based on diversity of population density, median age, urbanization of the population, projected average contact rates, and gross domestic product (GDP) as in Table 1 . We used a SEIR(L) epidemiological model (with parameter values from existing literature) and we used it to compute the near disease-free mathematical expressions of the effective reproduction numbers in terms of initial exponential growth of infection. While we made some specific choices on the structure of the compartments and the flow rates between them, we note that the mathematical methodology we applied is universally applicable to other forms of SIR and SEIR models, and not just to ours. This makes our results applicable to many other variants of compartments models for infectious disease transmission, not just for SARS-Cov-2 (for instance the next concern in the wake of the pandemic are the next flu seasons and the possible mitigation measures in case high;y infectious flu variants will make an appearance). We were able to highlight quantitative relationships between the inferred weekly effective reproduction numbers and the estimated weekly effective numbers based on mobility reduction (as captured by Google mobility index) and mask-wearing (as captured from existing data). Figures 2 and 3 show that there is a sharp drop in R ef f −data at the initial lockdown, then most countries have maintained their effective reproduction numbers around 1 by a combination of reduced mobility, mask-wearing and some additional other NPI's not. Figure 7 shows a mobility-induced decrease of R ef f −mobil with a further decrease when we add mask wearing levels in each region. Further, our modeling analyses provide direct illustrations of the effectiveness of other NPI's in controlling transmission in each region. Figures 7 and 8 show that effective reproduction numbers in all regions have been helped by population adherence and practice of NPI's over and above reductions in mobility and mask protection. This is also obvious by comparing the values of the effective reproductions numbers inferred from data versus (in Figure 7 most regions maintain their values in the neighborhood of 1), while mobility reduction alone (as illustrated from Google mobility reports and in Figure 4 ) has steadily decrease beyond May 2020 in all locations. There are several assumptions underlying our study. Clearly, the mobility reduction as reflected in Google mobility reports is used here as representative across each of the populations, however that may not be quite accurate and it depends on the percentage of cellphone usage in a region and whether or not that percentage can be considered representative of an average individual. At the same time, the pre-pandemic projected contact All rights reserved. No reuse allowed without permission. (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 1, 2022. ; https://doi.org/10.1101/2022.04. 29.22274485 doi: medRxiv preprint Figure 11 : Maximum-likelihood estimate of mask ef f for each region so that the fit of R ef f −mobile to R ef f −data is optimized. Curves in dark pink color are the new R ef f −mobilemask estimates, where for each country we use the deduced value of mask ef f listed. rates from [21] are themselves estimates, thus subject to further change or calibration, given the wealth of data from last year studies. Nevertheless, the overall ideas we followed are fairly straightforward and the quantification of the control on the pandemic via time series of effective reproduction numbers can be done as shown by using parameter values and data, without the need to model and fit SEIR model curves. Even if some of the assumptions/values in our analysis change, the methodology we propose is flexible, novel and easy to follow and implementing any new or updated piece of data available is straightforward. As future directions of research we are interested in All rights reserved. No reuse allowed without permission. (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 1, 2022. ; https://doi.org/10.1101/2022.04. 29.22274485 doi: medRxiv preprint using differing data sources for contact rates in some regions and differing estimates for mobility reductions to replace the Google mobility reports. We can also us differing S(E)IR models and more data on socio-economic and demographic factors that may lead to further insights into accounting for differences in pandemic evolution in diverse countries around the world. The Jacobian matrix near the disease-free equilibrium (DFE, which consists of S(0) = N and I(0) = 0) for the system of equations (2) is: Using the next generation matrix method around the DFE ([34]) we compute R 0 as the largest eigenvalue of the matrix F V −1 and we obtain its closed form expression: Further, using [55] , we can compute the eigenvalues of the reduced Jacobian above and find that there is one positive eigenvalue (responsible for the growth near the DFE) which can be derived in closed form: ρ = − γperasy + κperasy + γ − κ − γ − σ 2 + 2 γ 2 perasy 2 − 2 2 γκperasy 2 + 2 κ 2 perasy 2 − 2 2 γ 2 perasy + 4 2 γκperasy − 2 2 κ 2 perasy + 2 γ 2 − 2 2 γκ + 2 κ 2 4 + 2 γ 2 perasy − 2 γκperasy − 2 γperasyσ + 2 κperasyσ − 2 γ 2 + 2 γκ + 2 γσ − 2 κσ + 4βσ + γ 2 − 2γσ + σ 2 4 , which in turn can be solved for an expression of β as a function of the growth factor ρ near the DFE: β = β(ρ) := γperasyρ + γperasyσ − κperasyρ − κperasyσ − γρ − γσ + κρ + κσ + γρ + γσ + ρ 2 + ρσ σ . Finally we can estimate R 0 as a function of the growth factor near the DFE in each region using 9 as: R 0 (ρ) = γperasyρ + γperasyσ − κperasyρ − κperasyσ − γρ − γσ + κρ + κσ + γρ + γσ + ρ 2 + ρσ σ(( γperasy − κperasy − γ + κ + γ)) . (10) First several months of the pandemic, there was considerable debate on the effect of face masks on limiting the spread of the COVID-19 pandemic and whether to recommend the general public to use a face mask. Later, articles, scientific reports, and data proved the impact of face masks in altering the outcomes of peak hospitalization [5] . Notably, face masks are found to be useful in both preventing asymptomatic transmission and illness in healthy persons. We adapt our previously developed SEIRL model [56] for transmission of COVID-19 with the impact of public use of face masks. Moreover, varying efficacy and compliance of masks have an impact on the transmission dynamics and control of the COVID-19 pandemic [4] . All rights reserved. No reuse allowed without permission. (which was not certified by peer review) is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. The copyright holder for this preprint this version posted May 1, 2022. ; https://doi.org/10.1101/2022.04. 29.22274485 doi: medRxiv preprint A review [57] of observational studies estimates that surgical and comparable cloth masks are 67% effective in protecting the wearer. Some reports show that even a cotton T-shirt can block half of the inhaled aerosols and almost 80% of exhaled aerosols measuring 2µm across e.g. unpublished work by Linsey Marr, an environmental engineer at Virginia Tech in Blacksburg. We consider 50% efficacy in our model. The proportion of a population wearing face masks differs across countries/regions based on social norms, political reasons, the consequences of non-compliance e.g., fines. The results from a study surveying are different for example: • [58] shows 72% of U.S. adults say they either always wear a face mask or wear one often when going to public places. • Percentage of people who worn a face mask outside their home always is reported 93.9% in Italy and 12.1% in Sweden 12.1% by YouGov; Imperial College London [52] . We adapt our SEIRL model with the compliance of mask-wearing value denoted as compl in table 1 for each region. In contact transmissible diseases e.g. COVID-19, contact rate has an important role in epidemic models. Prem et al. [21] provide data-driven contact matrices in the home, work, school, and other locations for 152 countries of the world. We amalgamated the average contact rate in home, work, and other locations, then we computed the weighted average of a given projected contact matrix in a region with the corresponding proportions of 5-year age population groups in 2020 to determine one single average contact rate. Early transmission dynamics in wuhan, china, of novel coronavirus-infected pneumonia The reproductive number of covid-19 is higher compared to sars coronavirus Estimating in real time the efficacy of measures to control emerging communicable diseases A modelling framework to assess the likely effectiveness of facemasks in combination with 'lock-down' in managing the covid-19 pandemic To mask or not to mask: Modeling the potential for face mask use by the general public to curtail the covid-19 pandemic 5 Sars transmission, risk factors, and prevention in hong kong Risk factors for sars among persons without known contact with sars patients Modelling the covid-19 epidemic and implementation of population-wide interventions in italy 26 Early dynamics of transmission and control of covid-19: a mathematical modelling study Targeted adaptive isolation strategy for covid-19 pandemic 5 Safe blues: A method for estimation and control in the fight against covid-19 Effects of social distancing on the spreading of covid-19 inferred from mobile phone data Chopping the tail: how preventing superspreading can help to maintain covid-19 control Monitoring social distancing and sars-cov-2 transmission in brazil using cell phone mobility data Linking excess mortality to mobility data during the first wave of covid-19 in england and wales. SSM -Population Health The relationship between mobility and covid-19 pandemic: Daily evidence from an emerging country by causality analysis Public mobility data enables covid-19 forecasting and management at local and global scales Big data insight on global mobility during the covid-19 pandemic lockdown Mining google and apple mobility data: temporal anatomy for covid-19 social distancing Age-stratified transmission model of covid-19 in ontario with human mobility during pandemic's first wave Projecting social contact matrices in 152 countries using contact surveys and demographic data Inferring high-resolution human mixing patterns for disease modeling Social contacts and mixing patterns relevant to the spread of infectious diseases Quantifying population contact patterns in the united states during the covid-19 pandemic Changes in contact patterns shape the dynamics of the covid-19 outbreak in china Quantifying the impact of physical distance measures on the transmission of COVID-19 in the UK Evolving social contact patterns during the covid-19 crisis in luxembourg The differential impact of physical distancing strategies on social contacts relevant for the spread of covid-19 The effect of control strategies to reduce social mixing on outcomes of the covid-19 epidemic in wuhan, china: a modelling study for the Mathematical Modelling of Infectious Diseases Projecting contact matrices in 177 geographical regions: an update and comparison with empirical data for the covid-19 era A comprehensive estimation of country-level basic reproduction numbers r0 for covid-19: Regime regression can automatically estimate the end of the exponential phase in epidemic data What are the underlying transmission patterns of covid-19 outbreak? an age-specific social contact characterization Quantifying contact patterns in response to covid-19 public health measures in canada running title: Contact patterns during covid-19 in canada Mathematical modelling of covid-19 transmission and mitigation strategies in the population of ontario, canada Estimating the asymptomatic proportion of coronavirus disease 2019 (COVID-19) cases on board the diamond princess cruise ship, yokohama, japan, 2020 Estimation of the basic reproduction number, average incubation time, asymptomatic infection rate, and case fatality rate for COVID-19: Meta-analysis and sensitivity analysis Asymptomatic Infection and Asymptomatic Transmission Quantifying the role of social distancing, personal protection and case detection in mitigating COVID-19 outbreak in ontario Simulating preventative testing of sars-cov-2 in schools: policy implications A path out of covid-19 quarantine: an analysis of policy scenarios Impacts of covid-19 on public transport ridership in sweden: Analysis of ticket validations, sales and passenger counts Nära Var Tredje Bil Borta Från Stockholms Gator Critics Question Swedish Approach as Coronavirus Death Toll Reaches 1,000. The Guardian pid=1710000501&pickMembers%5B0%5D=1.1& pickMembers%5B1%5D=2.1&cubeTimeFrame.startYear=2016&cubeTimeFrame.endYear=2020&referencePeriods=20160101%2C20200101 Accessed The Institute for Health Metrics and Evaluation How Often Have You Worn a Face Mask Outside Your Home to Protect Yourself or Others from Coronavirus (COVID-19 Covid-19 in ontario: Modelling the pandemic by age groups incorporating preventative rapid-testing Covid-19 and use of non-traditional masks: how do various materials compare in reducing the risk of infection for mask wearers Estimating epidemic exponential growth rate and basic reproduction number Age-stratified transmission model of covid-19 in ontario with human mobility Physical distancing, face masks, and eye protection to prevent person-to-person transmission of sars-cov-2 and covid-19: a systematic review and meta-analysis Americans' Face Mask Usage Varies Greatly by Demographics Figure Captions 1 Transmission model in diagram form The yellow curve represents the outcome of equation (10), while the vertical pink dash line represents the start of nationwide social distancing orders in each region under consideration Lightest color patches signify biggest reductions in values of R ef f . In the second panel we see the cumulative incidence for the same regions Baseline rates were computed from Google mobility reports over a 6week interval of In each panel, the red vertical dash line represents the lockdown measures date and the pink curve shows the weekly effective reproduction numbers from the mobility affected by average contact rate in that region 12 7 A comparison of the effective reproduction number as obtained from the incidence data (yellow curves), to our theoretical estimate from Section 3.1 (solid pink curves) using mobility data. Beyond the date of mask mandate enactment in each region, we show the theoreticallyestimated incidence both with (solid pink) and without (dashed pink) the added effect of mask-wearing. The red and black dashed lines in each panel show the lockdown and mask-wearing mandate dates, respectively, in each region as publicly available 10 Effect of other NPI measures in Ontario and Argentina. Mandatory masks were introduced at dates represented by the vertical dashed line in each panel Maximum-likelihood estimate of mask ef f for each region so that the fit of R ef f −mobile to R ef f −data is optimized. Curves in dark pink color are the new R ef f −mobilemask estimates, where for each country we use the deduced value of mask ef f listed Not applicable. The first author's salary was supported from this grant as well. Not applicable. The data used and/or analyzed during the current study are derived from public domain resources. The This study does not require ethical approval. All data used or analyzed in this study are freely available in the public domain, and all methods are obtained directly from cited researchers. There is no direct interaction between the researchers and any individual or group in this study. The authors declare that they have no competing interests. Not applicable. All authors contributed equally to the design and implementation of the research, analysis of the results, and writing of the manuscript.Author details 1 Department of Mathematics & Statistics, University of Guelph, 50 Stone Road E., N1G 2W1, Guelph, Canada. 2 Modeling, Epidemiology and Data Science, Sanofi Pasteur, Toronto, Canada.