Computed values of the Bartlett test of sphericity which were greater than 0.5 and the KMO results with sig. <0.001 suggest that the data set was suitable for FA. Table 3 shows the Varimax rotated factor matrix consisting of the component factors, loadings of variables on each factor, and percentage of data variance explained by each factor. Four extracted factors described 76.20% of the data set variance. Factor 1 showed the highest factor loadings on Sr²⁺ (0.928) and Mg²⁺ (0.909), high factor loading on Ca²⁺ (0.749) and EC (0.696), and it explained 25.67% of the total variance (Figure 4(a)). Factor 2, which explained 18.76% of the total variance, exhibited high positive loadings on (0.825) and UV254 (0.824) and negative loading on F⁻ (−0.594). Factor 3 accounting for a similar amount of variance to Factor 2, showed high and positive loadings on both Na⁺ (0.791) and (0.761) (Figure 4(b)). Factor 4 explained 15.75% of the total variance, with major positive loadings on Cl⁻ (0.717) and K⁺ (0.686). Communalities were the results of sums of factors' squares, which meant that the proportion of variance can be explained by the extracted factors, and the rest is defined as a specific factor; both contribute to the variance of variables. Calculated results of communalities suggested that over than 65% of the variances for all variables except for and K⁺ can be explained by the four extracted factors, with EC (0.943) and Mg²⁺ (0.927) having the greatest communalities higher than 0.9.