Total: 1
Wastewater surveillance provides population-level evidence of infectious disease activity. As surveillance expands across multiple sites, repeated measurements must be analyzed jointly while accounting for differences among sites and the public health regions they represent. We develop a Bayesian spatio-temporal framework that decomposes wastewater measurements into a shared epidemic curve, spatially correlated regional deviations, persistent site effects, and within-site temporal dependence. This hierarchical decomposition accommodates unequal observation windows and sampling schedules, allowing regional trajectories to be compared while preserving site-level heterogeneity. A direct marginal spatial covariance specification further separates the magnitude of regional variation from the strength of association between neighboring regions, providing distinct interpretations of these two sources of spatial variation. We also introduce region-specific time shifts to distinguish regional leads and lags in epidemic timing from differences in signal levels and local fluctuations. Simulation studies demonstrate the recovery of these components and show that omitting true time shifts can lead to increased estimates of residual temporal variation. Application to Ontario SARS-CoV-2 wastewater surveillance data reveals persistent site heterogeneity, within-site temporal dependence, and spatial association among neighboring regions. Regional timing differences are more pronounced during the first year of Omicron than in the XBB-and-later period.