About this Journal Submit a Manuscript Table of Contents
Journal of Geological Research
Volume 2013 (2013), Article ID 309638, 9 pages
http://dx.doi.org/10.1155/2013/309638
Research Article

Fluoride in the Serra Geral Aquifer System: Source Evaluation Using Stable Isotopes and Principal Component Analysis

1Universidade Federal de Santa Catarina (UFSC), Campus Universitário, Trindade, 88.010-970 Florianópolis, SC, Brazil
2Universidade Federal do Rio Grande do Sul (UFRGS), Avenida Bento Gonçalves 9500, prédio 43126/103, 91501-970 Porto Alegre, RS, Brazil
3Centro de Pesquisas Geocronológicas (CPGeo), Instituto de Geociências (USP), Universidade de São Paulo (USP), Rua do Lago 562, Cidade Universitária, 05508-080 São Paulo, SP, Brazil
4Universidade do Estado de Santa Catarina, UDESC-FAED, 88.035-001 Florianópolis, SC, Brazil

Received 30 January 2013; Revised 25 April 2013; Accepted 29 April 2013

Academic Editor: Umberta Tinivella

Copyright © 2013 Arthur Schmidt Nanni et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Groundwater with anomalous fluoride content and water mixture patterns were studied in the fractured Serra Geral Aquifer System, a basaltic to rhyolitic geological unit, using a principal component analysis interpretation of groundwater chemical data from 309 deep wells distributed in the Rio Grande do Sul State, Southern Brazil. A four-component model that explains 81% of the total variance in the Principal Component Analysis is suggested. Six hydrochemical groups were identified. δ18O and δ2H were analyzed in 28 Serra Geral Aquifer System samples in order to identify stable isotopes patterns and make comparisons with data from the Guarani Aquifer System and meteoric waters. The results demonstrated a complex water mixture between the Serra Geral Aquifer System and the Guarani Aquifer System, with meteoric recharge and ascending water infiltration through an intensive tectonic fracturing.

1. Introduction

In the last several decades, the global water consumption has dramatically increased, especially for agriculture, water supply, and industrial uses. This paper examines the fluoride content in water from the southernmost region of the fractured subcontinental Serra Geral Aquifer System (SGAS), an important aquifer that supplies a large amount of water for an economic developed region in Southern Brazil.

Fluoride content in water is beneficial to human health and in a moderate concentration (0.7–1.2 mg/L) prevents dental cavities. When used in excess can be toxic causing human and animal dental and skeletal fluorosis, which has been detected in China [1, 2], India [3], Kenya [4], and Israel [5], among other countries. The drinking water limit recommended by the World Health Organization for fluoride is 1.5 mg/L [6]. In the study area, the SGAS fluoride average concentrations are around 0.24 mg/L, with a minimum value of 0.02 m/L and the highest at 3.03 mg/L.

The SGAS overlies the Guarani Aquifer System (GAS; [7]), which has been the focus of several recent studies due to its spatial extent and storage potential as a transboundary aquifer [8]. The area covered by the SGAS in Brazil, Uruguay, Argentina, and Paraguay is equivalent to 1,200,000 km2 (Figure 1). In Brazil, these groundwater resources should be efficiently managed to protect its water potential and quality. The climate in the study area ranges from subtropical to temperate, with precipitation average of 1,550 mm/year.

309638.fig.001
Figure 1: Location map showing the Serra Geral Formation (SGF) in South America (modified from [9]).

The SGAS reaches its maximum thickness of about 1,800 m along the central axis of the Parana Basin, located in the Sao Paulo State. In the Rio Grande do Sul State, this aquifer has a wedge shape, ranging from 1,000 m in the eastern area bordering the Santa Catarina State to thickness of 200 m in the western area adjacent to Argentina Border.

The aim of this study is to understand the spatial distribution and geologic control of fluoride content in the SGAS, identifying hydrochemical groups with high fluoride. Physicochemical analysis, geostatistics, and stable isotopes approaches were used along with a Geographic Information System (GIS) in order to identify the spatial distribution of the different water types.

1.1. Geologic and Hydrogeologic Framework

The study area is located in the Rio Grande do Sul State, Southern Brazil, and underlies 164,200 km2, extending from 27°S to 31°S and from 50°W to 57°W. The volcanic sequence has an average thickness of about 550 m and consists of basaltic to rhyolitic units of Mesozoic age; the last ones occur on the Serra Geral Formation (SGF) top sequence. The clayish soil thickness in the SGF ranges from few decimeters in the eastern area to several meters in the west and northwest. This soil controls water infiltration and, as a consequence, the SGAS recharge by meteoric waters [10].

Studies have suggested groundwater mixing processes between the SGAS and older confined sedimentary aquifer systems from the Parana Basin, including the GAS, through ascending groundwater infiltration [1116].

Distinct tectonic fracturing occurs in the study area, controlling geomorphological features and groundwater circulation. These structures are directly related to the opening of the South Atlantic Ocean, which is reflected in the major NE and NW fault and fracture directions. Neotectonics is considered to play an important hydrogeological role in the eastern region. Major tectonic systems define individual hydrochemical sectors, which are considered as individual hydrogeological blocks [15, 17, 18]. The tectonic block limited by the Terra de Areia-Posadas and Mata-Jaguari Fault Systems (Figure 2) is uplifted with respect to the neighboring blocks. In the north, the adjacent block is divided into smaller units with a gradual terrain lowering from east to west, influenced by NE normal faults, parallel to the Leao and Perimpo Fault Systems.

309638.fig.002
Figure 2: Sample points, hydrochemical groups distribution on SGAS, and major tectonic fault systems: , Mata-Jaguari, (2) Terra de Areia-Posadas, (3) Perimpo, and (4) Leao.

The identification of hydrochemical groups in the SGAS is an important tool for the definition of tectonic structures, supporting the fluoride ascending circulation hypothesis in the aquifer [22]. The meteoric recharge pattern as waters and also ascending infiltrations by waters with and in the GAS, resulting in a complex water mixture [15]. Several authors have pointed out that high fluoride content in groundwater can be related to ascending infiltrations waters [11, 15, 2224]. Such mixing and hydrochemical groups can both result from and waters [25].

2. Materials and Methods

2.1. Multivariate Statistics

Principal components analysis (PCA) was used to discriminate hydrochemical groups in the 309 physicochemical analyses of the SGAS in order to systematize the interpretation of a large number of variables through the maximization of differences between them [26].

The PCA was performed to 2, 3, and 4 principal components and normalized to a value of “1,” minimum eigenvalue accepted [27]. Subsequently, the outliers were cut by visual analysis in boxplots and dispersion charts for , , , , , , and , considering water classification. This statistical procedure restricts the analysis towards the fluoride focus, avoiding interferences from other parameters. The spatial distribution and interpretation of principal components defined by PCA were performed with a GIS environment.

2.2. Stable Isotopes

Stable isotopes were analyzed in the Geochronological Research Center of the Sao Paulo University. The analysis for oxygen was accomplished through the balanced isotopic reaction between water and CO2 over 24 hours.

The CO2 extraction was done in a vacuous line off-line and subsequently was read in a Thermo Finnigan Delta Plus/Advantage spectrometer with dual-inlet installed. The hydrogen was analyzed through the reaction of water with pure zinc at 500°C for liberation of the H2 gas.

Twenty seven samples were selected in order to establish stable isotope patterns in the SGAS. These results were considered together with eight additional ones belonging to the GAS (seven data from [28] and one from this study) and four from phreatic waters from an adjacent area [28]. This procedure is appropriate for understanding water mixture features between aquifers and to estimate their relative ages.

Unfiltered water samples were collected and stored in 100 mL polypropylene bottles to provide adequate volume for separate analyses of δ18O and δ2H. All data considered in this paper are expressed in ‰  with respect to Vienna Standard Mean Ocean Water (VSMOW) on the Delta scale: where is the 18O/16O or 2H/1H ratio of the sample and is the corresponding ratio in VSMOW. Analytical precisions were better than ±0.5‰  for δ18O and ±2‰  for δ2H.

The altitude and the distance from the ocean can influence stable isotopes patterns. The average altitude of the SGAS in the area is 500 m, and the spatial position ranges from the coast in the eastern sector to 800 km inland in the western sector.

Due to this, ISOHIS [29] isotope data from the precipitation sampling stations were selected both in Porto Alegre, Southern Brazil, located in coastal and low altitude condition (GNIP code 8396700) and Corrientes, Argentina, located in interior and low altitude (GNIP code 8716600) area showed to compare data obtained in this research.

3. Results and Discussion

3.1. Principal Components Analysis

The PCA results show the correspondence of each component to the analyzed parameters. The solution using the eigenvalue criterion results in four components that explain 81% of the total variance. Component 4 is highly weighted by fluoride (Table 1).

tab1
Table 1: Rotated matrix with four components and characteristic parameters. Data from [21].

Component 1 contains , , and , typical for calcium and magnesium bicarbonated waters with short to medium residence time. In Component 2, and prevail, indicating calcium and sodium-chlorinated and sulfate waters. Component 3 is distinguished by and , representing sodium bicarbonated waters. The followed by and defines Component 4, which comprises the high group linked to sodium-sulfated waters.

Based on these components, the solution for nine clusters was selected, and it shows the best approach to the hydrochemical groups. Subsequently, a cluster interpretation was performed based on central scores of each cluster. In order to assist the interpretation, each cluster was renamed according to the predominant ion (Table 2), defining hydrochemical groups. Thus, Component 1 has been renamed to CaMg, Component 2 to SO4, Component 3 to Na, and Component 4 to F.

tab2
Table 2: Cluster central scores, dominant components and hydrochemical groups. Data from [21].

Clusters 3 and 4 are composed predominantly by the F group. Clusters 1 and 7 comprise the CaMg group. Cluster 8 is composed by the Na group. Other clusters show prevalence of two components, resulting in the groups SO4Na, SO4F, and NaF. Cluster 5 comprises the Meteoric waters and other samples with no relationship to the four principal components.

3.2. Hydrochemistry

The geochemical data based on PCA results were plotted in a piper diagram (Figure 3), with the four renamed clusters defined by the PCA. The majority of wells (212 out of 309 wells) have a composition close to meteoric waters with nature and do not include high fluoride waters. Due to this, 212 samples are not represented in the diagrams in order to obtain a better fit by the other groups.

309638.fig.003
Figure 3: Hydrochemical groups distribution in the piper diagram, showing two fluoride hydrochemical types (I—sulfated and II—bicarbonated, [19]).

The CaMg group represents a predominant water type where appears in more than 50% of the wells (Figure 3).

The Na group shows their typical distribution in the piper diagram, representing waters.

SO4Na, SO4F, and NaF groups appear only in two wells each. Both wells compose distinct groups as shown in Figure 3.

F distribution in the piper diagram maintains an association with sulfated and bicarbonated groundwaters. These hydrochemical characteristics are also related to different groundwater sources or point out to mixture features during groundwater ascending infiltration, where ionic enrichment including increasing fluoride content can result from the interaction of ascending/descending water recharge due to longer residence time in the Na and CaMg hydrochemical groups [19]. Thus, the SO4F, NaF, F, and SO4Na hydrochemical groups have a closed relationship interaction with two or more fracture directions (Figure 4), where , , and fluoride can ascend.

309638.fig.004
Figure 4: Water flow interactions to different geological conditions.
3.3. Stable Isotopes

The 27 analyzed samples from SGAS are distributed among the six water groups defined in the PCA. For GAS waters, one sample was analyzed and aggregated other 7 results [20]. Another 4 results to phreatic waters [28] were aggregated as mentioned. The δ18O (‰) ranges from −5.0 to −4.2 (phreatic waters), from −5.7 to −4.2 (SGAS), and from −5.5 to −4.8 (GAS), while the δ2H (‰) varies from −27.4 to −20.8 (phreatic waters), from −37 to −24 (SGAS), and from −37 to −28 (GAS; Table 3). The SGAS shows a scattered distribution compared to the GAS (Figure 5). Nevertheless, the average of δ18O and δ2H values in both aquifer systems is similar, and the SGAS is located in an intermediate position but nearer to GAS when compared to phreatic waters. This reinforces the hypothesis of water circulation through the SGAS to the GAS, with minor influence of phreatic/meteoric waters in the sampled field.

tab3
Table 3: δ 18O and δ2H values of phreatic waters (a), GAS (a), and SGAS (b). Data from [28] (*) and [20] (**). TA is a code to Tres Arroios City (sampled by author). Fluoride content is exclusive to SGAS.
309638.fig.005
Figure 5: δ18O (‰) versus δ2H (‰) in GAS, SGAS, and phreatic aquifer.

The GAS data plotted in a map (Figure 6) show decreasing δ2H and δ18O values from east to west, reflecting the distance from the ocean and the groundwater flow direction through the GAS. The lowest δ18O and δ2H values in the GAS are related to the highest confinement and distance from direct recharge areas. In this sense, tectonic fractures can provide meteoric to medium residence time waters to the GAS across the Serra Geral Formation, indicating descending water recharge processes, corresponding to Condition 2 described into Figure 4.

fig6
Figure 6: Hydrogeological spatial trends of δ2H (a) and δ18O (b) in the GAS (data from [20]) and major tectonic fault systems: Mata-Jaguari, (2) Terra de Areia-Posadas, (3) Perimpo, and (4) Leao.

The SGAS δ18O and δ2H results plotted in the same map show a different arrangement from the GAS, with the lowest values in the northern, western, and eastern portions (Figure 7). These differences are related to the influence of tectonic fractures and, as a consequence, to the hydrochemical groups distribution [19].

fig7
Figure 7: Hydrogeological spatial trends of δ2H (a) and δ18O (b), hydrochemical groups in the SGAS and major tectonic fault systems: , Mata-Jaguari, (2) Terra de Areia-Posadas, (3) Perimpo, and (4) Leao.

In the northern region, the lowest δ18O and δ2H values coincide with a greater volcanic sequence thickness and thicker soils, providing for an increase in groundwater residence time and reduced meteoric water recharge. The prevalence of SO4F, NaF, and F hydrochemical groups in this portion confirms the long residence time and older relative age to this waters.

In the western sector of Figure 6, samples with the lowest δ18O and δ2H values are located along the Uruguay River Fault System in an NE-SW direction. In this area, the SGAS is thinner, and the intense tectonic fracturing favors to the groundwater ascension through the confined GAS.

In the eastern area, the lowest δ18O and δ2H values (samples 238 and 244) are related to the Leao and Terra de Areia-Posadas Fault Systems, which allow for groundwater ascension through the confined GAS.

The scatter diagram of fluoride versus δ18O and δ2H (Figure 8) reveals, in a broad sense, that lower isotope ratios (most negative values) are associated with higher fluoride content, reflecting water-rock interaction during a long residence time. Similar conditions were also described by [24] in the GAS in the Sao Paulo State, Southern Brazil. This behavior can be associated with ascending infiltration GAS waters.

fig8
Figure 8: Fluoride versus δ18O and δ2H and the SGAS hydrochemical groups.

4. Conclusions

In the Rio Grande do Sul State, extensive groundwater mixing processes operate between the SGAS and GAS through defined and intense tectonic fractures. Six hydrochemical groups were distinguished in the SGAS, three of them containing substantial fluoride contents in combination with and . Ascending recharge processes are considered to define ionic enrichment, especially in , , , and .

The comparison between δ18O and δ2H spatial distribution in the SGAS demonstrates a different scenario concerning the GAS. The differences are probably explained by the interference of different tectonic structures than those in the GAS which are responsible for the defined hydrogeological blocks. The combination of deep fault zones and high pressure and temperature conditions in a confined GAS supports the hypothesis of ascending groundwater infiltration.

The depleted δ18O and δ2H values suggest that fluoride enrichment in the SGAS is related to percolation through the GAS in confinement and long residence time conditions. The similar values of δ18O and δ2H in the two aquifer systems suggest, in some regions, that ascending recharge processes are more intense than meteoric recharge. In the Rio Grande do Sul State, this behavior increases from east to west, which is the groundwater flow direction in the GAS. Similar behavior is observed in the SGAS, except in some scarps and valleys, as is demonstrated by the distribution of the CaMg water group. The fluoride input to the groundwater is probably related to water-rock interaction in the GAS and in older aquifer systems below it.

Acknowledgments

The funding for this study was provided by the Brazilian National Research Council (CNPq), Rio Grande do Sul Research Foundation (Fapergs), and Rio Grande do Sul Environmental Agency (Fepam).

References

  1. N. Lin, J. Tang, and J. Bian, “Geochemical environment and health problems in China,” Environmental Geochemistry and Health, vol. 26, no. 1, pp. 81–88, 2004. View at Publisher · View at Google Scholar · View at Scopus
  2. W. Genxu and C. Guodong, “Fluoride distribution in water and the governing factors of environment in arid north-west China,” Journal of Arid Environments, vol. 49, no. 3, pp. 601–614, 2001. View at Publisher · View at Google Scholar · View at Scopus
  3. V. Vijaya Kumar, C. S. T. Sai, P. L. K. M. Rao, and C. S. Rao, “Studies on the distribution of fluoride in drinking water sources in Medchal Block, Ranga Reddy District, Andhra Pradesh, India,” Journal of Fluorine Chemistry, vol. 55, no. 3, pp. 229–236, 1991. View at Scopus
  4. W. K. Moturi, M. P. Tole, and T. C. Davies, “The contribution of drinking water towards dental fluorosis: a case study of Njoro division, Nakuru district, Kenya,” Environmental Geochemistry and Health, vol. 24, no. 2, pp. 123–130, 2002. View at Publisher · View at Google Scholar · View at Scopus
  5. U. Kafri, A. Arad, and L. Halicz, “Fluorine occurrence in groundwater in Israel and its significance,” Journal of Hydrology, vol. 106, no. 1-2, pp. 109–129, 1989. View at Scopus
  6. WHO, Fluorides, Environmental Health Criteria no. 227, United Nations Environmental Programme, International Labour Organization, World Health Organization, Geneva, Switzerland, 2002.
  7. H. C. N. S. Campos, “Modelación conceptual y matemática del acuífero guaraní, Cono Sur,” Acta Geológica Leopoldensia, vol. 23, no. 4, pp. 3–50, 2000.
  8. L. M. Aráujo, A. B. França, and P. E. Potter, “Hydrogeology of the Mercosul aquifer system in the Paraná and Chaco-Paraná Basins, South America, and comparison with the Navajo-Nugget aquifer system, USA,” Hydrogeology Journal, vol. 7, no. 3, pp. 317–336, 1999. View at Scopus
  9. A. Batezelli, A. R. Saad, V. J. Fulfaro, A. C. Corsi, P. M. B. Landim, and J. A. J. Perinotto, “Análise de bacia aplicada às unidades mesozoicas do triângulo mineiro (sudeste do Brasil): uma estratégia na prospecção de recursos hídricos subterrâneos,” Águas Subterrâneas, vol. 19, no. 1, pp. 61–73, 2005.
  10. P. A. R. Reginato and A. J. Strieder, “Integração de dados geológicos na prospecção de aqüíferos fraturados na Formação Serra Geral,” Revista Águas Subterrâneas, vol. 20, no. 1, pp. 1–14, 2006.
  11. M. Szikszay, J. Teissedre, U. Barner, and E. Matsui, “Geochemical and isotopic characteristics of spring and groundwater in the State of São Paulo, Brazil,” Journal of Hydrology, vol. 54, no. 1–3, pp. 23–32, 1981. View at Scopus
  12. G. Fraga, Origem de fluoreto em águas subterrâneas dos Sistemas Aqüíferos Botucatu e Serra Geral da Bacia do Paraná [Ph.D. thesis], Programa de Pós-Graduação em Recursos Minerais e Hidrogeologia, Universidade de São Paulo, 1992.
  13. C. V. Portela Filho, F. J. F. Ferreira, E. F. Rosa Filho, A. C. Buchmann, and S. P. Rostirolla, “Estudo preliminar da conexão entre os aquíferos Serra Geral e Guarani com base em dados aeromagnetométricos e hidroquímicos,” in XII Congresso Brasileiro de Águas Subterrâneas, vol. 1, pp. 74–86, Editora da ABAS, São Paulo, Brazil, 2002.
  14. C. V. Portela Filho, F. J. F. Ferreira, E. F. Rosa Filho, A. C. Buchmann, and S. Rostirolla, “Compartimentação magnética-estrutural do sistema aquífero Serra Geral e sua conectividade com o Sistema Aqüífero Guarani na região central do arco de Ponta Grossa (Bacia do Paraná),” Revista Brasileira de Geociências, vol. 3, pp. 369–381, 2005.
  15. J. L. F. Machado, Compartimentação espacial e arcabouço hidroestratigráfico do Sistema Aqüífero Guarani no Rio Grande do Sul [Ph.D. thesis], Programa de Pós-Graduação em Geologia—Área de Concentração em Geologia Sedimentar, Universidade do Vale do Rio Dos Sinos, São Leopoldo, Brazil, 2005.
  16. L. F. Scheibe and R. Hirata, “O contexto tectônico dos sistemas aqüíferos guarani e serra geral em santa catarina: uma revisão,” in XV Congresso Brasileiro de Águas Subterrâneas, pp. 1–14, Associação Brasileira de Águas Subterrâneas, Curitiba, Brazil, 2008.
  17. N. A. Lisboa, Fácies, estratificações hidrogeoquímicas e seus controladores geológicos, em unidades hidrogeológicas do Sistema Aqüífero Serra Geral, na Bacia do Paraná, Rio Grande do Sul [Ph.D. thesis], Curso de Pós-Graduação em Geociências, Universidade Federal do Rio Grande do Sul, 1996.
  18. N. A. Lisboa and E. Menegotto, “Diferenciações hidrogeoquímicas no sistema aqüífero Serra Geral no Rio Grande do Sul,” in XII Simpósio Brasileiro de Recursos Hídricos, pp. 489–496, 1997.
  19. A. Nanni, A. Roisenberg, M. P. C. Marimon, and A. P. Viero, “The hydrochemical facies and anomalous fluoride content in the Serra Geral aquifer system, southern Brazil: a GIS approach with tectonic and principal component analysis,” Environmental Geology, vol. 58, no. 6, pp. 1247–1255, 2009. View at Publisher · View at Google Scholar · View at Scopus
  20. A. A. Kimmelmann e Silva, A. C. Reboucas, and M. M. F. Santiago, “14C analyses of groundwater from the Botucatu Aquifer System in Brazil,” Radiocarbon, vol. 31, no. 3, pp. 926–933, 1989. View at Scopus
  21. A. Nanni, A. Roisenberg, J. M. G. Fachel, G. Mesquita, and C. Danieli, “Fluoride characterization by principal component analysis in the hydrochemical facies of Serra Geral Aquifer System in Southern Brazil,” The Anais da Academia Brasileira de Ciências, vol. 80, no. 4, pp. 693–701, 2008.
  22. O. A. B. Licht, A geoquímica multielementar na gestão ambiental—identificação e caracterização de províncias geoquímicas naturais, alterações antrópicas da paisagem, áreas favoráveis à prospecção mineral e regiões de risco para a saúde no estado do Paraná, Brasil [Ph.D. thesis], Pós-Graduação em Geologia, Área de Concentração—Geologia Ambiental, Setor de Ciências da Terra, Universidade Federal do Paraná, Curitiba, Brazil, 2001.
  23. A. V. L. Bittencourt, E. F. Rosa Filho, E. C. Hindi, and A. C. Buchmann Filho, “A influência dos basaltos e de misturas com águas de aquíferos sotopostos nas águas subterrâneas do Sistema Aqüífero Serra- Geral na bacia do rio piquiri, Paraná—BR,” Revista Águas Subterrâneas, vol. 17, pp. 67–75, 2003.
  24. O. Sracek and R. Hirata, “Geochemical and stable isotopic evolution of the Guarani Aquifer System in the state of São Paulo, Brazil,” Hydrogeology Journal, vol. 10, no. 6, pp. 643–655, 2002. View at Publisher · View at Google Scholar · View at Scopus
  25. E. F. Rosa Filho, J. X. Montano, and U. Duarte, “Estudo do movimento das águas subterrâneas do Sistema Aqüífero Guarani (SAG) através de isótopos, no Paraná, São Paulo e no Uruguai,” Revista Latino Americana de Hidrogeologia, pp. 109–121, 2007.
  26. A. O. Invernizzi and S. M. B. Oliveira, “Hydrochemical characterization of a watershed through factor analysis,” Revista de Águas Subterrâneas, vol. 18, pp. 67–77, 2004.
  27. H. F. Kaiser, “The varimax criterion for analytic rotation in factor analysis,” Psychometrika, vol. 23, no. 3, pp. 187–200, 1958. View at Publisher · View at Google Scholar · View at Scopus
  28. M. P. Marimom, Origem das Anomalias de Fluoreto nas Águas Subterrâneas utilizadas para abastecimento público da Formação Santa Maria, Estado do Rio Grande do Sul, Brasil [Ph.D. thesis], Curso de Pós-graduação em Geociências, Universidade Federal do Rio Grande do Sul, Porto Alegre, Brazil, 2006.
  29. IAEA—International Atomic Energy Agency, Isotope Hydrology Information System, the ISOHIS Database, 2004, http://isohis.iaea.org.