Next Article in Journal
Azima tetracantha Leaf Methanol Extract Inhibits Gastric Cancer Cell Proliferation through Induction of Redox Imbalance and Cytochrome C Release
Next Article in Special Issue
Removal of Copper, Nickel, and Zinc Ions from an Aqueous Solution through Electrochemical and Nanofiltration Membrane Processes
Previous Article in Journal
Wettability of Carbonate Reservoir Rocks: A Comparative Analysis
Previous Article in Special Issue
Prediction of the Long-Term Performance Based on the Seepage-Stress-Damage Coupling Theory: A Case in South-to-North Water Diversion Project in China
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Applied Sciences
Volume 12
Issue 1
10.3390/app12010132
applsci-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Groundwater Quality Monitoring Using In-Situ Measurements and Hybrid Machine Learning with Empirical Bayesian Kriging Interpolation Method

by
Delia B. Senoro
1,2,3,4,*,
Kevin Lawrence M. de Jesus
2,3,4,
Leonel C. Mendoza
4,5,
Enya Marie D. Apostol
4,5,
Katherine S. Escalona
4,6 and
Eduardo B. Chan
7
1
School of Civil, Environmental and Geological Engineering, Mapua University, Intramuros, Manila 1002, Philippines
2
School of Graduate Studies, Mapua University, Intramuros, Manila 1002, Philippines
3
School of Chemical, Biological, Materials Engineering and Sciences, Mapua University, Intramuros, Manila 1002, Philippines
4
Resiliency and Sustainable Development Center, Yuchengco Innovation Center, Mapua University, Manila 1002, Philippines
5
College of Teacher Education, Mindoro State University, Calapan 5200, Oriental Mindoro, Philippines
6
College of Arts and Sciences, Mindoro State University, Victoria 5205, Oriental Mindoro, Philippines
7
Dyson College of Arts and Sciences, Pace University, New York, NY 10038, USA
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(1), 132; https://doi.org/10.3390/app12010132
Submission received: 5 December 2021 / Revised: 17 December 2021 / Accepted: 17 December 2021 / Published: 23 December 2021
(This article belongs to the Special Issue Water Quality Modelling, Monitoring and Mitigation)
Browse Figures
Versions Notes

Abstract

:

Featured Application

In-Situ and Hybrid Machine Learning—Geostatistical Interpolation method for groundwater quality monitoring applications.

Abstract

This article discusses the assessment of groundwater quality using a hybrid technique that would aid in the convenience of groundwater (GW) quality monitoring. Twenty eight (28) GW samples representing 62 barangays in Calapan City, Oriental Mindoro, Philippines were analyzed for their physicochemical characteristics and heavy metal (HM) concentrations. The 28 GW samples were collected at suburban sites identified by the coordinates produced by Global Positioning System Montana 680. The analysis of heavy metal concentrations was conducted onsite using portable handheld X-Ray Fluorescence (pXRF) Spectrometry. Hybrid machine learning—geostatistical interpolation (MLGI) method, specific to neural network particle swarm optimization with Empirical Bayesian Kriging (NN-PSO+EBK), was employed for data integration, GW quality spatial assessment and monitoring. Spatial map of metals concentration was produced using the NN-PSO-EBK. Another, spot map was created for observed metals concentration and was compared to the spatial maps. Results showed that the created maps recorded significant results based on its MSEs with values such as 1.404 × 10−4, 5.42 × 10−5, 6.26 × 10−4, 3.7 × 10−6, 4.141 × 10−4 for Ba, Cu, Fe, Mn, Zn, respectively. Also, cross-validation of the observed and predicted values resulted to R values range within 0.934–0.994 which means almost accurate. Based on these results, it can be stated that the technique is efficient for groundwater quality monitoring. Utilization of this technique could be useful in regular and efficient GW quality monitoring.

1. Introduction

Water quality is associated with ecosystem preservation, economic growth and social development [1]. Groundwater (GW) quality is critical to the Philippines’ overall water resource; hence, monitoring should be given attention. Population expansion and the acceleration of modernization as well as industrialization have resulted in an increased water demand [2]. It is inevitable that quality of both surface water and GW is compromised in areas where economy is in transition and where there is increasing urbanization, industrialization, and agricultural activities [3]. However, due to scarcity of elemental laboratory instruments in the Philippines, the preparation requirements of the laboratory station-based instruments for samples, and the travel time from sampling sites to laboratory stations become a challenge to GW quality monitoring especially to elemental concentration analysis.
The population of the Philippines has rice in its regular daily meals. Rice fields and techniques to produce high and good quality yields are among the agricultural areas and programs, respectively, that are being supported by the Philippine government. Calapan City, Oriental Mindoro province in the Philippines, with GW as primary water source, is among the top producers of rice in the country. Rice fields require highly flattened area and good water quality as among the criteria to increase yields with good rice quality. Having the Philippines within the tropical storm belt, Calapan City experienced regular flooding as shown in Figure 1. This condition become a challenge to the water quality for rice fields and for domestic supply.
The Philippines is rich in natural resources such as metals and non-metallic minerals [5] with tropical climate where annual rainfall is high. However, anthropogenic activities, due to economic development, population positive growth rate and urbanization, caused unintentional adverse effects to the environment. When the pristine environment is disturbed, and minerals are exposed to oxygen and water, chemical reactions happen. Similar condition takes place during weathering. Having rich in mineral resources, high annual rainfall, large flattened agricultural areas, and GW as primary source of water supply, GW quality monitoring is important. Water quality that has elevated concentration of metals such as arsenic (As), copper (Cu), iron (Fe) nickel (Ni), manganese (Mn), lead (Pb), and other metals known for its toxicity characteristics would have acute and/or chronic adverse effects to human health [6,7,8].
Several technologies for measurement of the presence of metals in GW exists such as Inductively Coupled Plasma—Mass Spectrometry (ICP—MS), Inductively Coupled Plasma—Optical Emission Spectrometry (ICP—OES), and Atomic Absorption Spectrometry (AAS). These approaches are laboratory-based and require several days before result of analysis is available. This condition become not suitable for field monitoring and measurements [9] in sub-urban, rural areas, areas where access is a challenge and when time is a primary criterion in the analysis. These conditions require in-situ measurements, and accurate detection become a critical component in monitoring the GW quality. In-situ measurements provide observations on a rapid phase as well as covering the wider areas especially those with difficulty in access. This is in contrast of laboratory—based methods with significant limitations such as expensive instrumentation making limited availability, complex sample preparation and applicability in field conditions [10]. There are on-site detection and monitoring techniques such as electrochemical analysis [10,11], cyclic voltammetry (CV) [12], anodic stripping voltammetry (ASV) [13], square wave anodic stripping voltammetry (SWASV) [14], electro-chemical impedance spectroscopy (EIS) [15], electrochemiluminescence (ECL) [16] and the use of piezoelectric biosensors [17]. However, these techniques have drawbacks such as background noise control, unable to fulfill the current requirements for selectivity [9], detection limits of CV [12], insolubility of metals and the multiple peaks of ASV [13], complicated interferences and complex matrices in SWASV [14], EIS’s inability to identify different ions [15], frequent fouling of electrodes in the case of ECL [16], and only few enzymes are sensitive to heavy metals for the case of biosensors [17]. Hence, the use of portable x-ray fluorescence spectrometry (pXRF) technique in onsite metals detection and analysis is appropriate in rugged condition yet provides user of accurate and rapid analysis. This is in contrast of laboratory-based methods with significant limitations such as expensive instrumentation, complex sample preparation and applicability in field conditions [18]. Therefore, non-destructive analytical technique, such as relatively simple spectra line void of many interferences and rapid multi-element analyses [19], contributes significantly to the successful implementation of this study.
Concentration maps were frequently used tool for spatial monitoring. Spatial information in water resources are limited, and GW quality data can be obtained only through spot sampling. However, this procedure often requires extensive manpower and resources [10]. The issue in this practice and the determination of sample locations density influences the accuracy of the generated spatial maps [19]. The integration of in–situ measurements and GIS—based spatial interpolation techniques offer an improvement in the presentation and display of the status of GW quality in an area. The use of this integrated approach provides a clear and intelligent base—maps which can be utilized by researchers, policy makers, implementors for planning proposals [20] and creation of strategic programs.
Several studies on GW monitoring and assessment implemented using GIS—based approaches focused on different Southeast Asian countries such as the Philippines [1,21], Thailand [22], Malaysia [23], Singapore [24], Indonesia [25], Cambodia [26], Laos [27], and Vietnam [28]. This article illustrates the quantification and mapping of the concentrations of heavy metals such as Ba, Cu, Fe, Mn, and Zn in GW and presents the utilization of in-situ GW quality monitoring that uses a hybrid machine learning-geospatial interpolation technique and pXRF. This type of technique and analyses give prompt, accurate data and information on the current GW quality. This is to address the challenges encountered on-ground during sampling activities, the scarcity of instruments in the Philippines due to its price, and the complex samples preparation required by some laboratory-based instruments. Analyzing heavy metal concentrations in a faster, accurate and convenient method can help the researchers, authorities, water utility companies and local government units in making prompt decision, guidelines and strategic programs.

2. Materials and Methods

2.1. Description of the Study Area

The study area is Calapan City, in the province of Oriental Mindoro, Philippines. This is a third-class city and one of only two cities in the MIMAROPA region of the Philippines. It is the capital of the island province of Mindoro and located on the island’s northeastern shore. It has a population of about 150,000 people (about 25,000 households) and 62 barangays (the smallest local government unit) [29]. The city lies within 13°22′ N Latitude and 121°9′ E Longitude and Mindoro is located approximately 13°11′ N Latitude and 121°53′ E Longitude south of Mainland Luzon. The island of Mindoro is popularly known of rice production. Calapan City has an area of 217.30 square kilometers. Deep and shallow wells, in addition to piped water supply, are currently the primary sources of water in the city.

2.2. Collection/Treatment of Groundwater Samples

The GW samples were collected from twenty-eight (28) suburban deep and shallow well sites following the USEPA SESDPROC-301-R3/SESDPROC-111-R4 [30] as shown in Figure 2. The GW samples were collected using stainless steel sampler and polyethylene (PE) bottles. The PE bottles were thoroughly pre-washed with Type 1 water. Each PE bottle was carefully labeled, sealed, and placed temporarily in coolers for metals concentration detection. This is in preparation for the detection of the presence of metals concentration in all collected GW samples.

2.3. Physicochemical and Metal Concentrations Analysis

Temperature, pH, electric conductivity (EC), and total dissolved solids (TDS) of the samples were determined onsite using a multi-parameter water analyzer (HANNA HI 9811-5) with HI1285-5 probe (electrode) and HI7007, HI70031, HI70032, buffer solutions for calibration [31]. The HI7007, HI70031, HI70032 solution were used for pH, EC and TDS calibration, respectively. While HI700661 solution was used for cleaning the electrode. The physicochemical values detected in groundwater were compared to the permissible limits specified in the 2017 Philippine National Standards for Drinking Water (PNSDW) [32] and the WHO Drinking Water Guidelines [33,34]. These water parameters were used in the hybrid machine learning technique.
The heavy metal concentration analysis employed the use of portable handheld Olympus Vanta X-Ray Fluorescence Spectrometry. This pXRF is a rapid onsite accurate elemental analyzer that could be used for various environmental media including water [35,36,37,38,39]. The pXRF was set on geochem mode and recorded metals concentration in ppm (mg/L) detected from GW samples. Target metals were Ba, Cu, Fe, Mn and Zn.

2.4. Spatial Concentration Mapping Using Machine Learning Informed Empirical Bayesian Kriging (EBK) Method

Hybrid machine learning—geostatistical interpolation (MLGI) method was employed. The EBK technique was used to produce spatial concentration maps of the GW’s physicochemical properties and heavy metal concentrations. By sub-setting and replicating observed data, the EBK automates the most time-consuming aspects of constructing a viable kriging model. EBK provides a distribution of semi-variogram models and compensates for semi-variogram estimate uncertainty. The EBK is more realistic and superior to other current geostatistical modeling methods owing to its dependence on limited maximum likelihood estimation. This is in contrast to other existing kriging models that rely on weighted least squares estimation. EBK has many significant benefits, including a low need for interactive modeling, more accurate prediction of standard error and projection for small datasets as compared to other traditional kriging techniques, and exact prediction of substantially non-static data [40].
The Artificial Neural Network (ANN) approach is a subset of methods for artificial intelligence inspired by biological neurons. It is capable of quickly acquiring patterns and forecasting the result of a problem in a multi-dimensional environment. ANN models are trained using datasets [41] to show the efficacy. The training algorithm and the transfer function that was utilized in the model are two critical components of the ANN model. The Levenberg—Marquardt (LM) algorithm was chosen as the training algorithm since it is the quickest function for training a network, and the hyperbolic tangent sigmoid function was used as the transfer function because it is the recommended transfer function for rapid processes [42,43].
Particle Swarm Optimization (PSO) is a population-based stochastic optimization technique inspired by biological communities’ collaborative nature. The PSO is initiated using a community of randomly generated particles as solution options. It looks for global optima via iterations in which particles with their own velocity fly around the search space following the current optimum particles, which is the best approach for finding the best solution. The PSO was integrated to the ANN to determine the weights and biases which gives the minimum error [44].
This hybrid technique integrated to the EBK method generated the spatial concentration maps of the target study area. The Neural Network—Particle Swarm Optimization (NN-PSO) approach was applied to generate the spatial concentration maps of physicochemical parameters and HM concentrations.

3. Results

Subsequent sections elaborate the results of the study and in comparison, of the WHO and PNSDW guidelines.

3.1. Physicochemical Groundwater Parameters

The recorded physical and chemical properties of groundwater of the 28 sampling points are shown in Table 1, and in comparison to WHO (2017) and PNSDW (2017) guidelines. The detailed description of the sampling locations of the study area with the physicochemical properties of GW were exhibited in Table A1 of Appendix A.
The GW temperatures ranged from 26.2 to 33.6 degree Celsius which could lead to increased release rates of metals concentration especially within the water temperature of 30–35 degrees Celsius [47]. Furthermore, the study of Zhu, et al. in 2010 [48] attributed the high temperature in GW to the boom of urbanization that was also observed in the City of Calapan. The recorded pH of GW ranged from 6.7 to 8.8 which is within the pH range guidelines set by the WHO and PNSDW [49]. The release rates of metals were affected by a lower water pH. Lower pH of water means acidic water and known to be aggressive, enhancing the breakdown of Fe and Mn resulting in an unpleasant taste in water [50]. This condition could have adverse effects including heavy metal poisoning and toxicity [51,52,53]. The majority of the water samples are slightly basic which could be attributed to the existence of carbonates and bicarbonates [54]. The TDS and EC range recorded was 40–900 ppm and 100–1820 µS/cm, respectively. The TDS in GW found to be below WHO guidelines; however, beyond PNSDW guidelines. The TDS and EC had been found to have positive correlation [55]. The elevated EC of 1820 µS/cm has been attributed to inorganic chemicals in ionized form in water [54] such as metal elements.

3.2. Heavy Metal Concentrations

Presence of heavy metals were investigated in GW samples collected from the 28 sampling sites as indicated in Figure 2. Detected concentrations were compared to the existing maximum allowable limit of the WHO and the PNSDW 2017. These limits are enumerated in Table 2. The toxicants found in the GW samples are discussed in more detail in the subsections below.

3.2.1. Barium

All sampling locations observed Ba concentrations below permissible limits (Figure 3) of WHO, USEPA and PNSDW. The presence of Ba in GW has been attributed to the weathering of rocks such as igneous rocks, sandstone, shale, and coal [57,58].

3.2.2. Copper

The Cu concentrations (Figure 4) in GW samples from all sampling locations did not exceed the WHO guideline of 1.3 mg/L. The Cu in trace amounts in GW was associated to the kind of rock that forms the aquifer [59]. Another possible source of copper is the pipeline. Also, Cu concentrations at all sampling sites were within the acceptable limit of PNSDW (1 mg/L) and WHO guidelines.

3.2.3. Iron

Iron stains laundry and plumbing fixtures at concentrations more than 0.3 mg/L; it can also give metallic taste [60]. Hence, USEPA set an allowable concentration limit of 0.3 mg/L. The majority of the Fe in GW comes from minerals and sediments which may be in the form of particulate or dissolved [61]. The Fe concentration in each sampling location is presented in Figure 5. Sampling location 8 recorded an elevated Fe concentration. This was attributed to a longer residence time [62] which is associated to the type of subsurface (aquifer) that promotes longer residence time and creates opportunity for metals to react through chemical and physical weathering [63]. In addition, the area shown in Figure 2 illustrates the area of sampling point 8 of having lesser active wells. This condition also contributes to longer residence time of GW.

3.2.4. Manganese

Groundwater samples collected from all sampling locations did not exceed the WHO’s maximum permissible level for Mn concentration (Figure 6). The natural occurrence of Mn in GW can be influenced by several factors including TDS, GW level fluctuations, and the residence time. Agricultural operations and domestic wastewater are additional potential two sources of Mn that can adversely affect the GW quality [64,65].

3.2.5. Zinc

The highest concentration of Zn was recorded at Brgy. Gutad. However, this highest Zn concentration was within the WHO and PNSDW permissible limits. Several locations observed without Zn concentrations were at Brgy. Balingayan, Brgy. Maidlang, Brgy. Managpi, and Brgy. Personas. Zinc is naturally found in GW and the acidity affects the quality [64]; hence, it is important that monitoring is carried out. The acidity theory states that the higher the acidity (i.e., lower pH) of the water, the higher the Zn concentration. As observed in Table 1, GW samples from all sampling sites were slightly basic which explains the low concentration levels of Zn. The Zn concentration in each sampling location is presented in Figure 7.

3.3. Correlation Analysis

The correlations between the physicochemical characteristics were investigated using Pearson correlation analysis calculated through International Business Machine Statistical Package for Social Sciences (IBM SPSS). The r and p values were presented. The r value expresses the relationship between variables. The p value expresses the significance of the relationship. A lower p-value denotes statistical significance, whereas a higher p-value denotes the opposite. A negative correlation was found between pH and the other variables, while a positive correlation was found between the other parameters. At the 0.01 p level, all relationships were significant. The correlation matrix for the physicochemical parameters are presented in Table 3.
A substantial negative correlation was observed between pH and temperature (r = −0.665), pH and EC (r = −0.602), and pH and TDS (r = −0.603) which is similar to the findings of Abou Zakhem et al. in 2017 [66] and Sunkari and Abu in 2019 [67]. On the other hand, a substantial positive correlation was observed between temperature and EC (r = 0.657), temperature and TDS (r = 0.664), and EC and TDS (r = 0.995). This correlation values agreed to the findings of Wali et al. in 2021 [68].
Similar to physicochemical parameters, Pearson correlation analysis for the relationships between Ba, Cu, Fe, Mn, and Zn was also taken. Fe was positively correlated with Mn and Zn; and Mn was positively correlated with Zn. Positive substantial correlations between these metals indicated the same origin, are mutually dependent, and have similar transport characteristics [69]. The positive p with higher r values of this study illustrates relationship between metals; however, this relationship was not significant. The presence of these metals in GW is attributed to natural weathering of rocks. The correlation matrix for metals concentrations is shown in Table 4.

3.4. Spatial Concentration Mapping Using NN-PSO + EBK

The NN-PSO simulation was applied to accelerate the performance of the prediction capability of the EBK method. The simulation showed an excellent result as evident to the mean squared error (MSE) and correlation coefficient (R) values wherein the ideal value is 0 and 1, respectively. The NN-PSO simulation performed for the physicochemical parameters and heavy metal concentrations are presented in Table 5. Correlation plots of the R (validation) and R (testing) for the governing NN-PSO models of physicochemical parameters and heavy metals concentration are illustrated as Figure A1 of Appendix B.
The relationship between the number of neurons ranging from 1 to 30 and the corresponding AIC (Akaike Information Criterion) values obtained for the physicochemical parameters (temperature, pH, EC, and TDS) as well as the heavy metal concentrations are exhibited in Figure 8 and Figure 9, respectively. These figures represent the AIC values of all NN-PSO models for each hidden neuron that was simulated in this study. It was observed that the best models for the physicochemical parameters were determined from the 25, 29, 30, and 27 hidden neurons (HN) for temperature, pH, EC and TDS, respectively. The best models for Ba, Cu, Fe, Mn, and Zn were observed in 29, 29, 28, 29, 30 HN, respectively, for the heave metal concentrations.
The spatial concentration of the physicochemical parameters of GW in Calapan City, Oriental Mindoro was mapped using NN-PSO+EBK interpolation method. The highest temperature for GW recorded in the study area was 33.6 °C which was observed in Brgy. Parang. While the least temperature was observed in Brgy. Biga with recorded temperature of 26.2 °C. The highest pH for GW was observed at Brgy. Canubing I, with pH equal to 8.8. The lowest pH was detected at Brgy. Sto. Nino with pH of 6.7. The highest EC and TDS observed in Brgy. Ibaba West with EC and TDS value of 1820 μS/cm and 900 ppm, respectively. The least observed EC and TDS concentration was 100 μS/cm and 40 ppm, respectively which was recorded in Brgy. Sta. Rita. The spatial concentration of the physicochemical parameters of GW is shown in Figure 10.
The spatial concentration of the heavy metals of GW in Calapan City, Oriental Mindoro including Ba, Cu, Fe, Mn, and Zn was also mapped using the NN-PSO+EBK interpolation method. The heavy metal concentration maps generated using the NN-PSO+EBK method was presented in Figure 11.
The highest Ba concentration was measured in Brgy. Canubing I, where it was 7.9 times more than the background value for Ba measured in the research area. The average Ba concentration across all sample sites was 5.6 times more than the background value in the area of study. The highest concentrations of Cu were measured at several locations and recorded to be three times greater than the background concentration of copper. The mean concentration of Cu was found to be 2.1 times that of the background concentration. The Fe concentrations were found to be highest in Brgy. Gutad where it was found to be 6.1 times greater than the background value reported for the research region. Moreover, multiple sites were observed to exceed the WHO standards for Fe. These sites include Brgy. Camansihan, Brgy. Ibaba East, Brgy. Masipit, Brgy. Parang, Site 2 of Brgy. Personas, Brgy. San Vicente East, Brgy. Sta. Cruz, Brgy. Sta. Rita, and Brgy. Sto. Nino. Meanwhile, the mean concentration in Calapan City was just 0.5 percent more than the background level. Mn and Zn concentrations were highest in Brgy. Sto. Nino and Brgy. Gutad, respectively. However, these concentrations were still below the background concentration in the research region. The heavy metal concentration trend in the study area was observed to be Mn < Cu < Ba < Zn < Fe. Generally, these concentrations are within the WHO, USEPA and PNSDW limit except for Fe in several areas.

3.5. Cross Validation and Spot Sampling Evaluation Results

The predicted and observed values were compared to the NN-PSO+EBK method using the correctness measures to test the robustness of the predicted models. The results shown in Table 6 exhibit a robust and accurate result based on the R values close to 1. The cross-validation results suggested that all values provided more accurate spatial distribution for the study area. The cross-validation results are presented in Table 6.
A spot sampling analysis was performed using the data from the households in different barangays of Calapan City. A total of 21,559 households were utilized in the spot sampling analysis which is presented in Figure 12. The distribution of the households included in the spot sampling analysis per barangay is presented in Figure 13.
The spot sampling results was compared to the spatial concentration maps created in Figure 11. Table 7 exhibits the spot sampling comparison results for all heavy metals detected in the GW resources in Calapan City. The results showed that the created maps provided good results based on its MSE which is approaching zero when contrasted to the spot sampling values.
Considering each barangay, the MSEs for each element were also obtained as presented in Table A2 of Appendix C. Figure 14 presents the summary of the spot sampling comparison results in each barangay for all heavy metals considered in the study.

4. Discussion

Oriental Mindoro, an island province, is vulnerable to GW pollution and degradation due to natural and human activities. Due to structural disadvantages and characteristics such as smaller land area and population, insufficient natural resources, geographical distribution, and other global factors beyond domestic control, a small island economy is less resilient to the threat of GW deterioration and contamination than larger and more diverse economies [70].
Water plays a critical part in achieving the United Nations’ Sustainable Development Goals (SDGs). One of the problems that population has been experiencing is ensuring that everyone achieves SDG 6 (clean water and sanitation) which seeks to guarantee universal access to, and sustainable management of water. Continuous data integration and frequent monitoring remain to be critical components to achieving SDG 6 [71]. Hence, creating tools to aid in carrying out GW monitoring is significant. Tools such as the hybrid NN-PSO+EBK in making GW monitoring convenient to researchers and authorities are important.
Various heavy metals were detected at various sampling sites across Calapan City. The mean concentration of these metals in GW remained below the WHO and PNSDW acceptable levels. The recorded in-situ physicochemical characteristics were also compared to the WHO and PNSDW acceptable limits. The average GW temperature observed in the study area was 29.99 °C while the average GW pH observed was 7.69. Both figures are within the permissible range of WHO and PNSDW. One (1) sampling location exhibited pH value exceeding the maximum limit for pH of PNSDW. The average EC for the area of study was 560.36 μS/cm. This is within the permissible limit of the WHO. One (1) sampling location exhibited an EC value exceeding the WHO limit of 1500 μS/cm. The EC observation made in this location was categorized as Type II. The EC greater than 1500 μS/cm but less than 3000 μS/cm implies medium salts enrichment [72]. The average TDS is 277.15 ppm which was below the maximum allowable limit by the PNSDW and WHO. Though, TDS levels did not exceed the permissible limits set by WHO and PNSDW but were substantially lower or higher than the suggested TDS range of 600−1000 impairing palatability. Specifically, data of water samples from Ibaba West recorded TDS of 900. This number is at the high side of the limit which suggested the probability of impaired taste. On the other hand, the TDS concentrations recorded in Balingayan (60), Biga (50), Canubing (50), Comunal (90), Ibaba East Site 1 (160), Managpi (100), Personas (60), and Sta. Rita (40) were all significantly lower than the recommended range (600−1000 ppm) which may result to flat and insipid flavor.
The heavy metals concentration for GW was also observed in the study area. The Fe concentrations detected in multiple locations were above the WHO and PNSDW standards. However, the rest of the heavy metals detected were within the permissible limits. Water with elevated metals concentration has the potential to cause several public health issues. Health risks associated with elevated Fe in GW is probable. Pollutants entering the human body through drinking water have been shown in numerous studies to have detrimental health consequences for consumers. Micronutrients are essential in living organisms; however, elevated concentration adversely affects public health. Similar case with Mn which is necessary for humans; however, excessive quantities will have negative consequences. Neurological disorders, such as aberrant walking, ataxia, muscle hypotonicity, and a face devoid of lasting emotions, are frequently associated to Mn [64]. Dysfunction of liver was also reported [73]. Furthermore, excess Mn concentration has been demonstrated to produce neurotoxicity in infants receiving parenteral nourishment [74]. Excess Mn has been also linked to a lower level of IQ in children [64].
Meanwhile, asbestos-related cancer is believed to be caused by free radicals, which are produced by iron. Free radicals produced by iron can cause cancer by oxidizing DNA and causing DNA damage [75]. Additionally, elevated levels of Mn and Fe in drinking water have been associated to a decrease in birth weight in term-born infants [76]. Furthermore, since animals’ intestinal mucosa is highly porous, the fast absorption into the blood has been attributed to the Ba2+ ions which are rapidly absorbed from the gastrointestinal system and lungs. Moreover, it has been observed that Ba poisoning mostly affects the cardiovascular system; nevertheless, renal dysfunction has been documented as well [77].
The use of in-situ and hybrid machine learning—geostatistical methods are an integral part of data integration for GW quality monitoring. The impact of GW contamination in an island province had been a threat to public health especially when GW is used as primary source of domestic, agricultural and industrial water supply. Application of NN-PSO+EBK hybrid technique enables the establishment of spatial variability map of the contaminants that contributes to the depletion of GW quality. As a result, future undesired consequences could be avoided using this monitoring technique. The NN-PSO+EBK can offer periodic and long-term data that can be utilized for permanent monitoring of GW quality and risk assessments. Also, this tool can be utilized as early warning of GW quality for detrimental effects [78] by human activities and/or natural weathering.

5. Conclusions

An in-situ approach and hybrid MLGI, i.e., NN-PSO+EBK, was applied to assess and evaluate the GW quality in Calapan City, Oriental Mindoro, Philippines. Physicochemical characteristics and metals concentrations were detected onsite at various sampling locations. Generally, the physicochemical analysis of GW samples met the WHO and PNSDW guidelines. The average values for temperature, pH, EC, and TDS were within the permissible limits though few sampling locations exceeded the permissible limits of WHO and PNSDW. The pH of all samples was within the limits set by the PNSDW. Barangays Buhuan, Camansihan, Gutad, Ibaba East (Site 2), Ibaba West, Ilaya, Lazareto, Maidlang, Masipit, Nag-iba II, Pachoca, Parang, San Vicente East, Sta. Cruz, and Sto. Nino recorded elevated EC values. This was attributed to the addition of leachable salts. Also, the recorded TDS values suggested probable impaired palatability by having values significantly below the recommended range of 600−1000 ppm. Heavy metals analysis showed that only Fe detected in multiple GW samples had concentration above the WHO and PNSDW maximum permitted levels. This condition presents health concerns to the consumers. The record on Fe concentration in Brgy. Gutad GW samples were above WHO and PNSDW limit. Other GW samples recorded target metals concentration within the WHO and PNSDW permissible limits. The spot sampling analysis results showed that the generated maps by hybrid technique such as NN-PSO+EBK were reliable in describing the heavy metal concentration in the city of Calapan based on its MSE, R and AIC values.
This study is useful as a reference to providing techniques on gathering data for GW quality monitoring to help attain SDGs 6. It is suggested to conduct a study targeting other metals and regular monitoring of its concentration using this hybrid MLGI technique. Additionally, the regular monitoring is necessary to better understanding of the possible health consequences. Furthermore, a health risk assessment based on GW quality should be conducted. Another, preliminary interventions on GW quality control is necessary.

Author Contributions

Conceptualization, D.B.S.; methodology, K.L.M.d.J. and D.B.S.; software, K.L.M.d.J. and D.B.S.; validation, D.B.S., L.C.M., K.S.E. and E.M.D.A.; formal analysis, K.L.M.d.J. and D.B.S.; investigation, L.C.M., E.M.D.A., K.S.E. and D.B.S.; resources, D.B.S.; data curation, D.B.S., L.C.M., E.M.D.A., K.S.E. and E.B.C.; writing—original draft preparation, K.L.M.d.J., L.C.M. and E.M.D.A.; writing—review and editing, E.B.C., D.B.S. and K.L.M.d.J.; visualization, K.L.M.d.J.; supervision, D.B.S.; project administration, D.B.S.; funding acquisition, D.B.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Philippine Council for Health Research and Development of the Department of Science and Technology, Philippines.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data are contained in the manuscript.

Acknowledgments

This is to recognize the ‘in-kind’ support of Mapua University, Manila, Philippines, and the Calapan City local government units. In addition, the Philippines Mines and Geosciences Bureau by providing the hazard maps for MIMAROPA region.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table
Table A1. Observed Physicochemical Properties of Groundwater Samples.
Table A1. Observed Physicochemical Properties of Groundwater Samples.
Sampling No.BarangayLatitudeLongitudeTemp (°C)pHEC (µS/cm)TDS (ppm)
1Balingayan (Site 1)13.31903° N121.13432° E28.37.913060
2Balingayan (Site 2)13.32454° N121.13555° E28.58.413060
3Biga13.32791° N121.17312° E26.28.512050
4Buhuan13.31451° N121.22395°E29.27.6660320
5Camansihan13.33399° N121.22656° E30.77.81200590
6Canubing I13.35590° N121.14091° E27.18.813050
7Comunal13.31267° N121.16494° E27.27.920090
8Gutad13.35518° N121.25278° E32.47.1900440
9Ibaba East (Site 1)13.41517° N121.17836° E31.67.5350160
10Ibaba East (Site 2)13.41484° N121.17769° E32.57.5970480
11Ibaba West13.41478° N121.17676° E31.47.41820900
12Ilaya13.41181° N121.18548° E31.17.1780380
13Lazareto13.42972° N121.19940° E317990490
14Maidlang13.39711° N121.22727° E327.5600290
15Managpi13.32512° N121.19595° E28.17.5220100
16Masipit13.38917° N121.16190° E31.57.3570270
17Nag-iba II13.34643° N121.25301° E30.67.3820400
18Pachoca13.41061° N121.16840° E29.28.1690340
19Palhi13.37502° N121.20703° E30.27.4410200
20Panggalaan13.30148° N121.19908° E30.17.6140160
21Panggalaan13.30027° N121.20041° E28.48.3180150
22Parang13.40059° N121.21769° E33.67.7910450
23Personas (Site 1)13.30623° N121.14083° E28.78.314060
24Personas (Site 2)13.30930° N121.13945° E27.48.414060
25San Vicente East13.31045° N121.17980° E32.97.6750370
26Sta. Cruz13.31633° N121.23461° E29.77.3500240
27Sta. Rita13.35212° N121.13091° E29.77.710040
28Sto. Nino13.40712° N121.18545° E30.36.71140560

Appendix B

Applsci 12 00132 g0a1a 550Applsci 12 00132 g0a1b 550
Figure A1. R Value Plots for Validation and Testing Phase of the NN-PSO: (a) Temperature; (b) pH; (c) EC; (d) TDS; (e) Ba; (f) Cu; (g) Fe; (h) Mn; (i) Zn.
Figure A1. R Value Plots for Validation and Testing Phase of the NN-PSO: (a) Temperature; (b) pH; (c) EC; (d) TDS; (e) Ba; (f) Cu; (g) Fe; (h) Mn; (i) Zn.
Applsci 12 00132 g0a1aApplsci 12 00132 g0a1b

Appendix C

Table
Table A2. Spot Sampling Results in each Barangay of Calapan City.
Table A2. Spot Sampling Results in each Barangay of Calapan City.
Brgy. No.BarangayLatitudeLongitudeElev. (m)Average MSE
BaCuFeMnZn
1Balingayan13.3241° N121.1407° E13.20.00001510.00004120.00001040.00000070.0000138
2Balite13.4131° N121.1580° E7.60.00020070.00023810.00096900.00001100.0025765
3Batino13.3494° N121.2201° E7.80.00005780.00006240.00098730.00000280.0003248
4Bayanan I13.3679° N121.1685° E8.80.00055070.00002400.00011220.00000400.0002384
5Bayanan II13.3560° N121.1699° E12.50.00020420.00013880.00053540.00000420.0001622
6Biga13.3270° N121.1733° E14.10.00000800.00001080.00003960.00000060.0001686
7Bondoc13.3867° N121.2010° E165.50.00008440.00001330.00254900.00001100.0003627
8Bucayao13.3066° N121.1915° E17.30.00001460.00003160.00000470.00000390.0000776
9Buhuan13.3106° N121.1915° E12.80.00000690.00000810.00026980.00000030.0000034
10Bulusan13.4037° N121.2012° E28.40.00001270.00007300.00097660.00000100.0000702
11Calero13.4159° N121.1831° E9.40.00000790.00014240.00001260.00000200.0000081
12Camansihan13.3428° N121.2290° E7.10.00000420.00001830.00044830.00000070.0004709
13Camilmil13.4061° N121.1760° E8.40.00008270.00001290.00008280.00000590.0000046
14Canubing I13.3554° N121.1423° E8.60.00000880.00000430.00309140.00000550.0012961
15Canubing II13.3261° N121.1216° E12.90.00020420.00010330.00021910.00000180.0017159
16Comunal13.3075° N121.1606° E19.50.00000850.00006940.00005580.00000080.0000270
17Guinobatan13.3829° N121.1818° E9.70.00037230.00001990.00031310.00000730.0000374
18Gulod13.3433° N121.2073° E9.00.00045210.00001600.00074860.00000490.0002808
19Gutad13.3597° N121.2464° E7.20.00015610.00006870.00755500.00001080.0021064
20Ibaba East13.4149° N121.1788° E6.60.00000090.00000160.00002500.00000080.0000022
21Ibaba West13.4146° N121.1762° E5.90.00000270.00000180.00000100.00000010.0000021
22Ilaya13.4129° N121.1840° E8.70.00000400.00002320.00003720.00000020.0000020
23Lalud13.3993° N121.1739° E9.10.00036690.00001860.00022290.00000680.0000209
24Lazareto13.4286° N121.1995° E12.90.00002510.00003710.00003050.00000120.0000589
25Libis13.4149° N121.1847° E9.20.00000740.00010060.00000180.00000100.0000036
26Lumang Bayan13.4009° N121.1816° E6.80.00011560.00001130.00014950.00000420.0000077
27Mahal na Pangalan13.4082° N121.1502° E9.00.00042170.00024510.00579780.00000940.0056021
28Maidlang13.3883° N121.2339° E8.30.00004270.00001110.00174920.00000140.0000401
29Malad13.3396° N121.1588° E10.80.00010850.00001890.00025580.00000930.0001421
30Malamig13.3439° N121.1456° E10.70.00006400.00000770.00027690.00001100.0001567
31Managpi13.3282° N121.1997° E18.00.00003210.00000510.00004650.00000140.0000330
32Masipit13.3869° N121.1603° E6.60.00004840.00001620.00011010.00000020.0000723
33Nag-Iba I13.3400° N121.2721° E8.50.00008930.00003130.00326340.00001080.0010670
34Nag-Iba II13.3470° N121.2622° E11.00.00001540.00001900.00163470.00000600.0006289
35Navotas13.3739° N121.2486° E7.50.00043490.00001140.00651420.00000190.0018669
36Pachoca13.4108° N121.1677° E6.00.00004410.00002660.00004650.00000240.0001095
37Palhi13.3750° N121.2070° E18.60.00019920.00006330.00036300.00000460.0000934
38Panggalaan13.3012° N121.1990° E17.90.00001640.00002340.00004020.00000390.0000126
39Parang13.4035° N121.2182° E9.50.00001480.00002090.00012790.00000150.0000127
40Patas13.3452° N121.1222° E12.90.00092130.00012690.00058740.00000210.0000333
41Personas13.3083° N121.1438° E14.50.00000060.00003000.00002130.00000010.0000028
42Puting Tubig13.3470° N121.1887° E7.30.00031780.00007880.00032610.00000370.0001840
43San Antonio13.4259° N121.1956° E14.00.00000130.00013530.00005190.00000450.0000587
44San Rafael13.4216° N121.1911° E6.40.00000990.00021830.00001710.00000530.0000280
45San Vicente Central13.4120° N121.1787° E8.60.00000100.00000860.00000240.00000110.0000011
46San Vicente East13.4098° N121.1798° E9.30.00000090.00001060.00000440.00000140.0000004
47San Vicente North13.4138° N121.1785° E6.70.00000030.00000450.00000050.00000030.0000001
48San Vicente South13.4095° N121.1779° E7.60.00000720.00001150.00001230.00000310.0000018
49San Vicente West13.4124° N121.1765° E7.80.00000050.00000400.00000240.00000050.0000010
50Santa Cruz13.3169° N121.2370° E12.10.00012310.00001730.00066970.00000020.0001946
51Santa Isabel13.3654° N121.1577° E5.30.00039420.00002330.00055330.00000430.0000575
52Santa Maria Village13.4093° N121.1748° E6.30.00001540.00001150.00001670.00000400.0000019
53Santa Rita13.3489° N121.1303° E11.10.00001390.00003170.00001400.00000020.0000081
54Santo Niño13.4066° N121.1848° E8.20.00000040.00000750.00000570.00000040.0000003
55Sapul13.3651° N121.1885° E14.30.00084500.00001080.00039830.00001260.0001391
56Silonay13.3992° N121.2248° E8.80.00000040.00000040.00000110.00000000.0000007
57Suqui13.4177° N121.2040° E11.10.00000570.00004340.00017490.00000190.0000464
58Tawagan13.3712° N121.1448° E8.60.00006930.00001690.00190520.00000040.0012349
59Tawiran13.3950° N121.1680° E7.00.00021600.00002530.00026040.00000300.0001365
60Tibag13.4123° N121.1730° E8.20.00000180.00000400.00000270.00000060.0000014
61Wawa13.4025° N121.1453° E8.80.00046090.00024410.00858410.00000550.0063922

References

  1. Pineda-Pineda, J.J.; Martínez-Martínez, C.T.; Méndez-Bermúdez, J.A.; Muñoz-Rojas, J.; Sigarreta, J.M. Application of Bipartite Networks to the Study of Water Quality. Sustainability 2020, 12, 5143. [Google Scholar] [CrossRef]
  2. Mekonnen, M.M.; Hoekstra, A.Y. Four billion people facing severe water scarcity. Sci. Adv. 2016, 2, e1500323. [Google Scholar] [CrossRef] [Green Version]
  3. Jha, M.K.; Shekhar, A.; Jenifer, M.A. Assessing groundwater quality for drinking water supply using hybrid fuzzy-GIS-based water quality index. Water Res. 2020, 179, 115867. [Google Scholar] [CrossRef] [PubMed]
  4. Department of Environment and Natural Resources—Mines and Geosciences Bureau. Geohazard Maps. Available online: http://www.region4b.mgb.gov.ph/28-geohazard-maps/98-geohazard-maps (accessed on 25 October 2021).
  5. Philippine Statistics Authority. Compendium of Philippine Environment Statistics Component 2: Environmental Resources and their Use. Available online: https://psa.gov.ph/press-releases/id/163678 (accessed on 25 October 2021).
  6. Kumar, A.; Cabral-Pinto, M.; Kumar, A.; Kumar, M.; Dinis, P.A. Estimation of Risk to the Eco-Environment and Human Health of Using Heavy Metals in the Uttarakhand Himalaya, India. Appl. Sci. 2020, 10, 7078. [Google Scholar] [CrossRef]
  7. Kumar, A.; Jigyasu, D.K.; Kumar, A.; Subrahmanyam, G.; Mondal, R.; Shabnam, A.A.; Cabral-Pinto, M.M.S.; Malyan, S.K.; Chaturvedi, A.K.; Gupta, D.K.; et al. Nickel in terrestrial biota: Comprehensive review on contamination, toxicity, tolerance and its remediation approaches. Chemosphere 2021, 275, 129996. [Google Scholar] [CrossRef] [PubMed]
  8. Kumar, A.; Cabral Pinto, M.; Candeias, C.; Dinis, P.A. Baseline maps of potentially toxic elements in the soil of Garhwal Himalayas, India: Assessment of their eco-environmental and human health risks. Land Degrad. Dev. 2021, 10, 3984. [Google Scholar] [CrossRef]
  9. Hartzler, D.A.; Jain, J.C.; McIntyre, D.L. Development of a subsurface LIBS sensor for in situ groundwater quality monitoring with applications in CO 2 leak sensing in carbon sequestration. Sci. Rep. 2019, 9, 1–10. [Google Scholar]
  10. Bu, X.; Koide, K.; Carder, E.J.; Welch, C.J. Rapid analysis of residual palladium in pharmaceutical development using a catalysis-based fluorometric method. Org. Process Res. Dev. 2013, 17, 108–113. [Google Scholar] [CrossRef]
  11. Farghaly, O.A.; Hameed, R.A.; Abu-Nawwas, A.A.H. Analytical application using modern electrochemical techniques. Int. J. Electrochem. Sci. 2014, 9, 3287–3318. [Google Scholar]
  12. Kudr, J.; Richtera, L.; Nejdl, L.; Xhaxhiu, K.; Vitek, P.; Rutkay-Nedecky, B.; Hynek, D.; Kopel, P.; Adam, V.; Kizek, R. Improved electrochemical detection of zinc ions using electrode modified with electrochemically reduced graphene oxide. Materials 2016, 9, 31. [Google Scholar] [CrossRef]
  13. Estela, J.M.; Tomás, C.; Cladera, A.; Cerda, V. Potentiometric stripping analysis: A review. Crit. Rev. Anal. Chem. 1995, 25, 91–141. [Google Scholar] [CrossRef]
  14. Khadro, B.; Sikora, A.; Loir, A.S.; Errachid, A.; Garrelie, F.; Donnet, C.; Jaffrezic-Renault, N. Electrochemical performances of B doped and undoped diamond-like carbon (DLC) films deposited by femtosecond pulsed laser ablation for heavy metal detection using square wave anodic stripping voltammetric (SWASV) technique. Sens. Actuators B Chem. 2011, 155, 120–125. [Google Scholar] [CrossRef]
  15. Avuthu, S.G.R.; Narakathu, B.B.; Eshkeiti, A.; Emamian, S.; Bazuin, B.J.; Joyce, M.; Atashbar, M.Z. Detection of heavy metals using fully printed three electrode electrochemical sensor. In SENSORS; IEEE: New York, NY, USA, 2014; pp. 669–672. [Google Scholar]
  16. Cheng, L.; Liu, X.; Lei, J.; Ju, H. Low-potential electrochemiluminescent sensing based on surface unpassivation of CdTe quantum dots and competition of analyte cation to stabilizer. Anal. Chem. 2010, 82, 3359–3364. [Google Scholar] [CrossRef] [PubMed]
  17. Bansod, B.; Kumar, T.; Thakur, R.; Rana, S.; Singh, I. A review on various electrochemical techniques for heavy metal ions detection with different sensing platforms. Biosens. Bioelectron. 2017, 94, 443–455. [Google Scholar] [CrossRef]
  18. Borgese, L.; Zacco, A.; Bontempi, E.; Pellegatta, M.; Vigna, L.; Patrini, L.; Riboldi, L.; Rubino, F.M.; Depero, L.E. Use of total reflection X-ray fluorescence (TXRF) for the evaluation of heavy metal poisoning due to the improper use of a traditional ayurvedic drug. J. Pharm. Biomed. Anal. 2010, 52, 787–790. [Google Scholar] [CrossRef]
  19. De Jesus, K.; Senoro, D.B.; Dela Cruz, J.C.; Chan, E.B. A Hybrid Neural Network—Particle Swarm Optimization In-formed Spatial Interpolation Technique for Groundwater Quality Mapping in a Small Island Province of the Philippines. Toxics 2021, 9, 273. [Google Scholar] [CrossRef] [PubMed]
  20. El-Rawy, M.; Fathi, H.; Abdalla, F. Integration of remote sensing data and in situ measurements to monitor the water quality of the Ismailia Canal, Nile Delta, Egypt. Environ. Geochem. Health 2020, 42, 2101–2120. [Google Scholar] [CrossRef] [PubMed]
  21. Solis, K.L.B.; Macasieb, R.Q.; Parangat, R.C.; Resurreccion, A.C.; Ocon, J.D. Spatiotemporal Variation of Groundwater Arsenic in Pampanga, Philippines. Water 2020, 12, 2366. [Google Scholar] [CrossRef]
  22. Tiankao, W.; Chotpantarat, S. Risk assessment of arsenic from contaminated soils to shallow groundwater in Ong Phra Sub-District, Suphan Buri Province, Thailand. J. Hydrol. Reg. Stud. 2018, 19, 80–96. [Google Scholar] [CrossRef]
  23. Mogaji, K.A. Application of vulnerability modeling techniques in groundwater resources management: A comparative study. Appl. Water Sci. 2018, 8, 1–24. [Google Scholar] [CrossRef] [Green Version]
  24. Nistor, M.M.; Rahardjo, H.; Satyanaga, A.; Hao, K.Z.; Xiaosheng, Q.; Sham, A.W.L. Investigation of groundwater table distribution using borehole piezometer data interpolation: Case study of Singapore. Eng. Geol. 2020, 271, 105590. [Google Scholar] [CrossRef]
  25. Rusydi, A.F.; Onodera, S.I.; Saito, M.; Ioka, S.; Maria, R.; Ridwansyah, I.; Delinom, R.M. Vulnerability of groundwater to iron and manganese contamination in the coastal alluvial plain of a developing Indonesian city. SN Appl. Sci. 2021, 3, 1–12. [Google Scholar] [CrossRef]
  26. Lado, L.R.; Polya, D.; Winkel, L.; Berg, M.; Hegan, A. Modelling arsenic hazard in Cambodia: A geostatistical approach using ancillary data. Appl. Geochem. 2008, 23, 3010–3018. [Google Scholar] [CrossRef]
  27. Viossanges, M.; Pavelic, P.; Rebelo, L.M.; Lacombe, G.; Sotoukee, T. Regional mapping of groundwater resources in data-scarce regions: The case of Laos. Hydrology 2018, 5, 2. [Google Scholar] [CrossRef] [Green Version]
  28. Giang, P.Q.; Nakhapakorn, K.; Ussawarujikulchai, A. Effectiveness of different spatial interpolators in estimating heavy metal contamination in shallow groundwater: A case study of arsenic contamination in Hanoi, Vietnam. Environ. Nat. Resour. J. 2011, 9, 31–37. [Google Scholar]
  29. Agutaya, C.A.C.; Zamora, J.T. Developmental Projects in Calapan City, Philippines: Localization Perspectives. Am. J. Educ. Res. 2018, 6, 133–136. [Google Scholar]
  30. United States Environmental Protection Agency (U.S. E.P.A). Operating Procedure for In Situ Water Quality Monitoring (SESDPROC-111-R4). Available online: https://www.epa.gov/sites/default/files/2015-06/documents/Insitu-Water-Quality-Mon.pdf (accessed on 16 December 2021).
  31. Hanna Instruments. HI9811-5 Portable pH/EC/TDS/Temperature. Available online: https://www.hannainst.com/portable-ph-ec-tds-temperature-meter-hi9811-5.html (accessed on 16 December 2021).
  32. Magalona, M.L.; Peralta, M.M.; Lacsamana, M.S.; Sabularse, V.C.; Pelegrina, A.B.; De Guzman, C.C. Analysis of Inorganic Arsenic (As (III) and Total As) and Some Physicochemical Parameters in Groundwater Samples from Selected Areas in Bulacan, Batangas, and Laguna, Philippines. KIMIKA 2019, 30, 28–38. [Google Scholar] [CrossRef]
  33. Groover, K.D.; Izbicki, J.A. Elemental analysis using a handheld X-ray fluorescence spectrometer (No. 2016-3043). U.S. Geol. Surv. Fact Sheet 2016, 2015–3043. [Google Scholar] [CrossRef]
  34. Analytical Methods Committee AMCTB No. 89. Hand-held X-ray fluorescence spectrometry. Anal. Methods 2019, 11, 2498–2501. [Google Scholar] [CrossRef]
  35. Muramatsu, Y.; Izawa, R.; Nishioka, H.; Nogami, T. X-ray fluorescence analysis of dilute heavy-metals in water using a portable x-ray fluorescence spectrometer with the metal-adsorbent, tobermorite. X-sen Bunseki No Shinpo 2009, 40, 195–201. [Google Scholar]
  36. Pearson, D.; Weindorf, D.C.; Chakraborty, S.; Li, B.; Koch, J.; Van Deventer, P.; de Wet, J.; Kusi, N.Y. Analysis of metal-laden water via portable X-ray fluorescence spectrometry. J. Hydrol. 2018, 561, 267–276. [Google Scholar] [CrossRef]
  37. Zhou, S.; Yuan, Z.; Cheng, Q.; Zhang, Z.; Yang, J. Rapid in situ determination of heavy metal concentrations in polluted water via portable XRF: Using Cu and Pb as example. Environ. Pollut. 2018, 243, 1325–1333. [Google Scholar] [CrossRef]
  38. Pearson, D.; Chakraborty, S.; Duda, B.; Li, B.; Weindorf, D.C.; Deb, S.; Brevik, E.; Ray, D.P. Water analysis via portable X-ray fluorescence spectrometry. J. Hydrol. 2017, 544, 172–179. [Google Scholar] [CrossRef] [Green Version]
  39. Crocombe, R.A.; Leary, P.E.; Kammrath, B.W. Portable Spectroscopy and Spectrometry, Applications; John Wiley & Sons: Hoboken, NJ, USA, 2021; Volume 2, ISBN 978-111-963-640-3. [Google Scholar]
  40. Magesh, N.S.; Chandrasekar, N.; Elango, L. Occurrence and distribution of fluoride in the groundwater of the Tamiraparani River basin, South India: A geostatistical modeling approach. Environ. Earth Sci. 2016, 75, 1–16. [Google Scholar] [CrossRef]
  41. Shariati, M.; Mafipour, M.S.; Mehrabi, P.; Bahadori, A.; Zandi, Y.; Salih, M.N.; Nguyen, H.; Dou, J.; Song, X.; Poi-Ngian, S. Application of a hybrid artificial neural network-particle swarm optimization (ANN-PSO) model in behavior prediction of channel shear connectors embedded in normal and high-strength concrete. Appl. Sci. 2019, 9, 5534. [Google Scholar] [CrossRef] [Green Version]
  42. Kayabasi, E.; Ozturk, S.; Celik, E.; Kurt, H. Determination of cutting parameters for silicon wafer with a Diamond Wire Saw using an artificial neural network. Sol. Energy 2017, 149, 285–293. [Google Scholar] [CrossRef]
  43. Ozturk, S.; Kayabasi, E.; Celik, E.; Kurt, H. Determination of lapping parameters for silicon wafer using an artificial neural network. J. Mater. Sci.: Mater. Electron. 2018, 29, 260–270. [Google Scholar] [CrossRef]
  44. Lin, Y.H.; Hu, Y.C. Electrical energy management based on a hybrid artificial neural network-particle swarm optimization-integrated two-stage non-intrusive load monitoring process in smart homes. Processes 2018, 6, 236. [Google Scholar] [CrossRef] [Green Version]
  45. World Health Organization. Guidelines for Drinking-Water Quality, 3rd ed.; World Health Organization: Geneva, Switzerland, 2004; Volume 1, ISBN 92-4-154638-7. [Google Scholar]
  46. Department of Health—Food and Drug Administration. Administrative Order No. 2017-0010—Philippine National Standards for Drinking Water of 2017. Available online: https://www.fda.gov.ph/administrative-order-no-2017-0010-philippine-national-standards-for-drinking-water-of-2017 (accessed on 25 October 2021).
  47. Li, H.; Shi, A.; Li, M.; Zhang, X. Effect of pH, temperature, dissolved oxygen, and flow rate of overlying water on heavy metals release from storm sewer sediments. J. Chem. 2013, 2013, 434012. [Google Scholar] [CrossRef]
  48. Zhu, K.; Blum, P.; Ferguson, G.; Balke, K.D.; Bayer, P. The geothermal potential of urban heat islands. Environ. Res. Lett. 2010, 5, 044002. [Google Scholar] [CrossRef]
  49. Lomboy, M.; Riego de Dios, J.; Magtibay, B.; Quizon, R.; Molina, V.; Fadrilan-Camacho, V.; See, J.; Enoveso, A.; Barbosa, L.; Agravante, A. Updating national standards for drinking-water: A Philippine experience. J. Water Health 2017, 15, 288–295. [Google Scholar] [CrossRef] [PubMed]
  50. Edwards, K.A.; Classen, G.A.; Schroten, E.H.J. The Water Resource in Tropical Africa and Its Exploitation. Available online: https://www.ilri.org/publications/water-resource-tropical-africa-and-its-exploitation (accessed on 20 November 2021).
  51. Anyanwu, B.O.; Ezejiofor, A.N.; Igweze, Z.N.; Orisakwe, O.E. Heavy metal mixture exposure and effects in developing nations: An update. Toxics 2018, 6, 65. [Google Scholar] [CrossRef] [Green Version]
  52. Fenton, T.R.; Huang, T. Systematic review of the association between dietary acid load, alkaline water and cancer. BMJ open 2016, 6, e010438. [Google Scholar] [CrossRef] [PubMed]
  53. Jan, A.T.; Azam, M.; Siddiqui, K.; Ali, A.; Choi, I.; Haq, Q.M. Heavy metals and human health: Mechanistic insight into toxicity and counter defense system of antioxidants. Int. J. Mol. Sci. 2015, 16, 29592–29630. [Google Scholar] [CrossRef] [Green Version]
  54. Patil, V.T.; Patil, P.R. Physicochemical Analysis of Selected Groundwater Samples of Amalner Town in Jalgaon District, Maharashtra, India. E-J. Chem. 2010, 7, 111–116. [Google Scholar] [CrossRef] [Green Version]
  55. Oyem, H.H.; Oyem, I.M.; Ezeweali, D. Temperature, pH, electrical conductivity, total dissolved solids and chemical oxygen demand of groundwater in Boji-BojiAgbor/Owa area and immediate suburbs. Res. J. Environ. Sci. 2014, 8, 444. [Google Scholar] [CrossRef] [Green Version]
  56. Ramasamy, S. Barium in Drinking Water—Background document for development of WHO Guidelines for Drinking-water Quality. Available online: https://www.who.int/water_sanitation_health/water-quality/guidelines/chemicals/barium-background-jan17.pdf (accessed on 26 October 2021).
  57. Kroschwitz, J.I.; Howe-Grant, M. Kirk-Othmer Encyclopedia of Chemical Technology, 3rd ed.; Grayson, M., Eckroth, D., Eds.; John Wiley and Sons: New York, NY, USA, 1978; Volume 3, pp. 457–463. [Google Scholar]
  58. Miner, S. Air Pollution Aspects of Barium and Its Compounds; Contract No. Ph-22-68-25; Litton Systems, Inc.: Bethesda, MD, USA, 1969; p. 69. [Google Scholar]
  59. Pinto, V.M.; Hartmann, L.A. Flow-by-flow chemical stratigraphy and evolution of thirteen Serra Geral Group basalt flows from Vista Alegre, southernmost Brazil. An. Acad. Bras. Cienc. 2011, 83, 425–440. [Google Scholar] [CrossRef] [Green Version]
  60. Vieira, I.F.B.; Rolim Neto, F.C.; Carvalho, M.N.; Caldas, A.M.; Costa, R.C.A.; Silva, K.S.D.; Parahyba, R.D.B.V.; Pacheco, F.A.L.; Fernendes, L.F.S.; Pissarra, T.C.T. Water Security Assessment of Groundwater Quality in an Anthropized Rural Area from the Atlantic Forest Biome in Brazil. Water 2020, 12, 623. [Google Scholar] [CrossRef] [Green Version]
  61. Beyene, G.; Aberra, D.; Fufa, F. Evaluation of the suitability of groundwater for drinking and irrigation purposes in Jimma Zone of Oromia, Ethiopia. Groundw. Sustain. Dev. 2019, 9, 100216. [Google Scholar] [CrossRef]
  62. Carasek, F.L.; Baldissera, R.; Oliveira, J.V.; Scheibe, L.F.; Dal Magro, J. Quality of the groundwater of the Serra Geral aquifer system of Santa Catarina west region, Brazil. Groundw. Sustain. Dev. 2020, 10, 100346. [Google Scholar] [CrossRef]
  63. Maher, K. The dependence of chemical weathering rates on fluid residence time. Earth Planet. Sci. Lett. 2010, 294, 101–110. [Google Scholar] [CrossRef]
  64. Sholehhudin, M.; Azizah, R.; Sham, A.S.S.M.; Amiruddin, Z. Analysis of Heavy Metals (Cadmium, Chromium, Lead, Manganese, and Zinc) in Well Water in East Java Province, Indonesia. Malays. J. Med. Health Sci. 2021, 17, 146–153. [Google Scholar]
  65. Zhang, Z.; Xiao, C.; Adeyeye, O.; Yang, W.; Liang, X. Source and mobilization mechanism of iron, manganese and arsenic in groundwater of Shuangliao City, Northeast China. Water 2020, 12, 534. [Google Scholar] [CrossRef] [Green Version]
  66. Abou Zakhem, B.; Al-Charideh, A.; Kattaa, B. Using principal component analysis in the investigation of groundwater hydrochemistry of Upper Jezireh Basin, Syria. Hydrol. Sci. J. 2017, 62, 2266–2279. [Google Scholar]
  67. Sunkari, E.D.; Abu, M. Hydrochemistry with special reference to fluoride contamination in groundwater of the Bongo District, Upper East Region, Ghana. Sustain. Water. Resour. Manag. 2019, 5, 1803–1814. [Google Scholar] [CrossRef]
  68. Wali, S.U.; Alias, N.; Harun, S.B. Reevaluating the hydrochemistry of groundwater in basement complex aquifers of Kaduna Basin, NW Nigeria using multivariate statistical analysis. Environ. Earth Sci. 2021, 80, 1–25. [Google Scholar] [CrossRef]
  69. Senoro, D.B.; De Jesus, K.L.M.; Yanuaria, C.A.; Bonifacio, P.B.; Manuel, M.T.; Wang, B.N.; Kao, C.C.; Wu, T.N.; Ney, F.P.; Natal, P. Rapid site assessment in a small island of the Philippines contaminated with mine tailings using ground and areal technique: The environmental quality after twenty years. In IOP Conference Series: Earth and Environmental Science; IOP Publishing: Bristol, UK, 2019; Volume 351, p. 012022. [Google Scholar]
  70. Wright, N. Small Island Developing States, disaster risk management, disaster risk reduction, climate change adaptation and tourism. Background Paper for the Global Assessment Report on DRR 2013. 2013. Available online: https://www.preventionweb.net/english/hyogo/gar/2013/en/bgdocs/Wright,%20N.,%202013.pdf (accessed on 20 November 2021).
  71. Lopez-Maldonado, Y.; Batllori-Sampedro, E.; Binder, C.R.; Fath, B.D. Local groundwater balance model: Stakeholders’ efforts to address groundwater monitoring and literacy. Hydrol. Sci. J. 2017, 62, 2297–2312. [Google Scholar] [CrossRef]
  72. Galhardi, J.A.; Bonotto, D.M. Hydrogeochemical features of surface water and groundwater contaminated with acid mine drainage (AMD) in coal mining areas: A case study in southern Brazil. Environ. Sci. Pollut. 2016, 23, 18911–18927. [Google Scholar]
  73. Williams, M.; Todd, G.D.; Roney, N.; Crawford, J.; Coles, C.; McClure, P.R.; Garey, J.D.; Zaccaria, K.; Citra, M. Toxicological Profile for Manganese; Agency for Toxic Substances and Disease Registry (US): Atlanta, GA, USA, 2012. [Google Scholar]
  74. Erikson, K.M.; Thompson, K.; Aschner, J.; Aschner, M. Manganese neurotoxicity: A focus on the neonate. Pharmacol. Ther. 2007, 113, 369–377. [Google Scholar] [CrossRef] [Green Version]
  75. Bhasin, G.; Kauser, H.; Athar, M. Iron augments stage-I and stage-II tumor promotion in murine skin. Cancer Lett. 2002, 183, 113–122. [Google Scholar] [CrossRef]
  76. Grazuleviciene, R.; Nadisauskiene, R.; Buinauskiene, J.; Grazulevicius, T. Effects of elevated levels of manganese and iron in drinking water on birth outcomes. Pol. J. Environ. Stud. 2009, 18, 819–825. [Google Scholar]
  77. Kravchenko, J.; Darrah, T.H.; Miller, R.K.; Lyerly, H.K.; Vengosh, A. A review of the health impacts of barium from natural and anthropogenic exposure. Environ. Geochem. Health. 2014, 36, 797–814. [Google Scholar] [CrossRef] [PubMed]
  78. MacDonald, A.M.; Bonsor, H.C.; Ahmed, K.M.; Burgess, W.G.; Basharat, M.; Calow, R.C.; Dixit, A.; Foster, S.S.D.; Gopal, K.; Lapworth, D.J.; et al. Groundwater quality and depletion in the Indo-Gangetic Basin mapped from in situ observations. Nat. Geosci. 2016, 9, 762–766. [Google Scholar] [CrossRef]
Applsci 12 00132 g001 550
Figure 1. Flood Hazard Map of Calapan City [4].
Figure 1. Flood Hazard Map of Calapan City [4].
Applsci 12 00132 g001
Applsci 12 00132 g002 550
Figure 2. Map of the Study Area and Sampling Sites.
Figure 2. Map of the Study Area and Sampling Sites.
Applsci 12 00132 g002
Applsci 12 00132 g003 550
Figure 3. Concentration of Ba in GW samples.
Figure 3. Concentration of Ba in GW samples.
Applsci 12 00132 g003
Applsci 12 00132 g004 550
Figure 4. Concentration of Cu in GW samples.
Figure 4. Concentration of Cu in GW samples.
Applsci 12 00132 g004
Applsci 12 00132 g005 550
Figure 5. Concentration of Fe in GW samples.
Figure 5. Concentration of Fe in GW samples.
Applsci 12 00132 g005
Applsci 12 00132 g006 550
Figure 6. Concentration of Mn in GW samples.
Figure 6. Concentration of Mn in GW samples.
Applsci 12 00132 g006
Applsci 12 00132 g007 550
Figure 7. Concentration of Zn in GW samples.
Figure 7. Concentration of Zn in GW samples.
Applsci 12 00132 g007
Applsci 12 00132 g008 550
Figure 8. The AIC Values for Physicochemical Parameters.
Figure 8. The AIC Values for Physicochemical Parameters.
Applsci 12 00132 g008
Applsci 12 00132 g009 550
Figure 9. The AIC Values for Heavy Metals.
Figure 9. The AIC Values for Heavy Metals.
Applsci 12 00132 g009
Applsci 12 00132 g010 550
Figure 10. Physicochemical parameters map of Calapan City (a) Temperature, (b) pH, (c) EC, and (d) TDS.
Figure 10. Physicochemical parameters map of Calapan City (a) Temperature, (b) pH, (c) EC, and (d) TDS.
Applsci 12 00132 g010
Applsci 12 00132 g011a 550Applsci 12 00132 g011b 550
Figure 11. Heavy metal concentration map of Calapan City (a) Ba, (b) Cu, (c) Fe, (d) Mn, and (e) Zn.
Figure 11. Heavy metal concentration map of Calapan City (a) Ba, (b) Cu, (c) Fe, (d) Mn, and (e) Zn.
Applsci 12 00132 g011aApplsci 12 00132 g011b
Applsci 12 00132 g012 550
Figure 12. Concentration of Zn in GW samples.
Figure 12. Concentration of Zn in GW samples.
Applsci 12 00132 g012
Applsci 12 00132 g013 550
Figure 13. Distribution of the Number of Households included in the Spot Sampling Analysis.
Figure 13. Distribution of the Number of Households included in the Spot Sampling Analysis.
Applsci 12 00132 g013
Applsci 12 00132 g014a 550Applsci 12 00132 g014b 550
Figure 14. Summary of MSE per Barangay considering (a) Ba, (b) Cu, (c) Fe, (d) Mn, and (e) Zn.
Figure 14. Summary of MSE per Barangay considering (a) Ba, (b) Cu, (c) Fe, (d) Mn, and (e) Zn.
Applsci 12 00132 g014aApplsci 12 00132 g014b
Table
Table 1. The groundwater physical and chemical properties.
Table 1. The groundwater physical and chemical properties.
Sampling No.Temperature (°C)pHEC (µS/cm)TDS (ppm)
128.37.913060
228.58.413060
326.28.512050
429.27.6660320
530.77.81200590
627.18.813050
727.27.920090
832.47.1900440
931.67.5350160
1032.57.5970480
1131.47.41820900
1231.17.1780380
1331.07.0990490
1432.07.5600290
1528.17.5220100
1631.57.3570270
1730.67.3820400
1829.28.1690340
1930.27.4410200
2030.17.6140160
2128.48.3180150
2233.67.7910450
2328.78.314060
2427.48.414060
2532.97.6750370
2629.77.3500240
2729.77.710040
2830.36.71140560
WHO [45]30.06.5–8.54001000
PNSDW [46] 6.5–8.5-600
Table
Table 2. Permissible limits of metals in groundwater.
Table 2. Permissible limits of metals in groundwater.
Parameter, mg/LWHOUSEPA (2009)PNSDW 2017
Ba0.7 [56]2.000.70
Cu1.30 [45]1.301.00
Fe-0.301.00
Mn0.40 10.050.40
Zn-5.00-
1 Guidance value.
Table
Table 3. Correlation matrix of the physicochemical parameters.
Table 3. Correlation matrix of the physicochemical parameters.
ParameterTemppHECTDS
Temp1.000
pH−0.665 **1.000
EC0.657 **−0.602 **1.000
TDS0.664 **−0.603 **0.995 **1.000
** Correlation is significant at the 0.01 level (2-tailed).
Table
Table 4. Correlation matrix for the heavy metal concentration of groundwater samples.
Table 4. Correlation matrix for the heavy metal concentration of groundwater samples.
MetalBaCuFeMnZn
Ba1
Cu0.0551
Fe−0.136−0.0011
Mn−0.2350.0720.3201
Zn−0.0890.0130.870 **0.1901
** Correlation is significant at the 0.01 level (2-tailed).
Table
Table 5. NN-PSO Simulation Results.
Table 5. NN-PSO Simulation Results.
Hidden NeuronsNo. of ParticlesNo. of IterationsElapsed Time (sec)MSER
ValidationTesting
Temp2532000180.385590.012040.998780.99656
pH29102000119.160850.004340.990780.99039
EC3012000163.348820.001550.998510.99966
TDS2732000337.905920.000320.999250.99981
Ba2912000172.622022.44 × 10−60.988380.99286
Cu2932000366.128611.64 × 10−70.998260.99636
Fe2832000139.786930.000910.990100.99951
Mn2952000115.262851.34 × 10−70.982520.99725
Zn3012000142.606181.08 × 10−50.979450.99580
Table
Table 6. Cross Validation Results for NN-PSO+EBK simulation.
Table 6. Cross Validation Results for NN-PSO+EBK simulation.
CriteriaTemppHECTDSBaCuFeMnZn
R0.9890.9900.9340.9770.9940.9760.9850.9750.984
Table
Table 7. Spot Sampling Comparison Results.
Table 7. Spot Sampling Comparison Results.
CriteriaBaCuFeMnZn
MSE0.00014040.00005420.00062600.00000370.0004141
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Senoro, D.B.; de Jesus, K.L.M.; Mendoza, L.C.; Apostol, E.M.D.; Escalona, K.S.; Chan, E.B. Groundwater Quality Monitoring Using In-Situ Measurements and Hybrid Machine Learning with Empirical Bayesian Kriging Interpolation Method. Appl. Sci. 2022, 12, 132. https://doi.org/10.3390/app12010132

AMA Style

Senoro DB, de Jesus KLM, Mendoza LC, Apostol EMD, Escalona KS, Chan EB. Groundwater Quality Monitoring Using In-Situ Measurements and Hybrid Machine Learning with Empirical Bayesian Kriging Interpolation Method. Applied Sciences. 2022; 12(1):132. https://doi.org/10.3390/app12010132

Chicago/Turabian Style

Senoro, Delia B., Kevin Lawrence M. de Jesus, Leonel C. Mendoza, Enya Marie D. Apostol, Katherine S. Escalona, and Eduardo B. Chan. 2022. "Groundwater Quality Monitoring Using In-Situ Measurements and Hybrid Machine Learning with Empirical Bayesian Kriging Interpolation Method" Applied Sciences 12, no. 1: 132. https://doi.org/10.3390/app12010132

APA Style

Senoro, D. B., de Jesus, K. L. M., Mendoza, L. C., Apostol, E. M. D., Escalona, K. S., & Chan, E. B. (2022). Groundwater Quality Monitoring Using In-Situ Measurements and Hybrid Machine Learning with Empirical Bayesian Kriging Interpolation Method. Applied Sciences, 12(1), 132. https://doi.org/10.3390/app12010132

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