Next Article in Journal
Digital Twins: A Critical Discussion on Their Potential for Supporting Policy-Making and Planning in Urban Logistics
Next Article in Special Issue
Reliability of Relative Permeability Measurements for Heterogeneous Rocks Using Horizontal Core Flood Experiments
Previous Article in Journal
Assessing Walking and Cycling around Schools
Previous Article in Special Issue
Practical Application of the Galerkin Finite Element Method with a Mass Conservation Scheme under Dirichlet Boundary Conditions to Solve Groundwater Problems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Department of Nursing, Fooyin University, Kaohsiung City 831, Taiwan
2
Graduate Institute of Applied Geology, National Central University, Taoyuan City 320, Taiwan
3
Department of Water Resource and Environmental Engineering, Tamkang University, New Taipei City 251, Taiwan
4
Agricultural Engineering Research Center, Taoyuan City 320, Taiwan
*
Author to whom correspondence should be addressed.
Sustainability 2020, 12(24), 10608; https://0-doi-org.brum.beds.ac.uk/10.3390/su122410608
Submission received: 16 October 2020 / Revised: 20 November 2020 / Accepted: 23 November 2020 / Published: 18 December 2020

Abstract

:
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 specified 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 unconfined groundwater system. Thus, the quality of the water stored in shallow aquifers has been significantly affected 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, SO42−, Cl, Ca2+, Mg2+, Na+, and K+, shows a moderate correlation to land used for water conservation. Factor 2 (nitrate pollution), which includes the groundwater quality parameters of NO3-N and HCO3, 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 Fe3+ and Mn2+), which involves Mn2+ and Fe3+, 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.

1. 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.

2. Materials and Methods

2.1. 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 km2, 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.
Drilling and sequence stratigraphic studies for characterizing the subsurface geology and hydrogeology within a subsurface depth of 250 m were conducted from 1995 to 1998 during the “Integrated Taiwan Groundwater Monitoring Network Project”. Fifty hydrogeological investigation stations and 126 monitoring wells were established by the Water Resource Agency (WRA) [26]. The hydrogeological investigation stations each included one or more wells with screening at different depths for observation of groundwater levels and the monitoring of groundwater quality in the 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.

2.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 Table A1, Table A2 and Table 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+, Mg2+, Ca2+, Cl, SO42−, NO3-N, HCO3, TOC, As, Fe3+, and Mn2+ 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.

2.3. 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:
Z i , j = X i , j X i ¯ S i
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].

3. 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, Fe3+, Cl, Na+, EC, K+, Mg2+, As, SO42−, Mn2+, NO3-N, and TOC, with the coefficient of variation being above 1.0, followed by Ca2+ and HCO3. 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.
Factor 1 explains 28% of the total variance and shows strong positive absolute loading values for EC, Cl, SO42−, Ca2+, Mg2+, Na+, K+ and a moderate positive absolute loading value for land used for water conservation. Liu et al. [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, SO42−, Ca2+, Mg2+, 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.
Factor 2 explains 12% of the total variance, with a strong positive absolute loading value for HCO3, a strong negative loading values for NO3-N, strong negative absolute loading values for land used for fruit farming, moderate positive absolute loading values for human settlement, and moderate negative loading values for gravel. 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 HCO3 and NO3-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 NO3-N and low HCO3 concentrations. Jang and Chen [31] reported that the NO3-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 NO3-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 Mn2+ and Fe3+. The positive correlation between Mn2+ and Fe3+ 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 Mn2+ and Fe3+. This factor is related to reductive dissolution of Fe3+ and Mn2+, because the presence of Mn2+ and Fe3+ in the groundwater is a direct result of the reductive dissolution.
Factors 5 and 6 show very low absolute loading values for all groundwater quality parameters and so are not discussed here.

4. 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, SO42−, Ca2+, Mg2+, Na+, and K+ is moderately correlated to land used for water conservation. Factor 2 (nitrate pollution), which includes the groundwater quality parameters for NO3-N and HCO3, 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 Mn2+ and Fe3+), which is related to the concentration of Mn2+ and Fe3+, 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 NO3-N concentration closely correlated with land use for fruit farming and the thickness of the gravel in the unsaturated soil clearly identifies the that fruit farming is the major NO3-N source, and the thick layer of gravel present in the townships of Gaoshu, Yanpu, Changzih, Neipu, and Wanluan near the Central Mountain Range allows NO3-N to quickly move downward to deep aquifers. Following the identification of the NO3-N source, pollution control for these NO3-N contaminated regions could include public education to raise awareness of the use of fertilizer for fruit farming or change the land use for fruit farming.

Author Contributions

C.-P.L. and J.-S.C. conceived the study and its design, C.-H.W. performed the data analysis, T.-W.C. provided the data resources, C.-P.L., S.-W.W., T.-W.C. and J.-S.C. interpreted the results, C.-P.L. wrote the manuscript, S.-W.W. and J.-S.C. made a critical revision. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by Ministry of Science and Technology, Taiwan grant number MOST 107-2410-H-242-002 and The APC was funded Ministry of Science and Technology, Taiwan grant number MOST 109-2410-H-242-003.

Acknowledgments

The authors are grateful to the Ministry of Science and Technology, Taiwan for financial support of this research under contract MOST 107-2410-H-242-002.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1, Table A2 and Table 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.
Table A1. Data matrix for concentrations for 14 water quality parameters.
Table A1. Data matrix for concentrations for 14 water quality parameters.
WellsE.C.
(μ S/cm)
HCO3
(mg/L)
E. Coli
(CFU/100mL)
Cl
(mg/L)
SO42-
(mg/L)
NO3-N
(mg/L)
TOC
(mg/L)
As
(mg/L)
Fe3+
(mg/L)
Mn2+
(mg/L)
Ca2+
(mg/L)
Mg2+
(mg/L)
Na+
(mg/L)
K+
(mg/L)
#1993314.0450.51590.0480.550.0084.3050.600151.00022.97534.2754.503
#2542139.07027.759.50.0050.500.0040.4741.16045.60022.70035.5503.120
#3752274.8100109.51.0330.1220.720.017114.3352.61978.43337.13373.46715.908
#4700213.56317.569.41.9050.600.0010.0250.373104.25022.80021.15015.650
#510,385217.0316.658.050.0050.350.0103.0800.408115.00013.60018.4507.465
#635,541262.41416,562.51971.250.0260.300.0080.9771.311304.0001002.5007947.500318.625
#7498165.023818.31071.3100.350.0010.1020.106105.00015.65012.9001.800
#84700195.519819.040.10.0050.950.0096.1400.45076.80022.05022.4004.350
#9780307.7931.6102.8430.6040.670.0684.9290.246122.57124.25775.2432.490
#1041181.033.928.751.0050.500.0013.6850.29134.5508.34019.6001.385
#11634215.027584.959.352.8150.500.0010.4030.826118.00021.50039.40019.800
#12453133.520579.610.950.0050.800.0760.5650.13748.25018.20046.7002.930
#13515175.09,34517.072.9333.1900.400.0040.9030.16782.06714.86714.2331.513
#14405154.05018.951.555.7400.250.0000.1810.19459.75020.30023.6004.660
#153036184.57528.010.350.0050.350.0100.8820.30645.8008.01565.9001.970
#1643619.3285.335.252.6150.150.0030.4310.0299.2957.35012.8501.091
#17629171.01372.136.250.1050.400.0431.5300.10784.25016.70046.7001.535
#18869340.9231178.4108.40.2560.610.0634.8210.122111.00030.80045.2226.772
#19494177.060025.0103.80.4450.500.0243.7500.51786.15012.58519.3502.425
#20592171.5884.5151.252.6880.200.0000.0180.00295.20016.75014.0001.465
#2153454.070010.312.558.0950.150.0010.0620.00128.7008.43512.8502.030
#22475118.7223.564.51.0330.170.0000.0180.00155.03313.5336.1200.908
#2310,14576.03854.718.32.5350.200.0000.0130.00130.40011.4007.7750.764
#24816334.5526.71570.0600.580.0504.0500.094151.25018.60034.2502.518
#2516,72897.5314,450.016900.0600.250.0173.3001.545289.500833.5006435.000230.500
#26571169.55,50050.561.50.0050.700.0013.0250.60489.00017.25029.7505.305
#27515145.7332.4100.20.6500.130.0000.1960.00463.46717.5336.5270.944
#2832,787202.4816,402.52104.4250.0360.370.0240.3510.077260.325891.5006982.500279.900
#29660158.57538.81180.0250.500.0094.1350.550117.00019.85033.3504.400
#30571107.02505.5102.31.0900.300.0050.1100.00180.15013.90012.3501.920
#31436201.786.4324.7000.400.0000.0630.00466.80023.66711.0331.423
#32400261.129307.96.5330.0830.230.02211.5300.58664.7829.18612.7972.154
#33500124.3603.433.1331.7770.200.0000.0680.00245.20015.3006.6130.960
#3428138.2906.214.156.6000.080.0000.1070.00119.5506.9607.7151.373
#35606195.2833.9138.82.2820.360.0000.1350.00394.16017.94010.4481.630
#36591101.084.860.354.6000.150.0000.0410.00164.10014.9008.2001.130
#3710,333114.5104.5121.56.3350.250.0000.0320.00189.55016.8008.3601.850
#38724237.658.8120.40.7900.300.0071.2190.284141.20025.70027.2003.140
#3966980.06503.432.853.2300.200.0022.0452.56738.8508.11514.5001.013
#40592290.012026.32.20.0303.350.2910.6460.02920.20027.90078.25045.000
#41637102.0284.432.17.1350.200.0000.1430.00952.85015.7509.1201.225
#42587211.2904.2127.861.2920.180.0000.0340.00292.22018.46014.4801.980
#43502118.5507.293.850.8600.150.0011.9400.10074.35013.15014.6501.715
#44251113.03551.919.40.1850.500.0112.2300.49820.0507.16043.8501.169
#45369127.71003.91.60.0050.350.0010.3250.23942.1504.85530.3501.570
#46520169.51758.447.059.3800.200.0000.4810.013103.70018.3008.5701.840
Table A2. Data matrix for area percentages for 9 different land use categories.
Table A2. Data matrix for area percentages for 9 different land use categories.
WellsRice Field,%Dryland,
%
Fruit Farm,%Uncultivated Land,%Aquaculture,
%
Livestock,
%
Forest,
%
Water Conservation,
%
Human
Settlement,
%
#117122203100214
#216100025044
#3062000212027
#400325000445
#50222000124
#602205200228
#71466064068
#8517006300113
#9203200200183
#1005640000136
#112371101916088
#1223111101915088
#13154380940027
#1400640005325
#150140001568
#160189101005
#170050700711
#1835290100033
#194112110340125
#200123501501695
#21021000731
#220098000002
#2300000018073
#2400000003159
#250000000770
#260680930062
#2701859000004
#28004020086
#2983494220040
#300103002044364
#3106000001037
#3200550003040
#330099000000
#340377033001
#35521480000224
#36392129912323
#371374407107
#38134191254108
#3901641050523
#4061490200035
#411384031004
#423417231120032
#430963000179
#440091201104
#450441400012014
#460155203150010
Table A3. Data matrix for thickness for 4 different material types of unsaturated soil.
Table A3. Data matrix for thickness for 4 different material types of unsaturated soil.
WellsTypes of Unsaturated Soil
Depth (m)Gravel,%Coarse Sand,%Fine Sand,%Clay,%
#139.24940
#263.000315
#339.20314
#427.3183023
#524.018010
#624.50011
#736.20737
#830.90066
#955.60030
#1036.50045
#1129.40122
#1247.50200
#1333.80515
#1442.151630
#1553.00736
#1647.325440
#1735.83300
#1844.00302
#1938.20550
#20114.2173401
#2135.522070
#2271.544000
#2366.9130815
#2449.40136
#2525.33000
#2624.000016
#2778.3591231
#2826.00030
#2927.03500
#3071.9023115
#3162.031000
#32116.4130140
#3370.545000
#34108.963004
#3575.865000
#36112.030411
#37100.044000
#38111.6201041
#39149.531028
#4066.50003
#4195.894016
#42126.145702
#43117.8951213
#4451.828511
#45242.114400
#4677.737000
Table A4. Matrix of correlation coefficients for all variables.
Table A4. Matrix of correlation coefficients for all variables.
ECHCO3E. ColiClSO42NO3-NTOCAsFe3+Mn2+Ca2+Mg2+Na+K+Rice FieldDry-LandFruit FarmUncultivated
Land
Aqua-CultureLive-StockForestWater ConservationHuman
Settlement
GravelCoarse SandFine SandClay
EC1.000.09−0.100.910.91−0.16−0.08−0.03−0.050.180.750.910.920.92−0.16−0.22−0.24−0.040.31−0.15−0.060.26−0.13−0.15−0.22−0.01−0.13
HCO3 1.000.020.090.11−0.430.420.410.260.110.430.090.080.130.250.24−0.39−0.160.21−0.03−0.170.030.11−0.30−0.07−0.05−0.17
E. Coli 1.00−0.08−0.090.010.01−0.06−0.010.02−0.06−0.09−0.08−0.090.12−0.140.01−0.110.020.06−0.08−0.120.22−0.17−0.040.050.15
Cl 1.00 0.99−0.20−0.070.00−0.040.260.830.990.990.99−0.13−0.17−0.25−0.120.24−0.15−0.130.43−0.19−0.20−0.17−0.08−0.17
SO42 1.00−0.21−0.09−0.02−0.060.230.860.990.990.98−0.10−0.18−0.24−0.090.24−0.15−0.140.42−0.19−0.18−0.13−0.08−0.17
NO3-N 1.00−0.30−0.26−0.16−0.25−0.27−0.20−0.20−0.21−0.15−0.120.380.10−0.190.290.04−0.23−0.180.440.01−0.13−0.04
TOC 1.000.920.100.02−0.10−0.06−0.070.020.140.14−0.29−0.070.28−0.03−0.10−0.04−0.02−0.33−0.17−0.130.03
As 1.000.01−0.08−0.060.000.000.070.140.03−0.22−0.140.15−0.05−0.140.10−0.11−0.27−0.16−0.11−0.12
Fe3+ 1.000.590.00−0.03−0.04−0.02−0.060.44−0.17−0.07−0.07−0.040.21−0.060.07−0.16−0.050.00−0.02
Mn2+ 1.000.250.280.270.27−0.030.13−0.19−0.040.060.050.140.160.08−0.38−0.060.010.09
Ca2+ 1.000.840.830.820.05−0.09−0.36−0.040.280.01−0.230.49−0.13−0.23−0.08−0.15−0.21
Mg2+ 1.000.990.99−0.12−0.17−0.25−0.110.27−0.15−0.120.44−0.19−0.20−0.17−0.08−0.17
Na+ 1.000.99−0.12−0.18−0.24−0.110.27−0.15−0.130.43−0.19−0.20−0.17−0.07−0.16
K+ 1.00−0.12−0.17−0.26−0.110.30−0.14−0.130.40−0.20−0.22−0.18−0.09−0.15
Rice field 1.000.00−0.19−0.080.140.37−0.20−0.08−0.14−0.11−0.02−0.02−0.22
Dryland 1.00−0.35−0.07−0.01−0.050.09−0.240.12−0.01−0.01−0.18−0.20
Fruit farm 1.000.15−0.260.04−0.22−0.27−0.390.510.16−0.09−0.02
Uncultivated
land
1.00−0.150.04−0.09−0.150.200.13−0.01−0.010.28
Aquaculture 1.000.12−0.110.08−0.23−0.220.07−0.03−0.02
Livestock 1.00−0.11−0.09−0.17−0.030.00−0.16−0.14
Forest 1.00−0.120.22−0.070.220.110.20
Water conservation 1.00−0.21−0.260.05−0.18−0.08
Human Settlement 1.00−0.26−0.130.390.50
Gravel 1.00−0.07−0.21−0.26
Coarse sand 1.00−0.060.06
Fine sand 1.000.14
Clay 1.00

References

  1. Consuming Water Resources Database of Taiwan. Available online: http://wuss.wra.gov.tw (accessed on 30 June 2019).
  2. Agriculture Engineering Research Center. Survey, Analysis and Assessment of Groundwater Quality in Taiwan Areas in 2009; Water Resources Agency, Ministry of Economic Affairs, Executive Yuan: Taiwan, 2009; pp. 100–102. [Google Scholar]
  3. Agriculture Engineering Research Center. Survey, Analysis and Assessment of Groundwater Quality of Taiwan in 2010; Water Resources Agency, Ministry of Economic Affairs, Executive Yuan: Taiwan, 2010; pp. 151–158. [Google Scholar]
  4. Agriculture Engineering Research Center. Survey, Analysis and Assessment of Groundwater Quality of Taiwan in 2011; Water Resources Agency, Ministry of Economic Affairs, Executive Yuan: Taiwan, 2011; pp. 191–203. [Google Scholar]
  5. Agriculture Engineering Research Center. Survey, Analysis and Assessment of Groundwater Quality of Taiwan in 2012–2014 (1/3); Water Resources Agency, Ministry of Economic Affairs, Executive Yuan: Taiwan, 2012; pp. 85–101. [Google Scholar]
  6. Hsu, L.I.; Chen, G.S.; Lee, C.H.; Yang, T.Y.; Chen, Y.H.; Wang, Y.H.; Hsueh, Y.M.; Chiou, H.Y.; Wu, M.M.; Chen, C.J. Use of Arsenic-Induced Palmoplantar Hyperkeratosis and Skin Cancers to Predict Risk of Subsequent Internal Malignancy. Am. J. Epidemiol. 2013, 177, 202–212. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Liang, C.P.; Jang, C.S.; Chen, C.F.; Chen, J.S. Zonal Management of Multi-Purposes Groundwater Utilization Based on Water Quality and Impact on the Aquifer. Environ. Monit. Assess. 2016, 188, 431. [Google Scholar] [CrossRef] [PubMed]
  8. Liang, C.P.; Chien, Y.C.; Jang, C.S.; Chen, C.F.; Chen, J.S. Spatial Analysis of Human Health Risk Due to Arsenic Exposure through Drinking Groundwater in Taiwan’s Pingtung Plain. Int. J. Environ. Res. Public Health 2017, 14, 81. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Liang, C.P.; Chen, J.S.; Chien, Y.C.; Chen, C.F. Spatial Analysis of the Risk to Human Health from Exposure to Arsenic Contaminated Groundwater: A Kriging Approach. Sci. Total Environ. 2018, 627, 1048–1057. [Google Scholar] [CrossRef]
  10. Liang, C.P.; Hsu, W.S.; Chien, Y.C.; Wan, S.W.; Chen, J.S. The Combined Use of Groundwater Quality, Drawdown Index and Land Use to Establish a Multi-Purpose Groundwater Utilization Plan. Water Resour. Manag. 2019, 33, 4231–4247. [Google Scholar] [CrossRef]
  11. Bulka, C.M.; Jones, R.M.; Turyk, M.E.; Stayner, L.T.; Argos, M. Arsenic in Drinking Water and Prostate Cancer in Illinois Counties: An Ecologic Study. Environ. Res. 2016, 148, 450–456. [Google Scholar] [CrossRef] [Green Version]
  12. States, J.C. Disruption of Mitotic Progression by Arsenic. Biol. Trace Elem. Res. 2015, 166, 34–40. [Google Scholar] [CrossRef]
  13. Chen, C.; Jiang, X.; Zhang, Z. MicroRNA-155 Regulates Arsenite-Induced Malignant Transformation by Targeting Nrf2-Mediated Oxidative Damage in Human Bronchial Fepithelial Cells. Toxicol. Lett. 2017, 278, 38–47. [Google Scholar] [CrossRef]
  14. Liang, C.P.; Jang, C.S.; Liu, C.W.; Lin, K.H.; Lin, M.C. An Integrated GIS-Based Approach in Assessing Carcinogenic Risks via Food-Chain Exposure in Arsenic-Affected Groundwater Areas. Environ. Toxicol. 2010, 25, 113–123. [Google Scholar] [CrossRef]
  15. Liang, C.P.; Liu, C.W.; Jang, C.S.; Wang, S.W.; Lee, J.J. Assessing and Managing the Health Risk Due to Ingestion of Inorganic Arsenic from Fish and Shellfish Farmed in Blackfoot Disease Areas for General Taiwanese. J. Hazard. Mater. 2011, 186, 622–628. [Google Scholar] [CrossRef]
  16. Liang, C.P.; Jang, C.S.; Chen, J.S.; Wang, S.W.; Lee, J.J.; Liu, C.W. Probabilistic Health Risk Assessment for Ingestion of Seafood Farmed in Arsenic Contaminated Groundwater in Taiwan. Environ. Geochem. Health 2013, 35, 455–464. [Google Scholar] [CrossRef] [PubMed]
  17. Lerner, D.N.; Harris, R.C. The Relationship between Land Use and Groundwater Resources and Quality. Land Use Policy 2009, 26, 265–273. [Google Scholar] [CrossRef]
  18. Jang, C.S.; Lin, C.W.; Liang, C.P.; Chen, J.S. Developing a Reliable Model for Aquifer Vulnerability. Stoch. Environ. Res. Risk Assess. 2016, 30, 175–187. [Google Scholar] [CrossRef]
  19. Liang, C.P.; Hsu, S.Y.; Chen, J.S. An Analytical Model for Solute Transport in an Infiltration Tracer Test in Soil with a Shallow Groundwater Table. J. Hydrol. 2016, 540, 129–141. [Google Scholar] [CrossRef]
  20. Liu, C.W.; Lin, C.N.; Jang, C.S.; Ling, M.P.; Tsai, J.W. Assessing Nitrate Contamination and Its Potential Health Risk to Kinmen Residents. Environ. Geochem. Health 2011, 33, 503–514. [Google Scholar] [CrossRef]
  21. Wang, X. The Relationship of Land Use and Groundwater Quality: A Case Study of Rhode Island. Master’s Thesis, University of Rhode Island, Kingston, RI, USA, 1989. [Google Scholar]
  22. Eckhardt, D.V.A.; Stackelberg, P.E. Relation of Ground-Water Quality to Land Use on Long Island, New York. Ground Water 1995, 33, 1019–1031. [Google Scholar] [CrossRef]
  23. Jeong, C.H. Effect of Land Use and Urbanization on Hydrochemistry and Contamination of Groundwater from Taejon Area, Korea. J. Hydrol. 2001, 253, 194–210. [Google Scholar] [CrossRef]
  24. Ouyang, Y.; Zang, J.E.; Cui, L. Estimating Impact of Land Use on Groundwater Quality Using Trilinear Analysis. Environ. Monit. Assess. 2014, 186, 5353–5362. [Google Scholar] [CrossRef]
  25. Penha, A.M.; Chambel, A.; Murteira, M.; Morais, M. Influence of Different Land Use on Groundwater Quality in Southern Portugal. Environ. Earth Sci. 2016, 75, 622. [Google Scholar] [CrossRef]
  26. Water Resource Agency, Taiwan. Available online: https://www.wra.gov.tw/ (accessed on 30 June 2019).
  27. Heißerer, T.; Haslauer, C.; Bárdossy, A. Including Land Use Information for the Spatial Estimation of Groundwater Quality Parameters—1. Local Estimation Based on Neighbourhood Composition. J. Hydrol. 2016, 535, 688–698. [Google Scholar] [CrossRef]
  28. Central Geological Survey, Taiwan. Available online: https://www.moeacgs.gov.tw/ (accessed on 30 June 2019).
  29. Kaiser, H.F. The Varimax Criterion for Analytic Rotation in Factor Analysis. Psychometrika 1958, 23, 187–201. [Google Scholar] [CrossRef]
  30. Liu, C.W.; Lin, K.H.; Kuo, Y.M. Application of Factor Analysis in the Assessment of Groundwater Quality in a Blackfoot Disease Area in Taiwan. Sci. Total Environ. 2003, 313, 77–89. [Google Scholar] [CrossRef]
  31. Jang, C.S.; Chen, S.K. Integrating Indicator-Based Geostatistical Estimation and Aquifer Vulnerability of Nitrate-N for Establishing Groundwater Protection Zone. J. Hydrol. 2015, 523, 441–451. [Google Scholar] [CrossRef]
  32. Chen, W.F.; Liu, T.K. Dissolved Oxygen and Nitrate of Groundwater in Choshui Fan-Delta, Western Taiwan. Environ. Geol. 2003, 44, 731–737. [Google Scholar] [CrossRef]
Figure 1. Study area in the Pingtung Plain.
Figure 1. Study area in the Pingtung Plain.
Sustainability 12 10608 g001
Figure 2. Hydrogeological profile of the Pingtung Plain.
Figure 2. Hydrogeological profile of the Pingtung Plain.
Sustainability 12 10608 g002
Figure 3. Map of land use in the Pingtung Plain.
Figure 3. Map of land use in the Pingtung Plain.
Sustainability 12 10608 g003
Figure 4. Geographical visualization of factor scores for factor 1, the red area (high score) represent high seawater salinization in the coast.
Figure 4. Geographical visualization of factor scores for factor 1, the red area (high score) represent high seawater salinization in the coast.
Sustainability 12 10608 g004
Figure 5. Spatial distribution of the EC and water conservation.
Figure 5. Spatial distribution of the EC and water conservation.
Sustainability 12 10608 g005
Figure 6. Geographical visualization of factor score for factor 2, the green and white area (low score) represent high nitrate pollution.
Figure 6. Geographical visualization of factor score for factor 2, the green and white area (low score) represent high nitrate pollution.
Sustainability 12 10608 g006
Figure 7. Spatial distribution of the NO3-N, HCO3-, fruit farm, and thickness of gravel.
Figure 7. Spatial distribution of the NO3-N, HCO3-, fruit farm, and thickness of gravel.
Sustainability 12 10608 g007
Figure 8. Geographical visualization of factor score for factor 3, the red area (high score) represent high arsenic pollution.
Figure 8. Geographical visualization of factor score for factor 3, the red area (high score) represent high arsenic pollution.
Sustainability 12 10608 g008
Figure 9. Spatial distribution of As and TOC.
Figure 9. Spatial distribution of As and TOC.
Sustainability 12 10608 g009
Figure 10. Geographical visualization of factor score for factor 4, the red area (high score) represent high reductive dissolution of Fe3+ and Mn2+ in this area.
Figure 10. Geographical visualization of factor score for factor 4, the red area (high score) represent high reductive dissolution of Fe3+ and Mn2+ in this area.
Sustainability 12 10608 g010
Figure 11. Spatial distribution of Mn2+and Fe3+.
Figure 11. Spatial distribution of Mn2+and Fe3+.
Sustainability 12 10608 g011
Table 1. Descriptive statistics for 14 groundwater quality parameters, percentages of areas for 9 land use patterns and thicknesses of 4 unsaturated soil types from 46 observation wells in the Pingtung Plain.
Table 1. Descriptive statistics for 14 groundwater quality parameters, percentages of areas for 9 land use patterns and thicknesses of 4 unsaturated soil types from 46 observation wells in the Pingtung Plain.
ParameterAverage
(mg/L)
Maximum
(mg/L)
Minimum
(mg/L)
Standard Deviation
(-)
Coefficient of Variation
(-)
Exceed the Standard for DrinkingPercentage,%Exceed the Standard for AgriculturePercentage,%Exceed the Standard for AquaculturePercentage,%
EC3155.6035,541.13251.007486.012.4----1328----
HCO3170.89340.8919.3077.460.59----------
E. Coli503.079345.002.501608.823.21533----00
Cl1075.4516,562.501.853945.153.74949----
SO42185.862104.431.03467.962.53737----
NO3-N1.869.380.012.471.300--------
TOC0.443.350.080.481.100--------
As0.020.290.000.052.612264949
Fe3+4.08114.340.0116.774.1296337----
Mn2+0.372.620.000.601.629632044----
Ca2+88.51304.009.3063.230.7------------
Mg2+74.971002.504.86223.613.0------------
Na+488.157947.506.121778.953.6------------
K+22.13318.630.7668.913.1------------
Rice field4.5240.680.009.472.09
Dryland12.2661.560.0017.211.40
Fruit farm30.2199.330.0031.401.04
Uncultivated land0.515.180.001.212.36
Aquaculture6.7063.390.0012.941.93
Livestock2.3915.660.004.221.77
Forest2.6224.940.005.572.12
Water conservation7.2576.960.0016.302.25
Human Settlement20.4778.960.0021.491.05
Gravel14.4464.740.0019.071.32
Coarse sand4.2234.340.006.601.56
Fine sand2.2613.520.003.051.35
Clay4.0922.540.005.851.43
Table 2. The factor loading of varimax rotation.
Table 2. The factor loading of varimax rotation.
Common FactorEigenvaluesPercentages of Variance,%Cumulative Eigenvalues,%
17.422828
23.361240
32.50949
41.83756
51.55662
61.38567
Table 3. The factor loading of varimax rotation.
Table 3. The factor loading of varimax rotation.
ParameterFactor 1Factor 2Factor 3Factor 4Factor 5Factor 6
Water quality
EC0.9450.022−0.005−0.046−0.0290.007
HCO30.0740.6670.2940.1370.101−0.147
E. Coli−0.095−0.003−0.0200.0990.147−0.041
Cl0.9860.0310.0000.041−0.052−0.058
SO42−0.9830.053−0.0320.004−0.034−0.049
NO3-N−0.1530.669−0.212−0.1220.2230.017
TOC−0.0790.2130.9470.0500.0500.019
As−0.0270.1680.923−0.053−0.041−0.104
Fe3+−0.0760.1610.0150.858−0.076−0.098
Mn2+0.2440.089−0.0200.8330.0640.101
Ca2+0.8380.304−0.1170.0540.146−0.012
Mg2+0.9890.0300.0080.051−0.038−0.053
Na+0.9890.0240.0050.043−0.042−0.053
K+0.9850.0340.0860.050−0.024−0.041
Land use
Rice field−0.1460.3410.017−0.1060.612−0.219
Dryland−0.2200.389−0.0150.3720.013−0.180
Fruit farm−0.2500.742−0.148−0.090−0.135−0.039
Uncultivated land−0.067−0.077−0.123−0.0660.0440.803
Aquaculture0.3120.1460.322−0.1540.467−0.077
Livestock−0.113−0.162−0.0540.0840.8370.033
Forest−0.1260.006−0.1110.171−0.1380.012
Water conservation0.4150.0940.0610.051−0.123−0.050
Human settlement−0.1860.494−0.1670.033−0.2400.495
Geological material
Gravel−0.1760.494−0.289−0.249−0.123−0.177
Coarse sand−0.217−0.002−0.196−0.0850.085−0.029
Fine sand−0.0720.287−0.213−0.191−0.2390.017
Clay−0.140−0.0120.0610.045−0.1440.724
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liang, C.-P.; Wang, C.-H.; Wang, S.-W.; Chang, T.-W.; Chen, J.-S. Application of Factor Analysis for Characterizing the Relationships between Groundwater Quality and Land Use in Taiwan’s Pingtung Plain. Sustainability 2020, 12, 10608. https://0-doi-org.brum.beds.ac.uk/10.3390/su122410608

AMA Style

Liang C-P, Wang C-H, Wang S-W, Chang T-W, Chen J-S. Application of Factor Analysis for Characterizing the Relationships between Groundwater Quality and Land Use in Taiwan’s Pingtung Plain. Sustainability. 2020; 12(24):10608. https://0-doi-org.brum.beds.ac.uk/10.3390/su122410608

Chicago/Turabian Style

Liang, Ching-Ping, Chia-Hui Wang, Sheng-Wei Wang, Ta-Wei Chang, and Jui-Sheng Chen. 2020. "Application of Factor Analysis for Characterizing the Relationships between Groundwater Quality and Land Use in Taiwan’s Pingtung Plain" Sustainability 12, no. 24: 10608. https://0-doi-org.brum.beds.ac.uk/10.3390/su122410608

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop