The Uraniferous Province of Lagoa Real (Bahia, Brazil) is an area naturally enriched in U and Th, with the occurrence of Pb derived from the radioactive decay of U (238U e 235U) and Th (232Th), in the form of their stable isotopes 206Pb, 207Pb and 208Pb respectively, in addition to the presence of the natural isotope 204Pb. These mineralizations in uranium are the basis of the mining activity of this element, largely responsible for the economic development of the region. In these areas, the groundwater represents the main source of water for the population, since water scarcity is significant.
According to Licht (2001) , most of the natural or anthropic processes occurring on the Earth’s surface are characterized by associations of chemical elements (characteristic of the generating process) and geographically associated to their area of occurrence. However, these geochemical signatures may be masked by the superposition of other ones from different sources. Thus, the identifiable geochemical signal, e.g. in waters and sediments, would represent a sum of contributions originated from the natural environmental processes and human action. According to this author, when a complex geochemical signal is interpreted using a multidisciplinary approach, it is possible to decompose it, and, hence, natural or anthropogenic processes can be identified and isolated.
In order to investigate the correlation between a given element and the source that originated it (natural or anthropogenic), works from the 1960s onwards showed the applicability of the methodology of Pb isotope ratios in environmental studies (Nriagu, 1989) . Considering that physical or chemical processes do not affect the isotopic composition of Pb in the environment (Doe, 1970; Bollhöfer et al., 1999) , the Pb signature tends to be identical to the source that originated it. Therefore, the isotopic signatures of Pb of anthropogenic provenance are distinct from the isotopic signature of natural origin, and can be used as a tracer of the possible sources of Pb in the environment (Bollhöfer et al., 1999; Bollhöfer and Rosman, 2000; Bollhöfer and Rosman, 2001; Chen et al., 2005; Gioia et al., 2006; Gulson et al., 2007; Komárek et al., 2008; Gioia et al., 2010; Vecchia, 2015; Vecchia et al., 2015; Vecchia et al., 2017a and 2017b) .
In the area of the Gneissic-Granitic Complex, Cordani et al. (1992) obtained data on isotopic composition of Pb in granite rocks of the São Timóteo type (non-deformed granites) and in orthogneisses (deformed granites and associated with uraniferous mineralization). Subsequently, Iyer et al. (1999) demonstrated the applicability of Pb isotopes in the quantitative determination of U and Th mobility in different geological formations present in the São Francisco Craton. These authors studied several geological environments belonging to the São Francisco Craton, among which is the Lagoa Real Gneissic-Granitic Complex, using the isotopic data of Pb published by Cordani et al. (1992) .
According to Vecchia (2015) , Vecchia et al. (2015) and Vecchia et al. (2017a and 2017b) , the use of the Pb isotopic tool to study environmentally impacted areas requires knowledge of the Pb isotopic signatures of geographic samples as a determinant factor in the investigative process of environmental matrices. Therefore, this study uses the Pb isotopic data of geographic samples of the study area obtained by Cordani et al. (1992) along with the U, Th and Pb concentration data from Iyer et al., 1999 . The isotopic data of Pb produced by Cordani et al. (1992) , containing information on the preserved geological environment, were compared to the isotopic data of Pb from Vecchia et al. (2017b) , in order to diagnose influences of geogenic and/or anthropogenic sources in this area.
Statistical analyses allow the specification of the variables that best characterize the geogenic influences coming from the geological context and/or the anthropic activities in the environment. The multivariate statistical methods allow the extraction of useful information capable of reducing the dimensional representation of the data, organizing them in a structure that facilitates the visualization of the whole data set. Principal Component Analysis (PCA) is among the most widely used with this purpose. It compacts the data into a smaller number of variables (principal components) whose variance (and importance) is given by their eigenvalues. These components are mathematical functions that generate scores that have embedded the contributions (and importance) of each variable with its weight (loadings) (Manly, 2005; Mingoti, 2007; Zhao et al., 2011, Cavalcante et al., 2013; Sun et al., 2013; Siepak and Sojka, 2017) . It is important to emphasize that the results obtained with these techniques should, whenever possible, be confronted with information related to the research context, seeking an association on the possible processes in the study region (Vecchia, 2015; Vecchia et al., 2015; Vecchia et al., 2017a and 2017b) .
This work presents a statistical study based on the PCA, using Pb isotopic ratios, and U, Th and Pb concentrations from groundwater and sediments in the Lagoa Real Uranium Complex (Brazil). The main objective is to refine the information from the investigations carried out by Vecchia et al. (2017b) by analyzing the contributions of each variable in the diagnosis of possibly impacted areas in the Uraniferous Province area. This work aims to diagnose not only the water quality and sediments in the study area, but also to identify the areas possibly affected by the mining activity.
2. Study Area
The study area, Southeastern Brazil, is located in the semi-arid region characterized by a marked water deficit. It is considered the main active site of uranium in Brazil and South America. Pb isotopic signature of samples was determined in an area of 1200 km2 at Lagoa Real Uranium Complex, Bahia State, Brazil (Figure 1). This complex is geologically identified as a Gneissic-Granitic Complex with ages between of 1.72 and 1.77 Ga (Cordani et al. 1992) . Granite type São Timóteo (undeformed granites), orthogneiss (deformed granites and associated U-mineralization) and U-albitites (main hosts of the uranium mineralization) are the principal lithologies of the region and represents the geogenic sources (Vecchia et al., 2017b) .
The average annual rainfall is around 700 mm, with high variability in the spatial and temporal distribution (annual seasonality) (Vecchia, 2015) . According to Lamego Simões et al. (2003) , the climatological and hydrological characteristics,
Figure 1. Simplified geologic map with the mainlithologies representative of the geogenic sources in the Lagoa Real Uranium Complex: Albitites, Orthogneiss and São Timóteo Granites and the sample locations of sediments and groundwater. Modified by Vecchia et al. (2017b) .
associated to the conformation of the regional relief with flows to Atlantic slope, give rise to a hydrographic network in which intermittent drainage is common (Figure 1). Complementary hydrological data and geology of the study area can be found in Vecchia (2015) .
3. Experimental Methods
Sampling was carried out according to Vecchia et al. (2017b) . Groundwater and sediment sampling points were selected according to regional hydrology and geology, including points at the mining area and points located upstream and downstream thef U anomalies (contents higher than 1500 ppm). The deficit of surface water led to the choice of groundwater and sediment as study matrices. Sampling was carried out in the rainy season (December to January) because of the discontinuous flow of the rivers and the low precipitation level in the region (around 700 mm). The location of sampling was devised aiming to achieve a reliable geographical representation and to ensure a good spatial distribution of the data. Thirty samples were analyzed; 12 of groundwater (AP) and 18 of sediments (S) as shown in Figure 1. The sampling points were georeferenced, using GPS (Global Positioning System).
The sampling and preservation of water samples were carried out according to the Environmental Agency of the State of São Paulo directions— CETESB (2011) and Veridiana et al. (2008) . More details in Vecchia et al. (2017b) .
Approximately 1 kg of sediments samples, representative of the regional lithologies (Figure 1), were collected in depth from 50 cm to 80 cm, using a cylindrical collector. This material was transferred to low-density polyethylene (LDPE; NalgonÒ) containers. Subsequently, these samples were dried, disaggregated, homogenized and sieved. Fractions < 250 Mesh Tyler were sent to isotopic Pb and chemical analysis.
3.2. Pb Isotopic Geochemical Analytical Procedure
Pb isotopic ratios of groundwater and sediments were determined through Thermal Ionization Mass Spectrometry (TIMS), using a multi-collector Finnigan MAT 262 at the Geochronological Research Center of USP (University of São Paulo) according to the Vecchia et al. (2017b) .
3.3. Pb, U and Th Geochemistry Analytical Procedure
U, Th and Pb were analyzed in water and sediments samples by Inductively Couple Plasma Mass Spectrometry, ICP-MS (PerkinElmer - ELAN DRC-e) at the Center for Development of Nuclear Technology (CDTN), according to the Vecchia et al. (2017b) and procedures recommended by the U.S. EPA 2008.
3.4. Multivariate Analysis
Data were submitted to a Principal Component Analysis aiming to establish relations between U, Th and Pb contents and Pb isotopic ratios (206Pb/204Pb, 207Pb/204Pb, 208Pb/204Pb, 206Pb/207Pb, 208Pb/207Pbe 208Pb/206Pb) to the environment and their origin (geogenic or anthropogenic). The whole dataset (Tables 1-3 in Annex) consists of 12 sediment samples, 12 water samples (11 groundwater samples and 1 surface water sample) and 15 samples of Cordani et al. (1992) (07 São Timóteo Granite samples and 08 orthogneiss samples). MinitabÒ software, version 22.214.171.124, was used for statistical analysis.
4. Results and Discussion
Data (Tables 1-3 in Annex) were submitted to a Principal Component Analysis. This study was conducted in three steps. In the first one U, Th and Pb concentrations and Pb isotopic ratios of orthogneisses samples were analyzed together with data of groundwater and sediment samples. In the second step, only data of isotopic ratios from the same samples were used. In the third step, data of Pb isotopic ratios of São Timóteo granite were used together with data of groundwater and sediment samples. It is important to point out that these two types of rocks are characteristic of the geological formation of the study area and involve the genesis of uranium mineralization. Thus, they are important references of the geographic sources of Pb in the region. The use of these data aims to determine the type of rock that most influences the Pb isotopic signatures of groundwater and sediments and, therefore, which one better discriminates the environmental matrices investigated.
1) PCA using Pb isotopic ratios and U, Th and Pbconcentrations in sediments, groundwater and orthogneiss samples
Table 1 presents the loading values of each variable in each of the components as well as their respective eigenvalues. These data can be interpreted in terms of two main components, PC1 and PC2, which account for 81.2% of the data variability. The first component (PC1) with an eigenvalue of 4.6145 accounts for 51.3% of the data variance. It is mainly dominated by isotopic ratios206Pb/204Pb, 207Pb/204Pb, 206Pb/207Pb, 208Pb/206Pb and by U. This component can be interpreted as representing the contribution of this element and its decay to the geochemistry of sediments and water. The second component (PC2) with an eigenvalue of 2.6974 accounting for 30.0% of the data variance is dominated by the isotopic ratios 208Pb/204Pb, 208Pb/207Pb and by Th and Pb. This component can be understood as representing the contribution of Th and its decay in the general characterization of sediments and waters.
Table 1. Loading values of each variable in each of the principal components for samples of water, sediment and orthogneisses of the uraniferous province of Lagoa Real BA, for PCA carried out with all variables.
Figure 2 shows the score plot of the first two components (PC1 and PC2) for water, sediment and orthogneiss samples. The score of a sample represents the contribution of each variable measured and its respective weights to the characteristics of the sample. Sediment and rock samples are grouped around the origin of the diagram, showing little difference from each other, which shows the geogenic character of the sediments, since in the process of weathering they did not undergo enrichment, i.e., anthropic influence. Sample P19S stands out from this group. Its high score values of PC1 and PC2 suggests enrichment of U and Th whose increase in concentration leads to an increase in the score value for the score of the samples in each component. Also noteworthy is the sample P35S with positive score value of PC1 and negative score value of PC2. In this case, the sample is enriched in U rather than Th, being the U enrichment due to anthropic influence. The water samples stand out from the group with negative score values for both PC1 and PC2. This means that, in general, the values of the variables for this matrix are low, which is a reflex of the low solubility of U, Th and Pb species (and their isotopes) in water. Samples P22AP is enriched in these species while P18A is enriched in Th. Samples P22AP and P18AP stand out from the rest of the groundwater samples, resembling the sediment and orthogneiss group of samples, probably because their isotopic signatures are under the geogenic influence of the orthogneiss.
2) PCA using Pb isotopic ratios in sediments, groundwater and orthogneiss samples
The same procedure was applied to the same samples, removing U, Th and Pb concentrations, i.e., working only with the Pb isotopic ratios as variables. Table 2 shows the loading values for both components. Concerning the isotopic ratios there is no difference between the results obtained with the two approaches, i.e., calculations with and without U, Th and Pb concentrations. Figure 3 and Figure 4 reflect this result. The distributions of samples in both diagrams are the same.
Figure 2. Score plot of PC1 and PC2 for samples of water, sediment and orthogneisses of the uraniferous province of Lagoa Real BA, carried out with all variables.
Table 2. Loading values of each variable in each of the principal components for samples of water, sediment and orthogneisses of the uraniferous province of Lagoa Real (BA), for PCA carried without U, Th and Pb concentrations.
Figure 3. Score plot of PC1 and PC2for samples of water, sediment and orthogneisses of the uraniferous province of Lagoa Real BA, without U, Th and Pb concentrations, using Pb isotopic ratios.
Figure 4. Score plot of PC1 and PC2for samples of water, sediment and granite of the uraniferous province of Lagoa Real BA, using Pb isotopic ratios.
Exception is sample P19S for which the concentrations of these three elements are especially higher. This suggests that the isotopic ratios are sufficient to assess the characteristics of samples. An additional observation due to the absence of the concentration variable in this step is the increase in the percentage of explanation of the variability of the data of the principal components PC1 and PC2, evidencing the importance and efficacy of the isotopic ratios in the study of environmental samples.
3) PCA using Pb isotopic ratios data for sediments, groundwater and São Timóteo Granite samples.
For this matrix, only data of Pb isotopic ratios were used. They refer to the samples of São Timóteo type granite, groundwater and sediments. Table 3 shows the loading values of each variable in each of the main components.
These data can be interpreted in terms of two principal components, PC1 and PC2, which account for 97.8% of the data variability. The first component (PC1) with an eigenvalue of 3.9180 accounts for 65.3% of the data variance. It is mainly dominated by isotopic ratios206Pb/204Pb, 207Pb/204Pb, 206Pb/207Pb, 208Pb/206Pb the latter with negative loading value. This component can be interpreted as representing the contribution of U represented by its decay isotopes to the geochemistry of sediments. The second component (PC2) with an eigenvalue of 1.9517 accounting for 32.5% of the data variance is dominated by the isotopic ratios 208Pb/204Pb, 208Pb/207Pb, which are the decay products of Th. It represents the contribution of Th in the general characterization of sediments.
Figure 4 shows the score plot of the first two components (PC1 and PC2) for water, sediment and granite samples. In general, the same observations made with the previous PCA study (Figure 3) can be summarized here: no separation of sediment and rock samples leading to their geogenic character, separation of water samples, separation of samples P22AP, and P18AP, P19S and one sample (P35S) due to anthropogenic enrichment.
According to Vecchia et al. (2017b) there is a great variation of Pb isotopic
Table 3. Loading values of each variable in each of the principal components for samples of water, sediment and granite of the uraniferous province of Lagoa Real BA.
signatures for both sediment and groundwater samples related to the geogenic sources (representative of the regional lithologies). However, the variation for the sediments is significantly higher. Moreover, among the geogenic sources it was observed that the Pb isotopic signatures for sediments and groundwater, from the same geographical coordinate, indicated a natural enrichment in U and Th. The values of the majority of the Pb isotopic signatures for groundwater are considered low, in spite of being in accordance with the geogenic sources of the region, which are under the influence of granites. One sediment sample (P35S), collected in the U mining area, presented isotopic Pb signature typical of U anthropogenic sources (high 206Pb/204Pb ratio compared to 208Pb/204Pb ratio), which means U contamination. However, the PCA analysis proved that this was not the behavior evidenced by the P19S sample, which presented isotopic Pb signature in accordance with the geogenic sources of the region.
One of the major vulnerabilities of the nuclear area and/or mining activities to public opinion refers to frequent suspicions that nuclear and/or mining industry always generate(s) large environmental impacts, especially contamination of groundwater. These kinds of worries can find legal resonance even without the necessary pre-operational background data, and most of the decisions consider the analytical results of radionuclide concentrations from the operational phase as the main incriminating technical argument. The present study showed that the determination of the concentration of the contaminants is not enough to characterize environmental contamination. It also shows that the most effective analytical method to elucidate the role of anthropogenic interferences in the environment is the determination of isotopic ratios. It can even be used in environments with high concentrations of radionuclides.
This tool becomes even more relevant in regions where pre-operational data are scarce and incomplete, if not inexistent, or even in regions, such as Caetité (Bahia), where water sampling is very difficult due to low water availability. In addition, it has been proved that the Pb isotopic tool, besides being used in any type of environmental matrices, is still promising in the investigation of industrial environmental impacts from mining as well as from physical, chemical and thermal processes where the generation of NORM (Natural Occurring Radioactive Material) wastes is related.
Authors are thankful to CNEN for Doctoral scholarship; to CNPq and to FAPEMIG for financial support.
 Bollhöfer, A., & Rosman, K. J. R. (2000). Isotopic Source Signatures for Atmospheric Lead: The Southern Hemisphere. Geochimica et Cosmochimica Acta, 64, 3251-3262.
 Bollhöfer, A., & Rosman, K. J. R. (2001). Isotopic Source Signatures for Atmospheric Lead: The Northern Hemisphere. Geochimica et Cosmochimica Acta, 65, 1727-1740.
 Cavalcante, Y. L., Hauser-Davis, R. A., Saraiva, A. C. F., Brandão, I. L. S., Oliveira, T. F., & Silveira, A. M. (2013). Metal and Physico-Chemical Variations at a Hydroelectric Reservoir Analyzed by Multivariate Analyses and Artificial Neural Networks: Environmental Management and Policy/Decision-Making Tools. Science of the Total Environment, 442, 509-514.
 CETESB. Companhia de Tecnologia de Saneamento Ambiental. (2011). Guia nacional de coleta e preservação de amostras: água, sedimento, comunidades aquáticas e efluentes líquidos. 2. ed. São Paulo: CETESB; Brasília-DF: Agência Nacional de águas, ANA, Ministério do Meio Ambiente, 325 p.
 Chen, J., Tan, M., Li, Y.-L., Zhang, Y., Lu, W., Tong, Y., Li, Y., & Zhang, G. (2005). A Lead Isotopic Record of Shangai Atmospheric Led Emissions in Total Suspended Particles during the Period of Phasing out of Leaded Gasoline. Atmospheric Environment, 39, 1245-1253.
 Cordani, U. G., Iyer, S. S., Taylor, P. N., Kawashita, K., Sato, K., & Mcreath, I. (1992). Pb-Pb, Rb-Sr, and K-Ar Systematic of the Lagoa Real Uranium Province (South-Central Bahia, Brazil) and the Espinhaço Cycle (ca. 1.5-1.0 Ga). Journal of South American Earth Sciences, 5, 33-46.
 Gioia, S. M. C. L., Babinski, M., Weiss, D. J., & Kerr, A. A. F. S. (2010). Insights into the Dynamics and Sources of Atmospheric Lead and Particulate Matter in São Paulo, Brazil, from High Temporal Resolution Sampling. Atmospheric Research, 98, 478-485.
 Gioia, S. M. C. L., Pimentel, M. M., Tessler, M., Dantas, E. L., Campos, J. E. G., Guimarães, E. M., Maruoka, M. T. S., & Nascimento, E. L. C. (2006). Sources of Anthropogenic Lead in Sediments from an Artificial Lake in Brasília, Central Brazil. Science of the Total Environment, 356, 125-142.
 Gulson, B. L., Korsch, M. J., Dickson, B., Cohen, D., Mizon, K. J., & Davis, M. J. (2007). Comparison of Lead Isotopes with Source Apportionment Models, Including SOM, for Air Particulates. Science of the Total Environment, 381, 169-179.
 Iyer, S. S., Babinski, M., Marinho, M. M., Barbosa, J. S. F., Sato, I. M., & Salvador, V. L. (1999). Lead Isotope Evidence for Recent Uranium Mobility in Geological Formations of Brazil: Implications for Radioactive Waste Disposal. Applied Geochemistry, 14, 197-221.
 Komárek, M., Ettler, V., Chrastny, V., & Mihaljevic, M. (2008). Lead Isotopes in Environmental Sciences: A Review. Environment International, 34, 562-577.
 Lamego Simões, F. F. F., Fernandes, H. R. S. M., Franklin, M. R., Flexor, J. M., Fontes, S. L., Pereira, F. S. R., & Nascimento, F. M. F. (2003). Impactos de Mineração e Sustentabilidade no Semi-árido.Estudo de caso: Unidade de Concentração de Uranio—URA (Caetité, BA). Centro de Tecnologia Mineral—CETEM—Ministério da Ciência e Tecnologia. In: Simpósio Brasileiro de Recursos Hídricos, 15, 2003, Curitiba. Comunicaçãotécnica. Porto Alegre, ABRH, 18 p.
 Licht, O. A. B. (2001). Geoquímica Multielementar na Gestão Ambiental. Identificação e Caracterização de Províncias Geoquímicas Naturais, Alterações Antrópicas das Paisagens, áreas Favoráveis à Prospecção Mineral e Regiões de Risco no Estado do Paraná, Brasil. Tese de Doutorado. Curitiba: Universidade Federal do Paraná, 209 p.
 Siepak, M. & Sojka, M. (2017). Application of Multivariate Statistical Approach to Identify Trace Elements Sources in Surface Waters: A Case Study of Kowalskie and Stare Miasto Reservoirs, Poland. Environmental Monitoring and Assessment, 189, 364.
 Sun, C., Liu, J., Wang, Y., Sun, L., & Yu, H.-W. (2013). Multivariate and Geostatistical Analyses of the Spatial Distribution and Sources of Heavy Metals in Agricultural Soil in Dehui, Northeast China. Chemosphere, 92, 517-523.
 Vecchia, A. M. D. (2015). Investigações Isotópicas de Pb em águas Subterraneas e Sedimentos—Província Uranífera de Lagoa Real (Bahia). Tese de Doutorado. Centro de Desenvolvimento da Tecnologia Nuclear. Comissão Nacional de Energia Nuclear. Belo Horizonte, 222 p.
 Vecchia, A. M. D., Rodrigues, P. C. H., & Ladeira, A. C. Q. (2017a). Assinaturas isotópicas de Pb na avaliação do impacto ambiental causado pelo lançamento do efluente tratado pelo complexo mínero industrial do planalto de Poços de Caldas na represa de águas Claras e Ribeirão das Antas. Geonomos, 25, 14-23.
 Vecchia, A. M. D., Rodrigues, P. C. H., Ladeira, A. C. Q., & Rios, F. J. (2015). Interpretação de dados isotópicos de Pb em diferentes ambientes investigativos visando o diagnóstico de fontes geogênicas e/ou antrópicas. Geonomos, 22, 77-90.
 Vecchia, A. M. D., Rodrigues, P. C. H., Rios, F. J., & Ladeira, A. C. Q. (2017b). Investigations into Pb Isotope Signatures in Groundwater and Sediments in a Uranium-Mineralized Area. Brazilian Journal of Geology, 47, 147-158.
 Veridiana, M., Babinski, M., Ruiz, I., Sato, K., Souza, S., & Hirata, R. (2008). Analytical Procedures for Determining Pb and Sr Isotopic Compositions in Water by ID-TIMS. Química Nova, 31, 1836-1842.
 Zhao, J., Fu, G., Lei, K., & Li, Y.-W. (2011). Multivariate Analysis of Surface Water Quality in the Three Gorges Area of China and Implications for Water Management. Journal of Environmental Sciences, 23, 1460-1471.