Next Article in Journal
The Cost-Effectiveness of Lowering Permissible Noise Levels Around U.S. Airports
Next Article in Special Issue
A Meta-Analysis of Risk Factors for Post-Traumatic Stress Disorder (PTSD) in Adults and Children after Earthquakes
Previous Article in Journal
Correction: Nabeshima, Y., et al. Analysis of Japanese Articles about Suicides Involving Charcoal Burning or Hydrogen Sulfide Gas. Int. J. Environ. Res. Public Health 2016, 13, 1013
Previous Article in Special Issue
Integration of GIS, Electromagnetic and Electrical Methods in the Delimitation of Groundwater Polluted by Effluent Discharge (Salamanca, Spain): A Case Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hazard Ranking Method for Populations Exposed to Arsenic in Private Water Supplies: Relation to Bedrock Geology

1
Environmental Epidemiology Group, Centre for Radiation, Chemical and Environmental Hazards, Public Health England (PHE), Chilton, Oxfordshire OX11 0RQ, UK
2
London School of Hygiene and Tropical Medicine, Keppel Street, London WC1E 7HT, UK
3
Inorganic Geochemistry, Centre for Environmental Geochemistry, British Geological Survey, Nottingham NG12 5GG, UK
4
Statistics, Modelling and Economics Department, PHE, 61 Colindale Avenue, London NW9 5EQ, UK
5
Cornwall Council, Environmental Protection Team, Public Health and Protection, Camborne, Cornwall TR14 8SX, UK
6
School of Earth and Environmental Sciences, and Williamson Research Centre for Molecular Environmental Science, University of Manchester, Manchester M13 9PL, UK
*
Author to whom correspondence should be addressed.
Retired.
Int. J. Environ. Res. Public Health 2017, 14(12), 1490; https://doi.org/10.3390/ijerph14121490
Submission received: 22 September 2017 / Revised: 10 November 2017 / Accepted: 24 November 2017 / Published: 1 December 2017
(This article belongs to the Special Issue Natural Hazards and Public Health: A Systems Approach)

Abstract

:
Approximately one million people in the UK are served by private water supplies (PWS) where main municipal water supply system connection is not practical or where PWS is the preferred option. Chronic exposure to contaminants in PWS may have adverse effects on health. South West England is an area with elevated arsenic concentrations in groundwater and over 9000 domestic dwellings here are supplied by PWS. There remains uncertainty as to the extent of the population exposed to arsenic (As), and the factors predicting such exposure. We describe a hazard assessment model based on simplified geology with the potential to predict exposure to As in PWS. Households with a recorded PWS in Cornwall were recruited to take part in a water sampling programme from 2011 to 2013. Bedrock geologies were aggregated and classified into nine Simplified Bedrock Geological Categories (SBGC), plus a cross-cutting “mineralized” area. PWS were sampled by random selection within SBGCs and some 508 households volunteered for the study. Transformations of the data were explored to estimate the distribution of As concentrations for PWS by SBGC. Using the distribution per SBGC, we predict the proportion of dwellings that would be affected by high concentrations and rank the geologies according to hazard. Within most SBGCs, As concentrations were found to have log-normal distributions. Across these areas, the proportion of dwellings predicted to have drinking water over the prescribed concentration value (PCV) for As ranged from 0% to 20%. From these results, a pilot predictive model was developed calculating the proportion of PWS above the PCV for As and hazard ranking supports local decision making and prioritization. With further development and testing, this can help local authorities predict the number of dwellings that might fail the PCV for As, based on bedrock geology. The model presented here for Cornwall could be applied in areas with similar geologies. Application of the method requires independent validation and further groundwater-derived PWS sampling on other geological formations.

Graphical Abstract

1. Introduction

Domestic Private Water Supplies (PWS) provide water used for drinking and cooking to a proportion of the population varying according to the availability of alternative sources and distribution networks and householder attitudes. This proportion is estimated at approximately 1% in the UK [1]; from 10% to 30% in the United States [2], and over 90% in Bangladesh [3,4]. Groundwater (GW) is a common source for PWS, and its contamination by chemicals, of either geogenic or anthropogenic origin, can lead to preventable population exposure to health hazards [5,6,7,8]. In particular, inorganic arsenic is a recognised carcinogen and toxicant [9,10,11] commonly present in GW [12] and arguably the waterborne inorganic chemical with the highest global detrimental impact on public health [13].
Environmental Public Health Tracking (EPHT) can be defined as “the ongoing collection, integration, analysis and interpretation of data about environmental hazards, exposure to environmental hazards, and human health effects potentially related to exposure to environmental hazards. It includes dissemination of information learned from these data and implementation of strategies and actions to improve and protect public health” [14]. A key distinction between EPHT and traditional surveillance is the emphasis on data integration across hazard, exposure and health information systems, therefore an EPHT approach to arsenic (As) in PWS could have an initial goal to characterise and monitor population risk represented by As. This approach has been applied to public health monitoring of As in PWS in the USA [15,16,17,18,19]. In Louisiana, US, the EPHT approach has been used to prioritise parishes for further risk assessment and testing of As in PWS [16]. The EPHT approach could be applied to PWS in other geographical areas [20].
An EPHT programme for prevention of health effects deriving from exposure to As from PWS should consider information that can be obtained from the examination of a possible correlation between bedrock type and As concentration in GW-sourced PWS, as this could contribute to the estimation of population hazard, and related exposure and risk attributable to As. Geostatistical approaches to predicting GW arsenic concentrations have recently been widely implemented [21,22,23,24,25,26] and have been found to be useful particularly in shallow sediment/sedimentary-rock hosted aquifers. Ayotte et al. [27] suggested a possible geological source as a predictor of high As in New England PWS. Bedrock aquifers were reported to have higher As concentrations than unconsolidated aquifers. Arsenic concentrations in water from bedrock wells in variably calcareous metasedimentary rocks were significantly higher (p = 0.001) than water in bedrock wells in the combined differentiated metasedimentary and felsic igneous groups [27].
Whilst larger commercial and multi-use PWS in England are subject to risk assessment and testing for compliance to PWS Regulations (2016) [28] by local authorities, by 2016 68% of PWS had been risk assessed [29]. Single domestic PWS are not subject to risk assessment, unless requested by the owner or occupier. Water testing and treatment varies by area and in some regions potential As risk is unknown [29].
GW in areas of metalliferous mineralization is known to be particularly vulnerable to release of As. Metal sulphide minerals and their weathering products can release As and other trace elements as a result of mineral oxidation and desorption reactions. Release under anoxic GW conditions as a result of dissolution or desorption from metal oxides has also been well-documented [12]. Cornwall, in south west England, is known to have geologically enhanced As concentrations due to the presence of natural Cu-Sn-Zn-As mineralisation of rocks, which has been locally exacerbated by historical mining and ore processing activities. A water sampling programme of 497 small domestic PWS supplying 508 households was conducted in Cornwall during 2011–2013 and found that up to 35% of supplies had exceedances of one or more prescribed concentration or value (PCV) of a range of chemicals and that 20% of households had one or more exceedance of health-based values for drinking water (DW), including 12% above the manganese PCV [30,31]. Six percent of households failed the As PCV: the PCV and WHO guideline value [32] for As in PWS is 10 μg/L [32,33]. Biomonitoring of residents with PWS confirmed exposure to As through this exposure route, with high concentrations of As found in the urine, hair and toenails of some residents [34,35]. For PWS where DW was greater than 10 μg/L As, a strong correlation between DW As and urinary inorganic As was reported [34].
The confirmation of exposure to PWS As means that characterising As GW hazard, and hence population exposure in England, is an appropriate first step for a national programme of PWS tracking of chemicals in DW. This could subsequently address the related exposure hazard and health risk [36].
Hattis and Burmaster [37] argue the case for using distribution-based analysis in exposure assessments to examine variability and uncertainty. A distribution-based approach could be used to explore known relationships between geology and As concentrations in water, and then predict exposure in areas not systematically monitored.
The aim of this paper is to develop a methodology for a population-based hazard assessment and ranking of hazard for As in PWS used for drinking, based on geological bedrock classification. We estimate likely numbers of households at risk and how to best find those at risk based on bedrock geology.

2. Materials and Methods

2.1. Study Design

We conducted an observational, cross-sectional study of the population of Cornish PWS users. A model of exposure assessment predicts potential dwelling exposure through PWS based on the As distributions measured per simplified bedrock geological classification (SBGC). The method follows the model used to predict radon concentrations in dwellings in the UK by the National Radiological Protection Board (NRPB) [38], subsequently Public Health England (PHE) [39].

2.2. Study Population

The population was defined as users of DW originating from single domestic PWS in the county of Cornwall, in south west England, UK. Private well records used for access to users of PWS for domestic use were obtained by PHE in the form of a database compiled by the local authority, according to the provision of the PWS Regulations [28,33]. Residents of Cornwall were eligible; domestic use of the water was the primary focus, although in some instances it was used for other purposes (e.g., farms, dairy, food production, caravan parks or small rural businesses).

2.3. Rationale for Simplified Geological Classification

The rationale for the sampling design aimed to lead to an effective public health management approach, by integration of environmental (bedrock) and public health (concerning population groups) information for assessment of water results, as part of an EPHT proof of concept study. A SBGC was devised to approximate the geological features across areas sufficiently large to characterise potential population exposure to chemicals via DW, as described in Ander et al. [30]. The Cornwall area was divided into nine broad geological categories, defined on the basis of rock types (igneous/metamorphic/sedimentary), mineralogical characteristics (e.g., sulphide-bearing/slaty) and stratigraphy (age and therefore likely provenance). In addition to the nine categories, a further specific category was defined as a “mineralised zone” which is represented by the area 1 km around a mapped surface observation/record of Cornubian metalliferous mineralisation from the 1:50,000 scale British Geological Survey (BGS) digital map (available as a “polyline” Geographical Information System (GIS) file). The 1 km buffer zone provided the GIS “polygon” datafile, which was merged with the mapped solid geological information to provide the ten categorised geological groups. The SBGC was attributed to the centroid of each valid postcode in the datafile provided according to the geological classification.
The simplified classification was used for pre-sampling aggregation of the land area because there is a complex sequence of indurated and poorly permeable rocks. The sedimentary rock formations are variable but mainly argillaceous (fine-grained, clay-rich) types. Igneous rocks have compositions ranging from ultramafic to granitic, and may occupy large continuous areas, or be found as thinner layers intruded within (meta) sedimentary rocks. Mineralisation cross-cuts multiple sedimentary and/or igneous geological boundaries. Physical components can also effect GW flow and chemical composition. Weathering and rock structure (such as faulting and fracturing) can create preferential flow paths, facilitating mobilisation and transport of dissolved mass (even in rock units with otherwise very low bulk hydraulic conductivity values). Mineralisation is often associated with these structural features (upward movement of mineralising fluids are frequently associated with deep-seated fault structures). Structural features may create preferential routes and allow contact between moving water and mineralised ground [40].
These observations meant that a drinking-water sampling strategy for the area should include not only areas of known mineralisation, but also the geological formations represented. Figure 1 shows a map of the sampling area, divided into ten categories on the basis of geological types (basic intrusions are of limited spatial extent and were excluded).

2.4. Recruitment of Households Using PWS

The recruitment strategy aimed to reconstruct a population exposure profile for chemicals in water actually consumed by residents. The existing list of households known to use PWS was provided in late 2010 by Cornwall Council and used for two field campaigns, east Cornwall (2011) and west Cornwall (2013). Households were targeted where it was known that they had a GW source (borehole, spring capture or traditional large-diameter well) serving a PWS. Properties using public supplies, PWS surface water sources, or not using the PWS for domestic DW were excluded. A stratified, random design was used: records were classified into one of the SBGC, within which the records were randomised for sequence of contacting. Householders were encouraged to participate by the provision of the measured chemical data for their property, free-of-charge. Contact was made by letter and follow-up phone calls to ensure adequate recruitment. In east Cornwall, a local newspaper advertisement was also used by the county council. After recruitment for each phase, households that volunteered for the study had their PWS DW sampled in spring (March/April) in 2011 or 2013 and measured for As and a broad suite of 55 elements and other chemical parameters, including anions by ion chromatography, pH/alkalinity and conductivity [30].
From the local authority database, a total of 3095 records were retained that met the inclusion criteria. Thirty-three (1%) records were removed due to the use of surface water for the PWS rather than GW. Some records were incomplete, and these were clarified on contact with potential volunteer householders to ensure that the PWS met the study design criteria.
Contact details of residents were incomplete for 37% of the records. Where a valid address could not be established, the record could not be used (~30% of records). Over 2000 letters were sent inviting households to participate. Users of PWS without an appointment who approached the sampling teams in the field were also included, resulting in some convenience sampling where householders were available for water sampling and interviewing (~2%).
On recruitment, a household representative and/or PWS owner were asked a series of questions by telephone, about the nature of the water supply, water treatment methods, the number of people living in the property and the number of properties the supply served. Further details of the recruitment process are given in [30,31,41]. Further chemical data collection methods were detailed in Ander et al. [30]. Where possible, we collected paired GW and DW samples to show the effect of treatment, but this was possible for only 33% of the households. Using only GW results in the modelling would limit the stratification by SBGC, so we modelled distributions based on DW concentrations.
Where more than one property was using a PWS (shared supplies) the measured value was ascribed to all households which were using this water. For repeat or duplicate samples (<1%), we took the highest (worst case) As concentration, although there was little difference in the levels [30].

2.5. Transformations of Arsenic Distributions

Water chemistry data, however sampled, do not tend to follow a normal distribution. We intended to find a transformation of the As data collected for which the normal distribution approximation of the transformed data would be most appropriate for each SBGC. Traditionally, the following transformations are considered to “normalize” the data: cubic, square, identity, square root, log, inverse square root, inverse, inverse square and inverse cubic. The “ladder” command in Stata was used to explore the best fit transformation. The fit of the transformed data to normal distribution was then tested by Tukey’s Ladder of Powers [42] and the transformation for which the p-value from the test was the largest was selected as being the best transformation, subject to having at least 30 to 50 collected samples for the SBGC being considered. A graphical indication of the suitability of the normal distribution of each transformation was also employed, namely the quantile–quantile (Q–Q) plot, where deviations of the plotted points from a straight line being indicative that a normal distribution might not be suitable. When a variable is restricted to positive values and can vary over a wide range, it is often the situation that the logarithm of the data can be approximated by a normal distribution [43], in which case the data are said to follow a log-normal distribution. Concentrations of chemicals in environmental media often fit a log normal distribution [37,44] and have been used to model radon concentrations in indoor air [39]. This is also the case for As in, for example, DW [45,46], stream water [47], soil [48,49] and biomonitoring samples [50].
We explored the distribution of samples and calculated the geometric mean and standard deviation, examined the potential for fitting a log normal distribution and explored extremes. For SBGC with a small sample size (n < 30) or where there was no suitable transformation, we assumed log normality as that distribution was identified as the best fit for SBGC with a larger number of data points.

2.6. Geological Differences of Arsenic Distribution

To test the need for a model for As concentrations on each geology, we compared the log transformed distributions to assess whether they differed by SBGC. First, we tested the assumption of homoscedasticity (uniformity of variances) by using two variations of Levene’s test [51], both replacing the mean; one with the median and the other with the 10% trimmed mean. If the variances could be considered equal, then a t-test [43] for difference between means was applied, otherwise the Kruskal–Wallis rank sum test [43] was used to compare medians.

2.7. Hazard Assessment Model

Using the log transformed values, represented by the geometric mean and geometric standard deviation, we calculated the cumulative frequency distribution curve in MS Excel for PWS DW As on each SBGC. The probability (and proportion (%)) of dwellings that occurred in each As exposure category of interest was calculated; i.e., low (<1 μg/L), medium (1–5 μg/L), high (5–10 μg/L), and over the PCV (≥10 μg/L).
We compiled a hazard ranking tool in MS Excel to predict the number of dwellings in each exposure category, given the total number of dwellings with PWS on a SBGC. The exposure categories were as follows:
  • More than 10% of the dwellings over the PCV was highest risk-rank 1;
  • 5–10% of dwellings over the PCV-rank 2;
  • 1–5% of dwellings over the PCV-rank 3;
  • Less than 1% of dwellings over the PCV was the lowest risk-rank 4.
Rural land use was found to be a good indicator for the location of PWS, as 97% of PWS in Cornwall were in ‘rural’ areas according to the UK Office of National Statistics (ONS) Urban-Rural (2011) [52] land use classification for Census Output Areas (OAs) [31]. By overlaying the SBGC with the rural land use classifications and residential population layers in the GIS, we could predict the potential number of dwellings with PWS per SBGC. This was fed into the hazard ranking tool.

3. Results

3.1. Study Population and Arsenic Concentrations Measured

Of the Cornish population of 545,335 in 2014 [53], the Drinking Water Inspectorate (DWI) (2015) [1] estimated that 5.25% were served by 3811 PWS, considering all types of PWS, implying that approximately 28,630 people in Cornwall are drinking water from PWS [31].
From the ~2000 invitation letters sent out, around 25% of households took up the offer of water testing; 497 PWS were sampled in 2011–2013 across the SBGC (Figure 1). From the telephone survey, 516 dwellings were recruited and on interview were found to use the 497 PWS sampled. However, eight households did not use the water for drinking purposes. This left 508 domestic dwellings qualifying for inclusion in the study.
Measured As concentrations in DW ranged from 0.2 to 435 μg/L with a median of 0.37 μg/L, 25th percentile of 0.15 μg/L and 75th percentile of 1.46 μg/L (Table 1) [30]. Nearly 6% of sampled households utilised DW from PWS with As concentrations exceeding the PCV of 10 μg/L. A further 4% of households recorded PWS As concentrations between 5 and 10 μg/L, 30% of households had PWS As concentrations between 1–5 μg/L, and 60% were below 1 μg/L [31]. Descriptive statistics of the measured As concentrations per geological classification are shown in Table 1. The number of samples per SBGC ranged from 3 to 93 and on the mineralised zone, n = 140. The highest As concentrations were found on the Lower Carboniferous SBGC (group 4), with 20% of dwellings that were tested on this geology failing the PCV.

3.2. Transformations of Arsenic Distributions

The best fit transformation of the data to normality was identified using Tukey’s ladder of powers based on the largest p value and lowest chi-squared value (Table A1). For SBGC groups 2, 5 and 6 a logarithm was the best transformation. SBGC groups 3, 4, 8, 9 and 10 did not have the minimum number of observations to reliably find the transformation. For groups 1 and 7, normality was not appropriate for any of the transformations. For these and the SBGCs with insufficient observations (n < 30), we assumed a log transformation (Table A1) as this was consistent with other studies of As in environmental samples [45,46,47,48,49,50]. In exploring the fit of log normality to the distribution of As per SBGC, the visual Q–Q plots suggested that many of the groups could be transformed to normality using a log transformation, but variation at the upper bounds of the quartiles sometimes deviated from the lines (Table A1). However, as described in the methods (see Section 2.5) log normality was chosen as best fit, so in examining the distributions, we calculated the median, 25th and 75th percentile values, minimum, maximum, the geometric mean and geometric standard deviation which better described the variation per SBGC (Table 1).

3.3. Geological Differences of Arsenic Distribution

To compare distributions between SBGC using parametric statistics, we tested the assumptions upon which these are based. Using Levene’s robust test statistic for the equality of variances between the groups, we found that there was heteroscedasticity (p < 0.001). Therefore, we applied the Kruskal–Wallis rank sum test, which demonstrated that the As distributions were significantly different (p < 0.001) across SBGC. The medians ranged from 0.08 to 0.5 μg/L on different SBGC.

3.4. Hazard Assessment Model

Applying the log-normal distribution identified above as an appropriate mathematical method for describing the distribution of As in PWS in Cornwall, we modelled As levels per SBGC by using the geometric mean and geometric standard deviation. Using the cumulative frequency distribution, the probability (and proportion) of dwellings falling within the exposure categories of interest was obtained (Table 2). This was converted into the percentage of dwellings that were predicted to be in each exposure category. Using the proportion predicted to exceed the PCV, which ranged from 0 to 20%, we ranked the SBGC according to the risk of high As exposure. SBGC 4 was ranked highest (1), SBGC 1, 5 and 9 ranked second, SBGC 6 third and finally SBGC 2, 3, 7, 8 and 10 were lowest risk (rank 4) (Table 2).
The hazard ranking tool rendered in Microsoft Excel allows the user to input the number of PWS on a SBGC and the number of dwellings predicted to fail the PCV on each SBGC is calculated. For example, if 200 dwellings on each SBGC were identified to have PWS, 96 in total were predicted to fail the PCV for As (5%), with the most occurring on lower Carboniferous and volcanic geology (Group 4) (n = 41, 43% of the total dwellings failing the PCV) (Table S1).

4. Discussion

This paper describes a hazard assessment method to rank hazard according to potential for As concentrations in water by simplified geological classifications. Arsenic exposure varied across SBGCs. Dwellings on the Lower Carboniferous and Volcanics geology (Group 4) were associated with the highest potential for elevated As above the PCV. The method developed can be applied to assess and rank the geological hazards in other areas where As concentrations have been measured on different geologies. The design of the sampling was based on characterising exposure for a population-based risk assessment with public health perspectives.

4.1. Chemical Hazards in PWS

While the acute health effects of human consumption of As-contaminated water are well documented [54,55], the preventable public health burden from acute As exposure is mainly confined to historical incidents in the UK. Extensive delivery of public water supplies to most of the country and robust arrangements for monitoring and controlling DW standards in supplies [1,28,33,56] limits contemporary exposure to As. There is a wealth of evidence and published research on the carcinogenic properties of As and the effects of sustained uptake through this pathway [6,7,57,58,59,60,61].
The hazard assessment showed that As exposure varied across SBGCs with dwellings on the Lower Carboniferous and Volcanics geology (group 4) posing the highest hazard for As concentrations. SBGC alone does not predict GW As hazard; having As rich soil or mine spoil influencing borehole, GW or run off may contribute to elevated As concentrations in water among other factors (e.g., hydrology, well depth, rock structure, etc.) [27,61,62,63,64,65,66,67,68]. For this study, 6% of dwellings with PWS in Cornwall were previously reported by Ander et al. [30] and Crabbe et al. [31] to have elevated As above the PCV. Therefore, posing a risk of health impacts known to derive from As exposure, e.g., developing skin, lung, and/or bladder cancer, peripheral vascular disease, peripheral neurotoxicity and adverse foetal effects [10,69,70], skin disorders, cardiovascular disease, diabetes, and perinatal outcomes [6,7,57,71,72].
Arsenic is not the only chemical hazard potentially present in PWS. In Cornwall, Ander et al. [30] also reported exceedances of PCVs for pH, manganese, nitrate, aluminium, copper, iron, lead, nickel, fluoride, antimony and cadmium. However, in this paper we focused on As, owing to the incidence of exceedance above the PCV and it is the most concerning in terms of possible health impacts [4].

4.2. The Tracking Approach

The goal of an EPHT programme is to identify, monitor and reduce hazards and therefore lead to prevention of chronic disease. This is relevant also to hazards such as contaminants in PWS. Targeted sampling of wells providing PWS is an established practice as part of EPHT programmes in the USA, as it is a more efficient approach to characterisation of population hazard in relation to As [15]. For example, the Florida Department of Health (FDOH) conducts targeted private well testing in areas of concern, in partnership with the Florida Department of Environmental Protection (FDEP), as regulated by Florida Statute. Between 2007 and 2012, the FDOH sampled more than 11,000 wells for As and found approximately 1600 wells with As concentrations above 10 μg/L (this is targeted sampling and does not represent a random sample of wells) [73]. This contrasts with the design that would be usually chosen when characterising SBGC over a large area (e.g., by random sampling) [74].
As both EPHT and geology are complex frameworks involving several disciplines and activities, this study alone could not address all elements of either discipline. However, it is a starting point on which to build other aspects of the EPHT concepts. For example, completing the exposure model, health risk assessment, risk communication and risk management aspects.
One crucial aspect of building the overall EPHT approach to chemicals in PWS included the bedrock as one factor, similar to that applied by others [16,27,67], although in this case a simplified geology approach. Hossain et al. [66] showed the potential for educating local drillers to target safe aquifers on the basis of the color characteristics of sediments to predict risk associated with elevated As in Bangladesh. Ayotte et al. [67] developed a multivariable model for predicting GW As concentrations for exposure assessment. The model included geologic As sources, geochemical processes, hydrogeologic processes and land use as explanatory variables. Some individual bedrock geological units were associated with low As concentrations (e.g., gneiss), while the presence of As concentrations ≥5 μg/L in bedrock wells was associated primarily with local scale geological formations (Eliot, Kittery and Madrid Formations, with up to 13 fold greater odds). Arsenic in stream sediments increased the odds of elevated As concentrations in GW by 67% (p < 0.001). Geochemical processes were also shown to be associated with elevated As, e.g., Pleistocene marine inundation and wells within 3 km of intrusive igneous plutons. We observed similar effects with the “mineralisation zone” (Group 1) with elevated As in PWS (mean of 6.46 μg/L).
Katner et al. [16] also reported that poor water quality data and a large number of unregistered wells impeded risk determination of water As hazards in Louisiana, US. Even with limited data, they were able to rank parishes with a sample of 10 observations or more using the 90th percentile of GW As levels. Parishes were prioritised for monitoring and outreach based on a classification of either “potential hazard” of elevated arsenic and/or “high domestic well water use”. Three parishes within Louisiana were identified as high priority as they had both identified characteristics. Data were evaluated to identify patterns that might suggest additional areas for targeted outreach and monitoring. Through EPHT and water initiatives outreach, barriers to private well water testing were identified. Although not based on geology, the Katner et al. [16] study describes an alternative approach to prioritising areas for PWS risk or hazard assessment.
Overall, the PWS activities focused on chemicals, and in particular As, have provided an important proof of concept study for the development of EPHT in England, leading to increased awareness and acceptance of the EPHT approach.

4.3. Limitations

Based on the consideration that the population risk from As in PWS can be addressed by a dedicated EPHT programme, the study has some limitations. These include the following:
1. The majority of DW samples were subject to treatment (89%), being physical filtration (62%) or chemical treatment (47%) [30]. As participation in the study was voluntary, requiring participants active inclusion, supplies tested are probably more likely to be treated. The As concentrations used to build the model are in part a reflection of treatment adequacy rather than untreated borehole or GW concentrations. Although this is a consequence of the design that could be seen as a limitation from the point of view of a purely geological assessment, comparison of paired GW and DW samples in a subset (33%) indicated that As concentrations were very similar in most cases [30]. We chose such approach of a model built on exposure to As concentrations actually occurring for the population, as this is appropriate for an application to EPHT.
2. The relationship between geology and As in water will likely differ from Cornwall compared to other regions. However, Cornwall was previously reported with elevated As levels in DW, provided a high proportion of rural private water supplies and highly variable geology to test an EPHT approach to hazard assessment in private drinking water supplies according to geological classification [30,73,75,76,77,78]. Although the model may be valid in the area, it could easily be developed and applied to other geographical areas at a later stage.
3. The widespread mining activity in south west England resulted in mine spoil being left on the surface [62,63,64,65] allowing rainwater to carry As into the ground, thus geology is not the sole predictor of As [8,64]. Among other predictors, the local hydrogeology, a range of human activities, seasonality and rainfall may well play a role [67]. These could be considered later, as even in the absence of consideration of such other predictors, the effect of bedrock geology on hazard ranking is strong [68].
4. Geologies containing mineralisation may act as a confounder in the relationship between geology and As concentrations [8]. However, we accounted for this with the Mineralisation category (group 1), allowing for transportation of As in the environment by including a 1 km buffer. Where mineralization is mapped, alternative distance buffers could be explored, such as within 50 m to 500 m. Linear regression could be used to explore relationships between distance from a mineralised vein and transformed As concentrations in PWS.
5. Other than explicitly considering a Mineralisation category, we did not model spatial variation of As hazard within each bedrock geology. In radon mapping, PHE found that radon potential varied within areas shown on the geological maps as a single uniform geological unit [79] and the association of higher arsenic groundwaters with mineralization in this study similarly reflects the potential considerable lateral and vertical intra-unit variation of groundwater chemistry. In this study, we did not collect sufficient samples to allow us to plot the variation of As within each SBGC. Mapping on the basis of uniform geological groups for population level assessment was a first step in the analysis of these data. Further geostatistical analysis may provide insight into such variations, although ultimately the multi-scale complexities of hydrolithological heterogeneities might limit the cost-effective utility of more geologically detailed-based models.
6. Our model was limited by the number of samples in some SBGC and in those circumstances where a unit had fewer than 30 samples we assumed log normality for the reasons mentioned above, where in reality a different distribution may prevail. However, others who examined As distributions in water and soil have found the log-normal distribution to fit with a relatively small sample size (e.g., n = 33 [37], or n = 30 [47,48,79]). Log normality was a good fit for the majority of data points from Q–Q plots of the data. Therefore, the geometric means and geometric standard deviations better represented the central tendency of the distribution than the arithmetic mean and arithmetic standard deviation, so we used the former to model the distributions. Although the (geometric) mean (GM) and standard deviation were calculated from the collected samples for that geology, for those geologies with few samples, the confidence intervals (CI) of those quantities would be very wide (e.g., for group 8 GM 95% CIs are 0.006 to 1.103 μg/L) which in turn would greatly affect the proportions in each of the exposure categories and therefore the ranking. This means that the results for these geologies are less reliable than those with rather more samples.
7. Exposure to As in householders with PWS also depends upon the amount of water consumed, length of residency, seasonality of residency (e.g., permanent versus summer residents) as well as exposures through other routes (food, soil, air, dust) as explored in Cornwall and elsewhere in south west England [34,35,62,76,77]. Urine, hair and toenail samples taken from 127 households of the same cohort of PWS users indicated that As exposure via DW is likely a dominant source [34,35] as found by others [6,36,59,68,73,78,80,81,82].
The above considerations illustrate the limitations in applying this hazard assessment method as a basis for advice to regulators. We estimate likely dwelling numbers and proportions at risk and how to best find those at risk. However, some of the methods described will introduce bias, probably increasing the prediction of the scale of the risk (worst case estimate). Our results are being considered by Cornwall Council for application to their PWS management programme. Also, this paper reports on a methodology that might be applied in other regions, perhaps adapting hazard models according to the local geology. For example, the model could easily be applied to the neighbouring county of Devon, some 4580 km2 that share nine of the same SBGC geologies, covering ~70% of the land area.

4.4. Communication

The EPHT approach supports monitoring, surveillance and evaluation of interventions to minimise environmental exposures and effects on health. Ultimately, for small domestic unregulated PWS on higher risk geologies, owners should be encouraged to test for arsenic, and other chemicals with health based PCVs, in PWS and maintain or put in place effective water treatment systems [4]. Targeted interventions should address social vulnerability factors (e.g., income and education) that influence uptake of testing [17,18,19]. Considerable effort was made in communicating the results of the testing back to PWS users for this study and preparing public health and mitigation advice [31,34]. Fact sheets on the chronic health effects of As and other chemical exposures were devised and a programme of stakeholder and householder communications was essential to support the research [83].
Specific treatment can lower concentration of As and residents are advised to test regularly for As and regularly maintain any water filtration and treatment equipment. PHE produced advice for members of the public about As exposure [83,84]. Cornwall Council now require testing for As in PWS, as it is a known risk in the county. In practice, some, but not all, householders took action to lower their exposure to As in their PWS by maintenance of treatment systems or switching to bottled water [34].

4.5. Implications of This Work and Next Steps

Strengths of this approach allow for the identification and quantification of the populations affected and provides good evidence for allocating and prioritising precious public resources. For example, if a local authority were to test 500 PWS in a Lower Carboniferous and Volcanics geology area, the best estimate from this work is that they would find approximately 100 dwellings having DW over the PCV. If they tested 500 PWS, using this method of prioritisation by SBGC, say in hazard ranking groups 1 and 2, they would capture nearly half (45%) of all of the exceedances occurring in the area. Overlaying the high As risk SBGCs with another predictor of exposure, such as rural land uses (and/or population density), would prioritise areas where further exploration of As risk should occur.
Using a smaller distance unit for exposure related to distance from mineralisation (group 1) may provide a more accurate reflection of the relationship. It seems likely that a combination of distance methods and SBGC would provide more accurate indications of arsenic potential. Further development of the model would also need to consider the relationship with GW samples and geology, as well as additional consideration of other predictors of arsenic risk, e.g., sediment colour, well depth, treatment methods, rock structure, geochemical and hydrogeologic processes and land use as explanatory variables, alongside validation and testing in other areas. A model using geology and mineralogy (e.g., sulfide minerals) of the bedrock is an important first approach for tracking As in groundwater.
The model needs ground-truthing or independent validation by testing in areas not sampled but with similar geologies, such as in Devon, also in south west England. Taking the same approach, the methodology can be applied and expanded to other areas, by assessing the geological classifications and hazard according to those in each region. There is likely to be some commonality with similar geologies. Applying simulation techniques such as bootstrapping or Monte Carlo methods [39] to various distributions can test the model’s robustness. However, a further sampling campaign will be needed to test the estimated distributions. Alternatively, local authorities and regulatory bodies, such as the Drinking Water Inspectorate (DWI), can be encouraged to collate and provide data or co-ordinate data collection in order to test this approach.
This paper describes just one approach for hazard assessment as part of a risk assessment of the health effects of chemicals in PWS [6,85]. Exposure assessment also considers the time and location of exposure, route of entry, study population, magnitude of expected risk, relation between external exposure and biological response, and other contributing and modifying factors (e.g., social economic status, deprivation, occupation) [86]. Nieuwenhuijsen [87] defines a hierarchy of exposure data and surrogates based on the quality and type of available data, although the role of integration of several data sources for one valid exposure model is important. Risk assessments often make use of information that is cost-effective to collect and available within the confines of a study. Surrogate data for characterising exposures often use proxies for “residence or employment in the geographical area in reasonable proximity to the site where exposure can be assumed”. Here we use residence on a particular geology to assume a level of As exposure through PWS. Distance to a mineralised area may also approximate exposure as shown here for patterns of As within 1 km from a mineralised vein or area.
Local authorities may be able to use this method to contribute to the risk assessments required by legislation [28,33]. In the case of a supply provided to a single dwelling, the duty applies only where local authorities are requested to do so by the owner or occupier of a dwelling. The risk assessment must establish whether there is a risk of supplying water that would constitute a potential danger to human health [28]. With further development and testing of validity, the model can calculate the number of dwellings affected by elevated As, help prioritise areas for risk assessments, sampling and can help target areas for public information campaigns. As it stands, the limited sample sizes in some SBGCs limits the robustness of the model, but local authorities and other agencies will be able to test the model to build the rationale for its application.
Local authorities in the UK assess the hazard associated with PWS as required by national regulations [28,33], which was implemented from the EU Directive 98/83/EC [88]. Installation of chemical treatment can decrease As concentrations to levels below the PCV. However, the authority responsible for reviewing the need for risk assessments is the DWI.

5. Conclusions

For this study, the hazard assessment showed that As concentrations varied across SBGCs with dwellings on the Lower Carboniferous and Volcanic geology posing the highest hazard for elevated As above the PCV of 10 μg/L in PWS. Concentrations >10 μg/L increase the risk of developing skin, lung, bladder cancers and several other health impacts known to derive from As exposure.
Within each simplified geological area, concentrations of As in PWS were found or assumed to have log-normal distributions. Across the same areas, the proportion of dwellings predicted to have DW over the prescribed concentration value (PCV) for As ranged from 0% to 20%. From these results, a pilot predictive model was developed calculating the number and percentage of PWS falling within exposure level categories or above the PCV for As. Local authorities may use this to inform PWS risk assessments. The hazard rating will help local authorities to make an informed decision on the need to sample for As.
The geological and GW conditions that promote high arsenic concentrations are known and can help identify high-risk areas for public health [5,8]. Geology varies from surface to depth, therefore, characterising areas by bedrock geology evident in surface mapping also results in oversimplification. However, we found that utilising SBGC to characterise a specific region (Cornwall) for population level exposure groups and exposure to As was helpful to identify significant variation between geologies.
For future work, this paper highlights the potential for low-cost hazard assessment of As in PWS according to SBGC, which can be adjusted to specific regions. There is also the potential to model distributions of other chemicals in PWS and define a similar hazard assessment tool using SGBC along with other predictors, although this will require some adaptation of the model, alongside validation and testing in other areas. A model using mineralogy (e.g., sulfide minerals) of the bedrock would be an important approach for tracking As in groundwater in other areas.
This work can help to support the development of public health risk management, communication and other interventions for authorities in addressing the chemical quality of PWS. Ultimately, users of PWS should be encouraged to test for As and other chemicals if they are at risk, and maintain water treatment to lower consumption of potentially health threatening contaminants. Furthermore, As concentrations can change in any given well over the course of a few years so that regular monitoring is suggested in high-risk areas [8].

Supplementary Materials

The following are available online at www.mdpi.com/1660-4601/14/12/1490/s1, Table S1: Tool for calculating the number of dwellings exceeding the PCV in each exposure category.

Acknowledgments

The authors would like to thank all the householders who participated in the study. PHE colleagues are thanked for their support and input including the South West PHE Centre for field work co-ordination, recruitment of households and project management. Public health students also helped to recruit householders to the study. We thank the Cornwall Council Environmental Protection Team for access to the private water supplies data and for dissemination of results and advice to householders; Helen Taylor and Elliott Hamilton (BGS) for water sampling in the field and laboratory analyses; Stephen Robjohns, PHE Toxicology Group for toxicological and water advice; Zori Daraktchieva and Jon Miles (retired) of the PHE Radon Group for discussions about the application of the radon log-normal probability method are thanked. PHE designed the tracking proof of concept study, funded the study design, field collection and analysis within the Environmental Public Health Tracking programme. All other costs were met by the authors’ organisations. DM acknowledges the receipt of an NERC/CASE (with BGS) Ph.D. studentship (Contract No. GA/125/017, BUFI Ref: S204.2) PHE covered the costs to publish in open access.

Author Contributions

All authors contributed to the writing and reviewing of the paper. Giovanni S. Leonardi and Michael Studden conceived and designed the study; Helen Crabbe, Giovanni S. Leonardi, Rebecca Close and Tony Fletcher conceptualized the aims and objectives of this paper; Helen Crabbe prepared and analysed the data; Helen Crabbe and Giovanni S. Leonardi wrote the paper; Neville Q. Verlander contributed to statistical analysis and advice; Michael J. Watts, E. Louise Ander and Pauline L. Smedley contributed to geological advice and interpretation; Michael J. Watts, E. Louise Ander and Daniel R.S. Middleton collected water samples; Martin Gregory provided access to the original PWS (private water supplies) database and fed back results to PWS users; David A. Polya commented, edited and improved the manuscript ready for publication.

Conflicts of Interest

The authors declare no conflict of interest. PHE as the funding sponsors designed the study; funded the collection, analyses and interpretation of data; PHE funded the writing of the manuscript, and made the decision to publish the results.

Appendix A

Table A1. Summary of arsenic concentration distributions per Simplified Bedrock Geological Classification (SBGC); exploring distributions, log-normality and transformations.
Table A1. Summary of arsenic concentration distributions per Simplified Bedrock Geological Classification (SBGC); exploring distributions, log-normality and transformations.
Log-Transformed DataAssumption Made for Distribution
Geology/Rock Type (SBGC)Number of Properties SampledTransformation for Normal DistributionQ–Q PlotChi Square Test p Value *-
02 Granites and felsic intrusions93log Ijerph 14 01490 i0010.471Log normal
03 Bude and Crackington Formations11Insufficient observations Ijerph 14 01490 i0020.861log normality assumed
04 Lower Carboniferous and Volcanics19Insufficient observations Ijerph 14 01490 i0030.139log normality assumed
05 Middle amd Upper Devonian, with Tamar Gp73log Ijerph 14 01490 i0040.049Log normal
06 Meadfoot and Dartmouth Formations69log Ijerph 14 01490 i0050.256Log normal
07 Gramscatho Formation66none Ijerph 14 01490 i0060.009Log normality assumed
08 Minor Basic and ultrabasic intrusions3Insufficient observations Ijerph 14 01490 i007Insufficient observationslog normality assumed
09 Mylor Slate Formation18Insufficient observations Ijerph 14 01490 i0080.180log normality assumed
10 Lizard and Start Complexes16Insufficient observations Ijerph 14 01490 i0090.0004log normality assumed
01 Mineralised140none Ijerph 14 01490 i0100.017log normality assumed
Non-mineralised (groups 2–10)368none Ijerph 14 01490 i011<0.001n/a
All508none Ijerph 14 01490 i012<0.001n/a
* the null hypothesis is that the data fits the described distribution, so a p value of >0.05 suggests this is a good transformation.

References

  1. Drinking Water Inspectorate (DWI). Drinking Water 2014: Private Water Supplies in England, July 2015. Available online: http://dwi.defra.gov.uk/about/annual-report/2014/pws-eng.pdf (accessed on 7 June 2017).
  2. Bartholomay, R.C.; Carter, J.M.; Qi, S.L.; Squillace, P.J.; Rowe, G.L. Summary of Selected U.S Geological Survey Data on Domestic Well Water Quality for the Centers for Disease Control’s National Environmental Public Health Tracking Program; U.S. Geological Survey Scientific Investigations Report 2007-5213; CDC Stacks: Atlanta, GA, USA, 2007; p. 57.
  3. Flanagan, S.V.; Johnston, R.B.; Zheng, Y. Arsenic in tube well water in Bangladesh: Health and economic impacts and implications for arsenic mitigation. Bull. World Health Organ. 2012, 90, 839–846. [Google Scholar] [CrossRef] [PubMed]
  4. Frisbie, S.H.; Mitchell, E.J.; Mastera, L.J.; Maynard, D.M.; Yusuf, A.Z.; Siddiq, M.Y.; Ortega, R.; Dunn, R.K.; Westerman, D.S.; Bacquart, T.; et al. Public health strategies for western Bangladesh that address arsenic, manganese, uranium, and other toxic elements in drinking water. Environ. Health Perspect. 2009, 117, 410–416. [Google Scholar] [CrossRef] [PubMed]
  5. Polya, D.A.; Middleton, D.R.S. Arsenic in drinking water: Sources & human exposure routes. In Best Practice Guide on the Control of Arsenic in Drinking Water; Bhattacharya, P., Polya, D.A., Jovanovic, D., Eds.; IWA Publishing: London, UK, 2017; Chapter 1; ISBN 13-978-1-84339385-6. [Google Scholar]
  6. Chakraborti, D.; Das, B.; Rahman, M.M.; Nayak, B.; Pal, A.; Sengupta, M.K.; Ahamed, S.; Hossain, M.A.; Chowdhury, U.K.; Biswas, B.K.; et al. Arsenic in groundwater of the Kolkata Municipal Corporation (KMC), India: Critical review and modes of mitigation. Chemosphere 2017, 180, 437–447. [Google Scholar] [CrossRef] [PubMed]
  7. Chakraborti, D.; Rahman, M.M.; Das, B.; Chatterjee, A.; Das, D.; Nayak, B.; Pal, A.; Chowdhury, U.K.; Ahmed, S.; Biswas, B.K.; et al. Groundwater Arsenic Contamination and its Health Effects In India. Hydrogeol. J. 2017, 25, 1165–1181. [Google Scholar] [CrossRef]
  8. Nordstrom, D.K. Worldwide Occurrences of Arsenic in Ground Water. Science 2002, 296, 2143–2145. [Google Scholar] [CrossRef] [PubMed]
  9. NRC (National Research Council). Arsenic in Drinking Water: 2001 Update; National Academy Press: Washington, DC, USA, 2001. [Google Scholar]
  10. World Health Organisation (WHO). Environmental Health Criteria Document 224 Arsenic and Arsenic Compounds. In International Programme on Chemical Safety (IPCS); WHO Regional Office for Europe: Copenhagen, Denmark, 2001. [Google Scholar]
  11. International Agency for Research on Cancer (IARC). IARC Monographs on the Evaluation of Carcinogenic Risks to Humans; IARC: Lyon, France, 2012; Volume 100C, pp. 41–85. [Google Scholar]
  12. Smedley, P.L.; Kinniburgh, D.G. A review of the source, behaviour and distribution of arsenic in natural waters. Appl. Geochem. 2002, 17, 517–568. [Google Scholar] [CrossRef]
  13. Smith, A.H.; University of California at Berkeley, Berkeley, CA, USA. Personal communication, 2017.
  14. Centers for Disease Control and Prevention (US). National Environmental Public Health Tracking. 2017. Available online: http://www.cdc.gov/nceh/tracking/ (accessed on 14 June 2017).
  15. Jordan, M.; Duclos, C.; Kintziger, K.; Gray, A.; Bonometti, M.A. Using an environmental public health tracking biomonitoring study to validate safe water restoration efforts in Florida. J. Public Health Manag. Pract. 2015, 21 (Suppl. 2), S75–S79. [Google Scholar] [CrossRef] [PubMed]
  16. Katner, A.; Lackovic, M.; Streva, K.; Paul, V.; Trachtman, W.C. Evaluation of available data sources to prioritize parishes for arsenic monitoring and outreach related to private well drinking water. J. Public Health Manag. Pract. 2015, 21 (Suppl. 2), S93–S101. [Google Scholar] [CrossRef] [PubMed]
  17. Flanagan, S.V.; Spayd, S.E.; Procopio, N.A.; Chillrud, S.N.; Braman, S.; Zheng, Y. Arsenic in private well water part 1 of 3: Impact of the New Jersey Private Well Testing act on household testing and mitigation behavior. Sci. Total Environ. 2016, 562, 999–1009. [Google Scholar] [CrossRef] [PubMed]
  18. Flanagan, S.V.; Spayd, S.E.; Procopio, N.A.; Chillrud, S.N.; Braman, S.; Zheng, Y. Arsenic in private well water part 2 of 3: Who benefits the most from traditional testing promotion? Sci. Total Environ. 2016, 562, 1010–1018. [Google Scholar] [CrossRef] [PubMed]
  19. Flanagan, S.V.; Spayd, S.E.; Procopio, N.A.; Chillrud, S.N.; Braman, S.; Zheng, Y. Arsenic in private well water part 3 of 3: Socioeconomic vulnerability to exposure in Maine and New Jersey. Sci. Total Environ. 2016, 562, 1019–1030. [Google Scholar] [CrossRef] [PubMed]
  20. Behbod, B.; Lauriola, P.; Leonardi, G.; Crabbe, H.; Close, R.; Staatsen, B.; Knudsen, L.E.; Hoogh, K.; de Medina, S.; Semenza, J.C.; et al. Environmental and public health tracking to advance knowledge for planetary health. Eur. J. Public Health 2016, 26, 900. [Google Scholar] [CrossRef] [PubMed]
  21. Lado, L.R.; Polya, D.; Winkel, L.; Berg, M.; Hegan, A. Modelling arsenic hazard in Cambodia: A geostatistical approach using ancillary data. Appl. Geochem. 2008, 23, 3010–3018. [Google Scholar] [CrossRef]
  22. Winkel, L.; Berg, M.; Amini, M.; Hug, S.J.; Annette Johnson, C. Predicting groundwater arsenic contamination in Southeast Asia from surface parameters. Nat. Geosci. 2008, 1, 536–542. [Google Scholar] [CrossRef]
  23. Amini, M.; Abbaspour, K.M.; Berg, M.; Winkel, L.; Hug, S.F.; Hoehn, E.; Yang, H.; Johnson, A. Statistical Modeling of Global Geogenic Arsenic Contamination in Groundwater. Environ. Sci. Technol. 2008, 42, 3669–3675. [Google Scholar] [CrossRef] [PubMed]
  24. Rodríguez-Lado, L.; Sun, G.; Berg, M.; Zhang, Q.; Xue, H.; Zheng, Q.; Johnson, C.A. Groundwater arsenic contamination throughout China. Science 2013, 341, 866–868. [Google Scholar] [CrossRef] [PubMed]
  25. Sovann, C.; Polya, D.A. Improved groundwater geogenic hazard map for Cambodia. Environ. Chem. 2014, 11, 538–546. [Google Scholar] [CrossRef]
  26. Bretzler, A.; Berg, M.; Winkel, L.; Amini, M.; Rodriguez-Lado, L.; Sovann, C.; Polya, D.A.; Johnson, A. Geostatistical modelling of arsenic hazard in groundwaters. In Best Practice Guide on the Control of Arsenic in Drinking Water; Bhattacharya, P., Polya, D.A., Jovanovic, D., Eds.; IWA Publishing: London, UK, 2017; Chapter A2; ISBN 139-7-81-84339385-6. [Google Scholar]
  27. Ayotte, J.; Montgomery, S.; Flanagan, S.; Robinson, K. Arsenic in Groundwater in Eastern New England: Occurrence, Controls, and Human Health Implications. Environ. Sci. Technol. 2003, 37, 2075–2083. [Google Scholar] [CrossRef] [PubMed]
  28. The Private Water Supplies (England) Regulations 2016. Available online: http://www.legislation.gov.uk/uksi/2016/618/pdfs/uksi_20160618_en.pdf (accessed on 14 June 2017).
  29. Drinking Water Inspectorate (DWI). Drinking Water 2016: Private Water Supplies in England, July 2017. 2017. Available online: http://dwi.defra.gov.uk/about/annual-report/2017/pws-eng.pdf (accessed on 23 October 2017).
  30. Ander, E.L.; Watts, M.J.; Smedley, P.L.; Hamilton, E.M.; Close, R.; Crabbe, H.; Fletcher, T.; Studden, M.; Rimell, A.; Leonardi, G. Variability in private drinking water supply chemistry, and the impact of domestic treatment systems on water quality. Environ. Geochem. Health 2016. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Crabbe, H.; Close, R.; Rimmell, A.; Leonardi, G.; Watts, M.J.; Ander, E.L.; Hamilton, E.M.; Middleton, D.R.S.; Smedley, P.L.; Gregory, M.; et al. Estimating the population exposed to arsenic from groundwater-sourced private drinking water supplies in Cornwall, UK. In Best Practice Guide on the Control of Arsenic in Drinking Water; Bhattacharya, P., Polya, D.A., Jovanovic, D., Eds.; IWA Publishing: London, UK, 2017; Chapter A3; pp. 161–170. ISBN 139-7-81-84339385-6. [Google Scholar]
  32. World Health Organization (WHO). Guidelines for Drinking-Water Quality: Fourth Edition. 2011. Available online: http://whqlibdoc.who.int/publications/2011/9789241548151_eng.pdf?ua=1 (accessed on 14 June 2017).
  33. The Private Water Supplies Regulations 2009. Available online: http://dwi.defra.gov.uk/stakeholders/legislation/pwsregs2009.pdf (accessed on 14 June 2017).
  34. Middleton, D.R.S.; Watts, M.J.; Hamilton, E.M.; Ander, E.L.; Close, R.M.; Exley, K.S.; Crabbe, H.; Leonardi, G.S.; Fletcher, T.; Polya, D.A. Urinary arsenic profiles reveal exposures to inorganic arsenic from private drinking water supplies in Cornwall, UK. Sci. Rep. 2016. [Google Scholar] [CrossRef] [PubMed]
  35. Middleton, D.R.S.; Watts, M.J.; Hamilton, E.M.; Fletcher, T.; Leonardi, G.S.; Close, R.M.; Exley, K.S.; Crabbe, H.; Polya, D.A. Prolonged exposure to arsenic in UK private water supplies: Toenail, hair and drinking water concentrations. Environ. Sci. Process. Impacts 2016, 18, 562–574. [Google Scholar] [CrossRef] [PubMed]
  36. Leonardi, G.S.; Studden, M. Environmental Public Health Tracking Proof of Concept Study—Arsenic in Private Drinking Water Supplies (Hazard Tracking); Internal Report; Health Protection Agency: London, UK, 2010.
  37. Hattis, D.; Burmaster, D.E. Assessment of variability and uncertainty distributions for practical risk analyses. Risk Anal. 1994, 14, 713–730. [Google Scholar] [CrossRef]
  38. Wrixon, A.D.; Green, B.M.R.; Lomas, P.R.; Miles, J.C.H.; Cliff, K.D.; Francis, E.A.; Driscoll, C.M.H.; James, A.C.; O’Riordan, M.C.O. Natural Radiation Exposure in UK Dwellings; NRPB-R190; NRPB (National Radiological Protection Board): Oxfordshire, UK, 1988; p. 188. [Google Scholar]
  39. Daraktchieva, Z.; Miles, J.C.H.; McColl, N. Radon, the lognormal distribution and deviation from it. J. Radiol. Prot. 2014, 34, 183–190. [Google Scholar] [CrossRef] [PubMed]
  40. Domenico, P.A.; Schwartz, F.W. Physical and Chemical Hydrogeology, 2nd ed.; Wiley: New York, NY, USA, 1998; ISBN 978-0-471-59762-9. [Google Scholar]
  41. Close, R.; Studden, M.; Leonardi, G. Private Water Supply Study—Interim Report of Field Study Phase; Internal Report; Health Protection Agency (PHE): London, UK, 2012. [Google Scholar]
  42. Scott, D. Tukey’s Ladder of Powers Description. 2017. Available online: http://onlinestatbook.com/2/transformations/tukey.html (accessed on 16 June 2017).
  43. Armitage, P.; Berry, G. Statistical Methods in Medical Research, 3rd ed.; Blackwell: Oxford, UK, 1994. [Google Scholar]
  44. Ott, W.R. A physical explanation of the lognormality of pollutant concentrations. J. Air Waste Manag. Assoc. 1990, 40, 1378–1383. [Google Scholar] [CrossRef]
  45. Nicolli, H.B.; Suriano, J.M.; Peral, M.A.G.; Ferpozzi, L.H.; Baleani, O.A. Groundwater contamination with arsenic and other trace elements in an area of the pampa, province of Córdoba, Argentina. Environ. Geol. Water Sci. 1989, 14, 3–16. [Google Scholar] [CrossRef]
  46. Shumway, R.H.; Azari, R.S.; Kayhanian, M. Statistical approaches to estimating mean water quality concentrations with detection limits. Environ. Sci. Technol. 2002, 36, 3345–3353. [Google Scholar] [CrossRef] [PubMed]
  47. Lark, R.M.; Bearcock, J.M.; Ander, E.L. How does temporal variation affect the value of stream water as a medium for regional geochemical survey? J. Geochem. Explor. 2016, 169, 211–233. [Google Scholar] [CrossRef]
  48. British Standards Institute (BSI). Soil Quality: Guidance on the Determination of Background Values; BS EN ISO 19258:2011; International Organization for Standardization: Geneva, Switzerland, 2011. [Google Scholar]
  49. Chen, M.; Ma, L.Q.; Hoogeweg, C.G.; Harris, W.G. Arsenic background concentrations in Florida, USA Surface soils: Determination and Interpretation. Environ. Forensics 2001, 2, 117–126. [Google Scholar] [CrossRef]
  50. Samanta, G.; Sharma, R.; Roychowdhuryb, T.; Chakrabortic, D. Arsenic and other elements in hair, nails, and skin-scales of arsenic victims in West Bengal, India. Sci. Total Environ. 2004, 326, 33–47. [Google Scholar] [CrossRef] [PubMed]
  51. Levene, H. Contributions to Probability and Statistics: Essays in Honor of Harold Hotelling; Olkin, I., Ghurye, S.G., Hoeffding, W., Madow, W.G., Mann., H.B., Eds.; Stanford University Press: Redwood City, CA, USA, 1960; pp. 278–292. [Google Scholar]
  52. ONS Urban/Rural Classification. Available online: http://www.ons.gov.uk/ons/guide-method/geography/products/area-classifications/rural-urban-definition-and-la/rural-urban-definition--england-and-wales-/index.html (accessed on 15 June 2017).
  53. Office of National Statistics (ONS). Mid-Year Population Estimates for Cornwall. 2015. Available online: http://www.ons.gov.uk/ons/rel/pop-estimate/population-estimates-for-uk--england-and-wales--scotland-and-northern-ireland/mid-2014/stb---mid-2014-uk-population-estimates.html (accessed on 16 June 2017).
  54. Agency for Toxic Substances and Disease Registry (ATSDR). Toxicological Profile for Arsenic. US Department of Health and Human Services, Public Health Service: Washington, DC, USA, 2007. Available online: https://www.atsdr.cdc.gov/toxprofiles/tp.asp?id=22&tid=3 (accessed on 16 June 2017).
  55. United State Environmental Protection Agency (USEPA). Integrated Risk Information System (IRIS): Arsenic, Inorganic, CASRN 7440-38-2. Available online: https://cfpub.epa.gov/ncea/iris/iris_documents/documents/subst/0278_summary.pdf (accessed on 16 June 2017).
  56. Drinking Water Inspectorate (DWI). Drinking Water 2015: Private Water Supplies in England, July 2015. 2016. Available online: http://dwi.defra.gov.uk/about/annual-report/2015/pws-eng.pdf (accessed on 14 June 2017).
  57. Ayotte, J.D. Arsenic hazard and associated health risks: New England, USA aquifers. In Best Practice Guide on the Control of Arsenic in Drinking Water; Bhattacharya, P., Polya, D.A., Jovanovic, D., Eds.; IWA Publishing: London, UK, 2017; Chapter A1; ISBN 139-7-81-84339385-6. [Google Scholar]
  58. Abernathy, C.O.; Thomas, D.J.; Calderon, R.L. Health effects and risk assessment of arsenic. J. Nutr. 2003, 133, 1536S–1538S. [Google Scholar] [PubMed]
  59. Leonardi, G.; Vahter, M.; Clemens, F.; Goessler, W.; Gurzau, E.; Hemminki, K.; Hough, R.; Koppova, K.; Kumar, R.; Rudnai, P.; et al. Inorganic Arsenic and Basal Cell Carcinoma in Areas of Hungary, Romania, and Slovakia: A Case—Control Study. Environ. Health Perspect. 2012, 120, 721–726. [Google Scholar] [CrossRef] [PubMed]
  60. Smith, A.H.; Hopenhayn-Rich, C.; Bates, M.N.; Goeden, H.M.; Hertz-Picciotto, I.; Duggan, H.M.; Wood, R.; Kosnett, M.J.; Smith, M.T. Cancer risks from arsenic in drinking water. Environ. Health Perspect. 1992, 97, 259–267. [Google Scholar] [CrossRef] [PubMed]
  61. Smith, A.H.; Hopenhayn-Rich, C.; Bates, M.N.; Goeden, H.M.; Hertz-Picciotto, I.; Duggan, H.M. Health risk assessment of inorganic arsenic intake of Cambodia residents through groundwater drinking pathway. Water Res. 2010, 44, 5777–5788. [Google Scholar]
  62. Middleton, D.R.S.; Watts, M.J.; Beriro, D.J.; Hamilton, E.M.; Leonardi, G.S.; Fletcher, T.; Close, R.; Polya, D.A. Arsenic in residential soil and household dust in Cornwall, South West England: Potential human exposure and the influence of historical mining. Environ. Sci. Process. Impacts 2017, 19, 517–527. [Google Scholar] [CrossRef] [PubMed]
  63. Klinck, B.; Palumbo-Roe, B.; Cave, M.; Wragg, J. Arsenic Dispersal and Bioaccessibility in Mine Contaminated Soils: A Case Study from an Abandoned Arsenic Mine in Devon, UK, British Geological Survey (BGS) Research Report RR/04/003. Available online: http://nora.nerc.ac.uk/3681/1/RR04003.pdf (accessed on 16 June 2017).
  64. Palumbo-Roe, B.; Klinck, B. Bioaccessibility of arsenic in mine waste-contaminated soils: A case study from an abandoned arsenic mine in SW England (UK). J. Environ. Sci. Health 2007, 42, 1251–1261. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Camm, G.S.; Glass, H.J.; Bryce, D.W.; Butcher, A.R. Characterisation of a mining-related arsenic-contaminated site, Cornwall, UK. J. Geochem. Explor. 2004, 82, 1–15. [Google Scholar] [CrossRef]
  66. Hossain, M.; Bhattacharya, P.; Frape, S.K.; Jacks, G.; Islam, M.M.; Rahman, M.M.; von Brömssenae, M.; Hasand, M.A.; Ahmed, K.M. Sediment color tool for targeting arsenic-safe aquifers for the installation of shallow drinking water tubewells. Sci. Total Environ. 2014, 493, 615–625. [Google Scholar] [CrossRef] [PubMed]
  67. Ayotte, J.D.; Nolan, B.T.; Nuckols, J.R.; Canton, K.P.; Robinson, G.R.; Baris, D.; Hayes, L.; Karagas, M.; Bress, W.; Silverman, D.T.; et al. Modeling the probability of arsenic in groundwater in New England as a tool for exposure assessment. Environ. Sci. Technol. 2006, 40, 3578–3585. [Google Scholar] [CrossRef] [PubMed]
  68. Bondu, R.; Cloutier, V.; Rosa, E.; Benzaazoua, M. Mobility and speciation of geogenic arsenic in bedrock groundwater from the Canadian Shield in western Quebec, Canada. Sci. Total Environ. 2017, 574, 509–519. [Google Scholar] [CrossRef] [PubMed]
  69. World Health Organisation (WHO). Safety Evaluation of Certain Contaminants in Food. In Food Additives 2011: Geneva and Rome; Series: 63, FAO JECFA Monographs 8; WHO: Geneva, Switzerland, 2011. [Google Scholar]
  70. UK Teratology Information Service (UKTIS). Exposure to Arsenic in Pregnancy, (Accessed via TOXBASE for the UK) v1.1. 2012. Available online: http://www.medicinesinpregnancy.org/bumps/monographs/EXPOSURE-TO-ARSENIC-IN-PREGNANCY/ (accessed on 16 June 2017).
  71. Smith, A.H.; Steinmaus, C.M. Health effects of arsenic and Chromium in drinking Water: Recent human findings. Annu. Rev. Public Health 2009, 30, 107–122. [Google Scholar] [CrossRef] [PubMed]
  72. Hopenhayn, C. Arsenic in Drinking Water: Impact on human health. Elements 2006, 2, 103–107. [Google Scholar] [CrossRef]
  73. Johnson, L.R.; Farmer, J.G. Urinary arsenic concentrations and speciation in Cornwall residents. Environ. Geochem. Health 1989, 11, 39–44. [Google Scholar] [CrossRef] [PubMed]
  74. Florida Statutes 2016, Chapter 373, Part III: Water Resources, Regulation of Wells. Available online: http://law.justia.com/codes/florida/2016/title-xxviii/chapter-373/part-iii/ (accessed on 17 February 2017).
  75. Ander, L.; Lark, M.; Smedley, P.; Watts, M.; Hamilton, E.; Fletcher, T.; Crabbe, H.; Close, R.; Studden, M.; Leonardi, G. Testing how voluntary participation requirements in an environmental study affect the planned random sample design outcomes: Implications for the predictions of values and their uncertainty. EGU Gen. Assem. Conf. Abstr. 2015, 17, 15629. [Google Scholar]
  76. Smedley, P.L.; Allen, D. Baseline Report Series: 16, the Granites of South-West England; BGS Commissioned Report CR/04/255; British Geological Survey, NERC 2004. Available online: http://nora.nerc.ac.uk/3547/1/CR04255N.pdf (accessed on 16 June 2017).
  77. Xu, J.; Thornton, I. Arsenic in garden soils and vegetable crops in Cornwall, England: Implications for human health. Environ. Geochem. Health 1985, 7, 131–133. [Google Scholar] [CrossRef] [PubMed]
  78. Farago, M.E.; Thornton, I.; Kavanagh, P.; Elliott, P.; Leonardi, G.S. Health aspects of human exposure to high arsenic concentrations in soil in south-west England. In Arsenic: Exposure and Health Effects; Abernathy, C.O., Calderon, R.L., Chappell, W.R., Eds.; Springer: Dordrecht, The Netherlands, 1997. [Google Scholar]
  79. Miles, J.C.H.; Appleton, J.D. Mapping variation in radon potential both between and within geological units. J. Radiol. Prot. 2005, 25, 257–276. [Google Scholar] [CrossRef] [PubMed]
  80. Kavanagh, P.; Farago, M.E.; Thornton, I.; Goessler, W.; Kuehnelt, D.; Schlagenhaufen, C.; Irgolic, K.J. Urinary arsenic species in Devon and Cornwall residents, UK. A pilot study. Analyst 1998, 123, 27–29. [Google Scholar] [CrossRef] [PubMed]
  81. Button, M.; Jenkin, G.R.T.; Harrington, C.F.; Watts, M.J. Human toenails as biomarkers of exposure to elevated environmental arsenic. J. Environ. Monit. 2009, 11, 610–617. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Mondal, D.; Banerjee, M.; Kundu, M.; Banerjee, N.; Bhattacharya, U.; Giri, A.K.; Ganguli, B.; Sen Roy, S.; Polya, D.A. Comparison of drinking water, raw rice and cooking of rice as arsenic exposure routes in three contrasting areas of West Bengal, India. Environ. Geochem. Health 2010, 32, 463–477. [Google Scholar] [CrossRef] [PubMed]
  83. British Geological Survey (BGS). Arsenic: Are We Exposed? 2017. Available online: http://www.bgs.ac.uk/sciencefacilities/laboratories/geochemistry/igf/Biomonitoring/arsenicSW.html (accessed on 12 April 2017).
  84. PHE (Public Health England). Information for the Public on the Effects of Arsenic. 2016. Available online: https://www.gov.uk/government/uploads/system/uploads/attachment_data/file/562431/arsenic_general_information.pdf (accessed on 16 June 2017).
  85. Yassi, A.; Kjellstrom, T.; de Kok, T.; Guidotti, T. Risk Assessment Chapter in ‘Basic Environmental Health’, Chapter 3; Oxford University Press: Oxford, UK, 2001; pp. 104–142. [Google Scholar]
  86. Zeka, A. Exposure assessment for epidemiology. In Essentials of Environmental Epidemiology for Health Protection: A Handbook for Field Professionals; Kreis, I.A., Busby, A., Leonardi, G.L., Meara, J., Murray, V., Eds.; Oxford University Press: Oxford, UK, 2013; Chapter 8; pp. 68–76. [Google Scholar]
  87. Nieuwenhuijsen, M.J. Exposure Assessment in Occupational and Environmental Epidemiology; OUP: Oxford, UK, 2003. [Google Scholar]
  88. Council of the European Union. The Quality of Water Intended for Human Consumption. Council Directive 98/83/EC. 1998. Available online: http://ec.europa.eu/environment/water/water-drink/legislation_en.html (accessed on 16 June 2017).
Figure 1. Map of Cornwall showing Simplified Bedrock Geological Classifications (SBGC), location and levels of arsenic measured in Private Water Supplies (PWS) tested. Geological data reproduced with the permission of the British Geological Survey © Natural Environment Research Council (NERC).
Figure 1. Map of Cornwall showing Simplified Bedrock Geological Classifications (SBGC), location and levels of arsenic measured in Private Water Supplies (PWS) tested. Geological data reproduced with the permission of the British Geological Survey © Natural Environment Research Council (NERC).
Ijerph 14 01490 g001
Table 1. Descriptive statistics of arsenic concentrations (μg/L) in Private Water Supplies (PWS) in dwellings classified by Simplified Bedrock Geological Category (SBGC) measured in Cornwall, UK.
Table 1. Descriptive statistics of arsenic concentrations (μg/L) in Private Water Supplies (PWS) in dwellings classified by Simplified Bedrock Geological Category (SBGC) measured in Cornwall, UK.
Geology/Rock Type (SBGC)Number of Dwellings SampledPercentage of Total DwellingsArithmetic Mean AsMedian As *Minimum As $25% Tile75% TileMaximum AsGeometric Mean As Geometric Standard Deviation As% of Dwellings Measured at Levels: (As μg/L)
<11–55–10>10
01 Mineralised140286.460.710.020.222.672310.8851.805582769
02 Granites and felsic intrusions93181.240.50.020.261.1320.50.5751.214712711
03 Bude and Crackington Formations1120.120.090.020.040.170.370.0860.885100000
04 Lower Carboniferous and Volcanics19433.130.40.050.254.374351.2132.5425821021
05 Middle and Upper Devonian, with Tamar Gp73145.080.420.020.181.971780.6221.8196815115
06 Meadfoot and Dartmouth Formations69141.510.380.020.121.1914.20.4411.577682534
07 Gramscatho Formation66130.400.170.020.090.325.130.1751.10994330
08 Minor Basic and ultrabasic intrusions3<10.120.080.030.030.240.240.0831.040100000
09 Mylor Slate Formation1844.640.380.060.153.6125.60.6652.024726617
10 Lizard and Start Complexes1630.330.130.040.120.153.290.1410.97794600
Non-mineralised (groups 2–10)368723.630.310.020.130.924350.4111.637761744
All5081004.410.370.020.151.464350.5081.718711946
* Kruskal–Wallis test result comparing medians p < 0.001. $ Limit of detection is 0.02 μg/L.
Table 2. Data transformations; summary statistics of logged transformed values, probability of dwellings being in an exposure category and hazard ranking of Simplified Bedrock Geological Classifications (SBGC).
Table 2. Data transformations; summary statistics of logged transformed values, probability of dwellings being in an exposure category and hazard ranking of Simplified Bedrock Geological Classifications (SBGC).
Log Transformed Data% of Dwellings over the Geological Groups Predicted at Levels: (As μg/L)Hazard Ranking
Geology/Rock Type (SBGC)Number of Dwellings SampledGeometric MeanGeometric Standard Deviation<11–55–10>10
01 Mineralised1400.8851.80552.730.47.918.952
02 Granites and felsic intrusions930.5751.21467.628.72.810.934
03 Bude and Crackington Formations110.0860.88599.70.28004
04 Lower Carboniferous and Volcanics191.2132.54247.024.28.5420.31
05 Middle and Upper Devonian, with Tamar Gp730.6221.81960.327.16.256.342
06 Meadfoot and Dartmouth Formations690.4411.57769.824.03.792.383
07 Gramscatho Formation660.1751.10994.25.670.110.014
08 Minor Basic and ultrabasic intrusions30.0831.04099.20.730.070.044
09 Mylor Slate Formation180.6652.02458.026.16.929.032
10 Lizard and Start Complexes160.1410.97797.82.230.0104

Share and Cite

MDPI and ACS Style

Crabbe, H.; Fletcher, T.; Close, R.; Watts, M.J.; Ander, E.L.; Smedley, P.L.; Verlander, N.Q.; Gregory, M.; Middleton, D.R.S.; Polya, D.A.; et al. Hazard Ranking Method for Populations Exposed to Arsenic in Private Water Supplies: Relation to Bedrock Geology. Int. J. Environ. Res. Public Health 2017, 14, 1490. https://doi.org/10.3390/ijerph14121490

AMA Style

Crabbe H, Fletcher T, Close R, Watts MJ, Ander EL, Smedley PL, Verlander NQ, Gregory M, Middleton DRS, Polya DA, et al. Hazard Ranking Method for Populations Exposed to Arsenic in Private Water Supplies: Relation to Bedrock Geology. International Journal of Environmental Research and Public Health. 2017; 14(12):1490. https://doi.org/10.3390/ijerph14121490

Chicago/Turabian Style

Crabbe, Helen, Tony Fletcher, Rebecca Close, Michael J. Watts, E. Louise Ander, Pauline L. Smedley, Neville Q. Verlander, Martin Gregory, Daniel R. S. Middleton, David A. Polya, and et al. 2017. "Hazard Ranking Method for Populations Exposed to Arsenic in Private Water Supplies: Relation to Bedrock Geology" International Journal of Environmental Research and Public Health 14, no. 12: 1490. https://doi.org/10.3390/ijerph14121490

APA Style

Crabbe, H., Fletcher, T., Close, R., Watts, M. J., Ander, E. L., Smedley, P. L., Verlander, N. Q., Gregory, M., Middleton, D. R. S., Polya, D. A., Studden, M., & Leonardi, G. S. (2017). Hazard Ranking Method for Populations Exposed to Arsenic in Private Water Supplies: Relation to Bedrock Geology. International Journal of Environmental Research and Public Health, 14(12), 1490. https://doi.org/10.3390/ijerph14121490

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop