The primary objective of this study was to identify a universal wastewater biomarker for population normalization for SARS-CoV-2 wastewater-based epidemiology (WBE). A total of 2,624 wastewater samples (41 weeks) were collected weekly during May 2021- April 2022 from 64 wastewater facilities across Missouri, U.S. Three wastewater biomarkers, caffeine and its metabolite, paraxanthine, and pepper mild mottle virus (PMMoV), were compared for the population normalization effectiveness for wastewater SARS-CoV-2 surveillance. Paraxanthine had the lowest temporal variation and strongest relationship between population compared to caffeine and PMMoV. This result was confirmed by data from ten different Wisconsin WWTPs with gradients in population sizes, indicating paraxanthine is a promising biomarker of the real-time population across a large geographical region. The estimated real-time population was directly compared against the population patterns with human movement mobility data. Of the three biomarkers, population normalization by paraxanthine significantly strengthened the relationship between wastewater SARS-CoV-2 viral load and COVID-19 incidence rate the most (40 out of 61 sewersheds). Caffeine could be a promising population biomarker for regions where no significant exogenous caffeine sources (e.g., discharges from food industries) exist. In contrast, PMMoV showed the highest variability over time, and therefore reduced the strength of the relationship between sewage SARS-CoV-2 viral load and the COVID-19 incidence rate, as compared to wastewater data without population normalization and the population normalized by either recent Census population or the population estimated based on the number of residential connections and average household size for that municipality from the Census. Overall, the findings of this long-term surveillance study concluded that the paraxanthine has the best performance as a biomarker for population normalization for SARS-CoV-2 wastewater-based epidemiology.