Investigation and Application of a New Passive Sampling Technique

Jul 24, 2017 - The DGT measurement of these drugs was not influenced by the pH (4 to 9) and the ionic strength (0.001 M - 0.1 M) and unaffected by the...
3 downloads 10 Views 1MB Size
Subscriber access provided by UNIV OF NEWCASTLE

Article

Investigation and application of a new passive sampling technique for in-situ monitoring of illicit drugs in waste waters and rivers Changsheng Guo, Tingting Zhang, Song Hou, Jiapei Lv, Yuan Zhang, Fengchang Wu, Zhendong Hua, Wei Meng, Hao Zhang, and Jian Xu Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b00731 • Publication Date (Web): 24 Jul 2017 Downloaded from http://pubs.acs.org on July 24, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 30

Environmental Science & Technology

1

Investigation and application of a new passive sampling technique for in-situ monitoring of illicit

2

drugs in waste waters and rivers

3

Changsheng Guo†, Tingting Zhang‡, Song Hou†, Jiapei Lv†, Yuan Zhang†, Fengchang Wu†,

4

Zhendong Hua‡*, Wei Meng†, Hao Zhang§, Jian Xu†*

5

† State Key Laboratory of Environmental Criteria and Risk Assessment, Chinese Research

6

Academy of Environmental Sciences, 100012 Beijing, China

7

‡ National Narcotics Laboratory, Drug Intelligence and Forensic Center of the Ministry of Public

8

Security, Beijing 100193, China

9

§ Lancaster Environment Centre, Lancaster University, Lancaster LA1 4YQ, United Kingdom

1

ACS Paragon Plus Environment

Environmental Science & Technology 10

TOC art

11 12

Contents brief

13

DGT technique has been developed and applied for the first time as passive samplers to

14

quantitatively measure illicit drugs in situ in waste waters and rivers.

2

ACS Paragon Plus Environment

Page 2 of 30

Page 3 of 30 15

Environmental Science & Technology

ABSTRACT

16

Illicit drugs constitute a class of emerging contaminants that has been drawing significant

17

concern due to its potent pharmacological and biological activities. In this study, an in situ passive

18

sampling approach that uses diffusive gradients in thin films (DGT) was successfully tested for

19

measuring ketamine (KET), methamphetamine (METH), and amphetamine (AMP) in water. The

20

diffusion coefficients of KET, METH, and AMP in diffusive gel were (8.13±0.12)×10−6,

21

(8.55±0.14)×10−6, and (7.72±0.18)×10−6 cm2 s−1 at 22 °C, respectively. The capacities of an XAD

22

binding gel for KET, METH, and AMP were 92, 57, and 45 µg per binding gel disc, which were

23

suitable for long-term environmental monitoring. The DGT measurement of these drugs was not

24

influenced by the pH (4 to 9) and the ionic strength (0.001 M - 0.1 M) and unaffected by the water

25

flow, demonstrating the effectiveness of the XAD-based DGT for the in situ monitoring of illicit

26

drugs. DGT samplers were deployed in a WWTP influent and natural rivers in Beijing, China. The

27

ng L−1 levels of the drugs were high in the wastewater influent and low in river waters, with an

28

insignificant fluctuation during the seven-day monitoring. The DGT-measured concentrations were

29

comparable to the average concentrations determined by SPE method, which suggested that the

30

average data measured by DGT could be substituted for high-frequency grab sampling. This study

31

has demonstrated systematically for the first time that DGT is effective and accurate for monitoring

32

illicit drugs in wastewater and surface waters, and provides a powerful tool to investigating the

33

presence, transport, and environmental behaviors of these drugs in the aquatic ecosystem.

34

INTRODUCTION

35

Synthetic illicit drugs constitute a heterogeneous group of substances that has entered the drug

36

market around the world in the past quarter of a century.1,2 Although their use in some areas is lower

37

than that of traditional illicit drugs, such as cocaine, opium, heroin, and cannabis, they have the 3

ACS Paragon Plus Environment

Environmental Science & Technology 38

potential to spread in some populations, particularly among young people during social events,

39

which can result in significant problems. Data from surveys on populations indicate that each year

40

approximately 230 million people consume drugs of abuse worldwide.3 By the end of 2014, the

41

registered users in China has reached up to 2.9 million, two times more than the registered

42

population in 2008.4 However, the total number of drug users is approximately 5 times of the

43

number of the registered users. The actual population of drug users in China exceeded 14 million in

44

2014.4 The direct economic loss caused by drug abuse is 78.6 billion per year.4

45

A recently proposed approach for estimating the use of methamphetamine (METH) and other

46

illicit drugs in populations by measuring their metabolic residues in urban wastewater5–7 might

47

provide information on these synthetic drugs. The monitoring of illicit drug concentrations in

48

wastewater treatment plant (WWTP) influents is proposed to estimate the consumption of these

49

psychoactive substances at the community level.5,7,8

50

Limited studies were conducted to investigate the eco-toxicological impacts of drugs of abuse

51

at the chronic low levels in the aquatic environment.9–11 Similarly, reports of the ecotoxicological

52

impacts of illicit drugs at environmentally relevant concentrations have been few.12,13 Because of

53

their intrinsic psychoactive characteristics and possible adverse effects to the non-target biota, the

54

presence of drugs of abuse in the environment has been redeemed as an important issue requiring

55

much attention.

56

Many studies have been carried out to study the fluctuations in drug concentrations in the

57

aquatic environment during weekends and holidays using traditional discrete sampling methods

58

because of the known use of these substances for leisure.6,14 Nevertheless, inconsistent daily and

59

seasonal variations were reported by these studies, resulting in unclear conclusions. Few data of the

60

time-averaged concentrations of these drugs in the receiving river waters or WWTP influents are 4

ACS Paragon Plus Environment

Page 4 of 30

Page 5 of 30

Environmental Science & Technology

61

available. Grab and composite sampling, which are traditional water-sampling approaches, are

62

effective for measuring the presence of illicit drugs, however, these sampling approaches only

63

provide data at the sampling time and cannot reflect the circumstance changes, for instance, the

64

changes of the flow rate, chemical inputs, and the precipitation influence.11 The polar organic

65

chemical integrative sampler (POCIS) was developed to sample trace organic compounds, including

66

METH and amphetamine (AMP).11,15,16 However, this approach has limitations as well; for instance,

67

the environmental conditions such as water flow, the turbulence, and temperature can influence the

68

sampling rates and hence the concentration accuracy of the target compounds in water.

69

The diffusive gradients in thin films (DGT) technique was developed to measure labile

70

inorganic species and trace organic compounds in situ in aquatic environments.17 It has been

71

reported that DGT can measure 55 elements simultaneously at the detection limits from 0.001 to 1

72

ng mL-1 after a 24-h deployment.18 A total of 34 pharmaceuticals and pesticides could be measured

73

by DGT.19 Its robustness, simplicity, independence on flow to use, and apparent high stability

74

against outside interference make it attractive for measuring the time-integrated concentrations of

75

trace organics in waters, without tedious laboratory calibrations.17-19 DGT has been used recently for

76

both the qualitative and quantitative evaluations of antibiotics

77

and soil environment, which has provided alternatives for the DGT application to monitor organic

78

compounds in the aquatic environment.

20,21

and bisphenols

22

in aqueous

79

China Narcotics Bureau stated that METH and ketamine (KET) account for the largest share

80

(>99%) of synthetic illicit drugs in the Chinese drug consumption market.4 Du et al. 6 collected

81

water samples from 36 WWTPs in eighteen big cities covering major geographic areas of China to

82

determine the consumption of METH and KET. Wide occurrence of METH in the country was

83

observed, but there was no clear trend geographically. For KET, its usage was generally greater in 5

ACS Paragon Plus Environment

Environmental Science & Technology 84

the southern China than in the northern China. In previous studies, the active grab-sampling

85

technique followed by the solid-phase extraction (SPE) was used for quantifying the concentration

86

of drugs. One recent work using o-DGT to measure a diverse suite of polar organic compounds

87

including pharmaceutical drugs was carried out, in which the diffusion coefficients and the diffusive

88

boundary layers were measured, and the o-DGT was calibrated in the laboratory solution, without

89

considering its performance under wide range of environmental conditions.19 Although DGT is a

90

well-established technique in its fundamental theory and principle, it cannot be applied to illicit

91

drugs without systematic investigations if the technique is to be used for fully quantitative

92

measurement of in situ chemical concentrations. The illicit drug compounds are different in

93

chemical properties, structures and functions from antibiotics and other pharmaceutical compounds

94

tested by DGT. The aim of the current research was to develop DGT as a passive sampling

95

technique for accurately measuring the concentrations of illicit drugs, namely, KET, METH and

96

AMP, in situ. Different parameters and environmental conditions and their effects on DGT

97

performance were tested and investigated systematically. Finally, the developed DGT technique was

98

applied and validated in field conditions by the deployment of DGT devices in situ in a WWTP

99

influent and in the surface water of rivers.

100

EXPERIMENTAL SECTION

101

Chemicals and Reagents.

102

KET, AMP, METH, and their deuterated internal standards (KET-d4, AMP-d8, and METH-d8)

103

were obtained from Cerilliant Corporation (Round Rock, TX). The properties of the drugs are

104

provided in Supporting Information (SI) Table S1. Stock solutions of KET, AMP, or METH at 100

105

mg L−1 were prepared in methanol [MeOH, high-performance liquid chromatography grade] and

106

stored at −20 °C in the amber bottles. MeOH and acetonitrile (ACN) were obtained from Fisher 6

ACS Paragon Plus Environment

Page 6 of 30

Page 7 of 30

Environmental Science & Technology

107

Scientific (Poole, UK). Deionized water was prepared by a Milli-Q system (Millipore, MA, US).

108

SPE cartridges (Oasis MCX, 60 mg, 3 mL) were obtained from Waters Corporation (Milford, MA,

109

US). Polyethersulfone (PES) filter membrane was purchased from Pall Co., US. Mixed cellulose

110

ester (MCE), hydrophilic polytetrafluoroethylene (PTFE), and nylon (NL) filter membranes were

111

purchased from Shanghai Anpel Scientific Instrument Co., China. All membranes are 25 mm in

112

diameter and 0.45 µm in pore size.

113

DGT and Gel Preparation

114

Standard DGT devices for trace organics were purchased from DGT Research Ltd., UK. The

115

standard DGT consists of an XAD 18 agarose binding layer (0.5-mm thickness), an agarose

116

diffusive layer (0.8-mm thickness), and a PES filter membrane (0.14-mm thickness). The area of the

117

exposure window of the DGT plastic molding is 3.1 cm2. The preparation of binding and diffusive

118

gels followed the well-established procedures.20,22 To prepare the XAD layer, in 10 mL of hot

119

agarose solution 2 grams of XAD18 powder were added, and the solution was well mixed. The

120

solution was then poured into the space between two hot glasses, and cooled down at ambient

121

temperature to form the gel. The gels were put into the milli-Q water for hydration, and stored in the

122

salt solution (0.05 M NaCl).

123

Assessment of Adsorption by Sampling Apparatus and Materials.

124

Four different filters and two diffusive layers were soaked in 10 mL solutions containing 20 µg

125

L−1 KET, METH, or AMP. The solution was shaken for 12 h. DGT plastic moldings were soaked in

126

200 mL solutions of 20 µg L−1 KET, METH, or AMP because of their large size. The adsorbed mass

127

of the illicit drugs was calculated by the difference of drug concentrations before and after the

128

deployment of the testing materials.

129

Uptake Kinetics and Elution Efficiencies. 7

ACS Paragon Plus Environment

Environmental Science & Technology 130

The XAD binding gel discs were placed in 10 mL solution containing 20 µg L−1 KET, METH,

131

or AMP with a background NaCl concentration of 0.01 M, and shaken for different time (0.5 min to

132

24 h). The optimal elution conditions were obtained by calculating the elution efficiencies of the

133

illicit drugs by eluting the XAD layers in 10 mL of MeOH, ACN, 5% ammonia solution in MeOH,

134

or ethyl acetate for at least 1 d. Masses of target drugs in the eluents and the immersion solution

135

were used for elution efficiency calculation.

136

Diffusion Coefficient Measurements.

137

A diaphragm diffusion cell, which comprised two chambers that were linked by a 1.5-cm

138

circular window that housed a 0.80-mm thick agarose-based diffusive layer, was employed to

139

measure the diffusion coefficients of the drug compounds. The two chambers were washed first with

140

MeOH and then with Milli-Q water before use. The concentration gradient was the sole driving

141

force for the analyte transfer between the two chambers. The matrix-matched solutions in both

142

chambers were the same (50 mL, the same pH and ionic strength), except that the source chamber

143

contained 2 mg L−1 KET, METH, or AMP, while no illicit drugs were present in the receptor

144

chamber. The solutions in the two chambers were well stirred throughout the experiment. For at

145

least 3 h, 0.2 mL samples of solutions were taken at 15-min intervals from each chamber. The

146

diffusion coefficient D was calculated by the slope of the linear plot of the measured masses of each

147

drug in the receptor chamber versus time, as follows:

148

D = ( slope)

149

where A is the area of the window linking two chambers, C is the drugs concentration in the source

150

chamber, and ∆g is the thickness of the diffusive layer.

151

Calculation of DGT Measured Concentrations.

152

∆g , CA

(1)

After the deployment time (t, sec), the accumulated mass (M, in ng) of the analyte in the 8

ACS Paragon Plus Environment

Page 8 of 30

Page 9 of 30

Environmental Science & Technology

153

binding layer can be related to its solution concentration (i.e., CDGT, µg/L) according to Eq. (2).

154

CDGT =

155

where ∆g (cm) is the thickness of both the filters and the diffusive gel, D is the diffusion coefficient

156

(cm2 s−1), and A is the area (cm2) of the DGT’s exposure window.

(2)

M is calculated by Eq. (3) if it can be eluted with a known volume of eluent (Ve, L).

157

158

M ∆g , DAt

M=

Ce (Ve + Vgel ) fe

,

(3)

159

where Ce (µg/L) is the analyte concentration in the eluate; Vgel (L) is the volume of the binding gel,

160

which is typically 0.15 mL; and fe is the elution efficiency of the analyte from the binding gel.

161

Effects of Ionic Strength, pH, and Deployment Time

162

The effects of pH and ionic strength on DGT performance was tested by deploying the

163

standard DGT devices in different solutions: (1) 2.5 L of 20 µg L−1 KET, METH, or AMP solution

164

(pH 7) with varying NaCl levels (0.001- 0.5 M), and (2) 2.5 L solution with 0.01 M of NaCl and 20

165

µg L−1 KET, METH, or AMP at different pH values. The DGT devices were placed in 2.5 L

166

solutions (pH 7) containing 20 µg L−1 KET, METH, or AMP and 0.01 M of NaCl, and withdrawn at

167

various time intervals (12- 96 h) to evaluate the effect of deployment times. The relationship

168

between the DGT-accumulated mass and the thickness of diffusive gels was determined by

169

deploying the DGT devices equipped with different thickness diffusive gels (0.05 cm to 0.175 cm)

170

in the 2.5 L solution for 12 h containing 20 µg L−1 KET, METH, or AMP and 0.01 M of NaCl.

171

Capacity and Competition Effects.

172

The capacity of the XAD binding gels assembled in the DGT devices for the uptake of illicit

173

drugs were measured by immersing the devices in the solutions for 12 h containing 0.01 M of NaCl

174

and different levels of drugs: KET (0.02 mg L−1 to 25 mg L−1), METH (0.02 mg L−1 to 25 mg L−1),

175

or AMP (0.02 mg L−1 to 20 mg L−1). 9

ACS Paragon Plus Environment

Environmental Science & Technology 176

The competition effects among three illicit drugs in the solution were determined at various

177

solution concentrations. DGT device was placed for 12 h in different solutions containing 0.01 M of

178

NaCl: (1) KET or METH at 10 µg L−1, and the other two at 100 or 1000 µg L−1; (2) AMP at 10 µg

179

L−1, and KET and METH at 100 or 1000 µg L−1.

180

DGT Application and Validation in Field Conditions.

181

Four DGT devices were fixed together, with the exposure windows outward, to compare the

182

DGT measured concentrations and the SPE measured concentrations. The DGT devices, together

183

with temperature button data loggers (Maxim Integrated Products, United States) set to record

184

temperature every 2 h, were deployed in WWTP influents and urban rivers (located in Beijing,

185

China) for 7 days. The WWTP is in the south of Beijing City, which can treat 400,000 tons of

186

domestic wastewater per day. The WWTP covers an area of approximately 160 km2, serving

187

810,000 people. The treatment techniques of the WWTP include primary treatments and two

188

separate secondary treatment procedures, namely, traditional and reversed A2/O, each of which has a

189

domestic wastewater capacity of 200,000 tons per day. Surface water samples were collected in 15

190

urban river sites in Beijing, as shown in Figure S1. The Wenyu River and its tributaries (Sha River,

191

Beixiao River, Qing River, Ba River, and Tonghui River) cover an area of 2478 km2 and are

192

receiving rivers of most wastewater from the urban areas of Beijing. Wastewater samples were

193

collected at 10:00 and 16:00 each day. A previous study indicated that although the levels of METH

194

and AMP in a sewage pipeline varied significantly from 22:00 to 6:00, their concentrations in the

195

influents of several WWTPs in Beijing were relatively stable.23 All samples were transferred within

196

30 min to the laboratory, where the concentrations of drugs in the water samples (filtered by

197

Whatman GF/F filters) were determined with the established method (details in Supporting

198

Information) after the pre-concentration of 100 mL of water samples by SPE using an Oasis MCX 10

ACS Paragon Plus Environment

Page 10 of 30

Page 11 of 30

Environmental Science & Technology

199

cartridge (60 mg, 3 mL). The retrieved XAD binding gels from the DGT devices were soaked in 10

200

mL of MeOH for 24 h. Before instrumental analysis, the eluents from both the SPE cartridges and

201

DGT gels were passed through 0.22-µm membranes, followed by transferring to the HPLC sample

202

vials.

203

Instrumental Analysis.

204

The target compound analysis procedures by ultra-performance liquid chromatography tandem

205

mass spectrometry (UPLC-MS/MS, AB SCIEX, USA) are provided in the SI.

206

RESULTS AND DISCUSSION

207

Adsorption by DGT Moldings, Diffusive Layers, and Filter Membranes

208

The possible adsorption of drugs by the standard DGT moldings, diffusive gels, and filter

209

membranes was assessed. The results presented in Figure 1 shows minimal adsorption of drugs onto

210

the DGT moldings and two diffusive gels (agarose and PAM) (40%) and MCE (>20%). They were minimally absorbed on the

212

PES membranes (METH > AMP (Figure S3), which was in accordance with the order of the logKow values of

221

these compounds. The average uptake rates in the first half hour for KET, METH, and AMP were 11

ACS Paragon Plus Environment

Environmental Science & Technology 222

1.75, 1.45, and 1.03 ng cm−2 min−1, respectively, which were faster than the uptake rate required by

223

DGT devices to establish the concentration gradients. Hence, the DGT technique could be applied.

224

Diffusion Coefficients.

225

The time dependence of the masses of the three drugs diffused through an open-pore agarose

226

gel in the diffusion cell is shown in Figure S4. The masses of the drugs in the receptor chamber

227

increased linearly with time, whereas the drug concentrations in the source chamber were constant.

228

The slope of the mass plot in the receptor chamber versus time (Figure S4) allowed the calculation

229

of D values for the three drugs at the experimental temperature (22 °C). The obtained D values were

230

(8.13 ± 0.12) × 10−6, (8.55 ± 0.14) × 10−6, and (7.72 ± 0.18) × 10−6 cm2 s−1 for KET, METH, and

231

AMP, respectively. The diffusion coefficient at 25 °C can be calculated from measured D with the

232

following equation:17,20,22

233

log DT =

D (273 + T ) 1.37023(T − 25) + 8.36 ×10−4 (T − 25)2 + log 25 (4) 109 + T 298

234

where DT (cm2 s−1) is the analyte diffusion coefficient at temperature T (°C), and D25 (cm2 s−1) is the

235

analyte diffusion coefficient at 25 °C.

236

Elution Efficiencies and Capacity.

237

Four different eluents were tested to elute the drugs from the XAD binding gels that were

238

preloaded with the chemicals (Table S4). The average elution efficiencies of the three drugs were

239

more than 90% using 10 mL of MeOH (A) or ACN (B) in 6 h. For eluent C, that is, 10 mL of

240

MeOH with 5% ammonia solution, a similar performance was reported when eluting these illicit

241

drugs from mixed reverse-phase cation exchange cartridges.6,7 The elution efficiencies we obtained

242

for KET and METH were 89.5% and 85.6%, respectively. The use of ammonia solution had a small

243

effect on the elution efficiency of KET and METH, but a significant decrease for AMP, from 97% to

244

56-72%. When ethyl acetate (D) was used, the elution efficiencies were very low (0.05), and an insignificant difference was also observed between the weekend and

304

weekday measurements (p>0.05). Of the three studied compounds in the influent, METH and AMP

305

were ubiquitous with high concentrations. As shown in Figure 2, METH was detected in all samples

306

with concentrations ranging from 101.7 ng·L−1 to 181.1 ng·L−1, which implied that the abuse of

307

METH is popular in Beijing. AMP was also detected in all influents, with concentrations at tens of

308

ng·L−1 levels. In contrast to a previous research in which no KET was detected in the WWTPs in

309

Beijing,23,27 the use of KET among Beijing inhabitants was found in the current study, with

310

concentrations at several ng·L−1 levels (Figure 2).

311

To estimate the drug consumption at the community level, the concentrations of the target

312

compounds in the influent of WWTP were measured and used. 5,7,27 The amount of the drugs (grams)

313

discharged per day was obtained by using the average concentrations of the selected drugs in the 15

ACS Paragon Plus Environment

Environmental Science & Technology 314

untreated wastewater samples and the daily flow rates in the WWTP. The amount of drugs excreted

315

in the wastewater each person each day was estimated by dividing the total amount of drugs in the

316

influent by the number of population served by the WWTP. The value is generally normalized as

317

grams per day per 1000 people. The consumption data of each drug in the community were obtained

318

by the fraction of the excreted drugs in the urine samples proposed by Zuccato et al. 5 The daily

319

influent load of a drug residue that enters a certain WWTP can be calculated with the following

320

equation:

mg ) 1000 ⋅ day

Influent load of a chemical residue (

321

322

=

Concentration of chemical residue(ng ⋅ L−1 ) × Influent flow( L ⋅ d −1 ) 1 mg × 6 ( ). Population served 10 ng 1000

323

The mass loads of METH in sewage influents ranged from 50.2 mg (1000 inh)−1 d−1 to 89.3 mg

324

(1000 inh)−1 d−1 (Table 1). The high level of METH in the influent loading was in accordance with

325

the seizure data for illicit drugs released by Chinese authority, that the number of seized METH in

326

China rose eightfold from 2007 to 2011.34 As the metabolite of METH, AMP has been detected in

327

the influents in other three WWTPs in Beijing.27 In this study, the estimated dose of AMP varied

328

from 19.6 mg (1000 inh)−1 d−1 to 26.8 mg (1000 inh)−1 d−1. The AMP traces in sewage influents with

329

concentrations no more than 10% those of METH concentrations were deemed to have come from

330

the metabolism of METH (other than AMP consumption itself).5 Thus, the AMP in the wastewater

331

influents from most cities in China came from METH degradation, which indicated that AMP

332

consumption can be neglected in China. This finding is in contrast to the results in European

333

countries, where more AMP was found compared with METH.35-37 In the sewage samples AMP may

334

potentially come from the illicit use of METH and AMP, clinical use of METH and AMP, and

335

clinical use of AMP precursors.38 In China, neither AMP nor METH are clinically used; therefore, 16

ACS Paragon Plus Environment

Page 16 of 30

Page 17 of 30

Environmental Science & Technology

336

the presence of AMP in the influents are likely resulted from the illicit uses of both METH and AMP,

337

and the clinical use of selegiline, which is the AMP precursor used in China.23,27

338

The loads of KET in this study ranged from 2.6 mg (1000 inh)−1 d−1 to 3.7 mg (1000 inh)−1 d−1

339

(Table 1). The highest mean KET loads were found on Sundays. The concentrations of KET were

340

higher during weekends (mean 7.2 ng L−1) than weekdays (mean 5.8 ng L−1), which implied the

341

high usage of KET during weekends. Several early studies also reported the increased use of KET

342

during weekends.24,35,39 The comparison of the three drug concentrations with the levels in the

343

literature is provided in the SI.

344

In Situ Profiling of Illicit drugs in River Water in Beijing

345

DGT devices were deployed in a river system for monitoring in January of 2016 to extend the

346

application of the developed technique in this study. No report on the illicit drugs in the urban river

347

waters of Beijing exists. Fifteen sampling sites in six urban rivers in Beijing were selected, and

348

DGT devices were deployed for 7 days. Over these 7 days, composite water samples were also

349

collected in the same way as described earlier. The water temperature during sampling was in the

350

range of 3.3 ± 0.7 °C (Figure S10). Field blanks were analyzed and were too low to be quantifiable.

351

The concentrations of the three illicit drugs were generally at the ng L−1 level in the river waters.

352

The average DGT concentrations of the three illicit drugs (Figures 3a and S11) were close to the

353

average concentrations measured by SPE method, which demonstrated that the average

354

concentrations measured by DGT can be substituted for the high-frequency active sampling. The

355

masses of drugs in the DGT sampler linearly increased during the 7-day deployment, which

356

indicated that the devices could be deployed for at least 7 days to obtain reliable and accurate

357

measurements (Figure 3b). The concentrations were relatively stable during the DGT deployment,

358

as shown by the grab sample data in Figure 3b. METH was detected by the SPE and DGT methods 17

ACS Paragon Plus Environment

Environmental Science & Technology 359

in all sampling sites, with the highest concentrations ranging from 4.37 ng L−1 to 57.36 ng L−1.

360

These values were higher than those of the METH detected in other rivers, such as Jarama and

361

Manzanares Rivers (3.1 ng L−1 to 5.0 ng L−1) in Spain,40 River Lambro (2.1 ng L−1) in Italy,41 and

362

Yellow River (1.3 ng L−1), Songhua River (2.2 ng L−1), and Yangze River (1.9 ng L−1) in China,48

363

and lower than those in Taiwan (n.d.-917 ng L−1) 14,29 and downstream (62.6 ng L−1) of a WWTP at

364

Hasting, Nebraska,11 but comparable with that in Pear River (31.1 ng L−1).42 The high

365

concentrations and detection frequencies of METH in the surface water in Beijing were consistent

366

with the seizure data from the Chinese authorities, which showed that METH and KET are the two

367

major drugs with increasing abuse dosages in China.4 However, the measurements in the Beijing

368

surface water in this study showed very low concentrations of KET (n.d.-3.85 ng L−1) throughout

369

the sampling survey, which is consistent with the effluent data from the Beijing sewage treatment

370

plants.23,27 It is the first application of well tested DGT method for measuring the three illicit drugs

371

in situ in a WWTP and in river waters with good performance and convincing results.

372

Environmental Significance

373

Grab sampling is the sampling method currently used most widely due to its easy use without

374

special requirement on equipment. Results obtained from the grab sampling however, may be

375

misleading because of the combination of the high variations in concentration and low in sampling

376

frequency. In the rivers and WWTPs where the point sources and rains may result in significant

377

fluctuations of water flow, grab sampling for water quality monitoring should be paid more attention.

378

Grab sampling is also likely to miss discharging event, such as the consumption of illicit drugs at

379

entertainment venues, such as pubs, bars, and nightclubs, at any given time when the sample

380

collection usually is not carried out.

381

Passive sampling, which can measure a time-weighted concentration of contaminants, has been 18

ACS Paragon Plus Environment

Page 18 of 30

Page 19 of 30

Environmental Science & Technology

382

developed and applied to deal with the above concern. POCIS has been used to collect

383

pharmaceutical and illicit drugs in water, which is able to overcome the drawbacks of traditional

384

sampling and provide time-averaged result.

385

semi-quantitative analysis as its sampling rate is affected by environmental conditions such as water

386

flow.11,15 DGT is superior over other passive sampling approaches, particularly in terms of its steady

387

accumulating rates, which are unaffected by the flow rate of water. The sampling rates (RS/A) for

388

METH and AMP by POCIS were 5.3 and 6.3 mL d−1 cm−2 with a surface area of ∼41 cm2, 11 and the

389

RS/A values by DGT were 7.86 and 7.10 mL d−1 cm−2, respectively. The higher sampling rates of

390

DGT suggest a shorter deployment time for achieving the same detection limit. This study

391

demonstrates that DGT is a better sampling technique than conventional active sampling methods

392

for investigation on the behaviors of illicit drugs in the water environment. It is the first time we

393

could confidently recommend a simple, robust and reliable method for illicit drug measurements

394

and monitoring to the government authority. Although DGT was only tested for measuring three

395

specific illicit drugs in this study, it is probably able to measure other illicit drugs, which requires

396

further investigation.

397

Supporting Information. Sampling sites in Beijing urban rivers, details on analysis methods,

398

including water sample pre-treatment, instrumental analysis, QA/QCs, uptake kinetics of illicit

399

drugs by XAD gels, elution efficiencies of illicit drugs by various eluents, DGT capacity,

400

relationship between DGT measurement of illicit drugs and the deployment time or the thickness of

401

diffusion gels, effect of solution pH and ionic strength on the DGT performance, temperature during

402

DGT deployment, competing effects, and physicochemical properties of illicit drugs.

403

AUTHOR INFORMATION

404

Corresponding Author

11,15,16

However, POCIS is only effective for

19

ACS Paragon Plus Environment

Environmental Science & Technology

Tel:

86-10-84915103;

Fax:

86-10-84926073.

E-mail:

Page 20 of 30

405

*

406

[email protected] (H.Z.)

407

Notes

408

The authors declare no competing financial interest.

409

ACKNOWLEDGMENT

410

This work was funded by the National Science Foundation of China (No. 41673120), Special Fund

411

for Basic Scientific Research of Central Public Research Institutes (2016YSKY-030), and Beijing

412

Natural Science Foundation (8173058).

413

REFERENCES

414

1. Lu, L.; Fang, Y.; Wang, X., Drug Abuse in China: Past, Present and Future. Cell Mol Neurobiol.

415

2008, 28, (4), 479-490.

416

2. Liu, Z.; Lian, Z. H. I.; Zhao, C., Drug use and HIV/AIDS in China. Drug Alcohol Rev. 2006, 25,

417

(2), 173-175.

418

3. UNODC, World drug report. New York: UNODC 2013.

419

4. Office of China National Narcotic Control Commission. Annual Report on Drug Control in

420

China (Beijing, China) 2015.

421

5. Zuccato, E.; Chiabrando, C.; Castiglioni, S.; Bagnati, R.; Fanelli, R., Estimating community

422

drug abuse by wastewater analysis. Environ. Health Perspect. 2008, 116, (8), 1027-32.

423

6. Du, P.; Li, K.; Li, J.; Xu, Z.; Fu, X.; Yang, J.; Zhang, H.; Li, X., Methamphetamine and

424

ketamine use in major Chinese cities, a nationwide reconnaissance through sewage-based

425

epidemiology. Water Res. 2015, 84, 76-84.

426

7. Subedi, B.; Kannan, K., Mass loading and removal of select illicit drugs in two wastewater

427

treatment plants in New York state and estimation of illicit drug usage in communities through 20

ACS Paragon Plus Environment

[email protected]

(X.J.),

Page 21 of 30

Environmental Science & Technology

428

wastewater analysis. Environ. Sci. Technol. 2014, 48, (12), 6661-6670.

429

8. Bones, J.; Thomas, K. V.; Paull, B., Using environmental analytical data to estimate levels of

430

community consumption of illicit drugs and abused pharmaceuticals. J. Environ. Monit. 2007, 9, (7),

431

701-707.

432

9. Fent, K.; Weston, A. A.; Caminada, D., Ecotoxicology of human pharmaceuticals. Aquat.

433

Toxicol. 2006, 76, (2), 122-159.

434

10. Pounds, N.; Maclean, S.; Webley, M.; Pascoe, D.; Hutchinson, T., Acute and chronic effects of

435

ibuprofen in the mollusc Planorbis carinatus (Gastropoda: Planorbidae). Ecotoxicol Environ Saf.

436

2008, 70, (1), 47-52.

437

11. Bartelt-Hunt, S. L.; Snow, D. D.; Damon, T.; Shockley, J.; Hoagland, K., The occurrence of

438

illicit and therapeutic pharmaceuticals in wastewater effluent and surface waters in Nebraska.

439

Environ. Pollut. 2009, 157, (3), 786-791.

440

12. Evgenidou, E. N.; Konstantinou, I. K.; Lambropoulou, D. A., Occurrence and removal of

441

transformation products of PPCPs and illicit drugs in wastewaters: A review. Sci. Total Environ.

442

2015, 505, 905-926.

443

13. Mastroianni, N.; Bleda, M. J.; López de Alda, M.; Barceló, D., Occurrence of drugs of abuse in

444

surface water from four Spanish river basins: Spatial and temporal variations and environmental risk

445

assessment. J. Hazard. Mater. 2016, 316, 134-142.

446

14. Jiang, J.-J.; Lee, C.-L.; Fang, M.-D.; Tu, B.-W.; Liang, Y.-J., Impacts of emerging contaminants

447

on surrounding aquatic environment from a youth festival. Environ. Sci. Technol. 2014, 49, (2),

448

792-799.

449

15. Fedorova, G.; Randak, T.; Golovko, O.; Kodes, V.; Grabicova, K.; Grabic, R., A passive

450

sampling method for detecting analgesics, psycholeptics, antidepressants and illicit drugs in aquatic 21

ACS Paragon Plus Environment

Environmental Science & Technology 451

environments in the Czech Republic. Sci. Total Environ. 2014, 487, 681-687.

452

16. Harman, C.; Reid, M.; Thomas, K. V., In situ calibration of a passive sampling device for

453

selected illicit drugs and their metabolites in wastewater, and subsequent year-long assessment of

454

community drug usage. Environ. Sci. Technol. 2011, 45, (13), 5676-5682.

455

17. Davison, W.; Zhang, H., In situ speciation measurements of trace components in natural waters

456

using thin-film gels. Nature 1994, 367, 546-548.

457

18. Garmo, Ø. A.; Røyset, O.; Steinnes, E.; Flaten, T. P., Performance study of diffusive gradients in

458

thin films for 55 elements. Anal. Chem. 2003, 75 (14), 3573-3580.

459

19. Challis, J. K.; Hanson, M. L.; Wong, C. S., Development and calibration of an organic-diffusive

460

gradients in thin films aquatic passive sampler for a diverse suite of polar organic contaminants.

461

Anal. Chem. 2016, 88 (21), 10583-10591.

462

20. Chen, C.-E.; Zhang, H.; Ying, G.-G.; Jones, K. C., Evidence and recommendations to support

463

the use of a novel passive water sampler to quantify antibiotics in wastewaters. Environ. Sci.

464

Technol. 2013, 47, (23), 13587-13593.

465

21. Chen, C.-E.; Chen, W.; Ying, G.-G.; Jones, K. C.; Zhang, H., In situ measurement of solution

466

concentrations and fluxes of sulfonamides and trimethoprim antibiotics in soils using o-DGT.

467

Talanta 2015, 132, 902-908.

468

22. Zheng, J.-L.; Guan, D.-X.; Luo, J.; Zhang, H.; Davison, W.; Cui, X.-Y.; Wang, L.-H.; Ma, L. Q.,

469

Activated charcoal based diffusive gradients in thin films for in situ monitoring of bisphenols in

470

waters. Anal. Chem. 2014, 87, (1), 801-807.

471

23. Li, J.; Hou, L.; Du, P.; Yang, J.; Li, K.; Xu, Z.; Wang, C.; Zhang, H.; Li, X., Estimation of

472

amphetamine and methamphetamine uses in Beijing through sewage-based analysis. Sci. Total

473

Environ. 2014, 490, 724-732. 22

ACS Paragon Plus Environment

Page 22 of 30

Page 23 of 30

Environmental Science & Technology

474

24. Castiglioni, S.; Borsotti, A.; Senta, I.; Zuccato, E., Wastewater analysis to monitor spatial and

475

temporal patterns of use of two synthetic recreational drugs, ketamine and mephedrone, in Italy.

476

Environ. Sci. Technol. 2015, 49, (9), 5563-5570.

477

25. Chen, C.-E.; Zhang, H.; Jones, K. C., A novel passive water sampler for in situ sampling of

478

antibiotics. J. Environ. Monit. 2012, 14, (6), 1523-1530.

479

26. Togola, A.; Budzinski, H., Development of polar organic integrative samplers for analysis of

480

pharmaceuticals in aquatic systems. Anal. Chem. 2007, 79, (17), 6734-6741.

481

27. Khan, U.; van Nuijs, A. L.; Li, J.; Maho, W.; Du, P.; Li, K.; Hou, L.; Zhang, J.; Meng, X.; Li, X.,

482

Application of a sewage-based approach to assess the use of ten illicit drugs in four Chinese

483

megacities. Sci. Total Environ. 2014, 487, 710-721.

484

28. Lai, F. Y.; Bruno, R.; Leung, H. W.; Thai, P. K.; Ort, C.; Carter, S.; Thompson, K.; Lam, P. K.;

485

Mueller, J. F., Estimating daily and diurnal variations of illicit drug use in Hong Kong: a pilot study

486

of using wastewater analysis in an Asian metropolitan city. Forensic Sci. Int. 2013, 233, (1),

487

126-132.

488

29. Lin, A. Y.-C.; Wang, X.-H.; Lin, C.-F., Impact of wastewaters and hospital effluents on the

489

occurrence of controlled substances in surface waters. Chemosphere 2010, 81, (5), 562-570.

490

30. Kim, K. Y.; Lai, F. Y.; Kim, H.-Y.; Thai, P. K.; Mueller, J. F.; Oh, J.-E., The first application of

491

wastewater-based drug epidemiology in five South Korean cities. Sci. Total Environ. 2015, 524,

492

440-446.

493

31. Metcalfe, C.; Tindale, K.; Li, H.; Rodayan, A.; Yargeau, V., Illicit drugs in Canadian municipal

494

wastewater and estimates of community drug use. Environ. Pollut. 2010, 158, (10), 3179-3185.

495

32. Bijlsma, L.; Emke, E.; Hernández, F.; De Voogt, P., Investigation of drugs of abuse and relevant

496

metabolites in Dutch sewage water by liquid chromatography coupled to high resolution mass 23

ACS Paragon Plus Environment

Environmental Science & Technology 497

spectrometry. Chemosphere 2012, 89, (11), 1399-1406.

498

33. Huerta-Fontela, M.; Galceran, M. T.; Martin-Alonso, J.; Ventura, F., Occurrence of psychoactive

499

stimulatory drugs in wastewaters in north-eastern Spain. Sci. Total Environ. 2008, 397, (1), 31-40.

500

34. UNODC, Patterns and trends of amphetamine-type stimulants and other drugs: Asia and the

501

Pacific 2012. Vienna: UNODC (Global SMART Programme) 2012.

502

35. Andrés-Costa, M. J.; Rubio-López, N.; Suárez-Varela, M. M.; Pico, Y., Occurrence and removal

503

of drugs of abuse in wastewater treatment plants of Valencia (Spain). Environ. Pollut. 2014, 194,

504

152-162.

505

36. Vuori, E.; Happonen, M.; Gergov, M.; Nenonen, T.; Järvinen, A.; Ketola, R. A.; Vahala, R.,

506

Wastewater analysis reveals regional variability in exposure to abused drugs and opioids in Finland.

507

Sci. Total Environ. 2014, 487, 688-695.

508

37. Östman, M.; Fick, J.; Näsström, E.; Lindberg, R. H., A snapshot of illicit drug use in Sweden

509

acquired through sewage water analysis. Sci. Total Environ. 2014, 472, 862-871.

510

38. Khan, U.; Nicell, J. A., Sewer epidemiology mass balances for assessing the illicit use of

511

methamphetamine, amphetamine and tetrahydrocannabinol. Sci. Total Environ. 2012, 421, 144-162.

512

39. Terzic, S.; Senta, I.; Ahel, M., Illicit drugs in wastewater of the city of Zagreb (Croatia) –

513

Estimation of drug abuse in a transition country. Environ. Pollut. 2010, 158, (8), 2686-2693.

514

40. Mendoza, A.; Rodríguez-Gil, J. L.; González-Alonso, S.; Mastroianni, N.; López de Alda, M.;

515

Barceló, D.; Valcárcel, Y., Drugs of abuse and benzodiazepines in the Madrid Region (Central

516

Spain): Seasonal variation in river waters, occurrence in tap water and potential environmental and

517

human risk. Environ. Int. 2014, 70, 76-87.

518

41. Zuccato, E.; Castiglioni, S.; Bagnati, R.; Chiabrando, C.; Grassi, P.; Fanelli, R., Illicit drugs, a

519

novel group of environmental contaminants. Water Res. 2008, 42 (4–5), 961-968. 24

ACS Paragon Plus Environment

Page 24 of 30

Page 25 of 30

Environmental Science & Technology

520

42. Li, K.; Du, P.; Xu, Z.; Gao, T.; Li, X., Occurrence of illicit drugs in surface waters in China.

521

Environ. Pollut. 2016, 213, 395-402.

522

25

ACS Paragon Plus Environment

Environmental Science & Technology

Page 26 of 30

Table 1. Median or mean influent loads and concentrations of METH, AMP and KET from WWTPs in China and elsewhere

523 Location

METH

AMP

KET

Year

References

Load,

Concentration,

Load,

Concentration,

Load,

Concentration,

mg (1000 inh)-1 day-1

ng L-1

mg (1000 inh)-1 day-1

ng L-1

mg (1000 inh)-1 day-1

ng L-1

65.42±15.0

132.48±30.42

21.94±2.62

44.44±5.31

3.27±0.37

6.63±0.74

2015

Shanghai

31.2-35.2

n.r.

4.8-5.2

n.r.

2.5-2.9

13

2012

27

Guangzhou

33.5-48.4

n.r.

3.8-40.5

n.r.

33.6

89

2012

27

Shenzhen

121.7

n.r.

21.9

n.r.

38.7

500

2012

27

Beijing

8.6-1014.6

n.r.