Spatio-temporal Patterns and Source Identification of Water Pollution in Lake Taihu (china)

Various multivariate methods were used to analyze datasets of river water quality for 11 variables measured at 20 different sites surrounding Lake Taihu from 2006 to 2010 (13,200 observations), to determine temporal and spatial variations in river water quality and to identify potential pollution sources. Hierarchical cluster analysis (CA) grouped the 12 months into two periods (May to November, December to the next April) and the 20 sampling sites into two groups (A and B) based on similarities in river water quality characteristics. Discriminant analysis (DA) was important in data reduction because it used only three variables (water temperature, dissolved oxygen (DO) and five-day biochemical oxygen demand (BOD 5)) to correctly assign about 94% of the cases and five variables (petroleum, volatile phenol, dissolved oxygen, ammonium nitrogen and total phosphorus) to correctly assign >88.6% of the cases. In addition, principal component analysis (PCA) identified four potential pollution sources for Clusters A and B: industrial source (chemical-related, petroleum-related or N-related), domestic source, combination of point and non-point sources and natural source. The Cluster A area received more industrial and domestic pollution-related agricultural runoff, whereas Cluster B was mainly influenced by the combination of point and non-point sources. The results imply that comprehensive analysis by using multiple methods could be more effective for facilitating effective management for the Lake Taihu Watershed in the future.


Introduction
The water quality in lakes has recently become a matter of great concern, because of its negative effects on human health and its contribution to harmful algal blooms.Accompanying the rapid growth of human population and industries, lake water is being heavily polluted from various sources, like domestic sewage, industrial wastewater, stormwater, agricultural runoff and atmospheric depositions [1][2][3][4].Rivers, the natural pathways upstream of lakes, deliver most of the pollutant loadings into lakes [5,6].To prevent and control lake water pollution, it is important to understand the temporal and spatial variations in upstream river water quality through routine monitoring programs.Currently, most monitoring programs generate datasets with information including physical properties, Water 2016, 8, 86 2 of 11 aggregate organic constituents and nutrients and inorganic constituents, making it difficult to analyze and interpret the underlying spatio-temporal patterns of water pollution [7,8].Similarly, due to the complex interactions between natural processes and human activities, identifying potential pollution sources for river basins has also become a serious issue [7], which would eventually compromise the resilience and restoration of lakes.
Over the last decade, multivariate approaches have significantly advanced our understanding of the spatio-temporal patterns and pollution sources of river systems on the basis of water quality observations [8][9][10][11].For instance, cluster analysis (CA), discriminant analysis (DA) and principal component analysis (PCA) have been widely used to assess spatial or temporal variations of groundwater and surface water quality [12][13][14][15][16].However, reciprocal validation between any two multivariate analyses is still required to validate the results obtained by each approach [7].Although the spatial and temporal variations of surface water quality in reservoirs or lakes have been extensively studied [17][18][19], how to identify the pollution sources for watershed systems still remains elusive due to the lack of pollution census data [20][21][22].
In this study, on the basis of a dataset obtained from a five-year (2006-2010) water-monitoring program in Lake Taihu Watershed in China, both CA and DA were performed to identify the spatio-temporal patterns in river water quality.PCA analysis was further performed to identify the underlying pollution sources in different regions.In the end, this study discussed the implications of local watershed management.

Study Area and Monitoring Program
Lake Taihu, the third largest freshwater lake in China, is situated in southeast Jiangsu Province and the lowest reach of the Yangtze River Basin, China (30 ˝55 1 40 11 -31 ˝32 1 58 11 N; 119 ˝52 1 32 11 -120 ˝36 1 10 11 E; Figure 1).The lake has a surface area of 2338 km 2 with a maximum depth of 2.6 m and 1.9 m on average.The average water residence time is 309 days [23].Lake Taihu is one of the main drinking water sources for neighboring residents, albeit one of the most heavily-polluted freshwater lakes in China [24].Rapid economic expansion, population growth and industrial development surrounding Lake Taihu are accompanied with severe industrial pollution with a constant increase in domestic sewage and intensification of agricultural production [25,26], resulting in significant enrichment of organic substances and nutrients discharged into the surrounding riverine systems.This led to an overall decline in rivers' water quality and ecosystem services in the Lake Taihu Watershed [27].
Water quality data was observed monthly between January 2006 and December 2010 at 20 monitoring sites in the rivers surrounding Lake Taihu (Figure 1), where the upstream catchment area and annual mean streamflow were illustrated in the Supplementary Material, Table S1.The monitored rivers contribute approximately 70% of total discharge into Lake Taihu [28].Eleven water quality variables were selected for further analysis, i.e., water temperature (Temp), pH, conductivity (Cond), dissolved oxygen (DO), chemical oxygen demand (COD Mn ), five-day biochemical oxygen demand (BOD 5 ), ammonium nitrogen (NH 4 + -N), total phosphorus (TP), total petroleum hydrocarbons (Petro), volatile phenol (V-ArOH) and plumbum (Pb).Other variables (e.g., cadmium, cuprum, cyanide) were not considered in this study, because their concentrations are less than the detection limits.Selected variables refer to physical properties, organic constituents, nutrient constituents or biological properties of the rivers (

Multivariate Statistical Analysis
Spatio-temporal patterns in river water quality surrounding Lake Taihu were determined by CA using Ward's method with squared Euclidean distance [7,29,30].The linkage distance is expressed as (Dlink/Dmax) × 100, which represents the standardized quotient between the linkage distances for a particular case divided by maximal linkage distance [30].DA was used in backward stepwise mode to confirm the groups found by CA and to evaluate the spatio-temporal variations of the discriminant variables.In DA, the monitoring period or site variables were the clustering variables, while the parameters from originally-measured datasets were independent [11,31].Principal component analysis (PCA) is based on the assumption that there exists a bilinear model, which could explain the variance of observed water quality data by using less orthogonal variables, known as principal components [32].Once Kaiser's VARIMAX rotation is performed, factor Table 1.Univariate statistics and analytical methods in the analysis of the water quality variables for rivers surrounding Lake Taihu.SD and CV stand for standard deviation and coefficient of variation, respectively.Temp, temperature; Cond, conductivity; DO, dissolved oxygen; COD Mn , chemical oxygen demand; BOD 5 , chemical oxygen demand; TP, total phosphorus; Petro, total petroleum hydrocarbons; V-ArOH, volatile phenol.

Multivariate Statistical Analysis
Spatio-temporal patterns in river water quality surrounding Lake Taihu were determined by CA using Ward's method with squared Euclidean distance [7,29,30].The linkage distance is expressed as (D link /D max ) ˆ100, which represents the standardized quotient between the linkage distances for a particular case divided by maximal linkage distance [30].DA was used in backward stepwise mode to confirm the groups found by CA and to evaluate the spatio-temporal variations of the discriminant variables.In DA, the monitoring period or site variables were the clustering variables, while the parameters from originally-measured datasets were independent [11,31].Principal component analysis (PCA) is based on the assumption that there exists a bilinear model, which could explain the variance of observed water quality data by using less orthogonal variables, known as principal components [32].Once Kaiser's VARIMAX rotation is performed, factor loadings remain orthogonal and are no longer contributed to the explained maximum variance, while the scores become non-orthogonal.Using PCA, we could identify the unobservable latent pollution sources that affect the water quality of the upstream rivers.The details of multivariate analysis can be found in [30].

Data Pretreatment
The following data pretreatment methods were used for the river water quality dataset: (i) missing data were linearly interpolated based on the corresponding datasets; (ii) values below lower bounds were replaced with the mid-value between zero and the detection limit; (iii) the normality test was performed for each variable by analyzing its kurtosis and skewness, since most multivariate statistical analyses require the variables to be normally distributed.Unfortunately, skewness and kurtosis statistical test results suggest that most of the variables from the original data were severely deviated from the normal distribution at a confidence level of 95%.To satisfy the normality assumption, original data were logarithmically transformed, and further skewness and kurtosis values for the log-transformed data were significantly reduced and close to those for the normal distribution.All log-transformed variables were further z-scale standardized for CA and PCA analysis, while non-transformed original data were still used for DA analysis.

Temporal Similarity and Variation
Temporal hierarchical CA analysis generated a transient dendrogram and grouped the 12 months into either a two-cluster system at (D link /D max ) ˆ100 < 20 or a three-cluster system at (D link /D max ) ˆ100 < 16 (selecting the scale tree to (D link /D max ) ˆ100 option in order for the tree plot to be scaled to a standardized scale).The temporal difference in water quality between the two systems was significant (Figure 2a; p < 0.01).In the two-cluster system, Cluster 1 (first period) covered the time from May to November while Cluster 2 (second period) ranged from December to the next April.In three-cluster system, while Cluster 1 remains the same as the previous system, Cluster 2 was further divided into two clusters: Cluster 2 for April, December and March and Cluster 3 for January to February (Figure 2a).Such a grouping differed from the Lake Dianchi Watershed, which is located in Southwest China [33], because the length of the wet season in the Lake Taihu Watershed (May to October) is much longer than the Lake Dianchi Watershed (July to September).In both twoand three-cluster systems, the temporal variations from river water quality data from Lake Taihu were determined by both hydrological condition (i.e., wet or dry season) and water pollution characteristics (Figure 2b).For example, NH 4 + -N in the first period (May to November) was much higher than the other clusters (Figure 2c), due to the greater contribution of agricultural runoff in the wet season.However, the mean concentrations of BOD 5 and COD Mn in the first period were close to the other two periods (Figure 2c), because the relatively larger amount of domestic sewage was offset by the greater streamflow in the rivers.DA was then applied to evaluate the clusters' systems generated by the temporal CA method.The objectives of the DA were to test the significance of discriminant functions and to choose the most significant variables that contributed to the differences among clusters.For each discriminant function, the results of Wilks' lambda and chi-square analyses varied from 0.11 to 0.63 and 668 to 1972, respectively, with p < 0.001.This suggests that the temporal DA was reliable and effective [29,30].In the two-cluster scenario, DA produced two classification matrices (CMs) with 94% accuracy of classification using three discriminant variables: water temperature, DO and BOD 5 (Figure 2c).In the three-cluster scenario, DA produced three classification matrices (CMs) with 86.5% accuracy using two discriminant variables (i.e., COD Mn and NH 4 + -N), which were significantly different from each other among the three-cluster systems (Figure 2c).For instance, the average means of the Water 2016, 8, 86 5 of 11 three variables (BOD 5 , COD Mn and NH 4 + -N) in the third period were 3.4%, 8.4% and 12.5% higher than those of the second period (Figure 2c), respectively.The discrepancy between first and third period was similar to that between first and second period, whereas the gaps were much bigger for NH 4 + -N (Figure 2c).The two-cluster system divided the 12 months into wet and dry seasons, while the three-cluster system divided the 12 months into wet (May to November), moderate dry (March, April and December), severe dry seasons (January to February), according to local meteorological characteristics.Together, the backward stepwise DA results suggested that both two-and three-cluster systems explained the temporal similarities well.Water temperature, DO and BOD 5 were the three most significant variables in discriminating the water quality condition in different seasons in both systems.Although the coefficient of variation (CV; Table 1) could reflect the numerical variations of samples, it is inappropriate to be used to evaluate the temporal variation among different seasons in either system.
Water 2016, 8, x 5 of 10 moderate dry (March, April and December), severe dry seasons (January to February), according to local meteorological characteristics.Together, the backward stepwise DA results suggested that both two-and three-cluster systems explained the temporal similarities well.Water temperature, DO and BOD5 were the three most significant variables in discriminating the water quality condition in different seasons in both systems.Although the coefficient of variation (CV; Table 1) could reflect the numerical variations of samples, it is inappropriate to be used to evaluate the temporal variation among different seasons in either system.

Spatial Similarity and Variation
Through spatial CA analysis, we identified clusters of similar monitoring sites considering the effects of the temporal differences in spatial CA.Spatial similarity analysis was conducted for each individual temporal cluster, as well as all of the samples combined.Our results indicate that there were no significant differences among them, and further discussion will be focused on the spatial CA for all samples' combined data.Spatial CA produced two dendrograms with two clusters at (Dlink/Dmax) × 100 < 20 and three clusters at (Dlink/Dmax) × 100 < 8, respectively.In the two-cluster dendrogram, Cluster A covered the S2, S8, S13, S19 and S20 sites, while Cluster B covered the S1, S3 to S7, S9 to S12 and S14 to S18 sites (Figure 2b).In the three-cluster dendrogram, Cluster B from the previous two-cluster dendrogram was further split into two clusters: new Cluster B with the S4, S7,

Spatial Similarity and Variation
Through spatial CA analysis, we identified clusters of similar monitoring sites considering the effects of the temporal differences in spatial CA.Spatial similarity analysis was conducted for each individual temporal cluster, as well as all of the samples combined.Our results indicate that there were no significant differences among them, and further discussion will be focused on the spatial CA for all samples' combined data.Spatial CA produced two dendrograms with two clusters at Water 2016, 8, 86 6 of 11 (D link /D max ) ˆ100 < 20 and three clusters at (D link /D max ) ˆ100 < 8, respectively.In the two-cluster dendrogram, Cluster A covered the S2, S8, S13, S19 and S20 sites, while Cluster B covered the S1, S3 to S7, S9 to S12 and S14 to S18 sites (Figure 2b).In the three-cluster dendrogram, Cluster B from the previous two-cluster dendrogram was further split into two clusters: new Cluster B with the S4, S7, S9, S11, S15 and S18 sites and new Cluster C with the rest of them (Figure 2b).All classifications varied at a significance level of p < 0.01, which meets our expectation, since the sites in the same cluster shared similar natural backgrounds and had been affected by similar sources in a similar way.
The spatial DA was performed similarly as the temporal DA (Table 2).We performed Wilks' lambda and the chi-square analysis on each discriminant function.The values were within a range of 0.50 to 0.61 and 523 to 689, respectively, suggesting that the spatial DA had a similar discriminatory ability as the temporal DA.The spatial DA was performed using the original dataset with 11 variables after classifying into the two major groups (A and B) obtained from the spatial CA.Sites were the dependent variables, and the measured parameters were the independent variables.Backward mode discriminant functions successfully assigned >88.6% and >81.9% of the cases into the two-and three-cluster systems, respectively (Table 2).Moreover, the backward stepwise DA demonstrated that Petro, V-ArOH, DO, NH 4 + -N and TP were also significant discriminant variables for spatial variation (Figure 2d).The sites in Cluster A were situated in the highly developed area (i.e., Wujin District or Xishan District of Changzhou City) (Figure 3), where most of the industrial effluents and domestic sewage flow into the rivers directly.Most of the sites in Cluster B were located in the Tiaoxi River Basin (Figure 3).Tiaoxi River is the largest tributary of Lake Taihu, and it originates from mountainous area and moderately developed rural regions.The sites in Cluster C were located in northwestern and eastern Lake Taihu (Figure 3) in Yixing City and Wuzhong District of Suzhou City, where the major pollution sources include both point and non-point sources.

Identification of Potential Pollution Sources
Due to the similarity between two-and three-cluster systems (see Supplementary Material, Table S2), source identification of water pollution for the two-cluster system was only illustrated.Before conducting the PCA analysis, the Kaiser-Meyer-Olkin (KMO) and Bartlett's sphericity tests were performed on the parameter correlation matrix.The KMO results for Clusters A and B were 0.56 and 0.52, respectively, and Bartlett's sphericity results were 861 and 812 (p < 0.05), indicating that PCA could be used in dimensionality reduction.PCAs were applied to standardized log-transformed datasets (11 variables) to examine the differences between Clusters A and B and to identify the latent factors.PCA with VARIMAX rotation explained 75.6% and 67.0% of the total variance in Clusters A and B, respectively (Table 3).Such a performance of source identification was close to those for Lake Dianchi and Lake Chaohu in China [33,34].
flow into the rivers directly.Most of the sites in Cluster B were located in the Tiaoxi River Basin (Figure 3).Tiaoxi River is the largest tributary of Lake Taihu, and it originates from mountainous area and moderately developed rural regions.The sites in Cluster C were located in northwestern and eastern Lake Taihu (Figure 3) in Yixing City and Wuzhong District of Suzhou City, where the major pollution sources include both point and non-point sources.For Cluster A, the first varifactor (VF1), which explained 18.6% of the total variance, had only strong positive loadings on Petro and Pb, but a moderate loading on V-ArOH (Table 3).The element Pb is mainly from electronic manufacturing and chemical industries; V-ArOH is from paper-making and chemical industries; and Petro is from equipment manufacturing, metal smelting industries and chemical industries.VF1 represented chemical pollution, which is originated from industrial wastewater and discharged into the rivers.VF2 represented domestic pollution, explaining 17.09% of the total variance, and had strong positive loadings on TP, COD Mn and BOD 5 .VF3 could be interpreted as N-related industrial pollution, which accounted for 16.7% of the total variance and had strong positive loadings on conductivity (Cond) and NH 4 + -N.VF4 and VF5 explained 13.29% and 9.97% of the total variance, respectively.VF4 had strong positive loadings on water temperature, but strong negative loadings on DO, while VF5 had only strong positive loadings on pH.VF5 was attributed to the variability from the physicochemical source and represented natural sources impacted by seasonality.
According to the Pollution Source Census Survey of Jiangsu and Zhejiang provinces, V-ArOH, identified in VF1, is majorly generated by chemical manufactures in Lake Taihu area, and approximately 35% of V-ArOH is from the Wujin District of Changzhou City (Figure 3), a Cluster A catchment area.Major chemical manufactures in Wujin District were under tight regulation by local government, but the total V-ArOH discharge is still massive.Moreover, 89% of Pb discharge was from communication and electronic manufactures that are widely spread in the Class A catchment areas in Wujin District of Changzhou City and Xishan District of Wuxi City.Similarly, petroleum-related emissions are also from electronic manufacturing and chemical industries.
Since urbanization expansion was accompanied with population growth, domestic pollution became the primary source and is represented by VF2 (Table 3), with the major factors of TP, COD Mn and BOD 5 .High population density could lead to massive organic pollution without proper regulation.In the Class A catchment area, the population density reached 1148 persons/km 2 , and industrial wastewater discharge density was close to 90,000 tons/km 2 annually (Figure 3) Moreover, agricultural runoff of TN was 2.5 tons/km 2 annually.Therefore, domestic pollution in this area is becoming a serious issue and should be carefully considered within the pollution control plan in the future.
Further analysis of VF3 from historic statistical data suggests that nitrogen (N) emission was primarily from industrial wastewater.In the Cluster A catchment area, 55% of NH 4 + -N was from industrial wastewater, which is moderately correlated with V-ArOH, 31% from agricultural runoff and 14% from domestic sewage.Since the combined N emissions affect electronic conductivity, the potential pollution source VF3 could also be explained as N-related industrial pollution.For Cluster B, VF1 (accounting for 27.14% of the total variance) had strong positive loadings on COD Mn , TP, NH 3 + -N, conductivity and BOD 5 , which represent the combination of point and non-point sources.For instance, domestic wastewater discharge per area in Cluster B was only 34,000 tons/km 2 annually, which is less than half that in Cluster A (Figure 3).However, agricultural runoff of TN was up to 2.9 tons/km 2 annually, close to the intensity of that in Cluster A (Figure 3).As previously mentioned, VF2 explained 14.3% of the total variance and had positive loadings on water temperature, but strongly negative on DO.It represented natural sources impacted by seasonal change and hydrological conditions (Singh et al. [29]; Zhou et al., [26]).VF3 (13.77% of the total variance) was positively weighted by petroleum-related pollutions and had negative loadings on pH.Previously, we demonstrated that petroleum-related emissions were from multiple industries and represented the intensity of industrial development in a certain area.VF4 explained 11.7% of the total variance and had strong positive loadings on Pb and V-ArOH.Similar to the Cluster A area, VF4 was categorized as a chemical-related industrial pollution factor.Pb and V-ArOH discharges were mainly from electronic and chemical manufactures.However, it was still considered as an independent impact factor due to continuing economic growth and industrial development in this area.Analysis of the major factors and main pollution patterns in the highly-polluted area (Class A area) and the moderately-polluted area (Class B area) revealed that there were significant differences between these two groups.The Cluster A area was severely impacted by the heavy chemical industries.Recently, the Cluster A area performed better in controlling pollution under strong regulations and pollution control.Pb and V-ArOH levels were basically below the detection limits.Nevertheless, due to the huge total discharge amount, stronger regulations and pollution controls would be still required to reduce emission amounts from chemical manufacturing, electronic and communication manufactures, compared to current conditions.The most prominent source in the Cluster B area was domestic pollution, whereas industrial pollution, as an independent major factor, could not be ignored either.Moreover, since natural conditions are the second important pollution factor in this area, different seasons with distinct precipitation would result in different lake water quality.

Conclusions and Implications
Multivariate statistical methods were successfully applied to evaluate temporal and spatial variations in studying river water quality and identifying pollutions sources at river outlets surrounding Lake Taihu Watershed.Our results suggest that multiple methods are effective and compatible with each other and could be used in river water quality management in the future.Hierarchical CA clustered the 12 months into three periods and classified 16 sampling sites into two groups (A and B) based on the similarity of water quality characteristics.Both temporal and spatial DA analysis had the best performance with good discriminatory ability according to significance validation tests.They also identified several significant variables for discrimination among temporal or spatial groups.Analysis of temporal variation by DA required only three variables, but successfully assigned about 94% of the cases, and analysis of spatial variation required only five variables, but with more than 88.6% cases successfully assigned.In conclusion, the temporal and spatial similarities and differences could optimize monitoring programs with decreased monitoring frequency and a decreased number of sampling monitoring stations and monitoring variables, which could finally significantly reduce the subsequent costs.Moreover, PCA analysis identified four and five latent pollution sources for Clusters A and B, respectively, which are industrial source (chemical-related, petroleum-related or N-related), domestic source, a combination of point and non-point sources and natural source.Overall, our study provides important information in the understanding and characterization of pollution patterns surrounding the Lake Taihu area.However, how to accurately quantify the contributions of different pollution resources still remains elusive.Future study needs supplemental investigations on potential pollution resources and continuing monitoring of upstream rivers.Furthermore, comprehensive analysis in combination by using multiple methods could be more effective for facilitating effective management for the Lake Taihu Watershed in the future.

Figure 1 .
Figure 1.Study area and its monitoring sites.

Figure 1 .
Figure 1.Study area and its monitoring sites.

Figure 2 .
Figure 2. Temporal and spatial similarities of monitoring periods or sites produced by cluster analysis (CA).The values for water temperature, TP, Petro and V-ArOH are re-scaled by multiplying by 0.1, 10, 100 and 100.The units for concentration or temperature are mg/L and °C.(a) Dendrogram for temporal similarity of monitoring periods; (b) Dendrogram for spatial similarity of monitoring sites;(c) Differences in discriminant variables among two-cluster system; (d) Same as panel c but for three-cluster system.

Figure 2 .
Figure 2. Temporal and spatial similarities of monitoring periods or sites produced by cluster analysis (CA).The values for water temperature, TP, Petro and V-ArOH are re-scaled by multiplying by 0.1, 10, 100 and 100.The units for concentration or temperature are mg/L and ˝C.(a) Dendrogram for temporal similarity of monitoring periods; (b) Dendrogram for spatial similarity of monitoring sites; (c) Differences in discriminant variables among two-cluster system; (d) Same as panel c but for three-cluster system.

Figure 3 .
Figure 3. Spatial pattern of the grouping of monitoring sites and the corresponding emission rates per area of industrial source, domestic source and agricultural runoff.The grouping is determined by the three-cluster systems.The wastewater amounts of industrial and domestic pollution sources are extracted from the Pollution Source Census Survey of Jiangsu and Zhejiang provinces for the period pf 2007 to 2010, while the agricultural runoff of TN is obtained from Zhou et al. [31] and Hou et al. [32].The emission density of individual pollution sources in each cluster is shown in the inset of the top-left of the map.

Figure 3 .
Figure 3. Spatial pattern of the grouping of monitoring sites and the corresponding emission rates per area of industrial source, domestic source and agricultural runoff.The grouping is determined by the three-cluster systems.The wastewater amounts of industrial and domestic pollution sources are extracted from the Pollution Source Census Survey of Jiangsu and Zhejiang provinces for the period pf 2007 to 2010, while the agricultural runoff of TN is obtained from Zhou et al. [31] and Hou et al. [32].The emission density of individual pollution sources in each cluster is shown in the inset of the top-left of the map.

Table 1 .
Univariate statistics and analytical methods in the analysis of the water quality variables for rivers surrounding Lake Taihu.SD and CV stand for standard deviation and coefficient of variation, respectively.Temp, temperature; Cond, conductivity; DO, dissolved oxygen; CODMn, chemical oxygen demand; BOD5, chemical oxygen demand; TP, total phosphorus; Petro, total petroleum hydrocarbons; V-ArOH, volatile phenol.

Table 2 .
Classification matrices for backward discriminant analysis (DA) of temporal and spatial variations.

Table 3 .
Loadings of 11 measured variables on VARIMAX rotated factors of two clusters.