Health Risk Assessment of Household Drinking Water in a District in the UAE

: The quality of household drinking water in a community of 30 houses in a district in Abu Dhabi, United Arab Emirates (UAE) was assessed over a period of one year (January to November 2015). Standard analytical techniques were used to screen for water quality parameters and contaminants of concern. Water quality was evaluated in the 30 households at four sampling points: kitchen faucet, bathroom faucet, household water tank, and main water pipe. The sampling points were chosen to help identify the source when an elevated level of a particular contaminant is observed. Water quality data was interpreted by utilizing two main techniques: spatial variation analysis and multivariate statistical techniques. Initial analysis showed that many households had As, Cd, and Pb concentrations that were higher than the maximum allowable level set by UAE drinking water standards. In addition, the water main samples had the highest concentration of the heavy metals compared to other sampling points. Health risk assessment results indicated that approximately 30%, 55%, and 15% of the houses studied had a high, moderate, and low risk from the prolonged exposure to heavy metals, respectively. The analysis can help with planning a spatially focused sampling plan to conﬁrm the study ﬁndings and set an appropriate course of action. in Baniyas. The C3 cluster characterizes the western and the middle part of Baniyas, and shows a similarity in 98% of the water quality measurements. Phase-2 results (Figure 4) show that the ﬁrst cluster C1 includes around one-third of the studied houses (houses No. 1, 3, 6, 8, 17, 20, 21, 22, 23, 24, and 29). Their water quality parameters had a similarity of 95%. Also, the second cluster C2 includes 37% of the locations studied in Baniyas (houses No. 2, 7, 10, 11, 12, 16, 18, 19, 25, 28, and 30); however, its similarity was quite higher than C1 (97%). The third cluster C3 had the highest similarity between its water quality variables (98%). C3 includes 26% of the houses in Baniyas. Overall, phase-2 clusters had lower similarity between its clusters (87%).


Introduction
Access to safe and adequate drinking water is at the center of the development goals of any country. Concerns have been echoed on whether current practices in the water industry can successfully provide the intended water quality that is suitable for human consumption [1]. Standards and regulations for ensuring acceptable drinking water quality are issued by governments and special organizations [2]. Many local governments and municipalities maintain a drinking water distribution network by the continuous inspection of water mains, tanks, pipes, and junctions [3]. The United Arab Emirates is no exception; the Regulation and Supervision Bureau (RSB) in the Emirate of Abu Dhabi has issued the fourth edition of its water quality regulations [4]. Water quality in the distribution networks can be characterized by analyzing the representative parameters of the microbiological, physicochemical, and aesthetic aspects of water quality [5]. Regulations and standards help raise confidence in the quality of water supplied at the tap. However, the possibility of contamination incidents throughout the production and distribution processes remains the primary reason for consumer concerns regarding the safety of drinking water. The attribution of these incidents and factors influencing their frequency, severity, and consequences has been the subject of many investigations [6][7][8][9]. Of these incidents, contamination in the distribution network, including household plumbing, proves to be the most challenging to detect and control [10][11][12][13][14].
In addition to setting contaminants' maximum concentration levels, regulations also stipulate the frequency of water quality monitoring at various points in the water supply network to ensure compliance with standards [4]. The frequency of sampling usually depends on the population that is supplied water, and the volume of distributed water. However, sampling done at "supply points to consumers" usually means that the samples are taken before they enter the household. Several studies have reported increased concentrations of contaminants in the household tap due to factors related to plumbing, water temperature, and stagnation [13,15,16].
The United Arab Emirates (UAE) suffers from water scarcity where the availability of renewable water resources is insufficient to meet the country's water demands [17]. Nevertheless, the country has successfully provided safe drinking water in compliance with international standards by tapping into unconventional water resources (mainly desalination) as well as maintaining water networks and periodically monitoring the produced water quality. However, many factors could affect water quality within a drinking water distribution network, including the dissolution of pipe materials, and possible cracks and breaks due to aging and stress. Some districts in the UAE's large cities have drinking water networks that were designed for continuous water supply decades ago. However, high population growth and drinking water subsidies lead to increased water demand and system overload. The country continuously upgrades the water supply networks to maintain acceptable service levels; however, water quality monitoring is needed at the point-of-use to prioritize system improvement efforts.
To assess the drinking water quality supplied to consumers, investigations of water quality at the tap need to be conducted. These investigations typically have one purpose, which is the identification of contaminant levels that are higher than set standards [18]. This approach is inadequate, as it lacks the necessary attribution that is required for planning appropriate contaminant control measures. The attribution and source identification of contaminants requires the careful selection of sampling points and frequency. The resulting data are usually too complex for traditional statistical techniques [19]. Attempting to make sense of complex data, several studies have implemented multivariate statistical analysis techniques to identify similarities and common trends that could help in the better interpretation of sampling results [19][20][21][22][23].
In this study, a range of water quality parameters was monitored at various locations in several households in Baniyas, Abu Dhabi. Baniyas is an old district, where residents have raised some concerns regarding their drinking water quality. Since most of the water quality monitoring conducted by the distribution company is confined to sampling water in the distribution network, it would be of importance to investigate the water quality at the household level. The objectives of this study were to identify contaminants that do not comply with the set standards and explore the use of multivariate statistical techniques to interpret water quality data. Spatial (qualitative) analysis was used in this study to explore the spatial variation of contaminant concentrations at different households. Principal component analysis (PCA) and cluster analysis (CA) were used to group households with significant water quality issues to effectively plan future water quality monitoring. Health risk assessment was used to evaluate the risk related to each contaminant of concern. Finally, the implications of the analysis on planning for future water quality monitoring and enhancement are discussed.

Study Area
The district of Baniyas is located 37 km from the capital of the UAE, Abu Dhabi ( Figure 1). The population of Baniyas, according to a report published by the Urban Planning Council in 2011, is about 70,000 people who are a mixture of all segments of society in Abu Dhabi. The Abu Dhabi Water and Electricity Company (ADWEC) is responsible for providing Baniyas inhabitants with a drinking water supply from the Taweelah desalination plant. The quality of the water entering the distribution system is continuously monitored by ADWEC according to UAE drinking water standards. ADWEC is a wholly owned subsidiary of the Abu Dhabi Water and Electricity Authority (ADWEA). The company operates five desalination plants with approximately 938 MGD (Million Gallons per Day) capacity. In addition, the Abu Dhabi Distribution Company (ADDC) is responsible for distributing water to the consumers.

Sample Collection and Preparation
The sampling campaign aimed to be representative of the drinking water supplied to as large a proportion of the Baniyas population as possible. At the time of designing this campaign, there was no readily available data on the number of households and population in the supplied area. Thus, the district size and number of households was estimated based on satellite image analysis (around 3000 houses), which served as a proxy to the served population. Initially, 90 households were identified, which were mainly older houses that were suspected to have water quality issues. These households represented about 3% of the total houses in Baniyas, which provided good spatial coverage of the district. This coverage was considered acceptable for the purposes of our study, as this percentage is sufficient for identifying the main contaminants of concern. Consequently, 90 houses were contacted for voluntary participation in the study, of which only 30 households agreed to participate in the sampling campaign. The objectives of the research were explained to participating households. In addition, for participants who were interested, the outcomes of the study were shared.
Samples from four points in 30 households were collected on two phases (a total of 720 samples) ( Figure 2). Phase-1 sampling (winter season) was conducted during the months of January through April 2015, while phase-2 (summer season) samples were collected from August to November 2015. Grab samples were taken from the four sampling points. Samples were collected from household water tanks (WT) using a sampling rod and 500-mL sterile glass bottles. In addition, samples were collected from three other locations: the kitchen faucet (KF), a bathroom faucet (BF), and from the water main (WM) before it enters the household. Prior to sampling, sampling containers were disinfected with ethanol, rinsed once with ultra-high-quality water, and then rinsed twice with sampled water. Temperature, pH, and electric conductivity were measured onsite, and the sample

Sample Collection and Preparation
The sampling campaign aimed to be representative of the drinking water supplied to as large a proportion of the Baniyas population as possible. At the time of designing this campaign, there was no readily available data on the number of households and population in the supplied area. Thus, the district size and number of households was estimated based on satellite image analysis (around 3000 houses), which served as a proxy to the served population. Initially, 90 households were identified, which were mainly older houses that were suspected to have water quality issues. These households represented about 3% of the total houses in Baniyas, which provided good spatial coverage of the district. This coverage was considered acceptable for the purposes of our study, as this percentage is sufficient for identifying the main contaminants of concern. Consequently, 90 houses were contacted for voluntary participation in the study, of which only 30 households agreed to participate in the sampling campaign. The objectives of the research were explained to participating households. In addition, for participants who were interested, the outcomes of the study were shared.
Samples from four points in 30 households were collected on two phases (a total of 720 samples) ( Figure 2). Phase-1 sampling (winter season) was conducted during the months of January through April 2015, while phase-2 (summer season) samples were collected from August to November 2015. Grab samples were taken from the four sampling points. Samples were collected from household water tanks (WT) using a sampling rod and 500-mL sterile glass bottles. In addition, samples were collected from three other locations: the kitchen faucet (KF), a bathroom faucet (BF), and from the water main (WM) before it enters the household. Prior to sampling, sampling containers were disinfected with ethanol, rinsed once with ultra-high-quality water, and then rinsed twice with sampled water. Temperature, pH, and electric conductivity were measured onsite, and the sample was collected, stored, and transported to the laboratory according to standard methods for the Examination of Water and Wastewater [24]. was collected, stored, and transported to the laboratory according to standard methods for the Examination of Water and Wastewater [24].

Laboratory and Field Methods
In this study, water quality parameters were analyzed according to the Standard Methods for the Examination of Water and Wastewater by using field probes [25]. A multi-parameter probe (HANNA, HI2020, Hanna Instruments, Woonsocket, RI, USA) was used to measure the physicochemical properties of the samples: pH, temperature (Tem), and electric conductivity (EC). Other parameters were determined in the lab, such as: total dissolved solids (TDS), biochemical oxygen demand (BOD), and turbidity. In addition, the analysis included the determination of the heavy metals' concentration (As, Cd, Pb, Cr, Cu, Zn) and cationic composition of the samples (Mg 2+ , Ca 2+ , Ba 2+ , Na + , and K + ) using inductively coupled phase-mass spectrometry (ICP-AES-715). Moreover, ion exchange chromatography (IEC-1100, LevelOne, Dortmund, Germany) was used to determine the concentration of anions in the samples (Cl − , SO4 2− ).
Standard reference materials and procedural blanks were tested after each set of water samples (15 samples) to ensure that the accuracy of the measurements was maintained. The precision of the method was controlled by collecting triplicated samples randomly and testing the agreement between the readings. The results showed a satisfactory agreement between the two readings for the same sample.

Laboratory and Field Methods
In this study, water quality parameters were analyzed according to the Standard Methods for the Examination of Water and Wastewater by using field probes [25]. A multi-parameter probe (HANNA, HI2020, Hanna Instruments, Woonsocket, RI, USA) was used to measure the physicochemical properties of the samples: pH, temperature (Tem), and electric conductivity (EC). Other parameters were determined in the lab, such as: total dissolved solids (TDS), biochemical oxygen demand (BOD), and turbidity. In addition, the analysis included the determination of the heavy metals' concentration (As, Cd, Pb, Cr, Cu, Zn) and cationic composition of the samples (Mg 2+ , Ca 2+ , Ba 2+ , Na + , and K + ) using inductively coupled phase-mass spectrometry (ICP-AES-715). Moreover, ion exchange chromatography (IEC-1100, LevelOne, Dortmund, Germany) was used to determine the concentration of anions in the samples (Cl − , SO 4 2− ).
Standard reference materials and procedural blanks were tested after each set of water samples (15 samples) to ensure that the accuracy of the measurements was maintained. The precision of the method was controlled by collecting triplicated samples randomly and testing the agreement between the readings. The results showed a satisfactory agreement between the two readings for the same sample.

Data Analysis
Data were analyzed and interpreted using three major categories: qualitatively to study the spatial variation of each measurement, statistically to interpret the significance of each parameter, and finally by using a health risk analysis to assess the risk posed by each contaminant. The following sections will elaborate how these analyses carried out and the degree of their significance.

Spatial Analysis
Heavy metals and cations were analyzed qualitatively by producing maps using ArcGIS (ESRI, Redlands, CA, USA) to explore the spatial variation of each potential drinking water contaminant ( Figure 3). The purpose of this step is to identify areas that show a high level of concentration of a certain drinking water contaminant. In addition, the generated maps were used to detect the possible source of the contamination by comparing the variation of the contaminant concentration at each sampling point.

Multivariate Statistical Analysis
After laboratory tests were carried out, the analyzed data was compiled into one database. Descriptive statistics were used to explore the nature of the data and make initial interpretations about the water quality of the network and identify the contaminants of concern. The mean, median, minimum, maximum, and standard deviation values were calculated for each water quality parameter. Moreover, the concentration of each water quality parameter in the samples were compared with the maximum acceptable values set by UAE's drinking water standards [4] (Table 1).
In addition, multivariate statistical analysis techniques were used to understand and evaluate the variations in the water quality parameters from the water main, through the household water tank, to tap. A detailed multivariate statistical analysis was carried out using Minitab software v18 (Minitab Inc., State College, PA, USA). To reduce dimensionality and avoid biases resulting from different units of measurement, data was standardized (z-scale) for each parameter prior to the statistical analysis [19,26]. Multivariate linear Pearson's correlation coefficients between parameters were determined to indicate statistical dependence (Tables A1 and A2). Principal component analysis (PCA) was conducted to interpret the correlation between water quality parameters by reducing the dispersion matrix to a few interpretable linear combinations of the data [19,27]. Each linear combination will correspond to a principal component. The PCA method has been adopted in many water quality studies [19,28]. The method has proven its efficiency in segregating the water quality parameters and correlating them together. Minitab software v18 was used to conduct this analysis.
Cluster analysis (CA) was used to complement the multivariate analysis by creating groupings of similar water quality characteristics that have the same concentration patterns. Agglomerative hierarchical clustering (AHC) is one of the most common clustering methods. AHC has been used in the literature to classify and examine spatial variation among the sampling stations [22,29]. In this study, AHC has been carried out on a normalized dataset using Ward's method and Euclidean distances as a measure of similarity. This method uses the analysis of variance approach to evaluate the distances between clusters, attempting to minimize the sum of squares of any two clusters that can be formed at each step [30]. In this study, a cut-off level based on similarity between the measured distance was set to 97% or more.

Health Risk Assessment
Health risk assessment is a process of quantifying the occurrence likelihood of any variable that affects the health over a specified timescale. Basically, the risk assessment process involves: hazard identification, exposure assessment, dose response, and risk characterization [31,32]. The calculated magnitude of the risk level is the base of any health risk assessment, and is usually expressed in terms of a carcinogenic or non-carcinogenic health risk [33]. In this study, the degree of toxicity of the heavy metals detected in the water samples was evaluated to determine the level of health risk posed by the contaminants. The health risk assessment was carried out using two main indicators: hazard index (HI) and hazard quotient (HQ) as proposed by the United States (US) Environmental Protection Agency (EPA) [19,32,34]. The non-carcinogenic hazard or HQ is defined as the ratio between the human exposure to each heavy metal, which is also reported as average daily dose (ADD) and the reference dose (RfD). On the other hand, the HI is used to estimate the total non-carcinogenic risk. The HI can be determined by combining the individual HQs. Equations (1)-(3) were used to calculate ADD, HQ, and HI, respectively.
where ADD is the average daily dose (µg/kg/day), C wt is the heavy metal concentration (µg/L), IR is the water intake rate (2.0 L/day and 1.0 L/day for adults and children, respectively), EF is the exposure frequency (365 days/year), ED is the exposure duration (30 years and six years for adults and children, respectively), BW is the average body weight (76 kg and 30 kg for adults and children, respectively) [35], and AT is the average time (10,950 days/year and 2190 days/year for adults and children, respectively). RfD is the oral reference dose (µg/kg/day). In order to assess the risk associated with each heavy metal and their accumulative hazard in each house, this study relied on scales of chronic risk assessment set by previous studies [26,28], which commonly include the following categories:

Spatial Variation Analysis
The spatial variation of heavy metals and cations concentrations from the sampled locations was analyzed using ArcGIS interpolation tool and represented in Figures 3 and A1-A3. Based on the averaged concentrations, each individual parameter was evaluated and categorized into four groups based on meeting UAE standards. The most hazardous contaminants were presented in Figure 3, which includes As, Pb, and Cd. Other figures were developed for other parameters, but they did not show any risks due to their relatively low concentrations, which were below the limits set by UAE standards. In addition, the variation of each parameter concentration for the four sampling points was explored. Based on the results obtained, the WM samples had a high concentration of As, Pb and, Cd reaching up to 45 µg/L, 40 µg/L, and 7 µg/L, respectively ( Figure 3). Overall, the concentrations for other sampling points decreased as the water moves from the WM to WT and from the WT to the BF and KF. The sampling point KF had a relatively low concentration of heavy metals. While the spatial distribution of the sampling points was not sufficient for in-depth observations, it is clear from Figure 3 that the western part of Baniyas suffers from higher concentrations of heavy metals, while the northern and southern parts had fairly low values. The middle part showed a high variation in most of the maps due to the intense sampling points and the high variation between the readings.

Descriptive Statistics
The physicochemical water quality parameters (pH, temperature, electric conductivity (EC), total dissolved solids (TDS), turbidity, and biochemical oxygen demand (BOD) are statistically described and summarized in Table 1. The statistical values of each parameter were compared with UAE standards to obtain a general overview of the quality of the household drinking water. The reported pH for all of the samples was within the standard ranges. However, some of the samples detected quite low pH values of 6.33 (in phase-1) and 6.15 (in phase-2). However, they are considered close to standard limits (Table 1). For the EC values, the ranges were within the standards, as no samples detected EC values over 250. The values of the coefficient of variation reflected that the pH, temperature, EC, and TDS exhibited a minor variation, and were all below the standard limits. However, turbidity and BOD showed a high variation, reaching 63% and 55% respectively for phase-1, and 44% and 42%, respectively, for phase-2, although the values were below the standard limits. In general, the ranges of physical water quality parameters indicate that the quality of drinking water was acceptable with respect to its color and taste, and that its physicochemical quality complied with limits set by the Regulation and Supervision Bureau (RSB).
Heavy metals concentrations were determined for the 30 houses to assess the quality of the drinking water. The statistical description of all of the individual parameters of the heavy metals revealed a wide variation (Table 1). Based on the average concentrations of Cr, Cu, and Zn, these parameters were below the standard limits recommended by the RSB [3]. However, the maximum concentrations of As, Cd, and Pb that were detected were considerably higher than the recommended limits. Although the maximum values for these heavy metals exceeded the standard limits, the median values indicated that most of the samples had concentrations lower than set limits. Almost all of the heavy metals' concentrations exhibited a very high variation in their datasets, as the coefficient of variation ranges from 70% to 145% for phase-1 and 50% to 180% for phase-2. On the other hand, concentrations of cations parameters (Ca, Na, K, Mg, Ba) were far below the standard limits, and their variation was very low.

Principal Component Analysis
To explore the association between water quality parameters and contaminants, a principal component analysis (PCA) was conducted to reduce the original dataset to certain influencing factors [19,36]. To sidestep the measuring unit differences and numerical ranges, the dataset was standardized using the z-scale transformation technique. Then, the influential factors were defined as those encompassing at least 70% of the total variance. The results of the PCA for the two phases are presented in Table 2. The PCA results show that for phase-1, four principal components were obtained, accounting for 71% of the total variance.  Phase-2 PCA resulted in five principal components that encompassed 70% of the total variance.

Cluster Analysis
Cluster analysis was conducted for physicochemical parameters and heavy metal concentrations in the collected samples. The results show three clusters (C1, C2, and C3) that were formed from the analysis of both sampling phases. For samples taken in phase-1, the results of cluster analysis (Figure 4) show that for the first cluster C1, around 43% of the investigated houses (houses No. 6, 13, 14, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 29, and 30) had similarity in 98.9% of the water quality measurements. These houses are located in the northern and southern part of Baniyas, and their water distribution network could be forked from the same water main. The third cluster C3 only comprises houses No. 16 and 28, which represent the middle section of the study area; their water quality parameters had similarity in more than 99% of the measurements. The third cluster C3 encompasses 50% of the study locations in Baniyas. The C3 cluster characterizes the western and the middle part of Baniyas, and shows a similarity in 98% of the water quality measurements.
Health risk assessment was carried out to investigate the effect of heavy metals on adults and children. The HQ ranges for each heavy metal for adults and children are presented in Table 3; while the hazard index (HI) for each house is shown in Table 4. To avoid exaggerating the risk assessment, the analysis did not rely on the maximum concentration, as it could have been observed only once during the study period; the other concentrations (the mean and the minimum) could be more representative, and were used alongside the maximum in the risk assessment. Based on the categories of risk assessment, each water quality parameter was evaluated. For both adults and children, the HQ values for Cd, Cr, and Cu range from negligible (<0.1) to low (<1) chronic risk. In addition, the HQ for Zn is shown to be of negligible risk (<0.1) with no effect on the health of the adults and children. In contrast, the HQs for As and Pb range from low (<1) to high (>4.0), indicating their adverse health effects on the human body and potential non-carcinogenic concern via daily oral intake. The HI results revealed the overall risk of the heavy metals. The results show a significant hazard in approximately 30% of the houses studied in Baniyas (houses No. 5, 6, 15, 18, 20, 23, 25, 29, and 30), and categorized these houses at a very high risk (HI > 8). Moreover, about 55% of the houses experienced a moderate to high risk on human health for both adults and children. The remaining

Health Risk Assessment
Health risk assessment was carried out to investigate the effect of heavy metals on adults and children. The HQ ranges for each heavy metal for adults and children are presented in Table 3; while the hazard index (HI) for each house is shown in Table 4. To avoid exaggerating the risk assessment, the analysis did not rely on the maximum concentration, as it could have been observed only once during the study period; the other concentrations (the mean and the minimum) could be more representative, and were used alongside the maximum in the risk assessment. Based on the categories of risk assessment, each water quality parameter was evaluated. For both adults and children, the HQ values for Cd, Cr, and Cu range from negligible (<0.1) to low (<1) chronic risk. In addition, the HQ for Zn is shown to be of negligible risk (<0.1) with no effect on the health of the adults and children. In contrast, the HQs for As and Pb range from low (<1) to high (>4.0), indicating their adverse health effects on the human body and potential non-carcinogenic concern via daily oral intake.
The HI results revealed the overall risk of the heavy metals. The results show a significant hazard in approximately 30% of the houses studied in Baniyas (houses No. 5, 6, 15, 18, 20, 23, 25, 29, and 30), and categorized these houses at a very high risk (HI > 8). Moreover, about 55% of the houses experienced a moderate to high risk on human health for both adults and children. The remaining 15% of the houses had a low to moderate risk on human health. Overall, compared with adults, the higher HQ and HI values for children were observed, suggesting that children are more sensitive when exposed to heavy metals in water.

Discussion
This study examined the presence of six potentially toxic contaminants (arsenic, cadmium, lead, chromium, copper, and zinc) in household drinking water. Only two contaminants, arsenic and lead, present a concern for public health due to their frequent occurrence and high concentrations (Tables 1 and 3). Thus, this discussion will focus on the presence of lead and arsenic in drinking water. Many studies have investigated the effects of Pb on human health; the reports show that Pb has a severe impact on the hemopoietic, nervous, endocrine, and cardiovascular systems in the human body. For arsenic, it was shown that long-term exposure from drinking water can cause cancer in the skin, lungs, bladder, and kidneys [32,[37][38][39].
The spatial analysis of heavy metals (Figure 3) shows that the drinking water contamination is not regionally specific within the district of Baniyas; further analysis shows that the contamination is rather household-specific. While there were some differences in the principal components extracted from the two phases, there were some striking similarities, especially in the results for the heavy metals. Phase-1 PC4 and phase-2 PC2 both show the association between temperature and Pb concentration. The consistent patterns of some parameters (turbidity, cations, and heavy metals) in the PCA of the two phases suggest relatively strong associations between these parameters in Baniyas's drinking water. However, caution is warranted in the interpretation of the PCA for health consequences due to concerns of sampling adequacy for the collected data.
The results of clustering can help in planning for water quality monitoring in Baniyas. It may not be necessary to sample each house; sampling one house from each cluster can give the necessary information about the range of contamination in the cluster. This can reduce the cost and time of the study. It also helps and compliments the rapid assessment of the impact of any taken measures to enhance household water quality.
The presence of lead in drinking water is commonly attributed to plumbing components containing lead (mainly solders) that were previously used and are still available in the market [40]. These components, which are commonly known as 'end-of-line fittings', can particularly influence the water quality at the faucet. The presence of heavy metals (arsenic, cadmium, and lead) at high concentrations in the water samples may be a result of corrosion in the main pipes or household plumbing. However, the spatial and statistical analysis suggests that the source is from the water mains. This calls for a broader investigation of the water quality in the distribution mains.

Limitations of This Study
Some studies draw attention to the possibility of data skewness (bias) resulting from voluntary participation [40,41]. The premise of this claim is that volunteer participants will be willing to participate in the study when they are suspicious about the quality of their drinking water. This could result in most of the collected samples having high levels of contamination.
The multivariate statistical analysis would be more beneficial with a larger sample size [42] so that a detailed analysis of the contaminant in the drinking water samples could be carried out. In this study, the sample's size was limited due to several reasons. (a) There was no readily available data on the supplied area number of households and population at the time of designing the study. (b) Only 30% of the proposed number of households agreed to participate in this study. (c) Due to budget limitations, the choice had to be made between the number of contaminants that each sample was analyzed for, and the total number of households sampled.
The study did not examine the impact of the Baniyas distribution system's characteristics such as pipe materials, diameters, and network layout on drinking water contamination. It would be interesting to investigate the statistical relations between these characteristics and the presence of contaminants in drinking water.

Recommendations for Managing Contaminants in Drinking Water
Most environmental and health agencies advise consumers to flush the taps before consumption for two to three minutes [41,43,44]. This technique will help reduce the contaminant concentration (particularly lead); however, around 2000 L would be wasted annually per household, which is against the country's policy for reducing water consumption.
Other studies suggest using a tap made from a material that does not contain lead (lead-free taps), such as stainless steel taps. These taps could prevent lead dissolution in the plumbing system [41,45]. Yet, some studies reported the incidence of lead presence in water samples even when lead-free taps were used, which was attributed to in-line brass fittings within the household plumbing system.

Conclusions
In the present study, various water quality parameters (physicochemical parameters and heavy metals) were assessed for water samples collected at four sampling points from 30 houses in the Baniyas region over two seasons. The analyzed data was interpreted by utilizing two main techniques: spatial (qualitative) analysis and multivariate statistical techniques. The main findings are that for both phases, the arsenic, cadmium, and lead concentrations were higher than the maximum allowable levels set in UAE standards. In addition, the water main samples had the highest concentration of the heavy metals in comparison with the other sampling points. Overall, PCA revealed that some physical parameters could contribute to the elevated levels of heavy metals. Furthermore, a cluster analysis suggested that for the future rapid assessment of contamination status and evaluating the impact of water quality interventions, only one house from each cluster needed to be sampled/monitored. This could help in reducing the time and cost associated with water quality enhancement measures.
Health risk assessment results revealed that for As and Pb, the HQs ranged from low (<1) to high (>4.0) chronic risk, indicating their adverse health effects on the human body and potential chronic concern via daily oral intake. In addition, approximately 30% of the houses studied in Baniyas were categorized as houses at a very high risk. Overall, the results of the health risk assessment concluded that the HQ and HI values for children were observed to be higher than the adults, suggesting that children are more sensitive when exposed to heavy metals in water.
To this end, it is recommended to conduct further investigations into the concentrations of heavy metals to warrant intervention by the government to improve the water quality from the distribution mains and reduce health risks. Acknowledgments: The authors would like to thank Abu Dhabi Food Control Authority (ADFCA) for providing assistance in the microbial analysis done in this study. We also thank Mujahid Ali for helping with the water sampling and quality analysis. We are also grateful to the Departments of Geology and Chemistry at the United Arab Emirates University for their assistance in the physicochemical analysis of water samples.

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

Appendix A
Correlation Matrices have been produced in order to explore the relation between the physico-chemical parameters.

Appendix B
Appendix B represents the spatial variation of the heavy metals and cations have been analyzed using ArcGIS interpolation tool. Appendix B represents the spatial variation of the heavy metals and cations have been analyzed using ArcGIS interpolation tool.