Discriminant analysis (DA) and principal component analysis (PCA), as multivariate statistical techniques, are used to interpret large complex water quality data and assess their temporal and spatial variation in the basin of the Zohreh river. In this study, data sets of 16 water quality parameters collected from 1966 to 2013) in 4 stations (1554 observations for each parameter) were analyzed. PCA for data sets of Kheirabad, Poleflour, Chambostan and Dehmolla stations resulted in 4, 4, 4, and 3 latent factors accounting for 88. 985%, 93. 828%, 88. 648%, and 88. 68% of the total variance in water quality parameters, respectively. It is indicated that total dissolved solids (TDS), electrical conductivity (EC), chlorides (Cl− ), sodium (Na), sodium absorption ratio (SAR), and %Na were responsible for water quality variations which are mainly related to natural and anthropogenic pollution sources including climate effects, gypsum, and salt crystals in the supratidal of Zohreh river delta, fault zones of CHAMSHIR I and II, drainage of sugarcane fields, and domestic and industrial wastewaters discharge into the river. DA reduced the data set to only seven parameters (discharge, temperature, electrical conductivity, HCO3-, Cl-, %Na, and T-Hardness), affording more than 58. 5% correct assignations in temporal evaluations and describing responsible parameters for large variations in the quality of the Zohreh river.