🧬 Overview and Conceptual Motivation
Within the broad spectrum of mathematical approaches to epidemic analysis, the Gompertz Growth Epidemic Model occupies a distinctive position as a phenomenological framework. Rather than explicitly modeling interactions between susceptible and infectious individuals, this model focuses on the observed growth trajectory of cumulative cases. Its defining feature is that the relative growth rate of the epidemic decays exponentially over time. This makes the Gompertz model particularly effective for describing outbreaks that display early sub-exponential growth, rapid saturation, or strong slowing due to behavioral change, intervention, or population structure.
🏗️ Model Structure and Conceptual Flow
As a growth-based phenomenological model, the Gompertz formulation does not rely on a multi-compartment disease progression. Instead, it tracks the accumulation of reported cases over time through a single state variable.
Cumulative Incidence (C)
The total number of confirmed cases observed since the beginning of the outbreak.
Growth Process
The epidemic initially grows at a rate governed by an intrinsic growth parameter, but this rate is continuously dampened by an internal decay mechanism. This deceleration can represent depletion of susceptible clusters, behavioral adaptation, non-pharmaceutical interventions, or an effective carrying capacity of the population.
Flow structure:
C₀ → C(t) → K
Here, C₀ is the initial number of cases and K is the asymptotic final epidemic size. As time progresses, the effective growth rate decreases and eventually approaches zero.
🧮 Mathematical Formulation
The Gompertz epidemic model is commonly expressed as an ordinary differential equation describing the time evolution of cumulative cases.
Rate of change of cumulative cases
dC/dt = r · C · exp(−a · t)
An equivalent formulation that explicitly includes the final epidemic size K is:
dC/dt = r · C · [ ln(K) − ln(C) ]
In these equations, C denotes cumulative incidence at time t, r is the intrinsic initial growth rate, a is the exponential decay constant governing how rapidly growth slows, and K is the carrying capacity or final size of the epidemic.
📋 Table 1. Definition of Model Parameters
| Parameter | Symbol | Definition |
|---|---|---|
| Cumulative cases | C | Total number of reported infections at time t |
| Initial growth rate | r | Intrinsic early epidemic growth rate |
| Growth decay constant | a | Rate at which relative growth declines over time |
| Carrying capacity | K | Final epidemic size or saturation level |
| Time | t | Continuous time variable |
🌤️ Climatic Variable Integration and Environmental Sensitivity
Environmental conditions can influence epidemic growth by modifying transmission efficiency and host susceptibility. In the Gompertz framework, such effects are typically incorporated by allowing the growth parameter r to vary with time and climate.
A representative formulation is:
r(W) = r₍base₎ · [1 + η · sin(2π (t − φ) / 365)] · f(W)
In this expression, r₍base₎ is the baseline growth rate under optimal conditions, η represents the amplitude of seasonal forcing, and φ determines the seasonal phase. The function f(W) captures the response to a climatic variable W, such as temperature or absolute humidity. A common choice is
f(W) = exp(−k · |T − T₍opt₎|),
where T₍opt₎ denotes the optimal temperature for viral transmission efficiency. This structure allows seasonal acceleration or deceleration of epidemic growth to be represented within an otherwise simple growth model.
📊 Table 2. Realistic Parameter Ranges for Viral Epidemics
| Quantity | Typical Range | Interpretation |
|---|---|---|
| Initial growth rate (r) | 0.15 – 0.55 day⁻¹ | Early epidemic expansion speed |
| Growth decay constant (a) | 0.01 – 0.15 day⁻¹ | Strength of epidemic slowing |
| Final epidemic size (K) | 10³ – 10⁶ | Depends on population and outbreak scale |
| Initial doubling time | ln(2) / r | Early-phase growth metric |
🎯 Applicability and Limitations
Applicability
The Gompertz model is widely used for early-stage epidemic forecasting when biological parameters such as latent periods or contact rates are uncertain. It is particularly effective for describing sub-exponential growth in structured or behaviorally constrained populations and for generating rapid, top-down projections of peak timing and healthcare demand.
Key Assumptions and Weaknesses
The model does not explicitly represent biological mechanisms such as susceptibility, immunity, or recovery. It assumes a single epidemic wave with a fixed asymptote and therefore cannot naturally capture recurrent or multi-wave dynamics without re-estimation of parameters. In addition, the inherent asymmetry of the Gompertz curve, with peak growth occurring when cumulative cases reach approximately 37% of K, may not align with all epidemic trajectories.
📚 References
- Gompertz, B. (1825). On the Nature of the Function Expressive of the Law of Human Mortality. Philosophical Transactions of the Royal Society.
- Winsor, C. P. (1932). The Gompertz Curve as a Growth Curve. Proceedings of the National Academy of Sciences.
- Chowell, G., et al. (2016). Characterizing the reproduction number of epidemics with generalized growth models. Physics of Life Reviews.
- Ohnishi, T., et al. (2011). On the Gompertz Curve as a Model for an Infectious Disease. Journal of Theoretical Biology.
🧠 Analogy for Intuition
The Gompertz model can be compared to a car coasting uphill while running out of fuel. The initial growth rate r is the car’s starting speed. The decay constant a represents friction, gravity, and fuel depletion acting together to slow the vehicle. Eventually, the car comes to rest at a specific height, corresponding to the final epidemic size K. The internal mechanics of the engine need not be known to predict where the car will ultimately stop.