Next Article in Journal
Daily Fluctuations in the Isotope and Elemental Composition of Tap Water in Ljubljana, Slovenia
Previous Article in Journal
Ionic Runoff as a Way to Determine the Degree of Karst Denudation (Case Study Jasov Plateau, Slovak Karst, Slovakia)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of Rainfall Seasonality in Groundwater Chemistry at Western Region of São Paulo State—Brazil

Department of Chemistry and Biochemistry, School of Science and Technology—Sao Paulo State University (UNESP), Rua Roberto Simonsen 305, Presidente Prudente, SP 19060-900, Brazil
*
Author to whom correspondence should be addressed.
Water 2021, 13(11), 1450; https://doi.org/10.3390/w13111450
Submission received: 18 February 2021 / Revised: 24 April 2021 / Accepted: 19 May 2021 / Published: 21 May 2021
(This article belongs to the Section Water Quality and Contamination)

Abstract

:
The present study evaluated the spatiotemporal variation in concentration of cadmium, lead and copper ions in groundwater wells in the stratigraphic subdivision “Santo Anastácio” that belongs to the Bauru aquifer system in the western region of São Paulo State. Exploratory statistics methods were employed to investigate the response of the concentration of these metals in the aquifer through the pluviometric index of the region. The results show a direct dependence of the mean monthly flow of the metals in the groundwaters to the monthly rainfall flow. The observed behavior was cyclic with a gradual increase and decrease in the flow throughout time. Two groups of cyclic variation were identified. The seasonality of the mean monthly flow of Cd2+ and Pb2+ was inversely proportional to the magnitude of the pluviometric index of the region studied. Meanwhile, the seasonality of Cu2+ was directly correlated to the seasonable rainfall variability. These behaviors lead us to point out that cadmium and lead come from minerals present in the aquifer itself and the presence of copper in groundwater is associated with an anthropogenic action due to the region’s agricultural activity. The study helps us better comprehend the behavior of the whole groundwater system through a comparison with temporal hydrogeochemistry.

1. Introduction

The exploitation of groundwaters has been taking great proportions due to its various advantages, such as the quality of the water, the costs of the exploitation and the reduction in the treatment phases until it reaches the final consumer. Its use for public supply has gained significant relevance in areas where surface water is scarce. The increase in its consumption is mainly due to populational growth, climate changes, which cause a shortage of surface water, and the increase of surface water pollution [1,2,3,4].
Groundwater reserves are a dynamic resource subjected to qualitative and quantitative modifications as a result of various contamination sources (e.g., agricultural, industrial and domestic) [5,6,7]. The quality and characteristics of groundwaters are a function of natural processes (e.g., geology, groundwater flow direction, quality of groundwater recharge/water and rock interactions), anthropogenic activities (e.g., agricultural production, industrial growth, urbanization with an increase in groundwater exploitation) and atmospheric input [4,8,9,10]. The index of potentially toxic elements discard has been intensified especially due to irregular discard in natural reservoirs through anthropogenic actions, such as mining, industrialization, irregular landfill constructions, indiscriminate usage of pesticides and fertilizers in agriculture, domestic sewage, etc. The presence of these potentially toxic elements is responsible for adverse effects in the environment, causing damage to public health and the economy [11]. Apart from anthropogenic actions, the insertion of potentially toxic metals in aquatic systems occurs naturally through geochemical processes like weathering, the region’s geological structure and climate changes, in addition to the interaction of atmospheric particulate matter in the rain composition and its interactions when there is the percolation of this water in the soil, mixing into the existing aquifers throughout the leaching course in soil [12,13].
Facing the enormous importance of groundwater previously mentioned, countless researches towards identification and quantification of potentially toxic metals in hydrographic basins and aquifers have been conducted [14,15]. The seasonality of rainfall affects dynamically the conditions of groundwater quality and quantity, resulting in alterations in the redox conditions, in the concentrations of substances in situ and the levels of water in the aquifer. Studies of seasonal identification and spatial variation of the anthropogenic and natural effects aiming to improve and comprehend hydrogeochemical processes based on pluviometric indexes have been described in the literature [16,17,18,19,20]. In this present paper, we present the identification of seasonable characteristics and spatial variations in the natural and anthropogenic effects aiming to improve and comprehend hydrochemical processes of groundwater in the western region of São Paulo State. The study was conducted for 17 months following the variations in concentrations of lead, cadmium and copper in groundwaters.

2. Methods

2.1. Study Area

The study area is located in the west of the São Paulo State—Brazil. With an area of approximately 17,500 km2 and an estimated population of about 540,000 inhabitants (demographic density of 30.76 inhabitants/km2), where the Presidente Prudente municipality stands out as a major regional economic pole, the economy is based on agriculture, livestock, industry and commerce. The area belongs to the geomorphological province called Planalto Ocidental Paulista and the relief is characterized by a succession of smoothed hills composed of sandstone spikes. According to the Köppen classification [21], the climate in the Presidente Prudente region is Aw: mesothermal with hot summers and dry winters, annual rainfall, with an average of 130 mm (driest month in July, with an annual average of 39 mm), average temperature of the warmest month (February) of 25.5 °C and that of the coldest month (June) of 20 °C.
The hydrogeological study area is based on the Bauru Aquifer System highlighted as a green area on the map of the São Paulo State (Figure 1), composed of the Marilia, Adamantina, Araçatuba, Santo Anastácio and Caiuá aquifers [22]. The system behaves as a hydrogeological unit of regional extension, continuous, free and locally confined. Groundwater was collected from the Santo Anastácio aquifer. The Santo Anastácio Aquifer emerges in a narrow strip parallel to the Paraná River (yellow highlighted in Figure 1), but on a subsurface it advances in an easterly direction, extending for approximately 67,000 km2.
The Santo Anastácio Formation consists predominantly of very fine to medium-grained sandstones, rounded to sub-angular grains, reddish-brown in color, ferruginous and locally carbonate cementation. It is poor in sedimentary structures, with massive strata with a maximum thickness of 80 m, with incipient plane-parallel stratification. Regarding the classification from the point of view of chemical composition, the aquifer has a predominance of calcium or magnesium bicarbonated water. The aquifer transmissivity is around 50 to 100 m2/day.

2.2. Groundwater Sampling and Analysis

Groundwater samples were collected from 7 monitoring wells over a period of 17 months from August 2017 to December 2018. The sampling points were identified by W1, W2, W3, W4, W5, W6 and W7 (Figure 2). Samples were collected every two weeks throughout the collection months. Table 1 presents the data and characteristics of the locations where the collections were performed (points 1 to 7).
Samples were collected after pumping based on the International Standard MSZ EN ISO 5667–1. In the field, temperature and pH were measured in a Metrohm pH meter. The concentrations of copper, lead and cadmium in water were determined by differential pulse stripping voltammetry, according to the Metrohm procedure (VA Application Note nº V-86 version 01). All measurements were made in triplicate, the reproducibility of the analytical data was 5%. Results were statistically characterized and analyzed using Excel 2016 and Origin 2019 statistical software (OriginLab Corporation, Northampton, MA, USA).

3. Results and Discussion

The employment of exploratory statistics aimed to estimate the degree of variability assuming that the variables studied are affected by the pluviometric index. Thus, the discussion in the present paper is distributed in topics to each biweekly analyzed chemical component in groundwater.

3.1. Spatiotemporal Distribution and Exploratory Statistic of Cadmium in Groundwater

From the data obtained in the analysis of Cd2+ concentration, a basic descriptive statistic was performed (minimum, maximum, mean, standard deviation and median values) to understand the relation between the main parameters and the sampled wells in the research field. The values of mean monthly cadmium flow, among all wells, variated from 1.48 to 44 µg L−1. Table 2 presents the mean monthly values and their respective standard deviations of the obtained concentrations in the samples, as well as the values of minimum, maximum, mean, standard deviation and median in each well during the period studied.
The influence of the spatiotemporal behavior of the aquifer hydraulic charge on the mean monthly cadmium concentration was correlated to the parameter of mean monthly pluviometric level of the study region as if it was the recharge level or water volume level of the aquifer [23]. A Figure 3A shows a seasonable or cyclic behavior of the increasing and decreasing in the mean monthly flow of the mentioned metallic cation throughout time. With a statistical technique to aid the data interpretation, a trend curve was elaborated (short-dash-dot curve of Figure 3A) for the cadmium concentration in the aquifer in a temporal function in which the analysis occurred.
Through the model, we can observe that the cadmium concentration exhibits a crescent tendency in the dry season (April through September) and a decrescent tendency in the rainy season (October through March). The transformation of the concentration data and pluviometric index to the logarithmic base (see Figure 3B) demonstrates correlative and significative cycles in both parameters analyzed in temporal sampling function. The model indicates that the phenomenon has a periodicity, or it obeys a periodical function, where the mean logarithmic concentration of cadmium is inversely proportional to the magnitude of the monthly seasonable variability of rainfall. The cyclic variation in the cadmium mean monthly flow is indicative that the chemical element does not come from an external source and it is found in the aquifer’s area itself. During the rainy season, there is a dilution effect due to the increase in the aquifer’s water volume. Meanwhile, during the dry season, the low water level in the aquifer causes the opposite process, when the cadmium concentration increases.
Figure 4 presents a boxplot diagram [24] that provides a better representation of the observed data variation in each well for the cadmium concentration during the dry (Figure 4A) and rainy seasons (Figure 4B). The distinction between the rainy and dry seasons is very clear as it is shown in both boxplot diagrams. A general seasonal distinction is based on the only components from the hydric balance on a monthly scale as shown in Figure 3B. Analyzing the dataset of each well during the dry season (Figure 4A), well W3 presents the highest dispersion on the values of cadmium concentration in the interquartile range, which consists in a difference between the third and the first quartile. One of the concentration variability factors is the process of groundwater pumping during the driest seasons, influenced by the extremely low levels of the groundwater, resonating significantly in the chemical compound’s concentration in the aquifer [16]. Comparatively, the lowest dispersion in the interquartile range for the values of cadmium concentration was observed in the well W1. A factor for the lowest dispersion in the concentration of cadmium may be related to the location of the well W1 that is situated in a green area with a low population density. Analyzing the median lines, the wells W5 and W6 present a positive asymmetric distribution, indicating that the median is close to the first quartile or that the median value is lower than the mean value. Meanwhile, the wells W2, W3, W4 and W7 present a negative asymmetric distribution. It was verified for the well W1 that the mean and median values are coincidental demonstrating a symmetric distribution. It is important to point out that the median is the central tendency measure more appropriate when the data present asymmetric distribution since the arithmetic mean is influenced by the extreme values.
Contrastingly, cadmium concentration is lower in groundwater collected during the rainy season (Figure 4B) with a mean monthly below 10 μg L−1. As previously shown in Figure 3B, this result is due to the dilution of the concentration by groundwater recharge. However, the descriptive parameters in the boxplots present anomalous values (outliers) in almost every well studied (except in wells W1 and W2). The anomalous values are related to the first rainy month (October) in which the cadmium concentrations are still high in comparison to the dataset for the period, implying that the dilution factor is in the initial phase in the aquifer due to the dependence of the groundwater recharge rate [25].
Table 3 presents mean values for the monthly rainfall index (MRI), pH and cadmium concentration. The monthly mean values of pH were calculated based on the pH measurements of groundwater in each well studied.
To verify the dependence degree of the mean monthly flow to the rainfall index and the pH of the water collected, Pearson correlation coefficient was applied with paired Student’s t-test with a significance level of 5% on the correlation coefficient obtained, where H0: r = 0 and H1: r ≠ 0. As shown in Table 4, there is a moderate relation between the MRI parameters and the mean cadmium concentrations in groundwater.
Comparing the calculated t value (−5.352) to the critical value for the Student’s t-distribution (±2.131), t value is out of the region for the H0 hypothesis to be accepted. Thus, we can conclude that there is enough evidence to correlate the concentration parameters of Cd(II) to MRI. Concerning the mean monthly of Cd(II) and the mean pH, the Pearson coefficient obtained shows that there is a weak relation. Comparing the calculated t value (−0.379) to the critical value for the Student’s t-distribution (±2.131), t value is in the region for the H0 hypothesis to be accepted. Thus, we can conclude that there is not enough evidence suggesting that there is a relation between the cadmium concentration and the pH. The statistical result obtained through Pearson’s correlation indicates that the variability in groundwater pH is not caused by the variation in cadmium concentration in function of the region’s rainfall index.

3.2. Spatiotemporal Distribution and Exploratory Statistic of Lead in Groundwater

Table 5 presents the mean monthly values and their respective standard deviations of the obtained samples concentrations, as well as minimum, maximum, mean, standard deviation and median values in each well for study period.
The mean monthly lead flow variation (minimum = 1.39 µg L−1 and maximum = 44.9 µg L−1) in the wells was similar to what was observed for the cadmium concentrations. The concentration magnitude varied with the monthly period of the sampling and the geographical position of the wells studied. With the obtained data of lead concentrations, a spatiotemporal graph for the variation in total mean lead concentration was plotted in the logarithmic base (Figure 5B).
The profile of lead concentration variation in groundwater followed a periodicity model in function of time and the highest concentration magnitudes were observed during the dry season. Analyzing comparatively the lead spatiotemporal behavior, we can conclude that this element is found in the aquifer itself. Its dilution in the groundwater is due to increasing in hydric volume and consequently to the recharge flow (pluviometric level). For better visualization of what was verified before, Table 6 presents the lead mean concentration in wells W1 to W7 with their respective standard deviations, in the rainy season (October through March) and in the dry season (April through September). Evaluating the boxplot diagram, the highest dispersions and lead concentration distributions were observed during the dry season (Figure 6A). As observed before, well W1 presented the lowest variability of concentration in the sampling temporal scale. In general, the medians are very close to the concentration means, indicating that the observed concentrations in the dry season are a symmetrical normal distribution. For the rainy season, the mean monthly lead flow values (Figure 6B) were below 10 μg L−1. However, outlier concentrations were observed in wells W5, W6 and W7 and these outlier values correspond precisely to the beginning of the rainy season.
In the Pearson test to analyze the dependence degree of Pb2+ concentration to the mean pH of the collected water, with a correlation coefficient (ρ) of −0.0066, there is an indication that these two variables are not correlated. To confirm that the variability in lead concentration is not influenced by the pH or vice versa, a two-tailed Student’s t-test with a significance level of 5% was applied and considering the null and true hypothesis as r = 0 and r ≠ 0, respectively. Comparing the calculated t value (−0.26) to the critical value for the Student’s t-distribution (±2.131), the calculated t value is in the region for the H0 hypothesis to be accepted. Thus, we can conclude that there is not enough evidence to support that Pb2+ concentration and pH are correlated.

3.3. Spatio-Temporal Distribution and Exploratory Statistic of Copper in Groundwater

In the copper concentration spatiotemporal distribution study, it was identified an inverted behavior when compared to the cadmium and lead in groundwater. The mean monthly copper flow values, among all wells, variated from 4.8 to 2.479 µg L−1 (see Table 7). In addition, the standard deviation values are higher than the mean, indicating that the copper concentrations are distributed around a wide range of values throughout the entire period.
The copper mean monthly concentration through time and its relation of its logarithmic base to the pluviometric index are represented in Figure 7A,B, respectively. Analyzing Figure 7A, it can be observed that the copper concentration in the aquifer increases significantly during the month with higher rainfall intensity. After its concentration peak, the reduction of copper presence in groundwater is gradual and concomitantly with the pluviometric index. Its residence time or memory effect in the aquifer is around 5 months approximately. Through Figure 7B it is possible to observe that the copper concentration in groundwater is directly proportional to the pluviometric level. This behavior indicates strongly that the greater source of copper in groundwater is originated out of the aquifer. It is known that the recharge time can be around days to years depending directly on the hydrogeology properties of the aquifer and of the levels of direct recharge areas. Through the copper variability characterization in groundwater (see Figure 7B), the present study points out that the response of the aquifer’s recharge in function of rainfall regime is between 1–2 months, given that the magnitude of the maximum copper concentration is reached after the first period of rainfall. This result corroborates studies of the water table monitoring facing correlations to the pluviometric level in the region [26].
Applying the boxplot diagram (Figure 8) to the copper concentration variability for both pluviometric periods, it is observed more clearly the pluviometric influence when the higher copper concentrations happen during the rainy season (October through March) and its reduction during the dry season (April through September). Differently from the data obtained for cadmium and lead during the dry season, all wells presented outliers concentrations, which contributed to a higher asymmetric distribution of the data (Figure 8A). This adverse characteristic is related to the beginning of the dry season (April, May and June) when the copper concentration is still relatively high in the aquifer. Through Figure 8B that represents the rainy season, the concentrations magnitudes are about 5 times higher in comparison to the dry season mean. From this observation, it is presented in the item “Conclusion” an analysis of the possible causes of the increase in its concentration during the rainy season. Analyzing the dataset, wells W2, W5, W6 and W7 are the ones that present the higher concentration variability during the rainy season. With the mean monthly pluviometric index values, pH and Cu2+ concentration, it was verified the dependence degree of these parameters through Pearson’s test. As shown in Table 8, there is a strong correlation between the parameters MRI and copper mean concentrations in groundwater. This conclusion is confirmed through Student’s t-test with a significance level of 5%, where H0: r = 0 and H1: r ≠ 0. Comparing the calculated t value (8.743) to the critical value for the Student’s t-distribution (±2.131), the calculated t value is out of the region for the H0 hypothesis to be accepted. Thus, we can conclude that there is enough evidence that the parameters of Cu2+ concentration and mean monthly pluviometric index are correlated.
As for the relation between Cu2+ concentration and the groundwater pH, the Pearson coefficient value obtained shows that there is a moderated correlation. This conclusion is confirmed by the obtained t value (3.486) being higher than the critical value for the Student’s t-distribution (±2.131), indicating that the calculated t value is out of the region for the H0 hypothesis to be accepted. Thus, we can conclude that the pH of the aquifer is influenced by the presence and concentration of Cu2+.

3.4. Principal Component Analysis

A Principal Component Analysis (PCA) was employed to determine the discriminant functions in order to confirm the spatiotemporal variations of chemical elements in groundwater. Table 9 presents two principal components (PC) that obtained the most significant percentage of variability (totalizing 98.3% of the data variability) and the vectors’ values per parameter.
The criterion employed in this kind of analysis was the evaluation of the weight of the variable associated with a component. Weights above 0.7 are indicative of strong association; for values between 0.5–0.7, the variable is considered to be moderately associated; weights inferior to 0.5, the variable has a weak association to the component [27]. The first component (82.2%) is associated with a positive moderate variance for cadmium and lead during the dry season. In the same component, the copper variability presented a negative and moderate association for the same period. The second component (16.1%) was strongly associated solely with copper, which can be interpreted as the influence of the rainy season. According to a visual interpretation of the bidimensional representation of the two first principal components (Figure 9), the influence vectors are identified into two groups, as evidenced previously in the exploratory statistics analysis. The variable of copper concentration was close to CP2 axis, showing that the months of January and February (months of the highest pluviometric levels) present high correlations.
Contrastingly, cadmium and lead are more representatives to the CP1 axis, where we can observe that the vectors are equivalent and that the concentrations of these elements are correlative to the months of lowest pluviometric index (July through September). The acute angle of the vectors’ points to a high correlation between the two variables.

4. Conclusions

The results obtained through exploratory statistics in cadmium, lead and copper concentration variability provided some information about the processes that control the chemistry of groundwater in the study area. Cadmium and lead ions showed high concentrations in the low pluviometric index period, which can be reintroduced in groundwater through natural interaction water-mineral. During the dry season, the water table becomes lower and more constant [28], causing the phenomenon of concentration of these elements. In our research group, we have been studying the emanation of radon gas and its correlation to lead in soil in the region [29]. The presence of radon in the subsoil is an indication of alkaline igneous rocks and feldspathoids that present minerals with the content of uranium and that the lead in the aquifer is originated from radioisotopes generated in the series of radioactive decay of uranium.
The seasonality of copper in the aquifer was different and contrary to the results previously observed. The copper concentration increased exponentially in months that have the highest pluviometric level in the region and it decreases gradually during the months that follow the maximum peak in concentration. To analyze possible sources of copper in groundwater, reference materials about aquifer contamination correlated to economic activities and the land usage and occupation were consulted. Due to the observed characteristics, copper comes from a diffuse source originated by agricultural activity. The diffuse sources are characterized by presenting a wide area of contribution, deriving from activities that deposit pollutants sparsely, reaching bodies of water only in an intermittent way, associated with rainy periods [30,31]. The increment in copper concentration in groundwaters may have an origin in the leaching of agricultural chemicals through rainfall. The main economic base of the region is agribusiness, where it is most prominent the cultivation of sugar cane and yam. Agricultural inputs or subproducts employed with a corrective or nutritional purposed for both sugar cane and yam crops are copper salts-based [32]. Since we have no evidence of the contribution of the chemical elements of rain, further research should be done in the future to determine the composition of rain in the region and understand the contribution seasonal to underground water.

Author Contributions

Investigation, writing—original draft, L.A.M.; investigation, visualization, M.R.M.-S.; data curation, re-writing, statistical analysis, J.J.T.; conceptualization, methodology, validation, writing—review and editing, supervision, project administration, funding acquisition, M.F.S.T. All authors have read and agreed to the published version of the manuscript.

Funding

The authors acknowledge CEPID-FAPESP (2013/07296-2) and CNPq (301298/2017-3) for their financial support.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author: MFST, upon request.

Acknowledgments

The authors would like to acknowledge SJT.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cavazzana, G.H.; Lastoria, G.; Gabas, S.G. Surface-groundwater interaction in unconfined sedimentary aquifer system in the Brazil’s tropical wet region. RBRH-Rev. Bras. Recur. 2019, 24. [Google Scholar] [CrossRef] [Green Version]
  2. Machiwal, D.; Jha, M.K. Identifying sources of groundwater contamination in a hard-rock aquifer system using multivariate statistical analyses and GIS-based geostatistical modeling techniques. J. Hydrol. Reg. Stud. 2015, 4, 80–110. [Google Scholar] [CrossRef] [Green Version]
  3. Guo, X.R.; Zuo, R.; Meng, L.; Wang, J.S.; Teng, Y.G.; Liu, X.; Chen, M.H. Seasonal and Spatial Variability of Anthropogenic and Natural Factors Influencing Groundwater Quality Based on Source Apportionment. Int. J. Environ. Res. Public Health 2018, 15, 279. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Jayalakshmi, B.; Ramachandramoorthy, T.; Paulraj, A. Statistical interpretation on seasonal variations of groundwater quality in Ramanathapuram coastal tract, Tamil Nadu, India. Environ. Earth Sci. 2014, 72, 1271–1278. [Google Scholar] [CrossRef]
  5. Dubrovsky, N.M.; Burow, K.R.; Gronberg, J.A.M. Effects of two contrasting agricultural land uses on shallow groundwater quality in the San Joaquin Valley, California: Design and preliminary interpretation. In Models for Assessing and Monitoring Groundwater Quality; IAHS Press: Wallingford, UK, 1995; pp. 49–58. [Google Scholar]
  6. Frapporti, G.; Vriend, S.P.; Vangaans, P.F.M. Hydrogeochemistry of the Shallow Dutch Groundwater—Interpretation of the National Groundwater Quality Monitoring Network. Water Resour. Res. 1993, 29, 2993–3004. [Google Scholar] [CrossRef]
  7. Adimalla, N.; Li, P.Y.; Venkatayogi, S. Hydrogeochemical Evaluation of Groundwater Quality for Drinking and Irrigation Purposes and Integrated Interpretation with Water Quality Index Studies. Environ. Process. 2018, 5, 363–383. [Google Scholar] [CrossRef]
  8. Swamee, P.K.; Misra, A.K.; Poyil, R.P.; Choudhary, L.; Kaur, A.; Chintalapati, S. A linear and nonlinear approach for interpretation of groundwater quality based on water quality index. Environ. Dev. Sustain. 2018, 20, 2623–2642. [Google Scholar] [CrossRef]
  9. Huntington, T.G. Evidence for intensification of the global water cycle: Review and synthesis. J. Hydrol. 2006, 319, 83–95. [Google Scholar] [CrossRef]
  10. Oki, T.; Kanae, S. Global hydrological cycles and world water resources. Science 2006, 313, 1068–1072. [Google Scholar] [CrossRef] [Green Version]
  11. Miyazawa, M.; Giminez, S.M.N.; Yabe, M.J.S.; Oliveira, E.L.; Kamogawa, M.Y. Absorption and toxicity of copper and zinc in bean plants cultivated in soil treated with chicken manure. Water Air Soil Pollut. 2002, 138, 211–222. [Google Scholar] [CrossRef]
  12. Jalali, M. Salinization of groundwater in arid and semi-arid zones: An example from Tajarak, western Iran. Environ. Geol. 2007, 52, 1133–1149. [Google Scholar] [CrossRef]
  13. Subramani, T.; Rajmohan, N.; Elango, L. Groundwater geochemistry and identification of hydrogeochemical processes in a hard rock region, Southern India. Environ. Monit. Assess. 2010, 162, 123–137. [Google Scholar] [CrossRef]
  14. Nriagu, J.O. A Silent Epidemic of Environmental Metal Poisoning. Environ. Pollut. 1988, 50, 139–161. [Google Scholar] [CrossRef] [Green Version]
  15. Bodo, B.A. Heavy-Metals in Water and Suspended Particulates from an Urban Basin Impacting Lake-Ontario. Sci. Total Environ. 1989, 87–88, 329–344. [Google Scholar] [CrossRef]
  16. Sakakibara, K.; Tsujimura, M.; Song, X.F.; Zhang, J. Spatiotemporal variation of the surface water effect on the groundwater recharge in a low-precipitation region: Application of the multi-tracer approach to the Taihang Mountains, North China. J. Hydrol. 2017, 545, 132–144. [Google Scholar] [CrossRef]
  17. Rodriguez-Galiano, V.; Mendes, M.P.; Garcia-Soldado, M.J.; Chica-Olmo, M.; Ribeiro, L. Predictive modeling of groundwater nitrate pollution using Random Forest and multisource variables related to intrinsic and specific vulnerability: A case study in an agricultural setting (Southern Spain). Sci. Total Environ. 2014, 476, 189–206. [Google Scholar] [CrossRef]
  18. Baez-Cazull, S.E.; McGuire, J.T.; Cozzarelli, I.M.; Voytek, M.A. Determination of dominant biogeochemical processes in a contaminated aquifer-wetland system using multivariate statistical analysis. J. Environ. Qual. 2008, 37, 30–46. [Google Scholar] [CrossRef] [Green Version]
  19. Zheng, Y.; Chen, S.D.; Qin, H.P.; Jiao, J.J. Modeling the Spatial and Seasonal Variations of Groundwater Head in an Urbanized Area under Low Impact Development. Water 2018, 10, 803. [Google Scholar] [CrossRef] [Green Version]
  20. Shamsudduha, M.; Chandler, R.E.; Taylor, R.G.; Ahmed, K.M. Recent trends in groundwater levels in a highly seasonal hydrological system: The Ganges-Brahmaputra-Meghna Delta. Hydrol. Earth Syst. Sci. 2009, 13, 2373–2385. [Google Scholar] [CrossRef] [Green Version]
  21. Dubreuil, V.; Fante, K.P.; Planchon, O.; Neto, J.L.S. The types of annual climates in Brazil: An application of the classification of Koppen from 1961 to 2015. Confins 2018, 37. [Google Scholar] [CrossRef]
  22. Stradioto, M.R.; Teramoto, E.H.; Chang, H.K. Rock-solute reaction mass balance of water flowing within an aquifer system with geochemical stratification. Appl. Geochem. 2020, 123, 104784. [Google Scholar] [CrossRef]
  23. De Vries, J.J.; Simmers, I. Groundwater recharge: An overview of processes and challenges. Hydrogeol. J. 2002, 10, 5–17. [Google Scholar] [CrossRef]
  24. McBride, G.B. Using Statistical Methods for Water Quality Management: Issues, Problems and Solutions; John Wiley & Sons Inc.: Chichester, UK, 2005; pp. 1–315. [Google Scholar] [CrossRef]
  25. Wijngaard, R.R.; van der Perk, M.; van der Grift, B.; de Nijs, T.C.M.; Bierkens, M.F.P. The Impact of Climate Change on Metal Transport in a Lowland Catchment. Water Air Soil Pollut. 2017, 228, 107. [Google Scholar] [CrossRef] [Green Version]
  26. Goncalves, V.F.M.; Manzione, R.L. Groundwater Recharge Estimates at Bauru Aquifer System (Bas). Geo Uerj 2019. [Google Scholar] [CrossRef]
  27. Liu, C.W.; Lin, K.H.; Kuo, Y.M. Application of factor analysis in the assessment of groundwater quality in a blackfoot disease area in Taiwan. Sci. Total Environ. 2003, 313, 77–89. [Google Scholar] [CrossRef]
  28. Jiao, Y.F.; Liu, J.; Li, C.Z.; Qiu, Q.T.; Yu, X.Z.; Wang, W.; Zhang, G.J. Regulation storage calculation and capacity evaluation of the underground reservoir in the groundwater overdraft area of northern China. Environ. Earth Sci. 2020, 79, 18. [Google Scholar] [CrossRef]
  29. Moreira-Silva, M.R.; Saenz, C.A.T.; Nunes, J.O.R.; Godoy, M.; Teixeira, M.F.S. Evidence for a correlation between total lead concentrations in soils and the presence of geological faults. Environ. Chem. Lett. 2017, 15, 481–488. [Google Scholar] [CrossRef]
  30. Bendoricchio, G. Integrated management of water quality—The role of agricultural diffuse pollution sources. Water Sci. Technol. 1999, 39, vii. [Google Scholar]
  31. Albek, E.A.; Goncu, S.; Uygun, B.S.; Albek, M.; Avdan, Z.Y.; Gungor, O. Water quality monitoring with emphasis on estimation of point and diffuse pollution sources. Glob. Nest J. 2019, 21, 163–171. [Google Scholar] [CrossRef]
  32. Corbi, J.J.; Costa, C.G.; Gorni, G.R.; Colombo, V.; Rios, L. Environmental diagnosis of metals in streams near sugarcane cultivation areas: Current and historical analysis in the central region of the State of Sao Paulo. An. Acad. Bras. Cienc. 2018, 90, 2711–2719. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Hydrogeologic map of the Bauru Group in São Paulo state showing its outcropping lithostratigraphic units which correspond to homologous hydrogeologic units in the study area (modified from Ref. [22]). Highlighting the saturation and transmissivity of the Bauru Aquifer System in the portion of the study area.
Figure 1. Hydrogeologic map of the Bauru Group in São Paulo state showing its outcropping lithostratigraphic units which correspond to homologous hydrogeologic units in the study area (modified from Ref. [22]). Highlighting the saturation and transmissivity of the Bauru Aquifer System in the portion of the study area.
Water 13 01450 g001
Figure 2. Location of sampling wells in the study area.
Figure 2. Location of sampling wells in the study area.
Water 13 01450 g002
Figure 3. (A) Average temporal evolution of average cadmium concentration (n = 7) and the pluviometric levels (n = 39) in accordance with the seasonal period. The red short-dash-dot curve corresponds to the trend curve of the variability of cadmium concentration in the aquifer during the 17 months of study. (B) Normalized flow of cadmium concentration (black circle) and pluviometric index (red circle) throughout the sampling period. The results are presented as the mean value of 3 point.
Figure 3. (A) Average temporal evolution of average cadmium concentration (n = 7) and the pluviometric levels (n = 39) in accordance with the seasonal period. The red short-dash-dot curve corresponds to the trend curve of the variability of cadmium concentration in the aquifer during the 17 months of study. (B) Normalized flow of cadmium concentration (black circle) and pluviometric index (red circle) throughout the sampling period. The results are presented as the mean value of 3 point.
Water 13 01450 g003
Figure 4. Boxplot diagram for cadmium concentration in 7 wells during the 17 months of study. (A) Dry season and (B) Rainy season. Lines and crosses within boxes represent the median and mean, respectively. The asterisk corresponds to outlier (atypical values).
Figure 4. Boxplot diagram for cadmium concentration in 7 wells during the 17 months of study. (A) Dry season and (B) Rainy season. Lines and crosses within boxes represent the median and mean, respectively. The asterisk corresponds to outlier (atypical values).
Water 13 01450 g004
Figure 5. (A) Average temporal evolution of average lead concentration (n = 7) and the pluviometric levels (n = 39) in accordance with the seasonal period. The red short-dash-dot curve corresponds to the trend curve of the variability of cadmium concentration in the aquifer during the 17 months of study. (B) Normalized flow of lead concentration (black circle) and pluviometric index (red circle) throughout the sampling period.
Figure 5. (A) Average temporal evolution of average lead concentration (n = 7) and the pluviometric levels (n = 39) in accordance with the seasonal period. The red short-dash-dot curve corresponds to the trend curve of the variability of cadmium concentration in the aquifer during the 17 months of study. (B) Normalized flow of lead concentration (black circle) and pluviometric index (red circle) throughout the sampling period.
Water 13 01450 g005
Figure 6. Boxplot diagram for lead concentration in 7 wells during the 17 months of study. (A) Dry season and (B) Rainy season. Lines and crosses within boxes represent the median and mean, respectively. The asterisk corresponds to outlier (atypical values).
Figure 6. Boxplot diagram for lead concentration in 7 wells during the 17 months of study. (A) Dry season and (B) Rainy season. Lines and crosses within boxes represent the median and mean, respectively. The asterisk corresponds to outlier (atypical values).
Water 13 01450 g006
Figure 7. (A) Average temporal evolution of average copper concentration (n = 7) and the pluviometric levels (n = 39) in accordance with the seasonal period. The red short-dash-dot curve corresponds to the trend curve of the variability of cadmium concentration in the aquifer during the 17 months of study. (B) Normalized flow of copper concentration (black circle) and pluviometric index (red circle) throughout the sampling period. The results are presented as the mean value of 3 point.
Figure 7. (A) Average temporal evolution of average copper concentration (n = 7) and the pluviometric levels (n = 39) in accordance with the seasonal period. The red short-dash-dot curve corresponds to the trend curve of the variability of cadmium concentration in the aquifer during the 17 months of study. (B) Normalized flow of copper concentration (black circle) and pluviometric index (red circle) throughout the sampling period. The results are presented as the mean value of 3 point.
Water 13 01450 g007
Figure 8. Boxplot diagram for copper concentration in 7 wells during the 17 months of study. (A) Dry season and (B) Rainy season. Lines and crosses within boxes represent the median and mean, respectively. The asterisk corresponds to outlier (atypical values).
Figure 8. Boxplot diagram for copper concentration in 7 wells during the 17 months of study. (A) Dry season and (B) Rainy season. Lines and crosses within boxes represent the median and mean, respectively. The asterisk corresponds to outlier (atypical values).
Water 13 01450 g008
Figure 9. Graph of PC1 and PC2 weights for the three observed variables (Cu2+; Cd2+; Pb2+).
Figure 9. Graph of PC1 and PC2 weights for the three observed variables (Cu2+; Cd2+; Pb2+).
Water 13 01450 g009
Table 1. Characteristics of the areas where were acquired as sampling.
Table 1. Characteristics of the areas where were acquired as sampling.
PointW1W2W3W4W5W6W7
22°08′32.98″ S22°9′0.14″ S22°05′47.72″ S22°07′11.62″ S22°07′20.57″ S22°07′21.18″ S22°07′31.86″ S
51°27′49.79″ W51°25′23.60″ W51°25′42.95″ W51°24′50.49″ W51°24′30.98″ W51°24′29.77″ W51°24′42.25″ W
Altitude (m)435420449408446445422
Depth (m)120110175195170190172
UseIrrigationHuman consumptionHuman consumptionRecreationMaintenance and sanitationIrrigationPublic supply
Local aspectGreen area. Low population flow.Residential area. High vehicular and population flow.Green area.Residential area. High vehicular and population flow.Green area. High vehicular and population flow.Green area. High vehicular and population flow.Urbanized area.
Table 2. Mean monthly cadmium concentration (μg L−1) with standard deviation value determined for each studied well (W) during the 17 months.
Table 2. Mean monthly cadmium concentration (μg L−1) with standard deviation value determined for each studied well (W) during the 17 months.
MonthW1W2W3W4W5W6W7
Aug/XX12.7 ± 2.75.40 ± 3.2938.5 ± 1.224.1 ± 5.610.2 ± 0.315.3 ± 4.419.9 ± 5.2
Sep/XX9.93 ± 3.8418.8 ± 1.244.0 ± 6.426.6 ± 0.833.9 ± 5.630.1 ± 8.323.8 ± 1.1
Oct/XX6.63 ± 0.8310.2 ± 2.915.6 ± 4.023.5 ± 3.717.6 ± 3.711.8 ± 4.67.58 ± 2.99
Nov/XX4.48 ± 0.777.55 ± 0.393.24 ± 0.728.20 ± 0.375.76 ± 1.656.34 ± 0.753.48 ± 0.09
Dec/XX2.03 ± 1.402.22 ± 1.462.39 ± 0.147.09 ± 0.512.63 ± 1.745.33 ± 0.351.48 ± 0.16
Jan/YY1.75 ± 0.291.57 ± 0.511.87 ± 0.292.53 ± 0.733.74 ± 0.141.52 ± 0.151.99 ± 0.31
Feb/YY3.78 ± 0.172.80 ± 0.092.54 ± 0.453.41 ± 0.483.94 ± 0.112.18 ± 0.243.02 ± 0.07
Mar/YY4.71 ± 0.274.25 ± 0.673.22 ± 0.093.94 ± 0.065.46 ± 0.373.55 ± 0.045.08 ± 0.22
Apr/YY6.05 ± 0.535.65 ± 0.104.15 ± 0.675.06 ± 0.956.99 ± 0.023.97 ± 0.266.05 ± 0.40
May/YY8.80 ± 0.578.16 ± 0.356.03 ± 0.937.44 ± 1.198.55 ± 1.137.10 ± 0.308.10 ± 1.19
Jun/YY11.4 ± 0.19.81 ± 0.608.27 ± 0.3412.2 ± 4.411.8 ± 1.010.1 ± 0.110.9 ± 1.8
Jul/YY13.9 ± 0.212.7 ± 2.725.9 ± 3.818.9 ± 0.214.1 ± 0.313.6 ± 3.420.2 ± 2.4
Aug/YY17.2 ± 7.013.11 ± 9.133.3 ± 6.522.5 ± 3.716.3 ± 7.419.6 ± 6.424.1 ± 6.6
Sep/YY15.7 ± 5.915.7 ± 3.640.3 ± 5.622.7 ± 4.628.9 ± 7.532.0 ± 6.226.9 ± 6.0
Oct/YY7.10 ± 1.878.45 ± 2.7419.5 ± 5.620.2 ± 4.614.05 ± 4.914.6 ± 4.79.71 ± 4.0
Nov/YY3.98 ± 0.736.01 ± 1.845.72 ± 2.937.87 ± 0.494.50 ± 1.796.05 ± 1.104.91 ± 1.71
Dec/YY1.76 ± 0.872.12 ± 1.174.74 ± 2.734.65 ± 2.952.51 ± 1.013.99 ± 1.843.40 ± 2.23
Average *7.76 ± 4.967.91 ± 4.9915.2 ± 15.113.0 ± 8.711.2 ± 9.011.0 ± 9.110.6 ± 8.7
Minimum1.751.571.872.532.511.521.48
Median6.637.556.038.208.557.107.58
Maximum17.218.844.026.633.932.026.9
* Average value per well.
Table 3. Monthly rainfall index (MRI), mean monthly cadmium(II) concentration and mean monthly pH values determined by the average of all studied wells.
Table 3. Monthly rainfall index (MRI), mean monthly cadmium(II) concentration and mean monthly pH values determined by the average of all studied wells.
MonthMRI
(mm)
Mean [Cd] (μg L−1)Mean pH Value
Aug/XX3818.0 ± 10.97.97 ± 1.15
Sep/XX3726.7 ± 10.98.03 ± 1.08
Oct/XX12113.3 ± 6.07.93 ± 1.06
Nov/XX1295.58 ± 1.97.84 ± 1.31
Dec/XX1723.31 ± 2.077.92 ± 0.91
Jan/YY2222.14 ± 0.788.54 ± 1.16
Feb/YY1793.09 ± 0.658.43 ± 0.78
Mar/YY1514.32 ± 0.818.43 ± 1.30
Apr/YY805.42 ± 1.098.26 ± 1.34
May/YY847.74 ± 0.968.26 ± 1.08
Jun/YY5110.6 ± 1.368.17 ± 1.36
Jul/YY3917.0 ± 4.88.22 ± 1.08
Aug/YY3820.8 ± 6.68.24 ± 1.14
Sep/YY3726.0 ± 8.97.99 ± 1.33
Oct/YY12113.4 ± 5.27.68 ± 1.08
Nov/YY1295.58 ± 1.37.59 ± 1.31
Dec/YY1723.31 ± 1.27.67 ± 0.91
Table 4. Pearson correlation coefficient matrix for monthly rainfall index (MRI), cadmium(II) concentration and pH values for 17 months.
Table 4. Pearson correlation coefficient matrix for monthly rainfall index (MRI), cadmium(II) concentration and pH values for 17 months.
MRIMean CdMean pHCorrelationR2n
MRI1 12
Mean Cd−0.701611 moderate0.492317
Mean pH0.03187−0.096851weak0.009417
Table 5. Mean monthly lead(II) concentration (μg L−1) with standard deviation value determined for each studied well (W) during the 17 months.
Table 5. Mean monthly lead(II) concentration (μg L−1) with standard deviation value determined for each studied well (W) during the 17 months.
MonthW1W2W3W4W5W6W7
Aug/XX35.2 ± 9.627.2 ± 0.719.3 ± 3.422.9 ± 2.626.7 ± 5.624.2 ± 1.823.0 ± 3.2
Sep/XX18.7 ± 11.419.9 ± 0.237.7 ± 13.116.9 ± 3.128.2 ± 7.940.5 ± 5.523.4 ± 2.1
Oct/XX7.30 ± 0.5014.3 ± 3.518.18 ± 0.718.67 ± 0.1519.0 ± 4.413.5 ± 1.610.2 ± 1.9
Nov/XX4.05 ± 0.239.49 ± 0.696.23 ± 1.863.93 ± 0.238.12 ± 1.105.77 ± 0.476.27 ± 0.55
Dec/XX2.07 ± 0.593.50 ± 0.462.17 ± 0.201.46 ± 0.604.32 ± 0.593.39 ± 0.591.39 ± 0.13
Jan/YY2.77 ± 0.143.56 ± 1.102.06 ± 0.362.24 ± 0.253.54 ± 0.601.73 ± 0.753.39 ± 1.02
Feb/YY6.67 ± 1.725.00 ± 0.552.86 ± 0.173.49 ± 0.314.79 ± 0.482.89 ± 0.214.85 ± 0.47
Mar/YY8.50 ± 0.748.09 ± 0.403.49 ± 0.183.91 ± 0.106.08 ± 0.153.82 ± 0.065.36 ± 0.16
Apr/YY12.5 ± 4.810.3 ± 0.34.80 ± 0.134.82 ± 0.816.95 ± 0.114.23 ± 0.085.85 ± 0.13
May/YY16.8 ± 0.411.5 ± 0.77.41 ± 1.788.73 ± 2.3310.7 ± 2.86.52 ± 0.846.72 ± 0.53
Jun/YY17.4 ± 0.120.0 ± 3.411.0 ± 0.614.4 ± 0.2416.3 ± 0.89.05 ± 1.149.56 ± 2.70
Jul/YY23.0 ± 3.237.4 ± 5.3513.4 ± 2.116.4 ± 0.2944.9 ± 4.416. 3± 6.517.8 ± 1.89
Aug/YY33.7 ± 6.941.6 ± 8.518.6 ± 2.922.5 ± 1.6541.1 ± 0.937.8 ± 4.925.4 ± 5.3
Sep/YY23.4 ± 2.116.7 ± 4.435.5 ± 8.818.3 ± 2.532.7 ± 6.134.3 ± 9.029.4 ± 7.2
Oct/YY12.2 ± 1.110.8 ± 4.58.94 ± 1.1011.5 ± 3.813.5 ± 6.716.2 ± 6.412.2 ± 5.1
Nov/YY8.80 ± 2.756.77 ± 3.177.44 ± 1.777.09 ± 3.667.17 ± 1.666.92 ± 1.536.50 ± 2.12
Dec/YY2.62 ± 1.612.78 ± 0.603.04 ± 1.024.46 ± 3.713.24 ± 0.204.54 ± 2.213.36 ± 2.28
Average *13.9 ± 10.314.6 ± 11.510.1 ± 7.111.3 ± 10.816.3 ± 13.613.6 ± 12.911.24 ± 8.8
Minimum2.072.781.462.063.241.731.39
Median12.210.88.77.4410.76.926.72
Maximum35.241.622.937.744.940.529.4
* Average value per well.
Table 6. Mean total lead(II) concentration by well in dry and rainy season.
Table 6. Mean total lead(II) concentration by well in dry and rainy season.
WellDryRainy
W122.6 ± 8.16.10 ± 3.45
W223.1 ± 11.57.14 ± 3.89
W315.6 ± 6.35.19 ± 3.26
W418.5 ± 12.24.93 ± 2.74
W525.9 ± 13.87.75 ± 5.26
W620.8 ± 14.07.10 ± 5.68
W718.7 ± 10.76.41 ± 3.78
Table 7. Mean monthly copper(II) concentration (μg L−1) with standard deviation value determined for each studied well (W) during the 17 months.
Table 7. Mean monthly copper(II) concentration (μg L−1) with standard deviation value determined for each studied well (W) during the 17 months.
MonthW1W2W3W4W5W6W7
Aug/XX14.5 ± 1.224.6 ± 15.58.62 ± 2.049.03 ± 3.059.94 ± 3.326.72 ± 0.5235.2 ± 15.4
Sep/XX27.5 ± 4.232.7 ± 20.08.28 ± 7.0643.6 ± 2.044.80 ± 3.5918.0 ± 1.052.3 ± 42.3
Oct/XX63.5 ± 5.986.4 ± 46.146.7 ± 3.348.1 ± 3.223.3 ± 8.563.8 ± 43.287.9 ± 27.6
Nov/XX116.3 ± 27.0181.6 ± 21.594.9 ± 9.5115.1 ± 26.893.9 ± 46.6157.6 ± 1.6240.4 ± 20.3
Dec/XX356.9 ± 50.9333.0 ± 106.1516.4 ± 168.6637.0 ± 62.2344.0 ± 136.2249.9 ± 122.2569.6 ± 90.8
Jan/YY1031 ± 762479 ± 1671060 ± 73942.5 ± 79.7925.9 ± 131.41949 ± 761971 ± 14
Feb/YY534.2 ± 37.0867.6 ± 191.8601.1 ± 39.8506.6 ± 63.2553.8 ± 26.31207 ± 411791 ± 153
Mar/YY327.1 ± 150.0365.0 ± 43.8350.1 ± 36.6396.4 ± 16.8404.6 ± 0.4474.3 ± 5.12521.0 ± 17.1
Apr/YY152.0 ± 24.1265.1 ± 58.0247.5 ± 56.5338.9 ± 14.4337.7 ± 39.2336.5 ± 65.3433.7 ± 33.3
May/YY71.8 ± 17.0199.1 ± 14.8129.9 ± 62.7186.4 ± 62.2184.0 ± 76.1144.1 ± 9.4282.0 ± 98.2
Jun/YY34.3 ± 13.8106.6 ± 14.241.3 ± 3.096.6 ± 0.679.5 ± 16.827.6 ± 17.684.8 ± 2.98
Jul/YY20.4 ± 3.066.2 ± 27.827.0 ± 1.849.9 ± 9.422.4 ± 0.59.44 ± 1.5357.1 ± 7.3
Aug/YY14.5 ± 1.124.6 ± 12.68.62 ± 1.79.03 ± 2.510.0 ± 2.86.72 ± 0.4235.2 ± 12.6
Sep/YY52.3 ± 4.230.1 ± 12.225.9 ± 20.846.9 ± 4.96.53 ± 5.011.4 ± 7.911.8 ± 7.7
Oct/YY87.9 ± 2.874.9 ± 29.947.4 ± 19.495.4 ± 57.935.0 ± 13.643.6 ± 34.263.8 ± 43.2
Nov/YY240.4 ± 20.3148.0 ± 42.6105.0 ± 20.1171.2 ± 71.294.4 ± 27.5125.7 ± 45.6157.6 ± 1.6
Dec/YY569.6 ± 90.8295.0 ± 80.9426.7 ± 146.6610.1 ± 150.6405.2 ± 119.8222.0 ± 86.3399.9 ± 89.9
Average *218.5 ± 275.3328.2 ± 591.1220.3 ± 289.1253.1 ± 275.1207.9 ± 256.1297.2 ± 516.7399.7 ± 586.9
Minimum14.524.68.289.034.86.7211.8
Median87.914894.9115.193.9125.7157.6
Maximum103124791060942.5925.919491971
* Average value per well.
Table 8. Pearson correlation coefficient matrix for monthly rainfall index (MRI), copper(II) concentration and pH values for 17 months.
Table 8. Pearson correlation coefficient matrix for monthly rainfall index (MRI), copper(II) concentration and pH values for 17 months.
MRIMean CuMean pHCorrelationR2n
MRI1 12
Mean Cu0.802741 strong0.644417
Mean pH0.031870.588401weak0.346217
Table 9. Retention of data variability for the CP1 and CP2 axes and the relationships between cadmium, lead and copper concentrations.
Table 9. Retention of data variability for the CP1 and CP2 axes and the relationships between cadmium, lead and copper concentrations.
VariableCumulative: 98.3%
Coefficient of CP1 (82.2%)Coefficient of CP2 (16.1%)
[Cd2+]0.609250.34847
[Pb2+]0.606570.37380
[Cu2+]−0.510770.85956
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Maroubo, L.A.; Moreira-Silva, M.R.; Teixeira, J.J.; Teixeira, M.F.S. Influence of Rainfall Seasonality in Groundwater Chemistry at Western Region of São Paulo State—Brazil. Water 2021, 13, 1450. https://doi.org/10.3390/w13111450

AMA Style

Maroubo LA, Moreira-Silva MR, Teixeira JJ, Teixeira MFS. Influence of Rainfall Seasonality in Groundwater Chemistry at Western Region of São Paulo State—Brazil. Water. 2021; 13(11):1450. https://doi.org/10.3390/w13111450

Chicago/Turabian Style

Maroubo, Lais A., Marcos R. Moreira-Silva, José Jerônimo Teixeira, and Marcos F. S. Teixeira. 2021. "Influence of Rainfall Seasonality in Groundwater Chemistry at Western Region of São Paulo State—Brazil" Water 13, no. 11: 1450. https://doi.org/10.3390/w13111450

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

Article Metrics

Back to TopTop