Tracing geochemical sources and health risk assessment of uranium in groundwater of arid zone of India

Water quality degradation and metal contamination in groundwater are serious concerns in an arid region with scanty water resources. This study aimed at evaluating the source of uranium (U) and potential health risk assessment in groundwater of the arid region of western Rajasthan and northern Gujarat. The probable source of vanadium (V) and fluorine (F) was also identified. U and trace metal concentration, along with physicochemical characteristics were determined for 265 groundwater samples collected from groundwater of duricrusts and palaeochannels of western Rajasthan and northern Gujarat. The U concentration ranged between 0.6 and 260 μg L−1 with a mean value of 24 μg L−1, and 30% of samples surpassed the World Health Organization’s limit for U (30 μg L−1). Speciation results suggested that dissolution of primary U mineral, carnotite [K2(UO2)2(VO4)2·3H2O] governs the enrichment. Water–rock interaction and evaporation are found the major hydrogeochemical processes controlling U mineralization. Groundwater zones having high U concentrations are characterized by Na–Cl hydrogeochemical facies and high total dissolved solids. It is inferred from geochemical modelling and principal component analysis that silicate weathering, bicarbonate complexation, carnotite dissolution, and ion exchange are principal factors controlling major solute ion chemistry. The annual ingestion doses of U for all the age groups are found to be safe and below the permissible limit in all samples. The health risk assessment with trace elements manifested high carcinogenic risks for children.

www.nature.com/scientificreports/ F − and arsenic (As) contamination are reported in several studies [9][10][11] . While WHO (1999) has set a provisional guideline of 30 μg L −1 for U concentration in drinking water, and the Atomic Energy Regulatory Board (AERB) 12 in 2004 has set up a radiological limit of 60 μg L −1 , no such limit has yet been set by the Bureau of Indian Standards (BIS) 13 . The upper limit of concentration of heavy metals and physicochemical parameters in groundwater has also been limited by the WHO 14 and the BIS 13 . It is worthy to mention that the hydrogeochemical approach is essential to delineate the origin of U mineralization. Besides this, the role of hydrogeochemistry and geology have already been explored in the arid regions related to the evapotranspiration process in the tropical areas of Langer Heinrich (Southern Hemisphere), Klein Trekkopje (Namibia), Yeelirrie (Western Australia), Lake Maitland (southern Argentina) and Latin America [15][16][17][18][19] . The seasonal variations and the process of Uranium mobilisations in hard rock aquifers along with their health hazard indices in the south Indian aquifer indicated the redox reaction as a major governing factor 20,21 . However hydrogeochemical and U databases in groundwater of the palaeochannels in the chemical delta distributed in parts of Thar Desert of the western Rajasthan and northern Gujarat are meagre. Therefore, an attempt is made to create a geochemical database of the heavy metals and physicochemical parameters in the groundwater of the arid region of Rajasthan and Gujarat. The study becomes significant in itself because the groundwater serves as the main source for domestic and irrigation purposes in this region. The objectives of this study are as follows: (1) to identify the major hydogeochemical process associated with U mineralization in calcrete groundwaters along the palaeochannels of the Luni river system; (2) to identify the main source of U; and (3) to assess the spatial distribution and health risk assessment of U. Results of this study would be useful in framing policies to regulate water quality.

Material and methods
Study area and geology. The study area is bounded by latitudes 24°-29° N and longitudes 70°-76° E, covering the western part of Rajasthan and northern part of Gujarat, India. The sampling areas included the palaeodrainages of Luni, Khari bani, Jojori river, Sukri, Mithri, Sagi, and Guhiya rivers in Nagaur, Pali, Jalore, Sirohi, and Banaskantha districts of Rajasthan and Gujarat. Samples were collected from palaeochannels and groundwater from wells, tube-wells and taps. The sampling area falls in the arid and semi-arid region of the Thar desert along the carnotite mineralization area. The average rainfall in Pali district is 462 mm; whereas in Nagaur district, it is 310 mm. Barmer lies in the arid region with an average annual rainfall of 260 mm and potential evaporation of 1857 mm. In Jalore and Sirohi districts, the average rainfall is 280 mm and 606 mm. Luni is the major drainage system in western Rajasthan with many tributaries like Sukri, Sagi, Bandi, and Jawai. Guhiya river originates near Khariyaniv and Tharasani villages in Pali District in the hillocks. These palaeochannels drain through Malani Igneous Suite (MIS) especially Jalore granites and are located in the northwestern part of Rajasthan. These granites are adjacent to the paleochannels of late Tertiary to Quaternary settings; while a small portion belongs to Neoproterozoic, Paleoproterozoic, and Mesoproterozoic age 22 . The water table depths ranged from 50 to 150 m in the current study area 23 . In India, a favourable geological and climatic setup is evident in western Rajasthan and Northern Gujarat for the formation of surficial-type of U mineralization. The area under investigation is mostly covered with Quaternary formations concealing Sirohi, Erinpura, and MIS rocks. Calcrete occurrences along western Rajasthan and Northern Gujarat are well-known due to their arid climatic conditions, and these are indicators of palaeo-climatic conditions. The vast development of calcretes spreads over 0.32 million m 2 . Fluorite (CaF 2 ) and Barite (BaSO 4 ) occur in the form of veins in volcanic agglomerates. The Department of Mines and Geology, Government of Rajasthan has estimated fluorite reserves of 0.17 million tonnes in the Jalore area.
Groundwater sampling and chemical analysis. The total area was divided into a grid pattern of size 8 × 8 km 2 to cover the entire area of each district. The sampling location and places of U occurrences are shown in Fig. 1a. During the process of investigation, 265 samples were collected from the channels making the sample density of the order of one sample per square kilometre. The sample collection and storage were done according to the standard protocols described by ISO (2012). The groundwater samples were collected in plastic bottles after having been filtered through 0.45-μm millipore filter paper and acidified with 2-M HNO 3 (ultrapure merck) for cation and anion analysis. All the chemicals used for analysis were of analytical grade. pH was determined by pH-meter (ELICO, India) and conductivity by conductivity-meter. The pH and EC were analysed in the field. The acid-base titration method was used to analyse bicarbonate ions (HCO 3 − ), and sulphate ions (SO 4 2-) were determined by turbidimetry. Chloride ions (Cl − ) were determined using the volumetric method. F − was measured using the Fluoride electrode. The cations, calcium (Ca 2+ ), and magnesium (Mg 2+ ) were analysed using complexometric titration using ethylenediamine tetra acetic acid (EDTA.) The heavy metals zinc (Zn), copper (Cu), lead (Pb), lithium (Li), vanadium (V) and cobalt (Co) were detected using inductively coupled plasma atomic emission spectrometry. Three replicates were used to analyze all of the samples. As part of quality control, duplicate and standard checks were performed on every ten samples. In addition, a trace element standard reference material (SRM-1643f.) was examined. Supplementary material for methodology (Table S1a) detail the instrument's operating and optimized conditions, calibration, and QA/QC details.
The U content in the samples was measured by a LED Fluorimeter (Quantalase Enterprise Pvt. Ltd. model LF-2). It works on the principle of measurement of fluorescence of uranyl complex in the water sample. Here a pulsed LED UV light was used to excite uranyl species at 405 nm. The detection range of the instrument is 0.5-1000 ppb, with a detection limit of 0.2 ppb. Na + and K + were analysed using ELICO made flame photometer (Model No. 341). The accuracy and reliability of the water analysis were checked using the charge balance approach (error calculation) as follows:

Results
The present arid area is marked by carnotite mineralization where V and U are the prominent elements. Various geological factors influence the concentration of U in water. Enhanced concentration of bicarbonate in groundwater augments the formation of highly soluble uranyl carbonate complexes thereby resulting in the elevated U concentration in groundwater. U in groundwater in India is generally reported in the range of 1-100 µg L −1 . In the Indo-Gangetic plain, U has values ranging from 0.6 to 65.3 µg L −1 (Amritsar), 0.12-38.8 μg L −1 (Pathankot) and 0.12-38.8 μg L −1 in Gurdaspur 24 . Kumar et al. 25 have reported the value of U to range from 0.2 to 644 μg L −1 in Bhatinda (Punjab). In Rajasthan, Ganganagar district, the values of U varied from 2.5 to 171 μg L −1 . Whereas, the values in Sikar district were between 4 and 136 μg L −1 . Meanwhile, the occasional high occurrence of PO 4 3− levels (5.90 mg L −1 ) may be related to anthropogenic activities such as agriculture activities, which can also be an additional source of U in groundwater 26 . This contamination is due to the use of phosphatic fertilizers resulting in the additional contamination of groundwaters with phosphate 27 .
The statistical parameters of the physicochemical elements, U and heavy metals and standard ranges for defining the suitability of these elements in water for drinking purposes are listed in Table 1. The U concentration ranged between 0.6 and 260 μg L −1 with a mean value of 24.0 μg L −1 . U concentration in 30% of the 265 groundwater sources exceeded the provisional WHO (2012) health guidelines. Results show that 69.3% of samples have U values below 30 ppb; 27.3% of samples shows values between 30 and 60 ppb; 3% of the samples indicate values between 60 and 100 ppb. Moreover, 0.38% of samples have values greater than 100 ppb. The spatial distribution of U is shown in Fig. 1b. The highest value of U is observed in the piedmont zone of Sagi and Sukri rivers. It is worth noting that MIS suites of rocks are represented by rhyolite, volcanic agglomerate/breccia, andesite, and basalts. HCO 3 -is in moderate amount, ranging from 101 to 816 mg L −1 . SO 4 2− was increasing in the northeast (NE) direction; while K + , CO 3 2− , and PO 4 3− exhibited a similar distribution pattern, having low values. Spatial distribution is essential for evaluating spatiotemporal evaluation of different important parameters 28,29 . The spatial distribution of heavy metals and trace metals is indicated in Figs. S1 and S2. The spatial distribution of metals, trace elements and their origin is given in supplementary text l S3.
The higher salinity water is identified to have a high U concentration in the groundwater. Saline influence in the samples is calculated using the base exchange index (BEX). BEX can be used to identify refreshing or seawater mixing processes.
A BEX around zero indicates freshwater; while a negative BEX points towards saltwater intrusion; and alternatively, a positive value refers to refreshing. The BEX index is found to be negative in only 23% of the samples.  www.nature.com/scientificreports/ Weathering plots are used to estimate the relative contribution of geogenic source (evaporite dissolution, silicate weathering and carbonate dissolution) to its ionic load 32 , (supplementary text S3). Figure 2a-g reveal the weathering plots for Na + versus Cl − , Ca 2+ + Mg 2+ versus HCO 3 − , (Na + + K + ) versus TZ + , (Ca 2+ + Mg 2+ ) versus TZ + , (SO 4 + HCO 3 ) versus (Ca 2+ + Mg 2+ ; (Na + + K + ) versus (Cl − + SO 4 2 ) and (Na + -Cl − ) versus (Ca 2+ + Mg 2+ )-(SO 4 2−+ HCO 3 − ) 33 . Principal component analysis (PCA) of the data makes it evident that water-rock interaction has played a major role in determining the hydrochemical characteristics of groundwater in the arid zone. PCA was performed as the initial step for the cluster analysis (CA) procedure. The Kaiser-Meyer-Olkin value (0.65) and the results of Bartlett's test of sphericity (p < 0.001) showed that PCA could be used for data analysis. Four principal components were derived based on the Kaiser criteria of eigenvalue > 1. The percentage of variance for the four components corresponded to 25.66, 10.60, 7.95, and 7.51%, respectively.
The equilibrium state of the water concerning the mineral phase can be determined by calculating the Saturation Indices SI using analytical data. The SI of a mineral could be obtained from the following equation.

Discussion
Hydrochemical facies and water type. The piper plot in Fig. S3 highlights the hydrochemical facies of U < 30 ppb, 30 ppb < U < 60 ppb and U > 60 ppb in the group TDS1, TDS2, and TDS3. The correlation analysis indicated a strong correlation between Na, Cl, EC, and TDS reflecting saline water impact. Consequently, groundwater salinity may have originated from water-rock interaction, seawater intrusion, evaporation, deposition, and fossil saline water 35 . The major hydrochemical facies observed in the different TDS group samples are given in supplementary text S3. It is observed that most of the high U content samples belonged to Na-Cl hydrochemical facies.
Principal component analysis.  (Table S3). Rock source deduction revealed carbonate weathering due to the dissolution of calcite and dolomite minerals leading to the increase in pH and TDS in groundwater. As observed by Prasanna et al. 36 , the clay lens enhances ionic concentration in the adjacent groundwaters due to the process of ion exchange. The saline groundwater with the association of Ca-Mg-SO 4 was found to facilitate the U mobilization. Further, the redox conditions due to salinity variations also promote the enrichment of U in groundwater. The long residence time of groundwater augments the salinity, and high salinity generally favours the release of U due to the increase in residence time. Principal component 2 (PC2) showed an association between HCO 3 − , V and U (Table S3), showing bicarbonate complexation and carnotite dissolution, leading to uranium enrichment. Carnotite is a secondary U mineral observed in calcrete deposits. Carnotite is generally reported in a region with U-enriched groundwater, along with the SO 4 2− and CO 3 2− facies formation, especially in arid to semi-arid climates 27 . Other factors governing the carnotite formation process are pH, Eh, pCO 2 , availability of HCO 3 , U, K etc., in groundwater, or the rate of groundwater flow, alteration process, and stagnation of water, either by chemical or physical process. Groundwater is characterised by the dissolution of carnotite predominantly containing U, V, and K + (Table S3). As ground- www.nature.com/scientificreports/ water flows through a calcrete formation, it dissolves HCO 3 and adjoining clays exchange K + for hydrogen which leads to the increase of HCO 3 − and decrease of K and pH 27 . The dissolution of carnotite in the presence of CO 2 and water leads to the formation of K, U, V, and HCO 3 − Association of Pb, Ca 2+ and K + in the PCA 3, along with the negative correlation of pH represented 7.95% of the variance in PCA 3. The K + in the aquatic medium may be exchanged for H + in water, indicating the ion exchange process. The negative correlation of pH revealed that the cation exchange is the predominant process. PCA 4 explaining 7.51% of the variance indicated a positive correlation with heavy metals F − and Ni and a negative correlation of Li and Pb. The positive association between Ni and F − corresponded to their geogenic occurrence. Thus, PCA indicated that carnotite dissolution and rock water interaction including the evaporite dissolution are the processes responsible for the presence of U. Carnotite formation and speciation of U got augmented due to the leaching by meteoric waters from uranium-rich igneous rocks of the MIS. During their migration, the waters become more alkaline and saline due to progressive interactions and evaporation. Groundwater is intercepted by basement lithology and is directed towards the surface through soil suction/ capillary rise where ion exchange and evaporation increases, CO 2 degassing-pH decreases and change in water salinity, resulting in the breaking down of uranyl carbonate ions 20 . Uranyl carbonate is extremely stable in aqueous solutions, but it destabilizes due to evaporation. Decomplexation of uranyl carbonate in waters of higher salinities and evaporation converts into uranyl ion, which forms as nucleating seed over clay or hydroxide of iron and aluminium 37 . Subsequently, in the presence of V and K + ions, carnotite precipitation takes place in the pH range of 6.0-8.0. Given that higher V in groundwater is rare, V may be found present in alkali water of U-V mineralized regions. The higher V prohibits the movement of U 6+ in the pH between 4.0 and 8.0, leading to the thermodynamic stabilization of carnotite.
The predominant U species in groundwater was calculated by Eh-pH diagram. Stanley and Wilkin 38 argued that Eh-pH diagram constructed in the presence of hydroxide, sulphate, and carbonate ligands showed that calcrete waters were mostly characterised by neutral to alkaline and oxidising conditions favouring the U mobility through the formation of soluble anionic carbonate complexes, such as UO 2 (CO 3 ) 2 2− . Moreover, it is known that the mobility and absorption by phosphate, clay minerals organic materials, and iron oxides in an oxidising environment are governed by Eh and pH, While the adsorption dominates at low pH, the higher pH favours the formation of anionic uranyl carbonates. As early as 1978, Langmuir 39 described uranyl complexation and attributed it to pH and Eh as the governing factors of the process. In groundwater with higher Eh (oxidizing) and lesser pH (~ 5), the major forms of U are the uranyl ion UO 2 2+ as U 6+ . In the presence of F − , these waters lead to the formation of uranofluoride complexes. Subsequently, when the pH increases in these oxidising waters, the uranyl ion tries to accommodate the CO 3 2− to form the uranyl-carbonate complexes as UO 2 (CO 3 ) 2 2− . Correspondingly, at higher pH, it forms UO 2 (CO 3 ) 3 4−40 . The Eh-pH diagram for uranyl species (Fig. 3a) shows that the pH ranges from 6.9 to 9.2, and Eh ranges from 407 to 542 mV. The diagram shows that stability ranges for different forms of U concerning CO 3 2− and hydroxyl ion association with U, like U 4+ , UOH 3+ , UO 2 + , Uraninite and U (OH) −5 . The increase in pH shows the variation of uranyl form from UO 2 to UO 2 CO 3 , UO 2 (CO 3 ) 2 2− , UO 2 (CO 3 ) 3 4− . The availability of oxygen increases with Eh, subsequently, CO 3 2− increasing with pH. The samples in the UO 2 (CO 3 ) 2 2− and UO 2 (CO 3 ) 3 4− fields, reflecting the availability of CO 3 2− ions, and thus, increases uranyl carbonate species. The Eh-pH diagram of vanadium revealed the speciation of V. V mobility in natural waters mainly depends on its ability to form anion complexes. In the present study, HVO 4 2− and H 2 VO 4 − are the predominant V 5+ species in groundwater at pH ranges between 4.0 and 9.0 and indicates an oxidizing environment. The development of thermodynamic stable complexes of certain ligands with Vanadyl ions shows the persistence of species at pH > 6 (Breit and Wanty 1991) (Fig. 3b).
The order of dominance concerning the concentration of uranium carbonate species showed the trend UO 2 CO 3 < UO 2 (CO 3 ) 2 2− < UO 2 (CO 3 ) 3 4− reflecting the enhanced availability of carbonate ions (Fig. 3c). A parallel increase of CO 3 2− became evident concerning U ions in groundwater, leading to the complexation of uranylcarbonates. The concentration of uranyl sulphate species UO 2 (SO 4 ) 2 2− and UO 2 SO 4 were compared in groundwater samples (Fig. 3d). The ranges of sulphate species with U showed UO 2 SO 4 with a higher concentration than UO 2 (SO 4 ) 2 2− . While these species showed an increase in the concentration of these compounds. However, the samples are found to have been dispersed and not found to reflect a definite linear trend. The uranyl oxide (UO 2 + ) (Fig. 3e) also reflected a similar trend, ranging between 2.6 × 10 -20 and 9.3 × 10 -14 mol.
Saturation indices and water mineral equilibrium. The supersaturation of groundwater concerning calcium and magnesium bearing minerals suggested that there is a depletion of calcium and magnesium content due to their mineral precipitation in the groundwater system. In contrast, the samples were undersaturated with respect to evaporite minerals, such as anhydride, artinite, brucite, epsomite, fluorite, gypsum, and halite.
These results indicate that the mineralization is linked to the dissolution of evaporitic minerals (halite, anhydrite, and gypsum). www.nature.com/scientificreports/ All these minerals are complexes of UO 2 (Uraninite). The association of this complex with the hydroxyl molecule leads to schoepite, and with carbonates, it leads to the saturation of rutherfordine. The Ca(PO 4 ) 2 and KVO 4 association results in saturation states of autunite and carnotite, respectively. The calcrete observed inland could be mainly due to the carbonate rock precipitation in the phreatic zone of an aquifer 41 . The presence of calcrete in this zone facilitates the interaction between the laterally moving groundwater either by ion exchange or by chemical reaction process. The regions with drastic water level fluctuation, especially in the structurally disturbed hard rock aquifers favour calcrete deposits, and if the drainage is known to be rich in U, the carnotite minerals along these planes can be expected. This generally happens when the U from the granitic formation like MIS is released and combines with the V, derived from the mafic rocks. The MIS consists of complex acidic and is intrusive. V was also observed to be in Fe silicate minerals, like biotite and amphibole as a replacement for Fe in the crystal lattice 18 . An earlier investigation showed a similar occurrence of calcrete deposits in Australia having U associated with V 42 . Again, even much earlier than the year 1999, it was reported that the higher oxidation state of V facilitates the precipitation of U from the solution. The thermodynamic state of the carnotite composition is in near saturation conditions, even while other mineral compositions like uraninite is found undersaturated Fig. 3f. The stagnated water conditions favour evaporation, change in pCO 2 and pH, and subsequent saturation of waters with carbonate and calcium that ultimately help in the precipitation of calcium carbonates (calcrete).  Table 2. The greatest mean annual ingestion dose is observed in infants (0-6 mo) and (7-12 mo) were 49.30 and 56.34 µSvy −1 in comparison to all other age groups. Annual ingestion dose due to U in drinking water for males in age group 9-13, 14-18, and > 18 years ranged between 0.74 to 119.13; 1.01 to 161.40 and 0.22 to 35.11 µSvy −1 , respectively. For females, it varied between 0.65 to 104.24; 0.70 to 112.49 and 0.16 to 25.62 µSvy −1 in the same order of age groups ( Table 2). The annual effective dose due to the ingestion of U was found to be less than 100 µSvy −1 as per the European Commission (Table 2). Thus, the values are found within the radiological risk.
Radiological and chemical toxicity risk. Radiological risk assessment due to U ingestion was calculated by estimating the cancer risk. Table 2. There are three naturally occurring isotopes of U with half-lives of 2.4 × 10 5 , 7.0 × 10 8 , and 4.5 × 10 9 years. The cancer mortality risk coefficients Bq L −1 for the three isotopes of U 234 U, 235 U, 238 U were taken to be 6.1 × 10 −11 , 6.2 × 10 -11 , and 7.5 × 10 −10 , respectively. The average values for cancer mortality risk by U isotopes 234 U, 235 U, and 238 U are found 4.3 × 10 -11 , 2.09 × 10 -10 , and 3.52 × 10 -11 , respectively. Similarly, the mean cancer morbidity risk for 238 U, 235 U, and 234 U are 4.34 × 10 -10 , 5.71 × 10 -11 , and 5.66 × 10 -11 , respectively. U is a nephrotoxin that may result in kidney damage. Similarly, the chemotoxic dose was also computed for estimating health hazards. The lifetime average daily dose (LADD) found with a minimum and maximum value between 0.011 and 5.24 µg kg −1 d −1 , with a mean value of 0.47 µg kg −1 d −1 , which is much less than the reference dose (RfD) of 0.6 and 4.4 µg kg −1 d −1 prescribed by WHO 43 and AERB 12 , respectively. Also, the hazard quotient (HQ) shows arithmetic mean (AM), median and inter-quartile range (IQR) of 0.78, 0.64, and 0.71, respectively. Since HQ is less than unity, the analysed water samples could be used for drinking thus showing no considerable radiation risk to the population of the area under investigation (Table 2). Doses to various organs/tissues of an adult human are summarized in Table S4. Po-210 was seen to have been the major contributor to the dose to all organs, except for bone surface, for which Pb-210 is the source. The annual effective ingestion dose to the whole body was observed to have spanned between 0.114 and 41.791 µSv with an average of 4.710 µSv, which is of much less value than the recommended limit of 100 µSv 43 . Bone surfaces has the maximum share of dose (38%) Table 2. Radiotoxicity and chemical toxicity due to U (here U-238, U-235 and U-234 are isotopes of Uranium) in groundwater (DWI = Daily water intake). www.nature.com/scientificreports/ due to U and its daughters, followed by kidneys (14%), LLI Wall (12%), liver (5%), ULI Wall (5%), and small red marrow (4%). The dose coefficients calculated with the hair model were seen to be lower than ICRP's 44 biokinetic model because a considerable fraction of U in the blood (about 31.5%) is known to have excreted into the hair, and therefore resulting in no contribution to dose. This study elaborates the behaviour of uranium in groundwaters of complex terrain, especially in arid regions. The association of U to major ions and metals helps to understand the nature of reactions in groundwater of similar climatic conditions. The health risk assessment and identification of regions with higher U in groundwater helps the policymakers to manage the utility of the available water resources.
Classification of groundwater-based heavy metal index. The metal toxicity load (MTL) index was developed to find the toxicity of various metals in the groundwater samples (Table 3). From the results, it was found that 18.67% of sampling sites showed that toxicity of Ni was below their permissible toxicity load (825 mg L −1 ); while 81.33% of sites exceeded their permissible toxicity load. Similarly, 1.05% sampling sites exceeded the toxicity of Zn for permissible toxicity load, which is 136.95 mg L −1 . The 3.87% sampling sites went above toxicities of Pb for permissible toxicity load, which is 76.55 mg L −1 ; 27.11% sampling sites showed to have gone above toxicities of U for permissible toxicity load, which is 24.99 mg/L. Their percentage removal varied from 3.22% to 88.37%. Nonetheless, for F, V, Co, Cu, and Pb, their toxicity range showed to be below the permissible toxicity load. On the whole, the range of toxicity for metals was found above the permissible toxicity load (315. 19). The heavy metal pollution index (HPI) index was calculated for the metals of groundwater (Table 3). From the results of HPI, it was observed that HPI value was found to be less than 100, signifying less pollution in the study area.
Health risk assessment due to heavy metals. To understand the non-carcinogenic risk associated with heavy metals, an average daily dose of exposure CDI ing and CDI derm , Hazard Quotients, HQ ing and H Qderm, and Hazard Index (HI) corresponding to ingestion and dermal pathways were calculated separately for adults and children as tabulated in Table 4. The maximum daily dose in case of adults belonged to Li, i.e., 2.139 mg/kg/d in drinking water. The CDI ing in the case of adults followed the order Li . The HQ ing children followed the trend Li > Co > Pb > V > Ni > Cu > Zn. In the case of both adults and children, the HQs of Li contributed predominantly to the integrated HQ values due to higher concentrations among the metals considered. Since the average HI value for adults (0.965) showed to be less than unity and for the children (1.446) to be greater than unity; consequently, a significant non-carcinogenic risk was diagnosed for children. The mean carcinogenic risk was observed to be 7.62 for children, conversely, 5.08 for adults. The cancer index values for Ni and Pb were found higher than the ranges recommended by USEPA 45 of 1 × 10 -6 and 1 × 10 -4 for both adults and children (Table 4).

Conclusion
Hydrogeochemistry and geochemical speciation is carried out to investigate the mechanism responsible for the enrichment of U in the calcretes of the arid region of western Rajasthan. The non-carcinogenic and carcinogenic risks of the trace metals and uranium in the groundwater were also evaluated for health risk assessment. Carnotite dissolution is identified to be the primary cause of U enrichment in calcrete waters. High U concentration in groundwater presented a distinctive hydrogeochemical characteristic as follows: high TDS and alkalinity Table 3. Metal toxicity load of the groundwater following ATSDR (2017) relative to the toxicity level of metals accountable for human beings and computation of heavy metal pollution index (HPI) for groundwater of western Rajasthan and Northern Gujarat. www.nature.com/scientificreports/ along with high Na + and K + concentration, suggesting that weathering of geogenic source material, evaporation, and dissolution from mineral surfaces are the principal mechanism of U release. Based on the TDS values, the groundwater samples were divided into TDS < 1000 mg L −1 (Ca-MgHCO 3 type), 1000 < TDS < 2000 mg L −1 and TDS > 2000 mg L −1 (Na-Cl) type of water. High evaporation rates, typical of arid and semiarid climates generate saline groundwaters and alkaline pH, releasing U and F − from both the primary and the secondary material sources. Despite being a U-mineralized area, only 30% of samples exceeded their recommended WHO, since carnotite was seen to have been least affected by the dissolution in the low pCO 2 in desert soils and groundwater due to the absence or paucity of organic activity in the soil. Chemical speciation, computed using PHREEQCI of uranium/carbonate system, indicated that the predominant U species in alkaline conditions was UO 2 (CO 3 ) 3 4− . HVO 4 2− and H 2 VO 4 − were the predominant V 5+ species in groundwater in the completed study. Hydrogeochemical characteristics indicated that ionic load resulted from halite dissolution, silicate weathering, and ionic exchange. The metal toxicity load index suggested that the range were above the permissible toxicity load. The heavy metal pollution index obtained in this study is 56.34, signifying moderate pollution of heavy metals in the groundwater quality of the study area. The non-carcinogenic risks exhibited that Li, V, and Ni are the key pollutants affecting the health of human beings; conversely, carcinogenic risks inferred that Pb showed high health risks. However, the concentration of V, Li, Co, Cu, Zn, and Pb were found to be lower than those recommended by BIS 13 . 20% of samples surpassed their Ni limits. The aforesaid findings confirm the crucial role of hydrogeochemical processes in heavy metal enrichment in calcrete water. The results will be highly useful in understanding processes controlling heavy metal enrichment in arid waters.