📈 Conceptual Overview
Vector-borne infectious diseases such as Dengue, Zika, and Malaria require the simultaneous modeling of two biologically distinct populations: a vertebrate host and an arthropod vector. The Bailey–Dietz model extends the classical Ross–Macdonald framework by providing a clear system of ordinary differential equations that explicitly capture the bidirectional transmission cycle between humans and vectors.
This model is a foundational tool in mathematical epidemiology for understanding how infection is sustained through repeated host–vector–host interactions and how environmental drivers influence transmission intensity.
🏗️ Compartmental Structure and Transmission Flow
The Bailey–Dietz model tracks two interacting populations with different life histories and epidemiological roles.
Human Population (H)
The human population is divided into three epidemiological states:
• Susceptible (Sₕ): Individuals who can acquire infection
• Infected (Iₕ): Individuals capable of transmitting infection to vectors
• Recovered (Rₕ): Individuals with immunity
The total human population size Nₕ is typically assumed constant over short epidemic time scales.
Vector Population (V)
The vector population is divided into two states:
• Susceptible (Sᵥ): Vectors capable of acquiring infection
• Infected (Iᵥ): Vectors capable of transmitting infection
Vectors do not recover and remain infectious until death due to their short lifespan.
Transmission Cycle
• Human → Vector: A susceptible vector becomes infected after biting an infected human
• Vector → Human: An infected vector transmits the pathogen while biting a susceptible human
🧮 Mathematical Formulation
The dynamics of the Bailey–Dietz model are governed by the following system of ordinary differential equations.
Human Population Dynamics
dSₕ(t) / dt
= μₕ Nₕ − ( a · b · Iᵥ(t) / Nₕ ) Sₕ(t) − μₕ Sₕ(t)
dIₕ(t) / dt
= ( a · b · Iᵥ(t) / Nₕ ) Sₕ(t) − ( r + μₕ ) Iₕ(t)
dRₕ(t) / dt
= r Iₕ(t) − μₕ Rₕ(t)
Vector Population Dynamics
dSᵥ(t) / dt
= μᵥ Nᵥ − ( a · c · Iₕ(t) / Nₕ ) Sᵥ(t) − μᵥ Sᵥ(t)
dIᵥ(t) / dt
= ( a · c · Iₕ(t) / Nₕ ) Sᵥ(t) − μᵥ Iᵥ(t)
These equations explicitly link host and vector infection processes through biting-mediated transmission.
🌤️ Climatic Variable Integration: Temperature-Dependent Transmission
Transmission intensity in vector-borne diseases is strongly influenced by environmental temperature. The biting rate a and vector mortality μᵥ are frequently modeled as temperature-dependent functions.
A commonly used formulation for the biting rate is the Brière function:
a(T)
= c · T · ( T − Tₘᵢₙ ) · √( Tₘₐₓ − T )
where Tₘᵢₙ and Tₘₐₓ represent the lower and upper thermal limits for vector activity.
The Extrinsic Incubation Period (EIP), representing the time required for a pathogen to become transmissible inside the vector, is often modeled using a degree-day relationship:
EIP(T)
= DD / ( T − Tₜₕᵣₑₛₕₒₗd )
This formulation links temperature directly to the speed of pathogen development within the vector.
📋 Table 1. Parameter Definitions
| Parameter | Definition |
|---|---|
| Sₕ | Susceptible human population |
| Iₕ | Infected human population |
| Rₕ | Recovered human population |
| Nₕ | Total human population |
| Sᵥ | Susceptible vector population |
| Iᵥ | Infected vector population |
| Nᵥ | Total vector population |
| a | Vector biting rate |
| b | Transmission probability from vector to human |
| c | Transmission probability from human to vector |
| r | Human recovery rate |
| μₕ | Human natural mortality rate |
| μᵥ | Vector mortality rate |
| T | Ambient temperature |
| DD | Required degree-days for pathogen development |
| Tₜₕᵣₑₛₕₒₗd | Developmental temperature threshold |
📊 Table 2. Typical Parameter Ranges
| Parameter | Typical Range | Interpretation |
|---|---|---|
| a | 0.1 – 0.5 day⁻¹ | Average bites per vector per day |
| b | 0.2 – 0.8 | Vector-to-human transmission probability |
| c | 0.1 – 0.5 | Human-to-vector transmission probability |
| r | 0.05 – 0.2 day⁻¹ | Human infectious period of 5–20 days |
| μᵥ | 0.05 – 0.3 day⁻¹ | Vector lifespan of approximately 3–20 days |
| m | 1 – 100 | Vector-to-human ratio Nᵥ / Nₕ |
| Tₘᵢₙ | 10 – 18 °C | Lower thermal limit for vector activity |
| Tₘₐₓ | 32 – 40 °C | Upper thermal limit for vector activity |
🎯 Applicability and Limitations
Applicability
• Modeling Dengue, Yellow Fever, West Nile Virus, and similar arboviruses
• Evaluating vector control strategies such as insecticides or bed nets
• Studying seasonal transmission patterns driven by climate variability
Key Assumptions and Weaknesses
• Assumes constant host and vector populations
• Does not include explicit latent periods in its simplest form
• Ignores immature vector life stages such as larvae and pupae
• Assumes homogeneous mixing of hosts and vectors
Despite these limitations, the Bailey–Dietz model remains a central framework for understanding and controlling vector-borne epidemics.
📚 References
- Bailey, N. T. J. (1982). The Biomathematics of Malaria. Charles Griffin & Co Ltd.
- Dietz, K. (1975). Transmission and control of arbovirus diseases. Epidemiology.
- Ross, R. (1911). The Prevention of Malaria. John Murray.
- Macdonald, G. (1957). The Epidemiology and Control of Malaria. Oxford University Press.
- Mordecai, E. A., et al. (2013). Optimal temperature for malaria transmission is dramatically lower than previously predicted. Ecology Letters.
🔍 Analogy for Clarity
The Bailey–Dietz model can be visualized as a revolving door connecting two buildings—humans and vectors. For a pathogen to move, it must pass repeatedly through the door. The speed of the door depends on environmental conditions: when temperature is optimal, the door spins rapidly; when conditions are poor, transmission slows or stops entirely.