This paper introduces a time-since-recovery structured, multi-strain, multi-population model of avian influenza. Influenza A viruses infect many species of wild and domestic birds and are classified into two groups based on their ability to cause disease: low pathogenic avian influenza (LPAI) and high pathogenic avian influenza (HPAI). Prior infection with LPAI provides partial immunity towards HPAI. The model introduced in this paper structures LPAI-recovered birds (wild and domestic) with time-since-recovery and includes crossimmunity towards HPAI that can fade with time. The model has a unique disease-free equilibrium (DFE), unique LPAI-only and HPAIonly equilibria and at least one coexistence equilibrium. We compute the reproduction numbers of LPAI (R L ) and HPAI (R H ) and show that the DFE is locally asymptotically stable when R L < 1 and R H < 1. A unique LPAI-only (HPAI-only) equilibrium exists when R L > 1 (R H > 1) and it is locally asymptotically stable if HPAI (LPAI) cannot invade the equilibrium, that is, if the invasion numberR H L < 1 (R L H < 1). We show using numerical simulations that the ODE version of the model, which is obtained by discarding the time-since-recovery structures (making cross-immunity constant), can exhibit oscillations, and also that the pathogens LPAI and HPAI can coexist with sustained oscillations in both populations. Through simulations, we show that even if both populations (wild and domestic) are sinks when alone, LPAI and HPAI can persist in both populations combined. Thus, reducing the reproduction numbers of LPAI and HPAI in each population to below unity is not enough to eradicate the disease. The pathogens can continue to coexist in both populations unless transmission between the populations is reduced.
ARTICLE HISTORY