Hostname: page-component-848d4c4894-v5vhk Total loading time: 0 Render date: 2024-06-22T10:01:27.622Z Has data issue: false hasContentIssue false

Social divisions and risk perception drive divergent epidemics and large later waves

Published online by Cambridge University Press:  23 February 2023

Mallory J. Harris*
Affiliation:
Biology Department, Stanford University, Stanford, CA 94301, USA
Kimberly J. Cardenas
Affiliation:
Biology Department, Stanford University, Stanford, CA 94301, USA
Erin A. Mordecai
Affiliation:
Biology Department, Stanford University, Stanford, CA 94301, USA
*
*Corresponding author. E-mail: mharris9@stanford.edu

Abstract

During infectious disease outbreaks, individuals may adopt protective measures like vaccination and physical distancing in response to awareness of disease burden. Prior work showed how feedbacks between epidemic intensity and awareness-based behaviour shape disease dynamics. These models often overlook social divisions, where population subgroups may be disproportionately impacted by a disease and more responsive to the effects of disease within their group. We develop a compartmental model of disease transmission and awareness-based protective behaviour in a population split into two groups to explore the impacts of awareness separation (relatively greater in- vs. out-group awareness of epidemic severity) and mixing separation (relatively greater in- vs. out-group contact rates). Using simulations, we show that groups that are more separated in awareness have smaller differences in mortality. Fatigue (i.e. abandonment of protective measures over time) can drive additional infection waves that can even exceed the size of the initial wave, particularly if uniform awareness drives early protection in one group, leaving that group largely susceptible to future infection. Counterintuitively, vaccine or infection-acquired immunity that is more protective against transmission and mortality may indirectly lead to more infections by reducing perceived risk of infection and therefore vaccine uptake. Awareness-based protective behaviour, including awareness separation, can fundamentally alter disease dynamics.

Social media summary: Depending on group division, behaviour based on perceived risk can change epidemic dynamics & produce large later waves.

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
Copyright © The Author(s), 2023. Published by Cambridge University Press

Introduction

When an infectious disease causes substantial disease burden and death, people may perceive their risk of infection based on their awareness of the magnitude of disease-linked outcomes and respond by modifying their behaviour (An et al., Reference An, Hawley, Van Horn, Bacon, Yang and Resnicow2020; Cheok et al., Reference Cheok, Gatot, Sim, Ng, Tay, Howe and Koh2021; Gidengil et al., Reference Gidengil, Parker and Zikmund-Fisher2012; Ridenhour et al., Reference Ridenhour, Sarathchandra, Seamon, Brown, Leung, Johnson-Leon and Johnson-Leung2022; Yan et al., Reference Yan, Malik, Bayham, Fenichel, Couzens and Omer2021). In turn, protective behaviours like physical distancing, mask wearing and vaccination may suppress transmission, reducing peak and total infections and disease-linked mortality (Abaluck et al., Reference Abaluck, Kwong, Styczynski, Haque, Kabir, Bates-Jefferys and Mobarak2022; Toor et al., Reference Toor, Echeverria-Londono, Li, Abbas, Carter, Clapham and Gaythorpe2021; Yan et al., Reference Yan, Malik, Bayham, Fenichel, Couzens and Omer2021). Awareness-based behaviour describes protective measures that are adopted in response to epidemic intensity. Bidirectional feedback between protective behaviour and epidemic intensity can lead to unexpected and nonlinear dynamics, such as plateaus and oscillations in cases over time, if protective measures are abandoned over time (e.g. fatigue with non-pharmaceutical interventions may lead to a regular decline in adherence) or the strength of protection wanes (e.g. waning immunity from vaccination or infection) (Arthur et al., Reference Arthur, Jones, Bonds, Ram and Feldman2021; Eksin et al., Reference Eksin, Shamma and Weitz2017; Perra et al., Reference Perra, Balcan, Gonçalves and Vespignani2011; Weitz et al., Reference Weitz, Park, Eksin and Dushoff2020). Models that split the population into categories with respect to the disease (i.e. compartments) and mathematically define transition rates between different states are widely used to understand such complex epidemic dynamics. Compartmental models may incorporate awareness as a function of deaths or cases that reduces transmission evenly across the population (Arthur et al., Reference Arthur, Jones, Bonds, Ram and Feldman2021; Weitz et al., Reference Weitz, Park, Eksin and Dushoff2020). However, real populations are sharply divided in physical interactions, demography, ideology, education, housing and employment structures, and information access. These social divisions can impact the transmission of both pathogens and information within and between groups, altering epidemic dynamics. The impacts of such asymmetrically spreading disease and awareness in a highly divided population are not well understood (Acevedo-Garcia, Reference Acevedo-Garcia2000; Farmer, Reference Farmer1996; Grief & Miller, Reference Grief and Miller2017).

Populations may be subdivided based on an array of factors (e.g. race, ethnicity, age and geography), with marked differences in pathogen exposure and infection severity (Farmer, Reference Farmer1996; Greene et al., Reference Greene, Levin-Rector, Hadler and Fine2015; Li et al., Reference Li, Wang, Gao, Xu and Chen2016; Poteat et al., Reference Poteat, Millett, Nelson and Beyrer2020; Williams & Cooper, Reference Williams and Cooper2020; Zelner et al., Reference Zelner, Trangucci, Naraharisetti, Cao, Malosh, Broen and Delamater2020). Risk of pathogen introduction may vary between groups: high-income groups may encounter pathogens endemic to other regions through international travel, low-income groups may have heightened likelihood of exposure connected to poor housing quality and insufficient occupational protections, and certain regions and occupations experience greater risks of exposure to zoonotic illnesses (Benfer et al., Reference Benfer, Vlahov, Long, Walker-Wells, Pottenger, Gonsalves and Keene2021; Cubrich, Reference Cubrich2020; Dhewantara et al., Reference Dhewantara, Mamun, Zhang, Yin, Ding, Guo and Soares Magalhães2018; Greene et al., Reference Greene, Levin-Rector, Hadler and Fine2015; Pramasivan et al., Reference Pramasivan, Ngui, Jeyaprakasam, Liew, Low, Mohamed Hassan and Vythilingam2021). Once a pathogen is introduced, it may spread at different rates within groups based on factors like housing density and access to healthcare (Benfer et al., Reference Benfer, Vlahov, Long, Walker-Wells, Pottenger, Gonsalves and Keene2021; Poteat et al., Reference Poteat, Millett, Nelson and Beyrer2020; Quinn et al., Reference Quinn, Kumar, Freimuth, Musa, Casteneda-Angarita and Kidwell2011). Further, the severity of infection may vary directly with group identity owing to underlying biological differences (e.g. age or sex), as a function of co-morbidities especially prevalent in one group owing to underlying inequities (e.g. lung disease connected to environmental pollution or heart disease associated with factors driven by structural racism), or through heterogeneity in access to and quality of healthcare (Calvin et al., Reference Calvin, Winters, Wyatt, Williams, Henderson and Walker2003; Lane et al., Reference Lane, Morello-Frosch, Marshall and Apte2022; Li et al., Reference Li, Wang, Gao, Xu and Chen2016; Poteat et al., Reference Poteat, Millett, Nelson and Beyrer2020; Quinn et al., Reference Quinn, Kumar, Freimuth, Musa, Casteneda-Angarita and Kidwell2011; Takahashi et al., Reference Takahashi, Ellingson, Wong, Israelow, Lucas, Klein and Iwasaki2020; Williams & Cooper, Reference Williams and Cooper2020; Wu et al., Reference Wu, Nethery, Sabath, Braun and Dominici2020). Mixing, or between-group contact rates, can alter transmission dynamics. Physical barriers (e.g. geographic boundaries, schools, residential segregation and incarceration) and preferential contact with members of one's own group may reduce interactions and subsequent transmission between groups, a characteristic we describe as separated mixing (Arnold et al., Reference Arnold, Srinivasan, Rodriguez, Rydzak, Herzog, Gontu and Ferrari2022; Doherty et al., Reference Doherty, Schoenbach and Adimora2009; Greene et al., Reference Greene, Levin-Rector, Hadler and Fine2015; Harris et al., Reference Harris, Tessier-Lavigne and Mordecai2021; Rothenberg et al., Reference Rothenberg, Muth, Malone, Potterat and Woodhouse2005). Infectious disease models that account for differences in vulnerability within subgroups of a population and separated mixing can help to illustrate the emergence of health inequities and justify structural interventions to reduce these disparities (Jacquez et al., Reference Jacquez, Simon, Koopman, Sattenspiel and Perry1988; K. C. Ma et al., Reference Ma, Pei, Shaman, Dubrow and Chen2021; Richardson et al., Reference Richardson, Malik, Darity, Mullen, Morse, Malik and Jones2021; Zelner et al., Reference Zelner, Masters, Naraharisetti, Mojola, Chowkwanyun and Malosh2022). However, such models may miss an important behavioural dimension by failing to account for variation in awareness-based behaviour changes among groups.

Awareness and behavioural heterogeneity can significantly alter disease dynamics: for example, protective behaviour adoption based on disease status of social connections may slow pathogen transmission, while social clustering in vaccine exemptions may lead to outbreaks (Funk et al., Reference Funk, Gilad, Watkins and Jansen2009; Herrera-Diestra & Meyers, Reference Herrera-Diestra and Meyers2019; Omer et al., Reference Omer, Enger, Moulton, Halsey, Stokley and Salmon2008). Personal perception of disease severity may be influenced by population-level social norms and mass media, regardless of group identity. However, attitudes toward diseases and protective behaviours may also vary considerably between groups and correspond to actual risk and personal experiences of close social ties with the disease (Anthonj et al., Reference Anthonj, Diekkrüger, Borgemeister and Kistemann2019; Brug et al., Reference Brug, Aro, Oenema, de Zwart, Richardus and Bishop2004; Christensen et al., Reference Christensen, Pilling, Eyring, Dickerson, Sloan and Magnusson2020; Holtz et al., Reference Holtz, Zhao, Benzell, Cao, Rahimian and Yang2020; Oraby et al., Reference Oraby, Thampi and Bauch2014; Simione & Gnagnarella, Reference Simione and Gnagnarella2020). While prior awareness-based models have examined outcomes given different scales of information (i.e. local or global), we aim to characterize risk perception based on group-level information in a population split into two distinct and well-defined groups (Funk et al., Reference Funk, Salathé and Jansen2010). We define separated awareness as greater in- vs. out-group awareness of current epidemic conditions in a split population. We predict that, by producing behavioural responses more reflective of each group's risk, separated awareness may reduce differences between groups in disease burden that might otherwise occur (Steinegger et al., Reference Steinegger, Arola-Fernández, Granell, Gómez-Gardeñes and Arenas2022). Understanding the impacts of separation with respect to mixing and awareness on disease dynamics may be important for characterizing differences in epidemic burden and effectively intervening to mitigate population inequities (K. C. Ma et al., Reference Ma, Menkir, Kissler, Grad and Lipsitch2021; Richardson et al., Reference Richardson, Malik, Darity, Mullen, Morse, Malik and Jones2021; Steinegger et al., Reference Steinegger, Arola-Fernández, Granell, Gómez-Gardeñes and Arenas2022; Weston et al., Reference Weston, Hauck and Amlôt2018; Zelner et al., Reference Zelner, Masters, Naraharisetti, Mojola, Chowkwanyun and Malosh2022).

Here, we investigate the impacts of intergroup divisions on epidemic dynamics using an awareness-based model for transmission of an infectious disease, in which adoption of protective measures (either non-pharmaceutical interventions or vaccinations) is linked to recent epidemic conditions and mediated by awareness.

We ask:

  1. 1. How do separated awareness and mixing interact to affect differences between groups in epidemic dynamics?

  2. 2. How does fatigue interact with awareness separation to affect long-term epidemic dynamics?

  3. 3. When vaccines are introduced, how does immunity interact with awareness separation to affect long-term epidemic dynamics?

Methods

Non-pharmaceutical intervention model

We model disease transmission with awareness-based adoption of non-pharmaceutical interventions that reduce transmission rates. See Supplementary Figure S1 for a compartmental diagram for this model and Supplementary Table S1 for parameter definitions. We model disease transmission with a Susceptible–Infectious–Recovered–Deceased model, tracking the proportion of the population in each compartment through time. Susceptible individuals have never been infected or vaccinated. New infections arise through contact between susceptible and infected individuals, with the transmission coefficient β describing the rate at which the pathogen spreads. Individuals exit the infected compartment at per capita rate ρ, the inverse of infectious period 1/ρ and either recover or die. The fatality probability, or fraction of individuals exiting the infectious compartment who die, is μ (meaning that recovery after infection occurs with probability 1 − μ). In this model, recovered individuals have durable immunity and cannot be reinfected. The initial model does not include vaccine-derived immunity, an extension we consider below (Equation 3).

We further categorize the population based on whether they adopt behaviour that is Protective (P) or Unprotective (U). Compartment names contain two letters, the first indicating disease status and the second indicating behaviour (e.g. SU denotes Susceptible people with Unprotective behaviours). We track the behavioural status of Recovered and Deceased individuals (at the time of death), although they do not contribute directly to transmission. Protective measure efficacy against infection is determined by a scaling factor κ describing the degree to which the behaviour prevents infection (where κ = 0 corresponds to complete protection and κ = 1 corresponds to no protection). Protective measures affect the behaviour of both susceptible and infected individuals, so transmission rate is reduced by a factor of κ 2 in encounters where both parties have adopted protective measures. Living individuals can switch between protective and unprotective behaviour, and we assume that the rates of these behavioural transitions are independent of their own disease status. Unprotective individuals adopt protective behaviours based on awareness (α(t), Equation 2), or perceived epidemic intensity at a given point in time. Awareness is the product of disease-induced deaths over the past ℓ days (making ℓ a measure of memory) and a responsiveness constant θ. Protective behaviours are abandoned owing to fatigue at per capita rate ϕ.

To study the impact of social divisions, we further split the population into two groups of equal size, where group membership is fixed, and each group contains all epidemiological and behavioural compartments. The groups are labelled as a and b and indicated as a subscript in compartment names (e.g. SU a corresponds to the prevalence of Susceptible–Unprotective individuals in group a). We arbitrarily designate group a as having greater underlying vulnerability to infection or disease-linked mortality in all of the following scenarios. Parameters may vary between groups, as indicated by subscripts (e.g. θ a corresponds to responsiveness in group a). If parameters are equivalent for both groups, we exclude the subscript (e.g. θ = θ a = θ b).

Preferential within-group mixing is represented by homophily parameter h, corresponding to the proportion of contacts that are within-group. When h is 0.5, mixing is uniform, meaning that individuals are equally likely to contact members of their own group as members of the opposite group. As h approaches 1, mixing becomes increasingly separated, meaning that contacts are increasingly concentrated within groups. Similarly, we consider separation in awareness, $\epsilon$, or the relative weight of in-group vs. out-group awareness of deaths for protective behaviour.

The system of equations for group a is as follows (equations for group b can be derived symmetrically):

(1)$$\eqalign{& \mathop {SU_a}\limits^{.} = {-}\beta SU_a( {( h ) ( {IU_a + \kappa IP_a} ) + ( {1-h} ) ( {IU_b + \kappa IP_b} ) } ) -\theta SU_a\alpha _a( t ) + \phi SP_a \cr & \mathop {SP_a}\limits^{.} = {-}\beta \kappa SP_a( {( h ) ( {IU_a + \kappa IP_a} ) + ( {1-h} ) ( {IU_b + \kappa IP_b} ) } ) + \theta SU_a\alpha _a( t ) -\phi SP_a \cr & \mathop {IU_a}\limits^{.} = \beta SU_a( {( h ) ( {IU_a + \kappa IP_a} ) + ( {1-h} ) ( {IU_b + \kappa IP_b} ) } ) -\theta IU_a\alpha _a( t ) + ( {\phi -\rho } ) IP_a \cr & \mathop {IP_a}\limits^{.} = \beta \kappa SP_a( {( h ) ( {IU_a + \kappa IP_a} ) + ( {1-h} ) ( {IU_b + \kappa IP_b} ) } ) + \theta IU_a\alpha _a( t ) -( {\phi + \rho } ) IP_a\,\cr & \mathop {RU_a}\limits^{.} = ( {1-\mu } ) \rho IU_a-\theta RU_a\,\alpha _a( t ) + \phi RP_a \cr & \mathop {RP_a}\limits^{.} = ( {1-\mu } ) \rho IP_a + \theta RU_a\alpha _a( t ) -\phi RP_a\,\cr & \mathop {DU_a}\limits^{.} = \mu \rho IU_a\,\cr & \mathop {DP_a}\limits^{.} = \mu \rho IP_a\;} $$

where α a(t) is the awareness equation for group a:

(2)$$\alpha _a( t ) = \mathop \smallint \limits_{t-\ell }^t \left({( {\epsilon_a} ) \left({\mathop {DU_a}\limits^{.} + \mathop {DP_a}\limits^{.} } \right) + ( {1-\epsilon_a} ) \left({\mathop {DU_b}\limits^{.} + \mathop {DP_b}\limits^{.} } \right)} \right)\,{\rm d}t$$

Vaccination model

We develop an alternative model in which the awareness-based behaviour is vaccine uptake, rather than non-pharmaceutical interventions. See Supplementary Figure S2 for a compartmental diagram for this model and Supplementary Table S1 for parameter definitions. Here, the second letter of compartment names indicates immune status: Unprotective (U), Transmission and Mortality-reducing Immunity (T), or Mortality-reducing Immunity (M). This reflects our assumption that immunity initially reduces both transmission and mortality following infection or vaccination, and later wanes to reduce mortality but not infection.

As in the non-pharmaceutical intervention model, susceptible people without prior immunity (SU) may become infected and then recover or die according to baseline infection parameter values. Susceptible individuals may become vaccinated and transition directly to the recovered compartment, bypassing infection, at a rate dependent on awareness. There may be a lag between the beginning of the epidemic and vaccine introduction at time point t v (Supplementary Figures S13 and S14). To evaluate long-term immune effects of vaccination and infection on epidemic dynamics, we incorporate waning immunity by including distinct T and M compartments, as described above.

After vaccination or infection, individuals temporarily have complete protection from infection (RT). At per capita rate ω, they regain susceptibility to infection, this time with transmission and mortality-reducing immunity (i.e. ST). As in the non-pharmaceutical intervention model, transmission-reducing protection scales transmission rates for susceptible and infected individuals by a constant. Additionally, immunity reduces disease-linked mortality by scaling factor ζ. Transmission-reducing immunity is lost at per capita rate ϕ, while mortality-reducing immunity is retained over the course of the simulation, reflecting how neutralizing antibody production may decay over time while cellular immune responses are more durable (Siggins et al., Reference Siggins, Thwaites and Openshaw2021). Susceptible individuals with mortality-reducing immunity alone (SM) may regain transmission-reducing immunity via vaccination, which occurs based on the same awareness function as vaccination of people without immune protection.

The system of equations for this model in a population without groups is:

(3)$$\eqalign{& \mathop {SU}\limits^{.} = {-}{\rm \beta }SU( {IU + {\rm \kappa }IT + IM} ) -{\rm \theta }SU\mathop \smallint \limits_{t-\ell }^t \left({\mathop {DU}\limits^{.} + \mathop {DT}\limits^{.} + \mathop {DM}\limits^{.} } \right){\rm d}t\,\cr & \mathop {ST}\limits^{.} = {\rm \omega }RT-{\rm \beta \kappa }ST( {IU + {\rm \kappa }IT + IM} ) -\phi ST \cr & \mathop {SM}\limits^{.} = {-}{\rm \beta }SM( {IU + {\rm \kappa }IT + IM} ) -{\rm \theta }SM\mathop \smallint \limits_{t-\ell }^t \left({\mathop {DU}\limits^{.} + \mathop {DT}\limits^{.} + \mathop {DM}\limits^{.} } \right){\rm d}t + \phi ST\,\cr & \mathop {IU}\limits^{.} = {\rm \beta }SU( {IU + {\rm \kappa }IT + IM} ) -{\rm \rho }IU\,\cr & \mathop {IT}\limits^{.} = {\rm \beta \kappa }ST( {IU + {\rm \kappa }IT + IM} ) -{\rm \rho }IT \cr & \mathop {IM}\limits^{.} = {\rm \beta }SM( {IU + {\rm \kappa }IT + IM} ) -{\rm \rho }IM\,\cr & \mathop {RT}\limits^{.} = {\rm \rho }( {( {1-{\rm \mu }} ) IU + ( {1-{\rm \zeta \mu }} ) ( {IT + IM} ) } ) -{\rm \omega }RT + {\rm \theta }( {SU + SM} ) \mathop \smallint \limits_{t-\ell }^t \left({\mathop {DU}\limits^{.} + \mathop {DT}\limits^{.} + \mathop {DM}\limits^{.} } \right){\rm d}t\,\cr & \mathop {DU}\limits^{.} = ( {{\rm \mu \rho }} ) IU \cr & \mathop {DT}\limits^{.} = ( {{\rm \zeta \mu \rho }} ) IP\,\,\cr & \mathop {DM}\limits^{.} = ( {{\rm \zeta \mu \rho }} ) IM} $$

The equations for a split population with separated mixing and awareness can be derived following Equation (1).

Simulations

We ran simulations in R version 4.0.2, using the dede function in the deSolve package, which solves systems of differential equations (Soetaert et al., Reference Soetaert, Petzoldt and Setzer2010). The population begins as almost fully susceptible (S(0) ≈ 1), with a small initial infection prevalence (I(0)) to seed the outbreak and no protective behaviours. In the non-pharmaceutical intervention scenarios (scenario 1 and 2), the sole initial difference between groups is caused by introducing the pathogen into group a alone at prevalence I a(0) = 0.001. In the vaccination scenario (scenario 3), the pathogen is introduced in both groups at prevalence I(0) = 0.0005 and the fatality probability for group a is twice that of group b (μ a = 0.02 and μ b = 0.01). An interactive R Shiny app that allows users to simulate epidemics for the non-pharmaceutical intervention model across parameter values is available at https://mallory-harris.shinyapps.io/divided-disease/.

Results

1. Separated mixing and awareness

To understand how separation in awareness and mixing interact to alter short-term epidemic dynamics in a split population, we model awareness-based adoption of non-pharmaceutical interventions (Equation 1); all model parameters are defined in Supplementary Table S1 and a compartmental diagram is provided as Supplementary Figure S1. As described above, the pathogen is introduced in group a alone; all other parameters are equivalent between groups. To simplify short-term awareness-based behaviour, this scenario does not incorporate memory or fatigue (ℓ = 1 and ϕ = 0). First, we allow both mixing (h, which drives the contact and contagion process) and awareness ($\epsilon$, which drives protective behaviour adoption) to be either uniform (functioning like a single population) (0.5) or highly separated (0.99).

The groups experience identical epidemic dynamics regardless of awareness separation when mixing is uniform (Figure 1A, B), as the pathogen introduced into group a quickly spreads into group b and circulates evenly within and between groups. When groups mix separately, differences in epidemic dynamics between groups arise and depend on awareness separation (Figure 1C, D). Therefore, we focus the rest of our analyses on cases where mixing is separated to examine the impacts of awareness separation. When awareness is uniform, epidemic shape differs in both timing and magnitude between groups, increasing the peak size and total infections in the more vulnerable (earlier epidemic introduction) group a and decreasing both in group b (Figure 1C). Group a also has more cumulative deaths than group b under uniform awareness, while cumulative deaths across the full population (group a and group b combined) are approximately constant across different levels of awareness and mixing separation (Supplementary Figure S3).

Figure 1. Epidemic peaks are offset in time between groups when mixing is separated (C, D), and in magnitude when awareness is uniform but mixing is separated (C). Plots show the prevalence of infections over time in group a (pink) and group b (green) under four scenarios: awareness is uniform (A, C; $\epsilon = 0.5$) or separated (B, D; $\epsilon = 0.99$); mixing is uniform (A, B; h = 0.5) or separated (C, D; h = 0.99). We assume the pathogen is introduced only in group a at prevalence 0.001 and that all other parameters are equivalent between groups: transmission coefficient (β = 0.2), infectious period (1/ρ = 10), fatality probability (μ = 0.01), protective measure efficacy (κ = 0.3), responsiveness (θ = 100), memory (ℓ = 1), and fatigue (ϕ = 0). Lines overlap under uniform mixing (top row).

Awareness separation changes epidemic size in both groups by modulating how quickly protective behaviour arises relative to pathogen spread (Figure 2). Uniform awareness reduces total infections in group b, which adopts protective behaviour by observing mortality in group a at a point when infections within group b remain relatively low (Figures 1C and 2B, D, E). Meanwhile, uniform awareness causes group a to underestimate disease severity owing to the lack of early mortality in group b, leading to decreased early protective behaviour and a larger outbreak (Figures 1C and 2A, C, E). When awareness is separated, group b has little awareness of the emerging epidemic localized to group a, while group a responds to its relatively higher early disease burden with increased awareness, driving epidemic dynamics between the two groups to be similar in shape but delayed in time for group b (Figure 1D). Therefore, awareness separation reduces the differences between groups in epidemic shape (e.g. peak size, total infections), while mixing separation offsets them in time (Figure 1C, D, Supplementary Figures S4 and S5).

Figure 2. Separated awareness reduces between-group differences by reducing group b's awareness of the emerging epidemic and augmenting group a's response to the introduction of the pathogen. We initialize our model using the same parameters as Figure 1 with separated mixing (h = 0.99). We compare uniform awareness ($\epsilon = 0.5$; dashed lines) and separated awareness ($\epsilon = 0.99$; solid lines). At the top, we compare early time series (through t = 80) of (A) protective attitude prevalence in group a; (B) protective attitude prevalence in group b; (C) cumulative infections in group a; (D) cumulative infections in group b. (E) A phase portrait of protective attitude prevalence against cumulative infections in group a (pink) and group b (green). Points indicate values at t = 80, corresponding to the end of the time series in (A–D). Arrows indicate differences in protective attitude prevalence (grey) and cumulative infections (black) at t = 80 for separated vs. uniform awareness, with letters corresponding to time series panel labels.

Differences between groups in epidemic dynamics only arise at high levels of mixing separation (h > 0.9) but can occur at intermediate levels of awareness separation (Supplementary Figures S4 and S5; e.g. $\epsilon = 0.75$). Awareness separation also reduces differences between groups in severe outcomes when groups differ in their transmission coefficients, infectious periods, or fatality probabilities (Supplementary Figures S6–8).

2. Fatigue and awareness separation

We introduce memory and fatigue to examine the long-term impacts of separated awareness when awareness-driven protective behaviour is abandoned over time. Once again, the pathogen is introduced into group a alone and all other parameters are equivalent between groups. To maintain between-group differences, we assume separated mixing (h = 0.99).

In all cases, when protective behaviour wanes with fatigue, three distinct peaks emerge before transmission plateaus at low levels and declines gradually (Figure 3). The initial difference between groups with uniform awareness means that group b retains a relatively larger proportion of susceptible individuals who avoided infection in the first wave by rapidly adopting protective behaviours (Figures 1C and 3A). As a result, the second and third waves in group b exceed the first wave in peak and total infections (Figure 3A). Meanwhile, uniform awareness causes the second and third waves in group a to be smaller compared with separated awareness (Figure 3A vs. B). Under uniform awareness, the third wave in group a is considerably delayed, peaking around 800 days (vs. 450 days under separated awareness). At intermediate awareness separation ($\epsilon = 0.75)$, the first and second waves in group b are approximately equivalent in size (Supplementary Figure S10). As shown in the case without memory and fatigue (Figure 1), when both mixing and awareness are separated, the groups differ mainly in the timing of epidemic peaks rather than in their magnitude, before converging on a long and slow decline (i.e. shoulder; Figure 3B). In the full population, awareness separation may change infection prevalence over time but has no impact on cumulative deaths (Supplementary Figure S11).

Figure 3. Fatigue and long-term memory produce multiple epidemic peaks, which exceed the size of the initial peak in group b when uniform awareness and separated mixing leave that group with a high proportion of susceptible people following the first wave. We initialize the model with separated mixing (h = 0.99), long-term memory (ℓ = 30), and fatigue (ϕ = 0.02); all other parameters are the same as in Figure 1. We consider infections in group a (pink) and group b (green) over a longer time period (1000 days, compared with 200 days in Figure 1). The panels correspond to (A) uniform awareness ($\epsilon = 0.5$) and (B) separated awareness ($\epsilon = 0.99$).

3. Immunity and awareness separation

Next, we consider the implications of awareness-based vaccine uptake in a split population given waning immune protection against infection and durable protection against mortality (Equation 3, Supplementary Figure S2). We model immunity from prior infection as equivalent to immunity from vaccination. Unlike in the previous analyses, the pathogen is now introduced at the same prevalence in both populations simultaneously to ensure that groups a and b begin the post-vaccine period with similar levels of immunity. Group differences are driven by a fatality probability in group a that is twice that of group b. Again, we assume separated mixing (h = 0.99) to maintain distinct dynamics between the groups. We initiate vaccination at 200 days, after an initial large wave of infections. Our analyses focus on the period following the introduction of vaccines to understand how awareness separation modulates the impact of this protective measure across a period where infection is already well established in both populations but substantial proportions of the population remain susceptible.

After an initial large wave (displayed in Supplementary Figure S12), vaccination and waning immunity lead to damped cycles of infections and deaths (Figure 4). As was the case with the non-pharmaceutical intervention model (Figure 1), when awareness drives vaccination behaviour, separated awareness helps to reduce differences in mortality between groups (Figure 4D vs. C). Group a becomes vaccinated at a higher rate in response to the greater number of deaths observed in group a, an effect that is most notable during the second epidemic peak following vaccine introduction (Figure 4D). Therefore, group a also has fewer infections than group b in later waves under separated awareness (Figure 4B), while the two groups experience identical infection dynamics (despite the larger disparity in deaths) given uniform awareness (Figure 4A).

Figure 4. Waning immunity and awareness-based vaccination drive epidemic cycles; separated awareness reduces the disparity in deaths (C vs. D) as more vulnerable group a members become vaccinated at a higher rate. We consider infections (A, B) and deaths (C, D) in the post-vaccine period in group a (pink) and group b (green) where the fatality probability for group a is double that of group b (μ a = 0.02 and μ b = 0.01). The x-axis gives time since vaccination began (tv = 200). We compare uniform awareness ($\epsilon = 0.5$) (A, C) and separated awareness ($\epsilon = 0.99$) (B, D). Other parameter values are: β = 0.2 (transmission coefficient), κ = 0.05 (transmission-reducing immunity), ζ = 0.05 (mortality-reducing immunity), ω = ϕ = 0.01 (waning immunity), infectious period (1/ρ = 10), θ = 20 (responsiveness), ℓ = 30 (memory), h = 0.99 (separated mixing), I 0 = 0.0005 (initial infection prevalence). See Supplementary Figure S12 for a time series plot including the pre-vaccine period.

Because vaccination protects against infections and deaths, and recent deaths feed back to influence awareness-driven vaccine uptake, there is a potential tradeoff between immune protection from vaccines and epidemic dynamics. We explored this tradeoff by examining the effect of variation in immune protection on epidemic dynamics and their feedbacks on vaccine uptake rate, assuming that immune protection causes the same proportional reduction in transmission and mortality (κ = ζ). As expected, greater immune protection reduces the number of deaths by directly reducing the fatality probability. However, because of awareness-driven vaccine uptake, vaccination can produce diminishing returns at the population scale where doubling immune protection from death and infection only reduces total deaths by about one eighth owing to the compensatory reduction in vaccine uptake (Figure 5A), despite doubling individual protection for vaccinated people. Since a more effective immune response reduces mortality, the perceived risk associated with infection declines and fewer people become vaccinated (Figure 5B). The tradeoff between the direct impacts of immune protection on preventing infections and reduced uptake produces a nonlinear relationship between total infections and immune protection (Figure 5C). At low immune protection, infections remain approximately constant as immune protection improves. At higher levels of immune protection, reduced uptake with stronger protection leads to more infections (Figure 5C).

Figure 5. Greater immune protection (from vaccination and infection) leads to lower death rates (A), which in turn decreases vaccination rates (B) and increases infection rates (C); separated awareness reduces disparities in death rates (A) as groups are vaccinated at different rates proportional to their risks of death (B), creating differences in infection rates (C). We vary immune protection, defined as transmission-reducing immunity and mortality-reducing immunity, where both parameters are assigned the same values (κ = ζ). We assume immune protection is equivalent for vaccine- and infection-derived immunity. The x-axis is reversed because smaller values indicate stronger protection. We examine the impacts of stronger immune protection (lower values of κ and ζ) on total deaths (A), vaccinations (B) and infections (C) in the post-vaccine period (t = 200 through t = 2200). We consider the post-vaccine period to focus on the impacts of an awareness-based intervention administered under different levels of awareness separation. We compute each quantity for group a (pink) and group b (green) given uniform (dashed lines; $\epsilon = 0.5$) or separated (solid lines; $\epsilon = 0.99$) awareness. Other parameter values are the same as Figure 4.

Separated awareness drives differences between groups in vaccination behaviour – the higher-risk group a gets vaccinated at a higher rate in response to awareness of the greater cumulative mortality in that group (Figure 5B). This in turn increases differences in infections (group a experiences lower infection rates; Figure 5C) but decreases differences in mortality between groups (death rates are lower for group a but higher for group b than in the uniform awareness scenario; Figure 5A). Since group a is at a higher inherent risk of mortality given infection, separated awareness differentially promotes vaccination and reduces infection in this group, while uniform awareness causes group a to ignore its higher risk of mortality (Figure 5A, B, solid vs. dashed lines). Cumulative deaths increase especially quickly during the initial wave absent vaccination because the population lacks transmission- or mortality-reducing immunity. When vaccination begins earlier in the epidemic (prior to the initial peak around t = 100), separated awareness has greater potential to reduce the difference in cumulative deaths between the two groups (Supplementary Figures S13 and S14). Early vaccination may also reduce cumulative deaths and infections in each group (Supplementary Figure 14).

Discussion

Awareness separation and social divisions may interact to fundamentally alter disease dynamics, creating or erasing differences between groups in the timing and magnitude of epidemic peaks. Uniform awareness can exacerbate differences between population subgroups when the more vulnerable group (e.g. the group where the pathogen is introduced or the group with higher fatality probabilities) underestimates the in-group risk of disease and fails to adopt early protective measures (Figures 1 and 5). At the same time, the initially less-vulnerable group receives indirect protection from observing and responding to epidemic effects in the more vulnerable group, adopting protective measures that reduce their total and peak infections (Figures 1 and 5). However, when awareness-driven behaviour fades with fatigue, the relative disease burden may shift between groups such that the group that initially had fewer infections has relatively more infections in subsequent waves, especially when uniform awareness protects the initially less-vulnerable group during the first wave of infection (Figure 2). Awareness separation diminishes between-group differences in severe outcomes (Figures 1–5, Supplementary Figures S6–S8), but may do so by increasing differences in behaviour and infections (Figures 4 and 5, Supplementary Figure S8). For example, when the more vulnerable group has a higher rate of disease-linked mortality, awareness separation leads them to have higher vaccine uptake in response to their heightened perceived (and actual) risk, narrowing the difference in mortality (Figure 5). More broadly, awareness separation generally reduces differences in severe outcomes between groups by producing preferential uptake of preventative measures by the group with the greatest recent mortality, which is usually the group at greatest current risk.

In this model, greater awareness separation generally reduces differences in severe outcomes between groups, but the magnitude of these impacts may vary depending on disease properties (e.g. transmission coefficient) and behavioural and social processes (e.g. responsiveness to disease-linked mortality) (Supplementary Figures S9, S13 and S14). Outcomes may be further modulated by public health orders and the timing of different interventions. For example, there is greater potential for awareness separation to reduce between-group differences in mortality given earlier vaccine introduction (Supplementary Figures S13 and S14). The existing models could be modified to incorporate population-wide measures, particularly time-limited non-pharmaceutical intervention mandates, to study how social and behavioural processes may shift the optimal timing of interventions in the full population or either group (Ketcheson, Reference Ketcheson2021; Morris et al., Reference Morris, Rossine, Plotkin and Levin2021). Although this model and others assume that protective behaviour uptake is independent of disease status (Mehta & Rosenberg, Reference Mehta and Rosenberg2020; Smaldino & Jones, Reference Smaldino and Jones2021), the model could be modified to link behaviour with known disease status (e.g. accelerated uptake of or reduced fatigue with protective measures by people with symptomatic infections) (Eksin et al., Reference Eksin, Shamma and Weitz2017; Funk et al., Reference Funk, Gilad, Watkins and Jansen2009). To assess the robustness of our conclusions about the effects of awareness separation, the same scenarios could be evaluated across different models of awareness-based behaviour changes, including saturation at a certain threshold for deaths (Weitz et al., Reference Weitz, Park, Eksin and Dushoff2020), consideration of both lethal and non-lethal impacts of disease (e.g. hospitalizations and cases), or optimization to balance the benefits of protection against the costs of various measures (Arthur et al., Reference Arthur, Jones, Bonds, Ram and Feldman2021; Barrett et al., Reference Barrett, Bisset, Leidig, Marathe and Marathe2011; Eksin et al., Reference Eksin, Shamma and Weitz2017). The latter approach may clarify a point that is not addressed in our analysis: although awareness separation may reduce disparities in severe disease-linked outcomes, this phenomenon is not necessarily equitable or desirable. In fact, if self-protection is associated with significant costs, already-vulnerable populations may suffer compounding costs as they balance self-protection against significant disease risk without adequate support from a broader community that does not share their risks (Atchison et al., Reference Atchison, Bowman, Vrinten, Redd, Pristerà, Eaton and Ward2021; Barrett et al., Reference Barrett, Bisset, Leidig, Marathe and Marathe2011; Jay et al., Reference Jay, Bor, Nsoesie, Lipson, Jones, Galea and Raifman2020; Skinner-Dorkenoo et al., Reference Skinner-Dorkenoo, Sarmal, Rogbeer, André, Patel and Cha2022). Further, structural inequities often leave population subgroups that are vulnerable to larger, more severe outbreaks with reduced access to protective measures like health education, treatment, vaccination, and paid leave (Cardona et al., Reference Cardona, Felipe, Fischer, Sehgal and Schwartz2021; Christensen et al., Reference Christensen, Pilling, Eyring, Dickerson, Sloan and Magnusson2020; Clouston et al., Reference Clouston, Yukich and Anglewicz2015; Dryhurst et al., Reference Dryhurst, Schneider, Kerr, Freeman, Recchia, van der Bles and van der Linden2020; Heymann et al., Reference Heymann, Sprague, Earle, McCormack, Waisath and Raub2021; Poteat et al., Reference Poteat, Millett, Nelson and Beyrer2020; Ridenhour et al., Reference Ridenhour, Sarathchandra, Seamon, Brown, Leung, Johnson-Leon and Johnson-Leung2022; Simione & Gnagnarella, Reference Simione and Gnagnarella2020; Williams & Cooper, Reference Williams and Cooper2020). Resulting differences in rates of protective behaviour uptake and effectiveness can compound disparities between groups and reduce the protective impact of awareness separation for more-vulnerable groups.

Epidemics are complex phenomena that typically involve heterogeneous mixing among groups of people that differ in biological and social risk factors, dynamic evolution of host behaviour, pathogen infectiousness and immune evasion, and ever-changing epidemiological and policy responses to real and perceived risk. Despite this range of potential drivers, we show here that a simple model that captures two key social processes – awareness-driven protective behaviour in a split population that can be separated in mixing and awareness – can drive many of the complex dynamics observed in emerging epidemics like Covid-19. For example, when awareness is uniform and mixing is separated, the group in which the pathogen is introduced later can experience second and third waves that exceed the initial wave in size (Figure 3). This trend resembles one observed in the United States during the first year of the Covid-19 pandemic, where certain regions where the virus was introduced early (e.g. New York City metropolitan area) experienced a large early wave and relatively few infections over the rest of the year, while other regions (e.g. the southern United States) generally had small early waves and larger second and third waves. Many hypotheses have been introduced to explain this phenomenon (e.g. policy, seasonal climate factors, and population density) and several factors may have contributed to this pattern (Y. Ma et al., Reference Ma, Pei, Shaman, Dubrow and Chen2021; Sy et al., Reference Sy, White and Nichols2021). Yet in our model these dramatic differences among populations in epidemic waves occur despite the groups being identical in transmission rates and disease outcomes and are entirely due to awareness-driven behaviour with uniform awareness among groups (Figure 3). Although the current analysis does not examine causation, and observed trends during Covid-19 probably involved a confluence of drivers, we have demonstrated how a simple behavioural process can qualitatively reproduce complex epidemic dynamics observed in real populations. To understand the extent of awareness separation in real populations and the role of specific behavioural processes in observed trends, our model could be parameterized using a combination of epidemiological, survey, mobility and social media data (Chang et al., Reference Chang, Pierson, Koh, Gerardin, Redbird, Grusky and Leskovec2021; Shen et al., Reference Shen, Yao, Zhang, Evans, Cao and Zhang2021; Weitz et al., Reference Weitz, Park, Eksin and Dushoff2020).

Feedback between vaccine efficacy and awareness-based vaccine uptake can also produce the counterintuitive scenario where vaccines that cause a greater reduction in transmission and mortality lead to more cumulative infections, even as deaths are reduced (Figure 5). If, as we assume here, protective behaviour is driven by awareness of severe outcomes like mortality, awareness separation may reduce differences in deaths between groups while widening differences in cases (Figures 4 and 5). The potential for awareness separation in vaccine uptake to reduce between-group differences in mortality is greatest when vaccination is introduced earlier in the epidemic, indicating that intervention timing may have health equity implications (Supplementary Figures S13 and S14). Accounting for awareness-based adoption of protective behaviour is therefore critical for understanding complicated epidemic dynamics such as plateaus and cycles (Figures 3 and 4), accurately deploying protective measures and assessing their impact across different diseases and population subgroups (Arthur et al., Reference Arthur, Jones, Bonds, Ram and Feldman2021; Steinegger et al., Reference Steinegger, Arola-Fernández, Granell, Gómez-Gardeñes and Arenas2022; Weitz et al., Reference Weitz, Park, Eksin and Dushoff2020).

Here we have considered arbitrarily defined groups that can be separated in mixing and awareness but initially differ only in the timing of pathogen introduction (Figures 1–3), fatality probability (Figures 4 and 5, Supplementary Figure S8), pathogen transmission (Supplementary Figure S6) or infectious period (Supplementary Figure S7). Real social groupings may fall along a number of social, demographic and geographic lines, while the assumption of two distinct and identifiable groups may not fully capture relevant social dynamics. The most relevant groupings with respect to awareness and disease risk may also depend on the disease. For infectious diseases that are generally more prevalent and severe in children (e.g. pertussis and measles), risk may depend on age while awareness is split between parents of young children vs. adults without children or among parents with different sentiments towards childhood vaccination (Bhattacharyya & Bauch, Reference Bhattacharyya and Bauch2010). In the context of Covid-19, disease burden and attitudes toward preventative measures (e.g. masks and vaccines) have differed markedly across age, socioeconomic status and race, and over time, demonstrating how intersecting and imperfectly overlapping identities may interact to determine attitudes, protective behaviours and risk (Maroko et al., Reference Maroko, Nash and Pavilonis2020; Schulz et al., Reference Schulz, Mehdipanah, Chatters, Reyes, Neblett and Israel2020; van Holm et al., Reference van Holm, Wyczalkowski and Dantzler2020). Moreover, ideological and social factors that do not correspond directly to disease risk (e.g. political affiliation) may influence decision-making and cause the level of protective behaviour in certain subgroups to diverge sharply from their relative risk for severe disease, potentially overcoming the effects of awareness separation (Christensen et al., Reference Christensen, Pilling, Eyring, Dickerson, Sloan and Magnusson2020; Grossman et al., Reference Grossman, Kim, Rexer and Thirumurthy2020). This process could be incorporated into our model by splitting the population into additional groups with respect to a cultural contagion or (mis)information spread process and allowing protective measures to be adopted based on awareness or contact with protective in-group members and rejected through fatigue or aversion to protective measures displayed by the opposite group (Mehta & Rosenberg, Reference Mehta and Rosenberg2020; Smaldino & Jones, Reference Smaldino and Jones2021).

Although we assumed that awareness was directly proportional to recent mortality, external influences like partisanship (Christensen et al., Reference Christensen, Pilling, Eyring, Dickerson, Sloan and Magnusson2020; Grossman et al., Reference Grossman, Kim, Rexer and Thirumurthy2020), media coverage (Shanta & Biswas, Reference Shanta and Biswas2020), misinformation (Lee et al., Reference Lee, Choi and Britt2021) and policy (Yan et al., Reference Yan, Malik, Bayham, Fenichel, Couzens and Omer2021) may alter the perception of risk or the adoption of protective measures at both the individual and group level. Group identification and assessment of relative risk may be unclear or inaccurate based on uncertainty at the beginning of the outbreak, misinformation about risk factors, a gradient in risk (e.g. increasing risk with age), lack of data stratification or unobserved risk factors. Attitudes based on one disease may carry over to another disease even if risk factors differ. Relative risk across groups may also vary across time and space, potentially leading to inaccurate assessment based on prior conditions: for example, a mild initial epidemic wave can mislead a group into believing they are inherently more protected and thereby relaxing protective behaviours. Cognitive interventions that increase the accuracy of individual risk perception, especially in high-risk groups, may help to reduce between-group differences in disease burden (Sinclair, Hakimi, et al., Reference Sinclair, Hakimi, Stanley, Adcock and Samanez-Larkin2021; Sinclair, Stanley, et al., Reference Sinclair, Stanley, Hakimi, Cabeza, Adcock and Samanez-Larkin2021).

Our model may also be extended to other scenarios involving a transmission process and collective behaviour, particularly social contagions like the spread of rumours and trends. Additional parameter space may be explored via the R Shiny interactive app accompanying this project, which currently only incorporates the non-pharmaceutical intervention model (https://mallory-harris.shinyapps.io/divided-disease/). Considering awareness separation as a social process that may interact with mixing, fatigue, waning immunity, pathogen evolution and pharmaceutical and non-pharmaceutical interventions may help to explain how humans are affected by and respond to infectious diseases in the presence of social divisions.

Acknowledgments

The authors thank Ann Marie Abraham, Marissa Childs, Lisa Couper, Isabel Delwel, Caroline Glidden, Andy Huynh, James Holland Jones, Devin Kirk, Kelsey Lyberger, Maike Morrison, Nicole Nova, Noah Rosenberg, Alyson Singleton, Eloise Skinner for their helpful feedback.

Author contributions

MJH and EAM conceived of project and designed models; MJH conducted analyses; EAM and MJH interpreted simulations and wrote the manuscript. KJC led development of the R Shiny app and KJC, EAM, and MJH edited the app. EAM, MJH, and KJC revised the manuscript.

Financial support

MJH was supported by the Knight-Hennessy Scholars Program and the Achievement Rewards for College Scientists Scholarship. KJC was supported by the Biology Summer Undergraduate Research Program at Stanford and the Stanford King Center on Global Development. EAM was supported by the National Science Foundation (DEB-2011147), with support from the Fogarty International Center, the National Institute of General Medical Sciences (R35GM133439), the National Institute of Allergy and Infectious Diseases (R01AI168097 and R01AI102918), the Terman Award, and seed grants from the Stanford King Center on Global Development, Woods Institute for the Environment and Center for Innovation in Global Health.

Conflicts of interests

All authors declare that there are no competing interests.

Research transparency and reproducibility

Code used to conduct these analyses are available on Github at: https://github.com/mjharris95/divided-disease

An interactive R Shiny app based on the non-pharmaceutical intervention model is available at: https://mallory-harris.shinyapps.io/divided-disease/

Supplementary material

To view supplementary material for this article, please visit https://doi.org/10.1017/ehs.2023.2

References

Abaluck, J., Kwong, L. H., Styczynski, A., Haque, A., Kabir, Md. A., Bates-Jefferys, E., …, Mobarak, A. M. (2022). Impact of community masking on COVID-19: A cluster-randomized trial in Bangladesh. Science, 375(6577), eabi9069. https://doi.org/10.1126/science.abi9069CrossRefGoogle ScholarPubMed
Acevedo-Garcia, D. (2000). Residential segregation and the epidemiology of infectious diseases. Social Science & Medicine (1982), 51(8), 11431161. https://doi.org/10.1016/s0277-9536(00)00016-2CrossRefGoogle ScholarPubMed
An, L., Hawley, S., Van Horn, M. L., Bacon, E., Yang, P., & Resnicow, K. (2020). Development of a coronavirus social distance attitudes scale. Patient Education and Counseling. https://doi.org/10.1016/j.pec.2020.11.027Google ScholarPubMed
Anthonj, C., Diekkrüger, B., Borgemeister, C., & Kistemann, Thomas. (2019). Health risk perceptions and local knowledge of water-related infectious disease exposure among Kenyan wetland communities. International Journal of Hygiene and Environmental Health, 222(1), 3448. https://doi.org/10.1016/j.ijheh.2018.08.003CrossRefGoogle ScholarPubMed
Arnold, C. R. K., Srinivasan, S., Rodriguez, S., Rydzak, N., Herzog, C. M., Gontu, A., …, Ferrari, M. J. (2022). A longitudinal study of the impact of university student return to campus on the SARS-CoV-2 seroprevalence among the community members. Scientific Reports, 12(1), 8586. https://doi.org/10.1038/s41598-022-12499-5CrossRefGoogle ScholarPubMed
Arthur, R. F., Jones, J. H., Bonds, M. H., Ram, Y., & Feldman, M. W. (2021). Adaptive social contact rates induce complex dynamics during epidemics. PLOS Computational Biology, 17(2), e1008639. https://doi.org/10.1371/journal.pcbi.1008639CrossRefGoogle ScholarPubMed
Atchison, C., Bowman, L. R., Vrinten, C., Redd, R., Pristerà, P., Eaton, J., & Ward, H. (2021). Early perceptions and behavioural responses during the COVID-19 pandemic: A cross-sectional survey of UK adults. BMJ Open, 11(1), e043577. https://doi.org/10.1136/bmjopen-2020-043577CrossRefGoogle ScholarPubMed
Barrett, C., Bisset, K., Leidig, J., Marathe, A., & Marathe, M. (2011). Economic and social impact of influenza mitigation strategies by demographic class. Epidemics, 3(1), 1931. https://doi.org/10.1016/j.epidem.2010.11.002CrossRefGoogle ScholarPubMed
Benfer, E. A., Vlahov, D., Long, M. Y., Walker-Wells, E., Pottenger, J. L., Gonsalves, G., & Keene, D. E. (2021). Eviction, health inequity, and the spread of COVID-19: Housing policy as a primary pandemic mitigation strategy. Journal of Urban Health, 98(1), 112. https://doi.org/10.1007/s11524-020-00502-1CrossRefGoogle ScholarPubMed
Bhattacharyya, S., & Bauch, C. T. (2010). A game dynamic model for delayer strategies in vaccinating behaviour for pediatric infectious diseases. Journal of Theoretical Biology, 267(3), 276282. https://doi.org/10.1016/j.jtbi.2010.09.005CrossRefGoogle ScholarPubMed
Brug, J., Aro, A. R., Oenema, A., de Zwart, O., Richardus, J. H., & Bishop, G. D. (2004). SARS Risk perception, knowledge, precautions, and information sources, the Netherlands. Emerging Infectious Diseases, 10(8), 14861489. https://doi.org/10.3201/eid1008.040283CrossRefGoogle ScholarPubMed
Calvin, R., Winters, K., Wyatt, S. B., Williams, D. R., Henderson, F. C., & Walker, E. R. (2003). Racism and cardiovascular disease in African Americans. The American Journal of the Medical Sciences, 325(6), 315331. https://doi.org/10.1097/00000441-200306000-00003CrossRefGoogle Scholar
Cardona, S., Felipe, N., Fischer, K., Sehgal, N. J., & Schwartz, B. E. (2021). Vaccination disparity: Quantifying racial inequity in COVID-19 vaccine administration in Maryland. Journal of Urban Health: Bulletin of the New York Academy of Medicine, 98(4), 464468. https://doi.org/10.1007/s11524-021-00551-0CrossRefGoogle ScholarPubMed
Chang, S., Pierson, E., Koh, P. W., Gerardin, J., Redbird, B., Grusky, D., & Leskovec, J. (2021). Mobility network models of COVID-19 explain inequities and inform reopening. Nature, 589(7840), 8287. https://doi.org/10.1038/s41586-020-2923-3CrossRefGoogle ScholarPubMed
Cheok, G. J. W., Gatot, C., Sim, C. H. S., Ng, Y. H., Tay, K. X. K., Howe, T. S., & Koh, J. S. B. (2021). Appropriate attitude promotes mask wearing in spite of a significant experience of varying discomfort. Infection, Disease & Health, 26(2), 145151. https://doi.org/10.1016/j.idh.2021.01.002CrossRefGoogle Scholar
Christensen, S. R., Pilling, E. B., Eyring, J. B., Dickerson, G., Sloan, C. D., & Magnusson, B. M. (2020). Political and personal reactions to COVID-19 during initial weeks of social distancing in the United States. PLOS ONE, 15(9), e0239693. https://doi.org/10.1371/journal.pone.0239693CrossRefGoogle ScholarPubMed
Clouston, S. A. P., Yukich, J., & Anglewicz, P. (2015). Social inequalities in malaria knowledge, prevention and prevalence among children under 5 years old and women aged 15–49 in Madagascar. Malaria Journal, 14(1), 499. https://doi.org/10.1186/s12936-015-1010-yCrossRefGoogle Scholar
Cubrich, M. (2020). On the frontlines: Protecting low-wage workers during COVID-19. Psychological Trauma: Theory, Research, Practice and Policy, 12(S1), S186S187. https://doi.org/10.1037/tra0000721Google Scholar
Dhewantara, P. W., Mamun, A. A., Zhang, W.-Y., Yin, W.-W., Ding, F., Guo, D., …, Soares Magalhães, R. J. (2018). Epidemiological shift and geographical heterogeneity in the burden of leptospirosis in China. Infectious Diseases of Poverty, 7, 57. https://doi.org/10.1186/s40249-018-0435-2CrossRefGoogle ScholarPubMed
Doherty, I. A., Schoenbach, V. J., & Adimora, A. A. (2009). Sexual mixing patterns and heterosexual HIV transmission among African Americans in the southeastern United States. Journal of Acquired Immune Deficiency Syndromes (1999), 52(1), 114120. https://doi.org/10.1097/QAI.0b013e3181ab5e10CrossRefGoogle ScholarPubMed
Dryhurst, S., Schneider, C. R., Kerr, J., Freeman, A. L. J., Recchia, G., van der Bles, A. M., …, van der Linden, S. (2020). Risk perceptions of COVID-19 around the world. Journal of Risk Research, 23(7–8), 9941006. https://doi.org/10.1080/13669877.2020.1758193CrossRefGoogle Scholar
Eksin, C., Shamma, J. S., & Weitz, J. S. (2017). Disease dynamics in a stochastic network game: A little empathy goes a long way in averting outbreaks. Scientific Reports, 7(1), Article 1. https://doi.org/10.1038/srep44122CrossRefGoogle Scholar
Farmer, P. (1996). Social inequalities and emerging infectious diseases. Emerging Infectious Diseases, 2(4), 259269. https://doi.org/10.3201/eid0204.960402CrossRefGoogle ScholarPubMed
Funk, S., Gilad, E., Watkins, C., & Jansen, V. A. A. (2009). The spread of awareness and its impact on epidemic outbreaks. Proceedings of the National Academy of Sciences, 106(16), 68726877. https://doi.org/10.1073/pnas.0810762106CrossRefGoogle ScholarPubMed
Funk, S., Salathé, M., & Jansen, V. A. A. (2010). Modelling the influence of human behaviour on the spread of infectious diseases: A review. Journal of The Royal Society Interface, 7(50), 12471256. https://doi.org/10.1098/rsif.2010.0142CrossRefGoogle ScholarPubMed
Gidengil, C. A., Parker, A. M., & Zikmund-Fisher, B. J. (2012). Trends in risk perceptions and vaccination intentions: A longitudinal study of the first year of the H1N1 pandemic. American Journal of Public Health, 102(4), 672679. https://doi.org/10.2105/AJPH.2011.300407CrossRefGoogle ScholarPubMed
Greene, S. K., Levin-Rector, A., Hadler, J. L., & Fine, A. D. (2015). Disparities in reportable communicable disease incidence by census tract-level poverty, New York City, 2006–2013. American Journal of Public Health, 105(9), e27e34. https://doi.org/10.2105/AJPH.2015.302741CrossRefGoogle ScholarPubMed
Grief, S. N., & Miller, J. P. (2017). Infectious disease issues in underserved populations. Primary Care, 44(1), 6785. https://doi.org/10.1016/j.pop.2016.09.011CrossRefGoogle ScholarPubMed
Grossman, G., Kim, S., Rexer, J. M., & Thirumurthy, H. (2020). Political partisanship influences behavioral responses to governors’ recommendations for COVID-19 prevention in the United States. Proceedings of the National Academy of Sciences, 117(39), 2414424153. https://doi.org/10.1073/pnas.2007835117CrossRefGoogle ScholarPubMed
Harris, M., Tessier-Lavigne, E., & Mordecai, E. (2021). The interplay of policy, behavior, and socioeconomic conditions in early COVID-19 epidemiology in Georgia. Journal of the Georgia Public Health Association, 8(2). https://doi.org/10.20429/jgpha.2021.080204CrossRefGoogle Scholar
Herrera-Diestra, J. L., & Meyers, L. A. (2019). Local risk perception enhances epidemic control. PLOS ONE, 14(12), e0225576. https://doi.org/10.1371/journal.pone.0225576CrossRefGoogle ScholarPubMed
Heymann, J., Sprague, A., Earle, A., McCormack, M., Waisath, W., & Raub, A. (2021). US sick leave in global context: US eligibility rules widen inequalities despite readily available solutions. Health Affairs, 40(9), 15011509. https://doi.org/10.1377/hlthaff.2021.00731CrossRefGoogle ScholarPubMed
Holtz, D., Zhao, M., Benzell, S. G., Cao, C. Y., Rahimian, M. A., Yang, J., …, Aral, S. (2020). Interdependence and the cost of uncoordinated responses to COVID-19. Proceedings of the National Academy of Sciences, 117(33), 1983719843. https://doi.org/10.1073/pnas.2009522117CrossRefGoogle ScholarPubMed
Jacquez, J. A., Simon, C. P., Koopman, J., Sattenspiel, L., & Perry, T. (1988). Modeling and analyzing HIV transmission: The effect of contact patterns. Mathematical Biosciences, 92(2), 119199. https://doi.org/10.1016/0025-5564(88)90031-4CrossRefGoogle Scholar
Jay, J., Bor, J., Nsoesie, E. O., Lipson, S. K., Jones, D. K., Galea, S., & Raifman, J. (2020). Neighbourhood income and physical distancing during the COVID-19 pandemic in the United States. Nature Human Behaviour, 4(12), Article 12. https://doi.org/10.1038/s41562-020-00998-2CrossRefGoogle ScholarPubMed
Ketcheson, D. I. (2021). Optimal control of an SIR epidemic through finite-time non-pharmaceutical intervention. Journal of Mathematical Biology, 83(1), 7. https://doi.org/10.1007/s00285-021-01628-9CrossRefGoogle ScholarPubMed
Lane, H. M., Morello-Frosch, R., Marshall, J. D., & Apte, J. S. (2022). Historical redlining is associated with present-day air pollution disparities in U.S. Cities. Environmental Science & Technology Letters, 9(4), 345350. https://doi.org/10.1021/acs.estlett.1c01012CrossRefGoogle ScholarPubMed
Lee, J., Choi, J., & Britt, R. K. (2021). Social media as risk-attenuation and misinformation-amplification station: How social media interaction affects misperceptions about COVID-19. Health Communication, 0(0), 111. https://doi.org/10.1080/10410236.2021.1996920Google Scholar
Li, Z., Wang, P., Gao, G., Xu, C., & Chen, X. (2016). Age-period-cohort analysis of infectious disease mortality in urban–rural China, 1990–2010. International Journal for Equity in Health, 15(1), 55. https://doi.org/10.1186/s12939-016-0343-7CrossRefGoogle ScholarPubMed
Ma, K. C., Menkir, T. F., Kissler, S., Grad, Y. H., & Lipsitch, M. (2021). Modeling the impact of racial and ethnic disparities on COVID-19 epidemic dynamics. ELife, 10, e66601. https://doi.org/10.7554/eLife.66601CrossRefGoogle ScholarPubMed
Ma, Y., Pei, S., Shaman, J., Dubrow, R., & Chen, K. (2021). Role of meteorological factors in the transmission of SARS-CoV-2 in the United States. Nature Communications, 12(1), 3602. https://doi.org/10.1038/s41467-021-23866-7CrossRefGoogle ScholarPubMed
Maroko, A. R., Nash, D., & Pavilonis, B. T. (2020). COVID-19 and inequity: A comparative spatial analysis of New York City and Chicago hot spots. Journal of Urban Health, 97(4), 461. Complementary Index.CrossRefGoogle ScholarPubMed
Mehta, R. S., & Rosenberg, N. A. (2020). Modelling anti-vaccine sentiment as a cultural pathogen. Evolutionary Human Sciences, 2, e21. https://doi.org/10.1017/ehs.2020.17CrossRefGoogle Scholar
Morris, D. H., Rossine, F. W., Plotkin, J. B., & Levin, S. A. (2021). Optimal, near-optimal, and robust epidemic control. Communications Physics, 4(1), 78. https://doi.org/10.1038/s42005-021-00570-yCrossRefGoogle Scholar
Omer, S. B., Enger, K. S., Moulton, L. H., Halsey, N. A., Stokley, S., & Salmon, D. A. (2008). Geographic clustering of nonmedical exemptions to school immunization requirements and associations with geographic clustering of pertussis. American Journal of Epidemiology, 168(12), 13891396. https://doi.org/10.1093/aje/kwn263CrossRefGoogle ScholarPubMed
Oraby, T., Thampi, V., & Bauch, C. T. (2014). The influence of social norms on the dynamics of vaccinating behaviour for paediatric infectious diseases. Proceedings of the Royal Society B: Biological Sciences, 281(1780), 20133172. https://doi.org/10.1098/rspb.2013.3172CrossRefGoogle ScholarPubMed
Perra, N., Balcan, D., Gonçalves, B., & Vespignani, A. (2011). Towards a characterization of behavior-disease models. PLOS ONE, 6(8), e23084. https://doi.org/10.1371/journal.pone.0023084CrossRefGoogle ScholarPubMed
Poteat, T., Millett, G. A., Nelson, L. E., & Beyrer, C. (2020). Understanding COVID-19 risks and vulnerabilities among black communities in America: The lethal force of syndemics. Annals of Epidemiology, 47, 13. https://doi.org/10.1016/j.annepidem.2020.05.004CrossRefGoogle ScholarPubMed
Pramasivan, S., Ngui, R., Jeyaprakasam, N. K., Liew, J. W. K., Low, V. L., Mohamed Hassan, N., …, Vythilingam, I. (2021). Spatial distribution of Plasmodium knowlesi cases and their vectors in Johor, Malaysia: In light of human malaria elimination. Malaria Journal, 20(1), 426. https://doi.org/10.1186/s12936-021-03963-0CrossRefGoogle ScholarPubMed
Quinn, S. C., Kumar, S., Freimuth, V. S., Musa, D., Casteneda-Angarita, N., & Kidwell, K. (2011). Racial disparities in exposure, susceptibility, and access to health care in the US H1N1 influenza pandemic. American Journal of Public Health, 101(2), 285293. https://doi.org/10.2105/AJPH.2009.188029CrossRefGoogle ScholarPubMed
Richardson, E. T., Malik, M. M., Darity, W. A., Mullen, A. K., Morse, M. E., Malik, M., …, Jones, J. H. (2021). Reparations for Black American descendants of persons enslaved in the U.S. and their potential impact on SARS-CoV-2 transmission. Social Science & Medicine, 276, 113741. https://doi.org/10.1016/j.socscimed.2021.113741CrossRefGoogle ScholarPubMed
Ridenhour, B. J., Sarathchandra, D., Seamon, E., Brown, H., Leung, F.-Y., Johnson-Leon, M., …, Johnson-Leung, J. (2022). Effects of trust, risk perception, and health behavior on COVID-19 disease burden: Evidence from a multi-state US survey. PLOS ONE, 17(5), e0268302. https://doi.org/10.1371/journal.pone.0268302CrossRefGoogle ScholarPubMed
Rothenberg, R., Muth, S. Q., Malone, S., Potterat, J. J., & Woodhouse, D. E. (2005). Social and geographic distance in HIV risk. Sexually Transmitted Diseases, 32(8), 506512.CrossRefGoogle ScholarPubMed
Schulz, A. J., Mehdipanah, R., Chatters, L. M., Reyes, A. G., Neblett, E. W., & Israel, B. A. (2020). Moving health education and behavior upstream: Lessons From COVID-19 for addressing structural drivers of health inequities. Health Education & Behavior, 47(4), 519524. https://doi.org/10.1177/1090198120929985CrossRefGoogle ScholarPubMed
Shanta, S. S., & Biswas, Md. H. A. (2020). The impact of media awareness in controlling the spread of infectious diseases in terms of SIR model. Mathematical Modelling of Engineering Problems, 7(3), 368376. https://doi.org/10.18280/mmep.070306CrossRefGoogle Scholar
Shen, L., Yao, R., Zhang, W., Evans, R., Cao, G., & Zhang, Z. (2021). Emotional attitudes of Chinese citizens on social distancing during the COVID-19 outbreak: Analysis of social media data. JMIR Medical Informatics, 9(3), e27079. https://doi.org/10.2196/27079CrossRefGoogle ScholarPubMed
Siggins, M. K., Thwaites, R. S., & Openshaw, P. J. M. (2021). Durability of Immunity to SARS-CoV-2 and other respiratory viruses. Trends in Microbiology, 29(9), 862. https://doi.org/10.1016/j.tim.2021.06.009CrossRefGoogle ScholarPubMed
Simione, L., & Gnagnarella, C. (2020). Differences between health workers and general population in risk perception, behaviors, and psychological distress related to COVID-19 spread in Italy. Frontiers in Psychology, 11. https://www.frontiersin.org/article/10.3389/fpsyg.2020.02166CrossRefGoogle ScholarPubMed
Sinclair, A. H., Hakimi, S., Stanley, M. L., Adcock, R. A., & Samanez-Larkin, G. R. (2021). Pairing facts with imagined consequences improves pandemic-related risk perception. Proceedings of the National Academy of Sciences, 118(32), e2100970118. https://doi.org/10.1073/pnas.2100970118CrossRefGoogle ScholarPubMed
Sinclair, A. H., Stanley, M. L., Hakimi, S., Cabeza, R., Adcock, R. A., & Samanez-Larkin, G. R. (2021). Imagining a personalized scenario selectively increases perceived risk of viral transmission for older adults. Nature Aging, 1(8), 677683. https://doi.org/10.1038/s43587-021-00095-7CrossRefGoogle ScholarPubMed
Skinner-Dorkenoo, A. L., Sarmal, A., Rogbeer, K. G., André, C. J., Patel, B., & Cha, L. (2022). Highlighting COVID-19 racial disparities can reduce support for safety precautions among White U.S. residents. Social Science & Medicine, 301, 114951. https://doi.org/10.1016/j.socscimed.2022.114951CrossRefGoogle ScholarPubMed
Smaldino, P. E., & Jones, J. H. (2021). Coupled dynamics of behaviour and disease contagion among antagonistic groups. Evolutionary Human Sciences, 3, e28. https://doi.org/10.1017/ehs.2021.22CrossRefGoogle Scholar
Soetaert, K., Petzoldt, T., & Setzer, R. W. (2010). Solving differential equations in R: Package deSolve. Journal of Statistical Software, 33, 125. https://doi.org/10.18637/jss.v033.i09CrossRefGoogle Scholar
Steinegger, B., Arola-Fernández, L., Granell, C., Gómez-Gardeñes, J., & Arenas, A. (2022). Behavioural response to heterogeneous severity of COVID-19 explains temporal variation of cases among different age groups. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences, 380(2214), 20210119. https://doi.org/10.1098/rsta.2021.0119Google Scholar
Sy, K. T. L., White, L. F., & Nichols, B. E. (2021). Population density and basic reproductive number of COVID-19 across United States counties. PLOS ONE, 16(4), e0249271. https://doi.org/10.1371/journal.pone.0249271CrossRefGoogle ScholarPubMed
Takahashi, T., Ellingson, M. K., Wong, P., Israelow, B., Lucas, C., Klein, J., …, Iwasaki, A. (2020). Sex differences in immune responses that underlie COVID-19 disease outcomes. Nature, 588(7837), article 7837. https://doi.org/10.1038/s41586-020-2700-3CrossRefGoogle ScholarPubMed
Toor, J., Echeverria-Londono, S., Li, X., Abbas, K., Carter, E. D., Clapham, H. E., …, Gaythorpe, K. A. (2021). Lives saved with vaccination for 10 pathogens across 112 countries in a pre-COVID-19 world. ELife, 10, e67635. https://doi.org/10.7554/eLife.67635CrossRefGoogle Scholar
van Holm, E. J., Wyczalkowski, C. K., & Dantzler, P. A. (2020). Neighborhood conditions and the initial outbreak of COVID-19: The case of Louisiana. Journal of Public Health, 16. https://doi.org/10.1093/pubmed/fdaa147Google Scholar
Weitz, J. S., Park, S. W., Eksin, C., & Dushoff, J. (2020). Awareness-driven behavior changes can shift the shape of epidemics away from peaks and toward plateaus, shoulders, and oscillations. Proceedings of the National Academy of Sciences, 117(51), 3276432771. https://doi.org/10.1073/pnas.2009911117CrossRefGoogle ScholarPubMed
Weston, D., Hauck, K., & Amlôt, R. (2018). Infection prevention behaviour and infectious disease modelling: A review of the literature and recommendations for the future. BMC Public Health, 18(1), 336. https://doi.org/10.1186/s12889-018-5223-1CrossRefGoogle ScholarPubMed
Williams, D. R., & Cooper, L. A. (2020). COVID-19 and health equity – A new kind of ‘herd immunity’. JAMA, 323(24), 24782480. https://doi.org/10.1001/jama.2020.8051CrossRefGoogle ScholarPubMed
Wu, X., Nethery, R. C., Sabath, M. B., Braun, D., & Dominici, F. (2020). Air pollution and COVID-19 mortality in the United States: Strengths and limitations of an ecological regression analysis. Science Advances, 6(45), eabd4049. https://doi.org/10.1126/sciadv.abd4049CrossRefGoogle ScholarPubMed
Yan, Y., Malik, A. A., Bayham, J., Fenichel, E. P., Couzens, C., & Omer, S. B. (2021). Measuring voluntary and policy-induced social distancing behavior during the COVID-19 pandemic. Proceedings of the National Academy of Sciences, 118(16), e2008814118. https://doi.org/10.1073/pnas.2008814118CrossRefGoogle ScholarPubMed
Zelner, J., Masters, N. B., Naraharisetti, R., Mojola, S. A., Chowkwanyun, M., & Malosh, R. (2022). There are no equal opportunity infectors: Epidemiological modelers must rethink our approach to inequality in infection risk. PLOS Computational Biology, 18(2), e1009795. https://doi.org/10.1371/journal.pcbi.1009795CrossRefGoogle ScholarPubMed
Zelner, J., Trangucci, R., Naraharisetti, R., Cao, A., Malosh, R., Broen, K., …, Delamater, P. (2020). Racial disparities in coronavirus disease 2019 (COVID-19) mortality are driven by unequal infection risks. Clinical Infectious Diseases, 72(5), e88e95. https://doi.org/10.1093/cid/ciaa1723CrossRefGoogle Scholar
Figure 0

Figure 1. Epidemic peaks are offset in time between groups when mixing is separated (C, D), and in magnitude when awareness is uniform but mixing is separated (C). Plots show the prevalence of infections over time in group a (pink) and group b (green) under four scenarios: awareness is uniform (A, C; $\epsilon = 0.5$) or separated (B, D; $\epsilon = 0.99$); mixing is uniform (A, B; h = 0.5) or separated (C, D; h = 0.99). We assume the pathogen is introduced only in group a at prevalence 0.001 and that all other parameters are equivalent between groups: transmission coefficient (β = 0.2), infectious period (1/ρ = 10), fatality probability (μ = 0.01), protective measure efficacy (κ = 0.3), responsiveness (θ = 100), memory (ℓ = 1), and fatigue (ϕ = 0). Lines overlap under uniform mixing (top row).

Figure 1

Figure 2. Separated awareness reduces between-group differences by reducing group b's awareness of the emerging epidemic and augmenting group a's response to the introduction of the pathogen. We initialize our model using the same parameters as Figure 1 with separated mixing (h = 0.99). We compare uniform awareness ($\epsilon = 0.5$; dashed lines) and separated awareness ($\epsilon = 0.99$; solid lines). At the top, we compare early time series (through t = 80) of (A) protective attitude prevalence in group a; (B) protective attitude prevalence in group b; (C) cumulative infections in group a; (D) cumulative infections in group b. (E) A phase portrait of protective attitude prevalence against cumulative infections in group a (pink) and group b (green). Points indicate values at t = 80, corresponding to the end of the time series in (A–D). Arrows indicate differences in protective attitude prevalence (grey) and cumulative infections (black) at t = 80 for separated vs. uniform awareness, with letters corresponding to time series panel labels.

Figure 2

Figure 3. Fatigue and long-term memory produce multiple epidemic peaks, which exceed the size of the initial peak in group b when uniform awareness and separated mixing leave that group with a high proportion of susceptible people following the first wave. We initialize the model with separated mixing (h = 0.99), long-term memory (ℓ = 30), and fatigue (ϕ = 0.02); all other parameters are the same as in Figure 1. We consider infections in group a (pink) and group b (green) over a longer time period (1000 days, compared with 200 days in Figure 1). The panels correspond to (A) uniform awareness ($\epsilon = 0.5$) and (B) separated awareness ($\epsilon = 0.99$).

Figure 3

Figure 4. Waning immunity and awareness-based vaccination drive epidemic cycles; separated awareness reduces the disparity in deaths (C vs. D) as more vulnerable group a members become vaccinated at a higher rate. We consider infections (A, B) and deaths (C, D) in the post-vaccine period in group a (pink) and group b (green) where the fatality probability for group a is double that of group b (μa = 0.02 and μb = 0.01). The x-axis gives time since vaccination began (tv = 200). We compare uniform awareness ($\epsilon = 0.5$) (A, C) and separated awareness ($\epsilon = 0.99$) (B, D). Other parameter values are: β = 0.2 (transmission coefficient), κ = 0.05 (transmission-reducing immunity), ζ = 0.05 (mortality-reducing immunity), ω = ϕ = 0.01 (waning immunity), infectious period (1/ρ = 10), θ = 20 (responsiveness), ℓ = 30 (memory), h = 0.99 (separated mixing), I0 = 0.0005 (initial infection prevalence). See Supplementary Figure S12 for a time series plot including the pre-vaccine period.

Figure 4

Figure 5. Greater immune protection (from vaccination and infection) leads to lower death rates (A), which in turn decreases vaccination rates (B) and increases infection rates (C); separated awareness reduces disparities in death rates (A) as groups are vaccinated at different rates proportional to their risks of death (B), creating differences in infection rates (C). We vary immune protection, defined as transmission-reducing immunity and mortality-reducing immunity, where both parameters are assigned the same values (κ = ζ). We assume immune protection is equivalent for vaccine- and infection-derived immunity. The x-axis is reversed because smaller values indicate stronger protection. We examine the impacts of stronger immune protection (lower values of κ and ζ) on total deaths (A), vaccinations (B) and infections (C) in the post-vaccine period (t = 200 through t = 2200). We consider the post-vaccine period to focus on the impacts of an awareness-based intervention administered under different levels of awareness separation. We compute each quantity for group a (pink) and group b (green) given uniform (dashed lines; $\epsilon = 0.5$) or separated (solid lines; $\epsilon = 0.99$) awareness. Other parameter values are the same as Figure 4.

Supplementary material: File

Harris et al. supplementary material

Harris et al. supplementary material

Download Harris et al. supplementary material(File)
File 2 MB