Recursive state and parameter estimation of COVID-19 circulating variants dynamics

COVID-19 pandemic response with non-pharmaceutical interventions is an intrinsic control problem. Governments weigh social distancing policies to avoid overload in the health system without significant economic impact. The mutability of the SARS-CoV-2 virus, vaccination coverage, and mobility restriction measures change epidemic dynamics over time. A model-based control strategy requires reliable predictions to be efficient on a long-term basis. In this paper, a SEIR-based model is proposed considering dynamic feedback estimation. State and parameter estimations are performed on state estimators using augmented states. Three methods were implemented: constrained extended Kalman filter (CEKF), CEKF and smoother (CEKF & S), and moving horizon estimator (MHE). The parameters estimation was based on vaccine efficacy studies regarding transmissibility, severity of the disease, and lethality. Social distancing was assumed as a measured disturbance calculated using Google mobility data. Data from six federative units from Brazil were used to evaluate the proposed strategy. State and parameter estimations were performed from 1 October 2020 to 1 July 2021, during which Zeta and Gamma variants emerged. Simulation results showed that lethality increased between 11 and 30% for Zeta mutations and between 44 and 107% for Gamma mutations. In addition, transmissibility increased between 10 and 37% for the Zeta variant and between 43 and 119% for the Gamma variant. Furthermore, parameter estimation indicated temporal underreporting changes in hospitalized and deceased individuals. Overall, the estimation strategy showed to be suitable for dynamic feedback as simulation results presented an efficient detection and dynamic characterization of circulating variants.