Шаблоны LeoTheme для Joomla.
GavickPro Joomla шаблоны

banner hydrology

Original Research Article

Spatial Variation in Arsenic Concentration in the Groundwater of Nawalparasi District, Terai Nepal

Ishwar Chandra Yadav1*, Ningombam Linthoingambi Devi2

1State Key Laboratory of Organic Geochemistry, Guangzhou Institute of Geochemistry, Chinese Academy of Sciences, Guangzhou 510640, China
2Centre for Environmental Sciences, Central University of South Bihar, BIT Campus, Patna-800014, Bihar

*Corresponding author: Dr. Ishwar Chandra Yadav, State Key Laboratory of Organic Geochemistry, Guangzhou Institute of Geochemistry, Chinese Academy of Sciences, Guangzhou 510640, China, Tel: +86-15626134294;
Email: icyadav.bhu@gmail.com

Submitted: 06-23-2015 Accepted: 07-14-2015 Published: 07-15-2015

Download PDF






The spatial variability of arsenic (As) in the groundwater was studied in the three villages of Nawalparasi district, Nepal. The average concentration of As varied as low (0.16mgL-1) in Kasiya to high (1.34mgL-1) in Thulokunwar with all the tubewells exceeding WHO (0.01mgL-1) and NIS (0.05mgL-1) limits for drinking waters. The concentration of As was varied even in closely located tubewells suggesting the possibility of spatial variation due to some conditions such as latitude, longitude and depth of tubewell. The As level also varied inversely with the depth of tubewells.

Keywords: Kasiya; Panchgawa; Thulokunwar; Tubewells; Depth of Water Table; Lithology


The water quality is the major issue not only in Nepal but also in countries across the globe. Generally, groundwater is more susceptible to arsenic (As) contamination than the surfacewater, because groundwater interacts with a numbers of different kinds of minerals in aquifer. It is know that most of the As-related health problems reported so far is associated with consumption of As-contaminated groundwater over period of time. Nowadays, As-contamination of groundwater is of increasing concern due to its high toxicity and the widespread occurrence [1]. The alarming concentrations of As have been reported in South East Asia viz. Bangladesh, Vietnam, West Bengal–India, Nepal, Cambodia, Mongolia, China, Thailand, Pakistan and Taiwan [2-7], Central and South America (Argentina, Chile, Mexico), North America (USA and Canada) [8-12] and Australia [13]. Elevated As concentrations are also reported in European countries i.e. Hungary, Romania, Italy, Greece, Spain, Finland, and Germany [1,14,15].

Groundwater As is highly variable on the local as well as regional scale. The groundwater chemistry also varies both at the spatial as well as temporal scale [16,17]. Although As concentrations vary over distances ranging from meters to kilometers, the cause of the lateral heterogeneity yet remains to be explored. Regional dependency of As contamination accounts for the geologic divisions of the area but at the smaller scales, this pattern is more random. van Geen et al [18] have reported that in the smaller village with higher resolution (6000 tube wells in a 25 km2 area), tube wells having shallower depths are highly variable in As load, and do not follow the clear spatial pattern. It seems that the knowledge of spatial variability of As remains crucial not only for evaluating the biogeochemical-hydrological events that control groundwater concentrations of As but also for developing strategies for sustainable use of shallow groundwater that will reduce exposure to As [19].

In recent years, groundwater As-contamination in Terai region, Nepal has been the major issue [20]. Temporal variations in groundwater As-contamination of the Terai region necessitate mitigation efforts as there has been very little investigation on the spatial and temporal pattern of groundwater chemistry especially in the high-As aquifers. In the present study an attempt has been made to evaluate the spatial distribution of As-contamination in the groundwater of Nawalparasi district of southern Nepal.

Geological setting of the area

Tectonically, Nepal Himalayas can be classified into four different zones (i.e. Tethys Himalayas, Higher Himalayas, Lesser Himalayas, sub Himalayas and Terai plain) from north to south followed by Terai Plain bordering them on the south [21,22]. The Tethysian zone which lies in the northern part of western central Nepal consists of high mountain ranges crowned by cretaceous fossiliferous sediments. The higher Himalayan zone consists of crystalline rocks (schist and gneiss). Main control Thrust (MCT) marks the contact of crystalline zone to that of lesser Himalayan sedimentary pile. The contact of the crystalline with the overlying Tethys sediments is gradational. To the south of higher Himalayas, the lesser Himalaya zone which consists of meta-sediments with different grades of metamorphism. The sub-Himalayan zone, which lies to south of the lesser Himalayas and separated by MBT (Main Boundary Thrust) from it consist of sedimentary rocks containing silt, sand, shale and pebble and boulder beds. The pebbles and boulders are quartzites, sandstone, phyllites etc. These formations form the low lying hill ranges that rise from the Indo-Gangetic plain on the south. These hills are called Churia hills which is the Nepalese designation of Siwalik [22].

The Terai plain lies in the southern part of Nepal and is separated from the Churia formations by the Himalayan Frontal thrust. Nawalparasi lies in the Terai plain as the continuation of Indo-Gangetic plain and has gentle slope toward south from an elevation of 200-300m in the north to as low as 63m in the south near Indian border [21, 22]. The lithology of the Terai sedimentary basin belongs to Holocene alluvium that includes the present day alluvial deposits, channel sand and gravel deposits and outwash deposits [23].

Materials and Methods

Study area

Nawalparasi, a part of Lumbini zone (one of the 75 districts) is located in the western Terai region of southern Nepal. It lies 147 Km. west of the capital city Kathmandu, and surrounded by Chitwan district in the east, Tanahun in the north, Palpa and Rupandehi districts in the west, and the Indian state of Uttar Pradesh in the south. The present study was conducted in three villages namely Thulokunwar, Kasiya and Panchgawa of Nawalparasi district of Nepal (Figure 1) where the As-contamination of groundwater has variously been reported [20,24]. A total of 24 tubewells were selected (8 each) in Thulokunwar (T01, T02, T02, T03, T04, T05, T06, T07 and T08), Kasiya (K01, K02, K03, K04, K05, K06, K07, and K08) and Pachgawa (P01, P02, P03, P04, P05, P06, P07, and P08) for sampling of groundwater samples.

Collection of Water Samples

A total of 24 groundwater samples (250 ml) were collected from tubewells, varying in depth from 6 to 48.7m in September, 2011. Tubewells were flushed well before collecting the samples to remove the stagnant water. One set of sample was collected with 1 ml of conc. HNO3 for metal analysis while the other set was collected without any preservative and used for analysis of major anions. The samples were stored at 4°C until analysis. Depth to water table (DWT) of the tubewells was measured by inserting a steel measuring tape in tubewell.

Physico-Chemical Analysis

Water pH and oxidation reduction potential (ORP) were measured in-situ using Hanna made portable kit (HI 98121 waterproof pH/ORP/Temp). The electrical conductivity (EC) was recorded using Milwaukee sharp EC tester (M111751). SO42-, PO43- and NO3- were analyzed spectrophotometerically (Systronics Visiscan 167) while HCO3 - was analyzed titrimetrically. Water samples were digested with conc. HNO3 following standard methods [25]. Fe, Mn, Ni, Pb and Zn were quantified by Flame Atomic Absorption Spectrometer (Perkin Elmer AAnalyst 800). Total As was analyzed using hydride generator atomic absorption spectrometer (HG-AAS). As (III) was completely reduced to As (V) by adding 10 ml each of conc. HCl and 5% (w/v) reducing solutions (KI and ascorbic acid) in 50 ml volumetric flask and incubating it at room temperature for 45 min. The samples were diluted with 1.5% HCl to the mark and were analyzed by AAS-HG. The detail protocol about As analysis by GH-AAS is described elsewhere [12].

Results and Discussion

Characteristics of Groundwater Parameters

Descriptive statistics has been used to characterize the composition of groundwater (Table 1). Maximum pH (7.60) was recorded in Kasiya followed by Panchgawa (7.44) and Thulokunwar (7.41). The slightly alkaline pH values as recorded for all the sites were also well within the limit of Nepal drinking water quality standards (NDWQS). The EC values of groundwater

 hydro fig 4.1

Figure 1. Map of Nawalparasi district showing study area and sampling locations.

hydro table  4.1

Table 1. Descriptive statistics of groundwater quality parameters.

hydro fig 4.2

Figure 2. Location-wise concentration of total As in Nawalparasi.

hydro fig 4.3

Figure 3. As concentration in tubewells of Nawalparasi by depth.

varied from 450 to 1438 μScm-1, 449 to 624 μScm-1 and 391 to 571 μScm-1 in Thulokunwar, Kasiya and Panchgawa, respectively and were in compliance with the NDWQS limit of 1500 μScm-1. The groundwater in this area showed low ORP values and ranged from -124 to -39 mv, -99 to 97 mv and -129 to -39 mv in Thulokunwar, Kasiya and Panchgawa, respectively. Negative ORP value in groundwater is the indication of the reducing groundwater conditions which promotes As release from aquifer sediments [26]. SO4 2−, NO3- and PO43- are major constituents of groundwater, and their elevated concentration is indicative of anthropogenic pollution. The concentration of SO42− ranged from 0.06 to 1.14 mgL-1, 0.06 to 0.58 mgL-1 and 0.09 to 0.22 mgL-1 at Thulokunwar, Kasiya and Panchgawa, respectively and was well within the limit of NDWQS (250 mgL-1). Low concentration of SO4 2−, NO3- in the groundwater of present study may be due to sulfate reduction and nitrate reduction in groundwater [27,28]. The average concentration of PO4 3- was either very low or not detectable, and well within the limit of Bureau of Indian Standard (BIS) for drinking water. However, no limit is set for PO43- by NDWQS. The SO4 2-, NO3- and PO43- concentrationdetected in groundwater of this area is consistent with previous study from Nepal [12 17,24,29,30] and other part of world [31-33].

Bicarbonate (HCO3 -) represents major source of alkalinity in groundwater. It varied from 402.5 to 920 mgL-1, 402.5 to 575 mgL-1 and 402 to 575 mgL-1 in Thulokunwar, Kasiya and Panchgawa, respectively. The concentration of Fe was maximum in Thulokunwar (3.71 mgL-1) exceeding permissible limit of WHO (0.30 mgL-1) in all samples and within the limit of NDWQS (3 mgL-1) for Fe. Mn was very low at all the sites having its mean 0.21, 0.07 and 0.07 mgL-1 in Thulokunwar, Kasiya and Panchgawa, respectively. The presence of high concentration of Fe and low concentration of Mn indicates the reducing condition of the groundwater. High level of Pb was recorded at all the sites which ranged from 0.08 to 0.13 mgL-1, 0.09 to 0.16 mgL-1 and 0.13 to 0.16 mgL-1 in Thulokunwar, Kasiya and Panchgawa, respectively. Whereas Ni was present in the limit of WHO, Zn was, however, very low and was below the WHO limit of drinking water.

The average As-concentration was 10 to 100 times more than the WHO limit (0.01 mgL-1) and 25-30 times higher than Nepalese interim standard (NIS, 0.05 mgL-1) limit for As in drinking water (Figure 2). All the tubewells analyzed under present study contained As more or less, and none of the tubewell found to be free from As. While sampling in Nawalparasi, we found that the local people are forced to drink the As-contaminated water in absence of alternate option for drinking water. They are at high risk of several As-poisoning. Hence, the Department of Water Supply and Sewerage and Ramgram municipality, who manage the water supply in the area, are also advised to supply the treated tank water (As-free) in these villages to protect the health of local residents.

Higher concentration of As was recorded in Thulokunwar and Panchgawa in the east and west side, respectively as compared to Kasiya where high As was recorded in the middle. The level of As varied from 0.16 mgL-1 in Kasiya to 1.34 mgL-1 in Thulokunwar with majority of the tubewell containing >0.01 mgL-1 of As, indicating that distribution pattern of As is highly heterogeneous. The As-concentration varied even within the closely located tubewells (Figure 2). The differences in mean As level may probably be due to geospatial conditions especially to latitudinal and longitudinal differences. This may also possibly be due to the vertical distribution of aquiclude such as dyke, fault or some unconformity etc. in the groundwater aquifer limiting the movement of As from one tubewell to another. Our finding is in agreement with previous study reported across the world [16, 18, 34, 35]. High level of As was detected in Thulokunwar and Panchgawa located adjacent to third grade rivers (rivers whose water source areas lie in Siwalik range) i.e. Bhaluhi and Jharahi river. The possible reason for high As content in such villages may be the sedimentary rocks of Siwalik origin brought by these third grade rivers. These rivers are much smaller in size and have insignificant or no flow during dry periods. Flash flood during monsoon results in the deposition of transported sediments in the nearby area. However, the exact reason as to why the As levels vary greatly within the closely located tubewells, is yet to be ascertained.

As Concentration of Tubewells as a Function Of Depth

The level of As changed with the tubewell depth (Figure 3). Tubewells with lesser depth (≤ 20 meter) had elevated As concentrations in contrast to those with higher depth (>20 meter) having lower As. However, both the tubewell types (with lesser or higher depth) contained As levels higher than the WHO and NIS limits of As for drinking water. Whereas maximum As (1.18 mgL-1) was detected for the tubewell with the depth of 19.81m, tubewells with the depth of 27.4m, however, had the least As (0.14 mgL-1). In general, the level of As varied inversely with the depth of the tubewells, suggesting that As concentration depends more or less on the tubewell depth.

As Concentration Vs Depth of Water Table

DWT gradually increased from east to west direction from Thulokunwar to Panchgawa through Kasiya. It ranged as high as 0.52m below the ground level (bgl) at Thulokunwar in the east to as low as 5.55m bgl at Panchgawa in the west. High DWT at Panchgawa village lies towards Jharahi River. Water table also increased from north (2.52m bgl) to south (4.02m bgl). The annual fluctuation of water table was low (0.34m bgl) and high (1.58m bgl) at Panchgawa and Thulokunwar, respectively. The eastern part is having the high water level fluctuation (1.23m bgl) which may be due to sand/gravel bearing aquifers of Panchgawa. The deepest DWT was recorded at Panchgawa located at an altitude of 100m (7.5m bgl) while the shallow DWT was common for Thulokunwar located at an altitude of 96m. The DWT decreased from west to east and also from south to north.


The spatial distribution of As showed significant variation even in closely located tubewells in the groundwater of Nawalparasi. The difference in As concentration over space may probably be due to spatial conditions such as latitude, longitude and depth of tubewells. All the tubewells contained many fold higher As than the WHO limit. However, none of the tubewells was free from As contamination, suggesting that all are unsafe and there is an urgent need to educate the residents of Nawalparasi to avoid drinking such contaminated tubewell water forthwith. The Department of Water Supply and Sewerage and Ramgram municipality, who manage the water supply in the area, are also advised to supply the treated tank water (As-free) in these three villages.



1.Smedley P L, Kinniburgh DG. A review of the source, behavior and distribution of arsenic in natural waters. Appl Geochem. 2002, 17(5): 517–568.

2.Agusa T, Kunito T, Fujihara J, Kubota R, Minh TB et al. Contamination by arsenic and other trace elements in tubewell water and its risk assessment to humans in Hanoi, Vietnam Environ Pollut. 2006, 139(1): 95–106.

3.Berg M, Tran HC, Nguyen TC, Pham HV, Schertenleib R et al. Arsenic contamination of groundwater and drinking water in Vietnam: a human health threat. Environ Sci Technol. 2001, 35(13): 2621–2626.

4.Guo H, Wang Y. Geochemical characteristics of shallow groundwater in Datong basin, northwestern China. J Geochem Explor. 2005, 87(3): 109–120.

5.Harvey CF, Ashfaque KN, Yu W, Badruzzaman ABM, Ali MA et al. Groundwater dynamics and arsenic contamination in Bangladesh. Chem Geol. 2006, 228: 112–136.

6.Hossain MF. Arsenic contamination in Bangladesh an overview. Agr Ecosyst Environ. 2006, 113(2006): 1–16.

7.Nickson RT, McArthur JM, Shrestha B, Kyaw-Myint TO, Lowry D. Arsenic and other drinking water quality issues, Muzaffargarh District, Pakistan. Appl Geochem. 2005, 20(1): 55–68.

8.Haque S, Johannesson KH. Arsenic concentrations and speciation along a groundwater flow path: the Carrizo Sand aquifer, Texas, USA. Chem Geol. 2006, 228(1-3): 57–71.

9.Kim MJ, Nriagu J, Haack S. Arsenic species and chemistry in groundwater of southeast Michigan. Environ Pollut. 2002, 120(2): 379–390.

10.Peters SC, Blum JD, Karagas MR, Chamberlain CP, Sjostrom DJ. Sources and exposure of the New Hampshire population to arsenic in public and private drinking water supplies. Chem Geol. 2006, 228(1-3): 72–84.

11.Price RE, Pichler T. Abundance and mineralogical association of arsenic in the Suwannee Limestone (Florida): implications for arsenic release during water–rock interac-tion. Chem Geol. 2006, 228(1-3): 44–56.

12.Yadav IC, Devi NL, Mohan D, Shihua Q, Singh S. Assessment of groundwater quality with special reference to arsenic in Nawalparasi District, Nepal using multivariate statistical techniques. Environ Earth Sci. 2014, 72: 259–273.

13.Appleyard SJ, Angeloni J, Watkins R. Arsenic-rich groundwater in urban area experiencing drought and increasing population density Perth Australia. Appl Geochem. 2006, 21(1): 83–97.

14.Aiuppa A, Avino R, Brusca L, Caliro S, Chiodini G et al. Mineral control of arsenic content in thermal waters from volcano- hosted hydrothermal systems: insights from island of Ischia and Phlegrean fields (Campanian Volcanic Province, Italy). Chem Geol. 2006, 229(4): 313–330.

15.Nordstrom DK. Worldwide occurrences of arsenic in groundwater. Science. 2002, 296(5576): 2143–2145.

16.Singh SK, Singh CK, Kumar KS, Gupta R, Mukherjee S. Spatial– temporal monitoring of groundwater using multivariate statistical techniques in Bareilly district of Uttarpradesh, India. J Hydrol Hydromech. 2009, 57: 45–54.

17.Yadav IC, Devi NL, Singh S. Spatial and temporal variation in arsenic in the groundwater of upstream of Ganges River Basin, Nepal. Environ Earth Sci. 2015, 73:1265–1279.

18.van Geen A, Ahmed KM, Seddique AA, Shamsudduha M. Community wells to mitigate the arsenic crisis in Bangladesh. Bull. World Health Organization. 2003, 81(9): 632– 638.

19.Fendorf S, Michael HA, van Geen A. Spatial and temporal variations of groundwater arsenic in south and Southeast Asia. Science. 2010, 328: 1123-1127.

20.Yadav IC, Dhuldhaj UP, Mohan D, Singh S. Current status of groundwater arsenic and its impacts on health and mitigation measures in the Terai basin of Nepal: An overview. Environ Rev. 2011, 19: 55-67.

21.Hagen T, J B Auden. Reports on geological survey of Nepal Preliminary Reconnaissance. Denkscher Schweiz Natur Ges. 1969, 86:185.

22.Sharma CK. Geology of Nepal Himalaya and adjacent countries. Sangeeta Sharma. Kathmandu. 1990, 479.

23.Tibbitts GC, Ogilbee W, Sharma CK. Groundwater resource investigations in Lumbini zone, Western Terai, Nepal. Interim Administrative Report released by Capital project and Engineering Division, USAID Nepal, 1973.

24.Yadav IC, Singh S, Devi NL, Mohan D, Pahari M et al. Spatial distribution of arsenic in groundwater of Southern Nepal. Rev Environ Contamin Toxicol. 2012, 218: 125-140.

25.APHA. Standard methods for the examination of water and wastewater, 20th edn. American Public Health Association, American Water Works Association, Water Pollution Control Federation, Washington, DC. 1998.

26.Guo H, Yang SZ, Tang XH, Li Y, Shen ZL. Groundwater geochemistry and its implications for arsenic mobilization in shallow aquifers of the Hetao basin, Inner Mongolia. Sci Total Environ. 2008, 393(1): 131–144.

27.Kinniburgh DG, Smedley PL. Arsenic contamination of groundwater in Bangladesh. Summary. BGS Technical Report WC/00/19 Volume 1. British Geological Survey & Department of Public Health Engineering. 2001, 15.

28.Nickson RT, McArthur JM, Burgess WG, Ahmed KM, Ravenscroft P et al. Arsenic poisoning of Bangladesh groundwater. Nature. 1998, 395: 338.

29.Chapagain SK, Pandey VP, Shrestha S, Nakamura T, Kazama F. Assessment of deep groundwater quality in Kathmandu valley using multivariate statistical techniques. Water Air Soil Pollut. 2010, 210: 277–288.

30.Gurung JK, Hiroaki I, Khadka MS. Geological and geochemical examination of arsenic contamination in groundwater in the Holocene Terai Basin, Nepal. Environ Geol. 2005, 49: 98–113.

31.Ahmed KM, Bhattacharya P, Hasan MA, Akhter SH, Alam SMM et al. Arsenic enrichment in groundwater of the alluvial aquifers in Bangladesh: an overview. Appl Geochem. 2004, 19(2): 181–200.

32.van Geen A, Zheng Y, Cheng Z, Aziz Z, Horneman A et al. A transect of groundwater and sediment properties in Araihazar, Bangladesh: further evidence of decoupling between As and Fe mobilization. Chem Geol. 2006, 228(1-3): 85–96.

33.Reza AHMS, Jean JS, Lee MK, Yang HJ, Liu CC. Arsenic enrichment and mobilization in the Holocene alluvial aquifers of the Chapai–Nawabganj district, Bangladesh: a geochemical and statistical study. Appl Geochem. 2010, 25(8): 1280–1289.

34.Satpathy KK, Mohanty AK, Sahu G, Sarguru S, Sarkar SK et al. Spatio-temporal variation in physicochemical properties of coastal waters off Kalpakkam, southeast coast of India, during summer, pre-monsoon and post-monsoon period. Environ Monit Assess. 2011, 180(1-4): 41–62.

35.Guo H, Zhang B, Li Y, Berner Z, Tang X et al. Hydrogeological and biogeochemical constrains of arsenic mobilization in shallow aquifers from the Hetao basin, Inner Mongolia. Environ Pollut. 2011, 159(4): 876–883.

36.Government of Nepal. National Drinking Water Quality Standards, 2062 and National Drinking Water QualityStandard Implementation Guideline, 2062 Singh Durbar, Kathmandu, Nepal: Ministry of Land Reform Management, 2005.

37.WHO. Guideline for drinking-water quality, 3rd edn (Recommendations). World Health Organization, Geneva, 2004.

Cite this article: Yadav I C. Spatial Variation in Arsenic Concentration in the Groundwater of Nawalparasi District, Terai Nepal. J J Hydrology. 2015, 1(1): 004.

Contact Us:
TRAIL # 150 W
E-mail : info@jacobspublishers.com
Phone : 512-400-0398