Modelling Foot and Mouth Disease in the Context of Active Immigrants ()
1. Introduction
In our interconnected world, the intersection of global health concerns and human mobility has captured the attention of researchers, policymakers, and public health experts [1] . The complex patterns of infectious disease transmission, particularly in the context of active immigrants, present intricate challenges demanding deep comprehension and innovative solutions [2] . Foot and Mouth Disease (FMD) serves as a prime example, a highly contagious viral infection bearing substantial economic consequences for livestock industries [3] [4] . Simultaneously, the mobility of active immigrants seeking improved lives abroad further complicates disease spread [5] . This introduction explores the interplay of FMD, active immigrants, and the potent role of mathematical modeling in bridging gaps and shaping effective disease management strategies.
Foot and Mouth Disease (FMD) serves as a vivid illustration of the intricate links between animal well-being, economic strength, and global disease propagation [6] [7] . Affecting cloven-hoofed animals like cattle, pigs, sheep, and goats, FMD causes mouth and hoof blisters, inflicting both substantial suffering and economic setbacks [8] [9] . Rapid action is essential due to its rapid transmission, often involving large-scale culling to contain the outbreak [10] [11] . Notably, the impact of animal mobility, especially among active immigrants, in amplifying disease transmission has gained significant research attention.
Immigrants seeking improved lives and economic prospects can unintentionally facilitate the spread of diseases due to their close community ties and cross-border mobility [12] [13] . Recognizing their role is crucial for developing effective strategies to prevent and control diseases like FMD.
Mathematical modeling provides an advanced method to understand how diseases spread with active immigrants [14] . These models quantify interactions between populations, movement, and disease, revealing outbreak possibilities, high-risk areas, and intervention effects [15] [16] . They also fill research gaps by evaluating strategies like vaccination, culling, quarantine, treatment, surveillance, and border controls.
However, despite mathematical modeling’s potential, significant research gaps remain regarding the interplay between FMD and active immigrants. These gaps encompass more precise human movement data, inclusion of socio-cultural influences on transmission, and adaptable models for changing disease dynamics [17] . Addressing these requires teamwork among epidemiologists, mathematicians, social scientists, and policymakers to create comprehensive prevention and control strategies [18] .
In our study of Foot and Mouth Disease, active immigrants, and mathematical modeling, we aim to enhance our grasp of disease transmission dynamics. Through addressing research gaps and using mathematical models, we intend to offer valuable insights into global health security complexities, presenting creative strategies to reduce the impact of infectious diseases on animals and humans.
2. Materials and Methods
2.1. Model Formulation
In this section, we outline the mathematical representation of the FMD compartmental model using a system of differential equations. The model captures FMD transmission dynamics, categorized into five compartments: Susceptible
, Exposed
, Asymptomatic
, Symptomatic
and Recover
. So,
. The corresponding flow diagram is depicted in Figure 1 and the corresponding model system of differential equations is presented in (1).
2.2. Model Diagram
From Figure 1, the transition between compartments can now be expressed into five non-linear differential equations defined as follows:
(1)
where
is the force of infection.
· Animals are recruited at a constant rate
and are assumed to be susceptible to infection. Animals in all compartments suffer natural mortality at rate
and
is a disease induced dearth rate.
Figure 1. Model flow chart for the dynamics of Foot and Mouth disease with active immigrants.
·
stands for force of infection, where
is the effective contact rate and
is a modification parameter.
·
and
represent progression rates from exposed to symptomatic and asymptomatic classes respectively.
·
and
are recovery rates for symptomatic and asymptomatic animals, leading to temporary immunity.
·
and
are varying influx rates for asymptomatic and symptomatic animals from other regions, with
and
.
The state variables and model parameters for the given model system are presented in Table 1 and Table 2 respectively.
Table 1. The model variables and their representations
Table 2. Parameters and their descriptions.
2.3. Model Analysis
Non-Negativity and Boundness of Model System (1)
Theorem 1. For the model system (1), there exists a unique solution in
, however, the solution is always positive for all value of
and remains in
.
Proof. From the first equation of the model equations, we have
, integrating with initial condition
yields:
(2)
From the second equation of the model equations, we have:
,
integrating again with initial condition
yields:
(3)
In a similar way, the rest of the equations of the model Equation (1) with initial conditions,
,
, and
gives:
(4)
Hence, all the solutions of model Equation (1) are non negative for all
. ¨
Theorem 2. Let
be the unique solution of the model system (1) for all
, then the solution
is bounded above, that is,
where Ω is the feasible region defined as:
which is interior denoted by int(Ω) and given by:
Proof. Here, we prove that the solutions of model system (1) are bounded for all
. Biologically, the lowest possible value of each states of the model system (1) is zero. Next we determine the upper-bound of states. Given that,
:
substitute the values of differential equations from the model system (1), simplification gives:
(5)
By integration gives:
(6)
Since the total population
is positive for all
. It is well defined that:
Therefore,
, for all
.
This proves that all solutions of the FMD model with initial conditions in Ω remain in Ω for all
.¨
2.4. Disease-Free Equilibrium Point (DFE)
In this section, we present the disease-free equilibrium point of the model system (1). We first set the right-hand side of the model system (1) equal to zero as follows:
(7)
In the absence of disease in the population, we have
and the model (7) admits a trivial equilibrium point commonly known as the disease-free equilibrium point. Thus, the model (7) has a disease-free equilibrium E0 given by:
(8)
The disease-free equilibrium point in disease modeling is important as it is used when computing the basic reproduction number
.
2.5. Basic Reproduction Number ( R 0 )
Following the next generation matrix approach as used by [23] [24] [25] , the non-negative matrix
that denotes the generation of new infection in the population and non-singular matrix
that denotes the disease transfer of infected compartments evaluated at the disease-free E0 are given as follows:
(9)
Thus, using the next-generation matrix approach in (9), it can easily be verified that the model system (1) has basic reproduction number
given by:
(10)
The threshold value
is an important epidemiological quantity in disease modeling. It determines the strength of transmission potential of a disease to invade the community. It has been scientifically demonstrated that when
the disease dies out in the population. However, the disease persist when
.
2.6. Local Stability of Disease-Free Equilibrium (DFE)
Here, we establish the local stability of disease-free equilibrium using trace and determinant method and Routh Hurwitz’s stability criterion.
Theorem 3. The DFE of the model system (1) at DFE (E0) is locally asymptotically stable if
and unstable
.
Proof. In proving this theorem, the Jacobian matrix of model system (1) at DFE (E0) is given as:
We have
,
negative, while the other three eigenvalues can be computed through the following cubic equations. The remaining eigenvalues are obtained from:
where
,
and
. The characteristic polynomial for matrix G is given by
, Whereby;
,
,
,
.
Here, we used Routh-Hurwitz criterion, the Routh-Hurwitz Criterion is a mathematical tool used to assess the stability of equilibrium points in a dynamical system, particularly in the context of linear time-invariant systems. In the context of epidemiological models, if all the coefficients in the first column of the Routh array have the same sign, and that sign is positive, then the system is locally asymptotically stable. Conversely, if there is at least one sign change in the first column, the system is locally unstable. A polynomial
has roots with a negative real part if it satisfies the following necessary and sufficient conditions:
1)
, 2)
, 3)
and 4)
.
Then,
,
. Again,
,
.
It can be shown that:
The value of
is given by:
Thus, the condition
holds if
. Therefore, it implies that DFE (E0) is locally asymptotically stable when
and unstable otherwise.¨
2.7. Global Stability of DFE
Analysis of the global stability of the disease-free equilibrium solution of the model system (1) is done by using the approach by Castillo-Chavez et al. (2002). The model system can be written as:
(11)
From Equation (11),
and
indicates the number of uninfected and infected cattle respectively. F and B represent the functions of infected and uninfected individuals respectively. The disease-free equilibrium solutions is denoted by
:
For DFE to be globally asymptotically stable, the following conditions must be satisfied simultaneously:
For which, the Jacobian matrix A is a Metzler matrix and Ω is the region where the model system yields a biologically meaningful. Now, if the model system (11) satisfies the above conditions, then Theorem 4 holds.
Theorem 4. For
, the disease-free equilibrium point E0 of model system (7) is globally asymptotically stable whenever
.
Proof. Condition 1 in (12) is proved as follows:
At DFE,
gives:
It is clear that the solutions
, and
as
regardless of the value of
and
. Implying that E0 is global asymptotically stable.
Condition 2 in (12) is proved as follows:
At DFE
gives:
It is clear that matrix A is Metzler matrix. Then,
Since
, it is clear that
.
Hence,
is globally asymptotically stable.¨
2.8. Endemic Equilibrium Point
Endemic equilibrium point denoted by
is obtained by solving for
from the model system (14):
(14)
whereby
2.8.1. Bifurcation Analysis.
Bifurcation analysis was done by using Center Manifold theory, let
,
,
,
,
and
, by using vector notation
(where T denotes the transpose). Then, the model system (1) can be written in the form
, with
. Then,
(15)
Choose
as a bifurcation parameter. Consider when
, then solving for
using
gives:
The Jacobian matrix at DFE (E0) and bifurcation parameter
is given by:
where
,
,
,
,
and
.
The right eigenvector associated with the zero eigenvalue of the Jacobian matrix
,
is defined such that
at
given by;
,
,
,
, and
.
The left eigenvector
such that
yields:
and
.
Theorem 5. Consider the general system of ordinary differential equations with a parameter
such that:
(16)
such that
, whereby 0 is an equilibrium point of the system with the following conditions:
1)
is the linearization matrix of the system (3.33) around the equilibrium 0 with
evaluated at 0.
2) Zero is a simple eigenvalue of M, and the rest eigenvalue of M have negative real parts.
3) Matrix M has a right eigenvector w and a left eigenvector v corresponding to the zero eigenvalue.
Let
be the kth component of f and
(17)
The local dynamics of (3.33) around 0 are totally determined by a and b.
1)
,
. When
with
, 0 is locally asymptotically stable, and there exists a positive unstable equilibrium; when
, 0 is unstable and there exists a negative and locally asymptotically stable equilibrium.
2)
,
. When
with
, 0 is unstable; when
, 0 is locally asymptotically stable, and there exists a positive unstable equilibrium.
3)
,
. When
with
, 0 is unstable, and there exists a locally asymptotically stable negative equilibrium; when
, 0 is stable, and a positive unstable equilibrium appears.
4)
,
. When
changes from negative to positive, 0 changes its stability from stable to unstable. Correspondingly a negative unstable equilibrium becomes positive and locally asymptotically stable.
2.8.2. Bifurcation Coefficients
It is clear that
so, the consideration is on
and
.
Hence, the second order partial derivatives at disease-free equilibrium and at
are given by:
,
,
,
,
,
,
,
.
Then,
Hence,
.
Also, for the sign of b given as:
(18)
Therefore,
and
, hence the model system exhibit a forward bifurcation.
2.8.3. Global Stability of Endemic Equilibrium Point
To prove the global stability of EEP, we used quadratic Lyapunov function. The endemic equilibrium point is said to be globally asymptotically stable if its derivative is less than 0.
Theorem 6. The FMD has a unique endemic equilibrium point
for model system (1) that is globally asymptotically stable if
otherwise unstable.
Proof. Now, consider the quadratic Lyapunov function,
where
is the cattle population and
is the endemic equilibrium point. Then,
(19)
Clearly,
is a continuous and differentiable function. We have:
But,
(20)
And
(21)
Substitute (20) and (21) into
gives:
That is,
(22)
Therefore, it is clear that
which implies that the Endemic Equilibrium Point (
) is globally asymptotically stable.¨
2.9. Numerical Simulations of the Model
This section uses accessible parameter values from the literature as well as estimated ones to simulate the model system (1) numerically. Unless otherwise specified, simulation results will be based on parameter values from Table 2. The fourth-order Runge-Kutta method (RK4) is used to simulate the model system (1) in MATLAB software. This combination ensures accurate representation of the FMD transmission model dynamics. Our simulation aims to fill knowledge gaps about Foot and Mouth Disease and illustrate its spread in a population with active immigrants. The chosen initial conditions are:
,
,
,
, and
. These values are deliberately chosen to produce a desired model behavior.
Figure 2 describes the changes in animal population over time. The trend in the graph implies that active immigrants play a crucial role in shaping the dynamics of an FMD outbreak. The rapid decrease in susceptibility, combined with gradual changes in other compartments, suggests a limited disease establishment potential. This emphasizes the need for prompt interventions and highlights the complex interplay of disease spread within the context of active immigrants.
2.10. Results and Discussion
Our results can be summarized as follows: Figures 3-7, Figure 9 and Figure 10 demonstrate that, with an increase in infectious immigrants (
and
); the susceptible animal population gradually decreases to a steady state. This suggests
Figure 2. Evolution of animal population with time.
Figure 3. Effect of α1 on susceptible animal population.
Figure 4. Effect of α1 on exposed animal population.
a connection between immigrant influx and susceptible population decline. Infectious immigrants notably contribute to reducing susceptibility through disease transmission. Subsequent gradual decrease to its steady state shows ongoing disease dynamics and the establishment of a stable endemic equilibrium.
Figure 3 shows an initial decline within 10 days followed by a gradual rise with
Figure 5. Effect of α1 on asymptomatic animal population.
Figure 6. Effect of α1 on symptomatic animal population.
an increase in (
); this might reflect complex exposure-immigrant dynamics. The fluctuations possibly reveal interactions among infection introduction, host responses, and adaptation processes. Figure 7 is characterized by decreasing trend, followed by a sharp increase and gradual decrease before reaching the steady state, this pattern is shaped by the arrival of infectious immigrants (
), showcasing the dynamic interplay between disease introduction, population response, and
Figure 7. Effect of α2 on susceptible animal population.
Figure 8. Effect of α2 on exposed animal population.
Adaptation (see Figure 8). The steady state signifies ongoing disease dynamics and the formation of a stable endemic equilibrium.
In Figure 5, rising active immigrants (
) lead to initial increase, then gradual decrease due to possibly virulent immigrant strains causing higher symptomatic cases and reduce asymptomatic individuals, or new variants by infective immigrants. Figure 9 displays sharp increase followed by a gradual decrease towards a stable
Figure 9. Effect of α2 on asymptomatic animal population.
Figure 10. Effect of α2 on symptomatic animal population.
state with an increase in (
); this could be attributed to the fact that infectious immigrants can rapidly spread disease in local animals and shift animals to symptomatic states, reducing asymptomatic numbers. Some animals gain immunity, increasing asymptomatic or resistant individuals. Equilibrium shows disease-immune balance, highlighting dynamic adaptation.
Figure 10 portrays a sharp increase and gradual decrease to stable condition with an increase in (
); possible factors for this could include, a phase of local adaptation to the new agent, during this period, recovery could increase as the agent spreads, also, more animals acquire immunity through infections and recovery. The rapid rise in recoveries signifies initial immunity development to the new agent. Subsequent gradual increase and steady state show ongoing disease dynamics and the establishment of a stable endemic equilibrium.
The key observation is that the presence of infectious immigrants (
and
) is closely associated with increased susceptibility in animal populations. This implies that the local animal population becomes more susceptible to FMD when exposed to these infectious immigrants. This result underscores the complexity of disease transmission dynamics and has important implications for designing targeted interventions and management strategies. The findings provide valuable insights that contribute to a more comprehensive understanding of the factors influencing the establishment and spread of infectious diseases in populations.
3. Conclusions
In summary, the figures collectively highlight the intricate interplay of disease dynamics in animal populations, driven significantly by the influx of infectious immigrants (
and
). These immigrants significantly shape susceptibility and recovery populations, causing fluctuations and equilibrium. The observed patterns stress the need to account for immigrant influence in understanding disease transmission, population responses, and stability. These insights have implications for disease control and management strategies. In essence, the figures underscore the dynamic nature of disease processes and the significance of grasping the interaction between immigrant influx, host responses, and disease outcomes for effective management.
While these figures provide valuable insights, there are areas requiring further exploration. Specifically, researching the mechanisms of interactions between immigrants and local populations, understanding strain impacts and transmission dynamics, and investigating genetic, ecological, and epidemiological factors will refine our comprehension. Advancing understanding in these aspects will enhance precise disease management strategies. The study’s findings provide a foundation for shaping targeted and adaptive intervention strategies in the management of Foot and Mouth Disease. These recommendations emphasize the need for a multifaceted approach, combining early detection, preventive measures, adaptive management, and international collaboration to effectively address the complexities introduced by active immigrants in disease transmission dynamics. For the future, the study can be extended in a multidisciplinary approach, incorporating insights from epidemiology, immunology, economics, and genomics, which can provide a holistic understanding of the complexities associated with FMD and contribute to the global effort to control and mitigate the impact of this important livestock disease.
Acknowledgements
The authors would like to extend their appreciation to the supervisors for their valuable input and constructive suggestions that significantly enhanced the quality of this work.