Environmental Transmission - Practice

Overview

This app allows you to explore a model which allows for both direct transmission and transmission through an environmental stage. Read about the model in the “Model” tab. Then do the tasks described in the “What to do” tab.

Learning Objectives

The Model

Model Overview

This is a simple SIR model with an additional environmental compartment.

This model has the following compartments:

The processes being modeled are:

Model Implementation

The flow diagram for this model is shown in this figure:

Model Diagram

Model Diagram

The set of ordinary differential equations (ODE) for this model are:

\[ \begin{aligned} \dot S & = n - b_{I}IS -b_PPS - mS\\ \dot I & = b_{I}IS + b_PPS -gI - mI\\ \dot R & = gI - mR \\ \dot P & = qI - c P \end{aligned} \]

A comment on terminology

It’s tempting to use the letter E for the environmental stage. However, the letter E is generally used for the exposed stage (what I generally call pre-symptomatic), which is then called a SEIR model. Thus, I’m using the letter P here. You might see models that use E for exposed or environmental stages. In general, there are no naming rules for either model compartments or parameters, so you will see all kinds of letters used. They should all be explicitly specified by the authors such that there is no ambiguity.

What to do

The tasks below are described in a way that assumes that everything is in units of MONTHS (rate parameters, therefore, have units of inverse months). If any quantity is not given in those units, you need to convert it first.

Task 1

Set the model parameters to the following values. 1000 initially susceptible individuals, 1 initially infected host, no pathogen in the environment, simulation duration 1 year. We’ll first look at direct transmission. Set transmission rate of bI = 0.004, environmental transmission bP = 0. Set environmental pathogen shedding and clearance rates to 0. Assume that the duration of the infectious period is 15 days long (and a month has 30 days). Turn off births and deaths for now. Run the simulation. You should get the usual direct transmission dynamics and a single outbreak. You should be left with around 203 susceptibles at the end of the simulation. Remember that in a model like this, where the only flow from the susceptible class is outflow through becoming infected, everyone who is not in the susceptible compartment at the end has been infected. Thus the total number of infected during the outbreak is the different between susceptibles at the beginning and at the end. Use one of the equations you have learned about in the reproductive number apps to compute R0.

Record

Task 2

Now try various values for the rate of shedding into the environment, q, and environmental clearance c. Leave everything else unchanged. As those parameters move away from 0, what do you expect to see? Run the simulation and compare your expectations with the results. Anything surprising happening? Do you understand why you see what you see?

Record

Task 3

Now set q = 100 and c = 100 . Turn off (set to 0) direct transmission. Run to make sure you don’t get an outbreak. Now turn on environmental transmission, set it to bP = 0.004, the same value you had for direct transmission above. Leave everything else unchanged. Run the simulation. Compare results with those from task 1. You might or might not find the result surprising. Take a close look at the curves for I and P.

Record

Task 4

In the previous task, removal of pathogens in the environment happened at a rate that was so fast that it essentially tracked whatever was going on in the infected population. In such a case, one can make with is called a quasi-steady state approximation for the fast variable and assume that it’s change only depends on the change of the variable it tracks. There is some formal way of doing it, a quick and often (but not always) correct way is to set the right side of the equation of the fast variable to 0. Here, that is the environmental variable, so the equation becomes 0=qI-cP and from that P=qI/c. Since we had q=c in the previous task, I and P tracked each other (approximately). That is generally not the case. Explore a more general scenario by setting q = 100 and c = 50. Have direct transmission off and environmental transmission as in the previous task. Since decay is now slower, pathogen in the environment builds up to higher levels and produces an overall larger outbreak. Look at the maximum values for P and I. Compute the ratio and compare that with the equation. Keep exploring how things change for different values of q and c.

Record

Task 5

It is possible to compute the reproductive number for enviromental transmission. Recall that it is defined as the average number of infectious hosts (of 1 type) that are infected by one infectious host of the same type. Here, we need to go through the environment. First, we need to compute how many units of environmental pathogen an individual host produces before they recover, and then how many hosts are infected per unit of environmental pathogen. This can be done as follows: One infected host sheds pathogen for a duration of 1/g at a rate q, so the total environmental pathogen produced by that host is q/g. This amount of pathogen is around for an average time of 1/c, during which time it infects new hosts at rate bP for a total production of new infected hosts bP/c. To get the number of new infectious hosts produced by one infectious host, one needs to multiply these two quantities. Then, to get the basic reproductive number, one also needs to multiply by the inital number of susceptible to get R0=qbPS/(gc)). Use the values from the last task (q=100, c=50) to compute R0 for this simulation. Then use the final-size equation to compute the fraction of susceptible you expect to get for this R0 value and confirm with the simulation.

Record

Task 6

Use the settings from the previous task and alternate between starting with some non-zero number of infected individuals and no pathogen in the environment, or no infected individuals and some pathogen in the environment. Convince yourself that you always get an outbreak, and if it’s not clear why, look at the flow diagram for the model to understand which processes lead to more or less the same outbreak if you start with either non-zero infected individuals or non-zero levels of environmental contamination. Then also turn on direct transmission, such that you now have transmission through both pathways (environmental and direct) and explore how different values affect the outbreak.

Record

Task 7

Keep exploring by trying different parameters and transmission settings and see how they influence results. You can also turn on births/deaths and see what you get. As you continue your exploration, think about potential real infectious diseases that might be approximated by either one of the transmission types. You could try to find parameter values in the literature that could approximate that disease and run simulations mimicking a specific disease and setting.

Record

Further Information

References

Codeço, C T. 2001. Endemic and Epidemic Dynamics of Cholera: The Role of the Aquatic Reservoir. BMC Infectious Diseases 1: 1.