ABSTRACT
In this paper, we introduce the notion of practically susceptible population, which is a fraction of the biologically susceptible population. Assuming that the fraction depends on the severity of the epidemic and the public's level of precaution (as a response of the public to the epidemic), we propose a general framework model with the response level evolving with the epidemic. We firstly verify the well-posedness and confirm the disease's eventual vanishing for the framework model under the assumption that the basic reproduction number R 0 < 1 . For R 0 > 1 , we study how the behavioural response evolves with epidemics and how such an evolution impacts the disease dynamics. More specifically, when the precaution level is taken to be the instantaneous best response function in literature, we show that the endemic dynamic is convergence to the endemic equilibrium; while when the precaution level is the delayed best response, the endemic dynamic can be either convergence to the endemic equilibrium, or convergence to a positive periodic solution. Our derivation offers a justification/explanation for the best response used in some literature. By replacing "adopting the best response" with "adapting toward the best response", we also explore the adaptive long-term dynamics.
Subject(s)
Basic Reproduction Number , Communicable Diseases , Epidemics , Mathematical Concepts , Models, Biological , Humans , Basic Reproduction Number/statistics & numerical data , Epidemics/statistics & numerical data , Epidemics/prevention & control , Communicable Diseases/epidemiology , Communicable Diseases/transmission , Disease Susceptibility/epidemiology , Epidemiological Models , Biological Evolution , Computer SimulationABSTRACT
This paper examines the short-term or transient dynamics of SIR infectious disease models in patch environments. We employ reactivity of an equilibrium and amplification rates, concepts from ecology, to analyze how dispersals/travels between patches, spatial heterogeneity, and other disease-related parameters impact short-term dynamics. Our findings reveal that in certain scenarios, due to the impact of spatial heterogeneity and the dispersals, the short-term disease dynamics over a patch environment may disagree with the long-term disease dynamics that is typically reflected by the basic reproduction number. Such an inconsistence can mislead the public, public healthy agencies and governments when making public health policy and decisions, and hence, these findings are of practical importance.
Subject(s)
Communicable Diseases , Epidemiological Models , Humans , Models, Biological , Mathematical Concepts , Communicable Diseases/epidemiology , Ecology , Basic Reproduction Number , Population DynamicsABSTRACT
In this paper, a stoichiometric aquatic tri-trophic level model is proposed and analyzed, which incorporates the effect of light and phosphorus, as well as the fear effect in predator-prey interactions. The analysis of the model includes the dissipativity and the existence and stability of equilibria. The influence of environmental factors and fear effect on the dynamics of the system is particularly investigated. The key findings reveal that the coexistence of populations is positively influenced by an appropriate level of light intensity and/or the dissolved phosphorus input concentration; however, excessive levels of phosphorus input can disrupt the system, leading to chaotic behaviors. Furthermore, it is found that the fear effect can stabilize the system and promote the chances of population coexistence.
Subject(s)
Ecosystem , Models, Biological , Animals , Predatory Behavior , Population Dynamics , Fear , Phosphorus , Food ChainABSTRACT
We develop and analyze a mathematical model of oncolytic virotherapy in the treatment of melanoma. We begin with a special, local case of the model, in which we consider the dynamics of the tumour cells in the presence of an oncolytic virus at the primary tumour site. We then consider the more general regional model, in which we incorporate a linear network of lymph nodes through which the tumour cells and the oncolytic virus may spread. The modelling also considers the impact of hypoxia on the disease dynamics. The modelling takes into account both the effects of hypoxia on tumour growth and spreading, as well as the impact of hypoxia on oncolytic virotherapy as a treatment modality. We find that oxygen-rich environments are favourable for the use of adenoviruses as oncolytic agents, potentially suggesting the use of complementary external oxygenation as a key aspect of treatment. Furthermore, the delicate balance between a virus' infection capabilities and its oncolytic capabilities should be considered when engineering an oncolytic virus. If the virus is too potent at killing tumour cells while not being sufficiently effective at infecting them, the infected tumour cells are destroyed faster than they are able to infect additional tumour cells, leading less favourable clinical results. Numerical simulations are performed in order to support the analytic results and to further investigate the impact of various parameters on the outcomes of treatment. Our modelling provides further evidence indicating the importance of three key factors in treatment outcomes: tumour microenvironment oxygen concentration, viral infection rates, and viral oncolysis rates. The numerical results also provide some estimates on these key model parameters which may be useful in the engineering of oncolytic adenoviruses.
ABSTRACT
In this paper, we propose a delay differential equation model to describe the Wolbachia infection dynamics in mosquitoes in which the key factor of cytoplasmic incompactibility (CI) is incorporated in a more natural way than those in the literature. By analyzing the dynamics of the model, we are able to obtain some information on the impact of four important parameters: the competition capabilities of the wild mosquitoes and infected mosquitoes, the maternal transmission level and the CI level. The analytic results show that there are ranges of parameters that support competition exclusion principle, and there are also ranges of parameters that allow co-persistence for both wild and infected mosquitoes. These ranges account for the scenarios of failure of invasion, invasion and suppressing the wild mosquitoes, and invasion and replacing the wild mosquitoes. We also discuss some possible future problems both in mathematics and in modeling.
Subject(s)
Aedes , Wolbachia , Animals , Mathematical Concepts , Models, BiologicalABSTRACT
In this paper, we revisit the notion of infection force from a new angle which can offer a new perspective to motivate and justify some infection force functions. Our approach can not only explain many existing infection force functions in the literature, it can also motivate new forms of infection force functions, particularly infection forces depending on disease surveillance of the past. As a demonstration, we propose an SIRS model with delay. We comprehensively investigate the disease dynamics represented by this model, particularly focusing on the local bifurcation caused by the delay and another parameter that reflects the weight of the past epidemics in the infection force. We confirm Hopf bifurcations both theoretically and numerically. The results show that, depending on how recent the disease surveillance data are, their assigned weight may have a different impact on disease control measures.
Subject(s)
Communicable Diseases , Epidemics , Humans , Communicable Diseases/epidemiologyABSTRACT
In this paper, based on the classic Kermack-McKendrick SIR model, we propose an ordinary differential equation model to re-examine the COVID-19 epidemics in Wuhan where this disease initially broke out. The focus is on the impact of all those major non-pharmaceutical interventions (NPIs) implemented by the local public healthy authorities and government during the epidemics. We use the data publicly available and the nonlinear least-squares solver lsqnonlin built in MATLAB to estimate the model parameters. Then we explore the impact of those NPIs, particularly the timings of these interventions, on the epidemics. The results can help people review the responses to the outbreak of the COVID-19 in Wuhan, while the proposed model also offers a framework for studying epidemics of COVID-19 and/or other similar diseases in other places, and accordingly helping people better prepare for possible future outbreaks of similar diseases.
ABSTRACT
When perceiving a risk from predators, a prey may respond by reducing its reproduction and decreasing or increasing (depending on the species) its mobility. We formulate a patch model to investigate the aforementioned fear effect which is indirect, in contrast to the predation as a direct effect, of the predator on the prey population. We consider not only cost but also benefit of anti-predation response of the prey, and explore their trade-offs together as well as the impact of the fear effect mediated dispersals of the prey. In the case of constant response level, if there is no dispersal and for some given response functions, the model indicates the existence of an evolutionary stable strategy which is also a convergence stable strategy for the response level; and if there is dispersal, the analysis of the model shows that it will enhance the co-persistence of the prey on both patches. Considering the trait as another variable, we continue to study the evolution of anti-predation strategy for the model with dispersal, which leads to a three-dimensional system of ordinary differential equations. We perform some numerical simulations, which demonstrate global convergence to a positive equilibrium with the response level evolving towards a positive constant level, implying the existence of an optimal anti-predation response level.
Subject(s)
Adaptation, Physiological , Biological Evolution , Models, Biological , Predatory Behavior , Animals , Computer Simulation , Food Chain , Population DynamicsABSTRACT
Recent experimental study suggests that the engineered symbiotic bacteria Serratia AS1 may provide a novel, effective and sustainable biocontrol of malaria. These recombinant bacteria have been shown to be able to rapidly disseminate throughout mosquito populations and to efficiently inhibit development of malaria parasites in mosquitoes in controlled laboratory experiments. In this paper, we develop a climate-based malaria model which involves both vertical and horizontal transmissions of the engineered Serratia AS1 bacteria in mosquito population. We show that the dynamics of the model system is totally determined by the vector reproduction ratio [Formula: see text], and the basic reproduction ratio [Formula: see text]. If [Formula: see text], then the mosquito-free state is globally attractive. If [Formula: see text] and [Formula: see text], then the disease-free periodic solution is globally attractive. If [Formula: see text] and [Formula: see text], then the positive periodic solution is globally attractive. Numerically, we verify the obtained analytic result and evaluate the effects of releasing the engineered Serratia AS1 bacteria in field by conducting a case study for Douala, Cameroon. We find that ideally, by using Serratia AS1 alone, it takes at least 25 years to eliminate malaria from Douala. This implies that continued long-term investment is needed in the fight against malaria and confirms the necessity of integrating multiple control measures.
Subject(s)
Biological Control Agents , Malaria/prevention & control , Models, Biological , Serratia/physiology , Symbiosis/physiology , Animals , Anopheles/microbiology , Anopheles/parasitology , Bioengineering , Cameroon , Female , Humans , Malaria/transmission , Male , Mathematical Concepts , Mosquito Vectors/microbiology , Mosquito Vectors/parasitology , Serratia/geneticsABSTRACT
yme disease is transmitted via blacklegged ticks, the spatial spread of which is believed to be primarily via transport on white-tailed deer. In this paper, we develop a mathematical model to describe the spatial spread of blacklegged ticks due to deer dispersal. The model turns out to be a system of differential equations with a spatially non-local term accounting for the phenomenon that a questing female adult tick that attaches to a deer at one location may later drop to the ground, fully fed, at another location having been transported by the deer. We first justify the well-posedness of the model and analyze the stability of its steady states. We then explore the existence of traveling wave fronts connecting the extinction equilibrium with the positive equilibrium for the system. We derive an algebraic equation that determines a critical value c* which is at least a lower bound for the wave speed in the sense that, if cSubject(s)
Arachnid Vectors
, Deer/parasitology
, Ixodes
, Models, Biological
, Tick Infestations/veterinary
, Animal Distribution
, Animals
, Arachnid Vectors/microbiology
, Computer Simulation
, Geography
, Humans
, Ixodes/microbiology
, Lyme Disease/transmission
, Mathematical Concepts
, Population Growth
, Tick Infestations/parasitology
, United States
ABSTRACT
We propose and analyse a reaction-diffusion-advection predator-prey model in which we assume that predators move randomly but prey avoid predation by perceiving a repulsion along predator density gradient. Based on recent experimental evidence that anti-predator behaviors alone lead to a 40% reduction on prey reproduction rate, we also incorporate the cost of anti-predator responses into the local reaction terms in the model. Sufficient and necessary conditions of spatial pattern formation are obtained for various functional responses between prey and predators. By mathematical and numerical analyses, we find that small prey sensitivity to predation risk may lead to pattern formation if the Holling type â ¡ functional response or the Beddington-DeAngelis functional response is adopted while large cost of anti-predator behaviors homogenises the system by excluding pattern formation. However, the ratio-dependent functional response gives an opposite result where large predator-taxis may lead to pattern formation but small cost of anti-predator behaviors inhibits the emergence of spatial heterogeneous solutions.
Subject(s)
Models, Biological , Predatory Behavior , Animals , Computer Simulation , Ecosystem , Food Chain , Mathematical Concepts , Population DynamicsABSTRACT
Mosquito-borne diseases remain a significant threat to public health and economics. Since mosquitoes are quite sensitive to temperature, global warming may not only worsen the disease transmission case in current endemic areas but also facilitate mosquito population together with pathogens to establish in new regions. Therefore, understanding mosquito population dynamics under the impact of temperature is considerably important for making disease control policies. In this paper, we develop a stage-structured mosquito population model in the environment of a temperature-controlled experiment. The model turns out to be a system of periodic delay differential equations with periodic delays. We show that the basic reproduction number is a threshold parameter which determines whether the mosquito population goes to extinction or remains persistent. We then estimate the parameter values for Aedes aegypti, the mosquito that transmits dengue virus. We verify the analytic result by numerical simulations with the temperature data of Colombo, Sri Lanka where a dengue outbreak occurred in 2017.
Subject(s)
Models, Biological , Mosquito Vectors , Aedes/growth & development , Aedes/virology , Animals , Basic Reproduction Number , Computer Simulation , Dengue/epidemiology , Dengue/transmission , Dengue Virus/pathogenicity , Female , Global Warming , Humans , Mathematical Concepts , Mosquito Vectors/growth & development , Mosquito Vectors/virology , Population Dynamics/statistics & numerical data , TemperatureABSTRACT
In this paper, a very general model of impulsive delay differential equations in n-patches is rigorously derived to describe the impulsive control of population of a single species over n-patches. The model allows an age structure consisting of immatures and matures, and also considers mobility and culling of both matures and immatures. Conditions are obtained for extinction and persistence of the model system under three special scenarios: (1) without impulsive control; (2) with impulsive culling of the immatures only; and (3) with impulsive culling of the matures only, respectively. In the case of persistence, the persistence level is also estimated for the systems in the case of identical n patches, by relating the issue to the dynamics of multi-dimensional maps. Two illustrative examples and their numerical simulations are given to show the effectiveness of the results. Based on the theoretical results, some strategies of impulsive culling are provided to eradicate the population of a pest species.
Subject(s)
Models, Biological , Population Dynamics/statistics & numerical data , Age Factors , Animal Culling/statistics & numerical data , Animals , Computational Biology , Computer Simulation , Ecosystem , Extinction, Biological , Humans , Infection Control/statistics & numerical data , Mathematical Concepts , Pest Control/statistics & numerical data , Pest Control, Biological/statistics & numerical data , Stochastic ProcessesABSTRACT
Recent field experiments on vertebrates showed that the mere presence of a predator would cause a dramatic change of prey demography. Fear of predators increases the survival probability of prey, but leads to a cost of prey reproduction. Based on the experimental findings, we propose a predator-prey model with the cost of fear and adaptive avoidance of predators. Mathematical analyses show that the fear effect can interplay with maturation delay between juvenile prey and adult prey in determining the long-term population dynamics. A positive equilibrium may lose stability with an intermediate value of delay and regain stability if the delay is large. Numerical simulations show that both strong adaptation of adult prey and the large cost of fear have destabilizing effect while large population of predators has a stabilizing effect on the predator-prey interactions. Numerical simulations also imply that adult prey demonstrates stronger anti-predator behaviors if the population of predators is larger and shows weaker anti-predator behaviors if the cost of fear is larger.
Subject(s)
Fear , Models, Biological , Predatory Behavior , Adaptation, Physiological , Adolescent , Animals , Behavior, Animal , Demography , Humans , Population DynamicsABSTRACT
We propose a mathematical model for biocontrol of the invasive weed Fallopia japonica using one of its co-evolved natural enemies, the Japanese sap-sucking psyllid Aphalara itadori. This insect sucks the sap from the stems of the plant thereby weakening it. Its diet is highly specific to F. japonica. We consider a single isolated knotweed stand, the plant's size being described by time-dependent variables for total stem and rhizome biomass. It is the larvae of A. itadori that damage the plant most, so the insect population is described in terms of variables for the numbers of larvae and adults, using a stage-structured modelling approach. The dynamics of the model depends mainly on a parameter h, which measures how long it takes for an insect to handle (digest) one unit of F. japonica stem biomass. If h is too large, then the model does not have a positive equilibrium and the plant biomass and insect numbers both grow together without bound, though at a lower rate than if the insects were absent. If h is sufficiently small, then the model possesses a positive equilibrium which appears to be locally stable. The results based on our model imply that satisfactory long-term control of the knotweed F. japonica using the insect A. itadori is only possible if the insect is able to consume and digest knotweed biomass sufficiently quickly; if it cannot, then the insect can only slow down the growth which is still unbounded.
Subject(s)
Biological Control Agents , Fallopia japonica , Plant Weeds , Animals , Biomass , Fallopia japonica/growth & development , Fallopia japonica/parasitology , Hemiptera/physiology , Mathematical Concepts , Models, Biological , Plant Weeds/growth & development , Plant Weeds/parasitologyABSTRACT
While studies on pre-exposure prophylaxis (PrEP) and post-exposure prophylaxis (PEP) have demonstrated substantial advantages in controlling HIV transmission, the overall benefits of the programs with early initiation of antiretroviral therapy (ART) have not been fully understood and are still on debate. Here, we develop an immunity-based (CD4+ T cell count based) mathematical model to study the impacts of early treatment programs on HIV epidemics and the overall community-level immunity. The model is parametrized using the HIV prevalence data from South Africa and fully analyzed for stability of equilibria and infection persistence criteria. Using our model, we evaluate the effects of early treatment on the new infection transmission, disease death, basic reproduction number, HIV prevalence, and the community-level immunity. Our model predicts that the programs with early treatments significantly reduce the new infection transmission and increase the community-level immunity, but the treatments alone may not be enough to eliminate HIV epidemics. These findings, including the community-level immunity, might provide helpful information for proper implementation of HIV treatment programs.
Subject(s)
Antirheumatic Agents/therapeutic use , CD4-Positive T-Lymphocytes , Epidemics/statistics & numerical data , HIV Infections , Models, Theoretical , HIV Infections/epidemiology , HIV Infections/immunology , HIV Infections/therapy , HumansABSTRACT
A recent field manipulation on a terrestrial vertebrate showed that the fear of predators alone altered anti-predator defences to such an extent that it greatly reduced the reproduction of prey. Because fear can evidently affect the populations of terrestrial vertebrates, we proposed a predator-prey model incorporating the cost of fear into prey reproduction. Our mathematical analyses show that high levels of fear (or equivalently strong anti-predator responses) can stabilize the predator-prey system by excluding the existence of periodic solutions. However, relatively low levels of fear can induce multiple limit cycles via subcritical Hopf bifurcations, leading to a bi-stability phenomenon. Compared to classic predator-prey models which ignore the cost of fear where Hopf bifurcations are typically supercritical, Hopf bifurcations in our model can be both supercritical and subcritical by choosing different sets of parameters. We conducted numerical simulations to explore the relationships between fear effects and other biologically related parameters (e.g. birth/death rate of adult prey), which further demonstrate the impact that fear can have in predator-prey interactions. For example, we found that under the conditions of a Hopf bifurcation, an increase in the level of fear may alter the direction of Hopf bifurcation from supercritical to subcritical when the birth rate of prey increases accordingly. Our simulations also show that the prey is less sensitive in perceiving predation risk with increasing birth rate of prey or increasing death rate of predators, but demonstrate that animals will mount stronger anti-predator defences as the attack rate of predators increases.
Subject(s)
Behavior, Animal/physiology , Fear , Models, Biological , Predatory Behavior/physiology , Animals , Escape Reaction/physiologyABSTRACT
Many observational studies suggest that seasonal migratory birds play an important role in spreading Ixodes scapularis, a vector of Lyme disease, along their migratory flyways, and they are believed to be responsible for geographic range expansion of I. scapularis in Canada. However, the interplay between the dynamics of I. scapularis on land and migratory birds in the air is not well understood. In this study, we develop a periodic delay meta-population model which takes into consideration the local landscape for tick reproduction within patches and the times needed for ticks to be transported by birds between patches. Assuming that the tick population is endemic in the source region, we find that bird migration may boost an already established tick population at the subsequent region and thus increase the risk to humans, or bird migration may help ticks to establish in a region where the local landscape is not appropriate for ticks to survive in the absence of bird migration, imposing risks to public health. This theoretical study reveals that bird migration plays an important role in the geographic range expansion of I. scapularis, and therefore our findings may suggest some strategies for Lyme disease prevention and control.
Subject(s)
Animal Migration , Birds/parasitology , Ixodes/pathogenicity , Animals , Arachnid Vectors , Canada , Computer Simulation , Humans , Ixodes/microbiology , Lyme Disease/transmission , Mathematical Concepts , Models, Biological , Population Dynamics , SeasonsABSTRACT
This paper investigates consequences of vaccine implementation strategies for infectious diseases by a mathematical model. For an infectious disease, the degree of infection may vary widely among the individuals. Reports show that individuals belonging to certain groups possess considerably higher risk to infection. Incorporating this phenomenon into vaccination strategies, the host is categorized into different groups to measure the outcome of the vaccination. A mathematical model is proposed and analysed to evaluate this measure. Our results suggest that vaccinating a group with certain priority may lead to elimination of the disease effectively. The strategy is cost-effective as well.
Subject(s)
Communicable Diseases/epidemiology , Models, Biological , Vaccination , Basic Reproduction Number , Communicable Diseases/immunology , HumansABSTRACT
In this paper, we consider the evolutionary competition between budding and lytic viral release strategies, using a delay differential equation model with distributed delay. When antibody is not established, the dynamics of competition depends on the respective basic reproductive ratios of the two viruses. If the basic reproductive ratio of budding virus is greater than that of lytic virus and one, budding virus can survive. When antibody is established for both strains but the neutralization capacities are the same for both strains, consequence of the competition also depends only on the basic reproductive ratios of the budding and lytic viruses. Using two concrete forms of the viral production functions, we are also able to conclude that budding virus will outcompete if the rates of viral production, death rates of infected cells and neutralizing capacities of the antibodies are the same for budding and lytic viruses. In this case, budding strategy would have an evolutionary advantage. However, if the antibody neutralization capacity for the budding virus is larger than that for the lytic virus, the lytic virus can outcompete the budding virus provided that its reproductive ratio is very high. An explicit threshold is derived.