Department Water Resources and Drinking Water

Concentrations of naturally occurring fluoride in groundwater exceeding the WHO guideline of 1.5 mg/L have been detected in many parts of Pakistan. This may lead to dental or skeletal fluorosis and thereby poses a potential threat to public health. Utilizing a total of 5483 fluoride concentrations, comprising 2160 new measurements as well as those from other sources, we have applied machine learning techniques to predict the probability of fluoride in groundwater in Pakistan exceeding 1.5 mg/L at a 250 m spatial resolution. Climate, soil, lithology, topography, and land cover parameters were identified as effective predictors of high fluoride concentrations in groundwater. Excellent model performance was observed in a random forest model that achieved an Area Under the Curve (AUC) of 0.92 on test data that were not used in modeling. The highest probabilities of high fluoride concentrations in groundwater are predicted in the Thar Desert, Sargodha Division, and scattered along the Sulaiman Mountains. Applying the model predictions to the population density and accounting for groundwater usage in both rural and urban areas, we estimate that about 13 million people may be at risk of fluorosis due to consuming groundwater with fluoride concentrations >1.5 mg/L in Pakistan, which corresponds to ~6% of the total population. Both the fluoride prediction map and the health risk map can be used as important decision-making tools for authorities and water resource managers in the identification and mitigation of groundwater fluoride contamination.

The health of millions of people worldwide is negatively impacted by chronic exposure to elevated concentrations of geogenic fluoride in groundwater. Due to health effects including dental mottling and skeletal fluorosis, the World Health Organization maintains a maximum guideline of 1.5 mg/L in drinking water. As groundwater quality is not regularly tested in many areas, it is often unknown if the water in a given well or spring contains harmful levels of fluoride. Here we present a state-of-the-art global fluoride hazard map based on machine learning and over 400,000 fluoride measurements (10% of which >1.5 mg/L), which is then used to estimate the human population at risk. Hotspots indicated by the groundwater fluoride hazard map include parts of central Australia, western North America, eastern Brazil and many areas of Africa and Asia. Of the approximately 180 million people potentially affected worldwide, most reside in Asia (51–59% of total) and Africa (37-46% of total), with the latter representing 6.5% of the continent’s population. Africa also contains 14 of the top 20 affected countries in terms of population at risk. We also illuminate and discuss the key globally relevant hydrochemical and environmental factors related to fluoride accumulation.

Most people in Ghana have no or only basic access to safely managed water. Especially in rural areas, much of the population relies on groundwater for drinking, which can be contaminated with fluoride and lead to dental fluorosis. Children under the age of two are particularly susceptible to the adverse effects of fluoride and can retain 80-90% of a fluoride dose, compared to 60% in adults. Despite numerous local studies, no spatially continuous picture exists of the fluoride contamination across Ghana, nor is there any estimate of what proportion of the population is potentially exposed to unsafe fluoride levels. Here, we spatially model the probability of fluoride concentrations exceeding 1.0 mg/L in groundwater across Ghana to identify risk areas and estimate the number of children and adults exposed to unsafe fluoride levels in drinking water. We use a set of geospatial predictor variables with random forest modeling and evaluate the model performance through spatial cross-validation. We found that approximately 15% of the area of Ghana, mainly in the northeast, has a high probability of fluoride contamination. The total at-risk population is about 920,000 persons, or 3% of the population, with an estimated 240,000 children (0-9 years) in at-risk areas. In some districts, such as Karaga, Gushiegu, Tamale and Mion, 4 out of 10 children are potentially exposed to fluoride poisoning. Geology and high evapotranspiration are the main drivers of fluoride enrichment in groundwater. Consequently, climate change might put even greater pressure on the area's water resources. Our hazard maps should raise awareness and understanding of geogenic fluoride contamination in Ghana and can advise decision making at local levels to avoid or mitigate fluoride-related risks.

Naturally occurring, geogenic manganese (Mn) and iron (Fe) are frequently found dissolved in groundwater at concentrations that make the water difficult to use (deposits, unpleasant taste) or, in the case of Mn, a potential health hazard. Over 6000 groundwater measurements of Mn and Fe in Southeast Asia and Bangladesh were assembled and statistically examined with other physicochemical parameters. The machine learning methods random forest and generalized boosted regression modeling were used with spatially continuous environmental parameters (climate, geology, soil, topography) to model and map the probability of groundwater Mn > 400 μg/L and Fe > 0.3 mg/L for Southeast Asia and Bangladesh. The modeling indicated that drier climatic conditions are associated with a tendency of elevated Mn concentrations, whereas high Fe concentrations tend to be found in a more humid climate with elevated levels of soil organic carbon. The spatial distribution of Mn > 400 μg/L and Fe > 0.3 mg/L was compared and contrasted with that of the critical geogenic contaminant arsenic (As), confirming that high Fe concentrations are often associated with high As concentrations, whereas areas of high concentrations of Mn and As are frequently found adjacent to each other. The probability maps draw attention to areas prone to elevated concentrations of geogenic Mn and Fe in groundwater and can help direct efforts to mitigate their negative effects. The greatest Mn hazard is found in densely populated northwest Bangladesh and the Mekong, Red and Ma River Deltas of Cambodia and Vietnam. Widespread elevated Fe concentrations and their associated negative effects on water infrastructure pose challenges to water supply. The Mn and Fe prediction maps demonstrate the value of machine learning for the geospatial prediction modeling and mapping of groundwater contaminants as well as the potential for further constituents to be targeted by this novel approach.

While drinking water is known to create significant health risk in arsenic hazard areas, the role of exposure to arsenic through food intake is less well understood, including the impact of food trade. Using the best available datasets on crop production, irrigation, groundwater arsenic hazard, and international crop trade flows, we estimate that globally 17.2% of irrigated harvested area (or 45.2 million hectares) of 42 main crops are grown in arsenic hazard areas, contributing 19.7% of total irrigated crop production, or 418 million metric tons (MMT) per year of these crops by mass. Two-thirds of this area is dedicated to the major staple crops of rice, wheat, and maize (RWM) and produces 158 MMT per year of RWM, which is 8.0% of the total RWM production and 18% of irrigated production. More than 25% of RWM consumed in the South Asian countries of India, Pakistan, and Bangladesh, where both arsenic hazard and degree of groundwater irrigation are high, originate from arsenic hazard areas. Exposure to arsenic risk from crops also comes from international trade, with 10.6% of rice, 2.4% of wheat, and 4.1% of maize trade flows coming from production in hazard areas. Trade plays a critical role in redistributing risk, with the greatest exposure risk borne by countries with a high dependence on food imports, particularly in the Middle East and small island nations for which all arsenic risk in crops is imported. Intensifying climate variability and population growth may increase reliance on groundwater irrigation, including in arsenic hazard areas. Results show that RWM harvested area could increase by 54.1 million hectares (179% increase over current risk area), predominantly in South and Southeast Asia. This calls for the need to better understand the relative risk of arsenic exposure through food intake, considering the influence of growing trade and increased groundwater reliance for crop production.

Geogenic arsenic contamination in groundwaters poses a severe health risk to hundreds of millions of people globally. Notwithstanding the particular risks to exposed populations in the Indian sub-continent, at the time of writing, there was a paucity of geostatistically based models of the spatial distribution of groundwater hazard in India. In this study, we used logistic regression models of secondary groundwater arsenic data with research-informed secondary soil, climate and topographic variables as principal predictors generate hazard and risk maps of groundwater arsenic at a resolution of 1 km across Gujarat State. By combining models based on different arsenic concentrations, we have generated a pseudo-contour map of groundwater arsenic concentrations, which indicates greater arsenic hazard (> 10 μg/L) in the northwest, northeast and south-east parts of Kachchh District as well as northwest and southwest Banas Kantha District. The total number of people living in areas in Gujarat with groundwater arsenic concentration exceeding 10 μg/L is estimated to be around 122,000, of which we estimate approximately 49,000 people consume groundwater exceeding 10 µg/L. Using simple previously published dose–response relationships, this is estimated to have given rise to 700 (prevalence) cases of skin cancer and around 10 cases of premature avoidable mortality/annum from internal (lung, liver, bladder) cancers—that latter value is on the order of just 0.001% of internal cancers in Gujarat, reflecting the relative low groundwater arsenic hazard in Gujarat State.

Naturally occurring arsenic in groundwater affects millions of people worldwide. We created a global prediction map of groundwater arsenic exceeding 10 micrograms per liter using a random forest machine-learning model based on 11 geospatial environmental parameters and more than 50,000 aggregated data points of measured groundwater arsenic concentration. Our global prediction map includes known arsenic-affected areas and previously undocumented areas of concern. By combining the global arsenic prediction model with household groundwater-usage statistics, we estimate that 94 million to 220 million people are potentially exposed to high arsenic concentrations in groundwater, the vast majority (94%) being in Asia. Because groundwater is increasingly used to support growing populations and buffer against water scarcity due to changing climate, this work is important to raise awareness, identify areas for safe wells, and help prioritize testing.

Geogenic groundwater arsenic (As) contamination is pervasive in many aquifers in south and southeast Asia. It is feared that recent increases in groundwater abstractions could induce the migration of high-As groundwaters into previously As-safe aquifers. Here we study an As-contaminated aquifer in Van Phuc, Vietnam, located ~10 km southeast of Hanoi on the banks of the Red River, which is affected by large-scale groundwater abstraction. We used numerical model simulations to integrate the groundwater flow and biogeochemical reaction processes at the aquifer scale, constrained by detailed hydraulic, environmental tracer, hydrochemical and mineralogical data. Our simulations provide a mechanistic reconstruction of the anthropogenically induced spatiotemporal variations in groundwater flow and biogeochemical dynamics and determine the evolution of the migration rate and mass balance of As over several decades. We found that the riverbed–aquifer interface constitutes a biogeochemical reaction hotspot that acts as the main source of elevated As concentrations. We show that a sustained As release relies on regular replenishment of river muds rich in labile organic matter and reactive iron oxides and that pumping-induced groundwater flow may facilitate As migration over distances of several kilometres into adjacent aquifers.

Geogenic arsenic (As) contamination of groundwater poses a major threat to global health, particularly in Asia. To mitigate this exposure, groundwater is increasingly extracted from low-As Pleistocene aquifers. This, however, disturbs groundwater flow and potentially draws high-As groundwater into low-As aquifers.
Here we report a detailed characterisation of the Van Phuc aquifer in the Red River Delta region, Vietnam, where high-As groundwater from a Holocene aquifer is being drawn into a low-As Pleistocene aquifer. This study includes data from eight years (2010–2017) of groundwater observations to develop an understanding of the spatial and temporal evolution of the redox status and groundwater hydrochemistry.
Arsenic concentrations were highly variable (0.5-510 μg/L) over spatial scales of <200 m. Five hydro(geo)chemical zones (indicated as A to E) were identified in the aquifer, each associated with specific As mobilisation and retardation processes. At the riverbank (zone A), As is mobilised from freshly deposited sediments where Fe(III)-reducing conditions occur. Arsenic is then transported across the Holocene aquifer (zone B), where the vertical intrusion of evaporative water, likely enriched in dissolved organic matter, promotes methanogenic conditions and further release of As (zone C). In the redox transition zone at the boundary of the two aquifers (zone D), groundwater arsenic concentrations decrease by sorption and incorporations onto Fe(II) carbonates and Fe(II)/Fe(III) (oxyhydr)oxides under reducing conditions. The sorption/incorporation of As onto Fe(III) minerals at the redox transition and in the Mn(IV)-reducing Pleistocene aquifer (zone E) has consistently kept As concentrations below 10 μg/L for the studied period of 2010–2017, and the location of the redox transition zone does not appear to have propagated significantly. Yet, the largest temporal hydrochemical changes were found in the Pleistocene aquifer caused by groundwater advection from the Holocene aquifer. This is critical and calls for detailed investigations.

An index aquifer vulnerability study from western Canada (left) compared with a new logistic regression map of these vulnerability index values on the online GAP platform (right). The red colour shows areas with the highest vulnerability. The green areas are less vulnerable or adequately protected from surface contamination.

For about the past eight decades, high concentrations of naturally occurring fluoride have been detected in groundwater in different parts of India. The chronic consumption of fluoride in high concentrations is recognized to cause dental and skeletal fluorosis. We have used the random forest machine-learning algorithm to model a data set of 12 600 groundwater fluoride concentrations from throughout India along with spatially continuous predictor variables of predominantly geology, climate, and soil parameters. Despite only surface parameters being available to describe a subsurface phenomenon, this has produced a highly accurate prediction map of fluoride concentrations exceeding 1.5 mg/L at 1 km resolution throughout the country. The most affected areas are the northwestern states/territories of Delhi, Gujarat, Haryana, Punjab, and Rajasthan and the southern states of Andhra Pradesh, Karnataka, Tamil Nadu, and Telangana. The total number of people at risk of fluorosis due to fluoride in groundwater is predicted to be around 120 million, or 9% of the population. This number is based on rural populations and accounts for average rates of groundwater consumption from nonmanaged sources. The new fluoride hazard and risk maps can be used by authorities in conjunction with detailed groundwater utilization information to prioritize areas in need of mitigation measures.

It has been reported that persistent cyanide pollution occurs in artisanal small-scale gold mining (ASGM)-affected catchment areas in Burkina Faso. In the present study, the logistic regression method was employed to identify the factors that influence the spatial distribution of cyanide pollution as well as to predict the cyanide pollution map risk at catchment level. Soil samples were collected from two ASGM sites in the northern Zougnazagmiline ("North") site and southern Galgouli ("South") site parts of Burkina Faso, covering areas of 22 km² and 20 km², respectively. Free cyanide concentration in each sample was measured. It was shown that the spatial distribution of cyanide was solely controlled by the soil type in Zougnazagmiline and both the soil type and electric conductivity in Galgouli. On the other hand, the cyanidation zones within the two catchments were the places where the highest risk of cyanide pollution occurs, with probabilities of 0.8 and 1 in Zougnazagmiline and Galgouli, respectively. > 20% of the settled area in the Zougnazagmiline and 5% of that in Galgouli were exposed to cyanide pollution. Logistic regression was able to reliably predict cyanide contamination in areas affected by ASGM. The model could be useful for decision-makers to plan ASGM-site decontamination.

The semi-arid Sahel regions of West Africa rely heavily on groundwater from shallow to moderately deep (<100 m b.g.l.) crystalline bedrock aquifers for drinking water production. Groundwater quality may be affected by high geogenic arsenic (As) concentrations (>10 μg/L) stemming from the oxidation of sulphide minerals (pyrite, arsenopyrite) in mineralised zones. These aquifers are still little investigated, especially concerning groundwater residence times and the influence of the annual monsoon season on groundwater chemistry. To gain insights on the temporal aspects of As contamination, we have used isotope tracers (noble gases, ³H, stable water isotopes (²H, ¹⁸O)) and performed hydrochemical analyses on groundwater abstracted from tube wells and dug wells in a small study area in southwestern Burkina Faso. Results revealed a great variability in groundwater properties (e.g. redox conditions, As concentrations, water level, residence time) over spatial scales of only a few hundred metres, characteristic of the highly heterogeneous fractured underground. Elevated As levels are found in oxic groundwater of circum-neutral pH and show little relation with any of the measured parameters. Arsenic concentrations are relatively stable over the course of the year, with little effect seen by the monsoon. Groundwater residence time does not seem to have an influence on As concentrations, as elevated As can be found both in groundwater with short (<50 a) and long (>10³ a) residence times as indicated by ³He/⁴He ratios spanning three orders of magnitude. These results support the hypothesis that the proximity to mineralised zones is the most crucial factor controlling As concentrations in the observed redox/pH conditions. The existence of very old water portions with residence times >10³ years already at depths of <50 m b.g.l. is a new finding for the shallow fractured bedrock aquifers of Burkina Faso, suggesting that overexploitation of these relatively low-yielding aquifers may be an issue in the future.

This paper presents a first integrated survey on the occurrence and distribution of geogenic contaminants in groundwater resources of Western Amazonia in Peru. An increasing number of groundwater wells have been constructed for drinking water purposes in the last decades; however, the chemical quality of the groundwater resources in the Amazon region is poorly studied. We collected groundwater from the regions of Iquitos and Pucallpa to analyze the hydrochemical characteristics, including trace elements. The source aquifer of each well was determined by interpretation of the available geological information, which identified four different aquifer types with distinct hydrochemical properties. The majority of the wells in two of the aquifer types tap groundwater enriched in aluminum, arsenic, or manganese at levels harmful to human health. Holocene alluvial aquifers along the main Amazon tributaries with anoxic, near pH-neutral groundwater contained high concentrations of arsenic (up to 700 μg/L) and manganese (up to 4 mg/L). Around Iquitos, the acidic groundwater (4.2 ≤ pH ≤ 5.5) from unconfined aquifers composed of pure sand had dissolved aluminum concentrations of up to 3.3 mg/L. Groundwater from older or deeper aquifers generally was of good chemical quality. The high concentrations of toxic elements highlight the urgent need to assess the groundwater quality throughout Western Amazonia.

Arsenic-contaminated aquifers are currently estimated to affect ~150 million people around the world. However, the full extent of the problem remains elusive. This is also the case in Pakistan, where previous studies focused on isolated areas. Using a new data set of nearly 1200 groundwater quality samples throughout Pakistan, we have created state-of-the-art hazard and risk maps of arsenic-contaminated groundwater for thresholds of 10 and 50 μg/liter. Logistic regression analysis was used with 1000 iterations, where surface slope, geology, and soil parameters were major predictor variables. The hazard model indicates that much of the Indus Plain is likely to have elevated arsenic concentrations, although the rest of the country is mostly safe. Unlike other arsenic-contaminated areas of Asia, the arsenic release process in the arid Indus Plain appears to be dominated by elevated-pH dissolution, resulting from alkaline topsoil and extensive irrigation of unconfined aquifers, although pockets of reductive dissolution are also present. We estimate that approximately 50 million to 60 million people use groundwater within the area at risk, with hot spots around Lahore and Hyderabad. This number is alarmingly high and demonstrates the urgent need for verification and testing of all drinking water wells in the Indus Plain, followed by appropriate mitigation measures.

Arsenic contamination in groundwater from crystalline basement rocks in West Africa has only been documented in isolated areas and presents a serious health threat in a region already facing multiple challenges related to water quality and scarcity. We present a comprehensive dataset of arsenic concentrations from drinking water wells in rural Burkina Faso (n = 1498), of which 14.6% are above 10 μg/L. Included in this dataset are 269 new samples from regions where no published water quality data existed. We used multivariate logistic regression with arsenic measurements as calibration data and maps of geology and mineral deposits as independent predictor variables to create arsenic prediction models at concentration thresholds of 5, 10 and 50 μg/L. These hazard maps delineate areas vulnerable to groundwater arsenic contamination in Burkina Faso. Bedrock composed of schists and volcanic rocks of the Birimian formation, potentially harbouring arsenic-containing sulphide minerals, has the highest probability of yielding groundwater arsenic concentrations > 10 μg/L. Combined with population density estimates, the arsenic prediction models indicate that ~ 560,000 people are potentially exposed to arsenic-contaminated groundwater in Burkina Faso. The same arsenic-bearing geological formations that are positive predictors for elevated arsenic concentrations in Burkina Faso also exist in neighbouring countries such as Mali, Ghana and Ivory Coast. This study's results are thus of transboundary relevance and can act as a trigger for targeted water quality surveys and mitigation efforts.

Arsenic-contaminated groundwater used for drinking in China is a health threat that was first recognized in the 1960s. However, because of the sheer size of the country, millions of groundwater wells remain to be tested in order to determine the magnitude of the problem. We developed a statistical risk model that classifies safe and unsafe areas with respect to geogenic arsenic contamination in China, using the threshold of 10 micrograms per liter, the World Health Organization guideline and current Chinese standard for drinking water. We estimate that 19.6 million people are at risk of being affected by the consumption of arsenic-contaminated groundwater. Although the results must be confirmed with additional field measurements, our risk model identifies numerous arsenic-affected areas and highlights the potential magnitude of this health threat in China.

Arsenic contamination of shallow groundwater is among the biggest health threats in the developing world. Targeting uncontaminated deep aquifers is a popular mitigation option although its long-term impact remains unknown. Here we present the alarming results of a large-scale groundwater survey covering the entire Red River Delta and a unique probability model based on three-dimensional Quaternary geology. Our unprecedented dataset reveals that ∼7 million delta inhabitants use groundwater contaminated with toxic elements, including manganese, selenium, and barium. Depth-resolved probabilities and arsenic concentrations indicate drawdown of arsenic-enriched waters from Holocene aquifers to naturally uncontaminated Pleistocene aquifers as a result of > 100 years of groundwater abstraction. Vertical arsenic migration induced by large-scale pumping from deep aquifers has been discussed to occur elsewhere, but has never been shown to occur at the scale seen here. The present situation in the Red River Delta is a warning for other As-affected regions where groundwater is extensively pumped from uncontaminated aquifers underlying high arsenic aquifers or zones.

Contamination of groundwaters with geogenic arsenic poses a major health risk to millions of people. Although the main geochemical mechanisms of arsenic mobilization are well understood, the worldwide scale of affected regions is still unknown. In this study we used a large database of measured arsenic concentration in groundwaters (around 20,000 data points) from around the world as well as digital maps of physical characteristics such as soil, geology, climate, and elevation to model probability maps of global arsenic contamination. A novel rule-based statistical procedure was used to combine the physical data and expert knowledge to delineate two process regions for arsenic mobilization: “reducing” and “high-pH/oxidizing”. Arsenic concentrations were modeled in each region using regression analysis and adaptive neuro-fuzzy inferencing followed by Latin hypercube sampling for uncertainty propagation to produce probability maps. The derived global arsenic models could benefit from more accurate geologic information and aquifer chemical/physical information. Using some proxy surface information, however, the models explained 77% of arsenic variation in reducing regions and 68% of arsenic variation in high-pH/oxidizing regions. The probability maps based on the above models correspond well with the known contaminated regions around the world and delineate new untested areas that have a high probability of arsenic contamination. Notable among these regions are South East and North West of China in Asia, Central Australia, New Zealand, Northern Afghanistan, and Northern Mali and Zambia in Africa.

The use of groundwater with high fluoride concentrations poses a health threat to millions of people around the world. This study aims at providing a global overview of potentially fluoride-rich groundwaters by modeling fluoride concentration. A large database of worldwide fluoride concentrations as well as available information on related environmental factors such as soil properties, geological settings, and climatic and topographical information on a global scale have all been used in the model. The modeling approach combines geochemical knowledge with statistical methods to devise a rule-based statistical procedure, which divides the world into 8 different “process regions”. For each region a separate predictive model was constructed. The end result is a global probability map of fluoride concentration in the groundwater. Comparisons of the modeled and measured data indicate that 60−70% of the fluoride variation could be explained by the models in six process regions, while in two process regions only 30% of the variation in the measured data was explained. Furthermore, the global probability map corresponded well with fluorotic areas described in the international literature. Although the probability map should not replace fluoride testing, it can give a first indication of possible contamination and thus may support the planning process of new drinking water projects.

Arsenic contamination of groundwater resources threatens the health of millions of people worldwide, particularly in the densely populated river deltas of Southeast Asia. Although many arsenic-affected areas have been identified in recent years, a systematic evaluation of vulnerable areas remains to be carried out. Here we present maps pinpointing areas at risk of groundwater arsenic concentrations exceeding 10 g l^-1. These maps were produced by combining geological and surface soil parameters in a logistic regression model, calibrated with 1,756 aggregated and geo-referenced groundwater data points from the Bengal, Red River and Mekong deltas. We show that Holocene deltaic and organic-rich surface sediments are key indicators for arsenic risk areas and that the combination of surface parameters is a successful approach to predict groundwater arsenic contamination. Predictions are in good agreement with the known spatial distribution of arsenic contamination, and further indicate elevated risks in Sumatra and Myanmar, where no groundwater studies exist.