A Mathematical Description of the Dynamics of Coronavirus Disease 2019 (COVID-19): A Case Study of Brazil

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 ifR 0 < 1 and endemic equilibrium is locally asymptotically stable ifR 0 > 1 . Moreover, the disease-free equilibrium is globally asymptotically stable ifR 0 < 1 and endemic equilibrium is globally asymptotically stable ifR 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.

Mathematical modeling to study the impact of immigration on the dynamics of the COVID-19 pandemic: A case study for Venezuela

We propose two different mathematical models to study the effect of immigration on the COVID-19 pandemic. The first model does not consider immigration, whereas the second one does. Both mathematical models consider five different subpopulations: susceptible, exposed, infected, asymptomatic carriers, and recovered. We find the basic reproduction number R₀ using the next-generation matrix method for the mathematical model without immigration. This threshold parameter is paramount because it allows us to characterize the evolution of the disease and identify what parameters substantially affect the COVID-19 pandemic outcome. We focus on the Venezuelan scenario, where immigration and emigration have been important over recent years, particularly during the pandemic. We show that the estimation of the transmission rates of the SARS-CoV-2 are affected when the immigration of infected people is considered. This has an important consequence from a public health perspective because if the basic reproduction number is less than unity, we can expect that the SARS-CoV-2 would disappear. Thus, if the basic reproduction number is slightly above one, we can predict that some mild non-pharmaceutical interventions would be enough to decrease the number of infected people. The results show that the dynamics of the spread of SARS-CoV-2 through the population must consider immigration to obtain better insight into the outcomes and create awareness in the population regarding the population flow.

Recursive state and parameter estimation of COVID-19 circulating variants dynamics

COVID-19 pandemic response with non-pharmaceutical interventions is an intrinsic control problem. Governments weigh social distancing policies to avoid overload in the health system without significant economic impact. The mutability of the SARS-CoV-2 virus, vaccination coverage, and mobility restriction measures change epidemic dynamics over time. A model-based control strategy requires reliable predictions to be efficient on a long-term basis. In this paper, a SEIR-based model is proposed considering dynamic feedback estimation. State and parameter estimations are performed on state estimators using augmented states. Three methods were implemented: constrained extended Kalman filter (CEKF), CEKF and smoother (CEKF & S), and moving horizon estimator (MHE). The parameters estimation was based on vaccine efficacy studies regarding transmissibility, severity of the disease, and lethality. Social distancing was assumed as a measured disturbance calculated using Google mobility data. Data from six federative units from Brazil were used to evaluate the proposed strategy. State and parameter estimations were performed from 1 October 2020 to 1 July 2021, during which Zeta and Gamma variants emerged. Simulation results showed that lethality increased between 11 and 30% for Zeta mutations and between 44 and 107% for Gamma mutations. In addition, transmissibility increased between 10 and 37% for the Zeta variant and between 43 and 119% for the Gamma variant. Furthermore, parameter estimation indicated temporal underreporting changes in hospitalized and deceased individuals. Overall, the estimation strategy showed to be suitable for dynamic feedback as simulation results presented an efficient detection and dynamic characterization of circulating variants.

Analysis of COVID-19 in India using a vaccine epidemic model incorporating vaccine effectiveness and herd immunity

COVID-19 will be a continuous threat to human population despite having a few vaccines at hand until we reach the endemic state through natural herd immunity and total immunization through universal vaccination. However, the vaccine acts as a practical tool for reducing the massive public health problem and the emerging economic consequences that the continuing COVID -19 epidemic is causing worldwide, while the vaccine efficacy wanes. In this work, we propose and analyze an epidemic model of Susceptible-Exposed-Infected-Recovered-Vaccinated population taking into account the rate of vaccination and vaccine waning. The dynamics of the model has been investigated, and the condition for a disease-free endemic equilibrium state is obtained. Further, the analysis is extended to study the COVID-19 spread in India by considering the availability of vaccines and the related critical parameters such as vaccination rate, vaccine efficacy and waning of vaccine's impact on deciding the emerging fate of this epidemic. We have also discussed the conditions for herd immunity due to vaccinated individuals among the people. Our results highlight the importance of vaccines, the effectiveness of booster vaccination in protecting people from infection, and their importance in epidemic and pandemic modeling.

Computational and Mathematical Methods in Medicine Prediction of COVID-19 in BRICS Countries: An Integrated Deep Learning Model of CEEMDAN-R-ILSTM-Elman

Since the outbreak of COVID-19, BRICS countries have experienced different epidemic spread due to different health conditions, social isolation measures, vaccination rates, and other factors. A descriptive analysis is conducted for the spread of the epidemic in the BRICS countries. Considering the nonlinear and nonstationary characteristics of COVID-19 data, a principle of decomposition-reconstruction(R)-prediction-integration is proposed. Correspondingly, this paper constructs an integrated deep learning prediction model of CEEMDAN-R-ILSTM-Elman. Specifically, the prediction model is integrated by complete ensemble empirical mode decomposition (CEEMDAN), improved long-term and short-term memory network (ILSTM), and Elman neural network. First, the data is decomposed by adopting CEEMDAN. Then, by calculating the permutation entropy and average period, the decomposed eigenmode component IMFs are reconstructed into four sequences of high, medium, low level, and trend term. Thus, ILSTM and Elman algorithms are used for component sequence prediction, whose results are integrated as the final results. The ILSTM is established based on the LSTM model and the improved beetle antennae search algorithm (IBAS). The ILSTM mainly considers that the prediction accuracy of LSTM model is vulnerable to the influence of parameter selection. The IBAS with adaptive step size is used to automatically optimize the super parameters of LSTM model and to improve the modeling efficiency and prediction accuracy. Experimental results indicate that compared with other benchmark models, CEEMDAN-R-ILSTM-Elman integrated model predicts the number of newly confirmed cases of COVID-19 in BRICS countries with higher accuracy and lower error. Strict social policies have a greater impact on the infection rate and mortality rate of the epidemic. During July-August 2021, epidemic spread in BRICS countries will slow down, and the overall situation is still quite severe.

On mobility trends analysis of COVID-19 dissemination in Mexico City

This work presents a tool for forecasting the spread of the new coronavirus in Mexico City, which is based on a mathematical model with a metapopulation structure that uses Bayesian statistics and is inspired by a data-driven approach. The daily mobility of people in Mexico City is mathematically represented by an origin-destination matrix using the open mobility data from Google and the Transportation Mexican Survey. This matrix is incorporated in a compartmental model. We calibrate the model against borough-level incidence data collected between 27 February 2020 and 27 October 2020, while using Bayesian inference to estimate critical epidemiological characteristics associated with the coronavirus spread. Given that working with metapopulation models leads to rather high computational time consumption, and parameter estimation of these models may lead to high memory RAM consumption, we do a clustering analysis that is based on mobility trends to work on these clusters of borough separately instead of taken all of the boroughs together at once. This clustering analysis can be implemented in smaller or larger scales in different parts of the world. In addition, this clustering analysis is divided into the phases that the government of Mexico City has set up to restrict individual movement in the city. We also calculate the reproductive number in Mexico City using the next generation operator method and the inferred model parameters obtaining that this threshold is in the interval (1.2713, 1.3054). Our analysis of mobility trends can be helpful when making public health decisions.

Analysis of the second wave of COVID-19 in India based on SEIR model

India was under a grave threat from the second wave of the COVID-19 pandemic particularly in the beginning of May 2021. The situation appeared rather gloomy as the number of infected individuals/active cases had increased alarmingly during the months of May and June 2021 compared to the first wave peak. Indian government/state governments have been implementing various control measures such as lockdowns, setting up new hospitals, and putting travel restrictions at various stages to lighten the virus spread from the initial outbreak of the pandemic. Recently, we have studied the susceptible-exposed-infectious-removed (SEIR) dynamic modeling of the epidemic evolution of COVID-19 in India with the help of appropriate parameters quantifying the various governmental actions and the intensity of individual reactions. Our analysis had predicted the scenario of the first wave quite well. In this present article, we extend our analysis to estimate and analyze the number of infected individuals during the second wave of COVID-19 in India with the help of the above SEIR model. Our findings show that the people's individual effort along with governmental actions such as implementations of curfews and accelerated vaccine strategy are the most important factors to control the pandemic in the present situation and in the future.

Dynamical analysis of coronavirus disease with crowding effect, and vaccination: a study of third strain

Countries affected by the coronavirus epidemic have reported many infected cases and deaths based on world health statistics. The crowding factor, which we named "crowding effects," plays a significant role in spreading the diseases. However, the introduction of vaccines marks a turning point in the rate of spread of coronavirus infections. Modeling both effects is vastly essential as it directly impacts the overall population of the studied region. To determine the peak of the infection curve by considering the third strain, we develop a mathematical model (susceptible-infected-vaccinated-recovered) with reported cases from August 01, 2021, till August 29, 2021. The nonlinear incidence rate with the inclusion of both effects is the best approach to analyze the dynamics. The model's positivity, boundedness, existence, uniqueness, and stability (local and global) are addressed with the help of a reproduction number. In addition, the strength number and second derivative Lyapunov analysis are examined, and the model was found to be asymptotically stable. The suggested parameters efficiently control the active cases of the third strain in Pakistan. It was shown that a systematic vaccination program regulates the infection rate. However, the crowding effect reduces the impact of vaccination. The present results show that the model can be applied to other countries' data to predict the infection rate.

A dynamical map to describe COVID-19 epidemics

Nonlinear dynamics perspective is an interesting approach to describe COVID-19 epidemics, providing information to support strategic decisions. This paper proposes a dynamical map to describe COVID-19 epidemics based on the classical susceptible-exposed-infected-recovered (SEIR) differential model, incorporating vaccinated population. On this basis, the novel map represents COVID-19 discrete-time dynamics by adopting three populations: infected, cumulative infected and vaccinated. The map promotes a dynamical description based on algebraic equations with a reduced number of variables and, due to its simplicity, it is easier to perform parameter adjustments. In addition, the map description allows analytical calculations of useful information to evaluate the epidemic scenario, being important to support strategic decisions. In this regard, it should be pointed out the estimation of the number deaths, infection rate and the herd immunization point. Numerical simulations show the model capability to describe COVID-19 dynamics, capturing the main features of the epidemic evolution. Reported data from Germany, Italy and Brazil are of concern showing the map ability to describe different scenario patterns that include multi-wave pattern with bell shape and plateaus characteristics. The effect of vaccination is analyzed considering different campaign strategies, showing its importance to control the epidemics.

Experimental investigation on the susceptibility of minimal networks to a change in topology and number of oscillators

Understanding the global dynamical behavior of a network of coupled oscillators has been a topic of immense research in many fields of science and engineering. Various factors govern the resulting dynamical behavior of such networks, including the number of oscillators and their coupling schemes. Although these factors are seldom significant in large populations, a small change in them can drastically affect the global behavior in small populations. In this paper, we perform an experimental investigation on the effect of these factors on the coupled behavior of a minimal network of candle-flame oscillators. We observe that strongly coupled oscillators exhibit the global behavior of in-phase synchrony and amplitude death, irrespective of the number and the topology of oscillators. However, when they are weakly coupled, their global behavior exhibits the intermittent occurrence of multiple stable states in time. We report the experimental discovery of partial amplitude death in a network of candle-flame oscillators, in addition to the observation of other dynamical states including clustering, chimera, and weak chimera. We also show that closed-loop networks tend to hold global synchronization for longer duration as compared to open-loop networks.

COVID-19 dynamics considering the influence of hospital infrastructure: an investigation into Brazilian scenarios

COVID-19 dynamics is one of the most relevant subjects nowadays, and, in this regard, mathematical modeling and numerical simulations are of special interest. This paper describes COVID-19 dynamics based on a novel version of the susceptible-exposed-infectious-removed model. Removed population is split into recovered and death populations allowing a better comprehension of real situations. Besides, the total population is reduced based on the number of deaths. Hospital infrastructure is also included into the mathematical description allowing the consideration of collapse scenarios. Initially, a model verification is carried out calibrating system parameters with data from China outbreak that is considered a benchmark due the availability of data for the entire cycle. Afterward, Brazil outbreak is of concern, calibrating the model and developing numerical simulations. Results show several scenarios highlighting the importance of social isolation and hospital infrastructure. System dynamics has a strong sensitivity to transmission rate showing the importance of numerical simulations to guide public health decision strategies. Results also show that complex dynamical responses can emerge due to the oscillations of the transmission rate, being associated with distinct infection subsequent waves.

New approximations, and policy implications, from a delayed dynamic model of a fast pandemic

We study an SEIQR (Susceptible-Exposed-Infectious-Quarantined-Recovered) model due to Young et al. (2019) for an infectious disease, with time delays for latency and an asymptomatic phase. For fast pandemics where nobody has prior immunity and everyone has immunity after recovery, the SEIQR model decouples into two nonlinear delay differential equations (DDEs) with five parameters. One parameter is set to unity by scaling time. The simple subcase of perfect quarantining and zero self-recovery before quarantine, with two free parameters, is examined first. The method of multiple scales yields a hyperbolic tangent solution; and a long-wave (short delay) approximation yields a first order ordinary differential equation (ODE). With imperfect quarantining and nonzero self-recovery, the long-wave approximation is a second order ODE. These three approximations each capture the full outbreak, from infinitesimal initiation to final saturation. Low-dimensional dynamics in the DDEs is demonstrated using a six state non-delayed reduced order model obtained by Galerkin projection. Numerical solutions from the reduced order model match the DDE over a range of parameter choices and initial conditions. Finally, stability analysis and numerics show how a well executed temporary phase of social distancing can reduce the total number of people affected. The reduction can be by as much as half for a weak pandemic, and is smaller but still substantial for stronger pandemics. An explicit formula for the greatest possible reduction is given.