COVID-19: Development of a robust mathematical model and simulation package with consideration for ageing population and time delay for control action and resusceptibility

Trends of the COVID-19 dynamics in 2022 and 2023 vs. the population age, testing and vaccination levels

Introduction

The population, governments, and researchers show much less interest in the COVID-19 pandemic. However, many questions still need to be answered: why the much less vaccinated African continent has accumulated 15 times less deaths per capita than Europe? or why in 2023 the global value of the case fatality risk is almost twice higher than in 2022 and the UK figure is four times higher than the global one?

Methods

The averaged daily numbers of cases DCC and death DDC per million, case fatality risks DDC/DCC were calculated for 34 countries and regions with the use of John Hopkins University (JHU) datasets. Possible linear and non-linear correlations with the averaged daily numbers of tests per thousand DTC, median age of population A, and percentages of vaccinations VC and boosters BC were investigated.

Results

Strong correlations between age and DCC and DDC values were revealed. One-year increment in the median age yielded 39.8 increase in DCC values and 0.0799 DDC increase in 2022 (in 2023 these figures are 5.8 and 0.0263, respectively). With decreasing of testing level DTC, the case fatality risk can increase drastically. DCC and DDC values increase with increasing the percentages of fully vaccinated people and boosters, which definitely increase for greater A. After removing the influence of age, no correlations between vaccinations and DCC and DDC values were revealed.

Discussion

The presented analysis demonstrates that age is a pivot factor of visible (registered) part of the COVID-19 pandemic dynamics. Much younger Africa has registered less numbers of cases and death per capita due to many unregistered asymptomatic patients. Of great concern is the fact that COVID-19 mortality in 2023 in the UK is still at least 4 times higher than the global value caused by seasonal flu.

Implementation of a triangular probabilistic distribution for optimal parametrization of the SEIR model recovery rates with delay

Dynamical analysis of epidemiological models received significant attention after the global disaster of the Covid-19 pandemic. Although much of the attention has focused on choosing the appropriate model to describe and interpret the epidemiological data and on predicting the spread of the disease, working on uncovering the hidden dynamics by analyzing the predictive data and studying the algebraic properties of the models such as the invariant spaces, the center manifolds, and the Lyapunov functions in the general parameter space still have the potential to make a significant contribution to our understanding of the dynamics of spread. In this paper, in order to determine the recovery rate based on characteristic model parameters, the triangular probabilistic distribution is implemented in an epidemic delay differential equation with delays. By defining generic coefficients derived from field surveys and population characteristics, we propose a way to reduce the number of model parameters that need to be obtained during model fitting in order to find an effective way to handle the cultural and physiological diversity of societies and isolate the effect of counter measures for the epidemics.

Vaccination and isolation based control design of the COVID-19 pandemic based on adaptive neuro fuzzy inference system optimized with the genetic algorithm

The study of the COVID-19 pandemic is of pivotal importance due to its tremendous global impacts. This paper aims to control this disease using an optimal strategy comprising two methods: isolation and vaccination. In this regard, an optimized Adaptive Neuro-Fuzzy Inference System (ANFIS) is developed using the Genetic Algorithm (GA) to control the dynamic model of the COVID-19 termed SIDARTHE (Susceptible, Infected, Diagnosed, Ailing, Recognized, Threatened, Healed, and Extinct). The number of diagnosed and recognized people is reduced by isolation, and the number of susceptible people is reduced by vaccination. The GA generates optimal control efforts related to the random initial number of each chosen group as the input data for ANFIS to train Takagi-Sugeno (T-S) fuzzy structure coefficients. Also, three theorems are presented to indicate the positivity, boundedness, and existence of the solutions in the presence of the controller. The performance of the proposed system is evaluated through the mean squared error (MSE) and the root-mean-square error (RMSE). The simulation results show a significant decrease in the number of diagnosed, recognized, and susceptible individuals by employing the proposed controller, even with a 70% increase in transmissibility caused by various variants.

Mathematical Modeling and Numerical Simulation for the Outbreak of COVID-19 Involving Loss of Immunity and Quarantined Class

For the analysis of the recent deadly pandemic Sars-Cov-2, we constructed the mathematical model containing the whole population, partitioned into five different compartments, represented by SEIQR model. This current model especially contains the quarantined class and the factor of loss of immunity. Further, we discussed the stability of the SEIQR model (constructed on the basis of system of coupled differential equations). The basic reproduction that indicates the behavior of the disease is also estimated by the use of next-generation matrix method. Numerical simulation of this model is provided, the results are analyzed by theoretically strong numerical methods, and computationally known tool MATLAB Simulink is also used for visualization of the results. Validation of results by Simulink software and numerical methods shows that our model and adopted methodology are appropriate and accurate and could be used for further predictions on COVID-19. Our results suggest that the isolation of the active cases and strong immunization of patients or individuals play a major role to fight against the deadly Sars-Cov-2.

Comparing the dynamics of COVID-19 infection and mortality in the United States, India, and Brazil

This paper compares and contrasts the spread and impact of COVID-19 in the three countries most heavily impacted by the pandemic: the United States (US), India and Brazil. All three of these countries have a federal structure, in which the individual states have largely determined the response to the pandemic. Thus, we perform an extensive analysis of the individual states of these three countries to determine patterns of similarity within each. First, we analyse structural similarity and anomalies in the trajectories of cases and deaths as multivariate time series. Next, we study the lengths of the different waves of the virus outbreaks across the three countries and their states. Finally, we investigate suitable time offsets between cases and deaths as a function of the distinct outbreak waves. In all these analyses, we consistently reveal more characteristically distinct behaviour between US and Indian states, while Brazilian states exhibit less structure in their wave behaviour and changing progression between cases and deaths.

Uncertainty quantification in Covid-19 spread: Lockdown effects

We develop a Bayesian inference framework to quantify uncertainties in epidemiological models. We use SEIJR and SIJR models involving populations of susceptible, exposed, infective, diagnosed, dead and recovered individuals to infer from Covid-19 data rate constants, as well as their variations in response to lockdown measures. To account for confinement, we distinguish two susceptible populations at different risk: confined and unconfined. We show that transmission and recovery rates within them vary in response to facts, and that the diagnose rate is quite low, which leads to large amounts of undiagnosed infective individuals. A key unknown to predict the evolution of the epidemic is the fraction of the population affected by the virus, including asymptomatic subjects. Our study tracks its time evolution with quantified uncertainty from available official data, limited, however, by the data quality. We exemplify the technique with data from Spain, country in which late drastic lockdowns were enforced for months during the first wave of the current pandemic. In late actions and in the absence of other measures, spread is delayed but not stopped unless a large enough fraction of the population is confined until the asymptomatic population is depleted. To some extent, confinement can be replaced by strong distancing through masks in adequate circumstances.

Optimal COVID-19 lockdown strategies in an age-structured SEIR model of Northern Ireland

An age-structured SEIR model simulates the propagation of COVID-19 in the population of Northern Ireland. It is used to identify optimal timings of short-term lockdowns that enable long-term pandemic exit strategies by clearing the threshold for herd immunity or achieving time for vaccine development with minimal excess deaths.

Parameter identification in epidemiological models

We develop a Bayesian inference framework to estimate parameters in epidemiological models with quantified uncertainty. We consider SEIJR models involving populations of susceptible, exposed, infective, diagnosed, dead, and recovered individuals, and infer from data all the parameter models, in particular, the transmission and recovery rates. We exemplify the procedure using Covid-19 data from Spain since the onset of the current pandemic. Successive nonpharmaceutical actions such as lockdowns, distancing, and release of mobility restrictions, define stages in the data which are reflected in the parameter values. Tracking the evolution of the different populations with time, we can infer the evolution of the total fraction of subjects affected by the virus, including asymptomatic individuals. Using the resulting models as constraints in optimization problems for adequate costs we can gain insight on the parameter regimes in which the epidemic would remain controlled, even when migration effects are included.

Estimating a continuously varying offset between multivariate time series with application to COVID-19 in the United States

This paper introduces new methods to track the offset between two multivariate time series on a continuous basis. We then apply this framework to COVID-19 counts on a state-by-state basis in the United States to determine the progression from cases to deaths as a function of time. Across multiple approaches, we reveal an "up-down-up" pattern in the estimated offset between reported cases and deaths as the pandemic progresses. This analysis could be used to predict imminent increased load on a healthcare system and aid the allocation of additional resources in advance.

Estimation of COVID-19 recovery and decease periods in Canada using delay model

We derive a novel model escorted by large scale compartments, based on a set of coupled delay differential equations with extensive delays, in order to estimate the incubation, recovery and decease periods of COVID-19, and more generally any infectious disease. This is possible thanks to some optimization algorithms applied to publicly available database of confirmed corona cases, recovered cases and death toll. In this purpose, we separate (1) the total cases into 14 groups corresponding to 14 incubation periods, (2) the recovered cases into 406 groups corresponding to a combination of incubation and recovery periods, and (3) the death toll into 406 groups corresponding to a combination of incubation and decease periods. In this paper, we focus on recovery and decease periods and their correlation with the incubation period. The estimated mean recovery period we obtain is 22.14 days (95% Confidence Interval (CI) 22.00–22.27), and the 90th percentile is 28.91 days (95% CI 28.71–29.13), which is in agreement with statistical supported studies. The bimodal gamma distribution reveals that there are two groups of recovered individuals with a short recovery period, mean 21.02 days (95% CI 20.92–21.12), and a long recovery period, mean 38.88 days (95% CI 38.61–39.15). Our study shows that the characteristic of the decease period and the recovery period are alike. From the bivariate analysis, we observe a high probability domain for recovered individuals with respect to incubation and recovery periods. A similar domain is obtained for deaths analyzing bivariate distribution of incubation and decease periods.

Impact of time delays and environmental noise on the extinction of a population dynamics model

ABSTRACT

In this paper, we examine a population model with carrying capacity, time delay, and sources of additive and multiplicative environmental noise. We find that time delay, noise sources and their correlation induce regime shifts and transitions between the population survival state and the extinction state. To explore the transition mechanism between these two states, we analyzed the shift time to extinction, or the delayed extinction time, of populations. The main finding is that the extinction transition time as a function of the noise intensity shows a maximum, indicating the existence of an appropriate noise intensity leading to a maximal delayed extinction. This nonmonotonic behavior, with a maximum, is a signature of the noise-enhanced stability phenomenon, observed in many physical and complex metastable systems. In particular, this maximum increases (or decreases) as the cross-correlation intensity or the delay time in the death process increases. Furthermore, the signal-to-noise ratio as a function of noise intensity shows a maximum, which is a signature of the stochastic resonance phenomenon in the population dynamics model investigated in the presence of time delay and environmental noise.

Trends in COVID-19 prevalence and mortality: A year in review

This paper introduces new methods to study the changing dynamics of COVID-19 cases and deaths among the 50 worst-affected countries throughout 2020. First, we analyse the trajectories and turning points of rolling mortality rates to understand at which times the disease was most lethal. We demonstrate five characteristic classes of mortality rate trajectories and determine structural similarity in mortality trends over time. Next, we introduce a class of virulence matrices to study the evolution of COVID-19 cases and deaths on a global scale. Finally, we introduce three-way inconsistency analysis to determine anomalous countries with respect to three attributes: countries' COVID-19 cases, deaths and human development indices. We demonstrate the most anomalous countries across these three measures are Pakistan, the United States and the United Arab Emirates.

Nonlinear science against the COVID-19 pandemic

This special issue showcases recent uses of mathematical and nonlinear science methods in the study of different problems arising in the context of the COVID-19 pandemic. The sixteen original research papers included in this collection span a wide spectrum of studies including classical epidemiological models, new models accounting for COVID-19 specificities, non-pharmaceutical control measures, network models and other problems related to the pandemic.

Simulation of COVID-19 outbreaks via Graphical User Interface (GUI)

BACKGROUND

This research aimed to model the outbreak of COVID-19 in Malaysia and develop a GUI-based model.

DESIGN AND METHODS

The model is an improvement of the susceptible, infected, recovery, and death (SIRD) compartmental model.  The epidemiological parameters of the infection, recovery, and death rates were formulated as time dependent piecewise functions by incorporating the control measures of lockdown, social distancing, quarantine, lockdown lifting time and the percentage of people who abide by the rules. An improved SIRD model was solved via the 4th order Runge-Kutta (RK4) method and 14 unknown parameters were estimated by using Nelder-Mead algorithm and pattern-search technique. The publicly available data for COVID-19 outbreak in Malaysia was used to validate the performance of the model. The GUI-based SIRD model was developed to simulate the number of active cases of COVID-19 over time by considering movement control order (MCO) lifted date and the percentage of people who abide the rules.

RESULTS

The simulator showed that the improved SIRD model adequately fitted Malaysia COVID-19 data indicated by low values of root mean square error (RMSE) as compared to other existing models. The higher the percentage of people following the SOP, the lower the spread of disease. Another key point is that the later the lifting time after the lockdown, the lower the spread of disease.

CONCLUSION

These findings highlight the importance of the society to obey the intervention measures in preventing the spread of the COVID-19 disease.

COVID-19 modelling by time-varying transmission rate associated with mobility trend of driving via Apple Maps

Compartment-based infectious disease models that consider the transmission rate (or contact rate) as a constant during the course of an epidemic can be limiting regarding effective capture of the dynamics of infectious disease. This study proposed a novel approach based on a dynamic time-varying transmission rate with a control rate governing the speed of disease spread, which may be associated with the information related to infectious disease intervention. Integration of multiple sources of data with disease modelling has the potential to improve modelling performance. Taking the global mobility trend of vehicle driving available via Apple Maps as an example, this study explored different ways of processing the mobility trend data and investigated their relationship with the control rate. The proposed method was evaluated based on COVID-19 data from six European countries. The results suggest that the proposed model with dynamic transmission rate improved the performance of model fitting and forecasting during the early stage of the pandemic. Positive correlation has been found between the average daily change of mobility trend and control rate. The results encourage further development for incorporation of multiple resources into infectious disease modelling in the future.

Merits and Limitations of Mathematical Modeling and Computational Simulations in Mitigation of COVID-19 Pandemic: A Comprehensive Review

Mathematical models have assisted in describing the transmission and propagation dynamics of various viral diseases like MERS, measles, SARS, and Influenza; while the advanced computational technique is utilized in the epidemiology of viral diseases to examine and estimate the influences of interventions and vaccinations. In March 2020, the World Health Organization (WHO) has declared the COVID-19 as a global pandemic and the rate of morbidity and mortality triggers unprecedented public health crises throughout the world. The mathematical models can assist in improving the interventions, key transmission parameters, public health agencies, and countermeasures to mitigate this pandemic. Besides, the mathematical models were also used to examine the characteristics of epidemiological and the understanding of the complex transmission mechanism. Our literature study found that there were still some challenges in mathematical modeling for the case of ecology, genetics, microbiology, and pathology pose; also, some aspects like political and societal issues and cultural and ethical standards are hard to be characterized. Here, the recent mathematical models about COVID-19 and their prominent features, applications, limitations, and future perspective are discussed and reviewed. This review can assist in further improvement of mathematical models that will consider the current challenges of viral diseases.

Efficiency of communities and financial markets during the 2020 pandemic

This paper investigates the relationship between the spread of the COVID-19 pandemic, the state of community activity, and the financial index performance across 20 countries. First, we analyze which countries behaved similarly in 2020 with respect to one of three multivariate time series: daily COVID-19 cases, Apple mobility data, and national equity index price. Next, we study the trajectories of all three of these attributes in conjunction to determine which exhibited greater similarity. Finally, we investigate whether country financial indices or mobility data responded more quickly to surges in COVID-19 cases. Our results indicate that mobility data and national financial indices exhibited the most similarity in their trajectories, with financial indices responding quicker. This suggests that financial market participants may have interpreted and responded to COVID-19 data more efficiently than governments. Furthermore, results imply that efforts to study community mobility data as a leading indicator for financial market performance during the pandemic were misguided.

Detections and SIR simulations of the COVID-19 pandemic waves in Ukraine

Background. Unfortunately, the COVID-19 pandemic is still far from stabilizing. Of particular concern is the sharp increase in the number of diseases in June-July, September-October 2020 and February-March 2021. The causes and consequences of this sharp increase in the number of cases are still waiting for their researchers, but there is already an urgent need to assess the possible duration of the pandemic, the expected number of patients and deaths. Correct simulation of the infectious disease dynamics needs complicated mathematical models and many efforts for unknown parameters identification. Constant changes in the pandemic conditions (in particular, the peculiarities of quarantine and its violation, situations with testing and isolation of patients) cause various epidemic waves, lead to changes in the parameter values of the mathematical models.

Objective. In this article, pandemic waves in Ukraine will be detected, calculated and discussed. The estimations for durations and final sizes of the epidemic waves will be presented.

Methods. We propose a simple method for the epidemic waves detection based on the differentiation of the smoothed number of cases. We use the generalized SIR (susceptible-infected-removed) model for the dynamics of the epidemic waves. The known exact solution of the SIR differential equations and statistical approach were used. We will use different data sets for accumulated number of cases in order to compare the results of simulations and predictions.

Results. Nine pandemic waves were detected in Ukraine and corresponding optimal values of the SIR model parameters were identified. The number of cases and the number of patients spreading the infection versus time were calculated. In particular, the pandemic in Ukraine probably began in January 2020. If current trends continue, the end of the pandemic should be expected no earlier than in summer 2021.

Conclusions. The differentiation of the smoothed number of cases, the SIR model and statistical approach to the parameter identification are helpful to select COVID-19 pandemic waves and make some reliable estimations and predictions. The obtained information will be useful to regulate the quarantine activities, to predict the medical and economic consequences of the pandemic.

A simple but complex enough θ -SIR type model to be used with COVID-19 real data. Application to the case of Italy

Since the start of the COVID-19 pandemic in China many models have appeared in the literature, trying to simulate its dynamics. Focusing on modeling the biological and sociological mechanisms which influence the disease spread, the basic reference example is the SIR model. However, it is too simple to be able to model those mechanisms (including the three main types of control measures: social distancing, contact tracing and health system measures) to fit real data and to simulate possible future scenarios. A question, then, arises: how much and how do we need to complexify a SIR model? We develop a θ -SEIHQRD model, which may be the simplest one satisfying the mentioned requirements for arbitrary territories and can be simplified in particular cases. We show its very good performance in the Italian case and study different future scenarios.

Distribution of incubation periods of COVID-19 in the Canadian context

We propose a novel model based on a set of coupled delay differential equations with fourteen delays in order to accurately estimate the incubation period of COVID-19, employing publicly available data of confirmed corona cases. In this goal, we separate the total cases into fourteen groups for the corresponding fourteen incubation periods. The estimated mean incubation period we obtain is 6.74 days (95% Confidence Interval(CI): 6.35 to 7.13), and the 90th percentile is 11.64 days (95% CI: 11.22 to 12.17), corresponding to a good agreement with statistical supported studies. This model provides an almost zero-cost computational complexity to estimate the incubation period.

The Pandemic Holiday Blip in New York City

When it comes to pandemics, such as the one caused by the Coronavirus disease COVID-19, various issues and problems have arisen for the healthcare infrastructure and institutions. With increasing number of patients in need of urgent medical care and hospitalizations, the healthcare systems and regional hospitals may approach their maximum service capacity and may face shortage of various parameters, such as supplies including PPE, medications, therapeutic devices, ventilators, beds, and many more. The article at hand describes the development and framework of a simulation model that enables the modeling and evaluation of the COVID-19 pandemic progress. To achieve this, the model dynamically mimics and simulates the developments and time-dependent behavior of various crucial parameters of the pandemic, among others, the daily infection numbers and death rate. In addition, the model enables the simulation of single events and scenarios that occur outside of the regular pandemic developments as anomalies, such as holidays. Unlike traditional models, the proposed framework is based on factors and parameters closely derived from reality, such as the contact rate of individuals, which allows for a much more realistic representation. In addition, the real connection enables the assessment of effects of various influences regarding the development and progress of the pandemic, such as hospitalization numbers over time. All the aforementioned points are possible within the simulation framework and do not require awaiting the unfolding of the effects in reality. Thus, the model is capable of dynamically predicting how different scenarios turn out. The abilities of the model are demonstrated, illustrated, and proven in a specific case study that shows the impact of holidays, such as Passover and Easter in New York City when quarantine measures might have been ignored, and an increase in extended family gatherings temporarily occurred. As a result, the simulation showed significant impacts and disproportionate number of patients in need of medical care that could be potentially detrimental in reality. For example, compared to the previous trajectory of the pandemic, for a temporary increase of 50% in the contact rate of individuals, the model showed that the total number of cases would increase by 461 090, the maximum number of required hospitalizations would rise to 79 733, and the total number of fatalities would climb by 19 125 over 90 days. In addition to its function and proven capabilities, the model can and is furthermore planned to be adapted to other areas, not necessarily only metropolitan regions in order to expand the utilization of its predictive power. Such predictions could be used to derive regulatory measures and to test various policies for COVID-19 containment.

A dynamic model for the COVID-19 with direct and indirect transmission pathways

Two common transmission pathways for the spread of COVID-19 virus are direct and indirect. The direct pathway refers to the person-to-person transmission between susceptibles and infectious individuals. Infected individuals shed virus on the objects, and new infections arise through touching a contaminated surface; this refers to the indirect transmission pathway. We model the direct and indirect transmission pathways with a S A D O _I R ode model. Our proposal explicitly includes compartments for the contaminated objects, susceptible individuals, asymptomatic infectious, detected infectious, and recovered individuals. We compute the basic reproduction number and epidemic growth rate of the model and determine how these fundamental quantities relate to the transmission rate of the pathways. We further study the relationship between the rate of loss of immunity and the occurrence of backward bifurcation. An efficient statistical framework is introduced to estimate the parameters of the model. We show the performance of the model in the simulation scenarios and the real data from the COVID-19 daily cases in South Korea.

Visible and Real Sizes of New COVID-19 Pandemic Waves in Ukraine

Background. To simulate the COVID-19 pandemic dynamics, various data sets and different mathematical models can be used. In particular, previous simulations for Ukraine were based on smoothing of the dependence of the number of cases on time, classical and the generalized SIR (susceptible-infected-removed) models. Different simulation and comparison methods were based on official accumulated number of laboratory confirmed cases and the data reported by Johns Hopkins University. Since both datasets are incomplete (a very large percentage of infected persons are asymptomatic), the accuracy of calculations and predictions is limited. In this paper we will try to assess the degree of data incompleteness and correct the relevant forecasts. Objective. We aimed to estimate the real sizes of two new epidemic waves in Ukraine and compare them with visible dynamics based on the official number of laboratory confirmed cases. We also aimed to estimate the epidemic durations and final numbers of cases. Methods. In this study we use the generalized SIR model for the epidemic dynamics and its known exact solution. The known statistical approach is adopted in order to identify both the degree of data incompleteness and parameters of SIR model. Results. We have improved the method of estimating the unknown parameters of the generalized SIR model and calculated the optimal values ​​of the parameters. In particular, the visibility coefficients and the optimal values of the model parameters were estimated for two pandemic waves in Ukraine occurred in December 2020–March 2021. The real number of cases and the real number of patients spreading the infection versus time were calculated. Predictions of the real final sizes and durations of the pandemic in Ukraine are presented. If current trends continue, the end of the pandemic should be expected no earlier than in August 2022. Conclusions. New method of the unknown parameters identification for the generalized SIR model was proposed, which allows estimating the coefficients of data incompleteness as well. Its application for two pandemic waves in Ukraine has demonstrated that the real number of COVID-19 cases is approximately four times higher than those shown in official statistics. Probably, this situation is typical for other countries. The reassessments of the COVID-19 pandemic dynamics in other countries and clarification of world forecasts are necessary.

On a Discrete SEIR Epidemic Model with Two-Doses Delayed Feedback Vaccination Control on the Susceptible

A new discrete susceptible-exposed-infectious-recovered (SEIR) epidemic model is presented subject to a feedback vaccination effort involving two doses. Both vaccination doses, which are subject to a non-necessarily identical effectiveness, are administrated by respecting a certain mutual delay interval, and their immunity effect is registered after a certain delay since the second dose. The delays and the efficacies of the doses are parameters, which can be fixed in the model for each concrete experimentation. The disease-free equilibrium point is characterized as well as its stability properties, while it is seen that no endemic equilibrium point exists. The exposed subpopulation is supposed to be infective eventually, under a distinct transmission rate of that of the infectious subpopulation. Some simulation examples are presented by using disease parameterizations of the COVID-19 pandemic under vaccination efforts requiring two doses.

Mathematical Modeling of the Propagation of Covid-19 Pandemic Waves in the World

We develop a mathematical model of the coronavirus propagation in different countries (Brazil, India, US, Japan, Israel, Spain, Sweden), in the city of Moscow, and across the world. The pandemic spreads by a highly complex dynamics because it occurs in open nonhomogeneous systems where new infection foci erupt from time to time, triggering new transmission chains from infected to susceptible people. In general, statistical data collected as cumulative and epidemic curves are a superposition of many distinct local pandemic waves. In our modeling, we use the system of Feigenbaum’s discrete logistic equations (a logistic map) that describes the variation of the total number of infected over time. We show that this is the optimal model for the description of pandemic propagation in open nonhomogeneous systems with large errors in statistical data. We develop a procedure for isolating local waves, determining their model parameters, and predicting further evolution of each wave. We show that this model provides a good description of the statistical data and makes realistic forecasts. The forecast horizon depends on the degree of system closure and homogeneity. We calculate the start and end times of each wave, the peak, and the total number of infected in the current wave.

Sensitivity analysis of the infection transmissibility in the UK during the COVID-19 pandemic

The coronavirus (COVID-19) outbreak started in December 2019 and rapidly spread around the world affecting millions of people. With the growth of infection rate, many countries adopted different policies to control the spread of the disease. The UK implemented strict rules instructing individuals to stay at home except in some special circumstances starting from 23 March 2020. Accordingly, this study focuses on sensitivity analysis of transmissibility of the infection as the effects of removing restrictions, for example by returning different occupational groups to their normal working environment and its effect on the reproduction number in the UK. For this reason, available social contact matrices are adopted for the population of UK to account for the average number of contacts. Different scenarios are then considered to analyse the variability of total contacts on the reproduction number in the UK as a whole and each of its four nations. Our data-driven retrospective analysis shows that if more than 38.5% of UK working-age population return to their normal working environment, the reproduction number in the UK is expected to be higher than 1. However, analysis of each nation, separately, shows that local reproduction number in each nation may be different and requires more adequate analysis. Accordingly, we believe that using statistical methods and historical data can provide good estimation of local transmissibility and reproduction number in any region. As a consequence of this analysis, efforts to reduce the restrictions should be implemented locally via different control policies. It is important that these policies consider the social contacts, population density, and the occupational groups that are specific to each region.

Assessing the effects of time-dependent restrictions and control actions to flatten the curve of COVID-19 in Kazakhstan

This article presents the assessment of time-dependent national-level restrictions and control actions and their effects in fighting the COVID-19 pandemic. By analysing the transmission dynamics during the first wave of COVID-19 in the country, the effectiveness of the various levels of control actions taken to flatten the curve can be better quantified and understood. This in turn can help the relevant authorities to better plan for and control the subsequent waves of the pandemic. To achieve this, a deterministic population model for the pandemic is firstly developed to take into consideration the time-dependent characteristics of the model parameters, especially on the ever-evolving value of the reproduction number, which is one of the critical measures used to describe the transmission dynamics of this pandemic. The reproduction number alongside other key parameters of the model can then be estimated by fitting the model to real-world data using numerical optimisation techniques or by inducing ad-hoc control actions as recorded in the news platforms. In this article, the model is verified using a case study based on the data from the first wave of COVID-19 in the Republic of Kazakhstan. The model is fitted to provide estimates for two settings in simulations; time-invariant and time-varying (with bounded constraints) parameters. Finally, some forecasts are made using four scenarios with time-dependent control measures so as to determine which would reflect on the actual situations better.

On an SE(Is)(Ih)AR epidemic model with combined vaccination and antiviral controls for COVID-19 pandemic

In this paper, we study the nonnegativity and stability properties of the solutions of a newly proposed extended SEIR epidemic model, the so-called SE(Is)(Ih)AR epidemic model which might be of potential interest in the characterization and control of the COVID-19 pandemic evolution. The proposed model incorporates both asymptomatic infectious and hospitalized infectious subpopulations to the standard infectious subpopulation of the classical SEIR model. In parallel, it also incorporates feedback vaccination and antiviral treatment controls. The exposed subpopulation has three different transitions to the three kinds of infectious subpopulations under eventually different proportionality parameters. The existence of a unique disease-free equilibrium point and a unique endemic one is proved together with the calculation of their explicit components. Their local asymptotic stability properties and the attainability of the endemic equilibrium point are investigated based on the next generation matrix properties, the value of the basic reproduction number, and nonnegativity properties of the solution and its equilibrium states. The reproduction numbers in the presence of one or both controls is linked to the control-free reproduction number to emphasize that such a number decreases with the control gains. We also prove that, depending on the value of the basic reproduction number, only one of them is a global asymptotic attractor and that the solution has no limit cycles.

The spatial transmission of SARS-CoV-2 in China under the prevention and control measures at the early outbreak

Background

Since severe acute respiratory syndrome coronavirus, 2 (SARS-CoV-2) was firstly reported in Wuhan City, China in December 2019, Novel Coronavirus Disease 2019 (COVID-19) that is caused by SARS-CoV-2 is predominantly spread from person-to-person on worldwide scales. Now, COVID-19 is a non-traditional and major public health issue the world is facing, and the outbreak is a global pandemic. The strict prevention and control measures have mitigated the spread of SARS-CoV-2 and shown positive changes with important progress in China. But prevention and control tasks remain arduous for the world. The objective of this study is to discuss the difference of spatial transmission characteristics of COVID-19 in China at the early outbreak stage with resolute efforts. Simultaneously, the COVID-19 trend of China at the early time was described from the statistical perspective using a mathematical model to evaluate the effectiveness of the prevention and control measures.

Methods

In this study, the accumulated number of confirmed cases publicly reported by the National Health Committee of the People’s Republic of China (CNHC) from January 20 to February 11, 2020, were grouped into three partly overlapping regions: Chinese mainland including Hubei province, Hubei province alone, and the other 30 provincial-level regions on Chinese mainland excluding Hubei province, respectively. A generalized-growth model (GGM) was used to estimate the basic reproduction number to evaluate the transmissibility in different spatial locations. The prevention and control of COVID-19 in the early stage were analyzed based on the number of new cases of confirmed infections daily reported.

Results

Results indicated that the accumulated number of confirmed cases reported from January 20 to February 11, 2020, is well described by the GGM model with a larger correlation coefficient than 0.99. When the accumulated number of confirmed cases is well fitted by an exponential function, the basic reproduction number of COVID-19 of the 31 provincial-level regions on the Chinese mainland, Hubei province, and the other 30 provincial-level regions on the Chinese mainland excluding Hubei province, is 2.68, 6.46 and 2.18, respectively. The consecutive decline of the new confirmed cases indicated that the prevention and control measures taken by the Chinese government have contained the spread of SARS-CoV-2 in a short period.

Conclusions

The estimated basic reproduction number thorough GGM model can reflect the spatial difference of SARS-CoV-2 transmission in China at the early stage. The strict prevention and control measures of SARS-CoV-2 taken at the early outbreak can effectively reduce the new confirmed cases outside Hubei and have mitigated the spread and yielded positive results since February 2, 2020. The research results indicated that the outbreak of COVID-19 in China was sustaining localized at the early outbreak stage and has been gradually curbed by China’s resolute efforts.

An exploration of fractal-based prognostic model and comparative analysis for second wave of COVID-19 diffusion

The coronavirus disease 2019 (COVID-19) pandemic has fatalized 216 countries across the world and has claimed the lives of millions of people globally. Researches are being carried out worldwide by scientists to understand the nature of this catastrophic virus and find a potential vaccine for it. The most possible efforts have been taken to present this paper as a form of contribution to the understanding of this lethal virus in the first and second wave. This paper presents a unique technique for the methodical comparison of disastrous virus dissemination in two waves amid five most infested countries and the death rate of the virus in order to attain a clear view on the behaviour of the spread of the disease. For this study, the data set of the number of deaths per day and the number of infected cases per day of the most affected countries, the USA, Brazil, Russia, India, and the UK, have been considered in the first and second waves. The correlation fractal dimension has been estimated for the prescribed data sets of COVID-19, and the rate of death has been compared based on the correlation fractal dimension estimate curve. The statistical tool, analysis of variance, has also been used to support the performance of the proposed method. Further, the prediction of the daily death rate has been demonstrated through the autoregressive moving average model. In addition, this study also emphasis a feasible reconstruction of the death rate based on the fractal interpolation function. Subsequently, the normal probability plot is portrayed for the original data and the predicted data, derived through the fractal interpolation function to estimate the accuracy of the prediction. Finally, this paper neatly summarized with the comparison and prediction of epidemic curve of the first and second waves of COVID-19 pandemic to visualize the transmission rate in the both times.

Study on an SIHRS Model of COVID-19 Pandemic With Impulse and Time Delay Under Media Coverage

Media coverage plays an important role in prevention and control the spread of COVID-19 during the pandemic. In this paper, an SIHRS model of COVID-19 pandemic with impulse and time delay under media coverage is established. The positive and negative emotions of public are considered by the impact of confirmed cases and medical resources. In order to restrain the negative information of public, the factor of policies and regulations with impulse and time delay is introduced. Furthermore, the system model is simulated and verified by the reported data of COVID-19 pandemic in Wuhan. The main results are as follows: (1) When the implementation rate of the negative information generated by the confirmed cases gradually reduced to 0.4 times, the cumulative confirmed cases will be significantly reduced to about 37000, indicating that the popularization of pandemic related media information should be broad; (2) When the implementation rate affected by the amount of policies and regulations information gradually increases to 3 times, the cumulative confirmed cases will be significantly reduced to about 28000, indicating that the policies and regulations information should be continuously and incrementally reported; (3) When the inhibition rate of policies and regulation information on negative information gradually increases to 3 times, the cumulative confirmed cases will also be significantly reduced to about 27000 cases, indicating that the targeted policies and regulations information has a significant impact on inhibiting the corresponding negative emotions.

Parameter Estimation and Prediction of COVID-19 Epidemic Turning Point and Ending Time of a Case Study on SIR/SQAIR Epidemic Models

In this paper, the SIR epidemiological model for the COVID-19 with unknown parameters is considered in the first strategy. Three curves (S, I, and R) are fitted to the real data of South Korea, based on a detailed analysis of the actual data of South Korea, taken from the Korea Disease Control and Prevention Agency (KDCA). Using the least square method and minimizing the error between the fitted curve and the actual data, unknown parameters, like the transmission rate, recovery rate, and mortality rate, are estimated. The goodness of fit model is investigated with two criteria (SSE and RMSE), and the uncertainty range of the estimated parameters is also presented. Also, using the obtained determined model, the possible ending time and the turning point of the COVID-19 outbreak in the United States are predicted. Due to the lack of treatment and vaccine, in the next strategy, a new group called quarantined people is added to the proposed model. Also, a hidden state, including asymptomatic individuals, which is very common in COVID-19, is considered to make the model more realistic and closer to the real world. Then, the SIR model is developed into the SQAIR model. The delay in the recovery of the infected person is also considered as an unknown parameter. Like the previous steps, the possible ending time and the turning point in the United States are predicted. The model obtained in each strategy for South Korea is compared with the actual data from KDCA to prove the accuracy of the estimation of the parameters.

On a Controlled Se(Is)(Ih)(Iicu)AR Epidemic Model with Output Controllability Issues to Satisfy Hospital Constraints on Hospitalized Patients

An epidemic model, the so-called SE(Is)(Ih)(Iicu)AR epidemic model, is proposed which splits the infectious subpopulation of the classical SEIR (Susceptible-Exposed-Infectious-Recovered) model into four subpopulations, namely asymptomatic infectious and three categories of symptomatic infectious, namely slight infectious, non-intensive care infectious, and intensive care hospitalized infectious. The exposed subpopulation has four different transitions to each one of the four kinds of infectious subpopulations governed under eventually different proportionality parameters. The performed research relies on the problem of satisfying prescribed hospitalization constraints related to the number of patients via control interventions. There are four potential available controls which can be manipulated, namely the vaccination of the susceptible individuals, the treatment of the non-intensive care unit hospitalized patients, the treatment of the hospitalized patients at the intensive care unit, and the transmission rate which can be eventually updated via public interventions such as isolation of the infectious, rules of groups meetings, use of face masks, decrees of partial or total quarantines, and others. The patients staying at the non-intensive care unit and those staying at the intensive care unit are eventually, but not necessarily, managed as two different hospitalized subpopulations. The controls are designed based on output controllability issues in the sense that the levels of hospital admissions are constrained via prescribed maximum levels and the measurable outputs are defined by the hospitalized patients either under a joint consideration of the sum of both subpopulations or separately. In this second case, it is possible to target any of the two hospitalized subpopulations only or both of them considered as two different components of the output. Different algorithms are given to design the controls which guarantee, if possible, that the prescribed hospitalization constraints hold. If this were not possible, because the levels of serious infection are too high according to the hospital availability means, then the constraints are revised and modified accordingly so that the amended ones could be satisfied by a set of controls. The algorithms are tested through numerically worked examples under disease parameterizations of COVID-19.

European and US lockdowns and second waves during the COVID-19 pandemic

This paper investigates the lockdowns to contain the spread of the SARS-CoV-2 coronavirus in France, Germany, Italy, Spain, the UK and the US and also recent developments since these lockdowns have been relaxed. The analysis employs a two-stage SEIR model with different reproductive numbers pre- and post-lockdown. These parameters are estimated from data on the daily number of confirmed cases in a process that automatically detects the time at which the lockdown became effective. The model is evaluated by considering its predictive accuracy on current data and is then extended to a three-stage version to explore relaxations. The results show the extent to which each country was successful in reducing the reproductive number and demonstrate how the approach is able to model recent increases in the number of cases in all six countries, including the second peak in the US. The results also indicate that the current levels of relaxation in all five European countries could lead to significant second waves that last longer than the corresponding first waves. While there is uncertainty about the implications of these findings at this stage, they do suggest that a lot of vigilance is needed.

Intelligent computing with Levenberg-Marquardt artificial neural networks for nonlinear system of COVID-19 epidemic model for future generation disease control

The aim of this work is to design an intelligent computing paradigm through Levenberg-Marquardt artificial neural networks (LMANNs) for solving the mathematical model of Corona virus disease 19 (COVID-19) propagation via human to human interaction. The model is represented with systems of nonlinear ordinary differential equations represented with susceptible, exposed, symptomatic and infectious, super spreaders, infection but asymptomatic, hospitalized, recovery and fatality classes, and reference dataset of the COVID-19 model is generated by exploiting the strength of explicit Runge-Kutta numerical method for metropolitans of China and Pakistan including Wuhan, Karachi, Lahore, Rawalpindi and Faisalabad. The created dataset is arbitrary used for training, validation and testing processes for each cyclic update in Levenberg-Marquardt backpropagation for numerical treatment of the dynamics of COVID-19 model. The effectiveness and reliable performance of the design LMANNs are endorsed on the basis of assessments of achieved accuracy in terms of mean squared error based merit functions, error histograms and regression studies.

On Confinement and Quarantine Concerns on an SEIAR Epidemic Model with Simulated Parameterizations for the COVID-19 Pandemic

This paper firstly studies an SIR (susceptible-infectious-recovered) epidemic model without demography and with no disease mortality under both total and under partial quarantine of the susceptible subpopulation or of both the susceptible and the infectious ones in order to satisfy the hospital availability requirements on bed disposal and other necessary treatment means for the seriously infectious subpopulations. The seriously infectious individuals are assumed to be a part of the total infectious being described by a time-varying proportional function. A time-varying upper-bound of those seriously infected individuals has to be satisfied as objective by either a total confinement or partial quarantine intervention of the susceptible subpopulation. Afterwards, a new extended SEIR (susceptible-exposed-infectious-recovered) epidemic model, which is referred to as an SEIAR (susceptible-exposed-symptomatic infectious-asymptomatic infectious-recovered) epidemic model with demography and disease mortality is given and focused on so as to extend the above developed ideas on the SIR model. A proportionally gain in the model parameterization is assumed to distribute the transition from the exposed to the infectious into the two infectious individuals (namely, symptomatic and asymptomatic individuals). Such a model is evaluated under total or partial quarantines of all or of some of the subpopulations which have the effect of decreasing the number of contagions. Simulated numerical examples are also discussed related to model parameterizations of usefulness related to the current COVID-19 pandemic outbreaks.

Size and timescale of epidemics in the SIR framework

The most important features to assess the severity of an epidemic are its size and its timescale. We discuss these features in a systematic way in the context of SIR and SIR-type models. We investigate in detail how the size and timescale of the epidemic can be changed by acting on the parameters characterizing the model. Using these results and having as guideline the COVID-19 epidemic in Italy, we compare the efficiency of different containment strategies for contrasting an epidemic diffusion such as social distancing, lockdown, tracing, early detection and isolation.

Beyond just "flattening the curve": Optimal control of epidemics with purely non-pharmaceutical interventions

When effective medical treatment and vaccination are not available, non-pharmaceutical interventions such as social distancing, home quarantine and far-reaching shutdown of public life are the only available strategies to prevent the spread of epidemics. Based on an extended SEIR (susceptible-exposed-infectious-recovered) model and continuous-time optimal control theory, we compute the optimal non-pharmaceutical intervention strategy for the case that a vaccine is never found and complete containment (eradication of the epidemic) is impossible. In this case, the optimal control must meet competing requirements: First, the minimization of disease-related deaths, and, second, the establishment of a sufficient degree of natural immunity at the end of the measures, in order to exclude a second wave. Moreover, the socio-economic costs of the intervention shall be kept at a minimum. The numerically computed optimal control strategy is a single-intervention scenario that goes beyond heuristically motivated interventions and simple "flattening of the curve". Careful analysis of the computed control strategy reveals, however, that the obtained solution is in fact a tightrope walk close to the stability boundary of the system, where socio-economic costs and the risk of a new outbreak must be constantly balanced against one another. The model system is calibrated to reproduce the initial exponential growth phase of the COVID-19 pandemic in Germany.

Beyond just "flattening the curve": Optimal control of epidemics with purely non-pharmaceutical interventions

When effective medical treatment and vaccination are not available, non-pharmaceutical interventions such as social distancing, home quarantine and far-reaching shutdown of public life are the only available strategies to prevent the spread of epidemics. Based on an extended SEIR (susceptible-exposed-infectious-recovered) model and continuous-time optimal control theory, we compute the optimal non-pharmaceutical intervention strategy for the case that a vaccine is never found and complete containment (eradication of the epidemic) is impossible. In this case, the optimal control must meet competing requirements: First, the minimization of disease-related deaths, and, second, the establishment of a sufficient degree of natural immunity at the end of the measures, in order to exclude a second wave. Moreover, the socio-economic costs of the intervention shall be kept at a minimum. The numerically computed optimal control strategy is a single-intervention scenario that goes beyond heuristically motivated interventions and simple "flattening of the curve". Careful analysis of the computed control strategy reveals, however, that the obtained solution is in fact a tightrope walk close to the stability boundary of the system, where socio-economic costs and the risk of a new outbreak must be constantly balanced against one another. The model system is calibrated to reproduce the initial exponential growth phase of the COVID-19 pandemic in Germany.