Comparison of pharmaceutical, illicit drug, alcohol, nicotine and caffeine levels in wastewater with sale, seizure and consumption data for 8 European cities
BMC Public Health volume 16, Article number: 1035 (2016)
Monitoring the scale of pharmaceuticals, illicit and licit drugs consumption is important to assess the needs of law enforcement and public health, and provides more information about the different trends within different countries. Community drug use patterns are usually described by national surveys, sales and seizure data. Wastewater-based epidemiology (WBE) has been shown to be a reliable approach complementing such surveys.
This study aims to compare and correlate the consumption estimates of pharmaceuticals, illicit drugs, alcohol, nicotine and caffeine from wastewater analysis and other sources of information. Wastewater samples were collected in 2015 from 8 different European cities over a one week period, representing a population of approximately 5 million people. Published pharmaceutical sale, illicit drug seizure and alcohol, tobacco and caffeine use data were used for the comparison.
High agreement was found between wastewater and other data sources for pharmaceuticals and cocaine, whereas amphetamines, alcohol and caffeine showed a moderate correlation. methamphetamine and 3,4-methylenedioxymethamphetamine (MDMA) and nicotine did not correlate with other sources of data. Most of the poor correlations were explained as part of the uncertainties related with the use estimates and were improved with other complementary sources of data.
This work confirms the promising future of WBE as a complementary approach to obtain a more accurate picture of substance use situation within different communities. Our findings suggest further improvements to reduce the uncertainties associated with both sources of information in order to make the data more comparable.
The chemical analysis of the exogenous biomarkers of drug consumption in wastewater has been shown to be an interesting approach for studying drug use within a defined population . The approach, termed wastewater-based epidemiology (WBE), was first applied in 2008 to study illicit drug use in three Italian cities  and it has since been extended to a wider range of exogenous biomarkers specific to the use of alcohol [3–5], tobacco , caffeine  and selected pharmaceuticals [8, 9]. To date there has been very limited focus on evaluating the relationship between WBE data to other drug use measures.
The use of pharmaceuticals in society is usually determined by prescription data or self-reports. However, pharmaceutical prescription or sales data involve several sources of uncertainties such as the inappropriate use of a medication, household disposal through the toilet or sink, sales without prescription or illegal acquisition [10–12]. Self-reports on adherence are potentially subject to recall and reporter bias . An alternative could be to quantitatively measure pharmaceuticals in wastewater as an indicator of use [14–16].
Europe has a significant market for illicit drugs sustained by both domestic production and import from other regions . Different monitoring approaches have been used to investigate the scale of drug use. Police and custom’s seizures, population surveys, hospital admissions, drug-overdose deaths or treatment programs have been the most reliable data for drug use estimates [18–20]. Under-reporting is still apparent despite the use of a number of different sources of information to improve the level of accuracy of self-reporting . Seizure data may not give the correct picture and surveys suffer from reporting biases and low response rates . WBE has been introduced as a promising alternative or addition to these sources of data [23–25].
There is also a need to be able to accurately determine the consumption of licit drugs, such as alcohol, nicotine and caffeine. Alcohol consumption can be estimated from sales statistics but trafficking, illegal production, stock piling or waste can affect such data. Reid and co-workers for the first time applied WBE to obtain complementary information on community alcohol consumption by analysing the alcohol metabolite, ethyl sulfate (EtS) in wastewater . This approach has since been used in other studies to report the rates of alcohol consumption in populations [3, 26], but the value of these biomarkers has yet to be determined. Tobacco use is commonly deduced from tobacco sales statistics. Prevalence surveys and the aggregate between production and trade statistics can also yield the total consumption estimates . However, these data may bias the total consumption estimates (similar to alcohol) indicating that additional methods are required in order to improve tobacco use estimations. WBE has also been used to determine nicotine use through the measurement of its residues in wastewater at national level [6, 7, 28]. These data were compared with local tobacco sales statistics finding good agreement . Caffeine is present in coffee and tea as well as other products such as energy drinks, chocolate and certain medicines. Caffeine occurrence has been investigated to assess its patterns of use across different communities [7, 30].
Very few studies to date have correlated wastewater data with other sources of data. So far, pharmaceutical concentrations in wastewater have been related to the amounts used in a population by using prescription data [11, 31–33]. A Spanish study compared wastewater concentrations of 17 pharmaceuticals with annual sales data and found a good correlation . The same was true for studies from Australia , Italy , and Belgium . One study compared cocaine loads in wastewater in Oslo with cocaine measures from population survey data and drug tests among car drivers and found a good correlation . A recent study in Belgium compared WBE results with survey data collected simultaneously in a small community of around 30,000 inhabitants for alcohol and nicotine metabolites . No significant correlations were found, but the study was impaired by a typically low response to the surveys. Discrepancies between wastewater and national survey data have also been found in Australia .
The present study aims to compare WBE data with other sources of drug use data in order to assess its correlation. The study uses concentrations of pharmaceuticals, illicit drugs, alcohol, nicotine and caffeine in wastewater samples from different European cities and compares these with other sources of data for drug use. The data used for comparison were pharmaceutical prescriptions (morphine, methadone, atenolol, metoprolol, methylphenidate, carbamazepine, citalopram, paracetamol, oxazepam and diclofenac), police drug seizure data for illicit drugs (cocaine, amphetamine, and MDMA (ecstasy)), and consumption data for alcohol, tobacco and caffeine metabolites (EtS; cotinine, nicotine and hydroxycotinine; caffeine, paraxanthine, 1-methylxanthine, 7-methylxanthine, 1-methyluric acid and 1,7-dimethyluric acid).
Twenty-four hour composite inlet sewage samples were collected from 8 different European wastewater treatment plants (WWTP) over 7 consecutive days in March 2015. All samples were stored in plastic containers, immediately frozen at −20 °C to prevent degradation of the urinary metabolites and sent to Oslo and Milan within 24 h in cool boxes with dry ice or ice packs to keep the samples frozen. The average wastewater flow rate (litres per day (L/day)) through each of the treatment plants was recorded for each of the daily samples. Samples were collected every morning between 07:00 and 8:30 except for Zurich where it was performed at 00:00 (see sample details in Additional file 1: Table S1).
Extraction and analysis
The analysis of the different wastewater samples was performed using three previously validated methods [4, 7, 9]. Illicit drugs and pharmaceuticals were analysed by liquid chromatography coupled to quadrupole time-of-flight mass spectrometry (LC-QTOF) as described by Baz-Lomba and co-workers . Briefly, 50 ng of the isotopic labelled internal standards (ILIS) solution mix was spiked into wastewater (100 mL) and then extracted on a fully automatable solid phase extraction (SPE) (Horizon Technology, Salem, NH, USA) with HLB extraction disks (47 mm, I.D.; Horizon Technology, City, Country). 5 μL of the final eluent were injected into the LC-QTOF system. The compounds were chromatographically separated on a Waters Acquity UPLC system (Milford, MA, USA) fitted with a Acquity UPLC HSS C18 column (1.8 μm, 2.1 mm × 150 mm) (Waters, Milford, MA, USA). A Xevo G2-S Q-TOF mass spectrometer (Waters, Milford, MA USA) was used in positive ESI mode for acquisition using MSe, that allows both precursor and product ion data to be simultaneously acquired during a single run. Analyses were performed using an external reference (Lock-Spray™) with Leucine-enkephalin generating a reference ion in positive mode at m/z 556.2771 that was used for real-time mass corrections in order to maintain the mass accuracy.
EtS was analysed by ion-pair liquid chromatography coupled to tandem mass spectrometry (ion-pair LC-MS/MS) as described by Reid and co-workers . Briefly, an aliquot (1 ml) of each sample was spiked with deuterated internal standard (EtS-d5) at 50 ng/mL and centrifuged at 20,000 g for 10 mins. The supernatant was subsequently injected into the LC-MS/MS system (Waters, Milford, MA) equipped with an Acquity UPLC BEH C8 column (1.7 μm, 2.1 mm × 50 mm) (Waters, Milford, MA, USA). Mobile phases were 7 mM dihexylammonium acetate in water and methanol prepared by adding equimolar volumes of dihexylamine and acetic acid into the solvents.
The tobacco and coffee biomarkers were analysed by liquid chromatography-tandem mass spectrometry (LC–MS/MS) (API 5500 QqQ. Applied Biosystems-Sciex, Thornhill, Ontario, Canada) . Wastewater samples were filtered through 1.6 μm GF/A glass microfiber filters and 0.45 μm mixed cellulose membrane filters purchased from Whatman (Kent, UK). An aliquot (3 mL) was spiked with the mix labelled internal standard solution and loaded into the Oasis HLB SPE cartridges. For analysis, the extract was evaporated to dryness and reconstituted in 100 μL methanol-ultrapure water (20:80, v/v), centrifuged and transferred into glass vials for instrumental analysis. Finally, 1 μL of the final extract were injected in a LC–MS/MS system. Chromatographic separation was performed using a HPLC XTerra C18 column (3.5 μm, 1 mm × 100 mm) (Waters, Milford, MA, USA) and analytes were ionized using electrospray ionization in positive mode.
Drug consumption in the population connected to the respective WWTP was calculated by multiplying the measured concentration (Cww) of each compound (ng/L) by the average daily flow rate (qday) (m3/day) to obtain the daily mass loads (mg/day). Mass loads were then divided by the population served within the catchment area (Pww) to obtain the amount of the drug consumed per day per 1000 inhabitants expressed as population-normalised mass loads (Eq. 1). Correction factors (kCF) to account for the different excretion patterns were used only for pharmaceuticals. Moreover, concentration values that were < LOQ (method limit of quantification) were replaced by LOQ/2 if at least one day in the week had a concentration value above the LOQ .
Pharmaceuticals: Wastewater analysis vs sales data
The environmental population-normalised loads for the selected pharmaceuticals in wastewater were estimated for Oslo (Norway) from the per capita monthly sales data from 2012 and 2013 (Norwegian Drug Wholesales statistics, Norwegian Institute of Public Health). Norway has one of the most accurate national prescription and sales systems in the world and all data are electronically available . The possibility of gathering the sales data for the area of interest decreases the level of uncertainty. Other studies [31, 33] used national-based data and assumed that consumption pattern was the same for the whole national population.
VEAS, the Oslo WWTP, treats sewage for around 600,000 people of which Oslo contributes with the 70.5 % and Akershus with the 29.5 % (8 % from Asker and 21.5 % from Bærum). The monthly sales data (gathered for Oslo and Akershus) from 2012 and 2013 remain constant for all the pharmaceuticals during the course of the two years. For the calculation of the predicted loads the amount of defined daily doses (DDD) for each pharmaceutical was multiplied by the monthly turnover by dosage (TD) and divided by the number of days of each month. The measured loads for the pharmaceuticals in wastewater were in this case multiplied by the correction factors that took into account the different urinary excretion rates.
Illicit drugs: Wastewater analysis vs drug seizures
Measured population-normalised mass loads of benzoylecgonine, total amphetamines (as a sum of amphetamine and methamphetamine) and MDMA in wastewater were compared with the seizure data available from the 2015 European Drug Report . This report provided both the number and amount of seizures in 2013 (most recent data) for each substance except for Switzerland and the Netherlands. Bramness and co-workers  recommended the use of number of seizures for this kind of comparison since it provides a more accurate description of the drug situation rather than the amount seized which is more vulnerable to variations. The correlations for the three drugs are calculated based on the number of seizures in 2013 as recommended previously, however, further compared with the mean amount seized from 2010–2013 reported by the United Nations Office on Drugs and Crime (UNODC).
Legal stimulant drugs: Wastewater analysis vs sales data
Three legal stimulants were measured for comparison with sales and consumption. The alcohol consumption rates measured in wastewater were compared with the latest surveillance report on alcohol and health published by the WHO (World Health Organization, 2014) . The WHO reported consumption rates for the population aged 15 or older as L/year/person. Consumption rates were corrected for the entire population and then compared with the population-normalised loads of EtS measured in the wastewater.
The measured population-normalised mass loads of nicotine in wastewater were compared with tobacco sales data (presented as number of cigarettes per smoker per day) from 2014 as reported by the Tobacco Atlas (Euromonitor International Database) . These data were normalised by multiplying the corresponding correction factor for the percentage of population aged 15 or older. The mass loads for nicotine measured in wastewater were calculated as suggested by Castiglioni and co-workers  by summing up the loads of the two nicotine main metabolites, cotinine and hydroxycotinine.
The population-normalised mass loads of caffeine and its metabolites measured in wastewater were compared with the per capita coffee consumption data from 2013 referred to the total dry weight of coffee consumed rather than brewed volume (more sensitive to variations from market to market) . The correlation was performed with two approaches, by comparing i) the consumption data with the mass loads of caffeine, and ii) the sum of the six metabolites together. The interpretation of the wastewater results for caffeine needs to be carefully addressed since caffeine is present in many products, such as tea, energy drinks and chocolate.
Descriptive statistics for each data sources are presented as means. Spearman’s correlation coefficient (ρ) was used to assess the agreement between wastewater data and pharmaceutical sales data, police seizures and sales data, respectively for pharmaceuticals, illicit drugs and legal stimulants. P values less than 0.05 were considered statistically significant. All the analyses were performed in R, Version 3.2.2 .
Sales data for all of the targeted pharmaceuticals during the 24 months between 2012 and 2013 remained constant in Oslo. The coefficient of variation for sales varied between 7 % (oxazepam) and 14 % (carbamazepine). All of the selected pharmaceuticals were detected in all of the wastewater samples collected from Oslo. Paracetamol, metoprolol and carbamazepine showed the highest mean population-normalised mass loads (11.2, 6.5 and 2.4 g/day/1000 inh. respectively), while methadone and morphine were present at the lowest loads (19.0 and 94.2 mg/day/1000 inh. respectively). The day-to-day variability for each of the compounds over 7 days was relatively low with some exceptions (mean relative standard deviation (RSD) = 24 %, excluding methadone, morphine and methylphenidate). Methadone and morphine loads varied considerably during the 7 days without any pattern (RSD = 61 and 53 % respectively) while methylphenidate showed a high RSD (93 %), mainly due to the increased loads at the weekend.
A correction factor for the different urinary excretion rates was applied to make the data comparable with the pharmaceutical sales data for Oslo (Table 1). Sales data for all of the targeted pharmaceuticals for the 24 months between 2012 and 2013 remained constant in Oslo. The coefficient of variation for sales varied between 7 % (oxazepam) and 14 % (carbamazepine). Spearman’s rank correlation analysis was performed between the sales data in Oslo and the results obtained from the wastewater analysis for 10 different compounds showing a significant correlation (ρ = 0.85, P < 0.01) (Fig. 1).
Benzoylecgonine was present in all of the wastewater samples with population-normalised loads ranging between 151.5 (Oslo) and 672.6 mg/day/1000 inh. (Zurich) determined. With the exception of Oslo, all of the loads showed similar values for the cities studied. The population-normalised loads for amphetamine and methamphetamine found in wastewater from Oslo were similar (122.3 and 172.4 mg/day/1000 inh., respectively). Amphetamine was also measured in Bristol and Utrecht while methamphetamine was measured in Brussels, Milan, Zurich, Copenhagen, and at very low concentrations in some of the samples from Castellon. MDMA was detected in wastewater from all locations but not in every sample. The mean daily population-normalised loads for MDMA ranged from 5.4 mg/day/1000 inh. for Castellon to 121.5 mg/day/1000 inh. for Utrecht.
The measured wastewater loads for illicit drugs were compared with national drug seizure data  (Table 2). The lowest and highest number of cocaine seizures and quantity seized in 2013 was in Norway and Spain respectively. The Netherlands, Belgium and Italy also reported large amounts of cocaine while the United Kingdom reported the second highest number of seizures within the studied countries. The reported seizures for amphetamine were much higher in the United Kingdom. In general all the countries reported similar values except for Italy with very low number of seizures and amount seized reported. The Netherlands and Switzerland did not report data on amphetamine and methamphetamine. In general, few seizures were reported on methamphetamine. Norway presented the highest total number of seizures in 2013 (4,210 seizures) and mean amount seized from 2010–2013 (161 kg/year). For MDMA, the United Kingdom reported the highest number of seizures in 2013 (3,716 seizures) while The Netherlands the highest mean amount seized between 2010 and 2013 (1,750,867 tablets/year).
Correlation analysis was performed between wastewater data measured in 2015 and the number of seizures reported in 2013. The comparison for the illegal drugs was performed using the most recent data available (2015 European Drug Report with data from 2013 ). As for pharmaceuticals, there was a time mismatch between both datasets, however, seizure trends remained relatively stable during the last years. A significant positive correlation was found between the measured loads of benzoylecgonine and the number of seizures of cocaine reported by the different countries (ρ = 0.89, P = 0.03), while non-significant correlation was found between wastewater results and the number of seizures for amphetamines (ρ = 0.49, P = 0.36) and MDMA (ρ = 0.26, P = 0.66). Wastewater data were also compared with the mean amount seized from 2010–2013 reported by the UNODC, but non-significant correlation was found between the amount seized and the wastewater data for cocaine, amphetamines and MDMA.
Alcohol, nicotine and caffeine
Table 3 shows the measured population-normalised mass loads, sales and consumption data for the three legal drugs studied; alcohol, nicotine and caffeine. Copenhagen, Brussels and Oslo showed the highest EtS population-normalised mass loads in wastewater (7.7, 5.6 and 4.9 g/day/1000 inh. respectively) while Castellon and Milan showed the lowest level of alcohol use during the wastewater study period (both with 1.7 g/day/1000 inh). The reported WHO alcohol consumption rates were very similar for all the cities within a range of 11.9–14.1 L/year/person except for Italy and Norway, which showed lower levels (7.8 and 9.5 L/year/person respectively). Spearman’s rank correlation analysis showed, however no-correlation between the two data sources (ρ = 0.53, P = 0.18).
The measured load profiles of nicotine in wastewater varied among the 8 cities. The measurements in Oslo showed the highest mass loads level with 4.0 g/day/1000 inh. followed by Zurich, Brussels and Castellon with 3.6, 3.1 and 3.1 g/day/1000 inh. respectively. Milan and Utrecht showed the lowest population-normalised mass load rates (1.5 and 1.3 g/day/1000 inh. respectively). These data have been compared with tobacco consumption estimates provided by The Tobacco Atlas . The Tobacco Atlas reports Belgium as the country with the highest consumption rate of cigarettes (7.8 cigarettes/person/day) whereas Norway is reported as the country with the lowest rate (1.9 cigarettes/person/day). However, our analysis showed a non-significant correlation between the two data sources (ρ = −0.09, P = 0.84).
Oslo had the highest measured population-normalised loads for caffeine with a value of 42.9 g/day/1000 inh. followed by Copenhagen and Brussels with 39.4 and 25.8 g/day/1000 inh. respectively. By summing up all the caffeine metabolites Oslo still had a high mass load rates (129.7 g/day/1000 inh.), but Zurich and Copenhagen with 141.0 and 154.1 g/day/1000 inh. respectively were reported as the highest. There was no significant correlation between the measured loads of caffeine in wastewater and the caffeine consumption rates (ρ = 0.64, P = 0.09). The sum of all the caffeine metabolites measured in wastewater did not provide a better correlation (ρ = 0.33, P = 0.42).
This study compared the population-normalised mass loads of three different groups of compounds measured in wastewater with other drug use data sources. The pharmaceutical consumption was compared at local level in Oslo (Norway) while the comparison between wastewater data and other sources of data for the illicit and legal stimulant drugs was performed at international level. The wastewater data for the pharmaceuticals and also cocaine showed a high correlation with other data sources for use. For amphetamines, alcohol and caffeine a moderate, but non-significant correlation was observed, while wastewater data on MDMA and nicotine did not correlate with other sources of data.
The correlation between wastewater data for pharmaceuticals and sales data corroborates the findings of an earlier study finding good agreement between the sales data of metoprolol and cetirizine and measured loads of these drugs in wastewater in Oslo . In addition, the pharmaceuticals wholesale data used in the current study were from Oslo only and extracted only for the month of interest. Hence, the accuracy may be higher than national consumption data [31–33], which does take into consideration different local patterns of use , or prescription databases that may miss the contributions from hospital in-patients and sales of over-the-counter drugs . The current results seem to suggest that WBE data represent accurately the population drug use that is usually measured by sales data, confirming that wastewater data could be used with certainty to estimate drug use in a population.
Among the studied illicit drugs, cocaine measurements in wastewater showed the best agreement with seizure data. The wastewater results for cocaine further agree with the latest WBE multi-city study performed in Europe in 2014 in which United Kingdom, The Netherlands, Belgium, Switzerland and Spain reported the highest population-normalised loads in 2014 . In the same study, Utrecht and Castellon showed the lowest national population-normalised loads of cocaine measured in wastewater for their respective countries. These national spatial differences can explain some of the discrepancies between the WBE and the national seizure data.
For all the other comparisons in the present study the relationship between wastewater data and other sources of drug use data were less than optimal. Wastewater data presents intra-country and spatial variability and the patterns of drug use also vary over the course of the year , while seizures and consumption data represents an annual average within the whole country. WBE results have already shown different patterns of drug use within the same country [23, 45]. For many of our analysis we had very few observations making it possible for spurious measures to influence the analysis unduly. Moreover, the different sources of data may not compare because they measure different things (E.g. number of seizures and amount seized data did not agree completely, but together may improve the general picture of the drug situation in Europe describing the main points of entrance and possible trafficking distribution). In an Australian study the wastewater data was contrary to the drug use surveys but in agreement with prevalence data, police, health agencies and the media . The negative results should thus be interpreted carefully.
The sum of the amphetamine and methamphetamine measured mass loads in wastewater correlated only moderately well with the number of seizures. A contributing factor to the low agreement may be that amphetamines have a more defined “geographical drug-use patterns” than cocaine. While amphetamine is more usual in western Europe, methamphetamine use is mainly concentrated in northern Europe, Czech Republic and Slovakia . This is corroborated by previous WBE studies [23, 24, 38, 43, 46] and prevalence data .
No relationship was observed between WBE data for MDMA and the number of seizures. Despite this, wastewater results seem to be in agreement with other studies. For example, Utrecht and Bristol showed the highest population-normalised WBE loads of MDMA which is in harmony with survey data from other sources . MDMA is however a drug used more typically on special occasions and as a “party drug”, making it more sensitive to timing of WBE . Also, the MDMA market has changed dramatically during the last couple of years  and our poor results may have been worsened by the fact that we compared wastewater results from 2015 with seizures data from 2013. Furthermore, we only had six observations due to lack of seizure data from some countries .
Only a moderate correlation was found between EtS loads in wastewater and WHO sales data corroborating other studies finding similar discrepancies. A study from Santiago de Compostela and Milan estimated alcohol use on wastewater data and found results 39 and 68 % lower than the estimates reported by the WHO . A Belgian study found that WBE data for alcohol consumption were in agreement with values reported by the Belgian Scientific Institute of Public Health, but much lower than those reported by the WHO . Lastly, a multi-city international study of alcohol consumption for 20 different cities showed that a significant correlation was found for only 11 cities between the WHO estimations and the wastewater measurements . According to these studies the main limitations are related to suboptimal records over alcohol sales in many countries. The strict Norwegian government control over alcohol sales results in accurate sales statistics . Still, the unregistered consumption of alcohol increased to around 26 % in 2012 which would explain the high alcohol loads measured in wastewater according to what expected by WHO . In a previous study in Oslo, the alcohol use estimates obtained from WBE and those reported by the WHO were in good agreement . This indicates that WBE could be a suitable marker for alcohol consumption even with the observed moderate correlation with consumption data.
No relationship was found between wastewater loads for nicotine and cigarette sales data. Here Oslo appeared to be an outlier. The popularity of dipping tobacco in Norway has been shown to have risen between 2010–2015, reaching the same prevalence levels as cigarettes (10 % against the 13 % for cigarettes) . The levels of cotinine excreted from dipping tobacco users are in the same range than cigarettes  which would induce an overestimation when comparing wastewater with cigarettes consumption data. In addition, a large amount of dipping tobacco is disposed of into urinals and toilets resulting in increased nicotine residues in Oslo´s wastewater influent. Cigarette consumption may vary greatly within countries  and other studies have shown a better agreement when comparing WBE with smoking prevalence data at regional level [6, 7, 29, 52]. Nicotine metabolites have been proposed and used as quantitative biomarkers to estimate population size at local level  although the presented results and previous studies  demonstrate that the different tobacco culture-dependent and geographic behaviour may change between countries. Therefore, nicotine metabolites should not be used as a population biomarker when comparing data across different countries
Coffee consumption data and the loads of the main caffeine metabolites in wastewater showed a moderate correlation between the two data sources. There was, however, a decent relationship when the observations from Bristol were removed from the calculations. The high WBE loads of caffeine measured in Bristol could be related to the by far highest tea , chocolate  and caffeinated energy drinks consumption (also in Spain ). We compared the consumption data with a sum of all the caffeine metabolites resulting in an even lower agreement. Other studies have also found that metabolite excretion profiles of caffeine disagree with WBE . This may have different explanations, one being large variations in caffeine content between different coffee crops .
As described above, there are still a large number of limitations for both types of epidemiological data. A precise ecological approach would require data gathered from the same population and at the same time, but the greatly different time-scales involved in both approaches and the logistics of sampling make this a very challenging task.
Despite the relatively good agreement with the other sources of data, it is possible that wastewater data may be typified by low temporal representativeness and high spatial variability due to different drug use and availability trends from location to location and over time and therefore, WBE results need to be carefully interpreted. Thus far, coverage of most of the WBE studies performed is typically limited to a one-week sampling program, however, for the future more extensive random-stratified sampling schemes are recommended to increase the representativeness and decrease the level of uncertainty.
Furthermore, because WBE data is an aggregate of the consumption of all segments of the population, they do not provide an easy means of assessing the importance of sex, age or socio-economical information on the findings.
The strength of WBE is however in the ability to more easily carry out longitudinal studies. Data can be provided in a more continuous flow rather than the infrequent and perhaps sporadic data-points that may be generated from traditional surveys. WBE could therefore prove to be an excellent source of data for studying the effects of police operations, education programs or public heath campaigns aimed at reducing the use of a particular drug.
Wastewater analysis and conventional epidemiological indicators are not expected to be in exact agreement, but a degree of overlap should exist and be sufficient to demonstrate the complementary character of these approaches, as demonstrated by the case studies presented in this work.
A relatively good relationship between the WBE data and the other sources of information, especially for pharmaceuticals and cocaine, was found even with the limited number of cities in the present study. This study illustrates that WBE is in reality a snapshot of the drug situation in a determined location while other data often are more general, like annual prevalence or sales data, being less sensitive to quicker changes in the recent drug use patterns. Further research is needed to improve the accuracy on the consumption reports and the uncertainties associated with the WBE. The different nature of the datasets can lead to certain disagreement when comparing both sources of information. Despite this, this work has shown that WBE is currently providing complementary and valuable information to improve the general picture of the drug situation in Europe. The presented results although foreseeing a promising future, need to be interpreted carefully since they provide a rough idea about the level of agreement between WBE and other sources of data together with some interesting indicators on the current drug situation, but both sources are not exempt from uncertainties.
defined daily doses
isotopic labeled internal standards
litters per day
liquid chromatography-tandem mass spectrometry
liquid chromatography coupled to a quadrupole time-of-flight mass spectrometry
method limit of quantification
population served within the catchment area
wastewater daily flow rate
relative standard deviation
solid phase extraction
monthly turnover by dosage
United nations office on drugs and crime
World health organization
wastewater treatment plants
Thomas KV, Reid MJ. What else can the analysis of sewage for urinary biomarkers reveal about communities? Environ Sci Technol. 2011;45(18):7611–2.
Zuccato E, Chiabrando C, Castiglioni S, Bagnati R, Fanelli R. Estimating community drug abuse by wastewater analysis. Environ Health Perspect. 2008;116(8):1027–32.
Rodríguez-Álvarez T, Racamonde I, González-Mariño I, Borsotti A, Rodil R, Rodríguez I, Zuccato E, Quintana JB, Castiglioni S. Alcohol and cocaine co-consumption in two European cities assessed by wastewater analysis. Sci Total Environ. 2015;536:91–8.
Reid MJ, Langford KH, Mørland J, Thomas KV. Analysis and Interpretation of Specific Ethanol Metabolites, Ethyl Sulfate, and Ethyl Glucuronide in Sewage Effluent for the Quantitative Measurement of Regional Alcohol Consumption. Alcohol Clin Exp Res. 2011;35(9):1593–9.
Boogaerts T, Covaci A, Kinyua J, Neels H, van Nuijs AL. Spatial and temporal trends in alcohol consumption in Belgian cities: A wastewater-based approach. Drug Alcohol Depend. 2016.
Castiglioni S, Senta I, Borsotti A, Davoli E, Zuccato E. A novel approach for monitoring tobacco use in local communities by wastewater analysis. Tob Control. 2015;24(1):38–42.
Senta I, Gracia-Lor E, Borsotti A, Zuccato E, Castiglioni S. Wastewater analysis to monitor use of caffeine and nicotine and evaluation of their metabolites as biomarkers for population size assessment. Water Res. 2015;74:23–33.
Hernández F, Ibáñez M, Bade R, Bijlsma L, Sancho JV. Investigation of pharmaceuticals and illicit drugs in waters by liquid chromatography-high-resolution mass spectrometry. TrAC Trends Anal Chem. 2014;63:140–57.
Baz-Lomba JA, Reid MJ, Thomas KV. Target and suspect screening of psychoactive substances in sewage-based samples by UHPLC-QTOF. Analytica Chimica Acta. 2016;914:81–90.
Bound JP, Voulvoulis N. Household disposal of pharmaceuticals as a pathway for aquatic contamination in the United Kingdom. Environ Health Perspect. 2005;113(12):1705–11.
van Nuijs AL, Covaci A, Beyers H, Bervoets L, Blust R, Verpooten G, Neels H, Jorens PG. Do concentrations of pharmaceuticals in sewage reflect prescription figures? Environ Sci Pollut Res Int. 2015;22(12):9110–8.
Venhuis BJ, de Voogt P, Emke E, Causanilles A, Keizers PH. Success of rogue online pharmacies: sewage study of sildenafil in the Netherlands. BMJ. 2014;349(7965):g4317.
Brier MJ, Chambless D, Gross R, Su HI, DeMichele A, Mao JJ. Association between self-report adherence measures and oestrogen suppression among breast cancer survivors on aromatase inhibitors. Eur J Cancer. 2015;51(14):1890–6.
Ort C, Lawrence MG, Rieckermann J, Joss A. Sampling for pharmaceuticals and personal care products (PPCPs) and illicit drugs in wastewater systems: are your conclusions valid? A critical review. Environ Sci Technol. 2010;44(16):6024–35.
Osorio V, Marcé R, Pérez S, Ginebreda A, Cortina JL, Barceló D. Occurrence and modeling of pharmaceuticals on a sewage-impacted Mediterranean river and their dynamics under different hydrological conditions. Sci Total Environ. 2012;440:3–13.
Pereira AM, Silva LJ, Meisel LM, Lino CM, Pena A. Environmental impact of pharmaceuticals from Portuguese wastewaters: geographical and seasonal occurrence, removal and risk assessment. Environ Res. 2015;136:108–19.
EMCDDA. European Drug Report 2015: Trends and Developments. European Monitoring Centre for Drugs and Drug Addiction: Lisbon; 2015.
Kraus L, Augustin R, Frischer M, Kummler P, Uhl A, Wiessing L. Estimating prevalence of problem drug use at national level in countries of the European Union and Norway. Addiction. 2003;98(4):471–85.
Colpe LJ, Barker PR, Karg RS, Batts KR, Morton KB, Gfroerer JC, Stolzenberg SJ, Cunningham DB, First MB, Aldworth J. The National Survey on Drug Use and Health Mental Health Surveillance Study: calibration study design and field procedures. Int J Methods Psychiatr Res. 2010;19(S1):36–48.
Kroutil LA, Vorburger M, Aldworth J, Colliver JD. Estimated drug use based on direct questioning and open‐ended questions: responses in the 2006 National Survey on Drug Use and Health. Int J Methods Psychiatr Res. 2010;19(2):74–87.
Reid MJ, Langford KH, Grung M, Gjerde H, Amundsen EJ, Morland J, Thomas KV. Estimation of cocaine consumption in the community: a critical comparison of the results from three complimentary techniques. BMJ Open. 2012;2(6):e001637.
Castiglioni S, Borsotti A, Riva F, Zuccato E. Illicit drug consumption estimated by wastewater analysis in different districts of Milan: A case study. Drug Alcohol Rev. 2014;35(2):128-32.
Ort C, van Nuijs AL, Berset JD, Bijlsma L, Castiglioni S, Covaci A, de Voogt P, Emke E, Fatta-Kassinos D, Griffiths P, et al. Spatial differences and temporal changes in illicit drug use in Europe quantified by wastewater analysis. Addiction. 2014;109(8):1338–52.
Thomas KV, Bijlsma L, Castiglioni S, Covaci A, Emke E, Grabic R, Hernández F, Karolak S, Kasprzyk-Hordern B, Lindberg RH, et al. Comparing illicit drug use in 19 European cities through sewage analysis. Sci Total Environ. 2012;432:432–9.
Castiglioni S, Thomas KV, Kasprzyk-Hordern B, Vandam L, Griffiths P. Testing wastewater to detect illicit drugs: state of the art, potential and research needs. Sci Total Environ. 2014;487:613–20.
Rodríguez-Álvarez T, Rodil R, Cela R, Quintana JB. Ion-pair reversed-phase liquid chromatography–quadrupole-time-of-flight and triple-quadrupole–mass spectrometry determination of ethyl sulfate in wastewater for alcohol consumption tracing. J Chromatogr A. 2014;1328:35–42.
Guindon GE, Boisclair D. Past, current and future trends in tobacco use. 2003.
Lopes A, Silva N, Bronze MR, Ferreira J, Morais J. Analysis of cocaine and nicotine metabolites in wastewater by liquid chromatography–tandem mass spectrometry. Cross abuse index patterns on a major community. Sci Total Environ. 2014;487:673–80.
Rodríguez-Álvarez T, Rodil R, Rico M, Cela R, Quintana JB. Assessment of Local Tobacco Consumption by Liquid Chromatography–Tandem Mass Spectrometry Sewage Analysis of Nicotine and Its Metabolites, Cotinine and trans-3′-Hydroxycotinine, after Enzymatic Deconjugation. Anal Chem. 2014;86(20):10274–81.
Bueno MJ, Ucles S, Hernando MD, Davoli E, Fernandez-Alba AR. Evaluation of selected ubiquitous contaminants in the aquatic environment and their transformation products. A pilot study of their removal from a sewage treatment plant. Water Res. 2011;45(6):2331–41.
Carballa M, Omil F, Lema JM. Comparison of predicted and measured concentrations of selected pharmaceuticals, fragrances and hormones in Spanish sewage. Chemosphere. 2008;72(8):1118–23.
Lai FY, Ort C, Gartner C, Carter S, Prichard J, Kirkbride P, Bruno R, Hall W, Eaglesham G, Mueller JF. Refining the estimation of illicit drug consumptions from wastewater analysis: co-analysis of prescription pharmaceuticals and uncertainty assessment. Water Res. 2011;45(15):4437–48.
Verlicchi P, Al Aukidy M, Jelic A, Petrović M, Barceló D. Comparison of measured and predicted concentrations of selected pharmaceuticals in wastewater and surface water: A case study of a catchment area in the Po Valley (Italy). Sci Total Environ. 2014;470–471:844–54.
van Wel JH, Gracia-Lor E, van Nuijs AL, Kinyua J, Salvatore S, Castiglioni S, Bramness JG, Covaci A, Van Hal G. Investigation of agreement between wastewater-based epidemiology and survey data on alcohol and nicotine use in a community. Drug Alcohol Depend. 2016;162:170–5.
Tscharke BJ, Chen C, Gerber JP, White JM. Trends in stimulant use in Australia: A comparison of wastewater analysis and population surveys. Sci Total Environ. 2015;536:331–7.
Salvatore S, Bramness JG, Reid MJ, Thomas KV, Harman C, Roislien J. Wastewater-Based Epidemiology of Stimulant Drugs: Functional Data Analysis Compared to Traditional Statistical Methods. PLoS One. 2015;10(9):e0138669.
Norwegian Prescription Database. Available at: [http://www.norpd.no/]. Accessed April 2016.
Bramness JG, Reid MJ, Solvik KF, Vindenes V. Recent trends in the availability and use of amphetamine and methamphetamine in Norway. Forensic Sci Int. 2015;246:92–7.
Organization WH. Global status report on alcohol and health-2014. Geneva, Switzerland: World Health Organization; 2014. [http://apps.who.int/iris/bitstream/10665/112736/1/9789240692763_eng.pdf?ua=1]. Accessed April 2016.
Eriksen M, Mackay J, Schluger N, Gomeshtapeh FI, Drope J. The tobacco atlas, 5th edn. Georgia: The American Cancer Society, 2015. Available from: [http://3pk43x313ggr4cy0lh3tctjh.wpengine.netdnacdn.com/wp-content/uploads/2015/03/TA5_2015_WEB.pdf]. Accessed Sept 2016.
Caffeine (Coffee) Consumption By Country [http://www.caffeineinformer.com/caffeine-what-the-world-drinks]. Accessed April 2016.
Team RC. R: a language and environment for statistical computing, Vol 3.2.1. Vienna: R Foundation for Statistical Computing; 2013.
Reid MJ, Langford KH, Mørland J, Thomas KV. Quantitative assessment of time dependent drug-use trends by the analysis of drugs and related metabolites in raw sewage. Drug Alcohol Depend. 2011;119(3):179–86.
Ort C, Lawrence MG, Reungoat J, Eaglesham G, Carter S, Keller J. Determining the fraction of pharmaceutical residues in wastewater originating from a hospital. Water Res. 2010;44(2):605–15.
EMCDDA. Wastewater analysis and drugs: a European multi-city study. In: European Monitoring Centre for Drugs and Drug Addiction. 2015.
Kankaanpää A, Ariniemi K, Heinonen M, Kuoppasalmi K, Gunnar T. Use of illicit stimulant drugs in Finland: A wastewater study in ten major cities. Sci Total Environ. 2014;487:696–702.
Harman C, Reid M, Thomas KV. In situ calibration of a passive sampling device for selected illicit drugs and their metabolites in wastewater, and subsequent year-long assessment of community drug usage. Environ Sci Technol. 2011;45(13):5676–82.
Ryu Y, Barcelo D, Barron LP, Bijlsma L, Castiglioni S, de Voogt P, Emke E, Hernandez F, Lai FY, Lopes A, et al. Comparative measurement and quantitative risk assessment of alcohol consumption through wastewater-based epidemiology: An international study in 20 cities. Sci Total Environ. 2016.
Hordvin O. The Drug Situation in Norway 2012. Annual Report to the EMCDDA. Oslo: SIRUS; 2012.
Smoking habits, 2015 [https://www.ssb.no/en/helse/statistikker/royk/aar/2016-01-14#content]. Accessed April 2016.
Lee PN. Health risks related to dual use of cigarettes and snus – A systematic review. Regul Toxicol Pharmacol. 2014;69(1):125–34.
Mackuľak T, Birošová L, Grabic R, Škubák J, Bodík I. National monitoring of nicotine use in Czech and Slovak Republic based on wastewater analysis. Environ Sci Pollut Res. 2015;22(18):14000–6.
Chen C, Kostakis C, Gerber JP, Tscharke BJ, Irvine RJ, White JM. Towards finding a population biomarker for wastewater epidemiology studies. Sci Total Environ. 2014;487:621–8.
Where the world's biggest tea drinkers are [http://qz.com/168690/where-the-worlds-biggest-tea-drinkers-are/]. Accessed April 2016.
Alberts HC, Cidell J. Chocolate Consumption, Manufacturing, and Quality in Europe and North America. Econ Chocolate. 2015;119.
Zucconi S, Volpato C, Adinolfi F, Gandini E, Gentile E, Loi A, Fioriti L. Gathering consumption data on specific consumer groups of energy drinks. External Scientific Report for European Food Safety Authority. 2013;10(3):1–190.
Mitchell DC, Knight CA, Hockenberry J, Teplansky R, Hartman TJ. Beverage caffeine intakes in the U.S. Food Chem Toxicol. 2014;63:136–42.
Lienert J, Güdel K, Escher BI. Screening Method for Ecotoxicological Hazard Assessment of 42 Pharmaceuticals Considering Human Metabolism and Excretory Routes. Environ Sci Technol. 2007;41(12):4471–8.
ter Laak TL, Kooij PJ, Tolkamp H, Hofman J. Different compositions of pharmaceuticals in Dutch and Belgian rivers explained by consumption patterns and treatment efficiency. Environ Sci Pollut Res Int. 2014;21(22):12843–55.
Moffat AC, Osselton MD, Widdop B. Clarke’s analysis of drugs and poisons, vol. 3. London: Pharmaceutical press; 2011.
Faraj BA, Israili ZH, Perel JM, Jenkins ML, Holtzman SG, Cucinell SA, Dayton PG. Metabolism and disposition of methylphenidate-14C: studies in man and animals. J Pharmacol Exp Ther. 1974;191(3):535–47.
Shull HJ, Wilkinson GR, Johnson R, Schenker S. Normal disposition of oxazepam in acute viral hepatitis and cirrhosis. Ann Intern Med. 1976;84(4):420–5.
Jose Antonio Baz Lomba, Stefania Salvatore, Richard Bade, Erika Castrignanò, Ana Causanilles, Juliet Kinyua, Ann-Kathrin McCall, Pedram Ramin, Nikolaos I. Rousis, and Yeonsuk Ryu acknowledge the EU Marie-Skłodowska Curie Initial Training Network SEWPROF (Marie Curie-FP7-PEOPLE, grant number 317205) for their Early Stage Researcher grant and Emma Gracia-Lor for her Experienced Researcher grant. We thank the people and agencies who assisted in the collection of the wastewater samples, in particular Pia Ryrfors and colleagues at Vestfjorden Avløpselskap (VEAS, Oslo, Norway).
Availability of data and materials
Request for the data used in this study should be send to. Jose Antonio Baz Lomba (email: Joseantonio.email@example.com).
JABL: drafted the manuscript and performed part of the wastewater analysis; SSA performed the statistical analysis and contributed to the preparation of the manuscript; EGR and YRY performed part of the wastewater analysis and contributed to the preparation of the manuscript; RB, EC, AC, FH, BKH, JK, AKM, AvN, CO, BGP, PR, NIR and PdV organised the collection of the wastewater samplers and provided relevant data for WBE calculations. MR, SC, JGB and KVT advised on the execution of the study and contributed to the preparation of the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics and consent to participate
Wastewater samples information for the 8 European cities studied in 2015. The estimated population using the sewer network is showed in brackets (thousands of inhabitants within the catchment area). (DOCX 18 kb)
About this article
Cite this article
Baz-Lomba, J.A., Salvatore, S., Gracia-Lor, E. et al. Comparison of pharmaceutical, illicit drug, alcohol, nicotine and caffeine levels in wastewater with sale, seizure and consumption data for 8 European cities. BMC Public Health 16, 1035 (2016). https://doi.org/10.1186/s12889-016-3686-5
- Wastewater-based epidemiology
- Drug consumption, Correlation, Europe-wide study, Seizures, Sales statistics