Wells-Riley model

Last updated

The Wells-Riley model is a simple model of the airborne transmission of infectious diseases, [1] [2] developed by William F. Wells and Richard L. Riley for tuberculosis [3] and measles. [4]

Contents

Wells-Riley can also be applied to other diseases transmitted in the air, such as COVID-19. [5] [6] [7] The model describes the situation where one or more infected people are sharing a room (and so the room air) with other people who are susceptible to infection. It makes predictions for the probability that a susceptible person becomes infected. The prediction is that infection is more likely for small poorly ventilated rooms, and if the infected person is highly infectious.

The Wells-Riley is a highly simplified model of a very complex process, but does at least make predictions for how the probability of infection varies with things within our control, such as room ventilation.

Description of model

Wells-Riley assumes that the air contains doses of the infectious bacterium or virus, and that you become infected if you breathe one dose in, i.e., the probability a person becomes infected, , is given by [1] [2]

This dose is not a single bacterium or virus, but however many are needed to start an infection. These infectious doses are sometimes called 'quanta' - no relation to quantum physics. The doses are breathed out or otherwise emitted by the infectious person into the air of the room, such that in the room there is a concentration of these doses per unit volume of the air. If you are breathing in air at volume rate , then after a time in the room, the

The Wells-Riley then relies on standard Poisson statistics which predicts for the probability of infection

after a time in the room. This is just the Poisson statistics expression for the probability of one or more doses being inhaled, once we know the mean number.

So the prediction is that you are more likely to become infected if the concentration of infectious doses in the room air is high, or if you spend longer in the room. The concentration of doses will tend to be high in small, poorly ventilated rooms, and smaller in larger, better ventilated rooms.

Relation of the Wells-Riley model to epidemiology

Note that the Wells-Riley model approaches transmission of an airborne diseases as a physical transport problem, i.e., as the problem of how a virus or bacterium gets from one human body to another.[ citation needed ]

For transmission of COVID-19, for example, this would be how a virus breathed out by an infected person, can cross a room and be breathed in by a susceptible person. This is a different approach from that taken in the epidemiology of infectious diseases, which may gather information about who (e.g., nurses, factory workers) becomes infected, in what situations (e.g., the home, factories), and understand the spread of a disease in those terms - without considering how a virus or bacterium actually gets from one person to another.[ citation needed ]

However, the probability of infection predicted by the Wells-Riley model is close to the attack rate (also called secondary attack rate, and note that this 'rate' is a probability not a rate) in epidemiology. Compare the definition of in this page with the definition of the attack rate.

Mechanism of transmission of infectious diseases through the air

Wells-Riley is only applicable for transmission directly via the air, not via the susceptible person picking up the infectious agent from a surface (fomite transmission). [1] Because the model assumes the air is well mixed it does not account for the region within one or two metres of an infected person, having a higher concentration of the infectious agent. Breathing and speaking produce a cone of warm (body temperature), humid air that moves out into and dissipates into the room air over a distance of about one to two metres, [8] while a sneeze, with its much faster moving air, produces air movement up to metres away. [9] If the person breathing/speaking/sneezing is infected then an infectious agent such as tuberculosis bacterium or a respiratory virus is expected to be more concentrated in this cone of air, but the infectious agent can also (at least in some cases) spread into the room air. [10]

Assumptions made by the Wells-Riley model

Estimating the number of inhaled doses requires more assumptions. The assumptions made by the model are essentially: [2] [1]

Assumptions 4 to 6 mean that the concentration of doses in the room air, , is

Doses can be removed in three ways:

Assuming we can add the rates of these processes

for the lifetime of the infectious agent in air, the lifetime of a dose in the air before settling onto a surface or the floor, and the lifetime of the dose before it is removed by room ventilation or filtration. Then the concentration of doses is

If the susceptible person spends a time inside the room and inhales air at a rate (volume per unit time) then they inhale a volume and so a number of infectious doses

or

Prediction equation

Putting all this together, the Wells-Riley prediction for the probability of infection is

where:

The lifetime of room air is just one over the air changes per hour - one measure of how well ventilated a room is. Building standards recommend several air changes per hour, in which case will be tens of minutes.

Relation to carbon dioxide concentration in a room

The Wells-Riley model assumes that an infected person continuously breathes out infectious virus. They will also continuously breathe out carbon dioxide, and so excess carbon dioxide concentration has been proposed as a proxy infection risk. [11] [12] [13] In other words, the higher the carbon dioxide concentration in a room, the higher the risk of infection by an airborne disease. The excess concentration of carbon dioxide is that over the background level in the Earth's atmosphere, which is assumed to come from human respiration (in the absence of another source such as fire). [11] [12] [13] Then the excess concentration of carbon dioxide is [13]

for people each exhaling carbon dioxide at a rate . Carbon dioxide neither sediments out (it is a molecule) nor decays, leaving ventilation as the only process that removes it. In the second equality we used , i.e., the rate of production of carbon dioxide is the breathing rate (volume of air exhaled per second = volume of air inhaled per second) times the concentration of carbon dioxide in exhaled breath . [13] Note that this implies that we can estimate how well ventilated a room is [13] if we know how many people are in the room, and the room's volume, from

If for the virus ventilation is the dominant route for removal of the virus, i.e., , the Wells-Riley prediction for the infection probability is then

which predicts that the higher the room concentration of carbon dioxide, the higher the infection risk.

Application to COVID-19

Although originally developed for other diseases such as tuberculosis, Wells-Riley has been applied [5] [14] to try and understand transmission of COVID-19, notably for a superspreading event in a chorale rehearsal in Skagit Valley (USA). [5]

The Wells-Riley model is implemented as an interactive Google Sheets spreadsheet, [6] and interactive apps showing estimates of the probability of infection. [15] Even for the simple Wells-Riley model, the infection probability, , depends on seven parameters. The probability of becoming infected is predicted to increase with how infectious the person is ( - which may peak around the time of the onset of symptoms and is likely to vary hugely from one infectious person to another, [2] [16] how rapidly they are breathing (which for example will increase with exercise), the length of the time they are in the room, as well as the lifetime of the virus in the room air.

This lifetime can be reduced by both ventilation and by removing the virus by filtration. Large rooms also dilute the infectious agent and so reduce risk - although this assumes that the air is well mixed - a highly approximate assumption. [2] [8] [10] A study of a COVID-19 transmission event in a restaurant in Guangzhou, went beyond this well-mixed approximation, to show that a group of three tables shared air with each other, to a greater extent than with the remainder of the (poorly ventilated) restaurant. [17] One infected person on one of these tables (a few metres apart) infected people on the other two tables.

The COVID-19 pandemic has led to work on improving the Wells-Riley model to account for factors such as the virus being in droplets of varying size which have varying lifetimes, [18] and an improved model [18] also has an interactive app. [19]

Related Research Articles

<span class="mw-page-title-main">Instanton</span> Solitons in Euclidean spacetime

An instanton is a notion appearing in theoretical and mathematical physics. An instanton is a classical solution to equations of motion with a finite, non-zero action, either in quantum mechanics or in quantum field theory. More precisely, it is a solution to the equations of motion of the classical field theory on a Euclidean spacetime.

The Feynman–Kac formula, named after Richard Feynman and Mark Kac, establishes a link between parabolic partial differential equations and stochastic processes. In 1947, when Kac and Feynman were both faculty members at Cornell University, Kac attended a presentation of Feynman's and remarked that the two of them were working on the same thing from different directions. The Feynman–Kac formula resulted, which proves rigorously the real-valued case of Feynman's path integrals. The complex case, which occurs when a particle's spin is included, is still an open question.

In medicine, public health, and biology, transmission is the passing of a pathogen causing communicable disease from an infected host individual or group to a particular individual or group, regardless of whether the other individual was previously infected. The term strictly refers to the transmission of microorganisms directly from one individual to another by one or more of the following means:

In electronics, when describing a voltage or current step function, rise time is the time taken by a signal to change from a specified low value to a specified high value. These values may be expressed as ratios or, equivalently, as percentages with respect to a given reference value. In analog electronics and digital electronics, these percentages are commonly the 10% and 90% of the output step height: however, other values are commonly used. For applications in control theory, according to Levine, rise time is defined as "the time required for the response to rise from x% to y% of its final value", with 0% to 100% rise time common for underdamped second order systems, 5% to 95% for critically damped and 10% to 90% for overdamped ones. According to Orwiler, the term "rise time" applies to either positive or negative step response, even if a displayed negative excursion is popularly termed fall time.

Mathematical models can project how infectious diseases progress to show the likely outcome of an epidemic and help inform public health and plant health interventions. Models use basic assumptions or collected statistics along with mathematics to find parameters for various infectious diseases and use those parameters to calculate the effects of different interventions, like mass vaccination programs. The modelling can help decide which intervention(s) to avoid and which to trial, or can predict future growth patterns, etc.

Compartmental models are a very general modelling technique. They are often applied to the mathematical modelling of infectious diseases. The population is assigned to compartments with labels – for example, S, I, or R,. People may progress between compartments. The order of the labels usually shows the flow patterns between the compartments; for example SEIS means susceptible, exposed, infectious, then susceptible again.

Nondimensionalization is the partial or full removal of physical dimensions from an equation involving physical quantities by a suitable substitution of variables. This technique can simplify and parameterize problems where measured units are involved. It is closely related to dimensional analysis. In some physical systems, the term scaling is used interchangeably with nondimensionalization, in order to suggest that certain quantities are better measured relative to some appropriate unit. These units refer to quantities intrinsic to the system, rather than units such as SI units. Nondimensionalization is not the same as converting extensive quantities in an equation to intensive quantities, since the latter procedure results in variables that still carry units.

<span class="mw-page-title-main">Continuous stirred-tank reactor</span> Type of chemical reactor

The continuous stirred-tank reactor (CSTR), also known as vat- or backmix reactor, mixed flow reactor (MFR), or a continuous-flow stirred-tank reactor (CFSTR), is a common model for a chemical reactor in chemical engineering and environmental engineering. A CSTR often refers to a model used to estimate the key unit operation variables when using a continuous agitated-tank reactor to reach a specified output. The mathematical model works for all fluids: liquids, gases, and slurries.

A cyclostationary process is a signal having statistical properties that vary cyclically with time. A cyclostationary process can be viewed as multiple interleaved stationary processes. For example, the maximum daily temperature in New York City can be modeled as a cyclostationary process: the maximum temperature on July 21 is statistically different from the temperature on December 20; however, it is a reasonable approximation that the temperature on December 20 of different years has identical statistics. Thus, we can view the random process composed of daily maximum temperatures as 365 interleaved stationary processes, each of which takes on a new value once per year.

In statistics, the Kendall rank correlation coefficient, commonly referred to as Kendall's τ coefficient, is a statistic used to measure the ordinal association between two measured quantities. A τ test is a non-parametric hypothesis test for statistical dependence based on the τ coefficient. It is a measure of rank correlation: the similarity of the orderings of the data when ranked by each of the quantities. It is named after Maurice Kendall, who developed it in 1938, though Gustav Fechner had proposed a similar measure in the context of time series in 1897.

The alveolar gas equation is the method for calculating partial pressure of alveolar oxygen. The equation is used in assessing if the lungs are properly transferring oxygen into the blood. The alveolar air equation is not widely used in clinical medicine, probably because of the complicated appearance of its classic forms. The partial pressure of oxygen in the pulmonary alveoli is required to calculate both the alveolar-arterial gradient of oxygen and the amount of right-to-left cardiac shunt, which are both clinically useful quantities. However, it is not practical to take a sample of gas from the alveoli in order to directly measure the partial pressure of oxygen. The alveolar gas equation allows the calculation of the alveolar partial pressure of oxygen from data that is practically measurable. It was first characterized in 1946.

<span class="mw-page-title-main">Rotational diffusion</span> Mechanics concept

Rotational diffusion is the rotational movement which acts upon any object such as particles, molecules, atoms when present in a fluid, by random changes in their orientations. Whilst the directions and intensities of these changes are statistically random, they do not arise randomly and are instead the result of interactions between particles. One example occurs in colloids, where relatively large insoluble particles are suspended in a greater amount of fluid. The changes in orientation occur from collisions between the particle and the many molecules forming the fluid surrounding the particle, which each transfer kinetic energy to the particle, and as such can be considered random due to the varied speeds and amounts of fluid molecules incident on each individual particle at any given time.

<span class="mw-page-title-main">All-pay auction</span>

In economics and game theory, an all-pay auction is an auction in which every bidder must pay regardless of whether they win the prize, which is awarded to the highest bidder as in a conventional auction. As shown by Riley and Samuelson (1981), equilibrium bidding in an all pay auction with private information is revenue equivalent to bidding in a sealed high bid or open ascending price auction.

In mathematics — specifically, in stochastic analysis — Dynkin's formula is a theorem giving the expected value of any suitably smooth function applied to a Feller process at a stopping time. It may be seen as a stochastic generalization of the (second) fundamental theorem of calculus. It is named after the Russian mathematician Eugene Dynkin.

In actuarial science and applied probability, ruin theory uses mathematical models to describe an insurer's vulnerability to insolvency/ruin. In such models key quantities of interest are the probability of ruin, distribution of surplus immediately prior to ruin and deficit at time of ruin.

<span class="mw-page-title-main">Airborne transmission</span> Disease transmission by airborne particles

Airborne transmission or aerosol transmission is transmission of an infectious disease through small particles suspended in the air. Infectious diseases capable of airborne transmission include many of considerable importance both in human and veterinary medicine. The relevant infectious agent may be viruses, bacteria, or fungi, and they may be spread through breathing, talking, coughing, sneezing, raising of dust, spraying of liquids, flushing toilets, or any activities which generate aerosol particles or droplets.

The concept of a minimal infective dose (MID), also known as the infectious dose, has traditionally been used for infectious microorganisms that contaminate foods. MID was defined as the number of microorganisms ingested from which a pathology is observed in the consumer. For example, to cause gastrointestinal disorders, the food must contain more than 100,000 Salmonella per gram or 1000 per gram for salmonellosis. however, some viruses like DHBV( duck hepatitis B virus) need as low as 9.5 x 10(9) virus per milliliters to cause liver infections.To know the dose ingested, it is also necessary to know the mass of the portion. This may be calculated using the following formula:

Viral phylodynamics is the study of how epidemiological, immunological, and evolutionary processes act and potentially interact to shape viral phylogenies. Since the term was coined in 2004, research on viral phylodynamics has focused on transmission dynamics in an effort to shed light on how these dynamics impact viral genetic variation. Transmission dynamics can be considered at the level of cells within an infected host, individual hosts within a population, or entire populations of hosts.

The residence time of a fluid parcel is the total time that the parcel has spent inside a control volume (e.g.: a chemical reactor, a lake, a human body). The residence time of a set of parcels is quantified in terms of the frequency distribution of the residence time in the set, which is known as residence time distribution (RTD), or in terms of its average, known as mean residence time.

<span class="mw-page-title-main">Transmission of COVID-19</span> Mechanisms that spread coronavirus disease 2019

The transmission of COVID-19 is the passing of coronavirus disease 2019 from person to person. COVID-19 is mainly transmitted when people breathe in air contaminated by droplets/aerosols and small airborne particles containing the virus. Infected people exhale those particles as they breathe, talk, cough, sneeze, or sing. Transmission is more likely the closer people are. However, infection can occur over longer distances, particularly indoors.

References

  1. 1 2 3 4 Sze To GN, Chao CY (February 2010). "Review and comparison between the Wells-Riley and dose-response approaches to risk assessment of infectious respiratory diseases". Indoor Air. 20 (1): 2–16. doi:10.1111/j.1600-0668.2009.00621.x. PMC   7202094 . PMID   19874402.
  2. 1 2 3 4 5 Noakes CJ, Beggs CB, Sleigh PA, Kerr KG (October 2006). "Modelling the transmission of airborne infections in enclosed spaces". Epidemiology and Infection. 134 (5): 1082–91. doi:10.1017/S0950268806005875. PMC   2870476 . PMID   16476170.
  3. Riley RL (January 2001). "What nobody needs to know about airborne infection". American Journal of Respiratory and Critical Care Medicine. 163 (1): 7–8. doi:10.1164/ajrccm.163.1.hh11-00. PMID   11208616.
  4. Riley EC, Murphy G, Riley RL (May 1978). "Airborne spread of measles in a suburban elementary school". American Journal of Epidemiology. 107 (5): 421–32. doi: 10.1093/oxfordjournals.aje.a112560 . PMID   665658.
  5. 1 2 3 Miller SL, Nazaroff WW, Jimenez JL, Boerstra A, Buonanno G, Dancer SJ, et al. (March 2021). "Transmission of SARS-CoV-2 by inhalation of respiratory aerosol in the Skagit Valley Chorale superspreading event". Indoor Air. 31 (2): 314–323. doi:10.1111/ina.12751. PMC   7537089 . PMID   32979298.
  6. 1 2 Jimenez JL, Peng Z. "COVID-19 Aerosol Transmission Estimator".
  7. Peng Z, Rojas AP, Kropff E, Bahnfleth W, Buonanno G, Dancer S, Kurnitski J, Li Y, Loomans M, Marr L, Morawska L (2022-01-18). "Practical Indicators for Risk of Airborne Transmission in Shared Indoor Environments and Their Application to COVID-19 Outbreaks". Environmental Science & Technology. 56 (2): 1125–1137. doi:10.1021/acs.est.1c06531. hdl: 10261/259515 . ISSN   0013-936X. PMID   34985868. S2CID   245771265.
  8. 1 2 Yang F, Pahlavan AA, Mendez S, Abkarian M, Stone HA (2020-12-01). "Towards improved social distancing guidelines: Space and time dependence of virus transmission from speech-driven aerosol transport between two individuals". Physical Review Fluids. 5 (12): 122501. Bibcode:2020PhRvF...5l2501Y. doi: 10.1103/PhysRevFluids.5.122501 . ISSN   2469-990X.
  9. Bourouiba L (May 2020). "Turbulent Gas Clouds and Respiratory Pathogen Emissions: Potential Implications for Reducing Transmission of COVID-19". JAMA. 323 (18): 1837–1838. doi: 10.1001/jama.2020.4756 . PMID   32215590.
  10. 1 2 Jones NR, Qureshi ZU, Temple RJ, Larwood JP, Greenhalgh T, Bourouiba L (August 2020). "Two metres or one: what is the evidence for physical distancing in covid-19?". BMJ. 370: m3223. doi: 10.1136/bmj.m3223 . PMID   32843355.
  11. 1 2 Peng Z, Jimenez JL (2021-05-11). "Exhaled CO2 as a COVID-19 Infection Risk Proxy for Different Indoor Environments and Activities". Environmental Science & Technology Letters. 8 (5): 392–397. doi:10.1021/acs.estlett.1c00183. PMC   8043197 . PMID   37566374.
  12. 1 2 Burridge HC, Fan S, Jones RL, Noakes CJ, Linden PF (2021-09-28). "Predictive and retrospective modelling of airborne infection risk using monitored carbon dioxide". Indoor and Built Environment. 31 (5): 1363–1380. doi: 10.1177/1420326X211043564 . hdl: 10044/1/91129 . ISSN   1420-326X. S2CID   244244848.
  13. 1 2 3 4 5 Rudnick SN, Milton DK (2003). "Risk of indoor airborne infection transmission estimated from carbon dioxide concentration: Indoor airborne transmission of infectious diseases". Indoor Air. 13 (3): 237–245. doi: 10.1034/j.1600-0668.2003.00189.x . PMID   12950586.
  14. Evans MJ (2020-05-25). "Avoiding COVID-19: Aerosol Guidelines". medRxiv   10.1101/2020.05.21.20108894 .
  15. Wei-Haas M, Elliot K (11 August 2020). "Measure the risk of airborne COVID-19 in your office, classroom, or bus ride". National Geographic . Archived from the original on February 19, 2021.
  16. Tang JW, Bahnfleth WP, Bluyssen PM, Buonanno G, Jimenez JL, Kurnitski J, et al. (April 2021). "Dismantling myths on the airborne transmission of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2)". The Journal of Hospital Infection. 110: 89–96. doi:10.1016/j.jhin.2020.12.022. PMC   7805396 . PMID   33453351.
  17. Li Y, Qian H, Hang J, Chen X, Cheng P, Ling H, et al. (June 2021). "Probable airborne transmission of SARS-CoV-2 in a poorly ventilated restaurant". Building and Environment. 196: 107788. doi:10.1016/j.buildenv.2021.107788. PMC   7954773 . PMID   33746341.
  18. 1 2 Bazant MZ, Bush JW (2021-04-27). "A guideline to limit indoor airborne transmission of COVID-19". Proceedings of the National Academy of Sciences. 118 (17): e2018995118. doi: 10.1073/pnas.2018995118 . ISSN   0027-8424. PMC   8092463 . PMID   33858987.
  19. "COVID-19 Indoor Safety Guideline". indoor-covid-safety.herokuapp.com. Retrieved 2021-04-24.