Threshold dynamics of a reaction-advection-diffusion schistosomiasis epidemic model with seasonality and spatial heterogeneity.

Most water-borne disease models ignore the advection of water flows in order to simplify the mathematical analysis and numerical computation. However, advection can play an important role in determining the disease transmission dynamics. In this paper, we investigate the long-term dynamics of a periodic reaction-advection-diffusion schistosomiasis model and explore the joint impact of advection, seasonality and spatial heterogeneity on the transmission of the disease. We derive the basic reproduction number R 0 and show that the disease-free periodic solution is globally attractive when R 0 < 1 whereas there is a positive endemic periodic solution and the system is uniformly persistent in a special case when R 0 > 1 . Moreover, we find that R 0 is a decreasing function of the advection coefficients which offers insights into why schistosomiasis is more serious in regions with slow water flows.

Modelling phytoplankton-virus interactions: phytoplankton blooms and lytic virus transmission.

A dynamic reaction-diffusion model of four variables is proposed to describe the spread of lytic viruses among phytoplankton in a poorly mixed aquatic environment. The basic ecological reproductive index for phytoplankton invasion and the basic reproduction number for virus transmission are derived to characterize the phytoplankton growth and virus transmission dynamics. The theoretical and numerical results from the model show that the spread of lytic viruses effectively controls phytoplankton blooms. This validates the observations and experimental results of Emiliana huxleyi-lytic virus interactions. The studies also indicate that the lytic virus transmission cannot occur in a low-light or oligotrophic aquatic environment.

Markovian Approach for Exploring Competitive Diseases with Heterogeneity-Evidence from COVID-19 and Influenza in China.

Due to the complex interactions between multiple infectious diseases, the spreading of diseases in human bodies can vary when people are exposed to multiple sources of infection at the same time. Typically, there is heterogeneity in individuals' responses to diseases, and the transmission routes of different diseases also vary. Therefore, this paper proposes an SIS disease spreading model with individual heterogeneity and transmission route heterogeneity under the simultaneous action of two competitive infectious diseases. We derive the theoretical epidemic spreading threshold using quenched mean-field theory and perform numerical analysis under the Markovian method. Numerical results confirm the reliability of the theoretical threshold and show the inhibitory effect of the proportion of fully competitive individuals on epidemic spreading. The results also show that the diversity of disease transmission routes promotes disease spreading, and this effect gradually weakens when the epidemic spreading rate is high enough. Finally, we find a negative correlation between the theoretical spreading threshold and the average degree of the network. We demonstrate the practical application of the model by comparing simulation outputs to temporal trends of two competitive infectious diseases, COVID-19 and seasonal influenza in China.

Modeling the effect of imported malaria on the elimination programme in KwaZulu-Natal province of South Africa.

Introduction: with imported malaria cases in a given population, the question arises as to what extent the local cases are a consequence of the imports or not. We perform a modeling analysis for a specific area, in a region aspiring for malaria-free status. Methods: data on malaria cases over ten years is subjected to a compartmental model which is assumed to be operating close to the equilibrium state. Two of the parameters of the model are fitted to the decadal data. The other parameters in the model are sourced from the literature. The model is utilized to simulate the malaria prevalence with or without imported cases. Results: in any given year the annual average of 460 imported cases, resulted in an end-of-year season malaria prevalence of 257 local active infectious cases, whereas without the imports the malaria prevalence at the end of the season would have been fewer than 10 active infectious cases. We calculate the numerical value of the basic reproduction number for the model, which reveals the extent to which the disease is being eliminated from the population or not. Conclusion: without the imported cases, over the ten seasons of malaria, 2008-2018, the KwaZulu-Natal province would have been malaria-free over at least the last 7 years of the decade indicated. This simple methodology works well even in situations where data is limited.

Modelling the impact of precaution on disease dynamics and its evolution.

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.

Linking Spontaneous Behavioral Changes to Disease Transmission Dynamics: Behavior Change Includes Periodic Oscillation.

Behavior change significantly influences the transmission of diseases during outbreaks. To incorporate spontaneous preventive measures, we propose a model that integrates behavior change with disease transmission. The model represents behavior change through an imitation process, wherein players exclusively adopt the behavior associated with higher payoff. We find that relying solely on spontaneous behavior change is insufficient for eradicating the disease. The dynamics of behavior change are contingent on the basic reproduction number R a corresponding to the scenario where all players adopt non-pharmaceutical interventions (NPIs). When R a < 1 , partial adherence to NPIs remains consistently feasible. We can ensure that the disease stays at a low level or maintains minor fluctuations around a lower value by increasing sensitivity to perceived infection. In cases where oscillations occur, a further reduction in the maximum prevalence of infection over a cycle can be achieved by increasing the rate of behavior change. When R a > 1 , almost all players consistently adopt NPIs if they are highly sensitive to perceived infection. Further consideration of saturated recovery leads to saddle-node homoclinic and Bogdanov-Takens bifurcations, emphasizing the adverse impact of limited medical resources on controlling the scale of infection. Finally, we parameterize our model with COVID-19 data and Tokyo subway ridership, enabling us to illustrate the disease spread co-evolving with behavior change dynamics. We further demonstrate that an increase in sensitivity to perceived infection can accelerate the peak time and reduce the peak size of infection prevalence in the initial wave.

Dynamical analysis of a general delayed HBV infection model with capsids and adaptive immune response in presence of exposed infected hepatocytes.

The aim of this paper is to develop and investigate a novel mathematical model of the dynamical behaviors of chronic hepatitis B virus infection. The model includes exposed infected hepatocytes, intracellular HBV DNA-containing capsids, uses a general incidence function for viral infection covering a variety of special cases available in the literature, and describes the interaction of cytotoxic T lymphocytes that kill the infected hepatocytes and the magnitude of B-cells that send antibody immune defense to neutralize free virions. Further, one time delay is incorporated to account for actual capsids production. The other time delays are used to account for maturation of capsids and free viruses. We start with the analysis of the proposed model by establishing the local and global existence, uniqueness, non-negativity and boundedness of solutions. After defined the threshold parameters, we discuss the stability properties of all possible steady state constants by using the crafty Lyapunov functionals, the LaSalle's invariance principle and linearization methods. The impacts of the three time delays on the HBV infection transmission are discussed through local and global sensitivity analysis of the basic reproduction number and of the classes of infected states. Finally, an application is provided and numerical simulations are performed to illustrate and interpret the theoretical results obtained. It is suggested that, a good strategy to eradicate or to control HBV infection within a host should concentrate on any drugs that may prolong the values of the three delays.

Evaluation of age-structured vaccination strategies for curbing the disease spread.

Age structure is one of the crucial factors in characterizing the heterogeneous epidemic transmission. Vaccination is regarded as an effective control measure for prevention and control epidemics. Due to the shortage of vaccine capacity during the outbreak of epidemics, how to design vaccination policy has become an urgent issue in suppressing the disease transmission. In this paper, we make an effort to propose an age-structured SVEIHR model with the disease-caused death to take account of dynamics of age-related vaccination policy for better understanding disease spread and control. We present an explicit expression of the basic reproduction number R 0 , which determines whether or not the disease persists, and then establish the existence and stability of endemic equilibria under certain conditions. Numerical simulations are illustrated to show that the age-related vaccination policy has a tremendous influence on curbing the disease transmission. Especially, vaccination of people over 65 is better than for people aged 21-65 in terms of rapid eradication of the disease in Italy.

Generative Bayesian modeling to nowcast the effective reproduction number from line list data with missing symptom onset dates.

The time-varying effective reproduction number Rt is a widely used indicator of transmission dynamics during infectious disease outbreaks. Timely estimates of Rt can be obtained from reported cases counted by their date of symptom onset, which is generally closer to the time of infection than the date of report. Case counts by date of symptom onset are typically obtained from line list data, however these data can have missing information and are subject to right truncation. Previous methods have addressed these problems independently by first imputing missing onset dates, then adjusting truncated case counts, and finally estimating the effective reproduction number. This stepwise approach makes it difficult to propagate uncertainty and can introduce subtle biases during real-time estimation due to the continued impact of assumptions made in previous steps. In this work, we integrate imputation, truncation adjustment, and Rt estimation into a single generative Bayesian model, allowing direct joint inference of case counts and Rt from line list data with missing symptom onset dates. We then use this framework to compare the performance of nowcasting approaches with different stepwise and generative components on synthetic line list data for multiple outbreak scenarios and across different epidemic phases. We find that under reporting delays realistic for hospitalization data (50% of reports delayed by more than a week), intermediate smoothing, as is common practice in stepwise approaches, can bias nowcasts of case counts and Rt, which is avoided in a joint generative approach due to shared regularization of all model components. On incomplete line list data, a fully generative approach enables the quantification of uncertainty due to missing onset dates without the need for an initial multiple imputation step. In a real-world comparison using hospitalization line list data from the COVID-19 pandemic in Switzerland, we observe the same qualitative differences between approaches. The generative modeling components developed in this work have been integrated and further extended in the R package epinowcast, providing a flexible and interpretable tool for real-time surveillance.

Global analysis of an age-structured tuberculosis model with an application to Jiangsu, China.

Diagnostic delay for TB infected individuals and the lack of TB vaccines for adults are the main challenges to achieve the goals of WHO by 2050. In order to evaluate the impacts of diagnostic delay and vaccination for adults on prevalence of TB, we propose an age-structured model with latent age and infection age, and we incorporate Mycobacterium TB in the environment and vaccination into the model. Diagnostic delay is indicated by the age of infection before receiving treatment. The threshold dynamics are established in terms of the basic reproduction number R 0 . When R 0 < 1 , the disease-free equilibrium is globally asymptotically stable, which means that TB epidemic will die out; When R 0 = 1 , the disease-free equilibrium is globally attractive; there exists a unique endemic equilibrium and the endemic equilibrium is globally attractive when R 0 > 1 . We estimate that the basic reproduction number R 0 = 0.5320 (95% CI (0.3060, 0.7556)) in Jiangsu Province, which means that TB epidemic will die out. However, we find that the annual number of new TB cases by 2050 is 1,151 (95%CI: (138, 8,014)), which means that it is challenging to achieve the goal of WHO by 2050. To this end, we evaluate the possibility of achieving the goals of WHO if we start vaccinating adults and reduce diagnostic delay in 2025. Our results demonstrate that when the diagnostic delay is reduced from longer than four months to four months, or 20% adults are vaccinated, the goal of WHO in 2050 can be achieved, and 73,137 (95%CI: (23,906, 234,086)) and 54,828 (95%CI: (15,811, 206,468)) individuals will be prevented from being infected from 2025 to 2050, respectively. The modeling approaches and simulation results used in this work can help policymakers design control measures to reduce the prevalence of TB.

Modeling outbreaks of COVID-19 in China: The impact of vaccination and other control measures on curbing the epidemic.

This study aims to examine the development trend of COVID-19 in China and propose a model to assess the impacts of various prevention and control measures in combating the COVID-19 pandemic. Using COVID-19 cases reported by the National Health Commission of China from January 2, 2020, to January 2, 2022, we established a Susceptible-Exposed-Infected-Asymptomatic-Quarantined-Vaccinated-Hospitalized-Removed (SEIAQVHR) model to calculate the COVID-19 transmission rate and Rt effective reproduction number, and assess prevention and control measures. Additionally, we built a stochastic model to explore the development of the COVID-19 epidemic. We modeled the incidence trends in five outbreaks between 2020 and 2022. Some important features of the COVID-19 epidemic are mirrored in the estimates based on our SEIAQVHR model. Our model indicates that an infected index case entering the community has a 50%-60% chance to cause a COVID-19 outbreak. Wearing masks and getting vaccinated were the most effective measures among all the prevention and control measures. Specifically targeting asymptomatic individuals had no significant impact on the spread of COVID-19. By adjusting prevention and control parameters, we suggest that increasing the rates of effective vaccination and mask-wearing can significantly reduce COVID-19 cases in China. Our stochastic model analysis provides a useful tool for understanding the COVID-19 epidemic in China.

Stability analysis and numerical evaluations of a COVID-19 model with vaccination.

A novel (nonlinear) mathematical model for the transmission of Coronavirus 19 (COVID-19) with eight compartments and considering the impact of vaccination is examined in this manuscript. The qualitative behavior of the system such as the boundedness of solutions, the basic reproduction number, and the stability of the equilibrium points is investigated in detail. Some domestic real data collected from the Kerman University of Medical Science (KUMC) is used to estimate the parameters of the proposed model. We predict the dynamical behavior of the system through numerical simulations based on a combined spectral matrix collocation methodology. In this respect, we first linearize the nonlinear system of equations by the method of quasilinearization (QLM). Hence, the shifted version of Chebyshev polynomials of the second kind (SCPSK) is utilized along with the domain-splitting strategy to acquire the solutions of the system over a long time interval. The uniform convergence and upper bound estimation of the SCPSK bases are proved in a rigorous manner. Moreover, the technique of residual error functions is used to testify the accuracy of the QLM-SCPSK method. The presented numerical results justify the robustness and good accuracy of the QLM-SCPSK technique. The achieved numerical orders of convergence indicate that the QLM-SCSK algorithm has exponential rate of convergence. Using the linearization technique in one hand and the domain-splitting strategy on the other hand, enable us to predict the behaviour of similar disease problems with high accuracy and maximum efficiency on an arbitrary domain of interest.

Combining wastewater surveillance and case data in estimating the time-varying effective reproduction number.

Wastewater surveillance has been increasingly acknowledged as a useful tool for monitoring transmission dynamics of infections of public health concern, including the coronavirus disease (COVID-19). While a range of models have been proposed to estimate the time-varying effective reproduction number (Rt) utilizing clinical data, few have harnessed the viral concentration in wastewater samples to do so, leaving uncertainties about the potential precision gains with its use. In this study, we developed a Bayesian hierarchical model which simultaneously reconstructed the latent infection trajectory and estimated Rt. Focusing on the 2022 and early 2023 COVID-19 transmission trends in Singapore, where mass community wastewater surveillance has become routine, we performed estimations using a spectrum of data sources, including reported case counts, hospital admissions, deaths, and wastewater viral loads. We further explored the performance of our wastewater model across various scenarios with different sampling strategies. The results showed consistent estimates derived from models employing diverse data streams, while models incorporating more wastewater samples exhibited greater uncertainty and variation in the inferred Rts. Additionally, our analysis revealed prominent day-of-the-week effect in reported case counts and substantial temporal variations in ascertainment rates. In response to these findings, we advocate for a hybrid approach leveraging both clinical and wastewater surveillance data to account for changes in case-ascertainment rates. Furthermore, our study demonstrates the possibility of reducing sampling frequency or sample size without compromising estimation accuracy for Rt, highlighting the potential for optimizing resource allocation in surveillance efforts while maintaining robust insights into the transmission dynamics of infectious diseases.

Assessing the Time Evolution of COVID-19 Effective Reproduction Number in Brazil.

In this paper, we use a Bayesian method to estimate the effective reproduction number ( R ( t ) ), in the context of monitoring the time evolution of the COVID-19 pandemic in Brazil at different geographic levels. The focus of this study is to investigate the similarities between the trends in the evolution of such indicators at different subnational levels with the trends observed nationally. The underlying question addressed is whether national surveillance of such variables is enough to provide a picture of the epidemic evolution in the country or if it may hide important localized trends. This is particularly relevant in the scenario where health authorities use information obtained from such indicators in the design of non-pharmaceutical intervention policies to control the epidemic. A comparison between R ( t ) estimates and the moving average (MA) of daily reported infections is also presented, which is another commonly monitored variable. The analysis carried out in this paper is based on the data of confirmed infected cases provided by a public repository. The correlations between the time series of R ( t ) and MA in different geographic levels are assessed. Comparing national with subnational trends, higher degrees of correlation are found for the time series of R ( t ) estimates, compared to the MA time series. Nevertheless, differences between national and subnational trends are observed for both indicators, suggesting that local epidemiological surveillance would be more suitable as an input to the design of non-pharmaceutical intervention policies in Brazil, particularly for the least populated states.

A novel simulation-based analysis of a stochastic HIV model with the time delay using high order spectral collocation technique.

The economic impact of Human Immunodeficiency Virus (HIV) goes beyond individual levels and it has a significant influence on communities and nations worldwide. Studying the transmission patterns in HIV dynamics is crucial for understanding the tracking behavior and informing policymakers about the possible control of this viral infection. Various approaches have been adopted to explore how the virus interacts with the immune system. Models involving differential equations with delays have become prevalent across various scientific and technical domains over the past few decades. In this study, we present a novel mathematical model comprising a system of delay differential equations to describe the dynamics of intramural HIV infection. The model characterizes three distinct cell sub-populations and the HIV virus. By incorporating time delay between the viral entry into target cells and the subsequent production of new virions, our model provides a comprehensive understanding of the infection process. Our study focuses on investigating the stability of two crucial equilibrium states the infection-free and endemic equilibriums. To analyze the infection-free equilibrium, we utilize the LaSalle invariance principle. Further, we prove that if reproduction is less than unity, the disease free equilibrium is locally and globally asymptotically stable. To ensure numerical accuracy and preservation of essential properties from the continuous mathematical model, we use a spectral scheme having a higher-order accuracy. This scheme effectively captures the underlying dynamics and enables efficient numerical simulations.

Fractional epidemic model of coronavirus disease with vaccination and crowding effects.

Most of the countries in the world are affected by the coronavirus epidemic that put people in danger, with many infected cases and deaths. The crowding factor plays a significant role in the transmission of coronavirus disease. On the other hand, the vaccines of the covid-19 played a decisive role in the control of coronavirus infection. In this paper, a fractional order epidemic model (SIVR) of coronavirus disease is proposed by considering the effects of crowding and vaccination because the transmission of this infection is highly influenced by these two factors. The nonlinear incidence rate with the inclusion of these effects is a better approach to understand and analyse the dynamics of the model. The positivity and boundedness of the fractional order model is ensured by applying some standard results of Mittag Leffler function and Laplace transformation. The equilibrium points are described analytically. The existence and uniqueness of the non-integer order model is also confirmed by using results of the fixed-point theory. Stability analysis is carried out for the system at both the steady states by using Jacobian matrix theory, Routh-Hurwitz criterion and Volterra-type Lyapunov functions. Basic reproductive number is calculated by using next generation matrix. It is verified that disease-free equilibrium is locally asymptotically stable if R 0 < 1 and endemic equilibrium is locally asymptotically stable if R 0 > 1 . Moreover, the disease-free equilibrium is globally asymptotically stable if R 0 < 1 and endemic equilibrium is globally asymptotically stable if R 0 > 1 . The non-standard finite difference (NSFD) scheme is developed to approximate the solutions of the system. The simulated graphs are presented to show the key features of the NSFD approach. It is proved that non-standard finite difference approach preserves the positivity and boundedness properties of model. The simulated graphs show that the implementation of control strategies reduced the infected population and increase the recovered population. The impact of fractional order parameter α is described by the graphical templates. The future trends of the virus transmission are predicted under some control measures. The current work will be a value addition in the literature. The article is closed by some useful concluding remarks.

Transmission dynamics of a reaction-advection-diffusion dengue fever model with seasonal developmental durations and intrinsic incubation periods.

In this paper, we propose a reaction-advection-diffusion dengue fever model with seasonal developmental durations and intrinsic incubation periods. Firstly, we establish the well-posedness of the model. Secondly, we define the basic reproduction number ℜ 0 for this model and show that ℜ 0 is a threshold parameter: if ℜ 0 < 1 , then the disease-free periodic solution is globally attractive; if ℜ 0 > 1 , the system is uniformly persistent. Thirdly, we study the global attractivity of the positive steady state when the spatial environment is homogeneous and the advection of mosquitoes is ignored. As an example, we use the model to investigate the dengue fever transmission case in Guangdong Province, China, and explore the impact of model parameters on ℜ 0 . Our findings indicate that ignoring seasonality may underestimate ℜ 0 . Additionally, the spatial heterogeneity of transmission may increase the risk of disease transmission, while the increase of seasonal developmental durations, intrinsic incubation periods and advection rates can all reduce the risk of disease transmission.

Understanding HIV/AIDS dynamics: insights from CD4+T cells, antiretroviral treatment, and country-specific analysis.

In this article, we present a mathematical model for human immunodeficiency virus (HIV)/Acquired immune deficiency syndrome (AIDS), taking into account the number of CD4+T cells and antiretroviral treatment. This model is developed based on the susceptible, infected, treated, AIDS (SITA) framework, wherein the infected and treated compartments are divided based on the number of CD4+T cells. Additionally, we consider the possibility of treatment failure, which can exacerbate the condition of the treated individual. Initially, we analyze a simplified HIV/AIDS model without differentiation between the infected and treated classes. Our findings reveal that the global stability of the HIV/AIDS-free equilibrium point is contingent upon the basic reproduction number being less than one. Furthermore, a bifurcation analysis demonstrates that our simplified model consistently exhibits a transcritical bifurcation at a reproduction number equal to one. In the complete model, we elucidate how the control reproduction number determines the stability of the HIV/AIDS-free equilibrium point. To align our model with the empirical data, we estimate its parameters using prevalence data from the top four countries affected by HIV/AIDS, namely, Eswatini, Lesotho, Botswana, and South Africa. We employ numerical simulations and conduct elasticity and sensitivity analyses to examine how our model parameters influence the control reproduction number and the dynamics of each model compartment. Our findings reveal that each country displays distinct sensitivities to the model parameters, implying the need for tailored strategies depending on the target country. Autonomous simulations highlight the potential of case detection and condom use in reducing HIV/AIDS prevalence. Furthermore, we identify that the quality of condoms plays a crucial role: with higher quality condoms, a smaller proportion of infected individuals need to use them for the potential eradication of HIV/AIDS from the population. In our optimal control simulations, we assess population behavior when control interventions are treated as time-dependent variables. Our analysis demonstrates that a combination of condom use and case detection, as time-dependent variables, can significantly curtail the spread of HIV while maintaining an optimal cost of intervention. Moreover, our cost-effectiveness analysis indicates that the condom use intervention alone emerges as the most cost-effective strategy, followed by a combination of case detection and condom use, and finally, case detection as a standalone strategy.

SARS-CoV-2 Transmission in Alberta, British Columbia, and Ontario, Canada, January 2020-January 2022.

We estimated COVID-19 transmission potential and case burden by variant type in Alberta, British Columbia, and Ontario, Canada, during January 23, 2020-January 27, 2022; we also estimated the effectiveness of public health interventions to reduce transmission. We estimated time-varying reproduction number (Rt) over 7-day sliding windows and nonoverlapping time-windows determined by timing of policy changes. We calculated incidence rate ratios (IRRs) for each variant and compared rates to determine differences in burden among provinces. Rt corresponding with emergence of the Delta variant increased in all 3 provinces; British Columbia had the largest increase, 43.85% (95% credible interval [CrI] 40.71%-46.84%). Across the study period, IRR was highest for Omicron (8.74 [95% CrI 8.71-8.77]) and burden highest in Alberta (IRR 1.80 [95% CrI 1.79-1.81]). Initiating public health interventions was associated with lower Rt and relaxing restrictions and emergence of new variants associated with increases in Rt.

R 0 May Not Tell Us Everything: Transient Disease Dynamics of Some SIR Models Over Patchy Environments.

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.