Abstract
Dynamic models are used to assess the impact of three types of face masks– cloth masks, surgical/procedure masks and respirators– in controlling the COVID-19 pandemic in the United States. We showed that the pandemic would have failed to establish in the US if a nationwide mask mandate, based on using respirators with moderately-high compliance, had been implemented during the first two months of the pandemic. The other mask types would fail to prevent the pandemic from becoming established. When mask usage compliance is low to moderate, respirators are far more effective in reducing disease burden. Using data from the third wave, we showed that the epidemic could be eliminated in the US if at least 40% of the population consistently wore respirators in public. Surgical masks can also lead to elimination, but requires compliance of at least 55%. Daily COVID-19 mortality could be eliminated in the US by June or July 2021 if 95% of the population opted for either respirators or surgical masks from the beginning of the third wave. We showed that the prospect of effective control or elimination of the pandemic using mask-based strategy is greatly enhanced if combined with other nonpharmaceutical interventions that significantly reduce the baseline community transmission. This study further emphasizes the role of human behavior towards masking on COVID-19 burden, and highlights the urgent need to maintain a healthy stockpile of highly-effective respiratory protection, particularly respirators, to be made available to the general public in times of future outbreaks or pandemics of respiratory diseases that inflict severe public health and socio-economic burden on the population.
Author summary We developed and used dynamic models to assess the role of highly-effective face coverings on the control and mitigation of the COVID-19 pandemic in the US. The study indicates that implementing and sustaining mask mandates is useful in containing diseases like COVID-19. Additionally, the study suggests that prioritizing the use of respirators is more effective in combating the disease than using other mask types. Specifically, the COVID-19 pandemic would have been prevented from being established in the US if four in every five Americans started wearing respirators during the first two months of the pandemic. The study further shows that COVID-19 can be eliminated in the US if a universal masking strategy that emphasizes respirators, requiring only 23% compliance, is combined with other nonpharmaceutical interventions that can reduce community transmission by 20%. Furthermore, the daily COVID-19 death rate can be completely suppressed by June 2021 if 95% of the population consistently use respirators. The elimination will extend to January 2022 if cloth masks were adopted instead. We conclude that stockpiling and distributing highly-efficient face coverings, notably respirators, will be vital in effectively curtailing future epidemics and pandemics of respiratory diseases.
1 Introduction
Since its emergence at the end of December 2019, the novel coronavirus pandemic (COVID-19) has been causing devastating public health and socio-economic burden around the world. The COVID-19 pandemic has, as of March 27, 2021, accounted for over 126.5 million confirmed cases and over 2.77 million reported deaths globally [1, 2]. The United States has disproportionately contributed about 30 million reported cases, accounting for more than 25% of the total global burden, and over 548, 000 reported COVID-19 deaths [2, 3]. In fact, the US Centers for Disease Control and Prevention (CDC) estimated that the true number of COVID-19 cases in the US by the end of 2020 exceeds 80 million. This figure account for the fact that no wide-scale, random testing has been implemented to detect asymptomatic and presymptomatic cases, of which 4.1 million were hospitalized [4]. The estimated economic impact of COVID-19 in the US (measured in terms of aggregated mortality, morbidity, mental health conditions, and direct economic losses) is as high as 16 trillion USD [5–9]. Due to this enormous public health and socioeconomic burden, the use of all possible control and mitigation strategies is necessary to prevent further illness and death.
In response, both pharmaceutical interventions, including vaccines, monoclonal antibodies, and small-molecule drugs, as well as non-pharmaceutical interventions, such as social distancing, hand washing, and mask wearing, have been used to slow the spread and limit the impact of the pandemic. COVID-19 vaccines are now being deployed throughout the world, most widely in high income countries [10]. However, the rollout of these vaccines is being hampered by severe supply chain constraints even among high-income countries, leaving many individuals susceptible to infection and associated morbidity and mortality. This situation is especially pronounced in low- and middle-income countries (LMICs) that face real difficulties in competing for and acquiring scarce vaccine supplies. This disparity is projected to leave up to one quarter of the world’s population without COVID-19 vaccines until at least 2022 [11]. Even among those who ultimately have access, an estimated 32% of the global population is hesitant or unwilling to receive a COVID-19 vaccine, posing a significant challenge to the goal of achieving widespread vaccine-induced immunity against the pandemic [12]. Additionally, COVID-19 variants of concern are continuing to emerge that reduce the effectiveness of existing vaccines [13], which could drive the need for continuous development of reformulated vaccines. As such, while vaccines are a powerful tool, they represent only one tool in the armamentarium of all available control measures against the pandemic.
Of the non-pharmaceutical interventions, masking has been among the most contentious in specific countries, including the United States, due in part to politicization and conflicting early recommendations. Nevertheless, masking remains a crucial intervention, even in communities or populations with widespread vaccine deployment. Especially for those who are hesitant or unwilling to be vaccinated but who are willing to wear a mask, continued masking helps prevent disease transmission in otherwise susceptible individuals. Masking also presents a partial solution to variants of concern. Whereas variants have potential to evade vaccines, masks present a broad-spectrum physical barrier that is relatively insensitive to specific variants. With regard to variants of concern that demonstrate increased transmissibility, some public health experts have suggested that the use of higher-efficacy respiratory protection may be necessary to maintain similar levels of control compared to wild-type SARS-CoV-2.
As masking remains a vital strategy in the control of the COVID-19 pandemic, the efficacy of various mask options available to the general public should be considered in order to optimize protection against disease transmission. Broadly, several categories of respiratory protection exist. These include cloth masks, surgical or procedure masks, and respirators, which are defined as tight-fitting respiratory protective devices that meet or exceed the NIOSH N95 standard [14, 15]. Cloth masks and surgical or procedure masks are considered to be loose fitting face coverings and are designed to be used as source control. Their primary purpose is preventing the spread of disease from infected wearers to those around them, and not the protection of the wearer from others. Due to the unregulated and highly heterogeneous nature of cloth mask design and construction, efficacy estimates vary widely, ranging from 0−50% [14, 15]. Like cloth masks, surgical masks and procedure masks are also designed to be loose fitting and act as source control but are made from materials that have a theoretical filtration efficiency of up to 85% [16]. Due to their lack of a tight seal, the real world filtration efficiency of the surgical or procedure masks can be much lower [15]. In contrast to cloth and surgical masks, most respirators are designed to fit tightly to the face and represent true respiratory protection, protecting the wearer from respiratory hazards around them. Estimates for the efficiency of respirators filtering facepiece respirators are nearly 100% for charged biological particles such as respiratory aerosols [14–16]. In all cases, poor fit can impair efficiency, leading to decreased real-world fitted filtration efficiency ratings that are lower than theoretical maximums as determined solely by testing of the filtration material [17]. This effect is especially important when considering the differences between loose-fitting (cloth and surgical mask) and tight-fitting (respirator) devices. Loose-fitting devices are not designed to form a seal and thus enable significant filter bypass. This results in very low real-world fit factors in the low single-digits [17]. In comparison, respirators are able to form a tight seal and when properly fitted, are able to achieve fit factors at and above 100, which represents a 100-fold reduction in the concentration of particles inside the device compared to particles outside of the device [18]. Due to the highly-protective nature of properly fitting respirators, healthcare workers have used these devices for protection from SARS-CoV-2 throughout the pandemic.
Since the onset of the pandemic, the prevailing public health recommendations with regard to masking have been that cloth and surgical masks should be used by members of the general public partly due to being widely available. Recently, however, calls have been made to the CDC by public health and respiratory protection experts to change their public health guidance and to begin recommending the use of higher-efficacy respiratory protection, such as respirators, for individuals outside of healthcare roles [19]. These individuals could include those at high risk of acquiring SARS-CoV-2 infection or from populations that have suffered disproportionately, including essential workers such as those at grocery stores and meatpacking plants; people of specific races and ethnic background, including Black, Indigenous, Latino, and other people of color who have been heavily impacted by the pandemic; older individuals; and those with underlying medical conditions placing them at increased risk of severe illness and death, including cancer, chronic kidney disease, chronic obstructive pulmonary disease (COPD), Down syndrome, heart disease, transplant recipients, obesity, pregnancy, sickle cell disease, smoking, and diabetes [20]. The rationale for doing so stems from the theoretical improved protection that respirators could provide to those at high risk. As vaccination for many of these individuals is still months away even in high-income countries, and up to a year or more away for those in LMICs, the recommendation to use improved respiratory protection could play a role in lowering COVID-19 disease transmission and associated morbidity and mortality. As current variants of concern circulate and new concerning variants are likely to arise in the future, maximizing the available non-pharmaceutical interventions could also serve in a vaccine-sparing and therapeutic-sparing capacity by further limiting the ability of variants to emerge due to selective pressures against specific therapeutic monoclonal antibodies and vaccines.
The impact of the use of respiratory protection with efficacy far exceeding that of cloth and surgical masks for members of the general public is currently not known. Previous modeling studies have focused on the role of lower-efficacy face coverings, such as cloth and surgical masks, when adopted by a large majority of the population [14, 21–23]. What is not clear is the potential role that highly-efficacious respiratory protection could serve outside of healthcare and what levels of use in the population may be necessary to project significant changes compared with strategies employing lower-efficacy masks. As populations use a combination of varying categories of face coverings, it is also not known how combinations of different types of respiratory protection could impact disease transmission. Using mechanistic mathematical modeling, we studied the role of high-efficacy respiratory protection when compared with lower-efficacy solutions commonly suggested through public health messaging. This study clarifies the role that high-efficacy respiratory protection, such as respirators, could serve in decreasing COVID-19 disease transmission and speeding an end to the pandemic, as well as highlighting the implications of potential new strategies to fight future epidemics and pandemics.
2 Methods and Results
2.1 Formulation of Basic Model
We start with formulating a basic mathematical model for the transmission dynamics of COVID-19 in a population where a certain proportion habitually wear face masks (cloth, surgical or respirators) in public. The basic model considers the total population at time t, denoted by N(t), as a single group that is sub-divided into epidemiological compartments of individuals who are susceptible (S(t)), exposed (E(t)), pre-symptomatically-infectious (P(t)), symptomatically-infectious (I(t)), asymptomatically-infectious (A(t)), hospitalized (H(t)), and recovered R(t), so that N(t) = S(t) + E(t) + P (t) + I(t) + A(t) + H(t) + R(t). “Exposed” individuals are those who are already infected with the disease but are not yet infectious (i.e., they have latent infection).
The basic model is given by the following deterministic system of nonlinear differential equations for the rate of change of each of the epidemiological compartments (where a dot represents differentiation with respect to time):
In the basic model (2.1), β is the overall community transmission rate, 0 < cm < 1 is the proportion of individuals who wear faces masks in public, 0 < εm < 1 is the protective efficacy of face masks to prevent the transmission of infection to a susceptible individual and ηj (with j ∈ {P, I, A, H}) represents the modification parameter for the infectiousness of the pre-symptomatic (ηP), symptomatic (ηI), asymptomatic (ηA), and hospitalized (ηH) infectious individuals. The parameter, σE, is the progression rate of exposed (i.e., newly-infected individuals but not yet infectious) individuals to the pre-symptomatic class; σP is the rate at which pre-symptomatic individuals progress to either the symptomatic class or the asymptomatic class at the end of the pre-symptomatic period (i.e., after surviving the incubation period), with 0 < r < 1 as the proportion of the pre-symptomatic individuals who show clinical symptoms of the disease at the end of the pre-symptomatic period (the remaining proportion, 1 − r, do not show clinical symptoms at the end of the pre-symptomatic period, and are moved to the asymptomatically-infectious class). Symptomatic infectious individuals are hospitalized at a rate ϕI. Individuals in the symptomatic, asymptomatic and hospitalized infectious classes recover at a rate γI, γA and γH, respectively. Symptomatic and hospitalized individuals suffer disease-induced mortality at a rate δI and δH, respectively. The basic model forms the foundation for our estimation of unknown parameters, particularly the overall community transmission rate (β), the mortality rates (δI and δH), and the recovery rates (γI, γA and γH). The state variables and parameters of the basic model (2.1) are described in Table S1 of the Supplementary Material.
In the formulation of the basic model (2.1), it is assumed that the total population is well-mixed (i.e., every member of the community is equally likely to mix, and acquire infection, from every other member of the community) and the average waiting time in each epidemiological compartment is exponentially distributed [24]. It is also assumed that hospitalized individuals can transmit infection (at the rate ηHβ; this assumption can be relaxed, if relevant, by setting ηH = 0). Recovery from COVID-19 is assumed to confer permanent immunity [21–23, 25, 26]. Although a number of COVID-19 variants have been identified in Brazil, South Africa, UK and the US [27–30], they are not included in this model due to insufficient data on how they might interact with the predominant strain. We define a parameter, 0 ≤cr < 1, is the proportion of individuals in the community who observe other additional forms of non-pharmaceutical interventions for reducing commu-nity transmission, such as through social-distancing. The control reproduction number of the basic model (2.1), denoted by ℝc, is given by [31, 32] ℝc = ℝP + ℝI + ℝA + ℝH, where,and
represent, respectively, the constituent control reproduction numbers associated with disease transmission by the pre-symptomatic infec-tious (ℝP), symptomatically-infectious (ℝI), asymptomatically-infectious (ℝA), and hospitalized individuals (ℝH) [21–23]. The control reproduction number (ℝc) represents the average number of new COVID-19 cases generated by a typical infected individual introduced into a population where a certain proportion (0 < cm < 1) habitually wears face masks in public. The implication of this threshold quantity is that an outbreak will occur in the community if the value of ℝc exceeds unity, and no significant outbreak occurs otherwise [22, 23].
2.1.1 Calibration and parameter estimation of basic model
The basic model (2.1) contains 16 parameters, the values of 10 of which (cm, εm, ηP, ηI, ηA, ηH, σE, σP, r and ϕI) are known from the literature, while the values of the remaining six (β, γI, γA, γH, δI and δH) are unknown and will be obtained by fitting the basic model with observed data. Specifically, the basic model (2.1) will be fitted with real data to obtain realistic estimate of the associated unknown parameters. Specifically, we use the cumulative COVID-19 mortality data for the US (since this is more reliable than case data [22]) for three different periods (January 22 to June 30, 2020, July 1 to October 11, 2020, and October 12, 2020 to March 11, 2021), corresponding to the three observed COVID-19 pandemics waves in the US, to fit the single group model (2.1). In order to overcome the issue of avoidable errors introduced by fitting cumulative data to deterministic models [33], we also use new daily mortality data for the same period to check the quality of our fit. Based on the fitting, we obtain estimates for the community transmission rate (β), the mortality rates for symptomatic infectious and hospitalized individuals (δI and δH, respectively), and the recovery rates of symptomatic infectious, asymptomatic infectious, and hospitalized individuals (γI, γA, and γH, respectively) for each of the three pandemic waves. The fitting is achieved by using a simple nonlinear regression process, in which the sum of the squared differences between the observed cumulative mortality data for the US and the cumulative mortality from the model output (i.e., dt, where T is the time duration for the simulations) is minimized and a bootstrapping approach is used to compute the 95% credible intervals for the estimated parameters [34–36].
Time series plots illustrating the fit, and how well the output of the model matches the daily data for the same time period, are depicted in Figures 1 (a) and 1 (b), respectively. The estimated values of the parameters obtained from fitting the basic model are tabulated in Table 2.
(a) Time series plot of the observed/reported cumulative mortality data for the US (red dots) and predicted cumulative mortality for the US generated from the basic model (2.1). (b) Simulations of the basic model (2.1) using the fixed parameter values in Table 1 and the estimated parameter values in Table 2, illustrating the three COVID-19 waves observed in the US during the periods from January 20, 2020 to June 30, 2020 (first pandemic wave), July 1, 2020 to October 11, 2020 (second pandemic wave), and October 12, 2020 to March 11, 2021 (third pandemic wave). Dashed vertical blue lines demarcate the pandemic waves.
2.1.2 Results from the basic model
The basic model (2.1) is now simulated to assess the population-level impact of face masks usage in the community. For simulation purposes, we consider four categories of face masks, namely cloth mask (with estimated efficacy 30%), improved cloth masks or poorly-fitted surgical mask (with estimated efficacy of 50%), properly-fitted surgi-cal mask (with estimated efficacy of 70% [15, 16, 38]), and properly-fitted respirators (with estimated efficacy of 95%) [16, 38]. The impact of combining the face masks use strategy with an additional adoption of NPIs that target reduction in community transmission of COVID-19 (such as social-distancing, self-isolation, avoiding travels and/or public gatherings, avoiding indoor dining in restaurants, etc.), as measured by the parameter cr, will also be assessed. The simulations will be carried out using the fixed baseline parameter values presented in Table 1 and the estimated baseline parameter values in presented in Table 2.
Fixed baseline parameter values of the basic model (2.1). Although the baseline value of the mask compliance parameter (cm) is set at 0.50, we use different baseline values for the first and second pandemic waves. In particular, we set the baseline values for the first and second waves to cm ≈ 0.20 and cm ≈ 0.35, respectively.
Estimated (fitted) baseline parameter values and confidence intervals (CIs) for the basic model (2.1) using COVID-19 mortality data for the US for the period from (a) January 20, 2020 to June 30, 2020 (first pandemic wave), (b) July 1, 2020 to October 11, 2020 (second pandemic wave), and (c) October 12, 2020 to March 11, 2021 (third pandemic wave). The unit of each of the parameters is per day.
Percentage increase (denoted by a “+” sign) or reduction (denoted by a “−” sign) in cumulative mortality by April 30, 2021 (in comparison to the 570, 000 projected baseline value), if only cloth, improved cloth, surgical, or respirators (or equivalent) mask is used in the US, for various levels of reduction in community transmission due to the implementation of other NPIs implemented at the beginning of the third wave (as depicted in Figure 4).
The first set of simulations we carried out are for assessing the impact of early implementation of mask use strategy in the US Specifically, we simulated the scenario where only one of the aforementioned four categories of face masks is adopted in the community. For these simulations, we used the baseline values of the fixed and estimated parameters corresponding to the first wave of the pandemic (tabulated in Tables 1 and 2). The results obtained, depicted in Figure 2, show that the early implementation of universal face mask use strategy, with compliance of 80%, resulted in a dramatic reduction in the burden of the pandemic, particularly if high quality masks, notably surgical masks and respirators, are used. In particular, if 80% of Americans were wearing respirators starting from March 15, 2020 (which was almost 2 months after the index case on January 20, 2020, or four days after the WHO declared COVID-19 to be a global pandemic), up to 84% of the nearly 128, 000 baseline cumulative deaths recorded by June 30, 2020 would have been prevented (Figure 2 (a), green curve). On the other hand, if surgical masks were chosen instead, the reduction in the baseline cumulative mortality by June 30, 2020 would have been 43% (Figure 2(a), cyan curve). The reductions corresponding to the use of improved cloth mask or basic cloth masks only are 18% (Figure 2(a), magenta curve), or 5% (Figure 2(a), gold curve), respectively.
Simulations of the basic model (2.1) to assess the impact of wide adoption of face mask wearing policy (assuming 80% compliance) on the cumulative COVID-19 mortality ((a)-(c)) and daily COVID-19 mortality ((e)-(f)) in the US, starting from the onset of the first pandemic wave (January 20, 2020). Scenarios include cloth masks only (with estimated efficacy of 30%), improved cloth masks only (with estimated efficacy of 50%), surgical masks only (with estimated efficacy of 70%), or respirators only (with efficacy of 95%), The other parameter values used in the simulations are given in Table 1 and the first wave parameter values in Table 2 (a). The dotted vertical blue line depict the onset of the face mask use strategy, while the dashed vertical blue line depicts the end of the first pandemic wave: (a) and (d): March 15, 2020 (i.e., four days after the World Health Organization declared COVID-19 to be a global pandemic), (b) and (e): February 18, 2020 (i.e., one month after the US index case), and (c) and (f): February 2, 2020 (i.e., two weeks after the US index case).
Much higher reduction in the baseline cumulative mortality will be recorded if the universal masking strategy was implemented earlier than March 15, 2020. For example, if the universal masking strategy was started on February 18, 2020 (i.e., a month after the first reported case in the US), the pandemic would have hardly taken off if the (80%) masks-wearing proportion chose to wear respirators only (Figure 2 (b), green curve). On the other hand, choosing to wear cloth mask only or improved cloth mask only or (even) surgical masks only would have resulted in a sizable epidemic (Figure 2 (b), gold, magenta and cyan curves). Similar results are obtained if the mask use strategy was implemented two weeks after the index case in the US (Figure 2 (c)). Plots for daily mortality are depicted in Figures 2 (d)-(f), showing similar trends as the plots for the cumulative mortality (2 (a)-(c)). In all of these simulations, the COVID-19 pandemic will not have taken off in the US if the assumed 80% of the masks-wearing population chose to wear respirators (or equivalent) only. Furthermore, these simulations show that choosing properly-fitted surgical masks would also have made significant impact in reducing COVID-19 mortality, but certainly no where as effective as the case where respirators was preferred. The respirators and surgical masks are, of course, far more effective than using the improved or regular cloth mask.
Additional simulations of the basic model (2.1) are carried out to assess the combined impact of face masks usage and other NPIs. It should also be emphasized that, since the model (2.1) was parametrized using COVID-19 data from the onset of the outbreak in the US (i.e., using data from January 20, 2020 to mid March 2021), the combined effects of other NPIs that resulted in a reduction in community transmission of COVID-19 (such as social-distancing, travel restrictions, avoiding large gatherings and indoor dining in restaurants) are already embedded into the results/data. That is, our model parameter calibration includes some baseline level of the other
NPIs implemented in the community. To assess the impact of the additional levels of these other NPIs on the burden of the pandemic (from their baseline levels), we multiply the effective community transmission rate (β) by the fraction 1 − cr, where 0 < cr ≤ 1 is a measure of the additional reduction in the community transmission rate. Since our predictions are based on the estimated parameters for the third pandemic wave (i.e., the parameter values estimated using the data from October 12, 2020 to mid March 2021), the additional reduction in community transmission due to the implementation of other NPIs (modeled by the parameter cr) starts on October 12, 2020 and is maintained through the rest of the simulation (or prediction) period (for these, and future simulations in this paper, we used the fixed and estimated values of the parameters corresponding to the third wave of the pandemic). Figure 3 depicts a contour plot of the control reproduction number of the model (ℝc) as a function of face masks compliance (cm) and efficacy (εm) for various levels of increases in baseline values of reduction in community transmission due to other NPIs (cr). The figure indicates a maximum control reproduction number during the simulation period (October 12, 2020 to March 11, 2021) of ℝc = 1.6 (with a baseline value of ℝc = 1.2; signifying a continuing trajectory of disease spread in the US, since ℝc > 1 during the third wave. Further, this figure shows that ℝc decreases with increasing values of face masks compliance and efficacy. In the absence of additional reduction in community transmission due to the implementation of other NPIs (i.e., cr = 0), our simulations (Figure 3(a)) show that if individuals in the community choose to use cloth masks only (with estimated efficacy of 30%), the control reproduction number can never be brought to a value less than unity (hence, the disease will continue to spread in the population).
Contour plot of the control reproduction number (ℝc) of the basic model (2.1), as a function of the face masks compliance (cm) and efficacy (εm). (a) No reduction in community transmission from baseline due to the implementation of other NPI’s (i.e., cr = 0). (b) 10% reduction in community transmission due to implementation of other NPIs (i.e., cr = 0.1). (c) 20% reduction in community transmission due to the implementation of other NPIs (i.e., cr = 0.2). The other parameter values used in the simulations are as given in Tables 1 and 2 (c).
On the other hand, if the American populace choose to use improved cloth masks (with estimated efficacy of 50%) only, a compliance of at least 75% will be needed to bring the control reproduction number to a value less than one. The compliance level needed reduces to 55% if Americans choose, instead, to wear well-fitted surgical masks only. Finally, if respirators were adopted right from the beginning of the pandemic, this figure shows that only 40% compliance will be required to reduce the control reproduction number below unity. In other words, the pandemic can be effectively controlled (eliminated) if 40% of the American populace were wearing face masks of respirators quality right from the beginning of the pandemic outbreak. Furthermore, if the face masks use strategy is combined with other NPIs (i.e., if cr ≠ 0), our simulations show that the face masks compliance level needed to bring the control reproduction number to a value less than unity decreases significantly. For instance, if a 10% increase in reduction in community transmission due to other NPIs (from the baseline) can be achieved (i.e., if cr = 0.1), our simulations show that the reproduction number can be brought to a value less than one using the low-quality cloth masks if compliance is near 100% (Figure 3 (b)). For this scenario, the use of improved cloth masks only would require compliance of 60%, while the use of surgical masks and respirators (or equivalent) only would require compliance of only 45% and 32%, respectively, to bring the control reproduction number below one. For the case where the implementation of other NPIs can lead to a decrease in community transmission by 20%, our simulations show that the use of low-quality cloth masks can bring the reproduction number to a value less than one if at least 75% of Americans wear them (Figure 3 (c)). The required compliance reduces to 45% if improved cloth masks are chosen instead. On the other hand, only 32% and 23% compliance will be needed if Americans opted for surgical and respirators, respectively.
Additional simulations of the basic model are carried out to assess the impact of the face masks strategy on cumulative mortality for various levels of reduction in community transmission due to the implementation of other NPIs. Specifically, we simulated the scenario for COVID-19 dynamics in the US, starting from October 12, 2020, where only cloth, improved cloth masks, surgical masks, or respirators are used. We consider April 30, 2021 (the end of the first 100 days of the new administration in the US) as a reference point for projections. Figure 4 shows that, based on the baseline values of the fitted and fixed parameters in Tables 1-2, about 570, 000 cumulative deaths would be occur in the US by April 30, 2021 if only improved cloth masks are used (Figure 4 (a)-(d), blue curve). If implementing only a cloth mask strategy along with other NPIs that can reduce community transmission during the beginning of the third wave (October 12, 2020) by only 5% (respectively, 10%) (i.e., cr = 0.05 (respectively, cr = 0.10)), about 112, 000 (respectively, 32, 000) deaths above the projected 570, 000 baseline cumulative deaths by April 30, 2021 could be recorded (Figure 4 (a), magenta curve) (respectively, Figure 4 (a), gold curve). This represents approximately a 20% (respectively, 6%) increase in the projected baseline cumulative mortality by April 30, 2021. If the reduction in community transmission is increased to 15% (i.e., cr = 0.15), then about 64, 000 (11%) of the projected 570, 000 cumulative mortality by April 30, 2021 would have been averted (Figure 4 (a), cyan curve). About 30% (i.e., 170, 000) of the projected mortality by April 30, 2021 would be prevented if NPIs could reduce community transmission by 20% during the beginning of the third wave (Figure 4 (a), green curve). On the other hand, if a combination of respirators (or equivalent) are used, the projected cumulative mortality by April 30, 2021 reduces dramatically (Figure 4 (d)). Here, a 5% reduction in community transmission due to other NPIs would avert about 59% of the projected mortality (Figure 4 (b), magenta curve). A 20% reduction in community transmission (cr = 0.2) would avert about 60% of the projected mortality (Figure 4 (b), green curve). Similar results are obtained with respect to daily mortality (Figure S2 in the Supplementary Material).
Simulations of the basic model (2.1) to assess the combined impact of face mask use strategy and reduction in community transmission due to other NPIs (implemented at the beginning of the third wave) on cumulative COVID-19 mortality in the US (a) Only cloth masks are used. (b) Only improved cloth masks are used. (c) Only surgical masks are used. (d) Only respirator (or equivalent) are used. For the blue curves (baseline scenario) in Figures (a)-(d), the mask efficacy is maintained at 50% (i.e., εm = 0.5). The other parameter values used in the simulations are given in Tables 1 and 2 (c). The shaded magenta and green regions represent the beginning and end of the first and second pandemic waves in the US, respectively. The shaded yellow region represents the beginning of the third pandemic wave until March 11, 2021. Dashed vertical blue lines demarcate the three waves.
Additional simulations are carried out to assess the combined impact of using different kinds of face masks (cloth masks with an estimated efficacy of 30%, improved cloth masks with an estimated efficacy of 50%, surgical masks with an estimated efficacy of 70%, and respirators or equivalent with efficacy of approximately 95%) and other NPIs on the cumulative COVID-19 mortality. For these simulations, we choose April 30, 2021 as the reference time for comparison purposes. The results obtained, depicted in Figure 5, show that respirators are far more effective in curtailing COVID-19 mortality than using any of the other mask types, as expected. Further, the effectiveness of each mask type is enhanced if the face masks intervention is combined with increase in baseline levels of other NPIs. For instance, if Americans use cloth masks only with a compliance of 50%, from October 12, 2020, 750, 000 cumulative mortality will be recorded by April 30, 2021 (Figure 5 (a), magenta curve). This represents a 32% increase in the cumulative mortality recorded under the baseline scenario (570, 000), where the face masks compliance (cm) and increase in other NPIs (cr) are maintained at their baseline levels (Figure 5 (a), blue curve). On the other hand, if surgical masks are used from October 12, 2020, the cumulative mortality to be recorded on April 30, 2021 will be 335, 000 (Figure 5 (a), cyan curve), which represents a 41% decrease from the baseline scenario (depicted in Figure 5 (a), blue curve). Furthermore, if respirators are used from October 12, 2020, the cumulative mortality to be recorded on April 30, 2021 will be 240, 000 (Figure 5 (a), green curve), which represents a 58% decrease from the baseline scenario (depicted in Figure 5 (a), blue curve). Improved cloth masks correspond to the baseline scenario (Figure 5 (a), blue curve). If the face mask compliance (cm) is increased by 20% (i.e., if cm is now set to 0.7, instead of 0.5 above), then using cloth masks only from October 12, 2020 will result in a 12% increase in cumulative mortality by April 30, 2021 (Figure 5 (c), magenta curve). For this scenario, surgical masks will reduce the baseline cumulative mortality by 58% (Figure 5 (c), cyan curve), while respirators reduce the mortality by 60% (Figure 5 (c), green curve).
Simulations of the basic model (2.1) to assess the impact of masks type (i.e., cloth mask only, improved cloth masks only, surgical masks only, and respirators only), as a function of time, for various levels of face masks compliance (cm) and increasing levels of other NPIs from their baseline (cr), starting October 12, 2021. (a) various mask types with compliance set at 50% (i.e., cm = 0.5) and other NPIs maintained at baseline (i.e., cr = 0). (b) various mask types with compliance set at 60% (i.e., cm = 0.6) and other NPIs maintained at baseline (i.e., cr = 0). (c) various mask types with compliance set at 70% (i.e., cm = 0.6) and other NPIs maintained at baseline (i.e., cr = 0). (d) various mask types with compliance set at 50% (i.e., cm = 0.5) and other NPIs increased to 10% from their baseline (i.e., cr = 0.1). (e) various mask types with compliance set at 50% (i.e., cm = 0.5) and other NPIs increased to 20% from their baseline (i.e., cr = 0.2). (f) various face mask types with compliance set at 70% (i.e., cr = 0.7) and other NPIs increased to 20% from their baseline (i.e., cr = 0.2). Red dots represent observed cumulative mortality data. The other parameter values used in the simulations are as given in Table 1 and the third wave parameter values in Table 2 (c). The shaded magenta and green regions represent the beginning and end of the first and second pandemic waves in the US, respectively. The shaded yellow region represent the beginning of the third pandemic wave until March 11, 2021. Dashed vertical blue lines demarcate the three waves.
We further simulated the scenario where face mask compliance is kept at baseline (i.e., cm = 0.5) and NPIs increased by 10% and 20% from their baseline (i.e., cr = 0.1 and cr = 0.2), starting from October 12, 2020. In this scenario, cloth masks increase baseline cumulative mortality recorded by April 30, 2021 by 6%, when cr = 0.1 (Figure 5 (d), magenta) and reduce cumulative mortality by 30% when cr = 0.2 (Figure 5 (e), magenta). When cr = 0.1, surgical masks (Figure 5 (d), cyan curve) and respirators (Figure 5 (d), green curve) decrease the cumulative mortality by 54% and 59%, respectively. These percentage reductions due to surgical masks and respirators use increase to 58% (Figure 5 (e), cyan curve) and 60% 5 (e), green curve), respectively, if cr = 0.2. Finally, if face masks compliance is increased to 70% (i.e., cm = 0.7), while NPIs remain at 20% above their baseline levels (i.e., cr = 0.2), the use of cloth mask can reduce the baseline cumulative mortality recorded on April 30, 2021 by 46% (Figure 5 (f), magenta curve). Here, surgical masks (Figure 5 (f), cyan curve) and respirators (Figure 5 (f), green curve) reduce the baseline cumulative mortality recorded by April 30, 2021 by 60% and 61%, respectively.
2.2 Formulation of Extended Model: Effect of Behavior Changes With Respect to Masks Usage
Although the basic model (2.1) provides insights on the population-level impact of various types of face masks implemented as a sole intervention, or in combination with other NPIs, it does not allow for the assessment of the impact of changes in face mask adherence. For example, upon taking the oath of office in the middle of the 3rd wave on January 20, 2021, President Biden issued a mask mandate in all Federal buildings and public transportation systems [45, 46]. To account for changes in face mask usage, the basic model (2.1) is now extended and reformulated into a multigroup structure, which entails stratifying the total population into three groups: those who do not wear masks (Group 1); those who habitually wear other types of masks, except respirators (or equivalent), such as cloth masks, improved cloth/poorly-fitted surgical masks, or properly fitted surgical masks (Group 2); and those who habitually wear respirator or equivalent (Group 3) only in public. The resulting extended multigroup model also allows for back-and forth transitions between the three groups.
To develop the extended model, the total population at time t, denoted by N(t), is split into the total population of individuals who do not wear face masks in public (denoted by N1(t)), the total population of individuals who habitually wear other types of masks except respirators (or equivalent) only (denoted by N2(t)) and the total population of individuals who habitually wear the respirator (or equivalent) only (denoted by N3(t)). Each of these subgroups is further subdivided into seven mutually-exclusive epidemiological compartments of susceptible (S), exposed/latent (E), pre-symptomatic infectious (P), symptomatic infectious (I), asymptomatic infectious (A), hospitalized (H), and recovered R compartments, so that
The extended multigroup model also allows for the assessment of the impact of other NPIs implemented during the beginning of the third wave of the pandemic (as represented by the parameter cr, discussed earlier). The multigroup extended model is formulated via the same approach used in the formulation of the (single group) basic model (2.1). The equations for the rate of change of the state variables of the extended model are given by Equations (S.1)-(S.3) in the Supplementary Material. Parameters that represent change in face mask usage (i.e., the back-and-forth transitions between groups) are denoted by αij; i, j = 1, 2, 3 with i ≠ j, as presented in Table S2 of the Supplementary Material. Further, the state variables and parameters of the extended model are tabulated in Tables S3-S4, respectively). The remaining parameters of the extended model are the same as those used for the basic model (except for the indices, indicating the corresponding group the parameter is defined for). The extended model is also fitted using the cumulative mortality data for the US, in the Supplementary Material.
2.2.1 Results from the extended model
The extended model, given by Equations (S.1)-(S.3) in Supplementary Material, is now simulated to assess the population-level impact of change of masking behaviour (measured in terms of the parameters αij, i, j = 1, 2, 3 with i ≠ j) on the transmission dynamics and control of the disease, for various increases in the baseline coverage levels of other NPIs (cr) implemented after October 12, 2020. The results obtained, depicted in Figure 6, show that increases in the rate of behavior change of non-wearers to adopt to habitually wear respirators (α13) is more effective, in curtailing peak daily mortality and peak daily hospitalization, than corresponding increases in the rate at which non-wearers decide to wear surgical masks (α12). For instance, if 500, 000 non-wearers change their behaviour and begin to wear surgical masks only (i.e., α12 = 0.0015 per day, and α13 = 0), the baseline maximum peak daily mortality obtained when wearers do not abandon mask wearing from October 12, 2020 onward (5, 316 deaths, obtained from Figure 6 (a)) reduces to 3, 961 (which represents a 25.5% reduction in the baseline maximum peak daily mortality). However, if the 500, 000 non-wearers opt to wear respirators or equivalents, instead (so that, α13 = 0.0015 and α12 = 0), the baseline maximum peak daily mortality reduces to 3, 862 (corresponding to 27.4% reduction in the baseline maximum peak daily mortality). On the other hand, if the 500, 000 non-wearers choose to split 50-50 between the two masks (i.e., half wear surgical masks only and the remaining half wear respirators only), the baseline maximum peak daily mortality reduces to 3, 920 (which represents a 26.3% reduction in the baseline maximum peak daily mortality). It should be mentioned that lower reductions in peak daily mortality and hospitalizations are obtained if the comparison is between respirators and cloth or improved cloth masks, or between surgical masks and and cloth or improved cloth masks.
Heatmaps of the extended model (S.1)-(S.3), showing the peak daily mortality ((a)-(c)) and peak daily hospitalizations ((d)-(f)), as a function of the daily rate at which non-wearers opt to wear surgical masks (α12) or respirators (α13) for increasing baseline coverage levels of other NPIs (cr), starting October 12, 2020. (a) cr = 0, (b) c2 = 0.05, (c) cr = 0.1, (d) cr = 0, (e) cr = 0, 05 and (f) cr = 0.1. The values of the other parameters of the extended model used in these simulations are as given in Tables 1 and S6 in the Supplementary Material.
Figure 6 also shows that greater reductions in the peak daily mortality and the peak daily hospitalizations are achieved if the face masks use strategy is complemented with increasing baseline coverages of other NPIs (Figures 6 (c) and (f), in comparison to Figures 6 (a), (b), (d) and (e)). For instance, for the above case where 500, 000 non-wearers change their behavior and opt to wear one of the two masks types, combined with a 10% increase in baseline coverage of other NPIs (cr = 0.1), the maximum baseline daily mortality (5, 316) reduces to 2, 257 and 2, 183, respectively, if the non-wearers choose to wear surgical masks or respirators exclusively (6 (c)). This represents a 57.5% and 58.9% reduction in the baseline maximum peak mortality, respectively (recall that the corresponding reductions, under this scenario but with cr = 0, described above, were, respectively, 25.5% and 27.4%). Additionally, if the 500, 000 non-wearers are split equally between new surgical masks and respirator wearers, then the maximum baseline daily mortality (5, 316) reduces to 2, 233 (representing a 58% reduction in the baseline maximum peak mortality). Similar results are shown with respect to the baseline maximum daily hospitalization in Figures 6 (d)-(f).
We further simulated the extended model to compare the population-level effectiveness of the four mask types (cloth, improved cloth or poorly fitted surgical, properly fitted surgical, and respirators), for the case where all mask-wearers in the community choose to wear only one mask type, for various initial sizes of the total population of mask wearers (denoted by N2(0) + N3(0)). The results are depicted in Figure 7. This figure shows a devastating level of cumulative mortality (about 1.15 million deaths by April 30, 2021) under the worst-case scenario, where no non-wearer in the community begins to wear a face mask in public from the beginning of the third wave, October 12, 2020 (Figures 7 (a)-(c), purple curves). Similar devastating numbers are recorded for the daily mortality (with a peak daily mortality of about 9, 400 on January 4, 2021), under this scenario (Figures 7 (d)-(f)). For the case when 65% of the population were wearing face masks at the beginning of the third wave of the pandemic (i.e., N2(0) + N3(0) = 65% of the population), the cumulative mortality recorded on April 30, 2021 under the baseline scenario (i.e., using the baseline parameter values in Tables 1 and S5-S6 in the Supplementary Material), is 585, 000 (Figure 7 (a), blue curve). This represents a 49% reduction in cumulative mortality under the worst-case scenario.
Simulations of the extended model (S.1)-(S.3), given in Supplementary Material, showing the effects of mask types and initial size of total population of mask wearers (N2(0)+N3(0)) on the cumulative ((a)-(c) and daily ((d)-(f)) COVID-19 mortality. (a) and (d): total initial size of mask wearers is 65%; (b) and (e): total initial size of mask wearers is 80%; (c) and (f): total initial size of mask wearers is 95%. The values of the other parameters of the extended model used in the simulations are as given in Tables 1 and S5-S6 in the Supplementary Material. The shaded magenta and green regions represent the beginning and end of the first and second pandemic waves in the US, respectively. The shaded yellow region represents the beginning of the third pandemic wave until March 11, 2021. Dashed vertical blue lines demarcate the three waves.
If all mask wearers in the community opt for cloth masks only, from the beginning of the third wave of the pandemic (i.e., from October 12, 2020), the cumulative mortality recorded by April 30, 2021 is 650, 000 Figure 7 (a), magenta curve). Although this represents a 43% reduction in the worst-case cumulative mortality, it is higher than the cumulative mortality recorded under the baseline scenario. Using improved cloth or improperly fitted surgical masks only, on the other hand, results in cumulative mortality of 561, 000 by April 30, 2021 (Figure 7 (a), yellow curve), representing a 51% decrease from the worst-case cumulative mortality and a 4% decrease from the baseline cumulative mortality. For this scenario, there is a threshold time (around early September 2021) and cumulative mortality (about 612, 000) after which the recorded cumulative mortality with improved cloth mask exceeds the baseline cumulative mortality. (compare blue and yellow curves in Figure 7 (a)). Similarly, using properly fitted surgical masks only result in cumulative mortality of 488, 000 by April 30, 2021 (Figure 7 (a), cyan curve), representing a 58% decrease from the cumulative mortality recorded under the worst-case scenario and a 17% decrease from the baseline. However, if respirators are used, the projected cumulative mortality by April 30, 2021 is 418, 000 (Figure 7 (a), green curve), corresponding to a 64% decrease from the worst-case cumulative mortality scenario and a 29% decrease from the baseline cumulative mortality. Similar results are obtained with respect to the daily mortality (Figure 7 (d)).
Figure 7 further shows a marked decrease in cumulative mortality with increasing initial size of the population of individuals who opt to wear face masks from the beginning of the third wave of the pandemic (Figure 7 (b)). For instance, if the initial size of mask wearers is increased from 65% to 80%, the reductions in worst-case cumu-lative mortality is much higher than the reductions achieved when the initial population of mask wearers was 65% (compare (a) with (b)). Furthermore, unlike in the case where the initial population of mask wearers is 65%, the cumulative mortality recorded for all four mask types when the initial size of wearers is 80% are lower than the baseline cumulative mortality (Figure 7 (b)), with respirators much more effective than surgical and cloth masks. When 95% of members of the community chooses to wear a mask initially (i.e., from the very beginning of the third wave of the pandemic), our simulations show a dramatic decrease in cumulative mortality, in comparison to the baseline, regardless of the mask type used (Figure 7 (c)). Although the mask effectiveness is essentially the same, regardless of the mask type chosen, respirators are marginally better than the others, while properly fitted surgical masks are also marginally better than cloth masks (Figure 7 (c)). Figure 7 (f) further shows that if 95% of the US populace opted to wear respirators in public, from the beginning of the third pandemic wave, daily COVID-19 mortality recorded in the country will fall below one by mid June 2021. Same dramatic reduction in daily mortality will be achieved (under this scenario), but a month later, if surgical masks were preferred instead. However, if basic and improved cloth masks are preferred, such reduction can only be achieved around mid January 2022 and December 2021, respectively. In summary, the simulations in Figures 7 (a)-(c) show that respirators are always more effective, in reducing cumulative mortality (in comparison to either the worst-case or baseline scenario), than any of the other mask types considered in this study. Furthermore, these results show that the population-level effectiveness of the four mask types is essentially identical if the initial size of mask wearers is very high (e.g., > 85%). The results show similar pattern with respect to daily mortality (Figures 7 (d)-(f)).
The extended model is further simulated to assess the combined impact of the cumulative initial number of individuals who wear face masks and masking behavior on the spread of the pandemic. For these simulations, we take the reference point to be October 12, 2020 (i.e., the beginning of the third wave of the pandemic). Figure 8 shows that, for the case where 65% of Americans are wearing face masks by October 12, 2020, the projected cumulative mortality by April 30, 2021 for the scenario with no new mask wearers ((i.e., α12 = α13 = 0) and existing mask wearers do not abandon mask wearing (i.e., α21 = α31 = 0) exceeds the projected cumulative mortality for the baseline scenario (Figure 8 (a), gold curve; in comparison to blue curve) by about 4%. If, on the other hand, more Americans choose to wear masks after October 12, 2020 (in addition to the 65% initially) no new wearers abandon mask-wearing (i.e., α12 and α13 are kept at their baseline values, while α21 and α31 are set to zero), the projected cumulative mortality is decreased below the projected cumulative mortality at baseline, by about 26% (Figure 8 (a), green curve; in comparison to blue curve). For the case where the initial percentage of Americans who were wearing face masks by October 12, 2020 is 80%, our simulations show that even the case with no new mask wearers resulted in a decrease in projected baseline cumulative mortality, by about 31% when masks-wearers do not stop wearing masks (Figure 8 (b), magenta curve; in comparison to blue curve). This reduction increases to 43% if more Americans choose to wear masks, and no mask-wearers abandon their mask-wearing behavior (Figure 8 (b), green curve, in comparison to blue curve). Further reductions in projected cumulative mortality are recorded if the percentage of initial masks wearers by October 12, 2020 is 95% (Figure 8 (c)).
Simulations of the extended model (S.1)-(S.3), given in Supplementary Material. Figures (a)-(c): As-sessment of the impact of increases in the total initial size of mask wearers on cumulative mortality. (a) The initial population of mask wearers is 65%. (b) The initial population of mask wearers is 80%. (c) The initial population of mask wearers is 95%. The magenta curve (no behavior change) depicts the scenario in which there is no new masking and already masking individuals do not give up masking, while the green curve depicts the scenario in which new people opt for masking and existing masking individuals do not abandon masking. The magenta, cyan, and yellow regions denote the first pandemic wave, the second pandemic wave, and the portion of the third pandemic wave up to March 11, 2021, respectively, while the dashed vertical blue lines demarcate the fist and second wave and the portion of the third wave up to March 11, 2021. The values of the other parameters of the extended model are as given in Tables 1 and S5-S6 in Supplementary Material. The shaded magenta and green regions represent the beginning and end of the first and second pandemic waves, respectively. The shaded yellow region represents the beginning of the third pandemic wave until March 11, 2021. Dashed vertical blue lines demarcate the three waves.
The simulations in Figure 8 clearly show that there is a threshold level of initial size of population of mask wearers above which the projected baseline cumulative mortality by April 30,2021 will always be reduced even in the absence of new wearers (α12 = α13 = 0), while existing wearers remain wearers (i.e., α21 = α31 = 0, and α23 and α32 maintained at baseline). Larger reductions in the projected baseline cumulative mortality are recorded if more new wearers are attracted (and existing wearers do not opt out of wearing masks), as expected (green curves in Figure 8). Similar results for different levels of other NPI’s, as opposed to different initial sizes of mask-wearers (at the beginning of the third wave of the pandemic), are presented in the Supplementary material (Figure S4).
3 Discussion and Conclusions
With numerous classes of respiratory protection available, including cloth masks, surgical/procedure masks, and respirators, there has been little clarity on the widespread use of respirators among the general public and to what extent their use might curtail a pandemic driven by the respiratory transmission route. In this study, we formulated a basic dynamic epidemiology model for assessing the population-level impact of different types of face masks and other NPIs on curtailing and mitigating the burden of the COVID-19 pandemic in the US. The basic model, which takes the form of a deterministic system of nonlinear differential equations, is calibrated based on updated COVID-19 mortality data from the three waves of the COVID-19 pandemic in the US.
During the early stages of the COVID-19 pandemic in the US, masks were in limited supply. In particular, respirators were in exceptionally limited supply, and prioritized for frontline healthcare workers. We carried out numerical simulations of the basic model to assess the hypothetical scenario where masks are in abundant supply, and are deployed to the general public during the early stages of the pandemic. The simulations show that the early implementation of this masking strategy, particularly if respirators or well-fitted surgical masks are universally-adopted nationwide, would have greatly suppressed the burden (as measured in terms of cumulative mortality) and potential of the COVID-19 pandemic in the US In particular, there would have been no significant outbreak of the pandemic in the US if four out of every five Americans had chosen to mask up early and prioritized respirators or well-fitting surgical masks. Our study shows that the prospect of COVID-19 elimination in the US is greatly enhanced if a sizable proportion of the US population opt to wear high quality masks, notably respirators or tightly-fitted surgical/procedure masks, from the onset of the pandemic (at least 40% if respirators are used; and 55% if surgical/procedure masks are used). The compliance level needed to achieve such elimination is significantly higher (about 75%) if data for the first wave of the pandemic is used to calibrate the model, instead of the data for the third wave (see Figure S1 in Supplementary Material; this is because the baseline value of the control reproduction number of the model, ℝc ≈ 3.1, is much higher during the first wave, in comparison to its value during the third wave). In other words, our study supports the “hit hard, hit early” public health containment strategy, but with high quality masks. Thus, this study suggests that high-efficacy respirator protection should be stockpiled and made readily available to members of the public in times of future major outbreaks of respiratory diseases. This result is consistent with reports in [47–49], which emphasize the superiority of mask-based strategies using high quality masks. The mask compliance needed to eliminate the pandemic using respirators drops to almost one in five Americans if the masking strategy is complemented with other NPIs that can reduce community transmission by 20% from their baseline at the beginning of the third wave of the pandemic. This result is consistent with that reported by Hellewell et al. [50], which suggests that even minor decreases in community transmission of COVID-19 can have a significant impact on efforts aimed at containing the pandemic.
Our study shows that communities that prioritize the use of low efficacy face coverings, such as cloth masks, must complement this strategy with other community transmission-reducing NPIs (such as social-distancing, self-isolation, or even community lockdowns) at very high levels of effectiveness and coverage in order to effectively control the spread of the pandemic. On the other hand, communities that adopt high-efficacy respiratory protection would only need to, in the case of well-fitted surgical masks, complement the masking strategy with other NPIs at low to moderate levels of effectiveness and coverage. No complementing with other NPIs would be needed if respirators or equivalent, at moderate coverage, were used. Prioritising highly-efficient respiratory protection, particularly respirators, if adopted by 95% of the populace at the beginning of the third pandemic wave in the US, is estimated to eliminate the daily COVID-19 mortality by mid-June, 2021.
We showed that if the overall face masks compliance is moderate (e.g., 50% compliance), prioritizing low-efficacy face masks will not be effective in reducing the burden of the pandemic. In fact, the use of low-efficacy masks in this setting (with 50% compliance) increases the disease burden, in comparison to the baseline scenario. Opting, instead, for well-fitting surgical or respirators (in this 50% compliance scenario) will lead to a significant reduction in pandemic burden. If masks use compliance in the community is moderate (e.g., 70%), our study shows that surgical masks and respirators are almost equally as effective in reducing pandemic burden, with respirators being marginally more effective. However, the effectiveness of the face mask strategy, particularly if masks other than respirator or equivalent are prioritized, is greatly enhanced if the mask-based strategy is combined with other NPIs that reduce community transmission from its baseline level (e.g., NPIs that reduce baseline community trans-mission by at least 20%). Our simulations clearly support the supremacy of respirators over all other mask types in reducing and mitigating the burden of the COVID-19 pandemic in the US.
In order to assess the impact of human behavior, with regard to whether or not users choose to wear face masks in public, as well as relaxing mask mandates, we extended the basic model into a three-group model in which the total population is subdivided into the individuals who do not habitually wear face masks in public (Group 1), those who habitually wear cloth or surgical masks (Group 2), and those who habitually wear an respirators or equivalent (Group 3). This multigroup structure allows for the assessment of the impact of the initial sizes of the three subgroups and the adoption, adherence to, or relaxation of mass mandates on the burden of the pandemic.
Numerical simulations of the extended multi-group model demonstrate the importance of positive change of behavior, from non-masking to masking, in significantly reducing the burden (cumulative and daily mortality and hospitalizations) of the pandemic, as expected. This reinforces the urgent need for effective community outreach campaigns to improve mask adoption and opting for higher-efficacy options. This study shows that the effectiveness of one mask type over other types is dependent on the overall mask compliance level in the community. We showed that if the mask compliance level in the community is very high (e.g., 85%), then the four mask types (cloth, improved cloth, surgical, or respirator) are essentially equally as effective. For communities that can maintain high mask-use compliance, the type of mask chosen or adopted in the community may not matter because, for this high compliance setting, any of the four mask types will be very effective in significantly reducing the pandemic burden. This result is consistent with an observational study reported by Wang et al. [51], which indicated that universal household use of face masks was about 79% effective in curtailing the spread of COVID-19. As expected, our simulations also showed that more projected COVID-19 related deaths are averted as more individuals adopt mask wearing and while existing wearers do not opt out of wearing masks. This implies that continued adoption of, and strictly complying with, the face mask use strategy, especially if initial masking coverage in the community is high, is beneficial in effectively reducing and mitigating the burden of the COVID-19 pandemic in the US. This is particularly significant in most nations where vaccine availability continues to be highly-limited.
In summary, our study demonstrates that the prospect for the effective control or elimination of COVID-19 in the United States through modifying the public health masking strategy is promising both singly, if respirators are prioritized, or in combination with other NPIs that reduce baseline community transmission. Our results suggest that the COVID-19 pandemic could have been heavily suppressed or eliminated (i.e., the control reproduction number, ℝc, could have been reduced to a value less than one) if only around half of the US population opted to use respirators from the beginning of the pandemic, even in the absence of needing to complement masking with other NPIs. Such elimination can also be achieved if well-fitted surgical masks were opted for, but with compliance of at least 70% and combined with other NPIs that can reduce baseline community transmission by at least 20%. Our study gives a range of plausible NPI-based scenarios for effectively controlling, and mitigating the burden of, the COVID-19 pandemic in the US that would not necessarily require major economic disruptions, such as imposing community lockdowns or wide scale business and school closures. Some of the limitation of our modeling study is the assumption that the population is well-mixed, such that everyone is equally likely to mix with everyone else. Further, we assume mask efficacy as being constant. However, when not worn properly (i.e., if not tightly-fitted), the actual efficacy of the mask type may not be achieved. Nonetheless, modeling studies can provide insight into the transmission dynamics of the disease and how best to effectively combat its spread and mitigate its burden. As in the case of the prior two deadly coronavirus global outbreaks (SARS and MERS), our study suggests that COVID-19 is controllable using basic public health interventions, particularly using a mask-based strategy that prioritizes high-quality masks, particularly respirators at moderate to high levels of compliance.
Data Availability
All data used in this study is available publicly.
Supplementary Material: Formulation, Data Fitting and Parameter Estimation of of Extended Model
Author contributions
CNN and ABG designed, analysed, parametrized and simulated the dynamical models, in addition to drafting (and editing) the manuscript and supplementary material. JRK contributed in drafting the Introduction section and suggesting edits for the Abstract and Discussion sections of the manuscript. LM contributed in suggesting edits for the Abstract, Introduction and Discussion sections of the manuscript. MHB contributed in drafting the Discussion, in addition to suggesting edits for the rest of the manuscript.
Acknowledgments
One of the authors (ABG) acknowledge the support, in part, of the Simons Foundation (Award #585022) and the National Science Foundation (grant number: DMS-2052363). CNN acknowledges the support of the Simons Foundation (Award #627346).