Bayesian-based predictions of COVID-19 evolution in Texas using multispecies mixture-theoretic continuum models

Scalable computational algorithms for geospatial COVID-19 spread using high performance computing

A nonlinear partial differential equation (PDE) based compartmental model of COVID-19 provides a continuous trace of infection over space and time. Finer resolutions in the spatial discretization, the inclusion of additional model compartments and model stratifications based on clinically relevant categories contribute to an increase in the number of unknowns to the order of millions. We adopt a parallel scalable solver that permits faster solutions for these high fidelity models. The solver combines domain decomposition and algebraic multigrid preconditioners at multiple levels to achieve the desired strong and weak scalabilities. As a numerical illustration of this general methodology, a five-compartment susceptible-exposed-infected-recovered-deceased (SEIRD) model of COVID-19 is used to demonstrate the scalability and effectiveness of the proposed solver for a large geographical domain (Southern Ontario). It is possible to predict the infections for a period of three months for a system size of 186 million (using 3200 processes) within 12 hours saving months of computational effort needed for the conventional solvers.

A machine learning-based probabilistic computational framework for uncertainty quantification of actuation of clustered tensegrity structures

Clustered tensegrity structures integrated with continuous cables are lightweight, foldable, and deployable. Thus, they can be used as flexible manipulators or soft robots. The actuation process of such soft structure has high probabilistic sensitivity. It is essential to quantify the uncertainty of actuated responses of the tensegrity structures and to modulate their deformation accurately. In this work, we propose a comprehensive data-driven computational approach to study the uncertainty quantification (UQ) and probability propagation in clustered tensegrity structures, and we have developed a surrogate optimization model to control the flexible structure deformation. An example of clustered tensegrity beam subjected to a clustered actuation is presented to demonstrate the validity of the approach and its potential application. The three main novelties of the data-driven framework are: (1) The proposed model can avoid the difficulty of convergence in nonlinear Finite Element Analysis (FEA), by two machine learning methods, the Gauss Process Regression (GPR) and Neutral Network (NN). (2) A fast real-time prediction on uncertainty propagation can be achieved by the surrogate model, and (3) Optimization of the actuated deformation comes true by using both Sequence Quadratic Programming (SQP) and Bayesian optimization methods. The results have shown that the proposed data-driven computational approach is powerful and can be extended to other UQ models or alternative optimization objectives.

Uncertainty quantification in mechanistic epidemic models via cross-entropy approximate Bayesian computation

This paper proposes a data-driven approximate Bayesian computation framework for parameter estimation and uncertainty quantification of epidemic models, which incorporates two novelties: (i) the identification of the initial conditions by using plausible dynamic states that are compatible with observational data; (ii) learning of an informative prior distribution for the model parameters via the cross-entropy method. The new methodology's effectiveness is illustrated with the aid of actual data from the COVID-19 epidemic in Rio de Janeiro city in Brazil, employing an ordinary differential equation-based model with a generalized SEIR mechanistic structure that includes time-dependent transmission rate, asymptomatics, and hospitalizations. A minimization problem with two cost terms (number of hospitalizations and deaths) is formulated, and twelve parameters are identified. The calibrated model provides a consistent description of the available data, able to extrapolate forecasts over a few weeks, making the proposed methodology very appealing for real-time epidemic modeling.

Towards integration of time-resolved confocal microscopy of a 3D in vitro microfluidic platform with a hybrid multiscale model of tumor angiogenesis

The goal of this study is to calibrate a multiscale model of tumor angiogenesis with time-resolved data to allow for systematic testing of mathematical predictions of vascular sprouting. The multi-scale model consists of an agent-based description of tumor and endothelial cell dynamics coupled to a continuum model of vascular endothelial growth factor concentration. First, we calibrate ordinary differential equation models to time-resolved protein concentration data to estimate the rates of secretion and consumption of vascular endothelial growth factor by endothelial and tumor cells, respectively. These parameters are then input into the multiscale tumor angiogenesis model, and the remaining model parameters are then calibrated to time resolved confocal microscopy images obtained within a 3D vascularized microfluidic platform. The microfluidic platform mimics a functional blood vessel with a surrounding collagen matrix seeded with inflammatory breast cancer cells, which induce tumor angiogenesis. Once the multi-scale model is fully parameterized, we forecast the spatiotemporal distribution of vascular sprouts at future time points and directly compare the predictions to experimentally measured data. We assess the ability of our model to globally recapitulate angiogenic vasculature density, resulting in an average relative calibration error of 17.7% ± 6.3% and an average prediction error of 20.2% ± 4% and 21.7% ± 3.6% using one and four calibrated parameters, respectively. We then assess the model's ability to predict local vessel morphology (individualized vessel structure as opposed to global vascular density), initialized with the first time point and calibrated with two intermediate time points. In this study, we have rigorously calibrated a mechanism-based, multiscale, mathematical model of angiogenic sprouting to multimodal experimental data to make specific, testable predictions.

Modeling nonlocal behavior in epidemics via a reaction-diffusion system incorporating population movement along a network

The outbreak of COVID-19, beginning in 2019 and continuing through the time of writing, has led to renewed interest in the mathematical modeling of infectious disease. Recent works have focused on partial differential equation (PDE) models, particularly reaction-diffusion models, able to describe the progression of an epidemic in both space and time. These studies have shown generally promising results in describing and predicting COVID-19 progression. However, people often travel long distances in short periods of time, leading to nonlocal transmission of the disease. Such contagion dynamics are not well-represented by diffusion alone. In contrast, ordinary differential equation (ODE) models may easily account for this behavior by considering disparate regions as nodes in a network, with the edges defining nonlocal transmission. In this work, we attempt to combine these modeling paradigms via the introduction of a network structure within a reaction-diffusion PDE system. This is achieved through the definition of a population-transfer operator, which couples disjoint and potentially distant geographic regions, facilitating nonlocal population movement between them. We provide analytical results demonstrating that this operator does not disrupt the physical consistency or mathematical well-posedness of the system, and verify these results through numerical experiments. We then use this technique to simulate the COVID-19 epidemic in the Brazilian region of Rio de Janeiro, showcasing its ability to capture important nonlocal behaviors, while maintaining the advantages of a reaction-diffusion model for describing local dynamics.

Delay differential equations for the spatially resolved simulation of epidemics with specific application to COVID-19

In the wake of the 2020 COVID-19 epidemic, much work has been performed on the development of mathematical models for the simulation of the epidemic and of disease models generally. Most works follow the susceptible-infected-removed (SIR) compartmental framework, modeling the epidemic with a system of ordinary differential equations. Alternative formulations using a partial differential equation (PDE) to incorporate both spatial and temporal resolution have also been introduced, with their numerical results showing potentially powerful descriptive and predictive capacity. In the present work, we introduce a new variation to such models by using delay differential equations (DDEs). The dynamics of many infectious diseases, including COVID-19, exhibit delays due to incubation periods and related phenomena. Accordingly, DDE models allow for a natural representation of the problem dynamics, in addition to offering advantages in terms of computational time and modeling, as they eliminate the need for additional, difficult-to-estimate, compartments (such as exposed individuals) to incorporate time delays. In the present work, we introduce a DDE epidemic model in both an ordinary and partial differential equation framework. We present a series of mathematical results assessing the stability of the formulation. We then perform several numerical experiments, validating both the mathematical results and establishing model's ability to reproduce measured data on realistic problems.

Explainable Artificial Intelligence for COVID-19 Diagnosis Through Blood Test Variables

This work proposes an explainable artificial intelligence approach to help diagnose COVID-19 patients based on blood test and pathogen variables. Two glass-box models, logistic regression and explainable boosting machine, and two black-box models, random forest and support vector machine, were used to assess the disease diagnosis. Shapley additive explanations were used to explain predictions for the black-box models, while glass-box models feature importance brought insights into the most relevant features. All global explanations show the eosinophils and leukocytes, white blood cells are among the essential features to help diagnose the COVID-19. Moreover, the best model obtained an AUC of 0.87.

Modeling the Waves of Covid-19

The challenges with modeling the spread of Covid-19 are its power-type growth during the middle stages of the waves with the exponents depending on time, and that the saturation of the waves is mainly due to the protective measures and other restriction mechanisms working in the same direction. The two-phase solution we propose for modeling the total number of detected cases of Covid-19 describes the actual curves for many its waves and in many countries almost with the accuracy of physics laws. Bessel functions play the key role in our approach. The differential equations we obtain are of universal type and can be used in behavioral psychology, invasion ecology (transient processes), etc. The initial transmission rate and the intensity of the restriction mechanisms are the key parameters. This theory provides a convincing explanation of the surprising uniformity of the Covid-19 waves in many places, and can be used for forecasting the epidemic spread. For instance, the early projections for the 3rd wave in the USA appeared sufficiently exact. The Delta-waves (2021) in India, South Africa, UK, and the Netherlands are discussed at the end.

Bayesian inference of heterogeneous epidemic models: Application to COVID-19 spread accounting for long-term care facilities

We propose a high dimensional Bayesian inference framework for learning heterogeneous dynamics of a COVID-19 model, with a specific application to the dynamics and severity of COVID-19 inside and outside long-term care (LTC) facilities. We develop a heterogeneous compartmental model that accounts for the heterogeneity of the time-varying spread and severity of COVID-19 inside and outside LTC facilities, which is characterized by time-dependent stochastic processes and time-independent parameters in ∼ 1500 dimensions after discretization. To infer these parameters, we use reported data on the number of confirmed, hospitalized, and deceased cases with suitable post-processing in both a deterministic inversion approach with appropriate regularization as a first step, followed by Bayesian inversion with proper prior distributions. To address the curse of dimensionality and the ill-posedness of the high-dimensional inference problem, we propose use of a dimension-independent projected Stein variational gradient descent method, and demonstrate the intrinsic low-dimensionality of the inverse problem. We present inference results with quantified uncertainties for both New Jersey and Texas, which experienced different epidemic phases and patterns. Moreover, we also present forecasting and validation results based on the empirical posterior samples of our inference for the future trajectory of COVID-19.

Identifiability and predictability of integer- and fractional-order epidemiological models using physics-informed neural networks

We analyze a plurality of epidemiological models through the lens of physics-informed neural networks (PINNs) that enable us to identify time-dependent parameters and data-driven fractional differential operators. In particular, we consider several variations of the classical susceptible-infectious-removed (SIR) model by introducing more compartments and fractional-order and time-delay models. We report the results for the spread of COVID-19 in New York City, Rhode Island and Michigan states and Italy, by simultaneously inferring the unknown parameters and the unobserved dynamics. For integer-order and time-delay models, we fit the available data by identifying time-dependent parameters, which are represented by neural networks. In contrast, for fractional differential models, we fit the data by determining different time-dependent derivative orders for each compartment, which we represent by neural networks. We investigate the structural and practical identifiability of these unknown functions for different datasets, and quantify the uncertainty associated with neural networks and with control measures in forecasting the pandemic.

News Sentiment Informed Time-series Analyzing AI (SITALA) to curb the spread of COVID-19 in Houston

Coronavirus disease (COVID-19) has evolved into a pandemic with many unknowns. Houston, located in the Harris County of Texas, is becoming the next hotspot of this pandemic. With a severe decline in international and inter-state travel, a model at the county level is needed as opposed to the state or country level. Existing approaches have a few drawbacks. Firstly, the data used is the number of COVID-19 positive cases instead of positivity. The former is a function of the number of tests carried out while the number of tests normalizes the latter. Positivity gives a better picture of the spread of this pandemic as, with time, more tests are being administered. Positivity under 5% has been desired for the reopening of businesses to almost 100% capacity. Secondly, the data used by models like SEIRD (Susceptible, Exposed, Infectious, Recovered, and Deceased) lacks information about the sentiment of people concerning coronavirus. Thirdly, models that make use of social media posts might have too much noise and misinformation. On the other hand, news sentiment can capture long-term effects of hidden variables like public policy, opinions of local doctors, and disobedience of state-wide mandates. The present study introduces a new artificial intelligence (i.e., AI) model, viz., Sentiment Informed Time-series Analyzing AI (SITALA), trained on COVID-19 test positivity data and news sentiment from over 2750 news articles for Harris county. The news sentiment was obtained using IBM Watson Discovery News. SITALA is inspired by Google-Wavenet architecture and makes use of TensorFlow. The mean absolute error for the training dataset of 66 consecutive days is 2.76, and that for the test dataset of 22 consecutive days is 9.6. A cone of uncertainty is provided within which future COVID-19 test positivity has been shown to fall with high accuracy. The model predictions fare better than a published Bayesian-based SEIRD model. The model forecasts that in order to curb the spread of coronavirus in Houston, a sustained negative news sentiment (e.g., death count for COVID-19 will grow at an alarming rate in Houston if mask orders are not followed) will be desirable. Public policymakers may use SITALA to set the tone of the local policies and mandates.

Assessing the Spatio-temporal Spread of COVID-19 via Compartmental Models with Diffusion in Italy, USA, and Brazil

The outbreak of COVID-19 in 2020 has led to a surge in interest in the mathematical modeling of infectious diseases. Such models are usually defined as compartmental models, in which the population under study is divided into compartments based on qualitative characteristics, with different assumptions about the nature and rate of transfer across compartments. Though most commonly formulated as ordinary differential equation models, in which the compartments depend only on time, recent works have also focused on partial differential equation (PDE) models, incorporating the variation of an epidemic in space. Such research on PDE models within a Susceptible, Infected, Exposed, Recovered, and Deceased framework has led to promising results in reproducing COVID-19 contagion dynamics. In this paper, we assess the robustness of this modeling framework by considering different geometries over more extended periods than in other similar studies. We first validate our code by reproducing previously shown results for Lombardy, Italy. We then focus on the U.S. state of Georgia and on the Brazilian state of Rio de Janeiro, one of the most impacted areas in the world. Our results show good agreement with real-world epidemiological data in both time and space for all regions across major areas and across three different continents, suggesting that the modeling approach is both valid and robust.

Spatio-temporal predictive modeling framework for infectious disease spread

A novel predictive modeling framework for the spread of infectious diseases using high-dimensional partial differential equations is developed and implemented. A scalar function representing the infected population is defined on a high-dimensional space and its evolution over all the directions is described by a population balance equation (PBE). New infections are introduced among the susceptible population from a non-quarantined infected population based on their interaction, adherence to distancing norms, hygiene levels and any other societal interventions. Moreover, recovery, death, immunity and all aforementioned parameters are modeled on the high-dimensional space. To epitomize the capabilities and features of the above framework, prognostic estimates of Covid-19 spread using a six-dimensional (time, 2D space, infection severity, duration of infection, and population age) PBE is presented. Further, scenario analysis for different policy interventions and population behavior is presented, throwing more insights into the spatio-temporal spread of infections across duration of disease, infection severity and age of the population. These insights could be used for science-informed policy planning.

Least-Squares Finite Element Method for a Meso-Scale Model of the Spread of COVID-19

This paper investigates numerical properties of a flux-based finite element method for the discretization of a SEIQRD (susceptible-exposed-infected-quarantined-recovered-deceased) model for the spread of COVID-19. The model is largely based on the SEIRD (susceptible-exposed-infected-recovered-deceased) models developed in recent works, with additional extension by a quarantined compartment of the living population and the resulting first-order system of coupled PDEs is solved by a Least-Squares meso-scale method. We incorporate several data on political measures for the containment of the spread gathered during the course of the year 2020 and develop an indicator that influences the predictions calculated by the method. The numerical experiments conducted show a promising accuracy of predictions of the space-time behavior of the virus compared to the real disease spreading data.

A Bayesian Model of COVID-19 Cases Based on the Gompertz Curve

The COVID-19 pandemic has highlighted the need for finding mathematical models to forecast the evolution of the contagious disease and evaluate the success of particular policies in reducing infections. In this work, we perform Bayesian inference for a non-homogeneous Poisson process with an intensity function based on the Gompertz curve. We discuss the prior distribution of the parameter and we generate samples from the posterior distribution by using Markov Chain Monte Carlo (MCMC) methods. Finally, we illustrate our method analyzing real data associated with COVID-19 in a specific region located at the south of Spain.

Adaptive mesh refinement and coarsening for diffusion-reaction epidemiological models

The outbreak of COVID-19 in 2020 has led to a surge in the interest in the mathematical modeling of infectious diseases. Disease transmission may be modeled as compartmental models, in which the population under study is divided into compartments and has assumptions about the nature and time rate of transfer from one compartment to another. Usually, they are composed of a system of ordinary differential equations in time. A class of such models considers the Susceptible, Exposed, Infected, Recovered, and Deceased populations, the SEIRD model. However, these models do not always account for the movement of individuals from one region to another. In this work, we extend the formulation of SEIRD compartmental models to diffusion-reaction systems of partial differential equations to capture the continuous spatio-temporal dynamics of COVID-19. Since the virus spread is not only through diffusion, we introduce a source term to the equation system, representing exposed people who return from travel. We also add the possibility of anisotropic non-homogeneous diffusion. We implement the whole model in libMesh, an open finite element library that provides a framework for multiphysics, considering adaptive mesh refinement and coarsening. Therefore, the model can represent several spatial scales, adapting the resolution to the disease dynamics. We verify our model with standard SEIRD models and show several examples highlighting the present model's new capabilities.

Visualizing the invisible: The effect of asymptomatic transmission on the outbreak dynamics of COVID-19

Understanding the outbreak dynamics of the COVID-19 pandemic has important implications for successful containment and mitigation strategies. Recent studies suggest that the population prevalence of SARS-CoV-2 antibodies, a proxy for the number of asymptomatic cases, could be an order of magnitude larger than expected from the number of reported symptomatic cases. Knowing the precise prevalence and contagiousness of asymptomatic transmission is critical to estimate the overall dimension and pandemic potential of COVID-19. However, at this stage, the effect of the asymptomatic population, its size, and its outbreak dynamics remain largely unknown. Here we use reported symptomatic case data in conjunction with antibody seroprevalence studies, a mathematical epidemiology model, and a Bayesian framework to infer the epidemiological characteristics of COVID-19. Our model computes, in real time, the time-varying contact rate of the outbreak, and projects the temporal evolution and credible intervals of the effective reproduction number and the symptomatic, asymptomatic, and recovered populations. Our study quantifies the sensitivity of the outbreak dynamics of COVID-19 to three parameters: the effective reproduction number, the ratio between the symptomatic and asymptomatic populations, and the infectious periods of both groups. For nine distinct locations, our model estimates the fraction of the population that has been infected and recovered by Jun 15, 2020 to 24.15% (95% CI: 20.48%-28.14%) for Heinsberg (NRW, Germany), 2.40% (95% CI: 2.09%-2.76%) for Ada County (ID, USA), 46.19% (95% CI: 45.81%-46.60%) for New York City (NY, USA), 11.26% (95% CI: 7.21%-16.03%) for Santa Clara County (CA, USA), 3.09% (95% CI: 2.27%-4.03%) for Denmark, 12.35% (95% CI: 10.03%-15.18%) for Geneva Canton (Switzerland), 5.24% (95% CI: 4.84%-5.70%) for the Netherlands, 1.53% (95% CI: 0.76%-2.62%) for Rio Grande do Sul (Brazil), and 5.32% (95% CI: 4.77%-5.93%) for Belgium. Our method traces the initial outbreak date in Santa Clara County back to January 20, 2020 (95% CI: December 29, 2019-February 13, 2020). Our results could significantly change our understanding and management of the COVID-19 pandemic: A large asymptomatic population will make isolation, containment, and tracing of individual cases challenging. Instead, managing community transmission through increasing population awareness, promoting physical distancing, and encouraging behavioral changes could become more relevant.

Data-driven modeling of COVID-19-Lessons learned

Understanding the outbreak dynamics of COVID-19 through the lens of mathematical models is an elusive but significant goal. Within only half a year, the COVID-19 pandemic has resulted in more than 19 million reported cases across 188 countries with more than 700,000 deaths worldwide. Unlike any other disease in history, COVID-19 has generated an unprecedented volume of data, well documented, continuously updated, and broadly available to the general public. Yet, the precise role of mathematical modeling in providing quantitative insight into the COVID-19 pandemic remains a topic of ongoing debate. Here we discuss the lessons learned from six month of modeling COVID-19. We highlight the early success of classical models for infectious diseases and show why these models fail to predict the current outbreak dynamics of COVID-19. We illustrate how data-driven modeling can integrate classical epidemiology modeling and machine learning to infer critical disease parameters-in real time-from reported case data to make informed predictions and guide political decision making. We critically discuss questions that these models can and cannot answer and showcase controversial decisions around the early outbreak dynamics, outbreak control, and exit strategies. We anticipate that this summary will stimulate discussion within the modeling community and help provide guidelines for robust mathematical models to understand and manage the COVID-19 pandemic. EML webinar speakers, videos, and overviews are updated at https://imechanica.org/node/24098.

Is it safe to lift COVID-19 travel bans? The Newfoundland story

A key strategy to prevent a local outbreak during the COVID-19 pandemic is to restrict incoming travel. Once a region has successfully contained the disease, it becomes critical to decide when and how to reopen the borders. Here we explore the impact of border reopening for the example of Newfoundland and Labrador, a Canadian province that has enjoyed no new cases since late April, 2020. We combine a network epidemiology model with machine learning to infer parameters and predict the COVID-19 dynamics upon partial and total airport reopening, with perfect and imperfect quarantine conditions. Our study suggests that upon full reopening, every other day, a new COVID-19 case would enter the province. Under the current conditions, banning air travel from outside Canada is more efficient in managing the pandemic than fully reopening and quarantining 95% of the incoming population. Our study provides quantitative insights of the efficacy of travel restrictions and can inform political decision making in the controversy of reopening.

Visualizing the invisible: The effect of asymptomatic transmission on the outbreak dynamics of COVID-19

Understanding the outbreak dynamics of the COVID-19 pandemic has important implications for successful containment and mitigation strategies. Recent studies suggest that the population prevalence of SARS-CoV-2 antibodies, a proxy for the number of asymptomatic cases, could be an order of magnitude larger than expected from the number of reported symptomatic cases. Knowing the precise prevalence and contagiousness of asymptomatic transmission is critical to estimate the overall dimension and pandemic potential of COVID-19. However, at this stage, the effect of the asymptomatic population, its size, and its outbreak dynamics remain largely unknown. Here we use reported symptomatic case data in conjunction with antibody seroprevalence studies, a mathematical epidemiology model, and a Bayesian framework to infer the epidemiological characteristics of COVID-19. Our model computes, in real time, the time-varying contact rate of the outbreak, and projects the temporal evolution and credible intervals of the effective reproduction number and the symptomatic, asymptomatic, and recovered populations. Our study quantifies the sensitivity of the outbreak dynamics of COVID-19 to three parameters: the effective reproduction number, the ratio between the symptomatic and asymptomatic populations, and the infectious periods of both groups For nine distinct locations, our model estimates the fraction of the population that has been infected and recovered by Jun 15, 2020 to 24.15% (95% CI: 20.48%-28.14%) for Heinsberg (NRW, Germany), 2.40% (95% CI: 2.09%-2.76%) for Ada County (ID, USA), 46.19% (95% CI: 45.81%-46.60%) for New York City (NY, USA), 11.26% (95% CI: 7.21%-16.03%) for Santa Clara County (CA, USA), 3.09% (95% CI: 2.27%-4.03%) for Denmark, 12.35% (95% CI: 10.03%-15.18%) for Geneva Canton (Switzerland), 5.24% (95% CI: 4.84%-5.70%) for the Netherlands, 1.53% (95% CI: 0.76%-2.62%) for Rio Grande do Sul (Brazil), and 5.32% (95% CI: 4.77%-5.93%) for Belgium. Our method traces the initial outbreak date in Santa Clara County back to January 20, 2020 (95% CI: December 29, 2019 - February 13, 2020). Our results could significantly change our understanding and management of the COVID-19 pandemic: A large asymptomatic population will make isolation, containment, and tracing of individual cases challenging. Instead, managing community transmission through increasing population awareness, promoting physical distancing, and encouraging behavioral changes could become more relevant.

Is it safe to lift COVID-19 travel bans? The Newfoundland story

A key strategy to prevent a local outbreak during the COVID-19 pandemic is to restrict incoming travel. Once a region has successfully contained the disease, it becomes critical to decide when and how to reopen the borders. Here we explore the impact of border reopening for the example of Newfoundland and Labrador, a Canadian province that has enjoyed no new cases since late April, 2020. We combine a network epidemiology model with machine learning to infer parameters and predict the COVID-19 dynamics upon partial and total airport reopening, with perfect and imperfect quarantine conditions. Our study suggests that upon full reopening, every other day, a new COVID-19 case would enter the province. Under the current conditions, banning air travel from outside Canada is more efficient in managing the pandemic than fully reopening and quarantining 95% of the incoming population. Our study provides quantitative insights of the efficacy of travel restrictions and can inform political decision making in the controversy of reopening.