Maritsa River is one of the largest rivers flowing on Bulgarian territory. The quality of its waters is of substantial importance for irrigation, industrial, recreation and domestic use. Besides, part of the river is flowing on Turkish territory and the control and management of the Maritsa catchment is of mutual interst for the neighboring countires. Thus, performing interpretation and modeling of the river water quality is a major environmetric problem. Two multivariate statstical methods (Cluster analysis/CA/and Principal components analysis/PCA/) were applied for model assessment of the water quality of Maritsa River on Bulgarian territory. The study used long-term monitoring data from 21 sampling sites characterized by 8 surface water quality indicators. The application of CA to the indicators results in 3 significant clusters showing the impact of biological, anthropogenic and eutrophication sources. For further assessment of the monitoring data, PCA was implemented, which identified, again,three latent factors confirming, in principle, the clustering output. The latent factors were conditionally named "biologic", "anthropogenic" and "eutrophication" source. Their identification coinside correctly to the location of real pollution sources along the Maritsa River catchment. The linkage of the sampling sites along the river flow by CA identified four special patterns separated by specific tracers levels: biological and anthropogenic major impact for pattern 1, euthrophication major impact for pattern 2, background levels for pattern 3 and eutrophication and agricultural major impact for pattern 4. The apportionment models of the pollution determined the contribution of each one of identified pollution factors to the total concentration of each one of the water quality parameters. Thus, a better risk management of the surface water quality is achieved both on local and national level.