Evaluation of Water Social Service and Comprehensive Water Management Linked with Integrated River Evaluation

Various factors like climate change and population increase have limited water management evaluation. In South Korea particularly, although the management of water quality and water quantity has recently been integrated, a comprehensive policy has not yet been identified. This study, therefore, aims to propose a methodology for evaluating water social service for 18 basins near major water resources in South Korea. It aims to promote advanced water resource management, secure water equity, and improve inadequate policy implementation. In addition, it proposes a methodology for comprehensive water management evaluation linked with integrated river evaluation with respect to water quality and water quantity. Accordingly, contrary to the common assumption that the entire population has easy access to the supplied water, the status of water service was assessed objectively. The status of water management per sector was also visually represented, through which the vulnerabilities of water management could be intuitively diagnosed. Based on the possibility of utilizing the study results to determine the basic direction for water management, the methodology of this study has been proposed as a tool for establishing an efficient water management policy.


Introduction
Climate change is influencing physical and biological changes on Earth [1], where a portion of the global population still lacks access to sufficient water. Approximately only 2.5% of the total volume of water on Earth can be used by humans, and accordingly, the concern for water scarcity is continuously rising as 2.1 billion people around the globe do not have access to safe drinking water [2]. Further, water supply insufficiency is expected to reach 40% by 2030 [3]. Therefore, the emphasis is being placed on managing water, which is an essential element for sustaining life on Earth and has been recognized as a critical resource [4]. Hence, many countries around the world have adopted water management policies and established relevant policy measures. For effectively managing vulnerable water resources, the sustainability of appropriate water quality, water quantity, water availability, and water supply must be secured. Thus, the concept of Integrated Water Resource Management (IWRM) was first introduced in the 1970s [5] and has been utilized as a decision-making tool for water-related policies [6]. In particular, after the IWRM guidelines were publicly announced during the fifth World Water Forum held in 2009 [7], countries around the world are establishing and promoting IWRM related policies.
Along with management of the quality and quantity of water resources, a balance between water usage and population growth as well as the preservation and protection of water resources for maintaining their functionality and characteristics are required [8]. This is a concept of water social service, which is in line with the concept of water welfare.
Here, water social service refers to all the citizens living in a healthy, safe, and pleasant environment regardless of their social position or residential area. It is similar to the concepts of the human right to water [9,10], environmental equity [11], environmental fairness [12], and environmental justice [13][14][15]. In other words, it guarantees a better quality of life with respect to water resources, ensures equal distribution of environmental resources and services benefits, protection from environmental pollution, equal opportunities to participate in policy-making, and a fair distribution of policy outcomes.
Previous studies have proposed water stress index [16,17], water scarcity index [18], and water poverty index developed for sustainable water usage [19]. Furthermore, United States Environmental Protection Agency (EPA) [20] has proposed an evaluation method for social equity in terms of water environment service and water infrastructure; Hoffjan et al. TRansitions to the Urban water Services of Tomorrow (TRUST) [21] also proposed a method for evaluating the financial sustainability related to water supply and wastewater treatment using the Financial Sustainability Rating Tool. In a study on improving environmental convenience and equity, Comber et al. [22] reported a difference in the quality of service received according to the social status based on a network analysis of the accessibility to public goods services by groups with the same ethnicity and religion. The Organisation for Economic Co-operation and Development (OECD) [23] proposed using a principal component analysis method to objectively calculate the weights of comprehensive indicators, and Gómez-Limón and Sánchez-Fernandez [24] used a principal components analysis and analytic hierarchy process (AHP) to assess the sustainability of agriculture. Nguefack-Tsague et al. [25] used principal component analysis to analyze the importance of detailed indices that constitute the human development index (HDI). In South Korea, the management of water quantity and water quality, which was differentiated, has been integrated since 2019, and the "Basic Water Management Plan" is being established based on the IWRM. However, the basic water management plan at the national level still remains only a declaratory concept, and ambiguous water management standards have resulted in confusion in adopting water management policy measures. In addition, no evaluation system can simultaneously diagnose water quality and water quantity before implementing a water management system. In a related study, Hwang et al. [26] developed an integrated river evaluation method in terms of water quality and water quantity for the Paldang watershed encompassing the seven cities and guns surrounding the Paldang Lake located near Seoul (South Korea), thus proposing a methodology for simultaneously evaluating the water quantity and water quality of each watershed. However, most of the aforementioned previous studies evaluated water management conditions based on only one field, and also these studies have a difficulty which could not lead the result of the comprehensive evaluation with water-related components or sector-specific review contents. In other words, the comprehensive evaluation of water usage could result in the continuity of a stable human life. However, either the water management evaluation is biased, or the methodology of evaluating only certain fields could hinder the current situation diagnosis and future water usage planning.
This study, therefore, proposes a methodology and relevant results that could simultaneously consider multiple water management fields to determine comprehensive water management status, providing current water management information, intuitive problem recognition, future water management objectives, and planning to policy decision-maker and local residents. To this end, this paper aimed to evaluate the water social services while implementing a comprehensive evaluation of water management linked with water social service and considering the integrated river evaluation ratings proposed by Hwang et al. [26]. In addition, this paper proposes a methodology to devise a water management application plan applicable to regions having geographies similar to South Korea.

Study Area
The Paldang watershed, which is a major water resource in South Korea, was set as the study area, as shown in Figure 1. This water resource is used by 50% of South Korean population, which makes it a rare case. Hence, the central government of South Korea is enforcing very strict site regulations to ensure the security and stability of Paldang Lake. The inflow rivers that affect the water quantity and water quality of the Paldang Lake include the Bukhan River, Namhan River, and Gyeongan River, which are managed by the central government owing to the significance of the Paldang Lake. Moreover, the five cities and two counties of the Paldang watershed surrounding the Paldang Lake are under strict site regulations, including the water source conservation area, buffer zone, and special water preservation area. Thus, the majority of these consist of forest areas or typical rural areas. Therefore, this study examined 18 basins of the 3 inflow rivers within the special water preservation area and considers their significance (Table 1).

Data Collection
It is more efficient to conduct research only on the items with the highest relevance because of the vast data available on water social service. Previously, water social service has been defined as a concept of welfare, such as the human right to water, environmental equity, and environmental fairness. All the individuals are guaranteed access to water resources, environmental resources, and services benefits are equally shared, and they are fairly protected from environmental pollution. In a related study, Sullivan et al. [27] proposed the water poverty index for assessing the degree of water shortage by comprehensively examining human socioeconomic activities, environmental factors, and the physical availability of water. This was based on five components: the number of water resources, accessibility to water resources, water resource management capability, the efficiency of water usage, water quality, and ecological environment. Furthermore, the United States Department of Agriculture (USDA) has defined the concept of water management as controlling and moving water resources for minimizing the damage to human life and property and maximizing the efficiency and benefit of water usage. It further mentioned that water resource management occasionally involves changing common practices such as assigning water resources for different purposes or underground water recovery [28]. In the International Decade for Action "Water for Life" (2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015), proclaimed by the United Nation (UN), an integrated water resource management accepted worldwide was proposed to cope with the demands conflicting with the efficient, fair, and sustainable development and management of limited water resources [29].
Based on previous studies, the following 11 items were chosen as the research data in this study for selecting the indicators of water social service: the daily water consumption per person for the entire population, elderly population, beneficiaries of national basic livelihood, water supply amount, water usage amount, regional water treatment plant capacity, planned distribution of multi-regional water supply, water supply, and sewer service ratio, reclaimed water usage amount, cost recovery rate of water supply and sewage service, water and sewage price, and construction cost and maintenance cost of a water supply and sewage system in the 18 basins. The relevant data were collected from the National Water Management Information System (WAMIS) website [30]. The aforementioned data provided by WAMIS is collected by the local and central government. The South Korean government has determined the population aged 65 or older as an elderly population and classified the beneficiaries of national basic livelihood according to income standards. The South Korean government has been managing the elderly population and beneficiaries of national basic livelihood through the database in order to implement a reduction policy of sewage rent in terms of welfare policy. Naturally, the South Korean government has been collecting the water usage information of them, and the related data has been provided online for national operation and research purposes. The research period was set for eight years, from February 2009 to December 2016. For the data on integrated river evaluation rating with respect to water quantity and water quality for implementing a comprehensive evaluation of water management, the results of the study by Hwang et al. [26] were used.

Methodology
The research methodology is explained in Figure 2. First, the water social service index (WSSI) was selected to evaluate the water social service; relevant indicators were then calculated. The final indicators were selected after examining the adequacy of each indicator, and the equation for calculating the WSSI was developed. Here, to avoid the indicator factors influence each other, a weighted value was calculated and reflected in the WSSI equation. Subsequently, the ratings of water social services were divided into five categories using the indices deduced previously. Finally, the comprehensive evaluation of water management was performed by linking integrated river evaluation ratings and WSSI ratings. The detailed research procedures are further described in Sections 3.1-3.7; all the statistical analyses were performed using IBM SPSS Statistics 22 (New York, NY, USA).

Indicator Selection for Water Social Service Evaluation
In this study, the 15 indicators for four water social service categories shown in Table 2 were selected for evaluating the water social service. The water supply capacity category includes the three indicators as revenue water ratio, water supply amount, and water supply and demand ratio, which indicates the efficiency of the water supply capacity from the perspective of the water supplier. Where, "supply" means the distribution from water supply facilities to consumers. Further, "demand" means the required water quantity from predicted data of the national water supply plan. The water usage category consists of four indicators including the daily water usage per person, amount of replacement water, supply ratio of water system, and supply ratio of sewer treatment system. Where "usage" means actual use from the water consumer. These indicators imply the quantity and a variety of water usage being used from the user's perspective. The water service management category consists of the cost recovery rate of water supply and sewage as well as the self-reliance ratio of water supply and sewage as detailed indicators, which imply the management level of water services. Lastly, the water service equity category consists of the water supply and demand ratio and water usage amount of the elderly population as well as the beneficiaries of the national basic livelihood, through which the difference in the usage of water service according to socioeconomic characteristics is examined from the perspective of water justice. Rate of water charge to the unit cost of water supply production;a high indicator value implies that the water service management is adequate CRWS Cost recovery rate of sewage service Rate of sewage service charge to the cost of sewage treatment; a high indicator value implies that the water service management is adequate CRSS Self-reliance ratio of water supply Ratio of annual revenue to annual expenditure (construction + maintenance costs)/water charge; it represents the self-reliance of water supply service SRWS Self-reliance ratio of sewage service Ratio of annual revenue to annual expenditure (construction + maintenance costs)/water charge; it represents the self-reliance of sewage service

SRSS Water Service Equity
Water supply and demand ratio of the elderly population

Building Indicator's Data
The 15 indicators in the four categories for evaluating water social service can be calculated as follows.
Here, water demand amount (WDA, ton/year) is the product of population and water supply amount per person (WSAP, ton/year), whereas the WSA (ton/year) is the sum of regional water treatment plant capacity (ton/year) and planned distribution of multi-regional water supply (ton/year).

ARW (ton/year) = ASR (ton/year) + RWAS (ton/year)
(3) Calculating the indicators of water service management category Water service management category includes CRWS (%), CRSS (%), SRWS, and SRSS. The data from WAMIS [28] were used for CRWS and CRSS, whereas the SRWS and SRSS, which signify the independence of each service, were calculated using Equations (4) and (5), respectively. SRWS (non-dimension) = (construction costs + maintenance costs)/water charge (4) SRSS (non-dimension) = (construction costs + maintenance costs)/sewage charge (5) (4) Calculating the indicators of water service equity category Water service equity category includes WREP (%), WUEP (ton/year), WSDB (%), and WUSB (ton/year) in which each index was calculated using Equations (6)- (9). In Equations (6)-(9), water demand of the elderly population (WDEP ton/year) and water demand of the beneficiaries of national basic livelihood (WDBP, ton/year) is the product of the relevant population and the water supply amount per person. The amount of water supply to the elderly population (AWEP) and the amount of water supply to the beneficiaries of national basic livelihood (AWBP) was calculated as a ratio of the relevant population to the sum of regional water treatment plant capacity and planned distribution of multi-regional water supply, respectively. Furthermore, the data of daily water usage per elder (DWUE), daily water supply amount per elder (DWAE), daily water usage per beneficiary of national basic livelihood (DWUB), and water supply amount per beneficiary of national basic livelihood (DWAB) were taken from WAMIS [30].
WSDB (%) = DWUB/DWAB Here, WDEP and WDBP is the product of the relevant population and the WSA per person. The AWEP and AWBP were calculated as a ratio of the relevant population to the sum of regional water treatment plant capacity and planned distribution of multi-regional water supply, respectively. Here, the regional water treatment plant capacity means the water supply ability of water supply facilities built by a local government. The planned distribution of multi-regional water supply means the data of water supply project planned by a metropolitan government. This data has been offered on the internet.

Normality and Standardization of Each Indicator's Data
The raw data are assumed to have a normal distribution for statistical analyses; thus, the range of data has to be matched, or the data distribution needs to be modified to be similar for performing statistical analyses of the collected data. In addition, data must be standardized because different dimensions, scales, or units of the data collected for comparing the results of statistical analyses can generate errors when interpreting the results.
In this study, skewness was used to examine the normality and bias of the indices collected or calculated for each detailed indicator in Section 3.2. In this study, if the skewness of each detailed indicator was between −1 and 1, it was considered that there was no bias. If the skewness was below −1 or greater than 1, the normality was secured by converting with the common log and the natural log, respectively. Then, the standardization was performed to match the scale of the data.
For the standardization of data, a difference from the baseline, scale readjustment, Z-Score, categorical scale, percentage of annual difference over consecutive years, or cycle index can be used [31]. Considering the characteristics of this study, all index values must have the same directionality; however, inverse property of certain data can lead to errors in data interpretation. Thus, the methods mentioned above are not appropriate for standardization in this study. As an alternative, the minimum-to-maximum standardization method that can match the directionality is proposed in which the interpretation error, because of an inverse attribute of the data, can be prevented. Accordingly, the minimumto-maximum standardization method was adopted in this study, as shown in Figure 3, to standardize the index of each detailed indicator.

Indicator Selection and Reliability Review
The validity of a measurement tool must be verified because a valid and reliable measurement tool is an essential component of any study, and the interpretation of the findings depends on the validity of the measurement tool [32]. This study tried to use the previous study results to review the construct validity of whether 15 indicators are suitable as tools for evaluating water social services. However, few studies have been conducted in Korea to evaluate water social services, and no basic structure of data has been established for water social service evaluation. Therefore, this study was intended to determine the construct validity by conducting a factor analysis and remove indicators that cannot be interpreted or have low reliability. In general, factor analysis is divided into two categories: Exploratory Factor Analysis (EFA) and Confirmatory Factor Analysis (CFA), depending on the purpose. EFA is used by researchers when there is no definite information about any factors and number of factors, or when the underlying structure of the data is unknown, and CFA is used to hypothesize existing relationships between variables and demonstrate their relationship through factor analysis.
As mentioned earlier, this study judged that EFA is an appropriate method to review the construct validity for selected indicators, as it is difficult to find studies of water social service evaluation in South Korea. Principal Component Analysis (PCA) and Maximum Likelihood (ML) could be suggested as the EFA's factoring method. PCA and ML are both methods of summarizing or dimension reducing for multiple variables; however, factoring using PCA could lead to false results of larger factor loadings [33]. ML is introduced as a reasonable method to extract common factors under the assumption of the multivariate normal distribution as a method used when the normality of the variables is mostly satisfied [34]. For this reason, this study performed ML used as the factoring method of the EFA to review the construct validity of indicators.
Before conducting the EFA, the data's adequacy has to be measured. To this end, this study had reviewed the Kaiser-Meyer-Olkin (KMO) and Bartlett's test, and communality. In general, it is considered appropriate to conduct EFA when the KMO value is greater than 0.5, Bartlett's test has p (probability) less than 0.1, and communality is 0.5 or greater [35][36][37][38]. Further, the number of components was determined by extracting only the components with an eigenvalue greater than one by Keiser-Guttman method [39,40].
Factor rotation methods for factoring can be suggested with orthogonal rotation and oblique rotation, which ML typically was used the direct Oblimin method, one of the oblique rotation methods. In this study, the direct Oblimin method was used as a factor rotation method. Thus, this paper has reviewed the reliability of extracted indicators using Cronbach's alpha. Cronbach's alpha was used to review the internal consistency of indicators. Generally, confidence for application research is obtained when Cronbach's alpha is more significant than 0.8 [41].

Index Equation Design for Evaluating Water Social Service
This study proposed an additive function with the simplest form to calculate the WSSI using the indicator mentioned above's value. The result of using an additive function has the risk of overestimating when the value of each item is summed. To resolve this issue, an appropriate weighted value could be applied to each indicator; the additive function applied with a weighted value can be represented as shown in Equation (10), in which the sum of weights is one.
Here, W i denotes the weighted value of each indicator, and I i presents a standardized indicator value of each category.

Calculating the Weighted Value
Both, the ML and entropy weighted method (EWM), were performed for each indicator in this study to ensure that the weighted value does not allow the indicator factors to influence each other. It should be noted that there could be a compensatory issue owing to the cancellation effect of a low indicator value being diluted by a high indicator value. To prevent such an issue, this study was used a geometric mean to accumulate the indices. The weighted value with the smallest standard deviation was selected among all the values derived through the ML and EWM for calculating the WSSI and ensure that all indicators have an even influence on the outcomes.

Calculation of Weighted Value Using Principal Component Analysis Method
Factor analysis could be described as a technique for defining potential factor structures that describe the correlation structure between variables. The primary purpose of a factor analysis is to explore and identify the latent variables underlying the relationship between variables, i.e., the structure of the factors. This study had intended to calculate the weighted value using the characteristics of factor analysis. To this end, this study adopted the principal components analysis (PCA) method, which was reported by Hotelling [42].
This method, along with Factor analysis, is widely used to analyze qualitative indicators in the social and economic field [43][44][45].
The method of calculating weighted value using PCA method requires component matrix and eigenvalue, which in this study used the structure matrix and eigenvalues derived by ML for them. This study judged that the "structure matrix" is suitable for weighted value calculation because this matrix has elements which are a correlation coefficient that refers to the relationship between factors and variables.
The calculating weighted value using PCA method can be described as follows. When a matrix of k explanatory variables and n observed value is X, the linear combination Z can be expressed as follows: where A is an arbitrary symmetric matrix.
Considering the first linear combination, it can be expressed as follows: Here, when Lagrangian is applied to obtain the largest eigenvalue, it is shown below.
If the first eigenvalue is λ 1 , λ 1 will have the largest eigenvalue of (X X). The weighted value can be expressed as follows: Therefore, when λ 1 is identified, the weighted value can be calculated. Meanwhile, the size of X ij is determined by

Calculation of Weighted Value Using Entropy Weighted Methods
The entropy weighted methods (EWM) [38] is a method for determining the weighted value of each attribute. The method involves applying simple information theories to decision-making problems having a great number of attributes and alternatives. This particular method can be used for selecting the weighted value based on the attribute information of the evaluation items and for applying the weighted value in a vulnerability assessment. For calculating the weighted value, the attribute information of each evaluation item is extracted and normalized to calculate the entropy of each attribute information and the weighted value of each evaluation item.
The EWM has the advantage of being objective because the weighted value is calculated by considering the distribution of the secured data. Four basic methods used to measure information based on entropy include maximum entropy, joint entropy, conditional entropy, and trans-information. Maximum entropy is derived from statistical physics research and it is a technique that excludes assumptions about unknown facts and sees only known partial facts as a source of knowledge acquisition. This technique can function as a probability model with a uniform distribution for other values if the predefined constraints are met. Maximum entropy is defined as the entropy amount of the probability p(x n ) distribution of the discrete random variable X = {x n }. Maximum entropy H(X) refers to the uncertainty or information capacity of X, and limit entropy for hydrological random discrete variables is defined as Equation (11).
Generally, there is the maximum entropy, if the discrete probability variable X follows a uniform distribution. Further, having a uniform distribution of the maximum entropy model for unknown facts means that the model has maximum entropy for unknown contents. This is considered to be the same probability for two events when there is no information to clearly distinguish a given event; this means that uncertainty is maximized, resulting in maximum entropy.
For example, if y m = (m = 1, 2, · · · , N) is related to random x n , H(X|Y) represents the amount of X information that cannot be represented by the information in Y and it is expressed as Equation (12).
p(x n , y m )lnp(x n |y m ) (17) where the probability of combining p(x, y) : X = {x n } and Y = {y m } is p(x n |y m ) , and the conditional probability H(X|Y) of a given and is the conditional entropy of a given Y and X, which can be interpreted as a loss of information during the process of information transfer between X and Y [46]. The brief calculation process of the weighted value by EMW can be summarized as follows: Step 1: Constructing the attribute matrix using the data of indicators of each watershed.
Step 2: Normalizing the data of indicator.
Step4: Calculating weighted value for each indicator.

Assigning Ratings to Water Social Service Index
The rating of WSSI was divided into five sections. The calculated index values were assumed to have a Gaussian distribution in this study, and the ratings were given based on the standard deviation. The rating classification based on the standard deviation was determined as proposed in the study by Hwang et al. [26] in which 1 standard deviation (1 σ) was defined as the highest and the lowest ratings. Further, 16% on both left and right ends were defined as the highest rating of 1 and the lowest rating of 5, respectively. The section corresponding to the remaining 68% was divided as 0.5 standard deviation (0.5 σ), 67.5-84.0% was classified as rating 2, 32.5-67.5% was classified as rating 4, and 16.0-32.5% was classified as rating 5 (Table 3).

Comprehensive Evaluation of Water Management Based on Matrix Analysis
Another objective of this study is to perform a comprehensive evaluation of water management linked with evaluation results of water quantity and water quality in addition to evaluating water social service. Accordingly, the matrix method was used for the analysis in which the data are visualized by cross-linking two lists consisting of rows and columns. This method examines the correlation of data using a matrix consisting of two or more dimensions [47].
In this study, a comprehensive evaluation of water management was performed using the matrix shown in Figure 4. The x-axis represents the rating of integrated river evaluation, whereas the y-axis represents the rating of water social service. Number 3 on both the xand y-axis indicates Acceptable or rating 3 of the evaluation items and as the values increase along the x-axis and the y-axis, the rating becomes Poor; and as they decrease, the rating becomes Good. As shown in Figure 4, this matrix is a 2 × 2 matrix that can be divided into 4 sections; however, the matrix is divided into Area I, Area II, and Area III in this study. Area 1 indicates that both integrated river evaluation and water social service evaluations have good ratings where the result of a comprehensive evaluation of water management has a rating I. Two Area IIs indicate that either the integrated river evaluation or the water social service evaluation has a poor rating where the result of a comprehensive evaluation of water management has a rating II. Area III indicates that both integrated river evaluation and water social service evaluations have poor ratings where the result of a comprehensive evaluation of water management has a rating IV. Rating III is either not represented as an area in the matrix, and it is used when the value of both xand y-axis is 3; it is also the point of intersection in the center in Figure 4. This rating indicates that both integrated river evaluation and water social service evaluations have acceptable ratings where the result of a comprehensive evaluation of water management has a rating III. Thus, the comprehensive evaluation of water management has four ratings where smaller the rating number, better is the rating.

Suitability of Collected Data
As the calculation results of 15 indicators' data by Figure 3 during the study period cannot be fully presented due to the limited space, only the normalization and data conversion results of each detailed indicator' data are shown in Table 4. As mentioned in Section 3.3, skewness was analyzed to review the normality of each detailed indicator' data, and the results are shown in Table 4 ((a) Skewness before Conversion). The skewness of all detailed indicators' data is greater than 1 in Table 4 ((a) Skewness before Conversion), thus implying that the data is biased. To resolve this issue, the skewness of each detailed indicators' data was converted using a natural log to produce the results presented in Table 4 ((b) Skewness after Conversion). The results show that the skewness after conversion of all detailed indicators' data is between −1 and 1, which implies that the bias has been removed and the normality of indicators' data secured.

The Result of KMO, Bartlett's Unit Matrix Test, and Communality
The adequacy of data must be reviewed before EFA performing for extracting common factors. To this end, the KMO, Bartlett's test and communality for adequacy review of data were conducted, and the results are shown in Table 5. The KMO measured the adequacy to be 0.879, and the Bartlett's test resulted in a p-value less than 0.001. The data are appropriate for EFA when KMO is above 0.5 [48], Bartlett's test's probability value is below 0.1. The KMO adequacy measure and Bartlett's significance probability satisfied the criteria. The communality ranged between 0.844 and 1.000, however, the communalities of WUEP and WUSB were presented greater than 1.0.
If the communality is greater than 1.0 which is called the Haywood case [49]. When the Heywood case has occurred, the analysis results should not be accepted because the communality results are not good with the number of factors is not appropriate, the number of variables per factor is small, or the common factor model is not appropriate [50][51][52][53]. In this work, the indicators were removed and the suitability of the data was reviewed again.

Result of Explatory Factor Analysis
Two indicators with Haywood case were removed from the 15 indicators of Table 5, and a factor analysis was conducted again. KMO, Bartlett's test, and communality for 13 indicators were reviewed, KMO represent at 0.851 and Bartlett's p-value below 0.001 (Table 6). Further, the communality of all indicators ranged 0.736-0.998. Therefore, this study determined 13 indicators are suitable for factor analysis. The EFA of this study was employed the maximum likelihood as factoring method and adapted only the factor with an eigenvalue greater than 1.0 by the Kaiser-Guttman method. As the factoring result, the two factors were extracted that eigenvalue of factor 1 and factor 2 presented 10.500 and 1.606, respectively ( Table 7). The criterion for explained variance ratio and cumulative variance ratio of each factor must be specified at least what contribution should be to have practical significance. In this regard, Hair et al. [54] suggested that the social science field's cumulative explanatory variance should usually be 50-60%. The extraction sum of squared loading of two factors was shown to be 90.386% in Table 7; therefore, this study determined these factors have high explanatory power. The pattern matrix, structure matrix, Cronbach's alpha and Goodness of fit test for 13 indicators according to the factor analysis were shown in Table 8. Pattern matrix and structure matrix can be used for interpreting of factor analysis result [55]. Some previous studies emphasize structural matrix [56][57][58], or others insist that pattern matrix is essential for interpreting factor analysis result. The matrix selection criteria of factor analysis are unclear; however, this study had adopted a pattern matrix used in many prior studies. Furthermore, the indicators selected only elements representing values greater than 0.4 and utilized them for factor analysis. [59,60]. In Table 8, the indicators of the pattern matrix were represented in order of magnitude. WREP, WSDB, SRWS, and DWUP among the pattern matrix's Factor 1, CRWS, SUWS, RWR, WSA, SRST, and CRSS among the indicators of the pattern matrix's Factor 2 were shown to be less than 0.4. Furthermore, Chi-squared and p-value of Goodness of fit test were denoted at 765.973 and less than 0.001, respectively. Consequently, this study judged that factor analysis was conducted appropriately.
By comparing Factor 1 and Factor 2, the indicator was classified into two Factors: Factor 1 included seven indicators (CRWS, SUWS, RWR, WSA, SRST, CRSS, and ARW), and Factor 2 included six indicators (WREP, WSDB, SRWS, DWUP, and SRS). Thus, this study had denominated Factor 1 as "Water Infrastructure" and Factor 2 as "Water Safe" based on the characteristics of each Factor.
The Cronbach's alpha for each indicator was represented to determine the reliability of Factor in Table 8. The Cronbach's alpha of Factor 1 and Factor 2 appeared at 0.978 and 0.979, respectively. In general, confidence for application research is obtained when Cronbach's alpha is more significant than 0.8 [41]. As shown in Table 8, all indicators and the overall Cronbach's alpha appeared higher than 0.8, ensuring reliability without indicators that undermine reliability and maintaining internal consistency.

Calculation Results of the Weighted Value of Each Indicator
PCA method and EWM method were performed to calculate the weighted value of each indicator; the results were represented in Table 9. The weighted value and rank of PCA method showed that the weighted value for WREP was the highest at 0.0886, followed by the CRWS and WSDB at 0.0879 and 0.0875, respectively. Lastly, the indicator with the lowest weighted value, 0.0494, was the ARW. The weighted value and rank of EMW appeared differently from PCA method. DWUP's weighted value of EMW method was the highest at 0.0922, SRWS's weighted value was the lowest at 0.0475. As mentioned in Section 3.5.2, the weighted value with the smaller standard deviation between the PCA method and EWM method was adopted as the weighted value for WSSI estimation. Table 9 shows that the standard deviation of PCA method is 0.0113, whereas that of EWM method is 0.0150; consequently, the weighted value by PCA method was employed as the weighted value for WSSI estimation.

Water Social Service Index
Using the method explained in Section 3.5, an additive function was used to calculate WSSI as shown in Equation (12).
WI (Water Infrastructure) = W 1 I 1 + W 2 I 2 + W 3 I 3 + W 4 I 4 + W 5 I 5 + W 6 I 6 + W 7 I 7 (23) WW (Water Welfare) = W 8 I 8 + W 9 I 9 + W 10 I 10 + W 11 I 11 + W 12 I 12 + W 13 I 13 (24) Here, WI refers to the water infrastructure (factor 1), WW refers to the water welfare (factor 2). Furthermore, W 1 -W 13 represent the weighted values of each indicator in Table 9, whereas I 1 -I 13 represent the standardized indicators of each watershed. The WSSIs were calculated every year for each watershed using these values, and the WSSI ratings were given.

Water Social Service Index per Watershed
The mean WSSI during the study period is shown in Table 10. The average of mean water infrastructure, water welfare and WSSI for all study areas appeared at 0.2740, 0.2623, and 0.5309, respectively. The median of mean water infrastructure, water welfare, and WSSI for all study areas appeared at 0.3085, 0.2691, and 0.5728, respectively.   Figure 5 shows that the two indicators of water social service are higher when the area of the watershed is smaller, which could be attributed to the high efficiency of water-related infrastructure. In contrast, the index values are lower when the area of the watershed is larger, which indicates the low efficiency of water-related infrastructure. The study area where Paldang Lake is located is the tough sanction area. Although the study area is the metropolitan area of Korea, the Korean central government has applied very strong regulations to the study area because the Lake Paldang is a very important water resource. For this reason, the population density of study area is very low compared to other metropolitan areas, and local residents live in a large area with the low population density. Due to these characteristics of the study area, the Korean government has been also installing water-related infrastructure, however, in large areas, infrastructure efficiency is inevitably lower than in narrow areas. Although no direct numerical evidence is suggested in this study, the results of this study are relevant considering infrastructure efficiency, which is one of the general characteristics, because the efficiency of infrastructure might be low in regions with strict regulations and slow development, such as a case of this study area.

Evaluation Results of Water Social Service Index Rating
The ratings of water infrastructure, water welfare, and WSSI are presented during the study period in Table 11 and Figure 6. GJ and SD watersheds had a rating of 1 for water infrastructure and water welfare, which is considered very good, whereas BH, YH, and GJA had a rating of 5 for water infrastructure, and JJ, HC, and BH had a rating of 5, which is considered very poor. As shown in Table 11, the mean ratings of all the indicators for each watershed are similar without being biased. For example, all the mean ratings of GJ and SD, which were given rating of 1 for WSSI, are either 1, whereas all the mean ratings of HC, BH, YH, and GJA, which were given rating of five for WSSI, are either four or five. Specifically, it indicates that there is a substantial gap in the WSSI between the different watersheds. As the mean ratings are similar for all the WSSI, the watersheds with poor WSSI need improvements in all the sectors rather than in a specific sector. Figure 6 shows the changes in the annual rating of water infrastructure, water welfare, and WSSI. Most of the WSSI did not vary significantly except MG watershed. These results may correspond to the efficiency of water-related infrastructure mentioned previously in Section 4.5.2. However, the lack of change in WSSI over a long period also signifies a lack of awareness of water social service in each watershed. In particular, we have attributed this result to the indifference to the underprivileged of local government who are responsible for the installation and operation of water-related infrastructure. Another reason for the lack of change in the WSSI is low community awareness among the local residents. From the perspective of human right on water or water justice, a prolonged low WSSI signifies that the elderly population or the beneficiaries of national basic livelihood have been experiencing discrimination or limitations of water usage for the duration study period or longer. Particularly, the WSSI of JJ, GJA, and HC watersheds have remained at 4 or 5, which indicates that the residents in these respective regions have not been receiving universal welfare benefits in terms of water-related services for an extended period.
Accordingly, the water management evaluation based on the WSSI from the perspective of human right to water and water justice enables the status of universal water welfare to be assessed objectively, contrary to the vague expectations that all the users have access to the supplied water without any difficulty. Therefore, the results can be utilized as the basic materials for setting the future direction of advanced water management. Moreover, plans can be established for each sector depending on the results of water social service indicators, and thus a more efficient water management can be gradually designed.  Figure 7 shows the results of a comprehensive evaluation of water management of the study area based on a matrix analysis. The x-axis represents the rating of the integrated river evaluation, whereas the y-axis represents the rating of water social service. As previously mentioned in Section 2.2, the findings of the study by Hwang et al. [26] were used as the data for representing the x-axis. The results of comprehensive evaluation of water management are depicted as blue, black, and red arrows to represent the changes over time. The blue arrow represents the first movement, black arrow represents the following movement, and red arrow represents the last movement. As shown in Figure 7, results of the comprehensive evaluation of water management move horizontally with the flow of time. It can be inferred that most WSSI ratings are constant over time, such as Figure 6, and the horizontal movement is because of the changes in the ratings of integrated river evaluation. When the changes in the ratings of a comprehensive evaluation of water management are examined, there is no abrupt change in the comprehensive evaluation of water management of each watershed. In particular, the rating moved to an adjacent area of the 2 × 2 matrix or stayed within the area where the first comprehensive evaluation of water management result was given. Especially, in the 11 watersheds including GU, BG, ST, HC, YD, BH, YH, OS, MHy, GJA, and MG among the total 18 watersheds in this study, the water management level did not improve over time and reverted to the original rating.

Result of Comprehensive Evaluation of Water Management Based on Matrix Analysis
Furthermore, the WSSI of MG was one; however, their final integrated river evaluation rating was 3 or 4, which resulted in the comprehensive evaluation rating of water management to be rating II. The three watersheds BH, YH, and GJA have a poor rating for both the integrated river evaluation and the water social service thus having the comprehensive evaluation rating of water management placed in Rating IV. Hence, even if a watershed has a good water social service rating, the comprehensive evaluation of water management may be low owing to the poor rating of the integrated river evaluation. The watershed with a poor water social service rating also has a poor integrated river evaluation rating, thus having a low comprehensive evaluation rating of water management.
The comprehensive evaluation of water management using a matrix enables comprehensive monitoring of water management status of each watershed, and thus is an evaluation tool or a minimum indicator for water management of watersheds. The comprehensive evaluation of water management by matrix analysis enables the vulnerabilities in water management to be identified intuitively, unlike the other evaluation systems that are biased or concentrated on one side, thus allowing the policy makers or administrators of water-related matters to recognize the changes in water management status over a period. Moreover, the direction of complex water management can be set by visually representing the susceptible aspects of water management. The regular inspection results can be easily presented in a matrix, thus facilitating an immediate decision-making process.
However, the comprehensive evaluation of water management by matrix analysis only uses two axes; therefore, it has a limitation of using two criteria for evaluating the status of complex water management.

Conclusions
This study proposes a method to evaluate water social service. The water social service of each watershed in the study area was evaluated, and the vulnerable aspects in terms of the human right to water and water justice were identified. Moreover, the status of water social service over a period could be examined to address the low awareness of water welfare among the local governments and residents. In addition, a method of comprehensive evaluation of water management using a matrix analysis was proposed through which comprehensive monitoring of water management is possible. However, if some country has no source data or difficulty of using data, it might be challenging to produce results such as this study. This case could be explained as the limitations of this study though it is not only the fields considered in this study intended to judge water management comprehensively. Nevertheless, if the methodology proposed in this study is used using data that could be used depending on the circumstances or conditions, this paper believes that comprehensive water management evaluation of the watershed (or countries) will be possible.
The methodology proposed in this study has certain limitations, such as employing only two axes in matrix analysis; however, it can be expected to prevent discriminations in water usage against the disadvantaged group by evaluating the water social service to promote a society in which universal water welfare is secured. The comprehensive evaluation of water management using a matrix can help decision makers or administrators to recognize the changes in water management status over a period, and thus lead to efficient water management. Furthermore, the proposed methodology can promote the establishment of a systematic and comprehensive water management system at a national level, while being utilized as an alternative for water management in regions having a similar topology as South Korea.