The integrated monitoring system of water quality is eminently reliant on water quality trend data. This study aims to obtain water quality patterns related to land-use change over a periodic observation in the Opak sub-watershed, Indonesia, both from a seasonal and spatial point of view. Landsat image data from 2013 to 2020 and water quality data comprising 25 parameters were compiled and analyzed. This study observed that land use remarkably correlated to water quality, especially the building area representing the dense population and various anthropogenic activities, to pollute the water sources. Three types of pollutant sources were identified using principal component analysis (PCA), including domestic, industrial, and agricultural activities, which all influenced the variance in river water quality. The use of spatiotemporal-based and multivariate analysis was to interpret water quality trend data, which can help the stakeholders to monitor pollution and take control in the Opak sub-watershed. The results investigated 17 out of 25 water quality parameters, which showed an increasing trend from upstream to downstream during the observation time. The concentration of biological oxygen demand over five days (BOD5), chemical oxygen demand (COD), nitrite, sulfide, phenol, phosphate, oil and grease, lead, Escherichia coli (E. coli), and total coli, surpassed the water quality standard through spatial analysis.