Testing methods to estimate population size for wastewater treatment plants using census data: Implications for wastewater-based epidemiology

In wastewater-based epidemiology (WBE), wastewater loads are commonly reported as a per capita value. Census population counts are often used to obtain a population size to normalise wastewater loads. However, the methods used to calculate the population size of wastewater treatment plants (WWTPs) from census data are rarely reported in the WBE literature. This is problematic because the geographical extents of wastewater catchments and census area units rarely align perfectly with each other and exist at different spatial scales. This complicates efforts to estimate the number of people serviced by WWTPs in these census area units. This study compared four geospatial methods to combine wastewater catchment areas and census area units to calculate the census population size of wastewater treatment plants. These methods were applied nationally to WWTPs across New Zealand. Population estimates varied by up to 73 % between the methods, which could skew comparisons of per capita wastewater loads between sites. Variability in population estimates (relative standard deviation, RSD) was significantly higher in smaller catchments (r(s) = -0.727, P < .001), highlighting the importance of method selection in smaller sites. Census population estimates were broadly similar to those provided by wastewater operators, but significant variation was observed for some sites (ranging from 42 % lower to 78 % higher, RSD = 262 %). We present a widely applicable method to calculate population size from census, which involves disaggregating census area units by individual properties. The results reinforce the need for transparent reporting to maintain confidence in the comparison of WBE across sites and studies.