ABSTRACT
Ensuring the availability of high-quality water is a critical challenge. This study evaluates groundwater quality in the Palar basin, Tamil Nadu, India, during 1973, 1983, 1993, 2003, 2013, and 2022, focusing on pre- and post-monsoon seasons. Key irrigation indices – sodium absorption ratio (SAR), residual sodium carbonate (RSC), sodium percentage (Na%), permeability index (PI), magnesium hazard (MH), Kelly's ratio (KR), and potential salinity (PS) – were assessed using parameters such as calcium (Ca), magnesium (Mg), chlorine (Cl), potassium (K), and sulfate (SO4). pH, electrical conductivity (EC), SAR, Na%, RSC, MH, and water quality index (WQI) showed significant spatial and temporal variability. EC ranged from 250 to >3,000 μS/cm, with industrial hotspots like Ranipet exhibiting high levels. SAR was mostly ‘excellent’ (0–10) but exceeded 18 in 1993, indicating irrigation unsuitability. Na% values >60% degraded soil permeability, while RSC and MH exceeded thresholds near industrial zones, impairing soil and irrigation quality. WQI classified most samples as ‘none’ to ‘slight,’ but hotspots >450 posed consumption risks. Seasonal variations showed post-monsoon runoff effects, increasing Na and SO4. GIS mapping and statistical analyses highlighted contamination hotspots. Mitigation requires stricter regulations, improved wastewater treatment, groundwater recharge, and sustainable agriculture to ensure long-term water quality.
HIGHLIGHTS
Groundwater was analyzed from 1973 to 2022 during pre- and post-monsoon periods.
Indices, such as sodium absorption ratio, residual sodium carbonate, sodium percentage, permeability index, magnesium hazard, Kelly's ratio, and potential salinity, were calculated.
Geographic information system was used to map and compare water quality concentrations with standard values.
High WQI: Water quality index results show that groundwater quality is generally low.
Groundwater in the study area is largely safe for both drinking and irrigation.
INTRODUCTION
Water is essential for life, and groundwater is a critical source of freshwater, especially in arid and semi-arid regions. However, the over-exploitation of groundwater due to easy access and increasing energy availability has led to significant declines in groundwater levels over the past two decades. By 2025, global demand for drinking water is expected to rise by 44%, irrigation by 10%, and industrial water use by 81% (Burri et al. 2019). In India, issues related to groundwater quality are becoming increasingly severe, especially in water-scarce states like Tamil Nadu, where contamination from pollutants like fluoride poses major challenges. Studies conducted on the Palar River, one of the most contaminated rivers in India due to industrial pollution, highlight the water quality challenges in this region, with fluoride concentrations found between 0.18 and 0.22 mg/L, below the World Health Organization (WHO) 1.5 mg/L level (Kumar et al. 2019). Factors such as rapid population growth, industrialization, and pollution continue to exacerbate groundwater stress, threatening its availability for both drinking and agricultural purposes.
The Palar River has been heavily polluted by effluents from tanneries, particularly in Vellore District, where elevated chromium concentrations pose serious environmental and public health risks. This pollution has rendered the river unsuitable for consumption and agriculture, leading to widespread health issues, such as asthma, skin diseases, and stomach ailments, while thousands of acres of fertile land have been abandoned. In 2006, Ranipet, an industrial area along the river, was labeled one of the most polluted places in the world due to heavy contamination by salts and metals, including chromium (Pranavam et al. 2011).
Rapid urbanization in India has added complexity to the challenge of maintaining both the quantity and quality of groundwater. As a substantial portion of water demand is met by groundwater, urban growth and intensive agriculture continue to deteriorate water quality. Climate change further compounds these challenges by altering precipitation patterns, intensifying pressure on water resources (Awodumi & Akeasa 2017). Groundwater contamination, often due to wastewater discharge, poses significant risks to public health, and its detection is more challenging than surface water pollution (Zektser & Everett 2004; Singh & Kumar 2017).
To address these concerns, geographic information system (GIS) tools offer a powerful method for analyzing spatial patterns and trends in groundwater quality. GIS-based techniques help to monitor contamination and identify pollution sources, enabling informed decision-making for water resource management (Balakrishnan et al. 2011). In a study conducted in the flood plains of the upper Palar River, the analysis of 22 groundwater samples during the pre- and post-monsoon seasons revealed tannery effluents, solid wastes, and sewage as primary sources of pollution. Physico-chemical analysis indicated high alkali and sulfate (SO4) content, with statistical analysis supporting the connection to industrial pollution (Kuppuraj et al. 2012).
Similar studies in India have assessed groundwater quality using GIS and water quality indices (WQIs). Ram et al. (2021) mapped groundwater quality in Mahoba District, Uttar Pradesh, using ‘inverse distance weighted’ (IDW) interpolation on 43 samples. They found most of the area had acceptable quality, with some pockets of poor quality, highlighting the usefulness of WQI for assessing groundwater suitability. Dandge & Patil (2022) conducted a similar study in Bhokardan, Maharashtra, using GIS and remote sensing to analyze 12 water quality parameters. Their findings showed poor water quality in the pre-monsoon season, with improvements post-monsoon, emphasizing the importance of seasonal monitoring and WQI for groundwater management.
WQIs are essential tools for assessing water quality and are particularly useful for monitoring both drinking water and irrigation water quality. Various indices such as sodium absorption ratio (SAR), residual sodium carbonate (RSC), sodium percentage (Na%), permeability index (PI), and Kelly's ratio (KR) provide valuable insights into the suitability of water for irrigation purposes, and their application in the study of groundwater contamination is critical for understanding its impact on agriculture (Jafari et al. 2018; Khalid 2019). Regular monitoring and evaluation of these parameters are necessary for effective water resource management (Prathumratana et al. 2008). Al-Amry (2008) assessed groundwater quality in the Al-Salameh area, Yemen, using 25 samples. The results showed that 17 samples were suitable for drinking, while irrigation suitability was determined based on parameters like SAR, KR, and RSC, indicating general suitability for irrigation.
Despite the growing body of research on groundwater contamination in the region, significant gaps remain in understanding the full extent of pollution and its long-term effects on human health and the environment. This study aims to fill these gaps by assessing the quality of groundwater in the Palar River Basin, focusing on the impact of industrial pollution and the suitability of groundwater for drinking and irrigation purposes.
The primary objectives of this investigation are as follows:
1. To visualize the spatial distribution of irrigation indices using GIS techniques across the study area.
2. To compare spatial distribution maps of irrigation indices in different years to assess changes over time and to evaluate groundwater suitability for irrigation purposes based on the spatial analysis.
3. To support effective water resource management and sustainable agricultural practices by providing a comprehensive understanding of groundwater suitability trends.
STUDY AREA
The Palar Basin stands out as one of the most significant among Tamil Nadu's 17 basins. The area selected for this study is located between 12°37′19″ N–13°00′12″ N latitudes and 78°28′07″ E–79°34′51″ E longitudes with an areal extend of 10,163.151 km2 (Figure 1). The Palar Basin plays an essential role in Tamil Nadu's hydrology, agriculture, and water supply, making it highly relevant for assessing groundwater quality and irrigation suitability. As a transboundary basin shared with Karnataka and Andhra Pradesh, the Palar Basin supports a diverse range of agricultural activities and is a crucial source of groundwater for irrigation, drinking, and industrial use. Given its distinct geographical features – spanning the Eastern Ghats, plateau regions, and coastal plains – the basin experiences varied climatic and hydrological conditions, influencing groundwater recharge, quality, and availability.
By focusing on the Palar Basin, this study can reveal valuable insights into groundwater trends over time, providing a framework to manage and sustain water resources in a region where water demand is increasing and groundwater quality faces challenges from agricultural runoff and urbanization. The findings from this study will offer vital data to support informed decision-making for resource management, benefiting not only the Palar Basin but also contributing to broader water conservation and agricultural strategies across Tamil Nadu. The latitudes and longitudes of the sampling point in the study area are provided in Table 1.
Details of the study area
Well no . | District . | Taluk . | Village . | Latitude . | Longitude . | Subbasin . |
---|---|---|---|---|---|---|
1 | Kancheepuram | Madurandagam | Thennampattu | 12° 28′ | 79° 55′ 40″ | Kiliyar |
2 | Kancheepuram | Madurandagam | Velamur [ramapuram] | 12°28′40″ | 79° 46' | Kiliyar |
3 | Thiruvannamalai | Vandavasi | Tennangur | 12° 24′ 5″ | 79° 53′ 40″ | Kiliyar |
4 | Thiruvannamalai | Vandavasi | Salavedu | 12° 34′ | 79° 38′ | Kiliyar |
5 | Thiruvannamalai | Vandavasi | Peranamallur | 12° 28′ | 79° 46′ 32″ | Cheyyar |
6 | Thiruvannamalai | Vandavasi | Osur | 12° 34′ 17″ | 79° 25′ 59″ | Kiliyar |
7 | Thiruvannamalai | Thiruvannamalai | Nayadimangalam | 12° 26′ 15″ | 79° 40′ 10″ | Cheyyar |
8 | Thiruvannamalai | Polur | Kadaladi | 12° 24′ 18″ | 79° 6′ 6″ | Cheyyar |
9 | Thiruvannamalai | Polur | Mampattu | 12° 24′ 15″ | 78° 58′ 10″ | Cheyyar |
10 | Thiruvannamalai | Polur | Siruvallur | 12° 30′ 35″ | 79° 5′ 30″ | Cheyyar |
11 | Thiruvannamalai | Polur | Athuvambadi | 12° 27′ 45″ | 79° 1′ 20″ | Cheyyar |
12 | Thiruvannamalai | Polur | Melarani | 12° 35′ 50″ | 79° 7′ 35″ | Cheyyar |
13 | Thiruvannamalai | Chengam | Chengam Farm 2 | 12° 27′ 20″ | 79° 2′ 45″ | Cheyyar |
14 | Vellore | Katpadi | Vaduganthangal | 12° 16′ 20″ | 78° 43′ 10″ | Poiney (Ponnai) |
15 | Vellore | Gudiyatham | Arumbaruthi | 12° 58′ 10″ | 79° 4′ 10″ | Poiney (Ponnai) |
16 | Vellore | Gudiyatham | Panamadangi | 12° 58′ 5″ | 79° 12′ 25″ | Poiney (Ponnai) |
17 | Vellore | Vaniyambadi | Ambur | 13° 1′ 30″ | 79° 2′ 20″ | Upper Palar |
18 | Vellore | Vaniyambadi | Chengilikuppam | 12° 46′ | 78° 40' | Upper Palar |
19 | Vellore | Vaniyambadi | Kadavalam | 12° 42′ 40″ | 78° 39′ 30″ | Upper Palar |
20 | Vellore | Vaniyambadi | Agaramcheri | 12° 46′ 35″ | 78° 39′ 35″ | Agaramar |
21 | Kancheepuram | Uthiramerur | Malayankulam | 12° 53′ 40″ | 78° 53′ 50″ | Cheyyar |
22 | Thiruvannamalai | Vandavasi | Gangapuram | 12° 42′ 30″ | 79° 48′ 40″ | Cheyyar |
23 | Thiruvannamalai | Vandavasi | Sennavaram | 12° 34' | 79° 18' | Kiliyar |
24 | Thiruvannamalai | Vandavasi | Gengapuram | 12° 30′ 7″ | 79° 37′ 20″ | Cheyyar |
25 | Thiruvannamalai | Vandavasi | Nedungunam | 12° 33′ 47″ | 79° 19′ 20″ | Cheyyar |
26 | Thiruvannamalai | Polur | Padagam | 12° 28′ 7″ | 79° 23′ 13″ | Cheyyar |
27 | Thiruvannamalai | Polur | Chetpet | 12° 25′ 30″ | 79° 11′ 45″ | Cheyyar |
28 | Thiruvannamalai | Polur | Vadamathimangalam | 12° 27′ 45″ | 79° 21′ 15″ | Cheyyar |
29 | Thiruvannamalai | Polur | Venmani | 12° 34′ 50″ | 79° 11′ 45″ | Cheyyar |
30 | Thiruvannamalai | Cheyyar | Koolamandal | 12° 30′ 30″ | 79° 8′ 50″ | Cheyyar |
31 | Thiruvannamalai | Cheyyar | Echur | 12° 40′ 58″ | 79° 39′ 30″ | Kiliyar |
32 | Thiruvannamalai | Cheyyar | Koolamandal | 12° 33′ 18″ | 79° 34′ | Cheyyar |
33 | Thiruvannamalai | Cheyyar | Dusi | 12° 40′ 60″ | 79° 39′ 50″ | Cheyyar |
34 | Thiruvannamalai | Cheyyar | Korkai | 12° 46′ 22″ | 79° 41′ 10″ | Cheyyar |
35 | Thiruvannamalai | Cheyyar | Tirumpoondi | 12° 38′ 10″ | 79° 27' | Kiliyar |
36 | Thiruvannamalai | Cheyyar | Tirumani | 12° 34′ 30″ | 79° 33′ 35″ | Cheyyar |
37 | Vellore | Walajahpet | Lalapettai | 12° 36′ 28″ | 79° 21′ 5″ | Poiney (Ponnai) |
38 | Vellore | Vellore | Pennathur | 13° 0′ 20″ | 79° 18′ 10″ | Poiney (Ponnai) |
39 | Vellore | Vellore | Perumugai | 12° 50′ 40″ | 79° 6′ 40″ | Poiney (Ponnai) |
40 | Vellore | Vellore | Usur | 12° 56′ 33″ | 79° 11′ 49″ | Poiney (Ponnai) |
41 | Vellore | Vellore | Vallam | 12° 51′ 54″ | 79° 3′ 43″ | Cheyyar |
42 | Vellore | Vellore | Melvallam | 12° 45′ 15″ | 79° 9′ 18″ | Cheyyar |
43 | Vellore | Vellore | Usur | 12° 45′ 10″ | 79° 9′ 10″ | Poiney (Ponnai) |
44 | Thiruvannamalai | Chengam | Pachal | 12° 51′ 50″ | 79° 3′ 40″ | Cheyyar |
45 | Thiruvannamalai | Chengam | Periyakulam | 12° 17′ 38″ | 78° 57′ 26″ | Cheyyar |
46 | Thiruvannamalai | Cheyyar | Thandarai | 12° 19′ 30″ | 78° 59′ 50″ | Cheyyar |
47 | Thiruvannamalai | Cheyyar | Sumangali | 12° 39' | 79° 30' | Cheyyar |
48 | Thiruvannamalai | Cheyyar | Natteri | 12° 45' | 79° 34' | Vegavati |
49 | Vellore | Walajahpet | Gudimallur | 12° 49′ 10″ | 79° 30′ 5″ | Vegavati |
50 | Vellore | Vellore | Karadikudi | 12° 53′ 35″ | 79° 22′ 15″ | Agaramar |
51 | Vellore | Vellore | Pinnathurai | 12° 50′ 38″ | 78° 55′ 49″ | Agaramar |
52 | Vellore | Vellore | P.S. Mangalam | 12° 46′ 43″ | 78° 56′ 11″ | Poiney (Ponnai) |
53 | Vellore | Vellore | Poigai satyammangalam | 12° 55′ 5″ | 79° 3' | Poiney (Ponnai) |
54 | Vellore | Gudiyatham | Nellorepet | 12° 55′ 5″ | 79° 3' | Malattar |
55 | Vellore | Gudiyatham | Valathur | 12° 56′ 25″ | 78° 41′ 5″ | Kavundinyanadhi |
56 | Vellore | Gudiyatham | Erukkampattu | 12° 52′ 50″ | 78° 49′ 45″ | Malattar |
57 | Kancheepuram | Chengalpattu | Pattaravakkam | 12° 56′ 50″ | 78° 41′ 5″ | Lower Palar |
58 | Thiruvannamalai | Arni | Hariharan Colony | 12° 42′ 56″ | 80° 1′ 10″ | Cheyyar |
59 | Thiruvannamalai | Arni | Somanthangal | 12° 40′ | 79° 18' | Cheyyar |
60 | Thiruvannamalai | Arni | Devikapuram | 12° 42′ 50″ | 79° 12′ 45″ | Cheyyar |
61 | Thiruvannamalai | Arni | Vadugasathu | 12° 29′ 55″ | 79° 14′ 50″ | Cheyyar |
62 | Thiruvannamalai | Arni | Pulavanpadi | 12° 37′ | 79° 16′ 40″ | Cheyyar |
63 | Vellore | Vaniyambadi | Vaniyambadi | 12° 36′ 50″ | 79° 14′ 50″ | Upper Palar |
64 | Vellore | Vaniyambadi | Vellakuttai | 12° 37′ | 78° 38′ | Upper Palar |
65 | Vellore | Vaniyambadi | Thekkupattu | 12° 39' | 78° 41′ 40″ | Upper Palar |
66 | Vellore | Vaniyambadi | Dasiriappanur | 12° 38′ 55″ | 78° 33′ 35″ | Upper Palar |
67 | Kancheepuram | Thirukazhukkuntram | P V Kalathur | 12° 37′ 5″ | 78° 31' | Lower Palar |
68 | Kancheepuram | Thirukazhukkuntram | Panangattucherry | 12° 37′ 5″ | 79° 58′ 31″ | Lower Palar |
69 | Vellore | Arcot | Kalavai | 12° 30′ 10″ | 80° 4′ 55″ | Cheyyar |
70 | Vellore | Arcot | Pudupadi | 12° 46′ | 79° 23′ 5″ | Cheyyar |
71 | Vellore | Thiruppathur | Andiyappanur | 12° 51′ 55″ | 79° 22′ 25″ | Cheyyar |
72 | Kancheepuram | Sriperumbudur | Maduramangalam | 12° 22′ 10″ | 78° 42′ 15″ | Lower Palar |
73 | Kancheepuram | Uthiramerur | Walajabad | 12° 57′ | 79° 49' | Lower Palar |
74 | Thiruvannamalai | Cheyyar | Valavanur | 12° 47′ | 79° 49′ 50″ | Cheyyar |
75 | Vellore | Arcot | Melapalandai | 12° 48' | 79° 30' | Cheyyar |
76 | Vellore | Arcot | Varagur | 12° 39′ 30″ | 79° 23′ 42″ | Cheyyar |
77 | Vellore | Arcot | Varagurpattanam | 12° 48′ 40″ | 79° 14′ 15″ | Cheyyar |
78 | Kancheepuram | Sriperumbudur | Oragadam | 12° 48′ 50″ | 79° 14′ 50″ | Lower Palar |
79 | Kancheepuram | Uthiramerur | Thenneri | 12° 50′ 40″ | 79° 56′ 50″ | Lower Palar |
80 | Kancheepuram | Kancheepuram | Damal | 12° 52' | 79° 51′ 40″ | Vegavati |
81 | Kancheepuram | Kancheepuram | Kilambi | 12° 53′ | 79° 35′ 45″ | Vegavati |
82 | Kancheepuram | Kancheepuram | Magaral | 12° 51′ 30″ | 79° 39′ 20″ | Cheyyar |
83 | Vellore | Thiruppathur | Kothur | 12° 43′ | 79° 45′ 20″ | Upper Palar |
Well no . | District . | Taluk . | Village . | Latitude . | Longitude . | Subbasin . |
---|---|---|---|---|---|---|
1 | Kancheepuram | Madurandagam | Thennampattu | 12° 28′ | 79° 55′ 40″ | Kiliyar |
2 | Kancheepuram | Madurandagam | Velamur [ramapuram] | 12°28′40″ | 79° 46' | Kiliyar |
3 | Thiruvannamalai | Vandavasi | Tennangur | 12° 24′ 5″ | 79° 53′ 40″ | Kiliyar |
4 | Thiruvannamalai | Vandavasi | Salavedu | 12° 34′ | 79° 38′ | Kiliyar |
5 | Thiruvannamalai | Vandavasi | Peranamallur | 12° 28′ | 79° 46′ 32″ | Cheyyar |
6 | Thiruvannamalai | Vandavasi | Osur | 12° 34′ 17″ | 79° 25′ 59″ | Kiliyar |
7 | Thiruvannamalai | Thiruvannamalai | Nayadimangalam | 12° 26′ 15″ | 79° 40′ 10″ | Cheyyar |
8 | Thiruvannamalai | Polur | Kadaladi | 12° 24′ 18″ | 79° 6′ 6″ | Cheyyar |
9 | Thiruvannamalai | Polur | Mampattu | 12° 24′ 15″ | 78° 58′ 10″ | Cheyyar |
10 | Thiruvannamalai | Polur | Siruvallur | 12° 30′ 35″ | 79° 5′ 30″ | Cheyyar |
11 | Thiruvannamalai | Polur | Athuvambadi | 12° 27′ 45″ | 79° 1′ 20″ | Cheyyar |
12 | Thiruvannamalai | Polur | Melarani | 12° 35′ 50″ | 79° 7′ 35″ | Cheyyar |
13 | Thiruvannamalai | Chengam | Chengam Farm 2 | 12° 27′ 20″ | 79° 2′ 45″ | Cheyyar |
14 | Vellore | Katpadi | Vaduganthangal | 12° 16′ 20″ | 78° 43′ 10″ | Poiney (Ponnai) |
15 | Vellore | Gudiyatham | Arumbaruthi | 12° 58′ 10″ | 79° 4′ 10″ | Poiney (Ponnai) |
16 | Vellore | Gudiyatham | Panamadangi | 12° 58′ 5″ | 79° 12′ 25″ | Poiney (Ponnai) |
17 | Vellore | Vaniyambadi | Ambur | 13° 1′ 30″ | 79° 2′ 20″ | Upper Palar |
18 | Vellore | Vaniyambadi | Chengilikuppam | 12° 46′ | 78° 40' | Upper Palar |
19 | Vellore | Vaniyambadi | Kadavalam | 12° 42′ 40″ | 78° 39′ 30″ | Upper Palar |
20 | Vellore | Vaniyambadi | Agaramcheri | 12° 46′ 35″ | 78° 39′ 35″ | Agaramar |
21 | Kancheepuram | Uthiramerur | Malayankulam | 12° 53′ 40″ | 78° 53′ 50″ | Cheyyar |
22 | Thiruvannamalai | Vandavasi | Gangapuram | 12° 42′ 30″ | 79° 48′ 40″ | Cheyyar |
23 | Thiruvannamalai | Vandavasi | Sennavaram | 12° 34' | 79° 18' | Kiliyar |
24 | Thiruvannamalai | Vandavasi | Gengapuram | 12° 30′ 7″ | 79° 37′ 20″ | Cheyyar |
25 | Thiruvannamalai | Vandavasi | Nedungunam | 12° 33′ 47″ | 79° 19′ 20″ | Cheyyar |
26 | Thiruvannamalai | Polur | Padagam | 12° 28′ 7″ | 79° 23′ 13″ | Cheyyar |
27 | Thiruvannamalai | Polur | Chetpet | 12° 25′ 30″ | 79° 11′ 45″ | Cheyyar |
28 | Thiruvannamalai | Polur | Vadamathimangalam | 12° 27′ 45″ | 79° 21′ 15″ | Cheyyar |
29 | Thiruvannamalai | Polur | Venmani | 12° 34′ 50″ | 79° 11′ 45″ | Cheyyar |
30 | Thiruvannamalai | Cheyyar | Koolamandal | 12° 30′ 30″ | 79° 8′ 50″ | Cheyyar |
31 | Thiruvannamalai | Cheyyar | Echur | 12° 40′ 58″ | 79° 39′ 30″ | Kiliyar |
32 | Thiruvannamalai | Cheyyar | Koolamandal | 12° 33′ 18″ | 79° 34′ | Cheyyar |
33 | Thiruvannamalai | Cheyyar | Dusi | 12° 40′ 60″ | 79° 39′ 50″ | Cheyyar |
34 | Thiruvannamalai | Cheyyar | Korkai | 12° 46′ 22″ | 79° 41′ 10″ | Cheyyar |
35 | Thiruvannamalai | Cheyyar | Tirumpoondi | 12° 38′ 10″ | 79° 27' | Kiliyar |
36 | Thiruvannamalai | Cheyyar | Tirumani | 12° 34′ 30″ | 79° 33′ 35″ | Cheyyar |
37 | Vellore | Walajahpet | Lalapettai | 12° 36′ 28″ | 79° 21′ 5″ | Poiney (Ponnai) |
38 | Vellore | Vellore | Pennathur | 13° 0′ 20″ | 79° 18′ 10″ | Poiney (Ponnai) |
39 | Vellore | Vellore | Perumugai | 12° 50′ 40″ | 79° 6′ 40″ | Poiney (Ponnai) |
40 | Vellore | Vellore | Usur | 12° 56′ 33″ | 79° 11′ 49″ | Poiney (Ponnai) |
41 | Vellore | Vellore | Vallam | 12° 51′ 54″ | 79° 3′ 43″ | Cheyyar |
42 | Vellore | Vellore | Melvallam | 12° 45′ 15″ | 79° 9′ 18″ | Cheyyar |
43 | Vellore | Vellore | Usur | 12° 45′ 10″ | 79° 9′ 10″ | Poiney (Ponnai) |
44 | Thiruvannamalai | Chengam | Pachal | 12° 51′ 50″ | 79° 3′ 40″ | Cheyyar |
45 | Thiruvannamalai | Chengam | Periyakulam | 12° 17′ 38″ | 78° 57′ 26″ | Cheyyar |
46 | Thiruvannamalai | Cheyyar | Thandarai | 12° 19′ 30″ | 78° 59′ 50″ | Cheyyar |
47 | Thiruvannamalai | Cheyyar | Sumangali | 12° 39' | 79° 30' | Cheyyar |
48 | Thiruvannamalai | Cheyyar | Natteri | 12° 45' | 79° 34' | Vegavati |
49 | Vellore | Walajahpet | Gudimallur | 12° 49′ 10″ | 79° 30′ 5″ | Vegavati |
50 | Vellore | Vellore | Karadikudi | 12° 53′ 35″ | 79° 22′ 15″ | Agaramar |
51 | Vellore | Vellore | Pinnathurai | 12° 50′ 38″ | 78° 55′ 49″ | Agaramar |
52 | Vellore | Vellore | P.S. Mangalam | 12° 46′ 43″ | 78° 56′ 11″ | Poiney (Ponnai) |
53 | Vellore | Vellore | Poigai satyammangalam | 12° 55′ 5″ | 79° 3' | Poiney (Ponnai) |
54 | Vellore | Gudiyatham | Nellorepet | 12° 55′ 5″ | 79° 3' | Malattar |
55 | Vellore | Gudiyatham | Valathur | 12° 56′ 25″ | 78° 41′ 5″ | Kavundinyanadhi |
56 | Vellore | Gudiyatham | Erukkampattu | 12° 52′ 50″ | 78° 49′ 45″ | Malattar |
57 | Kancheepuram | Chengalpattu | Pattaravakkam | 12° 56′ 50″ | 78° 41′ 5″ | Lower Palar |
58 | Thiruvannamalai | Arni | Hariharan Colony | 12° 42′ 56″ | 80° 1′ 10″ | Cheyyar |
59 | Thiruvannamalai | Arni | Somanthangal | 12° 40′ | 79° 18' | Cheyyar |
60 | Thiruvannamalai | Arni | Devikapuram | 12° 42′ 50″ | 79° 12′ 45″ | Cheyyar |
61 | Thiruvannamalai | Arni | Vadugasathu | 12° 29′ 55″ | 79° 14′ 50″ | Cheyyar |
62 | Thiruvannamalai | Arni | Pulavanpadi | 12° 37′ | 79° 16′ 40″ | Cheyyar |
63 | Vellore | Vaniyambadi | Vaniyambadi | 12° 36′ 50″ | 79° 14′ 50″ | Upper Palar |
64 | Vellore | Vaniyambadi | Vellakuttai | 12° 37′ | 78° 38′ | Upper Palar |
65 | Vellore | Vaniyambadi | Thekkupattu | 12° 39' | 78° 41′ 40″ | Upper Palar |
66 | Vellore | Vaniyambadi | Dasiriappanur | 12° 38′ 55″ | 78° 33′ 35″ | Upper Palar |
67 | Kancheepuram | Thirukazhukkuntram | P V Kalathur | 12° 37′ 5″ | 78° 31' | Lower Palar |
68 | Kancheepuram | Thirukazhukkuntram | Panangattucherry | 12° 37′ 5″ | 79° 58′ 31″ | Lower Palar |
69 | Vellore | Arcot | Kalavai | 12° 30′ 10″ | 80° 4′ 55″ | Cheyyar |
70 | Vellore | Arcot | Pudupadi | 12° 46′ | 79° 23′ 5″ | Cheyyar |
71 | Vellore | Thiruppathur | Andiyappanur | 12° 51′ 55″ | 79° 22′ 25″ | Cheyyar |
72 | Kancheepuram | Sriperumbudur | Maduramangalam | 12° 22′ 10″ | 78° 42′ 15″ | Lower Palar |
73 | Kancheepuram | Uthiramerur | Walajabad | 12° 57′ | 79° 49' | Lower Palar |
74 | Thiruvannamalai | Cheyyar | Valavanur | 12° 47′ | 79° 49′ 50″ | Cheyyar |
75 | Vellore | Arcot | Melapalandai | 12° 48' | 79° 30' | Cheyyar |
76 | Vellore | Arcot | Varagur | 12° 39′ 30″ | 79° 23′ 42″ | Cheyyar |
77 | Vellore | Arcot | Varagurpattanam | 12° 48′ 40″ | 79° 14′ 15″ | Cheyyar |
78 | Kancheepuram | Sriperumbudur | Oragadam | 12° 48′ 50″ | 79° 14′ 50″ | Lower Palar |
79 | Kancheepuram | Uthiramerur | Thenneri | 12° 50′ 40″ | 79° 56′ 50″ | Lower Palar |
80 | Kancheepuram | Kancheepuram | Damal | 12° 52' | 79° 51′ 40″ | Vegavati |
81 | Kancheepuram | Kancheepuram | Kilambi | 12° 53′ | 79° 35′ 45″ | Vegavati |
82 | Kancheepuram | Kancheepuram | Magaral | 12° 51′ 30″ | 79° 39′ 20″ | Cheyyar |
83 | Vellore | Thiruppathur | Kothur | 12° 43′ | 79° 45′ 20″ | Upper Palar |
MATERIALS AND METHODS
Data
The investigation of the quality of groundwater of the Palar Basin was undertaken for the years 1993, 1983, 1973, 2003, 2013, and 2022. The data collected from Central Ground Water Board (CGWB), Chennai includes the information related to the water quality parameter (WQP) values for pre- and post-monsoon. The study took into account various physicochemical parameters, including pH, electrical conductivity (EC), and the concentrations of key ions such as calcium (Ca), magnesium (Mg), sodium (Na), potassium (K), chloride (Cl), and SO4. These parameters were utilized to compute the irrigation indices, which served as critical indicators for assessing water suitability for agricultural purposes.
Water quality parameters
The study examined a range of physicochemical parameters, including pH, EC, and concentrations of essential ions such as Ca, Mg, Na, K, Cl, and SO4. These parameters were used for assessment of drinking water quality and also used to calculate irrigation indices, which are vital indicators for evaluating the suitability of water for agricultural use. Table 2 presents the maximum permissible limits for water quality parameters and their impact on irrigation water.
Bureau of Indian Standards (BIS) maximum permissible limit (1998) of water quality parameters (Ravikumar et al. 2011)
Parameter . | Unit . | Max. permissible limit . | Effect on irrigation water . |
---|---|---|---|
pH | – | 6.5–8.5 | Optimal for most crops. pH outside this range can affect nutrient availability and soil structure |
Ca | mg/L | 200 | High levels improve soil structure but can lead to scaling |
Mg | mg/L | 100 | Excess Mg can contribute to salinity and affect soil permeability |
Na | mg/L | 200 | High Na causes soil dispersion, poor permeability, and reduced crop yield |
K | mg/L | 12 | Rarely problematic; an essential nutrient for plant growth. Levels above 10 mg/L can disrupt nutrient balance |
Cl | mg/L | 1,000 | High Cl can be toxic to sensitive crops, causing leaf burn |
SO4 | mg/L | 400 | High SO4 adds to overall salinity and affects sensitive crops |
Parameter . | Unit . | Max. permissible limit . | Effect on irrigation water . |
---|---|---|---|
pH | – | 6.5–8.5 | Optimal for most crops. pH outside this range can affect nutrient availability and soil structure |
Ca | mg/L | 200 | High levels improve soil structure but can lead to scaling |
Mg | mg/L | 100 | Excess Mg can contribute to salinity and affect soil permeability |
Na | mg/L | 200 | High Na causes soil dispersion, poor permeability, and reduced crop yield |
K | mg/L | 12 | Rarely problematic; an essential nutrient for plant growth. Levels above 10 mg/L can disrupt nutrient balance |
Cl | mg/L | 1,000 | High Cl can be toxic to sensitive crops, causing leaf burn |
SO4 | mg/L | 400 | High SO4 adds to overall salinity and affects sensitive crops |
Water quality using irrigation indices for irrigation water
The quality of water for irrigation has been assessed using various irrigation indices. It includes SAR, RSC, Na%, PI, magnesium hazard (MH), KR, and potential salinity (PS). The equations for calculating irrigation indices are given in Table 3, where all ion concentrations used in the equations are in meq/L. The groundwater categorization on irrigation indices is visible in Table 4. SAR, RSC, Na%, PI, MH, KR, PS were chosen for their ability to provide a comprehensive assessment of irrigation water quality and its impact on soil and plant health. These indices address critical factors such as sodium toxicity, soil permeability, salinity, and the balance of essential ions. SAR and Na% evaluate the risk of sodium-induced soil degradation, while RSC and PI focus on the alkalinity and permeability effects on soil. MH helps assess Mg toxicity risks, and KR indicates the ionic balance in water that influences plant growth. PS evaluates the salinity potential, which can affect plant roots and crop yield. Together, these indices offer a holistic understanding of water quality, ensuring a thorough evaluation of its suitability for irrigation in varying hydrogeological conditions.
Equations for the calculated irrigation indices
Irrigation indices . | Equation number . | References . |
---|---|---|
![]() | (1) | Raghunath (1987) |
![]() | (2) | Nageswara Rao et al. (2022) |
![]() | (3) | Todd (1980) |
![]() | (4) | Doneen (1964) |
![]() | (5) | Al-Shammiri et al. (2005), Joshi et al. (2009). |
![]() | (6) | Kelly (1940) |
![]() | (7) | Doneen (1964) |
Irrigation indices . | Equation number . | References . |
---|---|---|
![]() | (1) | Raghunath (1987) |
![]() | (2) | Nageswara Rao et al. (2022) |
![]() | (3) | Todd (1980) |
![]() | (4) | Doneen (1964) |
![]() | (5) | Al-Shammiri et al. (2005), Joshi et al. (2009). |
![]() | (6) | Kelly (1940) |
![]() | (7) | Doneen (1964) |
Groundwater categorization based on irrigation indices (Kelly 1940; Doneen 1964; Szabolcs 1964; Todd 1980; Wilcox 1955; Nageswara Rao et al. 2022)
Parameter/Indices . | Values . | Category . |
---|---|---|
SAR | <10 | Excellent |
10–18 | Good | |
18–26 | Doubtful | |
>26 | Unsuitable | |
RSC | <1.25 | Excellent |
1.25–2.5 | Good | |
>2.5 | Unsuitable | |
Na% | <20 | Excellent |
20–40 | Good | |
40–60 | Permissible | |
60–80 | Doubtful | |
>80 | Unsuitable | |
PI | >75 | Good |
25–75 | Slightly good | |
<25 | Unsuitable | |
MH | <50 | Suitable |
>50 | Unsuitable | |
KR | <1 | Suitable |
>1 | Unsuitable | |
PS | <3 | Safe |
>3 | Unsuitable |
Parameter/Indices . | Values . | Category . |
---|---|---|
SAR | <10 | Excellent |
10–18 | Good | |
18–26 | Doubtful | |
>26 | Unsuitable | |
RSC | <1.25 | Excellent |
1.25–2.5 | Good | |
>2.5 | Unsuitable | |
Na% | <20 | Excellent |
20–40 | Good | |
40–60 | Permissible | |
60–80 | Doubtful | |
>80 | Unsuitable | |
PI | >75 | Good |
25–75 | Slightly good | |
<25 | Unsuitable | |
MH | <50 | Suitable |
>50 | Unsuitable | |
KR | <1 | Suitable |
>1 | Unsuitable | |
PS | <3 | Safe |
>3 | Unsuitable |
The WQI is essential for determining overall water quality, while irrigation indices like Na%, SAR, PI, and MH are critical for evaluating groundwater's suitability for irrigation. These indices help assess soil salinity and sodicity, which impact crop productivity and soil health, making a combined analysis with WQI vital for sustainable agricultural practices (Shukla et al. 2023).
Water quality index
A composite index, known as the WQI, was employed to evaluate the suitability of groundwater for various purposes. The primary objective is to consolidate extensive and intricate water quality data into a single, comprehensive index (Ramakrishnaiah et al. 2009; Varol & Davraz 2015). The WQI is a valuable tool developed to simplify the assessment of overall water quality by combining various parameters into a single value, facilitating effective management and categorization of water sources for their suitability in different uses. Initially developed by Horton in 1965, the WQI has evolved over time, incorporating complex datasets to provide insights for governmental, public, and regulatory bodies, despite the absence of a universally recognized index method (Rana & Ganguly 2020). The standard values and unit weight values of irrigation indices are given in Table 5.
Standard values and unit weight values used for the calculation of WQI
Irrigation index . | Standard values (![]() | Unit weight values (![]() |
---|---|---|
SAR | 10 | 0.052356 |
RSC | 2.5 | 0.209424 |
Na% | 60 | 0.008726 |
PI | 25 | 0.020942 |
MH | 50 | 0.010471 |
KR | 1 | 0.52356 |
PS | 3 | 0.17452 |
Total weightage | 1 |
Irrigation index . | Standard values (![]() | Unit weight values (![]() |
---|---|---|
SAR | 10 | 0.052356 |
RSC | 2.5 | 0.209424 |
Na% | 60 | 0.008726 |
PI | 25 | 0.020942 |
MH | 50 | 0.010471 |
KR | 1 | 0.52356 |
PS | 3 | 0.17452 |
Total weightage | 1 |
Procedure for the determination of WQI (weighted arithmetic index method)
The WQI has been calculated in three steps using weighted arithmetic index methods (Ramakrishnaiah et al. 2009). The categorization of WQI in accordance with restriction in using groundwater for irrigation purpose is given in Table 6.
Categorization of the WQI for irrigation purposes (Kumar & Maurya 2023)
WQI value . | Restriction for using GW . |
---|---|
<150 | None |
150–300 | Slight |
300–450 | Moderate |
>450 | Severe |
WQI value . | Restriction for using GW . |
---|---|
<150 | None |
150–300 | Slight |
300–450 | Moderate |
>450 | Severe |
Quality rating or sub-index (Qn)
where
is the quality rating value for the nth water quality parameter;
is the actual value of the nth parameter;
is the ideal value of the parameter;
is the standard value of the nth parameter. Consider vi = 0 for all the samples except pH. For pH, the ideal value is set as vi = 7, corresponding to the neutral condition.
Unit weight value
Water quality index
Spatial modelling and interpolation through IDW
The regional distribution of groundwater quality metrics was interpolated using spatial modeling, an advanced capability of ArcGIS 10.4. The spatial resolution of the groundwater quality data used in this GIS analysis is approximately 5.7 km × 5.7 km, suggesting that the points are not evenly distributed across the study area. Finer resolution in well-dense areas improves mapping accuracy, while coarser resolution in sparse zones balances efficiency and data availability. Regions with closely spaced wells yield more accurate representations of water quality, whereas sparsely monitored areas have lower accuracy.
The IDW interpolation method was specifically utilized to create detailed spatial distribution maps of water quality parameters. This technique effectively interpolates data to assess observations between measured points. Numerous researchers have adopted the IDW method for interpolation in their studies, demonstrating its widespread application and reliability (Mueller et al. 2004). In assessing WQI, the IDW method within GIS modeling offers valuable spatial insights by interpolating parameters like temperature, dissolved oxygen (DO), nitrates (NO3), and total phosphorus across various depths. This integration with traditional methods allows for comprehensive spatial analysis of water quality, enhancing monitoring and understanding of pollution dynamics within water bodies (Vasistha & Ganguly 2022). In this study, IDW was selected for its effectiveness in estimating spatial variations in groundwater quality from scattered sample points. IDW assumes that closer points exert more influence, fitting well with our data distribution. Unlike Kriging, which requires a complex variogram model, IDW is computationally efficient and adaptable to irregular sampling patterns. Validated widely in groundwater studies, IDW aligns with established practices and provides reproducible, detailed maps. This method avoids the smoothing effect seen in Kriging, preserving small-scale quality variations and ensuring that local influences are accurately represented in our groundwater quality maps.
RESULTS AND DISCUSSIONS
pH
Electrical conductivity
Salt accumulation in groundwater poses a major environmental issue across several regions of Tamil Nadu, particularly those with substantial industrial activities. Industrial effluents often contain high concentrations of dissolved solids, heavy metals, and various contaminants that infiltrate groundwater, leading to increased EC and degraded water quality (Kanagaraj & Elango 2016). Over time, this contamination not only raises groundwater salinity but also threatens public health and agricultural sustainability. Extended use of high-EC water for irrigation can lead to soil salinization, adversely affecting crop productivity and soil structure (Richards 1954).
Sodium
Calcium
Magnesium
Potassium
Chlorine
Sulphate
Sodium absorption ratio
Residual sodium carbonate
Sodium percentage (Na %)
Permeability index
Magnesium hazard
Kelly's ratio
Potential salinity
Water quality index
The spatial distribution of WQI showed that the majority of water is coming under categories that have no restriction in pre- and post-monsoon periods except in 1993. An increase in water quality parameters like Ca, Mg, Na, and other parameters may contribute to the change in water quality. Contaminants dissolved in the runoff during monsoon reach the groundwater and at last, it joins with seawater resulting in the deterioration of water quality. Comparing the spatial distribution of the WQI for selected years (1973, 1983, 1993, 2003, and 2013) showed no drastic changes in the WQI up to 1993. In 1993, the greater part of the study area fell into the ‘slight restriction’ category. However, after 1993, it shifted to the ‘no restriction’ category. The western part of the region fell under the 'severe restriction' category during the pre-monsoon season of 1993.
Increased agricultural runoff, industrial discharges, or urbanization can introduce pollutants into groundwater, degrading its quality. Additionally, natural processes like mineral leaching and climate change impacts can alter groundwater composition, necessitating slight usage restrictions to protect human health and agricultural productivity. We have to take measures to control the decrease in the quality of water by strict enforcement of environmental regulations on industrial effluent discharge, improving the wastewater treatment facilities, and encouraging the construction of recharge structures in the desired areas.
From Table 7, it is evident that irrigation indices, such as RSC, MH, and PS, do not show significant seasonal variation (p > 0.05). Significant seasonal variations are observed in key irrigation indices such as SAR, Na%, PI, KR, and WQI (p < 0.05). These fluctuations indicate changes in water quality and soil interactions between the pre- and post-monsoon seasons. Specifically, the variations in sodium adsorption ratio, sodium percentage, and the balance of sodium with other cations may impact soil permeability and fertility, necessitating targeted adjustments in irrigation practices. The observed seasonal shifts in water quality emphasize the need for intensive management strategies during periods of significant change to ensure optimal soil and water health.
Significance levels of t-test for irrigation indices
Indices . | t-statistics . | p-value . | No significant difference . | Significant difference . |
---|---|---|---|---|
SAR | −0.30899 | 0.0026 | ✓ | |
RSC | 0.1826 | 0.8555 | ✓ | |
Na% | −2.4950 | 0.0143 | ✓ | |
PI | −2.2416 | 0.0272 | ✓ | |
MH | −0.3085 | 0.7584 | ✓ | |
KR | −3.0084 | 0.0033 | ✓ | |
PS | 0.2989 | 0.7657 | ✓ | |
WQI | −3.7821 | 0.003 | ✓ |
Indices . | t-statistics . | p-value . | No significant difference . | Significant difference . |
---|---|---|---|---|
SAR | −0.30899 | 0.0026 | ✓ | |
RSC | 0.1826 | 0.8555 | ✓ | |
Na% | −2.4950 | 0.0143 | ✓ | |
PI | −2.2416 | 0.0272 | ✓ | |
MH | −0.3085 | 0.7584 | ✓ | |
KR | −3.0084 | 0.0033 | ✓ | |
PS | 0.2989 | 0.7657 | ✓ | |
WQI | −3.7821 | 0.003 | ✓ |
LIMITATIONS OF STUDY
This research was designed to measure the spatial characterization of groundwater quality of the Palar basin for pre- and post-monsoon of the selected years. The study, relying on data from the CGWB, may be limited by its temporal resolution and spatial coverage, potentially affecting the granularity of the analysis. Although seasonal variations are considered, the complexity of seasonal factors influencing water quality, such as unrecorded weather events or extreme climate occurrences, may not be entirely captured. Additionally, the use of IDW interpolation for spatial analysis assumes uniformity in data distribution, which could introduce bias or inaccuracies in areas with sparse data points. Furthermore, the study may not fully address the impact of non-point source pollution, such as runoff from urban areas or diffuse agricultural runoff, which can significantly influence water quality in certain regions. On the contrary, it is important to recognize that assessing the quality of water in the context of climate change is the most effective way to understand the hydrological behavior of groundwater. Therefore, this research offers valuable insights into the spatial characteristics of water quality parameters of water resources, emphasizing that failure to maintain groundwater quality could have severe consequences for communities. As such, this study is highly relevant for water resources management, highlighting the need for effective measures to safeguard groundwater quality.
SUMMARY AND CONCLUSION
Groundwater quality and its relevance for irrigation purposes in the Palar Basin area have been evaluated. The water quality and pollution status of the Palar basin are critically important as they directly impact human health. Nearly 90% of diseases are caused by consuming contaminated water, and rivers serve as the primary water source. The Tamil Nadu government relies heavily on groundwater, using 80% of it for water supply. A change in water quality parameters over time was observed and irrigation indices were calculated.
Groundwater quality parameter . | Key findings . |
---|---|
pH | No significant variations; generally, within acceptable limits for irrigation |
EC | Ranged from 250 μS/cm to over 3,000 μS/cm, with elevated levels in industrial areas like Ranipet |
SAR | Predominantly ‘excellent’ (0–10), but exceeded 18 in specific southwest regions during 1993 pre-monsoon, indicating irrigation unsuitability |
Na% | Varied between 20 and 90%, with values >60% causing degradation of soil permeability |
PI | Shows spatial and seasonal variability influencing soil permeability and irrigation suitability |
RSC | Mostly within safe limits (<1.25 meq/L); localized areas showed values >2.5 meq/L, reducing soil productivity |
MH | Exceeded the 50% threshold in certain regions, especially near industrial and agricultural zones, making groundwater unsuitable for irrigation |
Ca | Concentrations up to 180 mg/L contributing to increased water hardness, affecting agricultural, industrial, and domestic use |
Mg | Concentrations up to 120 mg/L, contributing to water hardness |
Cl | Generally, below the permissible limit except in 1993, where it exceeded 2,000 mg/L |
SO4 | Increased during the post-monsoon period due to chemical cleaning agents, organic decomposition, industrial discharge, and mining activities |
WQI | Most samples were classified as ‘none’ to ‘slight’ restriction for use, with some hotspots having WQI >450 indicating severe restriction for irrigation |
Overall findings | Seasonal variations were observed, with post-monsoon surface runoff and dilution effects contributing to elevated sodium and SO4 levels. Localized Cl and SO4 increases pose minimal risks |
Groundwater quality parameter . | Key findings . |
---|---|
pH | No significant variations; generally, within acceptable limits for irrigation |
EC | Ranged from 250 μS/cm to over 3,000 μS/cm, with elevated levels in industrial areas like Ranipet |
SAR | Predominantly ‘excellent’ (0–10), but exceeded 18 in specific southwest regions during 1993 pre-monsoon, indicating irrigation unsuitability |
Na% | Varied between 20 and 90%, with values >60% causing degradation of soil permeability |
PI | Shows spatial and seasonal variability influencing soil permeability and irrigation suitability |
RSC | Mostly within safe limits (<1.25 meq/L); localized areas showed values >2.5 meq/L, reducing soil productivity |
MH | Exceeded the 50% threshold in certain regions, especially near industrial and agricultural zones, making groundwater unsuitable for irrigation |
Ca | Concentrations up to 180 mg/L contributing to increased water hardness, affecting agricultural, industrial, and domestic use |
Mg | Concentrations up to 120 mg/L, contributing to water hardness |
Cl | Generally, below the permissible limit except in 1993, where it exceeded 2,000 mg/L |
SO4 | Increased during the post-monsoon period due to chemical cleaning agents, organic decomposition, industrial discharge, and mining activities |
WQI | Most samples were classified as ‘none’ to ‘slight’ restriction for use, with some hotspots having WQI >450 indicating severe restriction for irrigation |
Overall findings | Seasonal variations were observed, with post-monsoon surface runoff and dilution effects contributing to elevated sodium and SO4 levels. Localized Cl and SO4 increases pose minimal risks |
In conclusion, the outcomes of the present study suggest that frequent and comprehensive groundwater analysis is essential to monitor the rate and types of pollution. Raising public awareness of the importance of preserving groundwater's highest standards of quality and purity is imperative.
AUTHORS’ CONTRIBUTIONS
This work was carried out in collaboration with all authors. All authors read and approved the final manuscript.
FUNDING
No funding was obtained for this study.
DATA AVAILABILITY STATEMENT
All relevant data are included in the paper or its Supplementary Information.
CONFLICT OF INTEREST
The authors declare there is no conflict.