Article Text
Abstract
OBJECTIVES To explore the use of routinely collected trihalomethane (THM) measurements for epidemiological studies. Recently there has been interest in the relation between byproducts of disinfection of public drinking water and certain adverse reproductive outcomes, including stillbirth, congenital malformations, and low birth weight.
METHOD Five years of THM readings (1992–6), collected for compliance with statutory limits, were analysed. One water company in the north west of England, divided into 288 water zones, provided 15 984 observations for statistical analysis. On average each zone was sampled 11.1 times a year. Five year, annual, monthly, and seasonal variation in THMs were examined as well as the variability within and between zones.
RESULTS Between 1992 and 1996 the total THM (TTHM) annual zone means were less than half the statutory concentration, at approximately 46 μg/l. Differences in annual water zone means were within 7%. Over the study period, the maximum water zone mean fell from 142.2 to 88.1 μg/l. Mean annual concentrations for individual THMs (μg/l) were 36.6, 8.0, and 2.8 for chloroform, bromodichloromethane (BDCM), and dibromochloromethane (DBCM) respectively. Bromoform data were not analysed, because a high proportion of the data were below the detection limit. The correlation between chloroform and TTHM was 0.98, between BDCM and TTHM 0.62, and between DBCM and TTHM −0.09. Between zone variation was larger than within zone variation for chloroform and BDCM, but not for DBCM. There was only little seasonal variation (<3%). Monthly variation was found although there were no consistent trends within years.
CONCLUSION In an area where the TTHM concentrations were less than half the statutory limit (48 μg/l) chloroform formed a high proportion of TTHM. The results of the correlation analysis suggest that TTHM concentrations provided a good indication of chloroform concentrations, a reasonable indication of BDCM concentrations, but no indication of DBCM. Zone means were similar over the years, but the maximum concentrations reduced considerably, which suggests that successful improvements in treatment have been made to reduce high TTHM concentrations in the area. For chloroform and BDCM, the main THMs, the component between water zones was greater than variation within water zones and explained most of the overall exposure variation. Variation between months and seasons was low and showed no clear trends within years. The results indicate that routinely collected data can be used to obtain exposure estimates for epidemiological studies at a small area level.
- chlorination byproducts
- exposure
- birth outcomes
- routine data
- trihalomethanes
Statistics from Altmetric.com
There is growing concern over the effects of disinfection byproducts of chlorination on human health. Early studies investigated links between disinfection byproducts of chlorination and increased risk of certain cancers, such as rectal, bladder, and colon cancer.1 More recent studies have focused on the possible effects of disinfection byproducts of chlorination on certain reproductive outcomes, including low birthweight, stillbirth, spontaneous abortion, and congenital abnormalities. The results from these studies have been mixed.2 3
Several of these reproductive epidemiological studies relied on crude exposure data. Some used comparisons between water source (groundv surface) and between water treatments (chlorinated vnon-chlorinated).4 5 Others have used routinely collected total trihalomethane (TTHM) measurements to form exposure categories.6-9 More recent studies have attempted to measure individual exposure with routinely collected THM measurements combined with a measure of ingestion such as number of glasses or water drunk per day.10-13 Misclassification of exposure as a result of crude exposure assessment may lead to biased risk estimates, which makes it more difficult to interpret the studies.14
Trihalomethanes (THMs), and many other disinfection byproducts of chlorination, are formed when chlorine reacts with naturally occurring organic substances present in water. They are generally the most common byproducts. It is thought that THM concentrations vary both temporally and spatially.15 16 In the water company under study much of the supply is upland surface water, which has higher concentrations of THM precursors than water drawn predominantly from aquifers. Surface water sources are more likely to produce chlorinated THMs, whereas treated water from ground water sources can produce higher concentrations of brominated THMs.
Water companies in the United Kingdom divide the area under their control into water supply zones, hydraulically defined in such a way that the quality of drinking water on entry to a zone would be the same for all consumers within that zone. These zones form the basis for monitoring the water quality of drinking water so they define populations of consumers to which routinely collected water data can be related. The water in each water zone is tested a minimum of four times a year to comply with the TTHM regulatory standard. Of all the disinfection byproducts, only THM concentrations are currently regulated, although they are not the only class of disinfection byproduct suggested as having adverse reproductive effects. Haloacetic acids, haloacetonitriles, and acetonitriles have also been implicated.3 17
To assess the effect of THM exposure concentrations on certain reproductive outcomes, routinely collected quarterly (seasonal) THM concentrations could be attached to a pregnancy trimester of women living in the water zone, but such exposure information would generally be based on only one measurement, and given the variation in exposure, that might be unreliable. Annual averages would provide more precise THM means, but would be based on perhaps four measurements and might not reflect exposure for any 3 month period within that year. A 5 year average THM concentration would be based on more measurements (between 20 and 60), but might be insensitive to timing of exposures relevant to birth outcomes.
The purpose of this study was to explore routinely collected THM data for use in epidemiological studies and to determine the best measure of THM exposure. The specific aims of this study were to examine; (a) the current exposure concentrations of TTHM and individual THMs in tap water in the north west area, (b) the correlation between TTHMs and individual THM concentrations, (c) the variability of total and individual THM concentrations between and within zones, seasons, months, and years. This will provide information to enable choice of the most appropriate routinely available exposure estimate for epidemiological studies, in particular reproductive health studies.
Methods
THE DATA
Routinely collected data on TTHM concentrations were obtained from a United Kingdom water company in the north of England for the period 1992–6. Information was also obtained on the geography of the supply zones, the type of water each receives and the population each supplies.
DATA COLLECTION
Four THMs are routinely measured by United Kingdom water companies. They are chloroform, bromodichloromethane (BDCM), dibromochloromethane (DBCM), and bromoform. Their concentrations in public drinking water are measured in samples from randomly selected taps in each water zone. The sum of the concentrations of the four individual THMs is the TTHM concentration. Under the Water Supply (Water Quality) (1989) Regulations, the water companies are required to collect four samples/year/water zone. The statutory limit for TTHMs is 100 μg/l as a rolling 3 month average. Exceedance of the limit requires monthly measurements to be taken until the TTHM concentrations decrease. In many zones monthly measurements are taken routinely for operational purposes.
WATER SUPPLY (OR COMPLIANCE) ZONES
Water zones differ in area, in the number of people each serves, and in the volume of water passing through each per day. Small, rural zones may serve 300 people 70 kl per day, whereas others may deliver as much as 16 000 kl per day to 50 000 people. The main water sources in the area are Thirlmere and Ullswater in the Lake District. Only 6% of zones are supplied exclusively by groundwater. The number of water zones may change. In north west England during 1992 there were 313 zones, in 1993 316, in 1994 317, in 1995 and 1996 323. Some of the new zones were created whole, whereas some were split from existing zones. It was decided to omit the data from all zones created after 1992, from all zones created from existing zones, and from the zones from which these were split. Only zones in which boundaries and source did not change over the study period and for which 5 years' worth of data were available were included. For the statistical analysis, this left 288 water zones supplied by 194 treatment works, serving about 6.8 million people.
LABORATORY ANALYSIS
Water samples were analysed by the water company with a capillary gas chromatograph with an electron captive detector. The detection limit was 2.5 μg/l for all THMs. The water company runs internal and external quality control tests and is accredited by the United Kingdom Accreditation Service (UKAS). The quality of the THM data is audited by the United Kingdom Drinking Water Inspectoarate.
STATISTICAL ANALYSIS
Descriptive analyses were carried out to calculate over time the zone means for total and individual THMs plus correlation coefficients between individual and total THMs. The distribution of theTTHMs could be best described with a normal distribution. However, the distributions of the individual THMs were positively skewed. No statistical analysis of the bromoform data was attempted as most of the readings (95%) were below the detection limit.
To find within and between zone variation the individual THM data were first transformed with Box-Cox transformations.18 This provided approximately normally distributed values. These were then analysed with a multilevel (variance components) model to estimate the mean and variance within and between zones for each THM on the transformed scale.19 This modelling was necessary to provide more reliable estimates of the variation within the zones with few data.
It is more meaningful to report estimates of the variation on the scale of the original data rather than on the scale of the transformed data. However, the variance may be a misleading summary on the original scale as the distributions were skewed. An alternative measure of the spread of the distributions within and between zones can be provided by the interquartile range (IQR). This is simply the difference between the upper and lower quartiles, giving the range of the middle half of the data, and makes no assumptions about the distribution of the data. Interquartile ranges within and between zones were therefore calculated for each transformed THM using the 25th and 75th percentiles of the cumulative normal distribution based on the appropriate estimates of mean and variance from the multilevel model. These IQRs were then back transformed onto the original scale of the data.
Results
TTHM CONCENTRATIONS
The range of TTHM concentrations across all water zones was 1.65–188.7 μg/l. The mean number of observations/water zone over 5 years was 55.5, or 11.1 a year, although this ranged from a minimum of four to a maximum of 45 a year.
Over the 5 years the mean annual zone TTHM concentrations were very similar, and all were less than half the statutory limit of 100 μg/l (table 1). The maximum zone TTHM mean fell by 38% over 5 years. Figure1 shows the distribution of TTHM zone means.
Means of the zone TTHM concentrations (μg/l) in the north west of England, 1992–6
Frequency distribution of 5 year total trihalomethane zone means, north west of England 1992–6.
INDIVIDUAL THM CONCENTRATIONS
Table 2 shows the zone means for individual THMs between 1992 and 1996. For chloroform, the means were similar over the years. For both BDCM and DBCM they rose slightly. The maximum recorded concentrations showed an increasing trend over the years. Most (95%) bromoform measurements were below the detection limit and no analysis of its concentrations was undertaken.
Characteristics of individual THM concentrations (μg/l) in the north west of England
CORRELATION BETWEEN YEARS
The mean annual zone TTHM concentration was plotted against the 5 year mean for each zone (fig 2) to assess the variation over the years. Overall the correlation between years was high (r⩾0.88). The correlation between the mean annual zone concentrations for the individual THMs was also high (chloroform r=0.88, BDCMr=0.87, DBCM r= 0.83).
Five year total trihalomethane zone means plotted against annual total trihalomethane zone means.
CORRELATION BETWEEN INDIVIDUAL THM AND TTHM
Figure 3 shows the correlations between annual zone means of the individual THM concentrations and annual zone means of the TTHM concentrations for 1995. Chloroform correlated well with the TTHM (r=0.98), but BDCM correlated less well with TTHMs (r=0.62). There was no correlation between TTHM and DBCM.
Zone means 1995. Scatter plot of individual THM zone means against total THM zone means.
VARIATION OF INDIVIDUAL THMS WITHIN AND BETWEEN ZONES
Table 3 shows the IQR between zones and a summary of the IQR for each distribution within a zone. For chloroform and BDCM the spread within zones is less than the spread between zones, particularly in later years. For DBCM, the mean of the IQRs within zones is greater than the IQR between zones for each year, but the medians of the IQRs within zones are lower. The IQRs within zones for chloroform and BDCM decrease throughout the 5 years, conversely the IQRs within zones for DBCM increase throughout the study period. For all three individual THMs, the ratio of the variation within zones to between zones increased throughout the study period.
Estimated interquartile ranges (IQRs) within and between zones for individual THMs (μg/l) in the north west of England
Discussion
In this study we found that concentrations of TTHM in the north west of England were on average just under half the statutory limit of 100 μg/l between 1992 and 1996. Water zone means remained at a similar level over the 5 years. Concentrations of TTHMs correlated well with chloroform, the predominant THM, less well with BDCM, but not with DBCM. The variation between zones was generally larger than the variation within zones for chloroform and BDCM concentrations, in particular in later years, but not for DBCM concentrations. Models incorporating year, season, or month showed that none of these explained more than 3% of the total variation in TTHMs.
Water zone mean concentrations of TTHM were very similar over the years, but the maximum zone means fell from 142.2 to 88.0 μg/l (32% reduction). This suggests that measures to reduce TTHM concentrations were well targeted towards the zones with the highest TTHM and have been successful in reducing high mean concentrations. Over the study period the variation within zones fell for chloroform and BDCM, and the variation between zones increased relative to the variation within zones. The fall in variability within each zone suggests that treatment methods became more efficient—for example, either at removing the THM precursor or at reducing THM formation in the network. The water company in this study introduced new technology into several treatment plants between 1992 and 1996 and this could have had some direct effect on stabilising the THM concentrations.
High variation between zones compared with variation within zones is desirable for geographical epidemiological studies. It is more likely that zone means are distinct from each other, and therefore there is less overlap in their distributions. In this case people living in different water zones would be exposed to distinctly different THM concentrations, which reduces the risk of exposure misclassification and bias in health risk estimates. Furthermore, when the variation within zones is small, zone means are likely to be more precise, even with a few measurements.
Concentrations of TTHM varied only little by season suggesting that the annual mean may be a good estimate for various seasons. This contrasts with some other studies that have shown that TTHM concentrations vary considerably throughout the year.15 16 20 We found that TTHM concentrations were generally lowest in the spring and highest at the end of summer and during the autumn, but there was a lack of consistency. Spring concentrations in some years were higher than autumn concentrations.21 An annual cycle would be consistent with the factors that affect the rate and quantity of THM formation; during the summer and autumn water temperatures are higher than in the winter and spring. Temperature is known to increase the rate of THM formation.15 Formation of THMs is also dependent on the quantity of organic precursors, and these are more abundant during the late summer and autumn, when there is increased leaching of organic matter from decaying vegetative matter. Improvements in treatment will also help to reduce the variation and this may explain our findings.
Our results showed that the correlation between yearly zone means and 5 year zone means was high, which suggests that there are areas where the THM concentrations are consistently high and consistently low in the region. This is a desirable feature for epidemiological studies as it reduces the potential for exposure misclassification.
In our study area the dominant THM was chloroform. In surface waters chloroform is often the most easily formed THM, followed by DCBM and DBCM. Bromoform is the least easily formed.15 This pattern of formation is reflected in their mean concentrations for 1995 and 1996, with chloroform existing at concentrations five times that of DCBM, 12 times that of DBCM, and 20 times that of bromoform. These ratios are consistent with those commonly found in chlorinated drinking water predominantly drawn from upland surface water sources. The high concentration of chloroform is consistent with the continued reaction of chlorine with organic matter and its production in preference to other THMs because of the low bromide content surface water.20
Concentrations of TTHM were correlated most strongly with chloroform, then BDCM. There was no correlation with DBCM. These findings suggest that TTHM concentrations most accurately reflect chloroform concentrations, particularly in areas where surface waters represent the dominant water source. However, given that TTHM is a good indicator only of chloroform and not of the other individual THM concentrations, data on other individual THM concentrations, and possibly other disinfection byproducts of chlorination, would need to be incorporated into epidemiological studies to fully explore any potential health effects of those byproducts. This is important as different THMs may be associated with different health outcomes.3 22 For example BDCM has been more strongly associated with spontaneous abortion than have other THMs.14
Over the study period the variations in chloroform and BDCM between zones were greater than the variations within zones. However, it should be pointed out that even in the year when variation between zones was at its highest (1996), a considerable proportion of the total variation was within zones. This variation within zones could be a result of real variability in THMs caused by changes in—for example, precursor concentrations, chlorine dose, residence time in the distribution system, water source, temperature, and pH.20 Limited sampling in a zone with high variability in THMs could lead to imprecise estimates and therefore exposure misclassification when water zone means are linked to an individual person within the zone. More information on variation within zones should be sought. Further statistical modelling with smoothing techniques and inclusion of information on—for example, residence time, pH, and temperature—could then be carried out to increase the precision of the water zone means.
The evidence presented here suggests that in this study area an annual water zone mean calculated from routinely collected health data may be a reasonable measure of exposure for a reproductive epidemiological study of those people resident within the zone. However, uncertainties remain. Concentrations of TTHMs were only strongly correlated with chloroform concentrations, and water zone means for other disinfection byproducts of chlorination should be used in the epidemiological study, if possible, to provide more accurate estimates. The situation in other areas of the country where the water supply characteristics are different may not be the same. Also these measures do not take into account personal exposure characteristics. For example, people can move between zones during pregnancy and there will be variation among pregnant women in patterns of water consumption, as well as showering, bathing, and swimming habits.23 However, it is possible that these factors are fairly similar on a population level by water zone and that the variation in water zone mean may be the main determinant of uptake of disinfection byproducts of chlorination in the population under study. Further work needs to be carried out to either confirm or refute this.
Acknowledgments
The Small Area Health Statistics Unit is funded by a grant from the Department of Health, Department of the Environment, Transport and the Regions, Health and Safety Executive, Scottish Office Home and Heath Department, Welsh Office, and Northern Ireland Department of Health and Social Services. The views expressed in this publication are those of the authors and not necessarily of the funding departments.