Groundwater Quality Evaluation of the Dawu Water Source Area Based on Water Quality Index (WQI): Comparison between Delphi Method and Multivariate Statistical Analysis Method

Water quality in the Dawu water source area is primarily studied through the traditional water quality assessment method that measures the chemical parameters of water against the existing standards, which is simple but cannot accurately describe the water quality. Although the water quality index (WQI) proposed by Horton is widely used for comprehensive water quality evaluation, parameter selection and weight determination are primarily based on the Delphi method, which is subjective and random. Moreover, in groundwater evaluation, the focus is primarily laid on general chemical parameters, such as Total Dissolved Solids, hydrogen ion concentration, Electrical Conductivity, and heavy metals, such as Hydrargyrum, Arsenic, and Chromium, with limited consideration for organic pollutants. In this study, WQI technology in combination with the entropy weight method was used to evaluate the groundwater environmental quality of the Dawu water source area, and the scientific results were analyzed by comparing the full index, Delphi, and multivariate statistical analysis methods. The results showed that the groundwater in the Dawu water source area generally had good quality and was potable and that the application of multivariate statistical analysis method was more suitable than the Delphi method in the index selection process.


Introduction
Groundwater is an important component of water resources and acts as a large natural "reservoir" that supplies water even in conditions of low surface water availability (i.e., dry season or interannual periods of prolonged drought). Using the WaterGAP Global Hydrology Model, global groundwater recharge was estimated at 12,666 km 3 /a during 1961-1990 under normal climate conditions, accounting for 32% of the total renewable water resources [1]. In 2017, China's total groundwater resources accounted for 28.89% of its total water resources, and groundwater supply accounted for 16.82% of its total water supply. However, groundwater pollution is a serious problem in China [2]. In 2019, among the 10,168 state-level groundwater quality monitoring sites in China, 14.4% were classified as category 1 to 3, 66.9% as category 4, and 18.8% as category 5. The above standard indices were based on Mn, Total hardness (TH), I, Total Dissolved Solids (TDS), velop WQI technology; the groundwater environment quality of Dawu water source area was evaluated, and the results were compared; the scientific basis of the multivariate statistical analysis method in the process of index selection was verified, and the groundwater pollution situation of the Dawu water source area was analyzed to provide a theoretical basis for sustainable development of the area.

Materials and Methods
The Dawu water source area is located at 36 • 30 31 -36 • 51 12 N, 118 • 02 16 -118 • 26 55 E, Linzi District, Zibo City, Shandong Province, China. It includes a part of Fengshui Town in Zhangdian District in the west, parts of the northeastern towns in Zichuan District in the south, and the western part of Qingzhou District in Weifang City in the east. The terrain slope gradually decreases from the south of the low hills to the northern piedmont slope plain. Figure 1 shows the location of the Dawu water source area.
Water 2021, 13, x FOR PEER REVIEW 3 of 17 In this study, the full index method, Delphi method, principal component analysis (PCA) approach combined with entropy weight assignment method were adopted to develop WQI technology; the groundwater environment quality of Dawu water source area was evaluated, and the results were compared; the scientific basis of the multivariate statistical analysis method in the process of index selection was verified, and the groundwater pollution situation of the Dawu water source area was analyzed to provide a theoretical basis for sustainable development of the area.

Materials and Methods
The Dawu water source area is located at 36°30′31′′-36°51′12′′ N, 118°02′16′′-118°26′55′′ E, Linzi District, Zibo City, Shandong Province, China. It includes a part of Fengshui Town in Zhangdian District in the west, parts of the northeastern towns in Zichuan District in the south, and the western part of Qingzhou District in Weifang City in the east. The terrain slope gradually decreases from the south of the low hills to the northern piedmont slope plain. Figure 1 shows the location of the Dawu water source area. This study was based on the 3D visualization information system construction of the water source of Zibo City from October 2017 to January 2018; 101 samples were collected with a sampling depth of 150-520 m. There were 122 detection indices, including the general chemical indices, such as metal ions, pH, and total hardness, and toxicology indices, This study was based on the 3D visualization information system construction of the water source of Zibo City from October 2017 to January 2018; 101 samples were collected with a sampling depth of 150-520 m. There were 122 detection indices, including the general chemical indices, such as metal ions, pH, and total hardness, and toxicology indices, such as nitrite, nitrate, heavy metals, organic pollutants, and pesticides. The sampling points were mining wells or pressure wells that often supply water, and the water was pumped for more than 10 mins before collection to ensure that the water samples obtained can reflect the real conditions of groundwater in the sampling points. Samples were collected in brown sampling bottles, which were pickled and washed with distilled water before sampling and washed with fresh water at least three times during sampling, and pH was measured on-site. After sampling, the container was immediately stoppered to avoid air gap, and the sample was promptly sent for analysis. Organic matter was detected in consonance with Standard examination methods for drinking water (GB/T5750) [22], primarily using Agilent6890N gas chromatograph (United States) and TU1810 UV-visible spectrophotometer (China). Groundwater samples were primarily karst water, which also accounted for pore water. In a preliminary analysis of the index data, none of the 37 pesticides were detected; eventually, the 33 organic compounds, such as dibromomethane, bromodichloromethane, and 1,2-dibromoethane, were also not detected. At the same time, referring to Standards for groundwater quality (GT14848-2017) [23], Standards for drinking water quality ( GB5749-2006) [24], The European Union Water Framework Directive (hereafter referred to as the "WFD") in 2000 [25] and the World Health Organization (WHO) Guidelines for drinking-water quality (the fourth edition) in July 2011 [26], the similar index and other indices not included in the scope of standards and guidelines were removed. Finally, 38 indices were determined as within the scope of this study. If the detected concentration was below the detection limit, it was set as the standard detection concentration for evaluation. Figure 2 shows the distribution of the sampling Wells. such as nitrite, nitrate, heavy metals, organic pollutants, and pesticides. The sampling points were mining wells or pressure wells that often supply water, and the water was pumped for more than 10 mins before collection to ensure that the water samples obtained can reflect the real conditions of groundwater in the sampling points. Samples were collected in brown sampling bottles, which were pickled and washed with distilled water before sampling and washed with fresh water at least three times during sampling, and pH was measured on-site. After sampling, the container was immediately stoppered to avoid air gap, and the sample was promptly sent for analysis. Organic matter was detected in consonance with Standard examination methods for drinking water (GB/T5750) [22], primarily using Agilent6890N gas chromatograph (United States) and TU1810 UVvisible spectrophotometer (China). Groundwater samples were primarily karst water, which also accounted for pore water. In a preliminary analysis of the index data, none of the 37 pesticides were detected; eventually, the 33 organic compounds, such as dibromomethane, bromodichloromethane, and 1,2-dibromoethane, were also not detected. At the same time, referring to Standards for groundwater quality (GT14848-2017) [23], Standards for drinking water quality ( GB5749-2006) [24], The European Union Water Framework Directive (hereafter referred to as the "WFD") in 2000 [25] and the World Health Organization (WHO) Guidelines for drinking-water quality (the fourth edition) in July 2011 [26], the similar index and other indices not included in the scope of standards and guidelines were removed. Finally, 38 indices were determined as within the scope of this study. If the detected concentration was below the detection limit, it was set as the standard detection concentration for evaluation. Figure 2 shows the distribution of the sampling Wells. In this study, the full index method, Delphi method, and PCA were adopted to evaluate the groundwater quality of the Dawu water source area. The research process is shown in the Figure 3 below. In this study, the full index method, Delphi method, and PCA were adopted to evaluate the groundwater quality of the Dawu water source area. The research process is shown in the Figure 3 below.
The Delphi method, also known as the expert survey method, was initiated and implemented by Dalkey and Helmer for RAND Corporation in 1946 [27]. It is primarily a feedback and anonymous consultation method; the general process is to collect opinions of experts on the problems to be predicted; then organize, statistically summarize, and give anonymous feedback to experts; then solicit opinions again, analyze and give feedback again, until a consistent opinion is obtained [28,29]. The Delphi method is a qualitative research method that can anonymously accept a large number of individual opinions from different sources and professional fields. Participants are allowed to openly provide criticism, which can thereby minimize social pressure or the influence of experts [27,30]. The standard Delphi method has five primary characteristics: the use of experts, anonymity, controlled feedback, multi-stage iteration, and the search for consensus through the statistical aggregation of group responses [31]. The Delphi method, also known as the expert survey method, was initiated and implemented by Dalkey and Helmer for RAND Corporation in 1946 [27]. It is primarily a feedback and anonymous consultation method; the general process is to collect opinions of experts on the problems to be predicted; then organize, statistically summarize, and give anonymous feedback to experts; then solicit opinions again, analyze and give feedback again, until a consistent opinion is obtained [28,29]. The Delphi method is a qualitative research method that can anonymously accept a large number of individual opinions from different sources and professional fields. Participants are allowed to openly provide criticism, which can thereby minimize social pressure or the influence of experts [27,30]. The standard Delphi method has five primary characteristics: the use of experts, anonymity, controlled feedback, multi-stage iteration, and the search for consensus through the statistical aggregation of group responses [31].
Multivariate statistical analysis was developed from classical statistics. It is a comprehensive analysis method that can analyze the statistical law of multiple objects and multiple indices under the condition that they are interrelated. Multivariate statistical techniques have been widely used in water quality data analysis [32][33][34][35][36][37]. These tools help simplify and organize large datasets to explain the observed relationships between several variables [38]. In this study, correlation analysis, PCA, and principal factor analysis (PFA) were used to select evaluation index parameters, and IBM SPSS Statistics 26 software (United States) was used for data visualization. Each water quality parameter contained important water quality information, and there may be different degrees of correlation between the parameters. Through comprehensive analysis of these correlations, we can obtain the water quality information from the data [35,36,[38][39][40]. Based on PCA, the water quality indices obtained from the same principal component had similar or the same sources contributing to pollution, while the water quality parameters within different principal components represent different sources of pollutants. Based on the correlation, PFA was carried out and factor load calculated. The contribution rate of water quality parameters in principal components can be determined through factor load calculation, to identify the primary pollution impact factors [38,[40][41][42].
The entropy weight method is based on Shannon entropy, proposed by Shannon in 1948 [43]. According to Shannon's theory, information is a function: the sender of the information chooses a piece of information from a set of possible choices, and the degree of Multivariate statistical analysis was developed from classical statistics. It is a comprehensive analysis method that can analyze the statistical law of multiple objects and multiple indices under the condition that they are interrelated. Multivariate statistical techniques have been widely used in water quality data analysis [32][33][34][35][36][37]. These tools help simplify and organize large datasets to explain the observed relationships between several variables [38]. In this study, correlation analysis, PCA, and principal factor analysis (PFA) were used to select evaluation index parameters, and IBM SPSS Statistics 26 software (United States) was used for data visualization. Each water quality parameter contained important water quality information, and there may be different degrees of correlation between the parameters. Through comprehensive analysis of these correlations, we can obtain the water quality information from the data [35,36,[38][39][40]. Based on PCA, the water quality indices obtained from the same principal component had similar or the same sources contributing to pollution, while the water quality parameters within different principal components represent different sources of pollutants. Based on the correlation, PFA was carried out and factor load calculated. The contribution rate of water quality parameters in principal components can be determined through factor load calculation, to identify the primary pollution impact factors [38,[40][41][42].
The entropy weight method is based on Shannon entropy, proposed by Shannon in 1948 [43]. According to Shannon's theory, information is a function: the sender of the information chooses a piece of information from a set of possible choices, and the degree of uncertainty associated with this choice behavior is central to this concept [44]. Consistent with the basic principles of information theory, information is a measure of the degree of system order, and entropy is a measure of the degree of system disorder. In line with the definition of information entropy, for a certain index, the entropy value can be used to determine the degree of dispersion of an index. A smaller information entropy value will lead to a greater degree of dispersion of an index, and a greater influence (i.e., weight) of the index on the comprehensive evaluation. If all the values of an index are equal, the index will not play any role in the comprehensive evaluation [43,[45][46][47]. Therefore, the information entropy tool can be used to calculate the weight of each index and provide a basis for the comprehensive evaluation of multiple indices. Li et al. [48] used the entropy weight fuzzy analysis method to evaluate the groundwater quality of industrial parks in Northwest China.
Indices can be divided into positive and negative indices. Positive indices are also known as efficiency criteria; a high original value of the index will result in a higher efficiency level, whereas a low original value of the index will result in a lower level. The inverse index is also known as the cost index. A smaller original value of the index will result in a higher level; a higher original value of the index will result in a lower level [45,49,50]. In the water quality index, except pH, all others are inverse indices. If the membership function is U w (x), the membership degree of the positive index is the inverse index membership degree is where x max and x min represent the maximum and minimum of the original data values of m samples of the index, respectively. In Equations (1) and (2), U w (x) ∈ [0,1]. The closer the membership function value is to 1, the better the water quality. On the basis of the membership degree of each index, the membership matrix, namely the decision matrix R w , was obtained as where r ij = U w (x); i = 1,2, . . . , m, j = 1,2, . . . , n; m is the number of samples; and n is the number of indices. The entropy method was used to determine the weight of each index, and the weight transformation was carried out to calculate the new weight of each index. The specific gravity of the index was transformed as E j represents the total contribution of all samples to the index, namely the entropy value, Here, the constant K = 1/ln (m), such that 0 ≤ E j ≤ 1 can be ensured; that is, E j has a maximum value of 1.
The reverse entropy was obtained as Finally, the index weight was calculated as The WQI is the simplest and most widely used index for assessing the overall quality of surface and groundwater water [17,33,51,52]. In this study, the WQI was [48] calculated as where w j is the weight, which ranges between 0 and 1, w 1 + w 2 + w 3 + . . . + w n = 1; q j is the quality level of each parameter, and n is the number of parameters used in WQI evaluation.
where C j is the concentration (mg/L) of each index in each water sample, C pH is the value of the pH in each water sample, and S j is the standard limit (mg/L) of each index in the Standards for drinking water quality (GB5749-2006).

Single Factor Evaluation
The research scope includes 38 indices, such as Na + , NH 4 + -N, Fe, TH, TDS, Hg, Lead (Pb), benzene, and xylene (total). First, all index data of 101 samples were summarized and preliminarily analyzed and compared with Standards for drinking water quality (GB5749-2006) (hereafter referred to as the "Standards"), the European Union Water Framework Directive (hereafter referred to as the "WFD") of 2000, and the WHO Guidelines for drinking-water quality (the fourth edition) of July 2011 (hereafter referred to as "Guidelines"). The results are shown in Table 1.
Among them, 16 indexes including NH 4 + -N, Fe, Cl − , SO 4 2− , pH, TH, TDS, chemical oxygen demand (permanganate index) (hereinafter referred to as "COD(Mn)"), Aluminum (Al), Manganese (Mn), Molybdenum (Mo), Argentum (Ag), 1,1-dichloroethylene, 1,1,2trichloroethane, chlorobenzene, and trichlorobenzene (total) are considered unnecessary or have not been established in the "Guidelines", while 1,2-dichloropropane is not specified in the "Standards". There are even more missing indicators in "WFD". The core of the water quality standards is the index and limit values. The index selection is primarily based on the water quality analysis to determine which pollutants are present in the water and their corresponding concentration levels. The limits are determined mainly based on toxicological data. The primary investigation requires significant time and manpower and material resources. Therefore, China's water quality standards are primarily based on the relevant standards adopted by the developed countries and organizations, such as WHO, Environmental Protection Agency (EPA) of the United States, and the European Union. Toxicological data can be used directly, but the contribution of different exposure routes and sources varies for countries and regions because of geographical factors, living standards, and habits, which require reassessment. Therefore, it is considered that the revised "Standards" should focus on the selection of indices. The indices that are unlikely to exceed the threshold (or where the apparent concentration is much lower than the healthbased guideline value) need not be included in the "Standards" for mandatory testing but can be used as reference indices for emergency management and response requirements. In addition, the "Standards" limits of F − , Nitrate nitrogen (NO 3 − -N), Nitrite nitrogen (NO 2 − -N), Hg, Nickel (Ni), tetrachloromethane, chloroethylene, and trichloroethylene vary from the "Guidelines" limits, of which 62.5% are stricter than the "Guidelines" limits. Therefore, the Chinese standard limits were adopted as the evaluation standards in this study, and for the values not included in the standard, the WHO guidelines were referred.
In accordance with the preliminary analysis of the data, 27 indices exceeded the limits, accounting for 71.05%. Each index of 101 sample wells was compared with Standards for groundwater quality (GT14848-2017), and 42 sample wells were classified as groundwater category 3 water bodies, 32 as groundwater category 4 water bodies, and 27 as groundwater category 5 water bodies, with no category 1 or category 2 water body. For the control standard, the exceeding rate of each sample index is shown in Table 2. The highest exceeding factor was NO 3 − -N, with an exceeding rate of 49.5%, followed by Fe and TH, with an exceeding rate of 42.57% and 31.68%, respectively. The exceeding rate of Cl − , TDS, NO 2 − -N, Mn, and 1,2-Trichloroethane was approximately 10.00%. The average exceeding rate of heavy metal pollutants was 2.48%, and that of organic pollutants was 2.50%. The concentrations of Na + , Cl − , SO 4 2− , TH, and TDS were primarily related to rock weathering, while the concentrations of NH 4 + -N, COD (Mn), F − , heavy metals, and organic matter were related to environmental pollution. In the single factor water quality evaluation method, if only one of all indices exceeded the standard value, the water quality is defined as exceeding the standards, which does not evaluate the groundwater quality accurately. This is because the impact (weight) of the index on the overall water quality is possibly minimal, and because of the different background water quality values, the water quality varies from place to place. Therefore, this study adopted the WQI method to comprehensively evaluate water quality.

Full Index Analysis
In accordance with Equations (1)-(3), the membership degree was calculated for 38 index data of 101 sample wells to obtain the decision matrix, In the light of Equations (4)- (7), the weights of 38 indexes were calculated as follows:

Delphi Method Analysis
In agreement with the Delphi method, we sent questionnaires to 60 experts in the water environment industry, and finally determined Na + , NH 4 + -N, Fe, Cl − , SO 4 2− , pH, TH, TDS, COD (Mn), Al, and Mn as the evaluation indices. To improve convenience for comparison, we also used the entropy value method to determine the weight in the Delphi method. According to Equations (1)-(3), the membership degree of 11 index data of 101 sample wells was calculated to obtain the decision matrix, In accordance with Equations (4)-(7), the weights of the 11 indexes are

Principal Component Analysis
The SPSS26 software was used for correlation analysis of the 38 indices, and the results showed that the concentration changes of these indicators were highly correlated. It means that these highly correlated index pollutants have the same or similar sources of pollution, and that they have the same change trends. So the concentration changes of one pollutant can be indicative the changes of other highly correlated pollutants [19]. We can screen the indicators through PCA.
Principal component analysis (PCA) and principal factor analysis (PFA) were carried out on 38 index data. PCA is a mathematical approach to dimensionality reduction. Through PCA, the original 38 indicators are recombined into several groups of new comprehensive indicators which are unrelated to each other to replace the original indicators. The information contained in each group of indicators is expressed by variance, that is, the greater the variance, the more information contained. Each set of indicators is called the principal component. The principal component 1 contains the most information, and then it decreases. In the process of principal component extraction, we select the principal components whose initial eigenvalue is greater than 1, namely, principal components 1 to 9. The total variance of the 9 principal components is presented in Table 3. The variance contribution rate of principal component 1 (PC1) was 29.188%, the variance contribution rate of principal component 2 (PC2) was 16.224%, and the accumulation variance contribution rate of the nine principal components was 85.356%. These values can better explain the characteristics of the 38 indexes, so 9 principal components were extracted. The composition matrix after rotation is shown in Table 4   In consonance with the PCA results, the indices with the highest correlation coefficient among the nine principal components were selected as evaluation indices, namely, Na + , Al, NO 3 − -N, Hg, tetrachloromethane, Mo, dichloromethane, xylene (total), and o-dichlorobenzene. In the light of Equations (1)-(3), the membership degree for nine index data of 101 sample wells was calculated to obtain the decision matrix, On the basis of Equations (4)-(7), the weights of the nine indices are w p = (0.095, 0.191, 0.120, 0.100, 0.090, 0.127, 0.100, 0.088, 0.087) Conforming to Equations (8) and (9), the WQI of 101 sample wells were calculated as follows:

Evaluation and Comparison
The groundwater was classified in agreement with the WQI evaluation results, as shown in Table 5 [48]. The WQI indices obtained by the full index method, Delphi method, and PCA analysis method were used for grading evaluation of 101 sample wells, as shown in Table 6. By the full index method, Class I sample wells accounted for 57.43%, Class II 21.78%, Class III 8.91%, Class IV 7.92%, and Class V 3.96%. In the Delphi method, Class I sample wells accounted for 28.71%, Class II 28.71%, Class III 17.82%, Class IV 8.91%, and Class V 15.84%. In the PCA method, Class I sample wells accounted for 74.26%, Class II 17.82%, Class III 5.94%, Class IV 0.99%, and Class V 0.99%. In all three methods, the average of Class I and II sample wells was 76.24%. Therefore, the groundwater quality in the Dawu water source area can be considered good and potable. Among the three methods, the full index method was the most accurate because it comprehensively considered the thirty-eight water quality indices and all the pollutants in groundwater. Therefore, the full index method was used as the benchmark for comparison. Figure 4 shows the comparison of WQI obtained by the three methods, with 101 sample well numbers on the abscissa and WQI values on the ordinate. The blue curve is the WQI calculated by the full index method; the orange curve is the Delphi method WQI, and the gray curve is the PCA method WQI. We line up the values at each point, which makes it easier to see the similarities between the various methods. As can be seen from the figure, the WQI value obtained by PCA method is more similar to that obtained by the full index method. A comparison between the Delphi method and the full index method by calculating shows that the same rate of WQI evaluation results of the Delphi method was 37.62%, which is worse than that of the full index method because the single-factor exceeding rate of the indices selected in the Delphi method was higher. The same rate of WQI evaluation results of the PCA method was 77.23%, which is better than that of the full index method. This is because the selected indices comprehensively consider the presence of organic pollutants; the exceeding rate of organic pollutants was very low, and therefore, the comprehensive evaluation results were good. The PCA method was more accurate than the Delphi method in terms of the same rate of evaluation results. Thus, we conclude that it is reasonable and scientific to apply the multivariate statistical analysis method in WQI technology.

Discussion
In this study, thirty-eight indices data from 101 sample wells in the Dawu water source area of China were analyzed, and the single factor method and WQI method were used to evaluate groundwater quality in this area. In the single factor water quality evaluation method, the weights of all indices are not considered. Even if only one of all indices exceed the standard and by a very small amount, the water quality is defined as exceeding the standard, and the evaluation results are one-sided. Therefore, this study adopted the WQI method for comprehensive evaluation.
In the index selection process of the WQI evaluation method, the full index method, the Delphi method, and the PCA method were adopted and the entropy method was used to determine the weight. Finally, the evaluation results were analyzed, and the scientific basis for the Delphi method and the PCA method were compared. Using the full index method, Class I sample wells accounted for 57.43%; using the Delphi method, these accounted for 28.71%; and using the PCA method, these accounted for 74.26%. Thus, from the three methods, the Class I and II sample wells had an average of 76.24%. Therefore, it may be concluded that the groundwater in the Dawu water source area generally has good quality and is potable. Selecting the full index method as the criterion, the same rate of WQI evaluation results of the Delphi method and the full index method was 37.62%, which was worse than that of the full index method. The same rate of PCA WQI evaluation was 77.23%, which was better than that of the full index method. Hence, the PCA method is more accurate than the Delphi method in terms of the same rate of evaluation results, and it is reasonable and scientific to apply the multivariate statistical analysis method and entropy weighting method in WQI technology.
Despite this, this study only carried out a comparative analysis on the detection data of the Dawu water source area; in the Delphi method, the selection of evaluation indices may consider more factors such as regional water quality background differences; more research is needed to determine whether the results of other groundwater or surface water assessments are consistent. At the same time, the large difference in the evaluation results between the full index method and the Delphi method also reminds scholars that they need to pay more attention to the rationality of the selection of indicators in the future research.

Discussion
In this study, thirty-eight indices data from 101 sample wells in the Dawu water source area of China were analyzed, and the single factor method and WQI method were used to evaluate groundwater quality in this area. In the single factor water quality evaluation method, the weights of all indices are not considered. Even if only one of all indices exceed the standard and by a very small amount, the water quality is defined as exceeding the standard, and the evaluation results are one-sided. Therefore, this study adopted the WQI method for comprehensive evaluation.
In the index selection process of the WQI evaluation method, the full index method, the Delphi method, and the PCA method were adopted and the entropy method was used to determine the weight. Finally, the evaluation results were analyzed, and the scientific basis for the Delphi method and the PCA method were compared. Using the full index method, Class I sample wells accounted for 57.43%; using the Delphi method, these accounted for 28.71%; and using the PCA method, these accounted for 74.26%. Thus, from the three methods, the Class I and II sample wells had an average of 76.24%. Therefore, it may be concluded that the groundwater in the Dawu water source area generally has good quality and is potable. Selecting the full index method as the criterion, the same rate of WQI evaluation results of the Delphi method and the full index method was 37.62%, which was worse than that of the full index method. The same rate of PCA WQI evaluation was 77.23%, which was better than that of the full index method. Hence, the PCA method is more accurate than the Delphi method in terms of the same rate of evaluation results, and it is reasonable and scientific to apply the multivariate statistical analysis method and entropy weighting method in WQI technology.
Despite this, this study only carried out a comparative analysis on the detection data of the Dawu water source area; in the Delphi method, the selection of evaluation indices may consider more factors such as regional water quality background differences; more research is needed to determine whether the results of other groundwater or surface water assessments are consistent. At the same time, the large difference in the evaluation results between the full index method and the Delphi method also reminds scholars that they need to pay more attention to the rationality of the selection of indicators in the future research.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.