Application of Factor Analysis for Characterizing the Relationships between Groundwater Quality and Land Use in Taiwan’s Pingtung Plain

: Although the average municipal water coverage in Taiwan is quite high, at 93.91%, only around half of the residents in the Pingtung Plain use tap water originating from the Taiwan Water Corporation to meet their needs. This means the exploitation of a substantial amount of groundwater as a source of water to meet drinking, agriculture, aquaculture, and industry requirements. Long-term groundwater quality surveys in Taiwan have revealed obvious contamination of the groundwater in several locations in the Pingtung Plain, with measured concentration levels of some groundwater quality parameters in excess of the permissible levels speciﬁed by the Taiwan Environmental Protection Administration. Clearly, establishing a sound plan for groundwater quality protection in this area is imperative for maximizing the protection of human health. The inappropriate use of hazardous chemicals and poor management of land use have allowed pollutants to permeate through unsaturated soil and ultimately reach the underlying shallow unconﬁned groundwater system. Thus, the quality of the water stored in shallow aquifers has been signiﬁcantly a ﬀ ected by land use. This study is designed to characterize the relationship between groundwater quality and land use in the Pingtung Plain. This goal is achieved by the application of factor analysis to characterize the measured concentrations of 14 groundwater quality parameters sampled from 46 observation wells, the area percentages for nine land use categories in the neighborhood of these 46 observation wells, and the thicknesses of four unsaturated types of soil based on core samples obtained during the establishment of 46 observation wells. The results show that a four-factor model can explain 56% of the total variance. Factor 1 (seawater salinization), which includes the groundwater quality parameters of EC, SO 42 − , Cl − , Ca 2 + , Mg 2 + , Na + , and K + , shows a moderate correlation to land used for water conservation. Factor 2 (nitrate pollution), which includes the groundwater quality parameters of NO 3 − -N and HCO 3 − , shows a strong correlation to land used for fruit farming and a moderate correlation to the thickness of the gravel comprising unsaturated soil. Factor 3 (arsenic pollution), which is composed of groundwater quality parameters of total organic carbon (TOC) and As, is very weakly affected by land use. Factor 4 (reductive dissolution of Fe 3 + and Mn 2 + ), which involves Mn 2 + and Fe 3 + , is weakly impacted by land use. Based on a geographic visualization of the scores for the four different factors and the patterns for land use, we can demarcate the areas where the groundwater in shallow unconfined aquifers is more vulnerable to being polluted by specific contaminants. We can then prioritize the areas where more intensive monitoring might be required, evaluate current land use practices, and adopt new measures to better prevent or control groundwater pollution.


Introduction
Groundwater is a vital and reliable component of water resource supplies around the world. In Taiwan, in areas where there is a shortage of surface water, residents depend upon groundwater to meet their basic daily water demands as well as for irrigation for agriculture and aquaculture. The Pingtung Plain in Southwestern Taiwan is an intensively productive agricultural area, mostly from the cultivation of crops and aquaculture. Approximately 50.5% of the land in the Pingtung Plain is used for agriculture and 5.5% for aquaculture, and agricultural activities have continued to intensify over the last few decades. The total groundwater used annually in the area is 24.3 million tons, with most being used for agriculture (63%), aquaculture (8%), industry and other demands (29%) [1]. Moreover, the residents of the Pingtung Plain are unusual, in that only around 50.8% use municipally supplied tap water even though the average tap water coverage in Taiwan is 93.91%.
Long-term investigations of groundwater quality in Taiwan have revealed obvious groundwater contamination, with the measured concentrations of some groundwater quality parameters at some monitoring wells in the Pingtung Plain in excess of the acceptable levels specified by the Taiwan Environmental Protection Administration [2][3][4][5]. The use of contaminated groundwater for drinking, irrigation, aquaculture, and industry can have potentially adverse effects on the health or activities of residents. Exposure to some hazardous chemicals is known to lead to a variety of acute and chronic health effects. For example, arsenic (As) is recognized as a toxicant and carcinogen which can cause Blackfoot disease as well as cancers of the liver, kidney, bladder, prostate, lymphoid tissue, skin, colon, lungs, and nasal cavity, ischemic heart disease, hyperpigmentation, hyperkeratosis, diabetes, and meningioma [6][7][8][9][10][11][12][13]. Arsenic can enter the food chain indirectly-for example, through the consumption of fish cultivated in As-affected groundwater or the ingestion of crops grown using As-contaminated groundwater for irrigation or directly by drinking As-contaminated water [14][15][16]. Nitrate is a naturally occurring form of nitrogen necessary for crop growth. However, decades of intensive farming in the same area can result in the leaching of excess nitrate from manures and fertilizers into shallow unconfined groundwater system [17][18][19]. Septic systems and cesspits are also sources of nitrates. All of this makes it one of the most common contaminants in groundwater worldwide. The ingestion of excess nitrate can cause blue baby syndrome, also known as methemoglobinemia, in infants, which can lead to brain damage and sometimes death (e.g., Liu et al. [20]). Both chloride and sodium are highly soluble chemical elements naturally found in groundwater, particularly in coastal areas. However, excess levels of chloride and sodium ions can not only adversely affect the taste of drinking water, rendering it unsuitable for drinking, but high concentrations of chloride and sodium ions in the groundwater can also have an inhibitive effect on the growth of crops. Iron and manganese are metals that occur naturally in soils, rocks, and minerals, but when groundwater comes into contact with these solid materials, they can be dissolved and, with their constituents, released into the water body. If groundwater containing high levels of iron and manganese is used for agriculture or aquaculture, it may inhibit the growth of cultivated crops or farm fish.
Achieving water security and availability for people worldwide is among the principal agenda for the Sustainable Development Goals adopted by all United Nations Member States. To achieve the goal of water security, establishing a sound plan for safe and sustainable groundwater quality management is imperative for protecting human health, ensuring food safety, maintaining access to clean and healthy water resources, and preserving the environment. Although they may be laborious, expensive, time-consuming, and complex, efforts for characterizing groundwater quality are required as a critical step in identifying good and safe water quality. The chemical composition of the groundwater in an aquifer is a direct consequence of the composition of the water that enters the aquifer and the interaction of this groundwater with the surrounding minerals deposited by various natural processes and/or anthropogenic activities. Such effects may significantly modify the various chemical compositions of the groundwater. Groundwater contamination originating from the incidental release of various types of wastes or the inappropriate management of materials used and products generated by industrial, agricultural, and public activities on the land surface is carried along with water infiltration through unsaturated soil, ultimately reaching the water in the underlying aquifer [17]. Groundwater in shallow unconfined aquifers is even more vulnerable to pollution by chemicals derived from anthropogenic activities. Moreover, the close connection between these shallow unconfined aquifers and the overlying land surfaces further supports the premise that groundwater quality is especially affected by land use.
Considering the obvious necessity to protect groundwater sources from contamination through human activities, studies on the relationship between groundwater quality and land use have gained increasing attention over the past couple of decades. One approach is simply to make a comparison between the contaminants in the groundwater and land uses in specific areas. For example, Wang [21] examined the relationship between groundwater quality and land use in Rhode Island in the United States using linear and multiple regression analysis. He found that the areas of residential land use were related to increased levels of chloride and sodium concentrations in well water. Nitrate concentrations were also closely related to increased residential land use. Eckhardt and Stackelberg [22] applied maximum-likelihood logistic regression analysis to investigate the relationship between groundwater quality and land use in Long Island, New York, in the United States. They considered five different study areas and several groundwater quality parameters. Their research results indicated that the presence of nitrate and boron in the groundwater was related to the use of undeveloped forested land. Jeong [23] analyzed the effect of land use and urbanization on groundwater contamination in the Taejon area in South Korea. They grouped groundwater quality parameters on the basis of land usage. Ouyang et al. [24] applied trilinear analysis to estimate the impact of land use on groundwater quality. They considered four different types of land use and five groundwater quality parameters. Their results demonstrated that nitrate and nitrite (NOx) were related to the presence of septic tanks and that total organic nitrogen (TON) was the dominant species in the groundwater beneath forested lands. Penha [25] investigated the effect of land use on groundwater quality in Southern Portugal by performing a comparison of 12 groundwater quality parameters and five major land uses. Groundwater beneath olive groves had high levels of electrical conductivity, calcium, potassium, sulfate, and phosphate. Dry crop land was correlated with the presence of calcium, magnesium, chloride, electrical conductivity, phosphate, and sulfate. Vineyard land was strongly correlated with high sulfate and phosphate levels.
The aforementioned studies deliver invaluable insights into the relationships between groundwater quality and land use. It should also be noted that the impact of land use on groundwater quality may also be significantly affected by the characteristics of unsaturated soil located between the land surface and the water table of a shallow unconfined aquifer. However, the existing studies on this topic are relatively sparse and even fewer include the factor of the material properties of the unsaturated soil in their investigations. The material properties of unsaturated soil will play an important role in determining aquifer vulnerability. Thus, this study aims to characterize the relationships between groundwater quality and land use in the Pingtung Plain by considering the material properties of the unsaturated soil. The study combines the data regarding groundwater quality parameters, land use patterns, and material properties of unsaturated soil. The factor analysis method is used to analyze the measured concentrations of 14 different water quality parameters for groundwater samples from 46 observation wells, the area percentages for nine different land use categories in the neighborhood of the 46 observation wells, and the thickness of four different types of unsaturated soil based on the core samples obtained during the establishment of the 46 observation wells. Specifically, the scores for the factors are mapped to demarcate the areas that are more vulnerable to land use contamination, prioritize the areas where more intensive monitoring of groundwater quality might be needed, evaluate current land use practices, and adopt new measures to better prevent or control potential sources of pollution.

Study Area
The Pingtung Plain lies in the southeast part of Taiwan and covers most of Pingtung County and a small part of Southeastern Kaohsiung City, as depicted in Figure 1. The plain is bound by the Taiwan Strait on the south, the foothills of the Central Mountain Range and river valleys on the north, the Fengshan Fault on the west, and the Chaozhou Fault on the east. It has a total area of 1270 km 2 , is approximately 60 km long from north to south and 20 km wide from west to east, and is divided into 30 townships. The total population is more than 870,000 persons distributed nonuniformly throughout the area. The Kaoping River is the largest river crossing the plain, flowing from the Central Mountain Range to the Taiwan Strait, with other shorter rivers such as the Tungkang River, Linbian River, and Shihwen River also passing through the plain. According to climate statistics for the period between 2010 and 2018, the Pingtung Plain receives an average annual precipitation ranging from 1160 to 3675 mm, with an average of 2428 mm. Most of the precipitation is concentrated between May and September.
Sustainability 2020, 12, x FOR PEER REVIEW 4 of 24 unsaturated soil based on the core samples obtained during the establishment of the 46 observation wells. Specifically, the scores for the factors are mapped to demarcate the areas that are more vulnerable to land use contamination, prioritize the areas where more intensive monitoring of groundwater quality might be needed, evaluate current land use practices, and adopt new measures to better prevent or control potential sources of pollution.

Study Area
The Pingtung Plain lies in the southeast part of Taiwan and covers most of Pingtung County and a small part of Southeastern Kaohsiung City, as depicted in Figure 1. The plain is bound by the Taiwan Strait on the south, the foothills of the Central Mountain Range and river valleys on the north, the Fengshan Fault on the west, and the Chaozhou Fault on the east. It has a total area of 1,270 km 2 , is approximately 60 km long from north to south and 20 km wide from west to east, and is divided into 30 townships. The total population is more than 870,000 persons distributed nonuniformly throughout the area. The Kaoping River is the largest river crossing the plain, flowing from the Central Mountain Range to the Taiwan Strait, with other shorter rivers such as the Tungkang River, Linbian River, and Shihwen River also passing through the plain. According to climate statistics for the period between 2010 and 2018, the Pingtung Plain receives an average annual precipitation ranging from 1,160 to 3,675 mm, with an average of 2,428 mm. Most of the precipitation is concentrated between May and September. different aquifers. The results of these drilling and sequence stratigraphic studies demonstrated that the geology underlying the plain is comprised of unconsolidated Late Pleistocene and Holocene age sediments and contains abundant groundwater. Most of the sediments consist of coastal and estuarine sand and mud, with abundant shallow marine and lagoon shells and foraminifers. The plain can be partitioned into proximal-fan and distal-fan areas. The deposits in the distal fan area can be further classified into eight overlapping sequences, including four marine sequences and four non-marine sequences. The non-marine sequences are composed of highly permeable coarse sediments and are thus classified as "aquifers", while the marine sequences which contain less permeable fine sediments are regarded as "aquitards". It should be noted that the aquitards are found mainly in the distal-fan area rather than in the proximal-fan area. Figure 2 shows the hydrogeological profile of the study area from west to east. There are four usable aquifers, labeled Aquifer 1, Aquifer 2, Aquifer 3, and Aquifer 4, from top to bottom, at depths of 0-70, 40-130, 90-180, and 160-250 m, respectively. The principal source of groundwater in the plain is from the infiltration of natural rainwater, which collects in the principal, ancient Quaternary reservoir. The proximal-fan area and the river valleys on the eastern and northern boundaries are the major regions for aquifer recharging. Groundwater flows from these regions to the western coastal area bordering the Taiwan Strait. The land use patterns depicted in Figure 3 are based on information obtained from the Land Use Investigation of Taiwan by the National Land Survey and Mapping Center of the Ministry of Interior, Taiwan. Approximately 50.5% of the area is used for agriculture and 5.5% for fishponds. During dry months or years, large amounts of groundwater are extracted to meet the water resource requirements for farmlands, fishponds, and households. This has led to an increase in the salinity of the groundwater, a reduction in the pollution diluting capability of the surface water, and an increase in the occurrence of severe land subsidence and seawater intrusion.
Ministry of Interior, Taiwan. Approximately 50.5% of the area is used for agriculture and 5.5% for fishponds. During dry months or years, large amounts of groundwater are extracted to meet the water resource requirements for farmlands, fishponds, and households. This has led to an increase in the salinity of the groundwater, a reduction in the pollution diluting capability of the surface water, and an increase in the occurrence of severe land subsidence and seawater intrusion.

Data Collection and Processing
Three types of data including groundwater quality parameters, land use categories, and material types of unsaturated soil are required to characterize the relationships between groundwater quality and land use. The data matrices of concentrations for 14 groundwater quality parameters, area percentages for 9 different land use categories, and thickness for 4 different material types of unsaturated soil are provided in Tables A1-A3 in Appendix A. After the "Integrated Groundwater Monitoring Network Project", long-term surveys of the groundwater quality in established groundwater observation wells have continued to be conducted by the Agricultural Engineering Research Center (AERC), with financial support from the Taiwan Water Resource Agency. The measurement data for the groundwater quality parameters considered in this study are mainly obtained from the annual reports published by the Water Resources Agency from 2014 to 2019. Only groundwater quality for Aquifer 1 (a shallow unconfined aquifer) is considered. We look at 14 groundwater quality parameters including electrical conductivity (EC), Escherichia coli (E. Coli), Na + , K + , Mg 2+ , Ca 2+ , Cl − , SO4 2− ,

Data Collection and Processing
Three types of data including groundwater quality parameters, land use categories, and material types of unsaturated soil are required to characterize the relationships between groundwater quality and land use. The data matrices of concentrations for 14 groundwater quality parameters, area percentages for 9 different land use categories, and thickness for 4 different material types of unsaturated soil are provided in Tables A1-A3 in Appendix A. After the "Integrated Groundwater Monitoring Network Project", long-term surveys of the groundwater quality in established groundwater observation wells have continued to be conducted by the Agricultural Engineering Research Center (AERC), with financial support from the Taiwan Water Resource Agency. The measurement data for the groundwater quality parameters considered in this study are mainly obtained from the annual reports published by the Water Resources Agency from 2014 to 2019. Only groundwater quality for Aquifer 1 (a shallow unconfined aquifer) is considered. We look at 14 groundwater quality parameters including electrical conductivity (EC), Escherichia coli (E. Coli), Na As, Fe 3+ , and Mn 2+ from 46 observation wells in the Pingtung Plain established by the Taiwan Water Resource Agency. Data for land use are obtained from the Land Use Investigation of Taiwan by the National Land Survey and Mapping Center of the Ministry of Interior, Taiwan. Because the data for land use patterns are categorical rather than numerical, they must be converted to facilitate the execution of factor analysis. The area percentages for different land use categories in the vicinity of the well have been used in several studies to facilitate the characterization of the relationship between land use and groundwater quality [27]. However, selection of the size of the area surrounding a well for calculating the area percentages of different land use categories is subjective. Wang [21] argued that a circular area around a well with a radius of 1000 ft was optimal to evaluate the relationship between groundwater quality and land use. Herein, we consider 9 land use categories including rice field, dryland, fruit farm, uncultivated land, aquaculture, livestock, forest, water conservation, and human settlement. Table A2 summarizes the area percentages for the 9 land use categories in the vicinity of the 46 observation wells. The type of geological material in unsaturated soil is often the main factor controlling the vertical movement of pollutants from the land surface as they seep into shallow unconfined aquifers. The thickness of the different types of geological material can be used to calculate the permeability of such unsaturated soil. The thicknesses of the various types of unsaturated soil are calculated from borehole data from observation wells that have been published by the Taiwan Central Geological Survey [28]. The soil types are classified into four groups including gravel, coarse sand, fine sand, and clay, corresponding to grain sizes of >2.0, 0.25-2.0, 0.063-0.25, and <0.063 mm, respectively.

Factor Analysis
Factor analysis (FA) is a multivariate statistical approach that is extensively used to describe the general relationships between several observed variables in terms of a potentially lower number of unobserved variables, which are called factors, with minimum loss of the original information. The values of the different variables should be standardized and normalized prior to FA to avoid the problem of "no commensurate units". The mean and standard deviation of the standardized variables are thus zero and unity, respectively. All variables are standardized by applying the following transformation: where Z i,j are the jth values of the ith standardized variables; X i,j are the jth observations of the ith variables, X i is the mean of the ith variable; S i is the standard deviation of the ith variable.
In this study, factor analysis is performed by using the principal component analysis (PCA) method. The procedure for FA is as follows. The correlation matrix, i.e., the array of correlation coefficients for all possible pairs of all standardized variables, is calculated. After the calculation of the correlation matrix, the appropriateness of the factor model is evaluated. The Kaiser-Meyer-Olkin (KMO) index is computed to measure the sampling adequacy that indicates the proportion of variance that is common variance, i.e., that which might be caused by underlying factors. A higher value generally indicates that FA analysis may be suitable. When the KMO value is greater than 0.5 and the significance (p value) is smaller than 0.001, it indicates that the data are suitable for FA. Subsequently, the eigenvalues and their corresponding eigenvectors (principal components) as well as the variance for the correlation matrix are obtained. The sequence of the factors corresponding to the magnitudes of the eigenvalues is defined. The values of the eigenvalues and variances associated with each individual factor are summed up to be expressed as a cumulative eigenvalue and percentage of variance, respectively. Although the factor matrix obtained in the extracted phase is indicative of the relationship between the factors and the individual variables, it is usually difficult to identify meaningful factors. Therefore, a process for rotation of the factor axis is commonly executed to yield a structure where the factors are clearly marked by high loadings for some variables and low loadings for others, thus facilitating the identification of meaningful factors. In this study, Kaiser's varimax rotation scheme is employed [29].

Results and Discussion
This study applies FA to characterize the relationship between the groundwater quality and land use, combining the measured concentrations for 14 groundwater quality parameters sampled from 46 observation wells belonging to the Taiwan WRA, the area percentages for nine land use categories in the vicinity of these 46 observation wells, and the thicknesses of four types of unsaturated soil according to the core samples obtained during the establishment of the 46 observation wells. Prior to investigation of the relationship between groundwater quality and land use by FA, a descriptive statistical analysis of all the collected data is routinely executed. Table 1 shows the descriptive statistics for the 14 groundwater quality parameters, area percentages for the nine land use patterns, and the thicknesses of the four types of unsaturated soil. As can be seen in Table 1, there is moderate to low variability (standard deviation) for the majority of the variables. The highest variabilities in the measured concentrations are obtained for the following groundwater quality parameters: E. Coli, Fe 3+ , Cl − , Na + , EC, K + , Mg 2+ , As, SO 4 2− , Mn 2+ , NO 3 − -N, and TOC, with the coefficient of variation being above 1.0, followed by Ca 2+ and HCO 3 − . Land use for water conservation and uncultivated land show the greatest variability of area percentages for land use patterns. The largest variability shown for thickness of unsaturated soil type is for coarse sand. Before performing factor analysis, KMO and Bartlett's test of sphericity are conducted to examine whether the data are suitable for FA. The computing software package IBM SPSS Statistics Version 22 (IBM Corp., Armonk, NY) is used to perform FA. A KMO value of 0.618 and p value of much less than 0.001 calculated in this study indicate that the data are suitable for FA. Then, the correlation matrix for all standardized variables is calculated, as shown in Table A4. Table 2 lists the eigenvalues, percentage of variances, cumulative eigenvalues, and cumulative percentage of variances associated with the first six factors. It can be seen that the first six factors account for 67% of the total variance of 27 variables.   Table 3 summarizes the individual factor loadings from the varimax rotation factor matrix for a four-factor model. Following Liu et al. [30], the terms "weak", "moderate", and "strong" are defined as corresponding to absolute loading values of <0.4, 0.4~0.6, and >0.6, respectively.   [30] applied FA to assess the groundwater quality in areas of Taiwan with Blackfoot disease. They identified this type of factor as duo to "the seawater salinization" because of the presence of EC, Cl − , SO 4 2− , Ca 2+ , Mg 2+ , Na + , and K + , which are indicative of the effects of seawater intrusion on groundwater quality. Figure 4 shows a geographic visualization of the scores for factor 1. High scores are found for observation wells in the townships of Sinyuan, Donggang, Linbian, and Fangliao, which are situated in coastal regions. Figure 5 depicts a geographic visualization of EC concentration of land used for water conservation.   Figure 6 shows a geographic visualization of the scores for factor 2. High negative absolute scores are found for observation wells in townships of Gaoshu, Yanpu, Changzih, Neipu, and Wanluan near the Central Mountain Range. Figure 7 shows a geographic visualization of HCO 3 − and NO 3 − -N concentration, land used for fruit farming, and the thickness of the gravel. It can be seen that land used for growing fruit has high NO 3 − -N and low HCO 3 − concentrations. Jang and Chen [31] reported that the NO 3 − -N concentration was highly correlated with land use for fruit farming. This result is also similar to those documented by Chen and Liu [32], who explained that the thick layer of gravel present in the proximal part of the Choushui River alluvial fan allows NO 3 − -N to quickly move downward to deep aquifers (more than 200 m down). This factor is related to "nitrate pollution". Factor 3 explains 9% of the total variance, with strong positive loading values for TOC and As. This close correlation between TOC and As was also reported by Liu et al. [30]. They explained that TOC is the major factor controlling the liberation of As. They found that As was adsorbed by the TOC near the redox boundary during oxidation and then liberated by the TOC by reductive dissolution/desorption or by ion exchange with seawater. The very low correlation between As and any land use patterns can be considered a direct consequence of the fact that As is a naturally occurring chemical contaminant. Figure 8 shows a geographic visualization of the scores for factor 3. High absolute score values are found at observation wells in the townships of Donggang, Linbian, Jiadong, and Nanzhou, which are located in the southwestern coastal regions. Figure 9 depicts the spatial distribution of As and TOC. Following Liu et al. [30], factor 3 is related to "arsenic pollution".    Factor 4 explains 7% of the total variance, with strong positive absolute loading values for Mn 2+ and Fe 3+ . The positive correlation between Mn 2+ and Fe 3+ can be explained by the fact that they are both derived from the reductive release of a mineral which typically exists in groundwater. Figure 10 shows a geographic visualization of the scores for factor 4. High positive absolute score values are found at observation wells in the townships of Daliao, Linyuan, which are located in the southwestern coastal regions. Figure 11 depicts a geographic visualization of Mn 2+ and Fe 3+ . This factor is related to reductive dissolution of Fe 3+ and Mn 2+ , because the presence of Mn 2+ and Fe 3+ in the groundwater is a direct result of the reductive dissolution.

Conclusions
In this study, factor analysis is applied to identify the relationship between groundwater quality and land use. The results show that a four-factor model can explain 56% of the total variance. Factor 1 (seawater salinization), which includes the groundwater quality parameters for EC, Cl − , SO 4 2− , Ca 2+ , Mg 2+ , Na + , and K + is moderately correlated to land used for water conservation. Factor 2 (nitrate pollution), which includes the groundwater quality parameters for NO 3 − -N and HCO 3 − , shows a close correlation to land used for fruit farming and the thickness of the gravel in the unsaturated soil. Factor 3 (arsenic pollution), which includes the groundwater quality parameters for TOC and As, is very weakly affected by land use patterns. Factor 4 (reductive dissolution of Mn 2+ and Fe 3+ ), which is related to the concentration of Mn 2+ and Fe 3+ , is also weakly impacted by land use patterns. For a sound plan for safe and sustainable groundwater quality management, identification of the source of contamination is the first priority. The

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A
Tables A1-A3 give the data matrices of concentrations for 14 groundwater quality parameters, area percentages for 9 different land use categories and thickness for 4 different material types of unsaturated soil. Table A4 give the matrix of correlation coefficients for all variables.