{"paper_id":"3e7a94ee-118a-4ab0-a79c-7857565eec98","body_text":"Statistical modeling of groundwater geochemistry in northeastern Brazil | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Statistical modeling of groundwater geochemistry in northeastern Brazil José Francisco de Francisco Oliveira Júnior, Alex Souza Moraes, and 7 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-1210837/v1 This work is licensed under a CC BY 4.0 License Status: Under Review Version 1 posted 4 You are reading this latest preprint version Abstract The objective is to characterize the hydrochemical and qualitative composition of the groundwater of the Beberibe aquifer located in the Northeast of Brazil (NEB) using geostatistical, geoprocessing, and multivariate analysis techniques. The Beberibe aquifer is located in the state of Pernambuco (PE), NEB, and is responsible for supplying water to the northern portion of the Recife Metropolitan Region (RMRE). The data collected came from 34 georeferenced mining and mineral water production units. From this, a Piper diagram was constructed and hydrogeological and hydrochemical spatial maps were generated. The data were submitted to descriptive statistics, geostatistics, and principal component analysis (PCA). The results showed that the artesian wells near the coast have a porous geological formation. In all artesian wells, the waters are of the mixed bicarbonate type and highly potable for human consumption. The coefficient of variation was medium to high for all hydrochemical variables. The variables showed a better fit to the Gaussian model, in the study of spatial dependence. In PCA, Cl was the one with the highest correlation with nitrate (NO 3 ). The anthropogenic influence on water quality at the center of RMRE, mainly on NO 3 concentrations, had a direct impact on water quality for human health. mineral water Beberibe aquifer geostatistics water quality Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Figure 8 1. Introduction Brazil has the largest freshwater reservoir in the world, equivalent to 14% of the total planet Earth, being a country of high relevance worldwide in the supply of mineral water (Tundisi, 2008; Kunrath et al., 2020). In recent years, bottled mineral water has grown significantly on a national scale and thus has become an attractive and profitable business in Brazil (DNPM, 2016; ABINAM, 2020 ). The worldwide consumption of mineral water reached 284.1 billion liters in 2016, while other beverages, for example, soft drinks reached 159.1 billion liters in the same period (Huang and Liu, 2017), due to these increases, the need for evaluations and analyzes of the quality of drinking water becomes essential and with high frequency (Shruti et al., 2020). The consumer bases his choice on extrinsic factors, such as brand and personal habits (Pacheco et al., 2018). Mineral waters are products differentiated by intrinsic factors, mainly by their aroma and flavor, which are classified according to the hydrochemical parameters that, in turn, result from the characteristics of the rocks and the aquifer (depth of the water table, as for circulation, confinement, and porosity) that the water flows until its capture and packaging (Dias and Bernardes, 2016; Guadayol et al., 2016; Pacheco et al., 2018). The quality of mineral water is highlighted by Nunes et al. (2018), who showed the importance of its determination, in which, the hydrochemical composition is the qualitative aspect of greatest importance to be investigated since it has a direct impact on human health. Bulia and Enzweiler ( 2018 ) evaluated the bottled water taken from 78 different sources in the state of São Paulo, Brazil, and found that the hydrochemical composition, in general, satisfactorily represented the groundwater in the aquifers corresponding to the bottled mineral water samples. Schwanz et al. (2016), highlighted the potential of concentrations of perfluorooctanoic acid (PFOA) classified in Group 2B (i.e., the possible presence of a carcinogen for humans), of bottled mineral water in Brazil. The statistical representation of the hydrochemical composition of mineral water is relatively representative. However, the complementary analysis of the spatial variability allows a careful assessment when investigating the distribution of mineral water in the regions and, thus, makes it possible to question the intrinsic and extrinsic factors that alter the hydrochemical composition of these waters. Studies like Gilabert-Alarcón et al. (2018), who evaluated the processes that affect the quality of groundwater, using hydrochemical tracers, as the main ionic composition and stable isotopes of water and nitrates, through the synthetic pollution index, allowed the spatial regionalization of groundwater quality. Concomitant with the use of geostatistical techniques that present satisfactory recent results in the spatialization of the hydrochemical composition of groundwater (Jang, 2015; Gu et al., 2017; Arshad et al., 2019 ; Bodrud-Doza et al., 2019 ). Bulia and Enzweiler ( 2018 ), highlighted that the groundwater sampling method based on bottled mineral water is a low-cost and viable alternative for determining the composition and hydrochemical quality of the water, without logistical and bureaucratic limitations with conventional sampling (laboratory) on a continental scale. Because of the above, the study aims to characterize the hydrochemical and qualitative composition of the groundwater of the Beberibe aquifer located in Northeastern Brazil via geostatistical, geoprocessing, and multivariate analysis techniques. 2. Material And Methods 2.1. Location and characterization of the study area The study was realized in the municipalities of Recife, Olinda, Paulista, São Lourenço da Mata, Camaragibe, and Paudalho, both municipalities are located in the state of Pernambuco (PE), Northeastern Brazil (NEB) (Figure 1 ). The municipality of Paudalho is located in the mesoregion of Zona da Mata, while Recife, Olinda, Paulista, São Lourenço da Mata, and Camaragibe are located in the mesoregion Metropolitana de Recife, the municipalities are located between the parallels 7°51'0.0'' S – 8°8'0.0'' S and between the meridians 35°11'0.0'' W – 34°54'0.0'' W. According to the Köppen-Geiger climate classification, the study region is in two climatic zones: “Am”―humid or sub-humid tropical climate and “Aw”―tropical climate with a dry winter season (Figure 1 ) - (Beck et al., 2018 ). [INSERT TO FIGURE 1 HERE] The research was conducted in the Beberibe aquifer, the main aquifer that supplies the northern part of the Metropolitan Region of Recife (RMRE), having been explored for over fifty years (CPRH, 2020). In the last decade, with the significant increase in the number of artesian wells, as well as the increase in industrial activities implemented, followed by population densification and the urbanization process, among others, there was a need to assess the quality of groundwater (CPRH, 2020). RMRE is located on a Deltaic Plain, which has distinct hydrogeological characteristics, composed of sediments of different origins: fluvial, marine, colluvial, and mangroves, which cover the coastal sedimentary basins of Pernambuco and Paraíba, respectively to the south and north of the region, separated by the structural divider Lineamento Pernambuco. Figure 2 illustrates the spatial distribution of hypsometry, to establish a detailed understanding of the Beberibe aquifer that contemplates the western portion of the Pernambuco coast, with a variation from 0 (coast of the state of Pernambuco) to 261 m of altitude. [INSERT TO FIGURE 2 HERE] 2.2. Data collect The in-situ data collection was carried out in 34 extraction units (i.e., artesian wells) that capture mineral water, georeferenced, based on the analyzes made available by the companies and following the rules of the Pernambuco Health Surveillance Agency (APEVISA, 2020). The representation of the spatial distribution of artesian wells is shown in Figure 1 . The record of the hydrochemical values were those made available on the label of each mineral water unit, according to electrical conductivity (EC, µS. Cm −1 ), hardness (Hn, mg. L −1 CaCO 3 ), hydrogen potential (pH, dimensionless), evaporative resistance (DR, mg. L −1 ), calcium (Ca, mg. L −1 ), magnesium (Mg, mg. L −1 ), sodium (Na, mg. L −1 ), potassium (K, mg. L −1 ), bicarbonate (HCO 3 , mg. L −1 ), sulfate (SO 4 , mg. L −1 ), chlorine (Cl, mg. L −1 ), iron (Fe, mg. L −1 ), manganese (Mn, mg. L −1 ), fluorine (F, mg. L −1 ), silicon (Si, mg. L −1 ), nitrate (NO 3 , mg. L −1 ), strontium (Sr, mg. L −1 ), zinc (Zn, mg. L −1 ) and barium (Ba, mg. L −1 ). Table 1 shows the maximum permitted values (MPVs) of the hydrochemical parameters of mineral water for human consumption, according to the Ministry of Health, Ordinance n°. 26, of January 15, 1999, and, of the Environmental Company of the State of São Paulo (CETESB) – (Ministério da Saúde, 1999; Ministério da Saúde, 2006; CETESB, 2016). Table 1 Hydrochemical parameters of water for human consumption and their maximum allowed values (MPVs) according to Brazilian legislation. Parameters Unit MPV Parameters Unit MPV 1 NO 3 mg. L −1 10 11 F mg. L −1 1.5 2 DR mg. L −1 1000 12 Zn mg. L −1 5 3 pH - 6 to 9 13 Fe mg. L −1 0.3 4 SO 4 mg. L −1 400 14 EC µS/cm 350 5 Na mg. L −1 200 15 Hn mg. L −1 CaCO 3 500 6 Ca mg. L −1 250 16 HCO 3 - - 7 Mg mg. L −1 100 17 Si - - 8 Cl mg. L −1 2 18 Sr - - 9 K mg. L −1 10 19 Ba mg. L −1 1 10 Mn mg. L −1 0.1 - - - 1 NO 3 : Nitrate; 2 DR: Evaporative resistance; 3 pH: Hydrogen potential; 4 SO 4 : Sulfate; 5 Na: Sodium; 6 Ca: Calcium; 7 Mg: Magnesium; 8 Cl: Chlorine; 9 K: Potassium; 10 Mn: Manganese; 11 F: Fluorine; 12 Zn: Zinc; 13 Fe: Iron; 14 EC: Electrical conductivity; 15 Hn: Hardness; 16 HCO 3 : Bicarbonate; 17 Si: Silicon; 18 Sr: Strontium; 19 Ba: Barium. 2.3. Piper's diagram The Piper diagram was constructed and applied in the classification and comparison of distinct groups of water, associated with dominant cations and anions, also called tri-linear diagrams, in this type of representation, the concentrations of each ion analyzed were given in percentages, as well, plotting the proportions of the main cations (Ca +2 , Mg +2 , Na + and K + ) and the main anions (HCO −3 , Cl − and SO −4 ) in two respective triangular diagrams, combining the information from the two triangles in a rhombus located between them, with the main types of chemical constituents of water (Piper, 1944; Anders et al., 2014 ; Yang et al., 2016). The proportions were plotted in triangular graphical representations and their scales, for the proportion of the variables, corresponding to 100%. The representations show the relative proportions of the main ions, but not their absolute concentrations. The procedure for plotting a point in Piper's diagram consisted of entering the percentage of certain cations and anions, in the triangular bases, and a line was drawn corresponding to the value of this percentage. The crossing of the lines drawn on the triangular bases for the cations and anions determined the point that is transcribed in the central diamond, where the hydrochemical classification of the water is. The software used for processing Piper's diagram was Grapher version 13. 2.4. Hydrogeological and hydrochemical map Based on the database provided by the Brazilian Institute of Geography and Statistics (IBGE), the hydrogeological (IBGE, 2018) and hydrochemical (IBGE, 2020) map of the study region was generated. The hydrogeological map of NEB cartographically represents the productivity of aquifers in this region, based on the values of flows and specific flows of 54,864 tubular wells. Such information is stored in a database, with the procedures used in the Map of Hydrogeological Domains/Subdomains of Brazil and the Hydrogeological Charters of Brazil to the Millionth (CPRM, 2007). The hydrochemical map of NEB's underground water sources brings together a collection of 10,478 hydrochemical analyzes, all from tubular wells and chemically homogeneous domains concerning potability, chemical facies, and the suitability of water for use in irrigation. The chemical reports were incorporated into a database (elaborated in Access), where they were classified in terms of potability, chemical facies, and the suitability of waters for use in irrigation, and represented in a shapefile format. These determinations were migrated and georeferenced in GeoMedia, where work was carried out to individualize chemically homogeneous zones, using geological, physiographic, and hydrogeological criteria, which allowed the demarcation of units that have similar characteristics within their limits. The data made available on the platform are in the shapefile format, with an established pre-classification, in which reclassification techniques by geoprocessing were applied to reclassify these shapefiles, for the area comprising the study. The shapefiles were processed using ArcGIS® software, version 10.2.2 (ESRI, Redlands, California, USA). 2.5. Statistical analysis The data obtained were subjected to analysis of descriptive statistics based on the mean, median, and coefficient of variation (CV, %), being classified as low when the CV < 12%; medium when 12% < CV < 24%, and high when CV > 24% (Warrick and Nielsen, 1980). Subsequently, the Kolmogorov-Smirnov normality test (KS test) was applied with a significance level set at p ≤ 0.01. The geostatistics analysis was performed based on the calculation of the classic semivariances (Eq. 1), which estimated the structure and the spatial dependence between the pairs of observations. $$\\text{γ}\\text{(h) = }\\frac{\\text{1}}{\\text{2N(h)}}\\text{ }{\\sum }_{\\text{i}\\text{ }\\text{=}\\text{ }\\text{1}}^{\\text{N(h)}}{\\left[\\text{Z}\\left({\\text{X}}_{\\text{i}}\\right)\\text{-}\\text{Z(}{\\text{X}}_{\\text{i}}\\text{ + h)}\\right]}^{\\text{2}}\\text{ }\\text{(}\\text{1}\\text{)}$$ where, γ(h) - is the experimental semi-variance estimator, obtained from the sampled values Z(X i ), Z(X i +h); N(h) - is the number of measured value pairs separated by the vector or delay distance; h - is the distance between sample pairs (i.e., it is the distance between two samples); Z(X i ) and Z(X i + h) - are the values of the i-th observation of the regionalized variable, collected at points X i and X i +h ( i = 1, ..., n), separated by the vector h. Spatial dependence was analyzed by adjusting the semivariogram based on the estimate of the semivariance using the GEO-EAS® program (England et al., 1989). The data were adjusted to the spherical, exponential, and Gaussian models (Eqs. 2, 3 , and 4), respectively, according to Deutsch et al. (1998). The spherical, exponential, and Gaussian models are referred to in the literature as transitive theoretical models and are more common for adjusting semivariograms (Gois et al., 2015). Spherical Model: $$\\text{γ}\\left(\\text{h}\\right)\\text{ =}\\left\\{\\begin{array}{c}\\begin{array}{c}\\\\ {\\text{C}}_{\\text{0}}\\text{ }\\text{+}\\text{ }\\text{C}\\text{ ·}\\end{array}\\\\ \\begin{array}{c}\\\\ {\\text{C}}_{\\text{0}}\\text{ }\\text{+}\\text{ }\\text{C}\\text{,}\\end{array}\\end{array}\\begin{array}{c}\\begin{array}{c}\\left[\\text{1.5}\\frac{\\text{h}}{\\text{a}}\\text{- 0.5}{\\left(\\frac{\\text{h}}{\\text{a}}\\right)}^{\\text{3}}\\right]\\text{,}\\end{array}\\\\ \\begin{array}{c}\\end{array}\\end{array}\\begin{array}{c}\\begin{array}{c}\\\\ \\text{ for}\\text{ }\\text{0 ≤ h ≤ }\\text{a}\\end{array}\\\\ \\begin{array}{c}\\\\ \\text{ }\\text{ }\\text{for}\\text{ }\\text{ }\\text{ }\\text{h > }\\text{a}\\end{array}\\end{array}\\right.\\text{ (2)}$$ Exponential Model: $$\\text{γ}\\left(\\text{h}\\right)\\text{ = }{\\text{C}}_{\\text{0}}\\text{ }\\text{+}\\text{ }\\text{C . }⌊\\text{1-}\\text{exp}\\left(\\text{-}\\frac{\\text{3h}}{\\text{a}}\\right)⌋\\text{ }$$ 3 Gaussian model: $$\\text{γ}\\left(\\text{h}\\right)\\text{ = }{\\text{C}}_{\\text{0}}\\text{ }\\text{+}\\text{ }\\text{C }\\text{·}\\left[\\text{1-}\\text{exp}\\left(\\text{-}\\frac{{\\text{(3h)}}^{\\text{2}}}{{\\text{a}}^{\\text{2}}}\\right)\\right]\\text{ (4)}$$ where, γ(h) - is the experimental semi-variance estimator; C 0 +C - is the sill (i.e., it is the nugget effect plus the variance dispersion, given by the acronyms C 0 and C, respectively); h - is the distance between sample pairs; a - is the range (m). The best models of the adjusted semivariograms were validated by the cross-validation of the Jack-Knifing test, in which the mean must be close to zero and the standard deviation close to 1 (e.g., Vauclin et al., 1983), the program used for that analysis was GEO-EAS® (England et al., 1989). The degree of spatial dependence (DSD) was classified according to Cambardella et al. ( 1994 ), which suggests strong dependence (St) (DSD < 25%); moderate dependence (Md) (DSD between 25 and 75%), and weak dependence (We): (DSD> 75%). For making the kriging maps, the Surfer 9 program (Golden Software, 2010) was used. For principal component analysis (PCA), 19 variables were admitted. Based on the principal components (PC), the covariance matrix was obtained to extract the eigenvalues that originate the eigenvectors. To identify the variables that presented correlation, the Kaiser criterion was used, considering the eigenvalues greater than 1.0, which generate components with a relevant amount of information contained in the original data (Kaiser, 1958). Pearson's correlation (r) was also performed for all variables, seeking to correlate with the PCA, to highlight the similarities between the variables. The program used for PCA and the correlation r was RStudio, version 3.6.1 (R Core Team, 2019). Linear regression models were established to justify and clarify the main hydrochemical characteristics that are correlated with nitrate (NO 3 ), a component that has a direct impact on human health. It was verified among the main variables, the one with the highest coefficient of determination (R 2 > 0.5) to reduce the main relevance and correlation with NO 3 . The analyzes were performed using the Minitab software, version 19. 3. Results And Discussion 3.1. Hydrogeological and hydrochemical characterization of the Beberibe aquifer Figure 3 shows spatially the hydrogeological characterization of the groundwater of the Beberibe aquifer. Most of the artesian wells in the mineral water production units are close to the coast, a densely populated region (Figures 1 and 3 ). The artesian wells close to the coast are of porous geological formation, with the lowest altitudes (Figure 2 ), with flow variations from medium to very high. From the western portion of the study region, there is a fissural geological formation, considering the higher altitudes, with greater rocky extension and flow variations from medium to very high (Figures 2 and 3 ). [INSERT TO FIGURE 3 HERE] There is a specific homogeneous flow occurring close to the coastal region, categorized as “very high productivity” (Figure 3 ). The obtained result corroborates with the studies carried out by Hu et al. (2017) who analyzed the hydrogeological configuration of groundwater for the entire Brazilian territory and highlighted good water retention capacity for recharging and discharging aquifers. The authors also pointed out that the aquifers close to the Atlantic Ocean play a role as a “wall or barrier”, because the water tables in these regions always remain at the same level as the sea surface at the limits of the coast, thus maintaining better retention, which will lead to homogenization of water flow. Although the original depth of the well is unknown, part of the water may come from shallower horizons (Pernambuco coast) and another part may come from the western region of the state, via deeper horizons, in superior topographic conditions. Such observations are in line with Silva et al. (2020), who studied the predominant hydrochemical mechanisms and processes that control the hydrochemistry of groundwater in a region of Acre, Northern Brazil, that regardless of the knowledge of the depth of the well, the origin of the water can come from different horizons, again highlighting the topographic conditions under the influence of a given region. Figure 4 highlights the hydrochemical characterization of the groundwater in the Beberibe aquifer. It appears that in all of the sampled points (artesian wells), the waters are of the mixed bicarbonate type (Bm) and high potability for human consumption (IBGE, 2020). However, one should not only consider the hydrochemical reactions that occur in the Beberibe aquifer, but also the anthropic actions that occur in RMRE. [INSERT TO FIGURE 4 HERE] The composition of groundwater can provide relevant information about interactions in aquifers since hydrochemistry is the result of specific reactions related to rock-water interactions that vary within an aquifer unit (Teramoto et al., 2020). Consequently, the interpretation of aquifer mineralogy and water chemistry data may indicate the transport of groundwater within and between aquifers. As recently highlighted by Teramoto et al. (2020), the hydrochemistry in the Guarani aquifer, due to its mixed bicarbonate composition, influences the geological formation of the porous type. This statement is confirmed in the present study since the geological composition is heterogeneous and presents an advanced state of decomposition, which implies a high variability of bicarbonates. 3.2. Piper's diagram The hydrochemical face is a term used to describe the percentage of the chemical composition of the water, resulting from the combined effect of the solution's kinetics, rock-water interactions, hydrogeological configurations, and sources of contamination. Piper's trilinear diagram (Piper, 1944) consists of a diamond and a pair of equilateral triangles, in which two triangles represent the anion and the cation, respectively, connected by a diamond-shaped diagram. Figure 5 shows the hydrochemical composition of the groundwater in the Beberibe aquifer. [INSERT TO FIGURE 5 HERE] The chemical composition of groundwater is mainly of the Na-K-Cl-SO 4 type (Figure 5 ), where Na-K corresponds to approximately 90% in the first triangle. The hydrochemical type Ca-Mg-Cl-SO 4 also appeared in Piper's diagram, to which, Ca-Mg corresponded to 20% of the composition. Na + is the dominant cation and Cl − is the dominant anion, while Ca 2+ , Mg 2+ , and HCO 3− are not predominant. The chemical facies Na-Cl indicate a low flow in the aquifer that allows the predominance of Na + as a result of the ion exchange process against Ca 2+ . There was a high concentration of Cl − in all samples, which suggests an evaporation process in the Beberibe aquifer. Sargazi et al. (2020), evaluated the quality of groundwater in the city of Zahedan, Sistan, and the province of Baluchestan, in Iran, and the results of the hydrochemical composition of the waters via the Piper diagram, are similar to the results obtained in the study, where the authors pointed out the predominance of Na-Cl-SO 4 and Ca-Mg-Cl-SO 4 , in the studied waters. Like Iqbal et al. (2018) and Pazand et al. (2018) studied the hydrochemical composition of groundwater in Abu Dhabi - United Arab Emirates and in the Ardestan basin in Central Iran, respectively, showed results similar to the study. 3.3. Descriptive statistics and geostatistics of the data Table 2 shows the temporal variability based on descriptive statistics (see section 2.5), for the hydrochemical parameters of the mineral water in the Beberibe aquifer. According to the criteria of Warrick and Nielsen (1980), the CV% presented values of medium too high for all the studied variables, in which the average variation (14.77%) was only for the hydrogen potential (pH), for the other variables had a high variation, CV > 29%. All parameters showed normal distribution, according to the KS test at the level of 1% probability. Table 2 Descriptive statistics of the hydrochemical parameters of mineral water. Variable Mean Median 20 Min 21 Max 22 CV 23 SD 24 A 25 K 1 NO 3 4.53 1.29 0.15 43.25 195.14 8.83 3.52 12.91 2 DR 66.17 57.90 20.00 124.98 42.13 27.87 0.61 -0.42 3 pH 5.31 5.41 4.10 7.70 14.77 0.78 1.32 2.77 4 SO 4 3.77 2.03 0.00 18.25 105.64 3.98 2.54 6.32 5 Na 9.21 8.69 5.00 19.94 29.44 2.71 2.04 6.59 6 Ca 0.90 0.61 0.00 2.66 77.37 0.70 0.96 0.12 7 Mg 1.44 1.00 0.40 5.30 77.84 1.12 2.04 4.03 8 Cl 12.27 11.97 0.12 26.51 40.00 4.91 0.29 3.34 9 K 4.29 4.65 0.00 10.45 81.45 3.50 0.17 -1.47 10 Mn 0.03 0.00 0.00 0.85 485.27 0.15 5.75 33.30 11 F 0.04 0.00 0.00 0.85 361.79 0.15 5.51 31.31 12 Zn 0.01 0.01 0.00 0.08 166.80 0.02 2.51 6.71 13 Fe 0.06 0.02 0.00 0.38 186.85 0.11 2.27 3.93 14 EC 86.88 78.70 40.00 197.00 41.64 36.18 1.68 3.14 15 Hn 9.68 7.13 4.86 30.50 62.85 6.08 1.82 3.36 16 HCO 3 8.94 8.48 0.00 22.87 71.13 6.36 0.40 -0.70 17 Si 13.43 13.99 0.00 43.66 66.33 8.91 0.89 2.71 18 Sr 0.01 0.01 0.00 0.05 143.13 0.01 1.85 3.28 19 Ba 0.04 0.01 0.00 0.23 138.00 0.06 1.62 2.78 1 NO 3 : Nitrate; 2 DR: Evaporative resistance; 3 pH: Hydrogen potential; 4 SO 4 : Sulfate; 5 Na: Sodium; 6 Ca: Calcium; 7 Mg: Magnesium; 8 Cl: Chlorine; 9 K: Potassium; 10 Mn: Manganese; 11 F: Fluorine; 12 Zn: Zinc; 13 Fe: Iron; 14 EC: Electrical conductivity; 15 Hn: Hardness; 16 HCO 3 : Bicarbonate; 17 Si: Silicon; 18 Sr: Strontium; 19 Ba: Barium; 20 Min: Minimum; 21 Max: Maximo; 22 CV: Coefficient of variation; 23 SD: Standard deviation; 24 A: Asymmetry; and 25 K: Kurtosis. The variability between the medium and high categories is mainly due to the hydrogeological and hydrochemical diversity, which in turn occurs in the geological formation of rocks and which has a direct impact on the chemistry of groundwater (IBGE, 2018; IBGE, 2020). The results obtained corroborate those mentioned by Yang et al. (2020), where they investigated the hydrochemical composition of groundwater to highlight its quality and assess its risk/suitability for irrigation purposes in the Ordos Basin, China, they identified high CV% for all variables studied, except for the pH, where the CV% presented a value referring to the medium variability. The highest CV% indicated that the ions were sensitive to the external environment, such as hydrological conditions, topography, and anthropogenic activities, as discussed by Yang et al. (2020). Another key factor is pH, which according to Liu et al. (2018), in a study of the dynamic characteristics of groundwater in the valley plain of the city of Lhasa (capital of the Autonomous Region of Tibet), the low CV% indicated that the pH of groundwater tends to be relatively stable, which differs from the present study, in which the pH was not stable and CV% is admitted with medium variability, according to the criterion of Warrick and Nielsen (1980) - (Table 2 ). The pH according to the standards established in Table 1 and compliance with the minimum values recorded in Table 2 , presented low standards, indicative of points of distribution of acidic mineral water for human health. High concentrations of NO 3 were observed in some mineral water distribution points, therefore, outside the limits established by Brazilian legislation (Table 1 ), probably caused by the socio-economic development of the study region, as highlighted earlier by Liu et al. (2018). The electrical conductivity (EC) showed high CV%, due to the oscillation of salinity in urban areas (natural effect due to the rock/soil factor, and anthropization), such variations are characteristics of low permeability and low hydraulic gradient, due to the dilution effect salts during the rainy season and evaporation in the dry season. The results similar to those obtained by Tlili-Zrelli et al. (2015), they studied the hydrogeochemistry of groundwater in an alluvial aquifer in northern Tunisia, where the zones of low permeability and low hydraulic gradient, occur due to the dilution effect of salts during the rainy season and evaporation during the dry season, which favors a high variability in the levels of salts diluted in the water due to the climate and the greater and/or less water availability. The analysis of the spatial variability for the hydrochemical parameters is shown in Table 3 . All the parameters evaluated, the Gaussian model was the one that best fit, which is indicated for studies of groundwater. The highlight for the coefficient R 2 with variations between 0.6 to 0.99. Table 3 Analysis of spatial variability followed by adjusted semivariograms for the hydrochemical parameters of mineral water. 20 Var 21 Mo 22 C 0 23 C 0 +C 24 a 25 R 2 26 C 0 27 DSD Jack-Knifing (C 0 +C) 28 m 29 SD 1 NO 3 Gau 2.000 10.500 25000.000 0.800 19.047 Ft -0.082 1.017 2 DR Gau 90.000 1396.000 4886.115 0.849 6.446 Ft 0.006 1.007 3 pH Gau 0.128 1.975 14440.107 0.933 6.481 Ft 0.070 0.998 4 SO 4 Gau 0.200 3.000 100.000 0.599 6.667 Ft -0.087 1.113 5 Na Gau 3.000 18.900 20000.000 0.728 15.873 Ft 0.012 1.044 6 Ca Gau 0.120 0.491 417.424 0.713 24.439 Ft -0.048 1.041 7 Mg Gau 0.136 1.292 3360.178 0.908 10.526 Ft -0.011 0.940 8 Cl Gau 5.000 28.000 13000.000 0.701 17.857 Ft 0.050 1.049 9 K Gau 4.000 22.000 5000.000 0.708 18.182 Ft 0.029 1.061 10 Mn Gau 1.50E-05 1.54E-04 16956.777 0.657 9.740 Ft 0.027 1.081 11 F Gau 1.00E-06 1.00E-02 1000.000 0.989 0.010 Ft -0.056 1.112 12 Zn Gau 2.00E-04 0.724 8861.172 0.935 0.027 Ft 0.049 1.032 13 Fe Gau 1.00E-6 0.002 169.741 0.917 4.234E-4 Ft 0.017 1.028 14 EC Gau 80.000 557.100 2468.172 0.961 14.360 Ft -0.057 1.070 15 Hn Gau 12.000 60.000 12500.000 0.790 20.000 Ft 0.017 1.071 16 HCO 3 Gau 25.000 110.000 4200.000 0.812 22.727 Ft 0.023 1.053 17 Si Gau 50.000 210.000 7000.000 0.905 23.810 Ft 0.049 1.004 18 Sr Gau 5.00E-06 7.30E-05 803.671 0.867 6.849 Ft -0.053 1.246 19 Ba Gau 8.00E-04 0.021 4588.202 0.630 3.810 Ft -0.010 1.080 1 NO 3 : Nitrate; 2 DR: Evaporative resistance; 3 pH: Hydrogen potential; 4 SO 4 : Sulfate; 5 Na: Sodium; 6 Ca: Calcium; 7 Mg: Magnesium; 8 Cl: Chlorine; 9 K: Potassium; 10 Mn: Manganese; 11 F: Fluorine; 12 Zn: Zinc; 13 Fe: Iron; 14 EC: Electrical conductivity; 15 Hn: Hardness; 16 HCO 3 : Bicarbonate; 17 Si: Silicon; 18 Sr: Strontium; 19 Ba: Barium; 20 Var: Variable; 21 Mo: Model; 22 C 0 : Nugget effect; 23 C 0 +C: Sill; 24 a : Range; 25 R 2 : Coefficient of determination; 26 C 0 /(C 0 +C): Degree of spatial dependence (%); 27 DSD: Class of the degree of spatial dependence; 28 m: Mean; 29 SD: Standard deviation; Gau: Gaussian; Ft: Strong. Yin et al. (2019), analyzed geostatistically the hydrochemical variations and the causes of NO 3 pollution of groundwater in a plain in southern Beijing, China. They highlighted the need for interactions between the transitive theoretical models to validate the best for each parameter. However, in the present study, cross-validation by Jack-Knifing (Vauclin et al., 1983) was applied, were all established semivariograms presented the mean close to 0 and the standard deviation close to 1, and thus used in kriging maps for each studied hydrochemical variable (Table 3 ). The DSD was strong for all the parameters analyzed, which indicates that the value of each parameter associated with a given location, was more similar to the value of a neighboring sample than to the rest of the location of the sample set, which reinforces the validation of the Gaussian model for the studied hydrochemical variables (Cambardella et al., 1994 ). Figure 6 illustrates the spatial distribution of the hydrochemical parameters of the mineral water sources in the Beberibe aquifer in the study area. Based on the semivariograms, kriging of the maps was performed, according to each model generated for the hydrochemical variables of the groundwater, implementing the nugget effect, level, and range of each model. [INSERT TO FIGURE 6 HERE] NO 3 concentrations occurred mainly on the coast of Pernambuco, with an emphasis on the municipalities of Recife and Olinda (Figure 6 ). In both municipalities, the NO 3 concentrations were higher, mainly because these municipalities have a high population density and have a high flow of tourists, vegetable cultivation near rivers (causing eutrophication), as well as the holding of major cultural events, such as Carnival in the city of Recife and Olinda, Brazil, considered cities that host the carnival festival with the largest flow of people in the world, with an average of 1 million people a day. Within this theme, one of the factors that imply high concentrations of NO 3 is caused by the high pollution generated by tourists and the consequent contamination of the main rivers, highlighting Capibaribe, Beberibe, and Tejipió. Also noteworthy is the fact that the city of Recife is among the 10 largest Brazilian metropolises, with a population of 1.555 million inhabitants, according to the IBGE - (IBGE, 2020). In this context, the rate of pollution generated in the urban centers of the metropolis, associated with the lack of sanitation in neighborhoods, provides for the contamination of rivers, which in turn contribute to recharge the aquifer, which implies the contamination of groundwater and the increase of NO 3 concentration. Therefore, there is still the spread of waterborne diseases, mainly diarrhea, typhoid fever, hepatitis A, leptospirosis, cholera, and intestinal infections, caused by bacteria. Concomitant with NO 3 concentrations, it is observed in the maps of the hydrochemical variables that the evaporative resistance (DR), sodium (Na), calcium (Ca), magnesium (Mg), chlorine (Cl), electrical conductivity (EC), hardness (Hn) and strontium (Sr), presented higher concentrations in the places with the highest concentration of NO 3 (Figure 6 ). From the coast that comprises the municipality of Recife and Olinda, it is where the greatest concentrations of these elements are found, among them are heavy metals and EC, which is associated with greater contractions of soluble solids in water, which in turn correlate with anthropic pollution. Corroborating the findings of the present study, Bodrud-Doza et al. ( 2016 ), carried out the characterization of groundwater quality based on water assessment indices, multivariate statistics, and geostatistics in central Bangladesh, and concluded that anthropogenic factors (e.g., industrial pollution, inadequate garbage discharge, inadequate sanitation) contribute groundwater quality. Bhuiyan et al. ( 2016 ) also assessed the groundwater quality of the Lakhimpur district in Bangladesh based on water quality indices, geostatistical methods, and multivariate analysis pointed out that a kriging method is an effective tool for decision-makers on groundwater quality management. It is noteworthy that the hydrogeological formation also contributes to NO 3 concentrations. It appears that the Pernambuco coast has a hydrogeological formation of the mixed bicarbonate type (Figure 4 ), which because of this, there is a high variability (> CV%) of the hydrochemical elements of groundwater (Table 2 ). This mixed diversity of bicarbonate can also be considered as one of the factors for this high concentration of NO 3 , as such regions are common of fissural geological formation, which in turn reflects and a greater tendency for contamination of these waters. According to the Ministério da Saúde (1999), Ministério da Saúde (2006), and CETESB (2016), concentrations of NO 3 , DR, Na, Ca, Mg, Cl, EC, Hn and Sr, only NO 3 and Cl did not comply with Brazilian legislation, with values higher than allowed (Table 1 ). However, to the west of the coast, values were observed within the limits allowed for NO 3 and Cl, which reinforces the anthropogenic impact in the coastal region (Figure 6 ). Corroborating the findings of Yin et al. (2019), who claim that industrial areas and intense urbanization promote greater groundwater pollution by NO 3 . The pH of groundwater on the coast met the acceptable standard for human health (pH between 6 and 9), with variations between acid to alkaline (pH between 6.4 and 7.6) (Ministério da Saúde, 1999; Ministério da Saúde, 2006; CETESB, 2016). However, the west region of the coast had a pH < 6 (Figure 6 ). Abu-alnaeem et al. ( 2018 ), evaluated the salinity and quality of groundwater in the Gaza coastal aquifer, Gaza Strip in Palestine, and highlighted that acidity and alkalinity close to neutralization, is derived from the dissolution of carbonate minerals in the form of bicarbonate (HCO 3− ). Corroborating the statement of the aforementioned authors, the coastal region of the present study, to the east, showed pH with variations between acid and alkaline, still because it presents mixed bicarbonate hydrogeological formation (Figure 4 ). The elements SO 4 , K, Mn, F, Zn, Fe, HCO 3 , Si, Sr and Ba were the elements that had the highest concentrations west of the coast (Figure 6 ). This study region was also the one with the lowest pH. It is noteworthy that only the elements SO 4 , F, Zn, and Ba, presented concentrations within the quality standards for human consumption (Ministério da Saúde, 1999; Ministério da Saúde, 2006; CETESB, 2016). The elements K, Mn, and Fe showed values above that allowed by current legislation. HCO 3 , Si, and Sr, on the other hand, do not have a defined quality standard. Hossain and Patra, (2020) carried out a contamination zoning and health risk assessment of trace elements in groundwater utilizing geostatistical modeling and reported that pH controls chemistry, which is a component that influences the hydrochemical standards of water and, thus, corroborates with the observations of the present study, in which the elements K, Mn, and Fe presented higher values for the region of lower pH. 3.4. Principal component analysis (PCA) The PCA presents the eigenvalue, the proportion, and the accumulation of the total variance for all variables studied (Table 4 ). The cumulative total variance for principal component 2 (PC2) was 49.61%, a significant result for the analyzed data set. According to the criteria established by Kaiser (1958), PC1 and PC2 presented their eigenvalues greater than 1, which indicated a significant information load and, thus, generate the biplot graphs. Table 4 Principal components of the hydrochemical parameters of mineral water. Components Eigenvalue Proportion (%) Cumulative (%) PC1 6.30 33.14 33.14 PC2 3.13 16.46 49.61 PC: Principal component. Yang et al. (2016) carried out the identification of hydrogeochemical processes and assessment of groundwater quality in the southeastern part of the Ordos basin in China, based on the PCA and obtained total cumulative variance for the PC2, of 59.26%, which is higher than the study. In Figure 7 , the PC’s biplot chart stands out and Pearson correlation of the hydrochemical parameters of the mineral water of the Beberibe aquifer in the study area, referring to the municipalities of Recife, Olinda, Paulista, Camaragibe, São Lourenço da Mata, and Paudalho. [INSERT TO FIGURE 7 HERE] For the PC's, the 19 studied variables (NO 3 ; DR; pH; SO 4 ; Na; Ca; Mg; Cl; K; Mn; F; Zn; Fe; EC; Hn; HCO 3 ; Si; Sr and Ba) were admitted and the correlated points refer to the artesian sample collection wells (A01 to A34). The variables DR, Na, Ca, Mg, Cl, EC, Hn, and Sr, showed a correlation with NO 3 , as observed in kriging maps (Figure 7 A). However, it is noteworthy that among the observed correlations, Cl was the one with the highest correlation with NO 3 , in which, the Cl projection line follows the NO 3 line, which characterizes this strong correlation (Figure 7 B). Similar results were found by Cao et al. ( 2020 ), where they studied the heterogeneity of an aquifer in a region of France, through PCA and observed a high correlation between Cl and NO 3 . The same results were found by Ayed et al. ( 2017 ), where they performed a hydrochemical characterization of groundwater using PCA for the southeastern region of Tunisia, and observed the same correlation between Cl and NO 3 . To clarify more precisely the correlations between the variables, Pearson's correlation was applied, where the correlation was significant between Cl and NO 3 (Figure 7 B). Corroborating the findings of this study, Ayed et al. ( 2017 ) also observed a strong correlation between Cl and NO 3 . 3.5. Linear regression analysis Based on the analysis of the kriging maps (Figure 6 ) and the PCA - (Figure 7 A), was performed a simple linear regression analysis of the electrical conductivity (EC), hardness (Hn), evaporative resistance (DR), magnesium (Mg), sodium (Na) and chlorine (Cl), with nitrate (NO 3 ), as they present a greater correlation with NO 3 in the previous analyzes (Figure 8 ). The R 2 criterion greater than 0.5 was established to represent only the most representative predictor variables with NO 3 . [INSERT TO FIGURE 8 HERE] Among the linear regressions, only NO 3 x DR showed R 2 < 0.5 (R 2 = 0.32), so the DR has no direct influence on NO 3 (Figure 8 C). It can be seen that Cl, showed a satisfactory R 2 with NO 3 (R 2 = 0.51) - (Figure 8 F), previously indicated by geostatistics based on kriging maps and PCA, and these two variables showed a high correlation between themselves. The results obtained corroborate those described by Tiwari et al. (2017), where they assessed the spatial distribution of groundwater quality and assessed regional suitability for drinking water in the Alwar region in India and found a strong correlation between NO 3 and Cl. The variables EC (R 2 = 0.68), Hn (R 2 = 0.57), Mg (R 2 = 0.67) and Na (R 2 = 0.56), showed satisfactory R 2 , respectively. Such variables have a direct correlation with NO 3 , thus being potential variables in the hydrochemical characterization of groundwater. 4. Conclusions Based on Brazilian legislation, there are strategic wells that have a direct impact on human health and, therefore, there is a need to reanalyze possible approaches for improvements that support health and quality of life. The composition and hydrochemical quality of groundwater in the metropolitan region of Recife and the coast of Pernambuco are directly influenced by hydrogeology, topographic conditions, and urban infrastructure (urban density). In geostatistical modeling, based on the Gaussian model, it stands out concerning the other transitive theoretical models in the hydrochemical characterization of the waters of the Beberibe aquifer and thus spatially displays the spatial dynamics of water composition and quality in the study area with greater accuracy. Multivariate analysis via PCA associated with Pearson's correlation is a promising technique for assessing the hydrochemistry of groundwater in the Beberibe aquifer together with the hydrogeological and hydrochemical maps of groundwater made it possible to efficiently characterize water quality standards for the coast of Pernambuco. The anthropogenic influence of the metropolitan center of Recife, has a direct impact on water quality for human health, especially on nitrate concentrations (NO 3 ). It is also noteworthy that the concentrations of chlorine (Cl) have a high correlation with NO 3 , being confirmed in the kriging maps and the PCA, Pearson's correlation, and in the linear regression. Other methodologies must be formulated and applied, to confirm or rule out the possibility of recharging the regional flow of groundwater. The methodology adopted in the study can be perfectly applied in other regions of the planet with conditions of social vulnerability and problems with water quality, with bias of monitoring and definition of public policies. Low-quality underground water should not be used for human consumption, as it becomes the main agent for the emergence of waterborne diseases, mainly diarrhea, typhoid fever, hepatitis A, leptospirosis, cholera, and intestinal infections, caused by bacteria, which it reaches from infants to the elderly, especially those who live in areas of social vulnerability and with low-quality water. Therefore, regular monitoring is essential to ensure good quality water and a healthy environment that provides quality of life for local families. Declarations 5. Acknowledgment To the Programa de Pós-Graduação em Engenharia Agrícola (PGEA), Programa de Pós-graduação em Química (PPGQ), and the Grupo de Pesquisa em Ambiência (GPESA) of the Federal Rural University of Pernambuco (UFRPE) for supporting the development of this research. The Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES - Finance Code 001) and the Fundação de Amparo à Ciência e Tecnologia do Estado de Pernambuco (FACEPE), for the financing of scholarships. The third and fourth authors are grateful for the Research Productivity grant from the Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) level 2 under process number 304060/2016-0 and 309681/2019-7, respectively. References ABINAM (2020) Associação Brasileira Indústria Águas Minerais. [online]. Available at:. [Accessed 15 July 2020] Abu-alnaeem MF, Yusoff I, Ng TF, Alias Y, Raksmey M (2018) Assessment of groundwater salinity and quality in Gaza coastal aquifer, Gaza Strip, Palestine: An integrated statistical, geostatistical and hydrogeochemical approaches study. Sci Total Environ 615:972–989. https://doi.org/10.1016/j.scitotenv.2017.09.320 Agência Pernambucana de Vigilância Sanitária – APEVISA (2020) Programa de Monitoramento e Fiscalização da Água Envasada (Pura). [online]. Available at:. [Accessed 11 January 2021] Anders R, Mendez GO, Futa K, Danskin WR (2014) A geochemical approach to determine sources and movement of saline groundwater in a coastal aquifer. Groundwater 52(5):756–768. https://doi.org/10.1111/gwat.12108 Arshad S, Mobeen M, Bashir S, Aziz T, Rehman A (2019) Spatializing Groundwater Quality Parameters and Their Impacts on Land Value in Khushab City, Punjab, Pakistan. International Journal of Economic and Environmental Geology 10(3):78–84. https://doi.org/10.46660/ijeeg.vol11.iss1.1919 Ayed B, Jmal I, Sahal S, Mokadem N, Saidi S, Boughariou E, Bouri S (2017) Hydrochemical characterization of groundwater using multivariate statistical analysis: the Maritime Djeffara shallow aquifer (Southeastern Tunisia). Environ Earth Sci 76(24):821. https://doi.org/10.1007/s12665-017-7168-6 Beck HE, Zimmermann NE, McVicar TR, Vergopolan N, Berg A, Wood EF (2018) Present and future Köppen-Geiger climate classification maps at 1-km resolution. Scientific Data 5:180214. https://doi.org/10.1038/sdata.2018.214 Bhuiyan MAH, Bodrud-Doza M, Islam AT, Rakib MA, Rahman MS, Ramanathan AL (2016) Assessment of groundwater quality of Lakshimpur district of Bangladesh using water quality indices, geostatistical methods, and multivariate analysis. Environ Earth Sci 75(12):1020. https://doi.org/10.1007/s12665-016-5823-y Bodrud-Doza M, Bhuiyan MAH, Islam SDU, Rahman MS, Haque MM, Fatema KJ, Rahman MA (2019) Hydrogeochemical investigation of groundwater in Dhaka City of Bangladesh using GIS and multivariate statistical techniques. Groundwater for Sustainable Development 8:226–244. https://doi.org/10.1016/j.gsd.2018.11.008 Bodrud-Doza M, Islam AT, Ahmed F, Das S, Saha N, Rahman MS (2016) Characterization of groundwater quality using water evaluation indices, multivariate statistics and geostatistics in central Bangladesh. Water Science 30(1):19–40. https://doi.org/10.1016/j.wsj.2016.05.001 Bulia IL, Enzweiler J (2018) The hydrogeochemistry of bottled mineral water in São Paulo state, Brazil. J Geochem Explor 188:43–54. https://doi.org/10.1016/j.gexplo.2018.01.007 Cambardella CA, Moorman TB, Novak JM, Parkin TB, Karlen DL, Turco RF, Konopka AE (1994) Field-scale variability of soil properties in central Iowa soils. Soil Science Society Amsterdam Journal 58:1501–1511. https://doi.org/10.2136/sssaj1994.03615995005800050033x Cao F, Jaunat J, Vergnaud-Ayraud V, Devau N, Labasque T, Guillou A, Ollivier P (2020) Heterogeneous behaviour of unconfined Chalk aquifers infer from combination of groundwater residence time, hydrochemistry and hydrodynamic tools. J Hydrol 581:124433. https://doi.org/10.1016/j.jhydrol.2019.124433 CETESB. (2016). QUALIDADE DAS ÁGUAS SUBTERRÂNEAS NO ESTADO DE SÃO PAULO . [online]. Available at: <https://cetesb.sp.gov.br/aguas-subterraneas/wp-content/uploads/sites/13/2013/11/Cetesb_QualidadeAguasSubterraneas2015_Web_20-07.pdf>. [Accessed 29 July 2020]. CPRH. (2020). RECURSOS HÍDRICOS SUBTERRÂNEOS. [online]. Available at: <http://www.cprh.pe.gov.br/downloads/25_Recursos_Hidricos_Subterraneos.pdf>. [Accessed 22 July 2020]. CPRM. (2007). Mapa de Domínios/Subdomínios Hidrogeológicos do Brasil 1:2.500.000 . [online]. Available at: <http://www.cprm.gov.br/publique/Hidrologia/Mapas-e-Publicacoes/Mapa-de-Dominios%7CSubdominios-Hidrogeologicos-do-Brasil-1%3A2.500.000-632.html>. [Accessed 1 July 2020]. Deutsch CV, Journel AG (1998) GSLIB Geostatistical Software Library and User’s Guide, Second edn. Oxford University Press, New York, p 369 England EJ, Sparks A, Robinson MD (1989) Geo—EAS (Geostatistical Environmental Assessment Software). Environmental Software 4(2):70–75. https://doi.org/10.1016/0266-9838(89)90033-6 Gilabert-Alarcón C, Daesslé LW, Salgado-Méndez SO, Pérez-Flores MA, Knöller K, Kretzschmar TG, Stumpp C (2018) Effects of reclaimed water discharge in the Maneadero coastal aquifer, Baja California, Mexico. Appl Geochem 92:121–139. https://doi.org/10.1016/j.apgeochem.2018.03.006 Gois G, Delgado RC, de Oliveira Júnior JF (2015) Modelos teóricos transitivos aplicados na interpolação espacial do Standardized Precipitation Index (SPI) para os episódios de El Niño forte no Estado do Tocantins, Brasil. Irriga 20(2):371–387. https://doi.org/10.15809/irriga.2015v20n2p371 GOLDEN SOFTWARE (2010) Surfer for windows version 9.0. Golden, Colorado, p 66 Gu X, Xiao Y, Yin S, Pan X, Niu Y, Shao J, Hao Q (2017) Natural and anthropogenic factors affecting the shallow groundwater quality in a typical irrigation area with reclaimed water, North China Plain. Environ Monit Assess 189(10):514. https://doi.org/10.1007/s10661-017-6229-3 Guadayol M, Cortina M, Guadayol JM, Caixach J (2016) Determination of dimethyl selenide and dimethyl sulphide compounds causing off-flavours in bottled mineral waters. Water Res 92:149–155. https://doi.org/10.1016/j.watres.2016.01.016 Hossain M, Patra PK (2020) Contamination zoning and health risk assessment of trace elements in groundwater through geostatistical modelling. Ecotoxicol Environ Saf 189:110038. https://doi.org/10.1016/j.ecoenv.2019.110038 Hu K, Awange JL, Forootan E, Goncalves RM, Fleming K (2017) Hydrogeological characterisation of groundwater over Brazil using remotely sensed and model products. Sci Total Environ 599:372–386. https://doi.org/10.1016/j.scitotenv.2017.04.188 Huang L, Liu Y (2017) Health information and consumer learning in the bottled water market. Int J Ind Organ 55:1–24. https://doi.org/10.1016/j.ijindorg.2017.08.002 IBGE (2018) Levantamento hidrogeologico. [online]. Available at: vN7H6tMg3R_3AI01ndqVTgh3mr000>. [Accessed 5 July 2020]. IBGE (2020) Hidroquímica subterrânea. [online]. Available at:. [Accessed 5 July 2020] Iqbal J, Nazzal Y, Howari F, Xavier C, Yousef A (2018) Hydrochemical processes determining the groundwater quality for irrigation use in an arid environment: The case of Liwa Aquifer, Abu Dhabi, United Arab Emirates. Groundwater for Sustainable Development 7:212–219. https://doi.org/10.1016/j.gsd.2018.06.004 Jang CS (2015) Geostatistical analysis for spatially characterizing hydrochemical features of springs in Taiwan. Environ Earth Sci 73(11):7517–7531. https://doi.org/10.1007/s12665-014-3924-z Kaiser HF (1958) The varimax criterion for analytic rotation in factor analysis. Psychometrika 23(3):187–200. https://doi.org/10.1007/BF02289233 Kunrath CCN, Patrocínio DC, Rodrigues MAS, Benvenuti T, Amado FDR (2020) Electrodialysis reversal as an alternative treatment for producing drinking water from brackish river water: A case study in the dry season, northeastern Brazil. Journal of Environmental Chemical Engineering 8(2):103719. https://doi.org/10.1016/j.jece.2020.103719 Liu J, Gao Z, Wang M, Li Y, Ma Y, Shi M, Zhang H (2018) Study on the dynamic characteristics of groundwater in the valley plain of Lhasa City. Environ Earth Sci 77(18):646. https://doi.org/10.1007/s12665-018-7833-4 Ministério da Saúde. (1999). PORTARIA Nº 26, DE 15 DE JANEIRO DE 1999 . [online]. Available at: <http://bvsms.saude.gov.br/bvs/saudelegis/anvisa/1999/prt0026_15_01_1999.html>. [Accessed 29 July 2020]. Ministério da Saúde. (2006). VIGILÂNCIA E CONTROLE DA QUALIDADE DA ÁGUA PARA CONSUMO HUMANO . [online]. Available at: <https://bvsms.saude.gov.br/bvs/publicacoes/vigilancia_controle_qualidade_agua.pdf>. [Accessed 29 July 2020]. Nunes LGDP, Oliveira MDV, de Souza AA, Lopes LDF, Dias PCES, Nogueira GB, Souza MAAD (2018) Water quality comparison between a supply network and household reservoirs in one of the oldest cities in Brazil. Int J Environ Health Res 29(2):173–180. https://doi.org/10.1080/09603123.2018.1531114 Pacheco MH, Esmerino EA, Capobiango CS, Cruz AG, Leddomado LS, Pimentel TC, de Freitas MQ (2018) Bottled mineral water: classic and temporal descriptive sensory analysis associated with liking. British Food Journal 120(7):1547–1560. https://doi.org/10.1108/BFJ-11-2017-0655 Pazand K, Khosravi D, Ghaderi MR, Rezvanianzadeh MR (2018) Identification of the hydrogeochemical processes and assessment of groundwater in a semi-arid region using major ion chemistry: A case study of Ardestan basin in Central Iran. Groundwater for Sustainable Development 6:245–254. https://doi.org/10.1016/j.gsd.2018.01.008 Piper AM (1944) A graphic procedure in the geochemical interpretation of water-analyses. Eos, Transactions American Geophysical Union 25(6):914–928. https://doi.org/10.1029/TR025i006p00914 R Core Team (2019) R: A language and environment for statistical computing. R Foundation for Statistical Computing Sargazi S, Mokhtari M, Ehrampoush MH, Almodaresi SA, Sargazi H, Sarhadi M (2020) The application of geographical information system (GIS) approach for assessment of groundwater quality of Zahedan city, Sistan and Baluchestan Province, Iran. Groundwater for Sustainable Development 100509. https://doi.org/10.1016/j.gsd.2020.100509 Schwanz TG, Llorca M, Farré M, Barceló D (2016) Perfluoroalkyl substances assessment in drinking waters from Brazil, France and Spain. Sci Total Environ 539:143–152. https://doi.org/10.1016/j.scitotenv.2015.08.034 Shruti VC, Pérez-Guevara F, Kutralam-Muniasamy G (2020) Metro station free drinking water fountain-A potential “microplastics hotspot” for human consumption. Environ Pollut 261:114227. https://doi.org/10.1016/j.envpol.2020.114227 Silva TA, Rosário FF, Silva Júnior GC (2020) Environmental isotopes and hydrochemical tracers applied to hydrogeological conceptual modeling of the southwest portion of the Amazon Aquifer System (Acre, Brazil). Appl Geochem 120:104670. https://doi.org/10.1016/j.apgeochem.2020.104670 Teramoto EH, Gonçalves RD, Chang HK (2020) Hydrochemistry of the Guarani Aquifer System modulated by mixing with underlying and overlying hydrostratigraphic units. Journal of Hydrology: Regional Studies 30:100713. https://doi.org/10.1016/j.ejrh.2020.100713 Tiwari K, Goyal R, Sarkar A (2017) GIS-based spatial distribution of groundwater quality and regional suitability evaluation for drinking water. Environmental Processes 4(3):645–662. https://doi.org/10.1007/s40710-017-0257-4 Tlili-Zrelli B, Gueddari M, Bouhlila R, Naceur M, O (2015) Groundwater hydrogeochemistry of Mateur alluvial aquifer (Northern Tunisia). J Hydrogeol Hydrol Eng 5(1). http://dx.doi.org/10.4172/2325-9647.1000128 Tundisi JG (2008) Water resources in the future: problems and solutions. Estudos avançados 22(63):7–16. https://doi.org/10.1590/S0103-40142008000200002 Vauclin M, Vieira SR, Vachaud G, Nielsen DR (1983) The Use of Cokriging with Limited Field Soil Observations 1. Soil Sci Soc Am J 47(2):175–184. 10.2136/sssaj1983.03615995004700020001x Warrick AW, Nielsen DR (1980) Spatial variability of soil physical properties in the field.. In: In: HILLEL D (ed) Applications of Soil Physics, vol 2. Academic, New York, pp 319–344 Yang Q, Li Z, Ma H, Wang L, Martín JD (2016) Identification of the hydrogeochemical processes and assessment of groundwater quality using classic integrated geochemical methods in the Southeastern part of Ordos basin, China. Environ Pollut 218:879–888. https://doi.org/10.1016/j.envpol.2016.08.017 Yang Q, Li Z, Ma H, Wang L, Martín JD (2016) Identification of the hydrogeochemical processes and assessment of groundwater quality using classic integrated geochemical methods in the Southeastern part of Ordos basin, China. Environ Pollut 218:879–888. https://doi.org/10.1016/j.envpol.2016.08.017 Yang Q, Li Z, Xie C, Liang J, Ma H (2020) Risk assessment of groundwater hydrochemistry for irrigation suitability in Ordos Basin, China. Nat Hazards 101:309–325. https://doi.org/10.1007/s11069-018-3451-4 Yin S, Xiao Y, Gu X, Hao Q, Liu H, Hao Z, Pei Q (2019) Geostatistical analysis of hydrochemical variations and nitrate pollution causes of groundwater in an alluvial fan plain. Acta Geophys 67(4):1191–1203. https://doi.org/10.1007/s11600-019-00302-5 Yin, S., Xiao, Y., Gu, X., Hao, Q., Liu, H., Hao, Z., ... & Pei, Q. (2019). Geostatistical analysis of hydrochemical variations and nitrate pollution causes of groundwater in an alluvial fan plain. Acta Geophysica, 67 (4), 1191-1203. https://doi.org/10.1007/s11600-019-00302-5 Cite Share Download PDF Status: Under Review Version 1 posted Reviews received at journal 09 Feb, 2022 Reviewers invited by journal 09 Feb, 2022 Editor assigned by journal 31 Dec, 2021 First submitted to journal 28 Dec, 2021 You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {\"props\":{\"pageProps\":{\"initialData\":{\"identity\":\"rs-1210837\",\"acceptedTermsAndConditions\":true,\"allowDirectSubmit\":false,\"archivedVersions\":[],\"articleType\":\"Research Article\",\"associatedPublications\":[],\"authors\":[{\"id\":82646666,\"identity\":\"ab1d7347-b6ee-40e7-9391-cff7f16c21bf\",\"order_by\":0,\"name\":\"José Francisco de Francisco Oliveira Júnior\",\"email\":\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA50lEQVRIiWNgGAWjYBACxgYGhgNAmgfEfgwWYmZuIFoLszEDgwGQYsSvBRmwSYO1MBDQwjwj9+CBj3sYZMylm49VF1T8ieZvB2r5UbENt8Nm5CUcnPGMgcdyzrG02zPOGOTOOMzYwNhz5jYeLTkGh3kOMPAY3Mgxu83bZpDbANTCzNhGlJb8b8UgLfNJ0JLDxgzSsoGglp43BgdngLTcOWYszXPGOHcjUMtBfH4xbM8x/vDhAIO9we3mh595KuRy550/fPDBjwo8WhrA1H8GBgkk0QM41QOBPJwlgUfVKBgFo2AUjGwAAIqIWMFOVRITAAAAAElFTkSuQmCC\",\"orcid\":\"https://orcid.org/0000-0002-6131-7605\",\"institution\":\"Federal University of Alagoas\",\"correspondingAuthor\":true,\"submittingAuthor\":false,\"prefix\":\"\",\"firstName\":\"José\",\"middleName\":\"Francisco de Francisco Oliveira\",\"lastName\":\"Júnior\",\"suffix\":\"\"},{\"id\":82646667,\"identity\":\"31bf35fb-9857-48fd-85d5-05e1db0b4cf7\",\"order_by\":1,\"name\":\"Alex Souza Moraes\",\"email\":\"\",\"orcid\":\"\",\"institution\":\"Universidade Federal Rural de Pernambuco\",\"correspondingAuthor\":false,\"submittingAuthor\":false,\"prefix\":\"\",\"firstName\":\"Alex\",\"middleName\":\"Souza\",\"lastName\":\"Moraes\",\"suffix\":\"\"},{\"id\":82646668,\"identity\":\"fc67ce93-1b9b-4ba9-862b-7a07e9351322\",\"order_by\":2,\"name\":\"Héliton Pandorfi\",\"email\":\"\",\"orcid\":\"\",\"institution\":\"Universidade Federal Rural de Pernambuco\",\"correspondingAuthor\":false,\"submittingAuthor\":false,\"prefix\":\"\",\"firstName\":\"Héliton\",\"middleName\":\"\",\"lastName\":\"Pandorfi\",\"suffix\":\"\"},{\"id\":82646669,\"identity\":\"b5c2e38b-49c8-4a9f-8f47-7914ec22af09\",\"order_by\":3,\"name\":\"Jhon Lennon Bezerra da Silva\",\"email\":\"\",\"orcid\":\"\",\"institution\":\"Universidade Federal Rural de Pernambuco\",\"correspondingAuthor\":false,\"submittingAuthor\":false,\"prefix\":\"\",\"firstName\":\"Jhon\",\"middleName\":\"Lennon Bezerra da\",\"lastName\":\"Silva\",\"suffix\":\"\"},{\"id\":82646670,\"identity\":\"6cb91c49-a21e-4f05-a1a8-85bd40a0e915\",\"order_by\":4,\"name\":\"Pedro Henrique Dias Batista\",\"email\":\"\",\"orcid\":\"\",\"institution\":\"Universidade Federal Rural de Pernambuco\",\"correspondingAuthor\":false,\"submittingAuthor\":false,\"prefix\":\"\",\"firstName\":\"Pedro\",\"middleName\":\"Henrique Dias\",\"lastName\":\"Batista\",\"suffix\":\"\"},{\"id\":82646671,\"identity\":\"bdf4afe0-5c0c-441d-af2e-69fd6ae47d7b\",\"order_by\":5,\"name\":\"Alexandre Maniçoba Rosa Ferraz Jardim\",\"email\":\"\",\"orcid\":\"\",\"institution\":\"Universidade Federal Rural de Pernambuco\",\"correspondingAuthor\":false,\"submittingAuthor\":false,\"prefix\":\"\",\"firstName\":\"Alexandre\",\"middleName\":\"Maniçoba Rosa Ferraz\",\"lastName\":\"Jardim\",\"suffix\":\"\"},{\"id\":82646672,\"identity\":\"076784c1-4f5b-49e3-bc12-cff6b547a928\",\"order_by\":6,\"name\":\"Gledson Luiz Pontes de Almeida\",\"email\":\"\",\"orcid\":\"\",\"institution\":\"Universidade Federal Rural de Pernambuco\",\"correspondingAuthor\":false,\"submittingAuthor\":false,\"prefix\":\"\",\"firstName\":\"Gledson\",\"middleName\":\"Luiz Pontes\",\"lastName\":\"de Almeida\",\"suffix\":\"\"},{\"id\":82646673,\"identity\":\"36bc9911-9f04-4257-9882-36111737b40b\",\"order_by\":7,\"name\":\"Taize Calvacante Santana\",\"email\":\"\",\"orcid\":\"\",\"institution\":\"Universidade Federal Rural de Pernambuco\",\"correspondingAuthor\":false,\"submittingAuthor\":false,\"prefix\":\"\",\"firstName\":\"Taize\",\"middleName\":\"Calvacante\",\"lastName\":\"Santana\",\"suffix\":\"\"},{\"id\":82646674,\"identity\":\"9d461dc1-beb3-482c-885e-404c06eb027f\",\"order_by\":8,\"name\":\"Marcio Mesquita\",\"email\":\"\",\"orcid\":\"\",\"institution\":\"Universidade Federal Rural de Pernambuco\",\"correspondingAuthor\":false,\"submittingAuthor\":false,\"prefix\":\"\",\"firstName\":\"Marcio\",\"middleName\":\"\",\"lastName\":\"Mesquita\",\"suffix\":\"\"}],\"badges\":[],\"createdAt\":\"2021-12-28 12:47:05\",\"currentVersionCode\":1,\"declarations\":\"\",\"doi\":\"10.21203/rs.3.rs-1210837/v1\",\"doiUrl\":\"https://doi.org/10.21203/rs.3.rs-1210837/v1\",\"draftVersion\":[],\"editorialEvents\":[],\"editorialNote\":\"\",\"failedWorkflow\":false,\"files\":[{\"id\":18138680,\"identity\":\"37d4c0ff-3ba0-4abb-84d1-a0afe8cef962\",\"added_by\":\"auto\",\"created_at\":\"2022-02-11 15:53:31\",\"extension\":\"png\",\"order_by\":1,\"title\":\"Figure 1\",\"display\":\"\",\"copyAsset\":false,\"role\":\"figure\",\"size\":221383,\"visible\":true,\"origin\":\"\",\"legend\":\"\\u003cp\\u003eSpatial map of location of the study area, municipalities of Recife (\\u003cstrong\\u003e1\\u003c/strong\\u003e), Olinda (\\u003cstrong\\u003e2\\u003c/strong\\u003e), Paulista (\\u003cstrong\\u003e3\\u003c/strong\\u003e), Camaragibe (\\u003cstrong\\u003e4\\u003c/strong\\u003e), São Lourenço da Mata (\\u003cstrong\\u003e5\\u003c/strong\\u003e), and Paudalho (\\u003cstrong\\u003e6\\u003c/strong\\u003e) and distribution of sampling points.\\u003c/p\\u003e\",\"description\":\"\",\"filename\":\"fig1.png\",\"url\":\"https://assets-eu.researchsquare.com/files/rs-1210837/v1/edb952532b8fbad33b18e5a5.png\"},{\"id\":18138567,\"identity\":\"348beb4e-aa03-4066-a57d-29d281f65a65\",\"added_by\":\"auto\",\"created_at\":\"2022-02-11 15:50:31\",\"extension\":\"png\",\"order_by\":2,\"title\":\"Figure 2\",\"display\":\"\",\"copyAsset\":false,\"role\":\"figure\",\"size\":321019,\"visible\":true,\"origin\":\"\",\"legend\":\"\\u003cp\\u003eSpatial altimetry (m) map of the study area, cities of Recife, Olinda, Paulista, Camaragibe, São Lourenço da Mata and Paudalho.\\u003c/p\\u003e\",\"description\":\"\",\"filename\":\"fig2.png\",\"url\":\"https://assets-eu.researchsquare.com/files/rs-1210837/v1/0e6bdf90e1ec7a716b28ce72.png\"},{\"id\":18138681,\"identity\":\"63c23758-42d2-4477-8119-8df8c0511b82\",\"added_by\":\"auto\",\"created_at\":\"2022-02-11 15:53:31\",\"extension\":\"png\",\"order_by\":3,\"title\":\"Figure 3\",\"display\":\"\",\"copyAsset\":false,\"role\":\"figure\",\"size\":245374,\"visible\":true,\"origin\":\"\",\"legend\":\"\\u003cp\\u003eSpatial map of the hydrogeological characterization of the groundwater of the Beberibe aquifer in the study area, cities of Recife, Olinda, Paulista, Camaragibe, São Lourenço da Mata, and Paudalho, and distribution of the sampling points.\\u003c/p\\u003e\",\"description\":\"\",\"filename\":\"fig3.png\",\"url\":\"https://assets-eu.researchsquare.com/files/rs-1210837/v1/df9ccc60d2bf260e6c7154f1.png\"},{\"id\":18138796,\"identity\":\"abf49882-3978-40db-834c-cf267f146422\",\"added_by\":\"auto\",\"created_at\":\"2022-02-11 15:56:31\",\"extension\":\"png\",\"order_by\":4,\"title\":\"Figure 4\",\"display\":\"\",\"copyAsset\":false,\"role\":\"figure\",\"size\":377583,\"visible\":true,\"origin\":\"\",\"legend\":\"\\u003cp\\u003eHydrochemical spatial map of the groundwater of the Beberibe aquifer in the study area, municipalities of Recife, Olinda, Paulista, Camaragibe, São Lourenço da Mata, and Paudalho, and distribution of the sampling points.\\u003c/p\\u003e\",\"description\":\"\",\"filename\":\"fig4.png\",\"url\":\"https://assets-eu.researchsquare.com/files/rs-1210837/v1/aa820aed50f7f668d3a2391d.png\"},{\"id\":18138570,\"identity\":\"a0af0902-3638-41bc-9756-ef998d02b15c\",\"added_by\":\"auto\",\"created_at\":\"2022-02-11 15:50:31\",\"extension\":\"png\",\"order_by\":5,\"title\":\"Figure 5\",\"display\":\"\",\"copyAsset\":false,\"role\":\"figure\",\"size\":108347,\"visible\":true,\"origin\":\"\",\"legend\":\"\\u003cp\\u003ePiper's diagram of the hydrochemical parameters of the mineral water of the Beberibe aquifer in the study area.\\u003c/p\\u003e\",\"description\":\"\",\"filename\":\"fig5.png\",\"url\":\"https://assets-eu.researchsquare.com/files/rs-1210837/v1/88e9585d295584b8c4ed783a.png\"},{\"id\":18138683,\"identity\":\"053a73a5-f7db-47b9-ae49-3c69372c2c40\",\"added_by\":\"auto\",\"created_at\":\"2022-02-11 15:53:32\",\"extension\":\"png\",\"order_by\":6,\"title\":\"Figure 6\",\"display\":\"\",\"copyAsset\":false,\"role\":\"figure\",\"size\":1498042,\"visible\":true,\"origin\":\"\",\"legend\":\"\\u003cp\\u003eSpatial distribution of the hydrochemical parameters of the mineral water from the Beberibe aquifer in the study area.\\u003c/p\\u003e\\u003cp\\u003eNO\\u003csub\\u003e3\\u003c/sub\\u003e: Nitrate; DR: Evaporative resistance; pH: Hydrogen potential; SO\\u003csub\\u003e4\\u003c/sub\\u003e: Sulfate; Na: Sodium; Ca: Calcium; Mg: Magnesium; Cl: Chlorine; K: Potassium; Mn: Manganese; F: Fluorine; Zn: Zinc; Fe: Iron; EC: Electrical conductivity; Hn: Hardness; HCO\\u003csub\\u003e3\\u003c/sub\\u003e: Bicarbonate; Si: Silicon; Sr: Strontium; Ba: Barium.\\u003c/p\\u003e\",\"description\":\"\",\"filename\":\"fig6.png\",\"url\":\"https://assets-eu.researchsquare.com/files/rs-1210837/v1/5f73fea702f1e60aed800e63.png\"},{\"id\":18138568,\"identity\":\"e0e049af-7468-48e5-af61-b131c11f4ef4\",\"added_by\":\"auto\",\"created_at\":\"2022-02-11 15:50:31\",\"extension\":\"png\",\"order_by\":7,\"title\":\"Figure 7\",\"display\":\"\",\"copyAsset\":false,\"role\":\"figure\",\"size\":187001,\"visible\":true,\"origin\":\"\",\"legend\":\"\\u003cp\\u003eBiplot plot of the principal components (A) and Pearson's correlation (B) of the hydrochemical parameters of the mineral water of the Beberibe aquifer in the study area.\\u003c/p\\u003e\",\"description\":\"\",\"filename\":\"fig7.png\",\"url\":\"https://assets-eu.researchsquare.com/files/rs-1210837/v1/7e1eb456068484c37cf70be8.png\"},{\"id\":18138573,\"identity\":\"81f7369b-15f6-48e3-9590-8370172bef1e\",\"added_by\":\"auto\",\"created_at\":\"2022-02-11 15:50:32\",\"extension\":\"png\",\"order_by\":8,\"title\":\"Figure 8\",\"display\":\"\",\"copyAsset\":false,\"role\":\"figure\",\"size\":72936,\"visible\":true,\"origin\":\"\",\"legend\":\"\\u003cp\\u003eLinear regression models for the variables with the highest correlation with nitrate (NO\\u003csub\\u003e3\\u003c/sub\\u003e); NO\\u003csub\\u003e3\\u003c/sub\\u003e x electrical conductivity (EC) (A); NO\\u003csub\\u003e3\\u003c/sub\\u003e x hardness (Hn) (B); NO\\u003csub\\u003e3\\u003c/sub\\u003e x evaporative resistance (DR) (C); NO\\u003csub\\u003e3\\u003c/sub\\u003e x magnesium (Mg) (D); NO\\u003csub\\u003e3\\u003c/sub\\u003e x sodium (Na) (E); NO\\u003csub\\u003e3\\u003c/sub\\u003e x chlorine (Cl) (F).\\u003c/p\\u003e\",\"description\":\"\",\"filename\":\"fig8.png\",\"url\":\"https://assets-eu.researchsquare.com/files/rs-1210837/v1/aa5c41e2adc7d19162b60b90.png\"},{\"id\":18138797,\"identity\":\"758add21-1899-4067-b4e5-82c02b89049f\",\"added_by\":\"auto\",\"created_at\":\"2022-02-11 15:56:36\",\"extension\":\"pdf\",\"order_by\":0,\"title\":\"\",\"display\":\"\",\"copyAsset\":false,\"role\":\"manuscript-pdf\",\"size\":1891412,\"visible\":true,\"origin\":\"\",\"legend\":\"\",\"description\":\"\",\"filename\":\"manuscript.pdf\",\"url\":\"https://assets-eu.researchsquare.com/files/rs-1210837/v1/0f43d2fd-3af5-41f3-89e4-21a67e0e34e0.pdf\"}],\"financialInterests\":\"\",\"formattedTitle\":\"\\u003cp\\u003eStatistical modeling of groundwater geochemistry in northeastern Brazil\\u003c/p\\u003e\",\"fulltext\":[{\"header\":\"1. Introduction\",\"content\":\"\\u003cp\\u003eBrazil has the largest freshwater reservoir in the world, equivalent to 14% of the total planet Earth, being a country of high relevance worldwide in the supply of mineral water (Tundisi, 2008; Kunrath et al., 2020). In recent years, bottled mineral water has grown significantly on a national scale and thus has become an attractive and profitable business in Brazil (DNPM, 2016; ABINAM, \\u003cspan citationid=\\\"CR1\\\" class=\\\"CitationRef\\\"\\u003e2020\\u003c/span\\u003e). The worldwide consumption of mineral water reached 284.1 billion liters in 2016, while other beverages, for example, soft drinks reached 159.1 billion liters in the same period (Huang and Liu, 2017), due to these increases, the need for evaluations and analyzes of the quality of drinking water becomes essential and with high frequency (Shruti et al., 2020).\\u003c/p\\u003e \\u003cp\\u003eThe consumer bases his choice on extrinsic factors, such as brand and personal habits (Pacheco et al., 2018). Mineral waters are products differentiated by intrinsic factors, mainly by their aroma and flavor, which are classified according to the hydrochemical parameters that, in turn, result from the characteristics of the rocks and the aquifer (depth of the water table, as for circulation, confinement, and porosity) that the water flows until its capture and packaging (Dias and Bernardes, 2016; Guadayol et al., 2016; Pacheco et al., 2018).\\u003c/p\\u003e \\u003cp\\u003eThe quality of mineral water is highlighted by Nunes et al. (2018), who showed the importance of its determination, in which, the hydrochemical composition is the qualitative aspect of greatest importance to be investigated since it has a direct impact on human health. Bulia and Enzweiler (\\u003cspan citationid=\\\"CR11\\\" class=\\\"CitationRef\\\"\\u003e2018\\u003c/span\\u003e) evaluated the bottled water taken from 78 different sources in the state of S\\u0026atilde;o Paulo, Brazil, and found that the hydrochemical composition, in general, satisfactorily represented the groundwater in the aquifers corresponding to the bottled mineral water samples. Schwanz et al. (2016), highlighted the potential of concentrations of perfluorooctanoic acid (PFOA) classified in Group 2B (i.e., the possible presence of a carcinogen for humans), of bottled mineral water in Brazil.\\u003c/p\\u003e \\u003cp\\u003eThe statistical representation of the hydrochemical composition of mineral water is relatively representative. However, the complementary analysis of the spatial variability allows a careful assessment when investigating the distribution of mineral water in the regions and, thus, makes it possible to question the intrinsic and extrinsic factors that alter the hydrochemical composition of these waters.\\u003c/p\\u003e \\u003cp\\u003eStudies like Gilabert-Alarc\\u0026oacute;n et al. (2018), who evaluated the processes that affect the quality of groundwater, using hydrochemical tracers, as the main ionic composition and stable isotopes of water and nitrates, through the synthetic pollution index, allowed the spatial regionalization of groundwater quality. Concomitant with the use of geostatistical techniques that present satisfactory recent results in the spatialization of the hydrochemical composition of groundwater (Jang, 2015; Gu et al., 2017; Arshad et al., \\u003cspan citationid=\\\"CR5\\\" class=\\\"CitationRef\\\"\\u003e2019\\u003c/span\\u003e; Bodrud-Doza et al., \\u003cspan citationid=\\\"CR9\\\" class=\\\"CitationRef\\\"\\u003e2019\\u003c/span\\u003e).\\u003c/p\\u003e \\u003cp\\u003eBulia and Enzweiler (\\u003cspan citationid=\\\"CR11\\\" class=\\\"CitationRef\\\"\\u003e2018\\u003c/span\\u003e), highlighted that the groundwater sampling method based on bottled mineral water is a low-cost and viable alternative for determining the composition and hydrochemical quality of the water, without logistical and bureaucratic limitations with conventional sampling (laboratory) on a continental scale.\\u003c/p\\u003e \\u003cp\\u003eBecause of the above, the study aims to characterize the hydrochemical and qualitative composition of the groundwater of the Beberibe aquifer located in Northeastern Brazil via geostatistical, geoprocessing, and multivariate analysis techniques.\\u003c/p\\u003e\"},{\"header\":\"2. Material And Methods\",\"content\":\"\\u003cdiv class=\\\"Section2\\\" id=\\\"Sec3\\\"\\u003e\\n \\u003ch2\\u003e2.1. Location and characterization of the study area\\u003c/h2\\u003e\\n \\u003cp\\u003eThe study was realized in the municipalities of Recife, Olinda, Paulista, S\\u0026atilde;o Louren\\u0026ccedil;o da Mata, Camaragibe, and Paudalho, both municipalities are located in the state of Pernambuco (PE), Northeastern Brazil (NEB) (Figure \\u003cspan class=\\\"InternalRef\\\"\\u003e1\\u003c/span\\u003e). The municipality of Paudalho is located in the mesoregion of Zona da Mata, while Recife, Olinda, Paulista, S\\u0026atilde;o Louren\\u0026ccedil;o da Mata, and Camaragibe are located in the mesoregion Metropolitana de Recife, the municipalities are located between the parallels 7\\u0026deg;51\\u0026apos;0.0\\u0026apos;\\u0026apos; S \\u0026ndash; 8\\u0026deg;8\\u0026apos;0.0\\u0026apos;\\u0026apos; S and between the meridians 35\\u0026deg;11\\u0026apos;0.0\\u0026apos;\\u0026apos; W \\u0026ndash; 34\\u0026deg;54\\u0026apos;0.0\\u0026apos;\\u0026apos; W. According to the K\\u0026ouml;ppen-Geiger climate classification, the study region is in two climatic zones: \\u0026ldquo;Am\\u0026rdquo;―humid or sub-humid tropical climate and \\u0026ldquo;Aw\\u0026rdquo;―tropical climate with a dry winter season (Figure \\u003cspan class=\\\"InternalRef\\\"\\u003e1\\u003c/span\\u003e) - (Beck et al., \\u003cspan class=\\\"CitationRef\\\"\\u003e2018\\u003c/span\\u003e).\\u003c/p\\u003e\\n \\u003cp\\u003e\\u003cstrong\\u003e[INSERT TO\\u003c/strong\\u003e FIGURE \\u003cspan class=\\\"InternalRef\\\"\\u003e1\\u003c/span\\u003e \\u003cstrong\\u003eHERE]\\u003c/strong\\u003e\\u003c/p\\u003e\\n \\u003cp\\u003eThe research was conducted in the Beberibe aquifer, the main aquifer that supplies the northern part of the Metropolitan Region of Recife (RMRE), having been explored for over fifty years (CPRH, 2020). In the last decade, with the significant increase in the number of artesian wells, as well as the increase in industrial activities implemented, followed by population densification and the urbanization process, among others, there was a need to assess the quality of groundwater (CPRH, 2020).\\u003c/p\\u003e\\n \\u003cp\\u003eRMRE is located on a Deltaic Plain, which has distinct hydrogeological characteristics, composed of sediments of different origins: fluvial, marine, colluvial, and mangroves, which cover the coastal sedimentary basins of Pernambuco and Para\\u0026iacute;ba, respectively to the south and north of the region, separated by the structural divider Lineamento Pernambuco. Figure \\u003cspan class=\\\"InternalRef\\\"\\u003e2\\u003c/span\\u003e illustrates the spatial distribution of hypsometry, to establish a detailed understanding of the Beberibe aquifer that contemplates the western portion of the Pernambuco coast, with a variation from 0 (coast of the state of Pernambuco) to 261 m of altitude.\\u003c/p\\u003e\\n \\u003cp\\u003e\\u003cstrong\\u003e[INSERT TO\\u003c/strong\\u003e FIGURE \\u003cspan class=\\\"InternalRef\\\"\\u003e2\\u003c/span\\u003e \\u003cstrong\\u003eHERE]\\u003c/strong\\u003e\\u003c/p\\u003e\\n\\u003c/div\\u003e\\n\\u003cdiv class=\\\"Section2\\\" id=\\\"Sec4\\\"\\u003e\\n \\u003ch2\\u003e2.2. Data collect\\u003c/h2\\u003e\\n \\u003cp\\u003eThe in-situ data collection was carried out in 34 extraction units (i.e., artesian wells) that capture mineral water, georeferenced, based on the analyzes made available by the companies and following the rules of the Pernambuco Health Surveillance Agency (APEVISA, 2020). The representation of the spatial distribution of artesian wells is shown in Figure \\u003cspan class=\\\"InternalRef\\\"\\u003e1\\u003c/span\\u003e.\\u003c/p\\u003e\\n \\u003cp\\u003eThe record of the hydrochemical values were those made available on the label of each mineral water unit, according to electrical conductivity (EC, \\u0026micro;S. Cm\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e), hardness (Hn, mg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e CaCO\\u003csub\\u003e3\\u003c/sub\\u003e), hydrogen potential (pH, dimensionless), evaporative resistance (DR, mg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e), calcium (Ca, mg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e), magnesium (Mg, mg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e), sodium (Na, mg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e), potassium (K, mg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e), bicarbonate (HCO\\u003csub\\u003e3\\u003c/sub\\u003e, mg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e), sulfate (SO\\u003csub\\u003e4\\u003c/sub\\u003e, mg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e), chlorine (Cl, mg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e), iron (Fe, mg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e), manganese (Mn, mg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e), fluorine (F, mg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e), silicon (Si, mg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e), nitrate (NO\\u003csub\\u003e3\\u003c/sub\\u003e, mg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e), strontium (Sr, mg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e), zinc (Zn, mg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e) and barium (Ba, mg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e).\\u003c/p\\u003e\\n \\u003cp\\u003eTable \\u003cspan class=\\\"InternalRef\\\"\\u003e1\\u003c/span\\u003e shows the maximum permitted values (MPVs) of the hydrochemical parameters of mineral water for human consumption, according to the Ministry of Health, Ordinance n\\u0026deg;. 26, of January 15, 1999, and, of the Environmental Company of the State of S\\u0026atilde;o Paulo (CETESB) \\u0026ndash; (Minist\\u0026eacute;rio da Sa\\u0026uacute;de, 1999; Minist\\u0026eacute;rio da Sa\\u0026uacute;de, 2006; CETESB, 2016).\\u003c/p\\u003e\\n \\u003cdiv class=\\\"gridtable\\\"\\u003e\\u003ctable border=\\\"1\\\" id=\\\"Tab1\\\"\\u003e\\n \\u003ccaption language=\\\"En\\\"\\u003e\\n \\u003cdiv class=\\\"CaptionNumber\\\"\\u003eTable 1\\u003c/div\\u003e\\n \\u003cdiv class=\\\"CaptionContent\\\"\\u003e\\n \\u003cp\\u003eHydrochemical parameters of water for human consumption and their maximum allowed values (MPVs) according to Brazilian legislation.\\u003c/p\\u003e\\n \\u003c/div\\u003e\\n \\u003c/caption\\u003e\\n \\u003ccolgroup cols=\\\"6\\\"\\u003e\\u003c/colgroup\\u003e\\n \\u003cthead\\u003e\\n \\u003ctr\\u003e\\n \\u003cth align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eParameters\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eUnit\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eMPV\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eParameters\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eUnit\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eMPV\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003c/tr\\u003e\\n \\u003c/thead\\u003e\\n \\u003ctbody\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e1\\u003c/sup\\u003eNO\\u003csub\\u003e3\\u003c/sub\\u003e\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003emg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e10\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e11\\u003c/sup\\u003eF\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003emg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e1.5\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e2\\u003c/sup\\u003eDR\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003emg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e1000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e12\\u003c/sup\\u003eZn\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003emg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e5\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e3\\u003c/sup\\u003epH\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e-\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e6 to 9\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e13\\u003c/sup\\u003eFe\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003emg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e0.3\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e4\\u003c/sup\\u003eSO\\u003csub\\u003e4\\u003c/sub\\u003e\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003emg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e400\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e14\\u003c/sup\\u003eEC\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u0026micro;S/cm\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e350\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e5\\u003c/sup\\u003eNa\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003emg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e200\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e15\\u003c/sup\\u003eHn\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003emg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e CaCO\\u003csub\\u003e3\\u003c/sub\\u003e\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e500\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e6\\u003c/sup\\u003eCa\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003emg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e250\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e16\\u003c/sup\\u003eHCO\\u003csub\\u003e3\\u003c/sub\\u003e\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e-\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e-\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e7\\u003c/sup\\u003eMg\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003emg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e100\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e17\\u003c/sup\\u003eSi\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e-\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e-\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e8\\u003c/sup\\u003eCl\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003emg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e2\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e18\\u003c/sup\\u003eSr\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e-\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e-\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e9\\u003c/sup\\u003eK\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003emg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e10\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e19\\u003c/sup\\u003eBa\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003emg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e1\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e10\\u003c/sup\\u003eMn\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003emg. L\\u003csup\\u003e\\u0026minus;1\\u003c/sup\\u003e\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e0.1\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e-\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e-\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e-\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003c/tbody\\u003e\\n \\u003ctfoot\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd colspan=\\\"6\\\"\\u003e\\u003csup\\u003e1\\u003c/sup\\u003eNO\\u003csub\\u003e3\\u003c/sub\\u003e: Nitrate; \\u003csup\\u003e2\\u003c/sup\\u003eDR: Evaporative resistance; \\u003csup\\u003e3\\u003c/sup\\u003epH: Hydrogen potential; \\u003csup\\u003e4\\u003c/sup\\u003eSO\\u003csub\\u003e4\\u003c/sub\\u003e: Sulfate; \\u003csup\\u003e5\\u003c/sup\\u003eNa: Sodium; \\u003csup\\u003e6\\u003c/sup\\u003eCa: Calcium; \\u003csup\\u003e7\\u003c/sup\\u003eMg: Magnesium; \\u003csup\\u003e8\\u003c/sup\\u003eCl: Chlorine; \\u003csup\\u003e9\\u003c/sup\\u003eK: Potassium; \\u003csup\\u003e10\\u003c/sup\\u003eMn: Manganese; \\u003csup\\u003e11\\u003c/sup\\u003eF: Fluorine; \\u003csup\\u003e12\\u003c/sup\\u003eZn: Zinc; \\u003csup\\u003e13\\u003c/sup\\u003eFe: Iron; \\u003csup\\u003e14\\u003c/sup\\u003eEC: Electrical conductivity; \\u003csup\\u003e15\\u003c/sup\\u003eHn: Hardness; \\u003csup\\u003e16\\u003c/sup\\u003eHCO\\u003csub\\u003e3\\u003c/sub\\u003e: Bicarbonate; \\u003csup\\u003e17\\u003c/sup\\u003eSi: Silicon; \\u003csup\\u003e18\\u003c/sup\\u003eSr: Strontium; \\u003csup\\u003e19\\u003c/sup\\u003eBa: Barium.\\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003c/tfoot\\u003e\\n \\u003c/table\\u003e\\n \\u003c/div\\u003e\\n\\u003c/div\\u003e\\n\\u003cdiv class=\\\"Section2\\\" id=\\\"Sec5\\\"\\u003e\\n \\u003ch2\\u003e2.3. Piper\\u0026apos;s diagram\\u003c/h2\\u003e\\n \\u003cp\\u003eThe Piper diagram was constructed and applied in the classification and comparison of distinct groups of water, associated with dominant cations and anions, also called tri-linear diagrams, in this type of representation, the concentrations of each ion analyzed were given in percentages, as well, plotting the proportions of the main cations (Ca\\u003csup\\u003e+2\\u003c/sup\\u003e, Mg\\u003csup\\u003e+2\\u003c/sup\\u003e, Na\\u003csup\\u003e+\\u003c/sup\\u003e and K\\u003csup\\u003e+\\u003c/sup\\u003e) and the main anions (HCO\\u003csup\\u003e\\u0026minus;3\\u003c/sup\\u003e, Cl\\u003csup\\u003e\\u0026minus;\\u003c/sup\\u003e and SO\\u003csup\\u003e\\u0026minus;4\\u003c/sup\\u003e) in two respective triangular diagrams, combining the information from the two triangles in a rhombus located between them, with the main types of chemical constituents of water (Piper, 1944; Anders et al., \\u003cspan class=\\\"CitationRef\\\"\\u003e2014\\u003c/span\\u003e; Yang et al., 2016).\\u003c/p\\u003e\\n \\u003cp\\u003eThe proportions were plotted in triangular graphical representations and their scales, for the proportion of the variables, corresponding to 100%. The representations show the relative proportions of the main ions, but not their absolute concentrations. The procedure for plotting a point in Piper\\u0026apos;s diagram consisted of entering the percentage of certain cations and anions, in the triangular bases, and a line was drawn corresponding to the value of this percentage. The crossing of the lines drawn on the triangular bases for the cations and anions determined the point that is transcribed in the central diamond, where the hydrochemical classification of the water is. The software used for processing Piper\\u0026apos;s diagram was Grapher version 13.\\u003c/p\\u003e\\n\\u003c/div\\u003e\\n\\u003cdiv class=\\\"Section2\\\" id=\\\"Sec6\\\"\\u003e\\n \\u003ch2\\u003e2.4. Hydrogeological and hydrochemical map\\u003c/h2\\u003e\\n \\u003cp\\u003eBased on the database provided by the Brazilian Institute of Geography and Statistics (IBGE), the hydrogeological (IBGE, 2018) and hydrochemical (IBGE, 2020) map of the study region was generated.\\u003c/p\\u003e\\n \\u003cp\\u003eThe hydrogeological map of NEB cartographically represents the productivity of aquifers in this region, based on the values of flows and specific flows of 54,864 tubular wells. Such information is stored in a database, with the procedures used in the Map of Hydrogeological Domains/Subdomains of Brazil and the Hydrogeological Charters of Brazil to the Millionth (CPRM, 2007).\\u003c/p\\u003e\\n \\u003cp\\u003eThe hydrochemical map of NEB\\u0026apos;s underground water sources brings together a collection of 10,478 hydrochemical analyzes, all from tubular wells and chemically homogeneous domains concerning potability, chemical facies, and the suitability of water for use in irrigation. The chemical reports were incorporated into a database (elaborated in Access), where they were classified in terms of potability, chemical facies, and the suitability of waters for use in irrigation, and represented in a shapefile format.\\u003c/p\\u003e\\n \\u003cp\\u003eThese determinations were migrated and georeferenced in GeoMedia, where work was carried out to individualize chemically homogeneous zones, using geological, physiographic, and hydrogeological criteria, which allowed the demarcation of units that have similar characteristics within their limits.\\u003c/p\\u003e\\n \\u003cp\\u003eThe data made available on the platform are in the shapefile format, with an established pre-classification, in which reclassification techniques by geoprocessing were applied to reclassify these shapefiles, for the area comprising the study. The shapefiles were processed using ArcGIS\\u0026reg; software, version 10.2.2 (ESRI, Redlands, California, USA).\\u003c/p\\u003e\\n\\u003c/div\\u003e\\n\\u003cdiv class=\\\"Section2\\\" id=\\\"Sec7\\\"\\u003e\\n \\u003ch2\\u003e2.5. Statistical analysis\\u003c/h2\\u003e\\n \\u003cp\\u003eThe data obtained were subjected to analysis of descriptive statistics based on the mean, median, and coefficient of variation (CV, %), being classified as low when the CV \\u0026lt; 12%; medium when 12% \\u0026lt; CV \\u0026lt; 24%, and high when CV \\u0026gt; 24% (Warrick and Nielsen, 1980). Subsequently, the Kolmogorov-Smirnov normality test (KS test) was applied with a significance level set at p \\u0026le; 0.01.\\u003c/p\\u003e\\n \\u003cp\\u003eThe geostatistics analysis was performed based on the calculation of the classic semivariances (Eq.\\u0026nbsp;1), which estimated the structure and the spatial dependence between the pairs of observations.\\u003c/p\\u003e\\n \\u003cdiv class=\\\"Equation\\\" id=\\\"Equa\\\"\\u003e\\n \\u003cdiv class=\\\"mathdisplay\\\" id=\\\"FileID_Equa\\\" name=\\\"EquationSource\\\"\\u003e$$\\\\text{\\u0026gamma;}\\\\text{(h) = }\\\\frac{\\\\text{1}}{\\\\text{2N(h)}}\\\\text{ }{\\\\sum }_{\\\\text{i}\\\\text{ }\\\\text{=}\\\\text{ }\\\\text{1}}^{\\\\text{N(h)}}{\\\\left[\\\\text{Z}\\\\left({\\\\text{X}}_{\\\\text{i}}\\\\right)\\\\text{-}\\\\text{Z(}{\\\\text{X}}_{\\\\text{i}}\\\\text{ + h)}\\\\right]}^{\\\\text{2}}\\\\text{ }\\\\text{(}\\\\text{1}\\\\text{)}$$\\u003c/div\\u003e\\n \\u003c/div\\u003e\\n \\u003cp\\u003ewhere,\\u003c/p\\u003e\\n \\u003cp\\u003e\\u0026gamma;(h) - is the experimental semi-variance estimator, obtained from the sampled values Z(X\\u003csub\\u003ei\\u003c/sub\\u003e), Z(X\\u003csub\\u003ei\\u003c/sub\\u003e+h);\\u003c/p\\u003e\\n \\u003cp\\u003eN(h) - is the number of measured value pairs separated by the vector or delay distance;\\u003c/p\\u003e\\n \\u003cp\\u003eh - is the distance between sample pairs (i.e., it is the distance between two samples);\\u003c/p\\u003e\\n \\u003cp\\u003eZ(X\\u003csub\\u003ei\\u003c/sub\\u003e) and Z(X\\u003csub\\u003ei\\u003c/sub\\u003e + h) - are the values of the \\u003cem\\u003ei-th\\u003c/em\\u003e observation of the regionalized variable, collected at points X\\u003csub\\u003ei\\u003c/sub\\u003e and X\\u003csub\\u003ei\\u003c/sub\\u003e+h (\\u003cem\\u003ei\\u003c/em\\u003e = 1, ..., n), separated by the vector h.\\u003c/p\\u003e\\n \\u003cp\\u003eSpatial dependence was analyzed by adjusting the semivariogram based on the estimate of the semivariance using the GEO-EAS\\u0026reg; program (England et al., 1989). The data were adjusted to the spherical, exponential, and Gaussian models (Eqs. 2, \\u003cspan class=\\\"InternalRef\\\"\\u003e3\\u003c/span\\u003e, and 4), respectively, according to Deutsch et al. (1998). The spherical, exponential, and Gaussian models are referred to in the literature as transitive theoretical models and are more common for adjusting semivariograms (Gois et al., 2015).\\u003c/p\\u003e\\n \\u003cp\\u003eSpherical Model:\\u003c/p\\u003e\\n \\u003cdiv class=\\\"Equation\\\" id=\\\"Equb\\\"\\u003e\\n \\u003cdiv class=\\\"mathdisplay\\\" id=\\\"FileID_Equb\\\" name=\\\"EquationSource\\\"\\u003e$$\\\\text{\\u0026gamma;}\\\\left(\\\\text{h}\\\\right)\\\\text{ =}\\\\left\\\\{\\\\begin{array}{c}\\\\begin{array}{c}\\\\\\\\ {\\\\text{C}}_{\\\\text{0}}\\\\text{ }\\\\text{+}\\\\text{ }\\\\text{C}\\\\text{ ·}\\\\end{array}\\\\\\\\ \\\\begin{array}{c}\\\\\\\\ {\\\\text{C}}_{\\\\text{0}}\\\\text{ }\\\\text{+}\\\\text{ }\\\\text{C}\\\\text{,}\\\\end{array}\\\\end{array}\\\\begin{array}{c}\\\\begin{array}{c}\\\\left[\\\\text{1.5}\\\\frac{\\\\text{h}}{\\\\text{a}}\\\\text{- 0.5}{\\\\left(\\\\frac{\\\\text{h}}{\\\\text{a}}\\\\right)}^{\\\\text{3}}\\\\right]\\\\text{,}\\\\end{array}\\\\\\\\ \\\\begin{array}{c}\\\\end{array}\\\\end{array}\\\\begin{array}{c}\\\\begin{array}{c}\\\\\\\\ \\\\text{ for}\\\\text{ }\\\\text{0 \\u0026le; h \\u0026le; }\\\\text{a}\\\\end{array}\\\\\\\\ \\\\begin{array}{c}\\\\\\\\ \\\\text{ }\\\\text{ }\\\\text{for}\\\\text{ }\\\\text{ }\\\\text{ }\\\\text{h \\u0026gt; }\\\\text{a}\\\\end{array}\\\\end{array}\\\\right.\\\\text{ (2)}$$\\u003c/div\\u003e\\n \\u003c/div\\u003e\\n \\u003cp\\u003eExponential Model:\\u003c/p\\u003e\\n \\u003cdiv class=\\\"Equation\\\" id=\\\"Equ1\\\"\\u003e\\n \\u003cdiv class=\\\"mathdisplay\\\" id=\\\"FileID_Equ1\\\" name=\\\"EquationSource\\\"\\u003e$$\\\\text{\\u0026gamma;}\\\\left(\\\\text{h}\\\\right)\\\\text{ = }{\\\\text{C}}_{\\\\text{0}}\\\\text{ }\\\\text{+}\\\\text{ }\\\\text{C . }\\u0026lfloor;\\\\text{1-}\\\\text{exp}\\\\left(\\\\text{-}\\\\frac{\\\\text{3h}}{\\\\text{a}}\\\\right)\\u0026rfloor;\\\\text{ }$$\\u003c/div\\u003e\\n \\u003cdiv class=\\\"EquationNumber\\\"\\u003e3\\u003c/div\\u003e\\n \\u003c/div\\u003e\\n \\u003cp\\u003eGaussian model:\\u003c/p\\u003e\\n \\u003cdiv class=\\\"Equation\\\" id=\\\"Equc\\\"\\u003e\\n \\u003cdiv class=\\\"mathdisplay\\\" id=\\\"FileID_Equc\\\" name=\\\"EquationSource\\\"\\u003e$$\\\\text{\\u0026gamma;}\\\\left(\\\\text{h}\\\\right)\\\\text{ = }{\\\\text{C}}_{\\\\text{0}}\\\\text{ }\\\\text{+}\\\\text{ }\\\\text{C }\\\\text{·}\\\\left[\\\\text{1-}\\\\text{exp}\\\\left(\\\\text{-}\\\\frac{{\\\\text{(3h)}}^{\\\\text{2}}}{{\\\\text{a}}^{\\\\text{2}}}\\\\right)\\\\right]\\\\text{ (4)}$$\\u003c/div\\u003e\\n \\u003c/div\\u003e\\n \\u003cp\\u003ewhere,\\u003c/p\\u003e\\n \\u003cp\\u003e\\u0026gamma;(h) - is the experimental semi-variance estimator;\\u003c/p\\u003e\\n \\u003cp\\u003eC\\u003csub\\u003e0\\u003c/sub\\u003e+C - is the sill (i.e., it is the nugget effect plus the variance dispersion, given by the acronyms C\\u003csub\\u003e0\\u003c/sub\\u003e and C, respectively);\\u003c/p\\u003e\\n \\u003cp\\u003eh - is the distance between sample pairs;\\u003c/p\\u003e\\n \\u003cp\\u003e\\u003cem\\u003ea\\u003c/em\\u003e - is the range (m).\\u003c/p\\u003e\\n \\u003cp\\u003eThe best models of the adjusted semivariograms were validated by the cross-validation of the Jack-Knifing test, in which the mean must be close to zero and the standard deviation close to 1 (e.g., Vauclin et al., 1983), the program used for that analysis was GEO-EAS\\u0026reg; (England et al., 1989).\\u003c/p\\u003e\\n \\u003cp\\u003eThe degree of spatial dependence (DSD) was classified according to Cambardella et al. (\\u003cspan class=\\\"CitationRef\\\"\\u003e1994\\u003c/span\\u003e), which suggests strong dependence (St) (DSD \\u0026lt; 25%); moderate dependence (Md) (DSD between 25 and 75%), and weak dependence (We): (DSD\\u0026gt; 75%). For making the kriging maps, the Surfer 9 program (Golden Software, 2010) was used.\\u003c/p\\u003e\\n \\u003cp\\u003eFor principal component analysis (PCA), 19 variables were admitted. Based on the principal components (PC), the covariance matrix was obtained to extract the eigenvalues that originate the eigenvectors. To identify the variables that presented correlation, the Kaiser criterion was used, considering the eigenvalues greater than 1.0, which generate components with a relevant amount of information contained in the original data (Kaiser, 1958).\\u003c/p\\u003e\\n \\u003cp\\u003ePearson\\u0026apos;s correlation (r) was also performed for all variables, seeking to correlate with the PCA, to highlight the similarities between the variables. The program used for PCA and the correlation r was RStudio, version 3.6.1 (R Core Team, 2019).\\u003c/p\\u003e\\n \\u003cp\\u003eLinear regression models were established to justify and clarify the main hydrochemical characteristics that are correlated with nitrate (NO\\u003csub\\u003e3\\u003c/sub\\u003e), a component that has a direct impact on human health. It was verified among the main variables, the one with the highest coefficient of determination (R\\u003csup\\u003e2\\u003c/sup\\u003e \\u0026gt; 0.5) to reduce the main relevance and correlation with NO\\u003csub\\u003e3\\u003c/sub\\u003e. The analyzes were performed using the Minitab software, version 19.\\u003c/p\\u003e\\n\\u003c/div\\u003e\"},{\"header\":\"3. Results And Discussion\",\"content\":\"\\u003cdiv class=\\\"Section2\\\" id=\\\"Sec9\\\"\\u003e\\n \\u003ch2\\u003e3.1. Hydrogeological and hydrochemical characterization of the Beberibe aquifer\\u003c/h2\\u003e\\n \\u003cp\\u003eFigure \\u003cspan class=\\\"InternalRef\\\"\\u003e3\\u003c/span\\u003e shows spatially the hydrogeological characterization of the groundwater of the Beberibe aquifer. Most of the artesian wells in the mineral water production units are close to the coast, a densely populated region (Figures \\u003cspan class=\\\"InternalRef\\\"\\u003e1\\u003c/span\\u003e and \\u003cspan class=\\\"InternalRef\\\"\\u003e3\\u003c/span\\u003e). The artesian wells close to the coast are of porous geological formation, with the lowest altitudes (Figure \\u003cspan class=\\\"InternalRef\\\"\\u003e2\\u003c/span\\u003e), with flow variations from medium to very high. From the western portion of the study region, there is a fissural geological formation, considering the higher altitudes, with greater rocky extension and flow variations from medium to very high (Figures \\u003cspan class=\\\"InternalRef\\\"\\u003e2\\u003c/span\\u003e and \\u003cspan class=\\\"InternalRef\\\"\\u003e3\\u003c/span\\u003e).\\u003c/p\\u003e\\n \\u003cp\\u003e\\u003cstrong\\u003e[INSERT TO\\u003c/strong\\u003e FIGURE \\u003cspan class=\\\"InternalRef\\\"\\u003e3\\u003c/span\\u003e \\u003cstrong\\u003eHERE]\\u003c/strong\\u003e\\u003c/p\\u003e\\n \\u003cp\\u003eThere is a specific homogeneous flow occurring close to the coastal region, categorized as \\u0026ldquo;very high productivity\\u0026rdquo; (Figure \\u003cspan class=\\\"InternalRef\\\"\\u003e3\\u003c/span\\u003e). The obtained result corroborates with the studies carried out by Hu et al. (2017) who analyzed the hydrogeological configuration of groundwater for the entire Brazilian territory and highlighted good water retention capacity for recharging and discharging aquifers. The authors also pointed out that the aquifers close to the Atlantic Ocean play a role as a \\u0026ldquo;wall or barrier\\u0026rdquo;, because the water tables in these regions always remain at the same level as the sea surface at the limits of the coast, thus maintaining better retention, which will lead to homogenization of water flow.\\u003c/p\\u003e\\n \\u003cp\\u003eAlthough the original depth of the well is unknown, part of the water may come from shallower horizons (Pernambuco coast) and another part may come from the western region of the state, via deeper horizons, in superior topographic conditions. Such observations are in line with Silva et al. (2020), who studied the predominant hydrochemical mechanisms and processes that control the hydrochemistry of groundwater in a region of Acre, Northern Brazil, that regardless of the knowledge of the depth of the well, the origin of the water can come from different horizons, again highlighting the topographic conditions under the influence of a given region.\\u003c/p\\u003e\\n \\u003cp\\u003eFigure \\u003cspan class=\\\"InternalRef\\\"\\u003e4\\u003c/span\\u003e highlights the hydrochemical characterization of the groundwater in the Beberibe aquifer. It appears that in all of the sampled points (artesian wells), the waters are of the mixed bicarbonate type (Bm) and high potability for human consumption (IBGE, 2020). However, one should not only consider the hydrochemical reactions that occur in the Beberibe aquifer, but also the anthropic actions that occur in RMRE.\\u003c/p\\u003e\\n \\u003cp\\u003e\\u003cstrong\\u003e[INSERT TO\\u003c/strong\\u003e FIGURE \\u003cspan class=\\\"InternalRef\\\"\\u003e4\\u003c/span\\u003e \\u003cstrong\\u003eHERE]\\u003c/strong\\u003e\\u003c/p\\u003e\\n \\u003cp\\u003eThe composition of groundwater can provide relevant information about interactions in aquifers since hydrochemistry is the result of specific reactions related to rock-water interactions that vary within an aquifer unit (Teramoto et al., 2020). Consequently, the interpretation of aquifer mineralogy and water chemistry data may indicate the transport of groundwater within and between aquifers.\\u003c/p\\u003e\\n \\u003cp\\u003eAs recently highlighted by Teramoto et al. (2020), the hydrochemistry in the Guarani aquifer, due to its mixed bicarbonate composition, influences the geological formation of the porous type. This statement is confirmed in the present study since the geological composition is heterogeneous and presents an advanced state of decomposition, which implies a high variability of bicarbonates.\\u003c/p\\u003e\\n\\u003c/div\\u003e\\n\\u003cdiv class=\\\"Section2\\\" id=\\\"Sec10\\\"\\u003e\\n \\u003ch2\\u003e3.2. Piper\\u0026apos;s diagram\\u003c/h2\\u003e\\n \\u003cp\\u003eThe hydrochemical face is a term used to describe the percentage of the chemical composition of the water, resulting from the combined effect of the solution\\u0026apos;s kinetics, rock-water interactions, hydrogeological configurations, and sources of contamination. Piper\\u0026apos;s trilinear diagram (Piper, 1944) consists of a diamond and a pair of equilateral triangles, in which two triangles represent the anion and the cation, respectively, connected by a diamond-shaped diagram. Figure \\u003cspan class=\\\"InternalRef\\\"\\u003e5\\u003c/span\\u003e shows the hydrochemical composition of the groundwater in the Beberibe aquifer.\\u003c/p\\u003e\\n \\u003cp\\u003e\\u003cstrong\\u003e[INSERT TO\\u003c/strong\\u003e FIGURE \\u003cspan class=\\\"InternalRef\\\"\\u003e5\\u003c/span\\u003e \\u003cstrong\\u003eHERE]\\u003c/strong\\u003e\\u003c/p\\u003e\\n \\u003cp\\u003eThe chemical composition of groundwater is mainly of the Na-K-Cl-SO\\u003csub\\u003e4\\u003c/sub\\u003e type (Figure \\u003cspan class=\\\"InternalRef\\\"\\u003e5\\u003c/span\\u003e), where Na-K corresponds to approximately 90% in the first triangle. The hydrochemical type Ca-Mg-Cl-SO\\u003csub\\u003e4\\u003c/sub\\u003e also appeared in Piper\\u0026apos;s diagram, to which, Ca-Mg corresponded to 20% of the composition. Na\\u003csup\\u003e+\\u003c/sup\\u003e is the dominant cation and Cl\\u003csup\\u003e\\u0026minus;\\u003c/sup\\u003e is the dominant anion, while Ca\\u003csup\\u003e2+\\u003c/sup\\u003e, Mg\\u003csup\\u003e2+\\u003c/sup\\u003e, and HCO\\u003csup\\u003e3\\u0026minus;\\u003c/sup\\u003e are not predominant. The chemical facies Na-Cl indicate a low flow in the aquifer that allows the predominance of Na\\u003csup\\u003e+\\u003c/sup\\u003e as a result of the ion exchange process against Ca\\u003csup\\u003e2+\\u003c/sup\\u003e. There was a high concentration of Cl\\u003csup\\u003e\\u0026minus;\\u003c/sup\\u003e in all samples, which suggests an evaporation process in the Beberibe aquifer.\\u003c/p\\u003e\\n \\u003cp\\u003eSargazi et al. (2020), evaluated the quality of groundwater in the city of Zahedan, Sistan, and the province of Baluchestan, in Iran, and the results of the hydrochemical composition of the waters via the Piper diagram, are similar to the results obtained in the study, where the authors pointed out the predominance of Na-Cl-SO\\u003csub\\u003e4\\u003c/sub\\u003e and Ca-Mg-Cl-SO\\u003csub\\u003e4\\u003c/sub\\u003e, in the studied waters. Like Iqbal et al. (2018) and Pazand et al. (2018) studied the hydrochemical composition of groundwater in Abu Dhabi - United Arab Emirates and in the Ardestan basin in Central Iran, respectively, showed results similar to the study.\\u003c/p\\u003e\\n\\u003c/div\\u003e\\n\\u003cdiv class=\\\"Section2\\\" id=\\\"Sec11\\\"\\u003e\\n \\u003ch2\\u003e3.3. Descriptive statistics and geostatistics of the data\\u003c/h2\\u003e\\n \\u003cp\\u003eTable \\u003cspan class=\\\"InternalRef\\\"\\u003e2\\u003c/span\\u003e shows the temporal variability based on descriptive statistics (see section 2.5), for the hydrochemical parameters of the mineral water in the Beberibe aquifer. According to the criteria of Warrick and Nielsen (1980), the CV% presented values of medium too high for all the studied variables, in which the average variation (14.77%) was only for the hydrogen potential (pH), for the other variables had a high variation, CV \\u0026gt; 29%. All parameters showed normal distribution, according to the KS test at the level of 1% probability.\\u003c/p\\u003e\\n \\u003cdiv class=\\\"gridtable\\\"\\u003e\\u003ctable border=\\\"1\\\" id=\\\"Tab2\\\"\\u003e\\n \\u003ccaption language=\\\"En\\\"\\u003e\\n \\u003cdiv class=\\\"CaptionNumber\\\"\\u003eTable 2\\u003c/div\\u003e\\n \\u003cdiv class=\\\"CaptionContent\\\"\\u003e\\n \\u003cp\\u003eDescriptive statistics of the hydrochemical parameters of mineral water.\\u003c/p\\u003e\\n \\u003c/div\\u003e\\n \\u003c/caption\\u003e\\n \\u003ccolgroup cols=\\\"9\\\"\\u003e\\u003c/colgroup\\u003e\\n \\u003cthead\\u003e\\n \\u003ctr\\u003e\\n \\u003cth align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eVariable\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eMean\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eMedian\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e20\\u003c/sup\\u003eMin\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e21\\u003c/sup\\u003eMax\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e22\\u003c/sup\\u003eCV\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e23\\u003c/sup\\u003eSD\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e24\\u003c/sup\\u003eA\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e25\\u003c/sup\\u003eK\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003c/tr\\u003e\\n \\u003c/thead\\u003e\\n \\u003ctbody\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e1\\u003c/sup\\u003eNO\\u003csub\\u003e3\\u003c/sub\\u003e\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e4.53\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.29\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.15\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e43.25\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e195.14\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e8.83\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e3.52\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e12.91\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e2\\u003c/sup\\u003eDR\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e66.17\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e57.90\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e20.00\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e124.98\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e42.13\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e27.87\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.61\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e-0.42\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e3\\u003c/sup\\u003epH\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e5.31\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e5.41\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e4.10\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e7.70\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e14.77\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.78\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.32\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e2.77\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e4\\u003c/sup\\u003eSO\\u003csub\\u003e4\\u003c/sub\\u003e\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e3.77\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e2.03\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.00\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e18.25\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e105.64\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e3.98\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e2.54\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e6.32\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e5\\u003c/sup\\u003eNa\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e9.21\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e8.69\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e5.00\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e19.94\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e29.44\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e2.71\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e2.04\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e6.59\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e6\\u003c/sup\\u003eCa\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.90\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.61\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.00\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e2.66\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e77.37\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.70\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.96\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.12\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e7\\u003c/sup\\u003eMg\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.44\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.00\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.40\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e5.30\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e77.84\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.12\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e2.04\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e4.03\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e8\\u003c/sup\\u003eCl\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e12.27\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e11.97\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.12\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e26.51\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e40.00\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e4.91\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.29\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e3.34\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e9\\u003c/sup\\u003eK\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e4.29\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e4.65\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.00\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e10.45\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e81.45\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e3.50\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.17\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e-1.47\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e10\\u003c/sup\\u003eMn\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.03\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.00\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.00\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.85\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e485.27\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.15\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e5.75\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e33.30\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e11\\u003c/sup\\u003eF\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.04\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.00\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.00\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.85\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e361.79\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.15\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e5.51\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e31.31\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e12\\u003c/sup\\u003eZn\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.01\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.01\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.00\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.08\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e166.80\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.02\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e2.51\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e6.71\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e13\\u003c/sup\\u003eFe\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.06\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.02\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.00\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.38\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e186.85\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.11\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e2.27\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e3.93\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e14\\u003c/sup\\u003eEC\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e86.88\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e78.70\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e40.00\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e197.00\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e41.64\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e36.18\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.68\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e3.14\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e15\\u003c/sup\\u003eHn\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e9.68\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e7.13\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e4.86\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e30.50\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e62.85\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e6.08\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.82\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e3.36\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e16\\u003c/sup\\u003eHCO\\u003csub\\u003e3\\u003c/sub\\u003e\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e8.94\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e8.48\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.00\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e22.87\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e71.13\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e6.36\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.40\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e-0.70\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e17\\u003c/sup\\u003eSi\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e13.43\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e13.99\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.00\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e43.66\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e66.33\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e8.91\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.89\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e2.71\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e18\\u003c/sup\\u003eSr\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.01\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.01\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.00\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.05\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e143.13\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.01\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.85\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e3.28\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e19\\u003c/sup\\u003eBa\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.04\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.01\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.00\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.23\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e138.00\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.06\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.62\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e2.78\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003c/tbody\\u003e\\n \\u003ctfoot\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd colspan=\\\"9\\\"\\u003e\\u003csup\\u003e1\\u003c/sup\\u003eNO\\u003csub\\u003e3\\u003c/sub\\u003e: Nitrate; \\u003csup\\u003e2\\u003c/sup\\u003eDR: Evaporative resistance; \\u003csup\\u003e3\\u003c/sup\\u003epH: Hydrogen potential; \\u003csup\\u003e4\\u003c/sup\\u003eSO\\u003csub\\u003e4\\u003c/sub\\u003e: Sulfate; \\u003csup\\u003e5\\u003c/sup\\u003eNa: Sodium; \\u003csup\\u003e6\\u003c/sup\\u003eCa: Calcium; \\u003csup\\u003e7\\u003c/sup\\u003eMg: Magnesium; \\u003csup\\u003e8\\u003c/sup\\u003eCl: Chlorine; \\u003csup\\u003e9\\u003c/sup\\u003eK: Potassium; \\u003csup\\u003e10\\u003c/sup\\u003eMn: Manganese; \\u003csup\\u003e11\\u003c/sup\\u003eF: Fluorine; \\u003csup\\u003e12\\u003c/sup\\u003eZn: Zinc; \\u003csup\\u003e13\\u003c/sup\\u003eFe: Iron; \\u003csup\\u003e14\\u003c/sup\\u003eEC: Electrical conductivity; \\u003csup\\u003e15\\u003c/sup\\u003eHn: Hardness; \\u003csup\\u003e16\\u003c/sup\\u003eHCO\\u003csub\\u003e3\\u003c/sub\\u003e: Bicarbonate; \\u003csup\\u003e17\\u003c/sup\\u003eSi: Silicon; \\u003csup\\u003e18\\u003c/sup\\u003eSr: Strontium; \\u003csup\\u003e19\\u003c/sup\\u003eBa: Barium; \\u003csup\\u003e20\\u003c/sup\\u003eMin: Minimum; \\u003csup\\u003e21\\u003c/sup\\u003eMax: Maximo; \\u003csup\\u003e22\\u003c/sup\\u003eCV: Coefficient of variation; \\u003csup\\u003e23\\u003c/sup\\u003eSD: Standard deviation; \\u003csup\\u003e24\\u003c/sup\\u003eA: Asymmetry; and \\u003csup\\u003e25\\u003c/sup\\u003eK: Kurtosis.\\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003c/tfoot\\u003e\\n \\u003c/table\\u003e\\n \\u003c/div\\u003e\\n \\u003cp\\u003eThe variability between the medium and high categories is mainly due to the hydrogeological and hydrochemical diversity, which in turn occurs in the geological formation of rocks and which has a direct impact on the chemistry of groundwater (IBGE, 2018; IBGE, 2020). The results obtained corroborate those mentioned by Yang et al. (2020), where they investigated the hydrochemical composition of groundwater to highlight its quality and assess its risk/suitability for irrigation purposes in the Ordos Basin, China, they identified high CV% for all variables studied, except for the pH, where the CV% presented a value referring to the medium variability.\\u003c/p\\u003e\\n \\u003cp\\u003eThe highest CV% indicated that the ions were sensitive to the external environment, such as hydrological conditions, topography, and anthropogenic activities, as discussed by Yang et al. (2020). Another key factor is pH, which according to Liu et al. (2018), in a study of the dynamic characteristics of groundwater in the valley plain of the city of Lhasa (capital of the Autonomous Region of Tibet), the low CV% indicated that the pH of groundwater tends to be relatively stable, which differs from the present study, in which the pH was not stable and CV% is admitted with medium variability, according to the criterion of Warrick and Nielsen (1980) - (Table \\u003cspan class=\\\"InternalRef\\\"\\u003e2\\u003c/span\\u003e).\\u003c/p\\u003e\\n \\u003cp\\u003eThe pH according to the standards established in Table \\u003cspan class=\\\"InternalRef\\\"\\u003e1\\u003c/span\\u003e and compliance with the minimum values recorded in Table \\u003cspan class=\\\"InternalRef\\\"\\u003e2\\u003c/span\\u003e, presented low standards, indicative of points of distribution of acidic mineral water for human health. High concentrations of NO\\u003csub\\u003e3\\u003c/sub\\u003e were observed in some mineral water distribution points, therefore, outside the limits established by Brazilian legislation (Table \\u003cspan class=\\\"InternalRef\\\"\\u003e1\\u003c/span\\u003e), probably caused by the socio-economic development of the study region, as highlighted earlier by Liu et al. (2018).\\u003c/p\\u003e\\n \\u003cp\\u003eThe electrical conductivity (EC) showed high CV%, due to the oscillation of salinity in urban areas (natural effect due to the rock/soil factor, and anthropization), such variations are characteristics of low permeability and low hydraulic gradient, due to the dilution effect salts during the rainy season and evaporation in the dry season. The results similar to those obtained by Tlili-Zrelli et al. (2015), they studied the hydrogeochemistry of groundwater in an alluvial aquifer in northern Tunisia, where the zones of low permeability and low hydraulic gradient, occur due to the dilution effect of salts during the rainy season and evaporation during the dry season, which favors a high variability in the levels of salts diluted in the water due to the climate and the greater and/or less water availability.\\u003c/p\\u003e\\n \\u003cp\\u003eThe analysis of the spatial variability for the hydrochemical parameters is shown in Table \\u003cspan class=\\\"InternalRef\\\"\\u003e3\\u003c/span\\u003e. All the parameters evaluated, the Gaussian model was the one that best fit, which is indicated for studies of groundwater. The highlight for the coefficient R\\u003csup\\u003e2\\u003c/sup\\u003e with variations between 0.6 to 0.99.\\u003c/p\\u003e\\u003ctable border=\\\"1\\\" id=\\\"Tab3\\\"\\u003e\\n \\u003ccaption language=\\\"En\\\"\\u003e\\n \\u003cdiv class=\\\"CaptionNumber\\\"\\u003eTable 3\\u003c/div\\u003e\\n \\u003cdiv class=\\\"CaptionContent\\\"\\u003e\\n \\u003cp\\u003eAnalysis of spatial variability followed by adjusted semivariograms for the hydrochemical parameters of mineral water.\\u003c/p\\u003e\\n \\u003c/div\\u003e\\n \\u003c/caption\\u003e\\n \\u003cthead\\u003e\\n \\u003ctr\\u003e\\n \\u003cth align=\\\"left\\\" rowspan=\\\"2\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e20\\u003c/sup\\u003eVar\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\" rowspan=\\\"2\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e21\\u003c/sup\\u003eMo\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\" rowspan=\\\"2\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e22\\u003c/sup\\u003eC\\u003csub\\u003e0\\u003c/sub\\u003e\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\" rowspan=\\\"2\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e23\\u003c/sup\\u003eC\\u003csub\\u003e0\\u003c/sub\\u003e+C\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\" rowspan=\\\"2\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e24\\u003c/sup\\u003e\\u003cem\\u003ea\\u003c/em\\u003e\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\" rowspan=\\\"2\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e25\\u003c/sup\\u003eR\\u003csup\\u003e2\\u003c/sup\\u003e\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e26\\u003c/sup\\u003eC\\u003csub\\u003e0\\u003c/sub\\u003e\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\" rowspan=\\\"2\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e27\\u003c/sup\\u003eDSD\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\" colspan=\\\"2\\\"\\u003e\\n \\u003cp\\u003eJack-Knifing\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003cth align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e(C\\u003csub\\u003e0\\u003c/sub\\u003e+C)\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e28\\u003c/sup\\u003em\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e29\\u003c/sup\\u003eSD\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003c/tr\\u003e\\n \\u003c/thead\\u003e\\n \\u003ctbody\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e1\\u003c/sup\\u003eNO\\u003csub\\u003e3\\u003c/sub\\u003e\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eGau\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e2.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e10.500\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e25000.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.800\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e19.047\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eFt\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e-0.082\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.017\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e2\\u003c/sup\\u003eDR\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eGau\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e90.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e1396.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e4886.115\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.849\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e6.446\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eFt\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.006\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.007\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e3\\u003c/sup\\u003epH\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eGau\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e0.128\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e1.975\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e14440.107\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.933\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e6.481\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eFt\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.070\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.998\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e4\\u003c/sup\\u003eSO\\u003csub\\u003e4\\u003c/sub\\u003e\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eGau\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e0.200\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e3.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e100.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.599\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e6.667\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eFt\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e-0.087\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.113\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e5\\u003c/sup\\u003eNa\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eGau\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e3.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e18.900\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e20000.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.728\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e15.873\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eFt\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.012\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.044\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e6\\u003c/sup\\u003eCa\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eGau\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e0.120\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e0.491\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e417.424\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.713\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e24.439\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eFt\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e-0.048\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.041\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e7\\u003c/sup\\u003eMg\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eGau\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e0.136\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e1.292\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e3360.178\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.908\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e10.526\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eFt\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e-0.011\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.940\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e8\\u003c/sup\\u003eCl\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eGau\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e5.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e28.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e13000.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.701\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e17.857\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eFt\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.050\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.049\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e9\\u003c/sup\\u003eK\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eGau\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e4.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e22.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e5000.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.708\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e18.182\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eFt\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.029\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.061\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e10\\u003c/sup\\u003eMn\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eGau\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e1.50E-05\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e1.54E-04\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e16956.777\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.657\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e9.740\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eFt\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.027\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.081\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e11\\u003c/sup\\u003eF\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eGau\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e1.00E-06\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e1.00E-02\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1000.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.989\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e0.010\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eFt\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e-0.056\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.112\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e12\\u003c/sup\\u003eZn\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eGau\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e2.00E-04\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e0.724\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e8861.172\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.935\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e0.027\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eFt\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.049\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.032\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e13\\u003c/sup\\u003eFe\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eGau\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e1.00E-6\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e0.002\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e169.741\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.917\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e4.234E-4\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eFt\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.017\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.028\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e14\\u003c/sup\\u003eEC\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eGau\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e80.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e557.100\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e2468.172\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.961\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e14.360\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eFt\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e-0.057\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.070\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e15\\u003c/sup\\u003eHn\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eGau\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e12.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e60.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e12500.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.790\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e20.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eFt\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.017\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.071\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e16\\u003c/sup\\u003eHCO\\u003csub\\u003e3\\u003c/sub\\u003e\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eGau\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e25.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e110.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e4200.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.812\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e22.727\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eFt\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.023\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.053\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e17\\u003c/sup\\u003eSi\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eGau\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e50.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e210.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e7000.000\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.905\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e23.810\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eFt\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.049\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.004\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e18\\u003c/sup\\u003eSr\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eGau\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e5.00E-06\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e7.30E-05\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e803.671\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.867\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e6.849\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eFt\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e-0.053\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.246\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e\\u003csup\\u003e19\\u003c/sup\\u003eBa\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eGau\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e8.00E-04\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e0.021\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e4588.202\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e0.630\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003e3.810\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eFt\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e-0.010\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e1.080\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003c/tbody\\u003e\\n \\u003ctfoot\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd colspan=\\\"10\\\"\\u003e\\u003csup\\u003e1\\u003c/sup\\u003eNO\\u003csub\\u003e3\\u003c/sub\\u003e: Nitrate; \\u003csup\\u003e2\\u003c/sup\\u003eDR: Evaporative resistance; \\u003csup\\u003e3\\u003c/sup\\u003epH: Hydrogen potential; \\u003csup\\u003e4\\u003c/sup\\u003eSO\\u003csub\\u003e4\\u003c/sub\\u003e: Sulfate; \\u003csup\\u003e5\\u003c/sup\\u003eNa: Sodium; \\u003csup\\u003e6\\u003c/sup\\u003eCa: Calcium; \\u003csup\\u003e7\\u003c/sup\\u003eMg: Magnesium; \\u003csup\\u003e8\\u003c/sup\\u003eCl: Chlorine; \\u003csup\\u003e9\\u003c/sup\\u003eK: Potassium; \\u003csup\\u003e10\\u003c/sup\\u003eMn: Manganese; \\u003csup\\u003e11\\u003c/sup\\u003eF: Fluorine; \\u003csup\\u003e12\\u003c/sup\\u003eZn: Zinc; \\u003csup\\u003e13\\u003c/sup\\u003eFe: Iron; \\u003csup\\u003e14\\u003c/sup\\u003eEC: Electrical conductivity; \\u003csup\\u003e15\\u003c/sup\\u003eHn: Hardness; \\u003csup\\u003e16\\u003c/sup\\u003eHCO\\u003csub\\u003e3\\u003c/sub\\u003e: Bicarbonate; \\u003csup\\u003e17\\u003c/sup\\u003eSi: Silicon; \\u003csup\\u003e18\\u003c/sup\\u003eSr: Strontium; \\u003csup\\u003e19\\u003c/sup\\u003eBa: Barium; \\u003csup\\u003e20\\u003c/sup\\u003eVar: Variable; \\u003csup\\u003e21\\u003c/sup\\u003eMo: Model; \\u003csup\\u003e22\\u003c/sup\\u003eC\\u003csub\\u003e0\\u003c/sub\\u003e: Nugget effect; \\u003csup\\u003e23\\u003c/sup\\u003eC\\u003csub\\u003e0\\u003c/sub\\u003e+C: Sill; \\u003csup\\u003e24\\u003c/sup\\u003e\\u003cem\\u003ea\\u003c/em\\u003e: Range; \\u003csup\\u003e25\\u003c/sup\\u003eR\\u003csup\\u003e2\\u003c/sup\\u003e: Coefficient of determination; \\u003csup\\u003e26\\u003c/sup\\u003eC\\u003csub\\u003e0\\u003c/sub\\u003e/(C\\u003csub\\u003e0\\u003c/sub\\u003e+C): Degree of spatial dependence (%); \\u003csup\\u003e27\\u003c/sup\\u003eDSD: Class of the degree of spatial dependence; \\u003csup\\u003e28\\u003c/sup\\u003em: Mean; \\u003csup\\u003e29\\u003c/sup\\u003eSD: Standard deviation; Gau: Gaussian; Ft: Strong.\\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003c/tfoot\\u003e\\n \\u003c/table\\u003e\\n \\u003cp\\u003eYin et al. (2019), analyzed geostatistically the hydrochemical variations and the causes of NO\\u003csub\\u003e3\\u003c/sub\\u003e pollution of groundwater in a plain in southern Beijing, China. They highlighted the need for interactions between the transitive theoretical models to validate the best for each parameter. However, in the present study, cross-validation by Jack-Knifing (Vauclin et al., 1983) was applied, were all established semivariograms presented the mean close to 0 and the standard deviation close to 1, and thus used in kriging maps for each studied hydrochemical variable (Table \\u003cspan class=\\\"InternalRef\\\"\\u003e3\\u003c/span\\u003e).\\u003c/p\\u003e\\n \\u003cp\\u003eThe DSD was strong for all the parameters analyzed, which indicates that the value of each parameter associated with a given location, was more similar to the value of a neighboring sample than to the rest of the location of the sample set, which reinforces the validation of the Gaussian model for the studied hydrochemical variables (Cambardella et al., \\u003cspan class=\\\"CitationRef\\\"\\u003e1994\\u003c/span\\u003e).\\u003c/p\\u003e\\n \\u003cp\\u003eFigure \\u003cspan class=\\\"InternalRef\\\"\\u003e6\\u003c/span\\u003e illustrates the spatial distribution of the hydrochemical parameters of the mineral water sources in the Beberibe aquifer in the study area. Based on the semivariograms, kriging of the maps was performed, according to each model generated for the hydrochemical variables of the groundwater, implementing the nugget effect, level, and range of each model.\\u003c/p\\u003e\\n \\u003cp\\u003e\\u003cstrong\\u003e[INSERT TO\\u003c/strong\\u003e FIGURE \\u003cspan class=\\\"InternalRef\\\"\\u003e6\\u003c/span\\u003e \\u003cstrong\\u003eHERE]\\u003c/strong\\u003e\\u003c/p\\u003e\\n \\u003cp\\u003eNO\\u003csub\\u003e3\\u003c/sub\\u003e concentrations occurred mainly on the coast of Pernambuco, with an emphasis on the municipalities of Recife and Olinda (Figure \\u003cspan class=\\\"InternalRef\\\"\\u003e6\\u003c/span\\u003e). In both municipalities, the NO\\u003csub\\u003e3\\u003c/sub\\u003e concentrations were higher, mainly because these municipalities have a high population density and have a high flow of tourists, vegetable cultivation near rivers (causing eutrophication), as well as the holding of major cultural events, such as Carnival in the city of Recife and Olinda, Brazil, considered cities that host the carnival festival with the largest flow of people in the world, with an average of 1 million people a day. Within this theme, one of the factors that imply high concentrations of NO\\u003csub\\u003e3\\u003c/sub\\u003e is caused by the high pollution generated by tourists and the consequent contamination of the main rivers, highlighting Capibaribe, Beberibe, and Tejipi\\u0026oacute;.\\u003c/p\\u003e\\n \\u003cp\\u003eAlso noteworthy is the fact that the city of Recife is among the 10 largest Brazilian metropolises, with a population of 1.555 million inhabitants, according to the IBGE - (IBGE, 2020). In this context, the rate of pollution generated in the urban centers of the metropolis, associated with the lack of sanitation in neighborhoods, provides for the contamination of rivers, which in turn contribute to recharge the aquifer, which implies the contamination of groundwater and the increase of NO\\u003csub\\u003e3\\u003c/sub\\u003e concentration. Therefore, there is still the spread of waterborne diseases, mainly diarrhea, typhoid fever, hepatitis A, leptospirosis, cholera, and intestinal infections, caused by bacteria.\\u003c/p\\u003e\\n \\u003cp\\u003eConcomitant with NO\\u003csub\\u003e3\\u003c/sub\\u003e concentrations, it is observed in the maps of the hydrochemical variables that the evaporative resistance (DR), sodium (Na), calcium (Ca), magnesium (Mg), chlorine (Cl), electrical conductivity (EC), hardness (Hn) and strontium (Sr), presented higher concentrations in the places with the highest concentration of NO\\u003csub\\u003e3\\u003c/sub\\u003e (Figure \\u003cspan class=\\\"InternalRef\\\"\\u003e6\\u003c/span\\u003e). From the coast that comprises the municipality of Recife and Olinda, it is where the greatest concentrations of these elements are found, among them are heavy metals and EC, which is associated with greater contractions of soluble solids in water, which in turn correlate with anthropic pollution.\\u003c/p\\u003e\\n \\u003cp\\u003eCorroborating the findings of the present study, Bodrud-Doza et al. (\\u003cspan class=\\\"CitationRef\\\"\\u003e2016\\u003c/span\\u003e), carried out the characterization of groundwater quality based on water assessment indices, multivariate statistics, and geostatistics in central Bangladesh, and concluded that anthropogenic factors (e.g., industrial pollution, inadequate garbage discharge, inadequate sanitation) contribute groundwater quality. Bhuiyan et al. (\\u003cspan class=\\\"CitationRef\\\"\\u003e2016\\u003c/span\\u003e) also assessed the groundwater quality of the Lakhimpur district in Bangladesh based on water quality indices, geostatistical methods, and multivariate analysis pointed out that a kriging method is an effective tool for decision-makers on groundwater quality management.\\u003c/p\\u003e\\n \\u003cp\\u003eIt is noteworthy that the hydrogeological formation also contributes to NO\\u003csub\\u003e3\\u003c/sub\\u003e concentrations. It appears that the Pernambuco coast has a hydrogeological formation of the mixed bicarbonate type (Figure \\u003cspan class=\\\"InternalRef\\\"\\u003e4\\u003c/span\\u003e), which because of this, there is a high variability (\\u0026gt; CV%) of the hydrochemical elements of groundwater (Table \\u003cspan class=\\\"InternalRef\\\"\\u003e2\\u003c/span\\u003e). This mixed diversity of bicarbonate can also be considered as one of the factors for this high concentration of NO\\u003csub\\u003e3\\u003c/sub\\u003e, as such regions are common of fissural geological formation, which in turn reflects and a greater tendency for contamination of these waters.\\u003c/p\\u003e\\n \\u003cp\\u003eAccording to the Minist\\u0026eacute;rio da Sa\\u0026uacute;de (1999), Minist\\u0026eacute;rio da Sa\\u0026uacute;de (2006), and CETESB (2016), concentrations of NO\\u003csub\\u003e3\\u003c/sub\\u003e, DR, Na, Ca, Mg, Cl, EC, Hn and Sr, only NO\\u003csub\\u003e3\\u003c/sub\\u003e and Cl did not comply with Brazilian legislation, with values higher than allowed (Table \\u003cspan class=\\\"InternalRef\\\"\\u003e1\\u003c/span\\u003e). However, to the west of the coast, values were observed within the limits allowed for NO\\u003csub\\u003e3\\u003c/sub\\u003e and Cl, which reinforces the anthropogenic impact in the coastal region (Figure \\u003cspan class=\\\"InternalRef\\\"\\u003e6\\u003c/span\\u003e). Corroborating the findings of Yin et al. (2019), who claim that industrial areas and intense urbanization promote greater groundwater pollution by NO\\u003csub\\u003e3\\u003c/sub\\u003e.\\u003c/p\\u003e\\n \\u003cp\\u003eThe pH of groundwater on the coast met the acceptable standard for human health (pH between 6 and 9), with variations between acid to alkaline (pH between 6.4 and 7.6) (Minist\\u0026eacute;rio da Sa\\u0026uacute;de, 1999; Minist\\u0026eacute;rio da Sa\\u0026uacute;de, 2006; CETESB, 2016). However, the west region of the coast had a pH \\u0026lt; 6 (Figure \\u003cspan class=\\\"InternalRef\\\"\\u003e6\\u003c/span\\u003e). Abu-alnaeem et al. (\\u003cspan class=\\\"CitationRef\\\"\\u003e2018\\u003c/span\\u003e), evaluated the salinity and quality of groundwater in the Gaza coastal aquifer, Gaza Strip in Palestine, and highlighted that acidity and alkalinity close to neutralization, is derived from the dissolution of carbonate minerals in the form of bicarbonate (HCO\\u003csup\\u003e3\\u0026minus;\\u003c/sup\\u003e). Corroborating the statement of the aforementioned authors, the coastal region of the present study, to the east, showed pH with variations between acid and alkaline, still because it presents mixed bicarbonate hydrogeological formation (Figure \\u003cspan class=\\\"InternalRef\\\"\\u003e4\\u003c/span\\u003e).\\u003c/p\\u003e\\n \\u003cp\\u003eThe elements SO\\u003csub\\u003e4\\u003c/sub\\u003e, K, Mn, F, Zn, Fe, HCO\\u003csub\\u003e3\\u003c/sub\\u003e, Si, Sr and Ba were the elements that had the highest concentrations west of the coast (Figure \\u003cspan class=\\\"InternalRef\\\"\\u003e6\\u003c/span\\u003e). This study region was also the one with the lowest pH. It is noteworthy that only the elements SO\\u003csub\\u003e4\\u003c/sub\\u003e, F, Zn, and Ba, presented concentrations within the quality standards for human consumption (Minist\\u0026eacute;rio da Sa\\u0026uacute;de, 1999; Minist\\u0026eacute;rio da Sa\\u0026uacute;de, 2006; CETESB, 2016). The elements K, Mn, and Fe showed values above that allowed by current legislation. HCO\\u003csub\\u003e3\\u003c/sub\\u003e, Si, and Sr, on the other hand, do not have a defined quality standard.\\u003c/p\\u003e\\n \\u003cp\\u003eHossain and Patra, (2020) carried out a contamination zoning and health risk assessment of trace elements in groundwater utilizing geostatistical modeling and reported that pH controls chemistry, which is a component that influences the hydrochemical standards of water and, thus, corroborates with the observations of the present study, in which the elements K, Mn, and Fe presented higher values for the region of lower pH.\\u003c/p\\u003e\\n\\u003c/div\\u003e\\n\\u003cdiv class=\\\"Section2\\\" id=\\\"Sec12\\\"\\u003e\\n \\u003ch2\\u003e3.4. Principal component analysis (PCA)\\u003c/h2\\u003e\\n \\u003cp\\u003eThe PCA presents the eigenvalue, the proportion, and the accumulation of the total variance for all variables studied (Table \\u003cspan class=\\\"InternalRef\\\"\\u003e4\\u003c/span\\u003e). The cumulative total variance for principal component 2 (PC2) was 49.61%, a significant result for the analyzed data set. According to the criteria established by Kaiser (1958), PC1 and PC2 presented their eigenvalues greater than 1, which indicated a significant information load and, thus, generate the biplot graphs.\\u003c/p\\u003e\\n \\u003cdiv class=\\\"gridtable\\\"\\u003e\\u003ctable border=\\\"1\\\" id=\\\"Tab4\\\"\\u003e\\n \\u003ccaption language=\\\"En\\\"\\u003e\\n \\u003cdiv class=\\\"CaptionNumber\\\"\\u003eTable 4\\u003c/div\\u003e\\n \\u003cdiv class=\\\"CaptionContent\\\"\\u003e\\n \\u003cp\\u003ePrincipal components of the hydrochemical parameters of mineral water.\\u003c/p\\u003e\\n \\u003c/div\\u003e\\n \\u003c/caption\\u003e\\n \\u003ccolgroup cols=\\\"4\\\"\\u003e\\u003c/colgroup\\u003e\\n \\u003cthead\\u003e\\n \\u003ctr\\u003e\\n \\u003cth align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eComponents\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eEigenvalue\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eProportion (%)\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003cth align=\\\"left\\\"\\u003e\\n \\u003cp\\u003eCumulative (%)\\u003c/p\\u003e\\n \\u003c/th\\u003e\\n \\u003c/tr\\u003e\\n \\u003c/thead\\u003e\\n \\u003ctbody\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003ePC1\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e6.30\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e33.14\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e33.14\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd align=\\\"left\\\"\\u003e\\n \\u003cp\\u003ePC2\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e3.13\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e16.46\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003ctd align=\\\"char\\\"\\u003e\\n \\u003cp\\u003e49.61\\u003c/p\\u003e\\n \\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003c/tbody\\u003e\\n \\u003ctfoot\\u003e\\n \\u003ctr\\u003e\\n \\u003ctd colspan=\\\"4\\\"\\u003ePC: Principal component.\\u003c/td\\u003e\\n \\u003c/tr\\u003e\\n \\u003c/tfoot\\u003e\\n \\u003c/table\\u003e\\n \\u003c/div\\u003e\\n \\u003cp\\u003eYang et al. (2016) carried out the identification of hydrogeochemical processes and assessment of groundwater quality in the southeastern part of the Ordos basin in China, based on the PCA and obtained total cumulative variance for the PC2, of 59.26%, which is higher than the study.\\u003c/p\\u003e\\n \\u003cp\\u003eIn Figure \\u003cspan class=\\\"InternalRef\\\"\\u003e7\\u003c/span\\u003e, the PC\\u0026rsquo;s biplot chart stands out and Pearson correlation of the hydrochemical parameters of the mineral water of the Beberibe aquifer in the study area, referring to the municipalities of Recife, Olinda, Paulista, Camaragibe, S\\u0026atilde;o Louren\\u0026ccedil;o da Mata, and Paudalho.\\u003c/p\\u003e\\n \\u003cp\\u003e\\u003cstrong\\u003e[INSERT TO\\u003c/strong\\u003e FIGURE \\u003cspan class=\\\"InternalRef\\\"\\u003e7\\u003c/span\\u003e \\u003cstrong\\u003eHERE]\\u003c/strong\\u003e\\u003c/p\\u003e\\n \\u003cp\\u003eFor the PC\\u0026apos;s, the 19 studied variables (NO\\u003csub\\u003e3\\u003c/sub\\u003e; DR; pH; SO\\u003csub\\u003e4\\u003c/sub\\u003e; Na; Ca; Mg; Cl; K; Mn; F; Zn; Fe; EC; Hn; HCO\\u003csub\\u003e3\\u003c/sub\\u003e; Si; Sr and Ba) were admitted and the correlated points refer to the artesian sample collection wells (A01 to A34). The variables DR, Na, Ca, Mg, Cl, EC, Hn, and Sr, showed a correlation with NO\\u003csub\\u003e3\\u003c/sub\\u003e, as observed in kriging maps (Figure \\u003cspan class=\\\"InternalRef\\\"\\u003e7\\u003c/span\\u003eA). However, it is noteworthy that among the observed correlations, Cl was the one with the highest correlation with NO\\u003csub\\u003e3\\u003c/sub\\u003e, in which, the Cl projection line follows the NO\\u003csub\\u003e3\\u003c/sub\\u003e line, which characterizes this strong correlation (Figure \\u003cspan class=\\\"InternalRef\\\"\\u003e7\\u003c/span\\u003eB).\\u003c/p\\u003e\\n \\u003cp\\u003eSimilar results were found by Cao et al. (\\u003cspan class=\\\"CitationRef\\\"\\u003e2020\\u003c/span\\u003e), where they studied the heterogeneity of an aquifer in a region of France, through PCA and observed a high correlation between Cl and NO\\u003csub\\u003e3\\u003c/sub\\u003e. The same results were found by Ayed et al. (\\u003cspan class=\\\"CitationRef\\\"\\u003e2017\\u003c/span\\u003e), where they performed a hydrochemical characterization of groundwater using PCA for the southeastern region of Tunisia, and observed the same correlation between Cl and NO\\u003csub\\u003e3\\u003c/sub\\u003e.\\u003c/p\\u003e\\n \\u003cp\\u003eTo clarify more precisely the correlations between the variables, Pearson\\u0026apos;s correlation was applied, where the correlation was significant between Cl and NO\\u003csub\\u003e3\\u003c/sub\\u003e (Figure \\u003cspan class=\\\"InternalRef\\\"\\u003e7\\u003c/span\\u003eB). Corroborating the findings of this study, Ayed et al. (\\u003cspan class=\\\"CitationRef\\\"\\u003e2017\\u003c/span\\u003e) also observed a strong correlation between Cl and NO\\u003csub\\u003e3\\u003c/sub\\u003e.\\u003c/p\\u003e\\n\\u003c/div\\u003e\\n\\u003cdiv class=\\\"Section2\\\" id=\\\"Sec13\\\"\\u003e\\n \\u003ch2\\u003e3.5. Linear regression analysis\\u003c/h2\\u003e\\n \\u003cp\\u003eBased on the analysis of the kriging maps (Figure \\u003cspan class=\\\"InternalRef\\\"\\u003e6\\u003c/span\\u003e) and the PCA - (Figure \\u003cspan class=\\\"InternalRef\\\"\\u003e7\\u003c/span\\u003eA), was performed a simple linear regression analysis of the electrical conductivity (EC), hardness (Hn), evaporative resistance (DR), magnesium (Mg), sodium (Na) and chlorine (Cl), with nitrate (NO\\u003csub\\u003e3\\u003c/sub\\u003e), as they present a greater correlation with NO\\u003csub\\u003e3\\u003c/sub\\u003e in the previous analyzes (Figure \\u003cspan class=\\\"InternalRef\\\"\\u003e8\\u003c/span\\u003e). The R\\u003csup\\u003e2\\u003c/sup\\u003e criterion greater than 0.5 was established to represent only the most representative predictor variables with NO\\u003csub\\u003e3\\u003c/sub\\u003e.\\u003c/p\\u003e\\n \\u003cp\\u003e\\u003cstrong\\u003e[INSERT TO\\u003c/strong\\u003e FIGURE \\u003cspan class=\\\"InternalRef\\\"\\u003e8\\u003c/span\\u003e \\u003cstrong\\u003eHERE]\\u003c/strong\\u003e\\u003c/p\\u003e\\n \\u003cp\\u003eAmong the linear regressions, only NO\\u003csub\\u003e3\\u003c/sub\\u003e x DR showed R\\u003csup\\u003e2\\u003c/sup\\u003e \\u0026lt; 0.5 (R\\u003csup\\u003e2\\u003c/sup\\u003e = 0.32), so the DR has no direct influence on NO\\u003csub\\u003e3\\u003c/sub\\u003e (Figure \\u003cspan class=\\\"InternalRef\\\"\\u003e8\\u003c/span\\u003eC). It can be seen that Cl, showed a satisfactory R\\u003csup\\u003e2\\u003c/sup\\u003e with NO\\u003csub\\u003e3\\u003c/sub\\u003e (R\\u003csup\\u003e2\\u003c/sup\\u003e = 0.51) - (Figure \\u003cspan class=\\\"InternalRef\\\"\\u003e8\\u003c/span\\u003eF), previously indicated by geostatistics based on kriging maps and PCA, and these two variables showed a high correlation between themselves. The results obtained corroborate those described by Tiwari et al. (2017), where they assessed the spatial distribution of groundwater quality and assessed regional suitability for drinking water in the Alwar region in India and found a strong correlation between NO\\u003csub\\u003e3\\u003c/sub\\u003e and Cl.\\u003c/p\\u003e\\n \\u003cp\\u003eThe variables EC (R\\u003csup\\u003e2\\u003c/sup\\u003e = 0.68), Hn (R\\u003csup\\u003e2\\u003c/sup\\u003e = 0.57), Mg (R\\u003csup\\u003e2\\u003c/sup\\u003e = 0.67) and Na (R\\u003csup\\u003e2\\u003c/sup\\u003e = 0.56), showed satisfactory R\\u003csup\\u003e2\\u003c/sup\\u003e, respectively. Such variables have a direct correlation with NO\\u003csub\\u003e3\\u003c/sub\\u003e, thus being potential variables in the hydrochemical characterization of groundwater.\\u003c/p\\u003e\\n\\u003c/div\\u003e\"},{\"header\":\"4. Conclusions\",\"content\":\"\\u003cp\\u003eBased on Brazilian legislation, there are strategic wells that have a direct impact on human health and, therefore, there is a need to reanalyze possible approaches for improvements that support health and quality of life. The composition and hydrochemical quality of groundwater in the metropolitan region of Recife and the coast of Pernambuco are directly influenced by hydrogeology, topographic conditions, and urban infrastructure (urban density).\\u003c/p\\u003e \\u003cp\\u003eIn geostatistical modeling, based on the Gaussian model, it stands out concerning the other transitive theoretical models in the hydrochemical characterization of the waters of the Beberibe aquifer and thus spatially displays the spatial dynamics of water composition and quality in the study area with greater accuracy. Multivariate analysis via PCA associated with Pearson's correlation is a promising technique for assessing the hydrochemistry of groundwater in the Beberibe aquifer together with the hydrogeological and hydrochemical maps of groundwater made it possible to efficiently characterize water quality standards for the coast of Pernambuco.\\u003c/p\\u003e \\u003cp\\u003eThe anthropogenic influence of the metropolitan center of Recife, has a direct impact on water quality for human health, especially on nitrate concentrations (NO\\u003csub\\u003e3\\u003c/sub\\u003e). It is also noteworthy that the concentrations of chlorine (Cl) have a high correlation with NO\\u003csub\\u003e3\\u003c/sub\\u003e, being confirmed in the kriging maps and the PCA, Pearson's correlation, and in the linear regression. Other methodologies must be formulated and applied, to confirm or rule out the possibility of recharging the regional flow of groundwater. The methodology adopted in the study can be perfectly applied in other regions of the planet with conditions of social vulnerability and problems with water quality, with bias of monitoring and definition of public policies.\\u003c/p\\u003e \\u003cp\\u003eLow-quality underground water should not be used for human consumption, as it becomes the main agent for the emergence of waterborne diseases, mainly diarrhea, typhoid fever, hepatitis A, leptospirosis, cholera, and intestinal infections, caused by bacteria, which it reaches from infants to the elderly, especially those who live in areas of social vulnerability and with low-quality water. Therefore, regular monitoring is essential to ensure good quality water and a healthy environment that provides quality of life for local families.\\u003c/p\\u003e\"},{\"header\":\"Declarations\",\"content\":\"\\u003cp\\u003e\\u003cstrong\\u003e5. Acknowledgment\\u003c/strong\\u003e\\u003c/p\\u003e\\n\\u003cp\\u003eTo the \\u003cem\\u003ePrograma de P\\u0026oacute;s-Gradua\\u0026ccedil;\\u0026atilde;o em Engenharia Agr\\u0026iacute;cola\\u003c/em\\u003e (PGEA), \\u003cem\\u003ePrograma de P\\u0026oacute;s-gradua\\u0026ccedil;\\u0026atilde;o em Qu\\u0026iacute;mica\\u003c/em\\u003e (PPGQ), and the \\u003cem\\u003eGrupo de Pesquisa em Ambi\\u0026ecirc;ncia\\u003c/em\\u003e (GPESA) of the Federal Rural University of Pernambuco (UFRPE) for supporting the development of this research. The \\u003cem\\u003eCoordena\\u0026ccedil;\\u0026atilde;o de Aperfei\\u0026ccedil;oamento de Pessoal de N\\u0026iacute;vel Superior\\u003c/em\\u003e (CAPES - Finance Code 001) and the \\u003cem\\u003eFunda\\u0026ccedil;\\u0026atilde;o de Amparo \\u0026agrave; Ci\\u0026ecirc;ncia e Tecnologia do Estado de Pernambuco\\u003c/em\\u003e (FACEPE), for the financing of scholarships. The third and fourth authors are grateful for the Research Productivity grant from the \\u003cem\\u003eConselho Nacional de Desenvolvimento Cient\\u0026iacute;fico e Tecnol\\u0026oacute;gico\\u003c/em\\u003e (CNPq) level 2 under process number 304060/2016-0 and 309681/2019-7, respectively.\\u003c/p\\u003e\"},{\"header\":\"References\",\"content\":\"\\u003col\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eABINAM (2020) \\u003cem\\u003eAssocia\\u0026ccedil;\\u0026atilde;o Brasileira Ind\\u0026uacute;stria \\u0026Aacute;guas Minerais.\\u003c/em\\u003e [online]. Available at:. [Accessed 15 July 2020]\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eAbu-alnaeem MF, Yusoff I, Ng TF, Alias Y, Raksmey M (2018) Assessment of groundwater salinity and quality in Gaza coastal aquifer, Gaza Strip, Palestine: An integrated statistical, geostatistical and hydrogeochemical approaches study. Sci Total Environ 615:972\\u0026ndash;989. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.scitotenv.2017.09.320\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eAg\\u0026ecirc;ncia Pernambucana de Vigil\\u0026acirc;ncia Sanit\\u0026aacute;ria \\u0026ndash; APEVISA (2020) Programa de Monitoramento e Fiscaliza\\u0026ccedil;\\u0026atilde;o da \\u0026Aacute;gua Envasada (Pura). [online]. Available at:. [Accessed 11 January 2021]\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eAnders R, Mendez GO, Futa K, Danskin WR (2014) A geochemical approach to determine sources and movement of saline groundwater in a coastal aquifer. Groundwater 52(5):756\\u0026ndash;768. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1111/gwat.12108\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eArshad S, Mobeen M, Bashir S, Aziz T, Rehman A (2019) Spatializing Groundwater Quality Parameters and Their Impacts on Land Value in Khushab City, Punjab, Pakistan. International Journal of Economic and Environmental Geology 10(3):78\\u0026ndash;84. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.46660/ijeeg.vol11.iss1.1919\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eAyed B, Jmal I, Sahal S, Mokadem N, Saidi S, Boughariou E, Bouri S (2017) Hydrochemical characterization of groundwater using multivariate statistical analysis: the Maritime Djeffara shallow aquifer (Southeastern Tunisia). Environ Earth Sci 76(24):821. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1007/s12665-017-7168-6\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eBeck HE, Zimmermann NE, McVicar TR, Vergopolan N, Berg A, Wood EF (2018) Present and future K\\u0026ouml;ppen-Geiger climate classification maps at 1-km resolution. Scientific Data 5:180214. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1038/sdata.2018.214\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eBhuiyan MAH, Bodrud-Doza M, Islam AT, Rakib MA, Rahman MS, Ramanathan AL (2016) Assessment of groundwater quality of Lakshimpur district of Bangladesh using water quality indices, geostatistical methods, and multivariate analysis. Environ Earth Sci 75(12):1020. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1007/s12665-016-5823-y\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eBodrud-Doza M, Bhuiyan MAH, Islam SDU, Rahman MS, Haque MM, Fatema KJ, Rahman MA (2019) Hydrogeochemical investigation of groundwater in Dhaka City of Bangladesh using GIS and multivariate statistical techniques. Groundwater for Sustainable Development 8:226\\u0026ndash;244. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.gsd.2018.11.008\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eBodrud-Doza M, Islam AT, Ahmed F, Das S, Saha N, Rahman MS (2016) Characterization of groundwater quality using water evaluation indices, multivariate statistics and geostatistics in central Bangladesh. Water Science 30(1):19\\u0026ndash;40. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.wsj.2016.05.001\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eBulia IL, Enzweiler J (2018) The hydrogeochemistry of bottled mineral water in S\\u0026atilde;o Paulo state, Brazil. J Geochem Explor 188:43\\u0026ndash;54. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.gexplo.2018.01.007\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eCambardella CA, Moorman TB, Novak JM, Parkin TB, Karlen DL, Turco RF, Konopka AE (1994) Field-scale variability of soil properties in central Iowa soils. Soil Science Society Amsterdam Journal 58:1501\\u0026ndash;1511. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.2136/sssaj1994.03615995005800050033x\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eCao F, Jaunat J, Vergnaud-Ayraud V, Devau N, Labasque T, Guillou A, Ollivier P (2020) Heterogeneous behaviour of unconfined Chalk aquifers infer from combination of groundwater residence time, hydrochemistry and hydrodynamic tools. J Hydrol 581:124433. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.jhydrol.2019.124433\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003eCETESB. (2016). \\u003cem\\u003eQUALIDADE DAS \\u0026Aacute;GUAS SUBTERR\\u0026Acirc;NEAS NO ESTADO DE S\\u0026Atilde;O PAULO\\u003c/em\\u003e. [online]. Available at: \\u0026lt;https://cetesb.sp.gov.br/aguas-subterraneas/wp-content/uploads/sites/13/2013/11/Cetesb_QualidadeAguasSubterraneas2015_Web_20-07.pdf\\u0026gt;. [Accessed 29 July 2020].\\u003c/li\\u003e\\n \\u003cli\\u003eCPRH. (2020). \\u003cem id=\\\"isPasted\\\"\\u003eRECURSOS H\\u0026Iacute;DRICOS SUBTERR\\u0026Acirc;NEOS.\\u003c/em\\u003e [online]. Available at: \\u0026lt;http://www.cprh.pe.gov.br/downloads/25_Recursos_Hidricos_Subterraneos.pdf\\u0026gt;. [Accessed 22 July 2020].\\u003c/li\\u003e\\n \\u003cli\\u003eCPRM. (2007). \\u003cem id=\\\"isPasted\\\"\\u003eMapa de Dom\\u0026iacute;nios/Subdom\\u0026iacute;nios Hidrogeol\\u0026oacute;gicos do Brasil 1:2.500.000\\u003c/em\\u003e. [online]. Available at: \\u0026lt;http://www.cprm.gov.br/publique/Hidrologia/Mapas-e-Publicacoes/Mapa-de-Dominios%7CSubdominios-Hidrogeologicos-do-Brasil-1%3A2.500.000-632.html\\u0026gt;. [Accessed 1 July 2020].\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eDeutsch CV, Journel AG (1998) GSLIB Geostatistical Software Library and User\\u0026rsquo;s Guide, Second edn. Oxford University Press, New York, p 369\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eEngland EJ, Sparks A, Robinson MD (1989) Geo\\u0026mdash;EAS (Geostatistical Environmental Assessment Software). Environmental Software 4(2):70\\u0026ndash;75. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/0266-9838(89)90033-6\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eGilabert-Alarc\\u0026oacute;n C, Daessl\\u0026eacute; LW, Salgado-M\\u0026eacute;ndez SO, P\\u0026eacute;rez-Flores MA, Kn\\u0026ouml;ller K, Kretzschmar TG, Stumpp C (2018) Effects of reclaimed water discharge in the Maneadero coastal aquifer, Baja California, Mexico. Appl Geochem 92:121\\u0026ndash;139. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.apgeochem.2018.03.006\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eGois G, Delgado RC, de Oliveira J\\u0026uacute;nior JF (2015) Modelos te\\u0026oacute;ricos transitivos aplicados na interpola\\u0026ccedil;\\u0026atilde;o espacial do Standardized Precipitation Index (SPI) para os epis\\u0026oacute;dios de El Ni\\u0026ntilde;o forte no Estado do Tocantins, Brasil. Irriga 20(2):371\\u0026ndash;387. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.15809/irriga.2015v20n2p371\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eGOLDEN SOFTWARE (2010) Surfer for windows version 9.0. Golden, Colorado, p 66\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eGu X, Xiao Y, Yin S, Pan X, Niu Y, Shao J, Hao Q (2017) Natural and anthropogenic factors affecting the shallow groundwater quality in a typical irrigation area with reclaimed water, North China Plain. Environ Monit Assess 189(10):514. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1007/s10661-017-6229-3\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eGuadayol M, Cortina M, Guadayol JM, Caixach J (2016) Determination of dimethyl selenide and dimethyl sulphide compounds causing off-flavours in bottled mineral waters. Water Res 92:149\\u0026ndash;155. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.watres.2016.01.016\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eHossain M, Patra PK (2020) Contamination zoning and health risk assessment of trace elements in groundwater through geostatistical modelling. Ecotoxicol Environ Saf 189:110038. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.ecoenv.2019.110038\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eHu K, Awange JL, Forootan E, Goncalves RM, Fleming K (2017) Hydrogeological characterisation of groundwater over Brazil using remotely sensed and model products. Sci Total Environ 599:372\\u0026ndash;386. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.scitotenv.2017.04.188\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eHuang L, Liu Y (2017) Health information and consumer learning in the bottled water market. Int J Ind Organ 55:1\\u0026ndash;24. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.ijindorg.2017.08.002\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eIBGE (2018) Levantamento hidrogeologico. [online]. Available at: \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003evN7H6tMg3R_3AI01ndqVTgh3mr000\\u0026gt;. [Accessed 5 July 2020].\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eIBGE (2020) Hidroqu\\u0026iacute;mica subterr\\u0026acirc;nea. [online]. Available at:. [Accessed 5 July 2020]\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eIqbal J, Nazzal Y, Howari F, Xavier C, Yousef A (2018) Hydrochemical processes determining the groundwater quality for irrigation use in an arid environment: The case of Liwa Aquifer, Abu Dhabi, United Arab Emirates. Groundwater for Sustainable Development 7:212\\u0026ndash;219. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.gsd.2018.06.004\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eJang CS (2015) Geostatistical analysis for spatially characterizing hydrochemical features of springs in Taiwan. Environ Earth Sci 73(11):7517\\u0026ndash;7531. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1007/s12665-014-3924-z\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eKaiser HF (1958) The varimax criterion for analytic rotation in factor analysis. Psychometrika 23(3):187\\u0026ndash;200. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1007/BF02289233\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eKunrath CCN, Patroc\\u0026iacute;nio DC, Rodrigues MAS, Benvenuti T, Amado FDR (2020) Electrodialysis reversal as an alternative treatment for producing drinking water from brackish river water: A case study in the dry season, northeastern Brazil. Journal of Environmental Chemical Engineering 8(2):103719. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.jece.2020.103719\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eLiu J, Gao Z, Wang M, Li Y, Ma Y, Shi M, Zhang H (2018) Study on the dynamic characteristics of groundwater in the valley plain of Lhasa City. Environ Earth Sci 77(18):646. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1007/s12665-018-7833-4\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003eMinist\\u0026eacute;rio da Sa\\u0026uacute;de. (1999).\\u0026nbsp;\\u003cem id=\\\"isPasted\\\"\\u003ePORTARIA N\\u0026ordm; 26, DE 15 DE JANEIRO DE 1999\\u003c/em\\u003e. [online]. Available at: \\u0026lt;http://bvsms.saude.gov.br/bvs/saudelegis/anvisa/1999/prt0026_15_01_1999.html\\u0026gt;. [Accessed 29 July 2020].\\u003c/li\\u003e\\n \\u003cli\\u003eMinist\\u0026eacute;rio da Sa\\u0026uacute;de. (2006). \\u003cem id=\\\"isPasted\\\"\\u003eVIGIL\\u0026Acirc;NCIA E CONTROLE DA QUALIDADE DA \\u0026Aacute;GUA PARA CONSUMO HUMANO\\u003c/em\\u003e. [online]. Available at: \\u0026lt;https://bvsms.saude.gov.br/bvs/publicacoes/vigilancia_controle_qualidade_agua.pdf\\u0026gt;. [Accessed 29 July 2020].\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eNunes LGDP, Oliveira MDV, de Souza AA, Lopes LDF, Dias PCES, Nogueira GB, Souza MAAD (2018) Water quality comparison between a supply network and household reservoirs in one of the oldest cities in Brazil. Int J Environ Health Res 29(2):173\\u0026ndash;180. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1080/09603123.2018.1531114\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003ePacheco MH, Esmerino EA, Capobiango CS, Cruz AG, Leddomado LS, Pimentel TC, de Freitas MQ (2018) Bottled mineral water: classic and temporal descriptive sensory analysis associated with liking. British Food Journal 120(7):1547\\u0026ndash;1560. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1108/BFJ-11-2017-0655\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003ePazand K, Khosravi D, Ghaderi MR, Rezvanianzadeh MR (2018) Identification of the hydrogeochemical processes and assessment of groundwater in a semi-arid region using major ion chemistry: A case study of Ardestan basin in Central Iran. Groundwater for Sustainable Development 6:245\\u0026ndash;254. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.gsd.2018.01.008\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003ePiper AM (1944) A graphic procedure in the geochemical interpretation of water-analyses. Eos, Transactions American Geophysical Union 25(6):914\\u0026ndash;928. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1029/TR025i006p00914\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eR Core Team (2019) R: A language and environment for statistical computing. R Foundation for Statistical Computing\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eSargazi S, Mokhtari M, Ehrampoush MH, Almodaresi SA, Sargazi H, Sarhadi M (2020) The application of geographical information system (GIS) approach for assessment of groundwater quality of Zahedan city, Sistan and Baluchestan Province, Iran. Groundwater for Sustainable Development 100509. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.gsd.2020.100509\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eSchwanz TG, Llorca M, Farr\\u0026eacute; M, Barcel\\u0026oacute; D (2016) Perfluoroalkyl substances assessment in drinking waters from Brazil, France and Spain. Sci Total Environ 539:143\\u0026ndash;152. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.scitotenv.2015.08.034\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eShruti VC, P\\u0026eacute;rez-Guevara F, Kutralam-Muniasamy G (2020) Metro station free drinking water fountain-A potential \\u0026ldquo;microplastics hotspot\\u0026rdquo; for human consumption. Environ Pollut 261:114227. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.envpol.2020.114227\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eSilva TA, Ros\\u0026aacute;rio FF, Silva J\\u0026uacute;nior GC (2020) Environmental isotopes and hydrochemical tracers applied to hydrogeological conceptual modeling of the southwest portion of the Amazon Aquifer System (Acre, Brazil). Appl Geochem 120:104670. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.apgeochem.2020.104670\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eTeramoto EH, Gon\\u0026ccedil;alves RD, Chang HK (2020) Hydrochemistry of the Guarani Aquifer System modulated by mixing with underlying and overlying hydrostratigraphic units. Journal of Hydrology: Regional Studies 30:100713. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.ejrh.2020.100713\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eTiwari K, Goyal R, Sarkar A (2017) GIS-based spatial distribution of groundwater quality and regional suitability evaluation for drinking water. Environmental Processes 4(3):645\\u0026ndash;662. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1007/s40710-017-0257-4\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eTlili-Zrelli B, Gueddari M, Bouhlila R, Naceur M, O (2015) Groundwater hydrogeochemistry of Mateur alluvial aquifer (Northern Tunisia). J Hydrogeol Hydrol Eng 5(1). \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttp://dx.doi.org/10.4172/2325-9647.1000128\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eTundisi JG (2008) Water resources in the future: problems and solutions. Estudos avan\\u0026ccedil;ados 22(63):7\\u0026ndash;16. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1590/S0103-40142008000200002\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eVauclin M, Vieira SR, Vachaud G, Nielsen DR (1983) The Use of Cokriging with Limited Field Soil Observations 1. Soil Sci Soc Am J 47(2):175\\u0026ndash;184. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003e10.2136/sssaj1983.03615995004700020001x\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eWarrick AW, Nielsen DR (1980) Spatial variability of soil physical properties in the field.. In: In: HILLEL D (ed) Applications of Soil Physics, vol 2. Academic, New York, pp 319\\u0026ndash;344\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eYang Q, Li Z, Ma H, Wang L, Mart\\u0026iacute;n JD (2016) Identification of the hydrogeochemical processes and assessment of groundwater quality using classic integrated geochemical methods in the Southeastern part of Ordos basin, China. Environ Pollut 218:879\\u0026ndash;888. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.envpol.2016.08.017\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eYang Q, Li Z, Ma H, Wang L, Mart\\u0026iacute;n JD (2016) Identification of the hydrogeochemical processes and assessment of groundwater quality using classic integrated geochemical methods in the Southeastern part of Ordos basin, China. Environ Pollut 218:879\\u0026ndash;888. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1016/j.envpol.2016.08.017\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eYang Q, Li Z, Xie C, Liang J, Ma H (2020) Risk assessment of groundwater hydrochemistry for irrigation suitability in Ordos Basin, China. Nat Hazards 101:309\\u0026ndash;325. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1007/s11069-018-3451-4\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eYin S, Xiao Y, Gu X, Hao Q, Liu H, Hao Z, Pei Q (2019) Geostatistical analysis of hydrochemical variations and nitrate pollution causes of groundwater in an alluvial fan plain. Acta Geophys 67(4):1191\\u0026ndash;1203. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1007/s11600-019-00302-5\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n \\u003cli\\u003e\\u003cspan\\u003eYin, S., Xiao, Y., Gu, X., Hao, Q., Liu, H., Hao, Z., ... \\u0026amp; Pei, Q. (2019). Geostatistical analysis of hydrochemical variations and nitrate pollution causes of groundwater in an alluvial fan plain. \\u003cem\\u003eActa Geophysica, 67\\u003c/em\\u003e(4), 1191-1203. \\u003cspan class=\\\"ExternalRef\\\"\\u003e\\u003cspan class=\\\"RefSource\\\"\\u003ehttps://doi.org/10.1007/s11600-019-00302-5\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/span\\u003e\\u003c/li\\u003e\\n\\u003c/ol\\u003e\"}],\"fulltextSource\":\"\",\"fullText\":\"\",\"funders\":[],\"hasAdminPriorityOnWorkflow\":false,\"hasManuscriptDocX\":true,\"hasOptedInToPreprint\":true,\"hasPassedJournalQc\":\"\",\"hasAnyPriority\":false,\"hideJournal\":false,\"highlight\":\"\",\"institution\":\"\",\"isAcceptedByJournal\":false,\"isAuthorSuppliedPdf\":false,\"isDeskRejected\":\"\",\"isHiddenFromSearch\":false,\"isInQc\":false,\"isInWorkflow\":true,\"isPdf\":false,\"isPdfUpToDate\":true,\"isWithdrawnOrRetracted\":false,\"journal\":{\"display\":true,\"email\":\"info@researchsquare.com\",\"identity\":\"environmental-earth-sciences\",\"isNatureJournal\":false,\"hasQc\":true,\"allowDirectSubmit\":false,\"externalIdentity\":\"enge\",\"sideBox\":\"Learn more about [Environmental Earth Sciences](https://www.springer.com/journal/12665)\",\"snPcode\":\"12665\",\"submissionUrl\":\"https://submission.nature.com/new-submission/12665/3\",\"title\":\"Environmental Earth Sciences\",\"twitterHandle\":\"\",\"acdcEnabled\":true,\"dfaEnabled\":true,\"editorialSystem\":\"em\",\"reportingPortfolio\":\"Springer Hybrid\",\"inReviewEnabled\":true,\"inReviewRevisionsEnabled\":false},\"keywords\":\"mineral water, Beberibe aquifer, geostatistics, water quality\",\"lastPublishedDoi\":\"10.21203/rs.3.rs-1210837/v1\",\"lastPublishedDoiUrl\":\"https://doi.org/10.21203/rs.3.rs-1210837/v1\",\"license\":{\"name\":\"CC BY 4.0\",\"url\":\"https://creativecommons.org/licenses/by/4.0/\"},\"manuscriptAbstract\":\"\\u003cp\\u003eThe objective is to characterize the hydrochemical and qualitative composition of the groundwater of the Beberibe aquifer located in the Northeast of Brazil (NEB) using geostatistical, geoprocessing, and multivariate analysis techniques. The Beberibe aquifer is located in the state of Pernambuco (PE), NEB, and is responsible for supplying water to the northern portion of the Recife Metropolitan Region (RMRE). The data collected came from 34 georeferenced mining and mineral water production units. From this, a Piper diagram was constructed and hydrogeological and hydrochemical spatial maps were generated. The data were submitted to descriptive statistics, geostatistics, and principal component analysis (PCA). The results showed that the artesian wells near the coast have a porous geological formation. In all artesian wells, the waters are of the mixed bicarbonate type and highly potable for human consumption. The coefficient of variation was medium to high for all hydrochemical variables. The variables showed a better fit to the Gaussian model, in the study of spatial dependence. In PCA, Cl was the one with the highest correlation with nitrate (NO\\u003csub\\u003e3\\u003c/sub\\u003e). The anthropogenic influence on water quality at the center of RMRE, mainly on NO\\u003csub\\u003e3\\u003c/sub\\u003e concentrations, had a direct impact on water quality for human health.\\u003c/p\\u003e\",\"manuscriptTitle\":\"Statistical modeling of groundwater geochemistry in northeastern Brazil\",\"msid\":\"\",\"msnumber\":\"\",\"nonDraftVersions\":[{\"code\":1,\"date\":\"2022-02-11 15:50:29\",\"doi\":\"10.21203/rs.3.rs-1210837/v1\",\"editorialEvents\":[{\"type\":\"communityComments\",\"content\":0},{\"type\":\"editorInvitedReview\",\"content\":\"\",\"date\":\"2022-02-10T01:54:11+00:00\",\"index\":0,\"fulltext\":\"\"},{\"type\":\"reviewersInvited\",\"content\":\"\",\"date\":\"2022-02-09T13:59:53+00:00\",\"index\":\"\",\"fulltext\":\"\"},{\"type\":\"editorAssigned\",\"content\":\"\",\"date\":\"2021-12-31T13:12:07+00:00\",\"index\":\"\",\"fulltext\":\"\"},{\"type\":\"submitted\",\"content\":\"Environmental Earth Sciences\",\"date\":\"2021-12-28T07:46:32+00:00\",\"index\":\"\",\"fulltext\":\"\"}],\"status\":\"published\",\"journal\":{\"display\":true,\"email\":\"info@researchsquare.com\",\"identity\":\"environmental-earth-sciences\",\"isNatureJournal\":false,\"hasQc\":true,\"allowDirectSubmit\":false,\"externalIdentity\":\"enge\",\"sideBox\":\"Learn more about [Environmental Earth Sciences](https://www.springer.com/journal/12665)\",\"snPcode\":\"12665\",\"submissionUrl\":\"https://submission.nature.com/new-submission/12665/3\",\"title\":\"Environmental Earth Sciences\",\"twitterHandle\":\"\",\"acdcEnabled\":true,\"dfaEnabled\":true,\"editorialSystem\":\"em\",\"reportingPortfolio\":\"Springer Hybrid\",\"inReviewEnabled\":true,\"inReviewRevisionsEnabled\":false}}],\"origin\":\"\",\"ownerIdentity\":\"6d999557-18ea-4b76-aa5e-e6d577a15c71\",\"owner\":[],\"postedDate\":\"February 11th, 2022\",\"published\":true,\"recentEditorialEvents\":[],\"rejectedJournal\":[],\"revision\":\"\",\"amendment\":\"\",\"status\":\"under-review\",\"subjectAreas\":[],\"tags\":[],\"updatedAt\":\"2022-08-19T04:19:34+00:00\",\"versionOfRecord\":[],\"versionCreatedAt\":\"2022-02-11 15:50:29\",\"video\":\"\",\"vorDoi\":\"\",\"vorDoiUrl\":\"\",\"workflowStages\":[]},\"version\":\"v1\",\"identity\":\"rs-1210837\",\"journalConfig\":\"researchsquare\"},\"__N_SSP\":true},\"page\":\"/article/[identity]/[[...version]]\",\"query\":{\"redirect\":\"/article/rs-1210837\",\"identity\":\"rs-1210837\",\"version\":[\"v1\"]},\"buildId\":\"cBFmMYwuxLRRLfASyISRj\",\"isFallback\":false,\"isExperimentalCompile\":false,\"dynamicIds\":[84888],\"gssp\":true,\"scriptLoader\":[]}","source_license":"CC-BY-4.0","license_restricted":false}