Reaction-diffusion spatial modeling of COVID-19: Greece and Andalusia as case examples

On modelling airborne infection risk

Airborne infection risk analysis is usually performed for enclosed spaces where susceptible individuals are exposed to infectious airborne respiratory droplets by inhalation. It is usually based on exponential, dose-response models of which a widely used variant is the Wells-Riley (WR) model. We revisit this infection-risk estimate and extend it to the population level. We use an epidemiological model where the mode of pathogen transmission, airborne or contact, is explicitly considered. We illustrate the link between epidemiological models and the WR and the Gammaitoni and Nucci models. We argue that airborne infection quanta are, up to an overall density, airborne infectious respiratory droplets modified by a parameter that depends on biological properties of the pathogen, physical properties of the droplet and behavioural properties of the individual. We calculate the time-dependent risk of being infected for two scenarios. We show how the epidemic infection risk depends on the viral latent period and the event time, the time infection occurs. Infection risk follows the dynamics of the infected population. As the latent period decreases, infection risk increases. The longer a susceptible is present in the epidemic, the higher its risk of infection for equal exposure time to the pathogen is.

Agent-based modeling to estimate the impact of lockdown scenarios and events on a pandemic exemplified on SARS-CoV-2

In actual pandemic situations like COVID-19, it is important to understand the influence of single mitigation measures as well as combinations to create most dynamic impact for lockdown scenarios. Therefore we created an agent-based model (ABM) to simulate the spread of SARS-CoV-2 in an abstract city model with several types of places and agents. In comparison to infection numbers in Germany our ABM could be shown to behave similarly during the first wave. In our model, we implemented the possibility to test the effectiveness of mitigation measures and lockdown scenarios on the course of the pandemic. In this context, we focused on parameters of local events as possible mitigation measures and ran simulations, including varying size, duration, frequency and the proportion of events. The majority of changes to single event parameters, with the exception of frequency, showed only a small influence on the overall course of the pandemic. By applying different lockdown scenarios in our simulations, we could observe drastic changes in the number of infections per day. Depending on the lockdown strategy, we even observed a delayed peak in infection numbers of the second wave. As an advantage of the developed ABM, it is possible to analyze the individual risk of single agents during the pandemic. In contrast to standard or adjusted ODEs, we observed a 21% (with masks) / 48% (without masks) increased risk for single reappearing participants on local events, with a linearly increasing risk based on the length of the events.

A rigorous theoretical and numerical analysis of a nonlinear reaction-diffusion epidemic model pertaining dynamics of COVID-19

The spatial movement of the human population from one region to another and the existence of super-spreaders are the main factors that enhanced the disease incidence. Super-spreaders refer to the individuals having transmitting ability to multiple pathogens. In this article, an epidemic model with spatial and temporal effects is formulated to analyze the impact of some preventing measures of COVID-19. The model is developed using six nonlinear partial differential equations. The infectious individuals are sub-divided into symptomatic, asymptomatic and super-spreader classes. In this study, we focused on the rigorous qualitative analysis of the reaction-diffusion model. The fundamental mathematical properties of the proposed COVID-19 epidemic model such as boundedness, positivity, and invariant region of the problem solution are derived, which ensure the validity of the proposed model. The model equilibria and its stability analysis for both local and global cases have been presented. The normalized sensitivity analysis of the model is carried out in order to observe the crucial factors in the transmission of infection. Furthermore, an efficient numerical scheme is applied to solve the proposed model and detailed simulation are performed. Based on the graphical observation, diffusion in the context of confined public gatherings is observed to significantly inhibit the spread of infection when compared to the absence of diffusion. This is especially important in scenarios where super-spreaders may play a major role in transmission. The impact of some non-pharmaceutical interventions are illustrated graphically with and without diffusion. We believe that the present investigation will be beneficial in understanding the complex dynamics and control of COVID-19 under various non-pharmaceutical interventions.

Vaccination compartmental epidemiological models for the delta and omicron SARS-CoV-2 variants

We explore the inclusion of vaccination in compartmental epidemiological models concerning the delta and omicron variants of the SARS-CoV-2 virus that caused the COVID-19 pandemic. We expand on our earlier compartmental-model work by incorporating vaccinated populations. We present two classes of models that differ depending on the immunological properties of the variant. The first one is for the delta variant, where we do not follow the dynamics of the vaccinated individuals since infections of vaccinated individuals were rare. The second one for the far more contagious omicron variant incorporates the evolution of the infections within the vaccinated cohort. We explore comparisons with available data involving two possible classes of counts, fatalities and hospitalizations. We present our results for two regions, Andalusia and Switzerland (including the Principality of Liechtenstein), where the necessary data are available. In the majority of the considered cases, the models are found to yield good agreement with the data and have a reasonable predictive capability beyond their training window, rendering them potentially useful tools for the interpretation of the COVID-19 and further pandemic waves, and for the design of intervention strategies during these waves.

Parameter estimation for contact tracing in graph-based models

We adopt a maximum-likelihood framework to estimate parameters of a stochastic susceptible-infected-recovered (SIR) model with contact tracing on a rooted random tree. Given the number of detectees per index case, our estimator allows to determine the degree distribution of the random tree as well as the tracing probability. Since we do not discover all infectees via contact tracing, this estimation is non-trivial. To keep things simple and stable, we develop an approximation suited for realistic situations (contract tracing probability small, or the probability for the detection of index cases small). In this approximation, the only epidemiological parameter entering the estimator is R0. The estimator is tested in a simulation study and is furthermore applied to COVID-19 contact tracing data from India. The simulation study underlines the efficiency of the method. For the empirical COVID-19 data, we compare different degree distributions and perform a sensitivity analysis. We find that particularly a power-law and a negative binomial degree distribution fit the data well and that the tracing probability is rather large. The sensitivity analysis shows no strong dependency of the estimates on the reproduction number. Finally, we discuss the relevance of our findings.

The Impact of Spring Festival Travel on Epidemic Spreading in China

The large population movement during the Spring Festival travel in China can considerably accelerate the spread of epidemics, especially after the relaxation of strict control measures against COVID-19. This study aims to assess the impact of population migration in Spring Festival holiday on epidemic spread under different scenarios. Using inter-city population movement data, we construct the population flow network during the non-holiday time as well as the Spring Festival holiday. We build a large-scale metapopulation model to simulate the epidemic spread among 371 Chinese cities. We analyze the impact of Spring Festival travel on the peak timing and peak magnitude nationally and in each city. Assuming an R0 (basic reproduction number) of 15 and the initial conditions as the reported COVID-19 infections on 17 December 2022, model simulations indicate that the Spring Festival travel can substantially increase the national peak magnitude of infection. The infection peaks arrive at most cities 1-4 days earlier as compared to those of the non-holiday time. While peak infections in certain large cities, such as Beijing and Shanghai, are decreased due to the massive migration of people to smaller cities during the pre-Spring Festival period, peak infections increase significantly in small- or medium-sized cities. For a less transmissible disease (R0 = 5), infection peaks in large cities are delayed until after the Spring Festival. Small- or medium-sized cities may experience a larger infection due to the large-scale population migration from metropolitan areas. The increased disease burden may impose considerable strain on the healthcare systems in these resource-limited areas. For a less transmissible disease, particular attention needs to be paid to outbreaks in large cities when people resume work after holidays.

The Role of Mobility in the Dynamics of the COVID-19 Epidemic in Andalusia

Metapopulation models have been a popular tool for the study of epidemic spread over a network of highly populated nodes (cities, provinces, countries) and have been extensively used in the context of the ongoing COVID-19 pandemic. In the present work, we revisit such a model, bearing a particular case example in mind, namely that of the region of Andalusia in Spain during the period of the summer-fall of 2020 (i.e., between the first and second pandemic waves). Our aim is to consider the possibility of incorporation of mobility across the province nodes focusing on mobile-phone time-dependent data, but also discussing the comparison for our case example with a gravity model, as well as with the dynamics in the absence of mobility. Our main finding is that mobility is key toward a quantitative understanding of the emergence of the second wave of the pandemic and that the most accurate way to capture it involves dynamic (rather than static) inclusion of time-dependent mobility matrices based on cell-phone data. Alternatives bearing no mobility are unable to capture the trends revealed by the data in the context of the metapopulation model considered herein.

A numerical study of COVID-19 epidemic model with vaccination and diffusion

The coronavirus infectious disease (or COVID-19) is a severe respiratory illness. Although the infection incidence decreased significantly, still it remains a major panic for human health and the global economy. The spatial movement of the population from one region to another remains one of the major causes of the spread of the infection. In the literature, most of the COVID-19 models have been constructed with only temporal effects. In this paper, a vaccinated spatio-temporal COVID-19 mathematical model is developed to study the impact of vaccines and other interventions on the disease dynamics in a spatially heterogeneous environment. Initially, some of the basic mathematical properties including existence, uniqueness, positivity, and boundedness of the diffusive vaccinated models are analyzed. The model equilibria and the basic reproductive number are presented. Further, based upon the uniform and non-uniform initial conditions, the spatio-temporal COVID-19 mathematical model is solved numerically using finite difference operator-splitting scheme. Furthermore, detailed simulation results are presented in order to visualize the impact of vaccination and other model key parameters with and without diffusion on the pandemic incidence. The obtained results reveal that the suggested intervention with diffusion has a significant impact on the disease dynamics and its control.

Backcasting COVID-19: a physics-informed estimate for early case incidence

It is widely accepted that the number of reported cases during the first stages of the COVID-19 pandemic severely underestimates the number of actual cases. We leverage delay embedding theorems of Whitney and Takens and use Gaussian process regression to estimate the number of cases during the first 2020 wave based on the second wave of the epidemic in several European countries, South Korea and Brazil. We assume that the second wave was more accurately monitored, even though we acknowledge that behavioural changes occurred during the pandemic and region- (or country-) specific monitoring protocols evolved. We then construct a manifold diffeomorphic to that of the implied original dynamical system, using fatalities or hospitalizations only. Finally, we restrict the diffeomorphism to the reported cases coordinate of the dynamical system. Our main finding is that in the European countries studied, the actual cases are under-reported by as much as 50%. On the other hand, in South Korea-which had a proactive mitigation approach-a far smaller discrepancy between the actual and reported cases is predicted, with an approximately 18% predicted underestimation. We believe that our backcasting framework is applicable to other epidemic outbreaks where (due to limited or poor quality data) there is uncertainty around the actual cases.

Modelling of spatial infection spread through heterogeneous population: from lattice to partial differential equation models

We present a simple model for the spread of an infection that incorporates spatial variability in population density. Starting from first-principle considerations, we explore how a novel partial differential equation with state-dependent diffusion can be obtained. This model exhibits higher infection rates in the areas of higher population density-a feature that we argue to be consistent with epidemiological observations. The model also exhibits an infection wave, the speed of which varies with population density. In addition, we demonstrate the possibility that an infection can 'jump' (i.e. tunnel) across areas of low population density towards areas of high population density. We briefly touch upon the data reported for coronavirus spread in the Canadian province of Nova Scotia as a case example with a number of qualitatively similar features as our model. Lastly, we propose a number of generalizations of the model towards future studies.

The COVID-19 pandemic as inspiration to reconsider epidemic models: A novel approach to spatially homogeneous epidemic spread modeling

Epidemic spread models are useful tools to study the spread and the effectiveness of the interventions at a population level, to an epidemic. The workhorse of spatially homogeneous class models is the SIR-type ones comprising ordinary differential equations for the unknown state variables. The transition between different states is expressed through rate functions. Inspired by -but not restricted to- features of the COVID-19 pandemic, a new framework for modeling a disease spread is proposed. The main concept refers to the assignment of properties to each individual person as regards his response to the disease. A multidimensional distribution of these properties represents the whole population. The temporal evolution of this distribution is the only dependent variable of the problem. All other variables can be extracted by post-processing of this distribution. It is noteworthy that the new concept allows an improved consideration of vaccination modeling because it recognizes vaccination as a modifier of individuals response to the disease and not as a means for individuals to totally defeat the disease. At the heart of the new approach is an infection age model engaging a sharp cut-off. This model is analyzed in detail, and it is shown to admit self-similar solutions. A hierarchy of models based on the new approach, from a generalized one to a specific one with three dominant properties, is derived. The latter is implemented as an example and indicative results are presented and discussed. It appears that the new framework is general and versatile enough to simulate disease spread processes and to predict the evolution of several variables of the population during this spread.

Data-driven reduced-order modeling of spatiotemporal chaos with neural ordinary differential equations

Dissipative partial differential equations that exhibit chaotic dynamics tend to evolve to attractors that exist on finite-dimensional manifolds. We present a data-driven reduced-order modeling method that capitalizes on this fact by finding a coordinate representation for this manifold and then a system of ordinary differential equations (ODEs) describing the dynamics in this coordinate system. The manifold coordinates are discovered using an undercomplete autoencoder-a neural network (NN) that reduces and then expands dimension. Then, the ODE, in these coordinates, is determined by a NN using the neural ODE framework. Both of these steps only require snapshots of data to learn a model, and the data can be widely and/or unevenly spaced. Time-derivative information is not needed. We apply this framework to the Kuramoto-Sivashinsky equation for domain sizes that exhibit chaotic dynamics with again estimated manifold dimensions ranging from 8 to 28. With this system, we find that dimension reduction improves performance relative to predictions in the ambient space, where artifacts arise. Then, with the low-dimensional model, we vary the training data spacing and find excellent short- and long-time statistical recreation of the true dynamics for widely spaced data (spacing of ∼ 0.7 Lyapunov times). We end by comparing performance with various degrees of dimension reduction and find a "sweet spot" in terms of performance vs dimension.

Modeling epidemic flow with fluid dynamics

In this paper, a new mathematical model based on partial differential equations is proposed to study the spatial spread of infectious diseases. The model incorporates fluid dynamics theory and represents the epidemic spread as a fluid motion generated through the interaction between the susceptible and infected hosts. At the macroscopic level, the spread of the infection is modeled as an inviscid flow described by the Euler equation. Nontrivial numerical methods from computational fluid dynamics (CFD) are applied to investigate the model. In particular, a fifth-order weighted essentially non-oscillatory (WENO) scheme is employed for the spatial discretization. As an application, this mathematical and computational framework is used in a simulation study for the COVID-19 outbreak in Wuhan, China. The simulation results match the reported data for the cumulative cases with high accuracy and generate new insight into the complex spatial dynamics of COVID-19.

An integrated simulation framework for the prevention and mitigation of pandemics caused by airborne pathogens

In this work, we developed an integrated simulation framework for pandemic prevention and mitigation of pandemics caused by airborne pathogens, incorporating three sub-models, namely the spatial model, the mobility model, and the propagation model, to create a realistic simulation environment for the evaluation of the effectiveness of different countermeasures on the epidemic dynamics. The spatial model converts images of real cities obtained from Google Maps into undirected weighted graphs that capture the spatial arrangement of the streets utilized next for the mobility of individuals. The mobility model implements a stochastic agent-based approach, developed to assign specific routes to individuals moving in the city, through the use of stochastic processes, utilizing the weights of the underlying graph to deploy shortest path algorithms. The propagation model implements both the epidemiological model and the physical substance of the transmission of an airborne pathogen (in our approach, we investigate the transmission parameters of SARS-CoV-2). The deployment of a set of countermeasures was investigated in reducing the spread of the pathogen, where, through a series of repetitive simulation experiments, we evaluated the effectiveness of each countermeasure in pandemic prevention.