skip to main content


Title: The doubling time analysis for modified infectious disease Richards model with applications to COVID-19 pandemic

In the absence of reliable information about transmission mechanisms for emerging infectious diseases, simple phenomenological models could provide a starting point to assess the potential outcomes of unfolding public health emergencies, particularly when the epidemiological characteristics of the disease are poorly understood or subject to substantial uncertainty. In this study, we employ the modified Richards model to analyze the growth of an epidemic in terms of 1) the number of times cumulative cases double until the epidemic peaks and 2) the rate at which the intervals between consecutive doubling times increase during the early ascending stage of the outbreak. Our theoretical analysis of doubling times is combined with rigorous numerical simulations and uncertainty quantification using synthetic and real data for COVID-19 pandemic. The doubling-time approach allows to employ early epidemic data to differentiate between the most dangerous threats, which double in size many times over the intervals that are nearly invariant, and the least transmissible diseases, which double in size only a few times with doubling periods rapidly growing.

 
more » « less
Award ID(s):
2011622
NSF-PAR ID:
10351008
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Mathematical Biosciences and Engineering
Volume:
19
Issue:
3
ISSN:
1551-0018
Page Range / eLocation ID:
3242 to 3268
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Simple mathematical tools are needed to quantify the threat posed by emerging and re-emerging infectious disease outbreaks using minimal data capturing the outbreak trajectory. Here we use mathematical analysis, simulation and COVID-19 epidemic data to demonstrate a novel approach to numerically and mathematically characterize the rate at which the doubling time of an epidemic is changing over time. For this purpose, we analyze the dynamics of epidemic doubling times during the initial epidemic stage, defined as the sequence of times at which the cumulative incidence doubles. We introduce new methodology to characterize epidemic threats by analyzing the evolution of epidemics as a function of (1) the number of times the epidemic doubles until the epidemic peak is reached and (2) the rate at which the doubling times increase. In our doubling-time approach, the most dangerous epidemic threats double in size many times and the doubling times change at a relatively low rate (e.g., doubling times remain nearly invariant) whereas the least transmissible threats double in size only a few times and the doubling times rapidly increases in the period of emergence. We derive analytical formulas and test and illustrate our methodology using synthetic and COVID-19 epidemic data. Our mathematical analysis demonstrates that the series of epidemic doubling times increase approximately according to an exponential function with a rate that quantifies the rate of change of the doubling times. Our analytic results are in excellent agreement with numerical results. Our methodology offers a simple and intuitive approach that relies on minimal outbreak trajectory data to characterize the threat posed by emerging and re-emerging infectious diseases. 
    more » « less
  2. Abstract Background

    No versatile web app exists that allows epidemiologists and managers around the world to comprehensively analyze the impacts of COVID-19 mitigation. Thehttp://covid-webapp.numerusinc.com/web app presented here fills this gap.

    Methods

    Our web app uses a model that explicitly identifies susceptible, contact, latent, asymptomatic, symptomatic and recovered classes of individuals, and a parallel set of response classes, subject to lower pathogen-contact rates. The user inputs a CSV file of incidence and, if of interest, mortality rate data. A default set of parameters is available that can be overwritten through input or online entry, and a user-selected subset of these can be fitted to the model using maximum-likelihood estimation (MLE). Model fitting and forecasting intervals are specifiable and changes to parameters allow counterfactual and forecasting scenarios. Confidence or credible intervals can be generated using stochastic simulations, based on MLE values, or on an inputted CSV file containing Markov chain Monte Carlo (MCMC) estimates of one or more parameters.

    Results

    We illustrate the use of our web app in extracting social distancing, social relaxation, surveillance or virulence switching functions (i.e., time varying drivers) from the incidence and mortality rates of COVID-19 epidemics in Israel, South Africa, and England. The Israeli outbreak exhibits four distinct phases: initial outbreak, social distancing, social relaxation, and a second wave mitigation phase. An MCMC projection of this latter phase suggests the Israeli epidemic will continue to produce into late November an average of around 1500 new case per day, unless the population practices social-relaxation measures at least 5-fold below the level in August, which itself is 4-fold below the level at the start of July. Our analysis of the relatively late South African outbreak that became the world’s fifth largest COVID-19 epidemic in July revealed that the decline through late July and early August was characterised by a social distancing driver operating at more than twice the per-capita applicable-disease-class (pc-adc) rate of the social relaxation driver. Our analysis of the relatively early English outbreak, identified a more than 2-fold improvement in surveillance over the course of the epidemic. It also identified a pc-adc social distancing rate in early August that, though nearly four times the pc-adc social relaxation rate, appeared to barely contain a second wave that would break out if social distancing was further relaxed.

    Conclusion

    Our web app provides policy makers and health officers who have no epidemiological modelling or computer coding expertise with an invaluable tool for assessing the impacts of different outbreak mitigation policies and measures. This includes an ability to generate an epidemic-suppression or curve-flattening index that measures the intensity with which behavioural responses suppress or flatten the epidemic curve in the region under consideration.

     
    more » « less
  3. Abstract Background

    Vector-borne diseases (VBDs) are important contributors to the global burden of infectious diseases due to their epidemic potential, which can result in significant population and economic impacts. Oropouche fever, caused by Oropouche virus (OROV), is an understudied zoonotic VBD febrile illness reported in Central and South America. The epidemic potential and areas of likely OROV spread remain unexplored, limiting capacities to improve epidemiological surveillance.

    Methods

    To better understand the capacity for spread of OROV, we developed spatial epidemiology models using human outbreaks as OROV transmission-locality data, coupled with high-resolution satellite-derived vegetation phenology. Data were integrated using hypervolume modeling to infer likely areas of OROV transmission and emergence across the Americas.

    Results

    Models based on one-support vector machine hypervolumes consistently predicted risk areas for OROV transmission across the tropics of Latin America despite the inclusion of different parameters such as different study areas and environmental predictors. Models estimate that up to 5 million people are at risk of exposure to OROV. Nevertheless, the limited epidemiological data available generates uncertainty in projections. For example, some outbreaks have occurred under climatic conditions outside those where most transmission events occur. The distribution models also revealed that landscape variation, expressed as vegetation loss, is linked to OROV outbreaks.

    Conclusions

    Hotspots of OROV transmission risk were detected along the tropics of South America. Vegetation loss might be a driver of Oropouche fever emergence. Modeling based on hypervolumes in spatial epidemiology might be considered an exploratory tool for analyzing data-limited emerging infectious diseases for which little understanding exists on their sylvatic cycles. OROV transmission risk maps can be used to improve surveillance, investigate OROV ecology and epidemiology, and inform early detection.

     
    more » « less
  4. Abstract

    Infectious diseases are a major threat for biodiversity conservation and can exert strong influence on wildlife population dynamics. Understanding the mechanisms driving infection rates and epidemic outcomes requires empirical data on the evolutionary trajectory of pathogens and host selective processes. Phylodynamics is a robust framework to understand the interaction of pathogen evolutionary processes with epidemiological dynamics, providing a powerful tool to evaluate disease control strategies. Tasmanian devils have been threatened by a fatal transmissible cancer, devil facial tumour disease (DFTD), for more than two decades. Here we employ a phylodynamic approach using tumour mitochondrial genomes to assess the role of tumour genetic diversity in epidemiological and population dynamics in a devil population subject to 12 years of intensive monitoring, since the beginning of the epidemic outbreak. DFTD molecular clock estimates of disease introduction mirrored observed estimates in the field, and DFTD genetic diversity was positively correlated with estimates of devil population size. However, prevalence and force of infection were the lowest when devil population size and tumour genetic diversity was the highest. This could be due to either differential virulence or transmissibility in tumour lineages or the development of host defence strategies against infection. Our results support the view that evolutionary processes and epidemiological trade‐offs can drive host‐pathogen coexistence, even when disease‐induced mortality is extremely high. We highlight the importance of integrating pathogen and population evolutionary interactions to better understand long‐term epidemic dynamics and evaluating disease control strategies.

     
    more » « less
  5. Abstract. Methane emissions from boreal and arctic wetlands, lakes, and rivers areexpected to increase in response to warming and associated permafrost thaw.However, the lack of appropriate land cover datasets for scalingfield-measured methane emissions to circumpolar scales has contributed to alarge uncertainty for our understanding of present-day and future methaneemissions. Here we present the Boreal–Arctic Wetland and Lake Dataset(BAWLD), a land cover dataset based on an expert assessment, extrapolatedusing random forest modelling from available spatial datasets of climate,topography, soils, permafrost conditions, vegetation, wetlands, and surfacewater extents and dynamics. In BAWLD, we estimate the fractional coverage offive wetland, seven lake, and three river classes within 0.5 × 0.5∘ grid cells that cover the northern boreal and tundra biomes(17 % of the global land surface). Land cover classes were defined usingcriteria that ensured distinct methane emissions among classes, as indicatedby a co-developed comprehensive dataset of methane flux observations. InBAWLD, wetlands occupied 3.2 × 106 km2 (14 % of domain)with a 95 % confidence interval between 2.8 and 3.8 × 106 km2. Bog, fen, and permafrost bog were the most abundant wetlandclasses, covering ∼ 28 % each of the total wetland area,while the highest-methane-emitting marsh and tundra wetland classes occupied5 % and 12 %, respectively. Lakes, defined to include all lentic open-waterecosystems regardless of size, covered 1.4 × 106 km2(6 % of domain). Low-methane-emitting large lakes (>10 km2) and glacial lakes jointly represented 78 % of the total lakearea, while high-emitting peatland and yedoma lakes covered 18 % and 4 %,respectively. Small (<0.1 km2) glacial, peatland, and yedomalakes combined covered 17 % of the total lake area but contributeddisproportionally to the overall spatial uncertainty in lake area with a95 % confidence interval between 0.15 and 0.38 × 106 km2. Rivers and streams were estimated to cover 0.12  × 106 km2 (0.5 % of domain), of which 8 % was associated withhigh-methane-emitting headwaters that drain organic-rich landscapes.Distinct combinations of spatially co-occurring wetland and lake classeswere identified across the BAWLD domain, allowing for the mapping of“wetscapes” that have characteristic methane emission magnitudes andsensitivities to climate change at regional scales. With BAWLD, we provide adataset which avoids double-accounting of wetland, lake, and river extentsand which includes confidence intervals for each land cover class. As such,BAWLD will be suitable for many hydrological and biogeochemical modellingand upscaling efforts for the northern boreal and arctic region, inparticular those aimed at improving assessments of current and futuremethane emissions. Data are freely available athttps://doi.org/10.18739/A2C824F9X (Olefeldt et al., 2021). 
    more » « less