The modified GWR models had higher R2 and reflected the actual spatial features.
A manual variable excluding-selecting method is explored to avoid multicollinearity.
Influences of the dominant indicator on water quality varied with space and seasons.
Protection policies need consider site-specific conditions and seasonal variations.