skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Forecasting Doubling Time of SARS-CoV-2 using Hawkes and SQUIDER Models
The rate of spread of an emerging epidemic is frequently characterized via the doubling time, which is the time it takes for the number of cases to double. This paper explores different ways to estimate doubling time, and investigates the estimation of doubling time in relationship to parameters in the HawkesN model and the SQUIDER (Susceptible, Quarantine, Undetected Infected, Infected, Dead, Exposed, Recovered) model. We observe an approximately exponential relationship between the productivity parameter κ in the HawkesN model and doubling time. We also evaluate the performance of the models in forecasting doubling times and compare to empirical doubling times using daily reported statewide totals for SARS-CoV-2 infections in California, and find that the HawkesN model forecasts doubling times more accurately, with 3.6% smaller root mean squared errors in Spring 2020, 79.4% smaller root mean squared errors in Autumn 2020, and 5.4% smaller root mean squared errors in Summer 2021. The HawkesN and SQUIDER models appear to forecast daily rate doubling times accurately at most times, though the SQUIDER forecasts of daily rate doubling times are far more volatile and thus occasionally have much larger errors, particularly in Fall 2020.  more » « less
Award ID(s):
2124433
PAR ID:
10552357
Author(s) / Creator(s):
; ;
Publisher / Repository:
European Society of Medicine
Date Published:
Journal Name:
Medical Research Archives
Volume:
12
Issue:
3
ISSN:
2375-1916
Page Range / eLocation ID:
1-19
Subject(s) / Keyword(s):
Keywords: Contagious diseases, epidemics, Hawkes model, Point process, SARS-Cov-2, Self-exciting.
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Heat storage change (HSC) is a crucial component of lake's thermal energy budget. Conventional temperature profile based models of HSC require location specific parameters such as lakebed topography. Based on the half‐order time‐derivative formula of heat fluxes, an analytical model was formulated for estimating HSC from water surface temperature and solar radiation without using geography dependent parameters. The proposed model was tested against field measurements at Poyang Lake, a shallow inland lake, which has pronounced seasonal variations in water level and lake area. Our analysis indicates that the model accurately simulates diurnal HSC with a coefficient of determination of 0.94 and a root mean squared error (RMSE) of 77.5 ± 21.6 Wm−2for the study period. Larger nighttime RMSE (75.0 ± 26.8 Wm−2) than the daytime value (55.1 ± 19.7 W m−2) is attributable to larger measurement errors of nighttime turbulent fluxes. The estimation of HSC independent of temperature profile and lake‐specific parameters by the proposed model facilitates remote sensing monitoring the HSC of global water bodies. 
    more » « less
  2. Abstract Water quality and freshwater ecosystems are affected by river discharge and temperature. Models are frequently used to estimate river temperature on large spatial and temporal scales due to limited observations of discharge and temperature. In this study, we use physically based river routing and temperature models to simulate daily discharge and river temperature for rivers in 138 basins in Alaska, including the entire Yukon River basin, from 1990–2021. The river temperature model was optimized for ice free months using a surrogate‐based model optimization method, improving model performance at uncalibrated river gages. A common statistical model relating local air and water temperature was used as a benchmark. The physically based river temperature model exhibited superior performance compared to the benchmark statistical model after optimization, suggesting river temperature model optimization could become more routine. The river temperature model demonstrated high sensitivity to air temperature and model parameterization, and lower sensitivity to discharge. Validation of the models showed a Kling‐Gupta Efficiency of 0.46 for daily river discharge and a root mean square error of 2.04°C for daily river temperature, improving on the non‐optimized physical model and the benchmark statistical model, which had root mean square errors of 3.24 and 2.97°C, respectively. The simulation shows that rivers in northern Alaska have higher maximum summer temperatures and more variability than rivers in the Central and Southern regions. Furthermore, this framework can be readily adapted for use across models and regions. 
    more » « less
  3. null (Ed.)
    This paper describes an implementation of the Combined Hybrid-Parallel Prediction (CHyPP) approach of Wikner et al. (2020) on a low-resolution atmospheric global circulation model (AGCM). This approach combines a physics-based numerical model of a dynamical system (e.g., the atmosphere) with a computationally efficient type of machine learning (ML) called reservoir computing (RC) to construct a hybrid model. This hybrid atmospheric model produces more accurate forecasts of most atmospheric state variables than the host AGCM for the first 7-8 forecast days, and for even longer times for the temperature and humidity near the earth's surface. It also produces more accurate forecasts than a purely ML-based model, or a model that combines linear regression, rather than ML, with the AGCM. The potential of the approach for climate research is demonstrated by a 10-year long hybrid model simulation of the atmospheric general circulation, which shows that the hybrid model can simulate the general circulation with substantially smaller systematic errors and more realistic variability than the host AGCM. 
    more » « less
  4. Abstract. Near-surface air temperature (Ta) is a key variable in global climatestudies. A global gridded dataset of daily maximum and minimum Ta (Tmax⁡ and Tmin⁡) is particularly valuable and critically needed inthe scientific and policy communities but is still not available. In this paper, we developed a global dataset of daily Tmax⁡ and Tmin⁡at 1 km resolution over land across 50∘ S–79∘ N from 2003 to 2020 through the combined use of ground-station-basedTa measurements and satellite observations (i.e., digital elevation model and land surface temperature) via a state-of-the-artstatistical method named Spatially Varying Coefficient Models with Sign Preservation (SVCM-SP). The root mean square errors in our estimates rangedfrom 1.20 to 2.44 ∘C for Tmax⁡ and 1.69 to 2.39 ∘C for Tmin⁡. We found that the accuracies were affectedprimarily by land cover types, elevation ranges, and climate backgrounds. Our dataset correctly represents a negative relationship betweenTa and elevation and a positive relationship between Ta and land surface temperature; it captured spatial and temporalpatterns of Ta realistically. This global 1 km gridded daily Tmax⁡ and Tmin⁡ dataset is the first of its kind, and weexpect it to be of great value to global studies such as the urban heat island phenomenon, hydrological modeling, and epidemic forecasting. The data havebeen published by Iowa State University at https://doi.org/10.25380/iastate.c.6005185 (Zhang and Zhou, 2022). 
    more » « less
  5. Empirical-statistical downscaling (ESD) can be a computationally advantageous alternative to dynamical downscaling in representing a high-resolution regional climate. Two distinct strategies of ESD are employed here to reconstruct near-surface winds in a region of rugged terrain. ESD is used to reconstruct the innermost grid of a multiply nested mesoscale model framework for regional climate downscaling. An analog ensemble (AnEn) and a convolutional neural network (CNN) are compared in their ability to represent near-surface winds in the innermost grid in lieu of dynamical downscaling. Downscaling for a 30 year climatology of 10 m April winds is performed for southern MO, USA. Five years of training suffices for producing low mean absolute error and bias for both ESD techniques. However, root-mean-squared error is not significantly reduced by either scheme. In the case of the AnEn, this is due to a minority of cases not producing a satisfactory representation of high-resolution wind, accentuating the root-mean-squared error in spite of a small mean absolute error. Homogeneous comparison shows that the AnEn produces smaller errors than the CNN. Though further tuning may improve results, the ESD techniques considered here show that they can produce a reliable, computationally inexpensive method for reconstructing high-resolution 10 m winds over complex terrain. 
    more » « less