Time series modelling for wastewater-based epidemiology of COVID-19: A nationwide study in 40 wastewater treatment plants of Belgium, February 2021 to June 2022.

BACKGROUND: Wastewater-based epidemiology (WBE) has been implemented to monitor surges of COVID-19. Yet, multiple factors impede the usefulness of WBE and quantitative adjustment may be required. AIM: We aimed to model the relationship between WBE data and incident COVID-19 cases, while adjusting for confounders and autocorrelation. METHODS: This nationwide WBE study includes data from 40 wastewater treatment plants (WWTPs) in Belgium (02/2021-06/2022). We applied ARIMA-based modelling to assess the effect of daily flow rate, pepper mild mottle virus (PMMoV) concentration, a measure of human faeces in wastewater, and variants (alpha, delta, and omicron strains) on SARS-CoV-2 RNA levels in wastewater. Secondly, adjusted WBE metrics at different lag times were used to predict incident COVID-19 cases. Model selection was based on AICc minimization. RESULTS: In 33/40 WWTPs, RNA levels were best explained by incident cases, flow rate, and PMMoV. Flow rate and PMMoV were associated with -13.0 % (95 % prediction interval: -26.1 to +0.2 %) and +13.0 % (95 % prediction interval: +5.1 to +21.0 %) change in RNA levels per SD increase, respectively. In 38/40 WWTPs, variants did not explain variability in RNA levels independent of cases. Furthermore, our study shows that RNA levels can lead incident cases by at least one week in 15/40 WWTPs. The median population size of leading WWTPs was 85.1 % larger than that of non­leading WWTPs. In 17/40 WWTPs, however, RNA levels did not lead or explain incident cases in addition to autocorrelation. CONCLUSION: This study provides quantitative insights into key determinants of WBE, including the effects of wastewater flow rate, PMMoV, and variants. Substantial inter-WWTP variability was observed in terms of explaining incident cases. These findings are of practical importance to WBE practitioners and show that the early-warning potential of WBE is WWTP-specific and needs validation.

SARS-CoV-2 Surveillance in Belgian Wastewaters.

Wastewater-based surveillance was conducted by the national public health authority to monitor SARS-CoV-2 circulation in the Belgian population. Over 5 million inhabitants representing 45% of the Belgian population were monitored throughout 42 wastewater treatment plants for 15 months comprising three major virus waves. During the entire period, a high correlation was observed between the daily new COVID-19 cases and the SARS-CoV-2 concentration in wastewater corrected for rain impact and covered population size. Three alerting indicators were included in the weekly epidemiological assessment: High Circulation, Fast Increase, and Increasing Trend. These indicators were computed on normalized concentrations per individual treatment plant to allow for a comparison with a reference period as well as between analyses performed by distinct laboratories. When the indicators were not corrected for rain impact, rainy events caused an underestimation of the indicators. Despite this negative impact, the indicators permitted us to effectively monitor the evolution of the fourth virus wave and were considered complementary and valuable information to conventional epidemiological indicators in the weekly wastewater reports communicated to the National Risk Assessment Group.