In recent years, mental disorders have become serious health problems. Monitoring the usage of psychotropic drugs can reflect the prevalence of mental disorders in the population, and wastewater analysis is an effective method for monitoring the use of psychoactive pharmaceuticals. In our study, we collected 58 samples from 12 wastewater treatment plants (WWTPs) in a southern city in China. We then analyzed the concentrations and loads of 18 frequently used psychoactive drugs and corresponding metabolites. The most abundant compounds in wastewater were diazepam, sulpiride, clozapine, venlafaxine, and quetiapine, with a detection frequency close to 100%. Sulpiride had the highest average load, followed by venlafaxine, clozapine, quetiapine, and diazepam. The higher load of diazepam in suburban WWTPs than in urban WWTPs reflects the large numbers of older people living in suburban areas in most of China. The lower loads of most target compounds in suburban WWTPs than in urban WWTPs indicate a greater likelihood of mental disorders in urban people than in suburban people. A correlation analysis of the target compounds revealed a correlation between the occurrence of venlafaxine, sulpiride, and clozapine. Overall, wastewater analysis is expected to be a significant tool for monitoring the consumption of psychoactive pharmaceuticals.

  • Timing and regional distribution of 18 psychoactive drugs and metabolites in a Chinese city.

  • The loads of most target compounds were lower in suburban WWTPs than in urban WWTPs.

  • Schools showed high loads of diazepam, indicating that diazepam may be abused among students.

  • Correlation was observed between venlafaxine, sulpiride, and clozapine loads.

In recent years, mental disorders have become a serious health problem and are receiving increasing attention. As a result, psychoactive drugs have become the most common prescription drugs globally, with over 30 billion doses prescribed daily in 2008 (Subedi & Kannan 2015). One meta-analysis study revealed that 25% of respondents had suffered from a mental disorder in the past year, and 29.2% had experienced a common mental disorder during their lifetimes (Steel et al. 2014). In China, which is one of the fastest developing countries, the production and use of psychoactive drugs are also increasing every year (Wang et al. 2015). The Chinese Epidemiological Survey on Mental Disorders reported the prevalence of mental disorders in China as 9.32% and the lifetime prevalence of schizophrenia and other mental disorders as 0.61%. These data, which confirm that mental disorders have become a major public health and social problem (Wittchen et al. 2011), emphasize the importance of monitoring and investigating substance dependence and mental disorders in the Chinese population.

Monitoring the use of psychotropic drugs can reflect the mental disorders in a population. In fact, most psychoactive drugs are metabolized in the human body and excreted either unchanged or as metabolites or conjugates (Calisto & Esteves 2009). The parent compounds and their metabolites then enter wastewater and are treated in wastewater treatment plants (WWTPs). This means that wastewater analysis is an effective method for monitoring psychoactive pharmaceutical use. Furthermore, wastewater analysis can also provide real-time, objective, and continuous information about and comparisons of the differences in psychoactive pharmaceutical consumption trends across time and space. Currently, wastewater analysis plays an important role in the monitoring of many psychoactive pharmaceuticals and illicit drugs, and it is currently undergoing comprehensive worldwide promotion.

Many studies have incorporated wastewater analysis to monitor the use of psychoactive drugs, but these have been conducted mainly in Europe, America, and Australia. Subedi & Kannan (2015) studied the occurrence of psychoactive drugs in wastewater in New York State, USA. Causanilles et al. (2017) researched the occurrence of pharmaceuticals in wastewater in Costa Rica. Kosjek et al. (2012) researched the occurrence of benzodiazepines in wastewater in Slovenia. Gurke et al. (2015) studied the occurrence of frequently prescribed pharmaceuticals and corresponding metabolites in wastewater in Germany. Several articles have been published regarding the occurrence of psychoactive drugs in wastewater in China (Zhou et al. 2010; Yuan et al. 2013; Sun et al. 2014; Sui et al. 2015; Wang et al. 2017). Nevertheless, most articles have focused only on Beijing (the capital of China), or only on a few psychoactive pharmaceuticals. In fact, the patterns of use of some psychoactive drugs differ from region to region (Du et al. 2019). In addition, information on regional and time variations in pharmaceuticals in wastewater is scarce.

In our study, we chose 18 frequently prescribed psychoactive drugs and their corresponding metabolites to investigate the occurrence of psychoactive pharmaceuticals and their corresponding metabolites in four suburban WWTPs, seven urban WWTPs, and one school WWTP on working days and weekends in a city in China. The selected psychoactive drugs included two antidepressants (venlafaxine and paroxetine), eight antipsychotic drugs (quetiapine, 7-hydroxyquetiapine, aripiprazole, dehydroaripiprazole, clozapine, risperidone, OH-risperidone, and sulpiride), and eight sedative-hypnotic drugs (promethazine, alprazolam, diazepam, nordiazepam, estazolam, flurazepam, midazolam, and hydroxymidazolam).

Chemicals and reagents

Promethazine, venlafaxine, quetiapine, 7-hydroxyquetiapine, alprazolam, aripiprazole, dehydroaripiprazole, clozapine, diazepam, nordiazepam, estazolam, flurazepam, midazolam, hydroxymidazolam, risperidone, OH-risperidone, paroxetine, sulpiride, and diazepam-d5 were purchased from Cerilliant (Round Rock, TX, USA). The mixed-stock standard solution (1 μg/mL) and internal standard solution (10 μg/mL) were diluted with methanol. Working standard solutions at concentrations of 100 and 5 ng/mL in methanol were prepared by serial dilutions of the stock standard solutions. A working internal standard solution at 20 ng/mL was generated by dilution in methanol.

High-performance liquid chromatography-grade methanol and acetonitrile were obtained from Chem Service (USA). Hydrochloric acid (analytical grade) was obtained from Yonghua Chemistry Co. (Jiangsu, China). Formic acid (98%) and ammonium acetate were obtained from Fluka (Buchs, Switzerland). An ammonium hydroxide solution (50%) was purchased from Aladdin Chemistry Co. (Shanghai, China). Ultrapure water was purified using a Milli-Q water system (Millipore, MA, USA). Glass microfiber filter papers (GF/C, 1.2 μm) were acquired from Whatman (Kent, UK).

Wastewater sample collection and treatment

The blank wastewater used for optimizing and validating the method was obtained from the inlet of a WWTP located in a small town in China. The blank sample was acidified with hydrochloric acid to pH 2, transported to the laboratory, and stored in the dark at −20 °C until analysis.

Authentic wastewater samples were collected every 2 hours from various WWTPs using automatic sampling devices. Twelve wastewater samples from the same day were combined and stored in 600 mL polyethylene terephthalate (PET) bottles as a 24 h composite sample. In our study, all 12 WWTPs were sampled for 4–6 days on working days and weekends (a total of 58 wastewater samples). In these WWTPs, Sj is located near some schools and mainly reflects the situation of the schools; ZHB, XS, SS, and CJZ are located in suburban districts; and DDH, KX, HC, WN, ZJB, XCP, and JXZ are located in the central part of the urban area. Further details about the sampling are included in Supplementary Table S1. The service population was provided by the corresponding WWTPs. After collection, the wastewater samples were acidified with hydrochloric acid to pH 2 and then transported to the laboratory for analysis. Samples were stored at −20 °C until pretreatment.

Sample pretreatment

Sample pretreatment followed the procedure described in a previous paper (Yuan et al. 2020), with minor modifications. Briefly, 50 mL of wastewater was filtered through glass microfiber filter papers, and an internal standard (100 μL of 20 ng/mL) was added to samples prior to solid-phase extraction (SPE). Oasis MCX cartridges were conditioned with methanol (4 mL) and acidified water (4 mL, pH 2 with hydrochloric acid). The wastewater samples were then loaded onto the cartridges at a flow rate of 2 mL/min. Before elution, the cartridges were washed with water (4 mL) and methanol (4 mL) at a flow rate of 2 mL/min and dried under vacuum. Analytes were eluted from the cartridges using 4 mL of 3% ammonium hydroxide in methanol at a flow rate of 1 mL/min. The eluate was then evaporated to dryness at 50 °C and redissolved in 200 μL mobile phase A/mobile phase B (60/40, v/v). Analyses were performed by injecting 20 μL into the ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS) system.

Instrumentation analysis

Chromatographic separation was performed using an AcquityTM Ultra Performance LC (Waters, USA) with an Allure PFPP column (100 mm × 2.1 mm, 5 μm) at a flow rate of 0.5 mL/min. Twenty mmol/L ammonium acetate, 5% acetonitrile, and 0.1% formic acid in water was mobile phase A. Acetonitrile was mobile phase B. The gradient was programed as follows: 0–1 min, from 40 to 45% B; 1–6.5 min, from 45 to 90% B; 6.5–12 min, 90%B; 12–12.5 min, from 90%B to 40%B; 12.5–15 min, 40%B.

Mass spectrometry analysis was performed using a Sciex 6500 Plus Q-trap quadrupole mass spectrometer (AB Sciex, USA). The analytes were analyzed in positive mode using multiple reaction monitoring (MRM) mode. The MRM parameters are shown in Supplementary Table S2. The data were evaluated using MultiQuant 3.0.2.

Method validation

Method validation included determinations of selectivity, linearity, limits of detection (LODs), lower limits of quantitation (LLOQs), extraction recovery, matrix effect, precision, and accuracy.

Selectivity was evaluated by analyzing six different blank wastewater samples. This evaluation was intended to demonstrate the potential interference of endogenous substances with the analyte or IS signals.

Linearity was determined by linear regression with 1/x weighting. The LODs and LLOQs were defined as the concentrations giving a signal-to-noise ratio (S/N) of at least 3 and 10, respectively. The LLOQs were also required to exhibit a relative standard deviation (RSD) of less than 20% and an accuracy ranging between 80 and 120%.

The method extraction recovery, matrix effect, accuracy, and precision were calculated by analyzing three concentrations of samples (5, 20, and 80 ng/L) with six replicates of each sample. The accuracy was expressed as the percentage ratio of the measured concentration to the nominal concentration. The precision was evaluated by the RSD. The extraction recovery was calculated by dividing the pre-extraction spiked sample areas by the post-extraction spiked sample areas. The matrix effect was calculated by dividing the post-extraction spiked sample areas by the standard solution.

Load estimation

The per capita load of compounds (Loadi,j) from each WWTP was estimated using the following equation:
where Ci,j is the residue concentrations of target compounds in WWTP j, Q is the influent flow of WWTPs, P is the population served by the specific WWTP, and n (1–12) is the number of WWTPs used in this study.
Averagei is the average load of the compounds in 12 WWTPs.

UPLC-MS/MS method optimization

To achieve the best peak shape, retention time, separation, and sensitivity, different high-performance liquid chromatography columns and different mobile phases were tested. Two different high-performance liquid chromatography columns, namely the Allure PFPP column (100 mm × 2.1 mm, 5 μm) and the Waters Acquity TM UPLC HSS T 3 column (100 mm × 2.1 mm, 1.8 μm), were evaluated. Different mobile phases, including water, methanol, or acetonitrile as solvents and the addition of formic acid, ammonium formate, acetic acid, and ammonium acetate, among others, as modifiers, were also evaluated. Ultimately, an Allure PFPP column and a mobile phase consisting of acetonitrile and water (20 mmol/L ammonium acetate, 5% acetonitrile, and 0.1% formic acid) were used. Figure 1 shows the chromatographic separation of the target compounds.
Figure 1

Total ion chromatogram of a mixture of target compounds in wastewater.

Figure 1

Total ion chromatogram of a mixture of target compounds in wastewater.

Close modal

Method validation

The results of method validation are shown in Supplementary Table S3. The correlation coefficients (r2) of linearity were all greater than 0.99. The LODs ranged from 0.2 to 2 ng/L and the LLOQs ranged from 0.5 to 5 ng/L. The accuracy and precision of the target compounds ranged from 88.2 to 111.6% and from 1.1 to 12.5%, respectively. The recovery and matrix effects of the target compounds ranged from 40.1 to 105.7% and from 60.3 to 99.7%, respectively.

Occurrence of targets in wastewater samples

The most abundant compounds in wastewater were diazepam, sulpiride, clozapine, venlafaxine, and quetiapine, which had a detection frequency close to 100%. The mean concentration of diazepam ranged from 1.0 ± 0.1 to 8.7 ± 0.6 ng/L, sulpiride from 9 ± 14.5 to 72.3 ± 14.0 ng/L, clozapine from 1.6 ± 3.6 to 24.5 ± 8.1 ng/L, venlafaxine from 3.2 ± 5.2 to 56.3 ± 18.6 ng/L, and quetiapine from 1.2 ± 1.5 to 17.8 ± 8.4 ng/L (Table 1). Nordiazepam is one of the metabolites of diazepam, but the nordiazepam concentration was below the LLOQ for most WWTPs (Table 1). The low frequency of nordiazepam may be related to its low excretion rate (Causanilles et al. 2017). By contrast, 7-hydroxyquetiapine, which is the major active metabolite of quetiapine (Hasselstrom & Linnet 2006), was detected in most WWTPs with mean concentrations ranging from 0.2 ± 0.3 to 1.3 ± 0.3 ng/L (Table 1). The detection frequency of estazolam was close to 75%, and its concentration was lower than 3 ng/L. Promethazine, alprazolam, and paroxetine had low detection frequencies. Other compounds, including aripiprazole, dehydroaripiprazole, flurazepam, midazolam, hydroxymidazolam, risperidone, and OH-Resperidone, were not detected in all WWTPs.

Table 1

Concentrations (ng/L) of target compounds at the WWTPs

WWTPAlprazolamClozapineDiazepamNordiazepamEstazolamVenlafaxinePromethazineParoxetineSulpirideQuetiapine7-hydroxyquetiapine
ZHB <LOD <LOD 5.3 <LOD <LOD <LOQ <LOD <LOD 1.7 <LOQ <LOD 
 < LOD < LOD 5.7 < LOD < LOD 1.3 < LOD < LOD 3.0 < LOQ < LOD 
 <LOD <LOD 6.4 <LOD <LOD 1.1 <LOD <LOD 2.9 <LOQ <LOD 
 <LOD <LOD 17.1 <LOQ <LOD <LOD <LOD <LOD 2.7 <LOQ <LOD 
 0.6 8.1 1.4 <LOD <LOD 12.5 <LOD <LOD 34.9 3.8 <LOD 
Mean 0.1 ± 0.3 1.6 ± 3.6 7.2 ± 5.9 <LOQ <LOD 3.2 ± 5.2 <LOD <LOD 9 ± 14.5 1.2 ± 1.5 <LOD 
HC <LOD 10.8 3.0 <LOD 1.4 24.7 <LOD <LOD 77.3 7.4 1.3 
 < LOD 20.3 2.2 < LOD 1.3 31.2 < LOD < LOD 51.5 7.3 0.8 
 <LOD 26.1 2.1 <LOD 1.9 37.7 <LOQ <LOD 53.5 10.4 1.5 
 <LOD 24.8 1.8 <LOD 1.6 39.3 <LOD <LOQ 51.0 10.7 1.0 
 <LOD 15.8 2.5 <LOD 1.3 33.0 <LOD <LOD 35.9 5.7 0.7 
Mean <LOD 19.6 ± 5.7 2.4 ± 0.4 <LOD 1.5 ± 0.2 33.2 ± 5.8 <LOD <LOD 53.8 ± 14.9 8.3 ± 2.2 1.1 ± 0.3 
JXZ <LOD 7.6 1.9 <LOD <LOQ 22.3 <LOD <LOD 20.5 11.8 0.9 
 < LOD 9.4 3.0 < LOD 1.0 22.4 < LOD < LOD 38.0 24.4 1.3 
 <LOD 20.9 3.0 <LOD 1.2 24.7 <LOD <LOD 55.3 23.3 1.1 
 <LOD 20.5 3.1 <LOD <LOQ 25.4 <LOD <LOD 61.1 23.6 1.6 
 <LOD 15.7 1.9 <LOD <LOQ 23.6 <LOD <LOD 33.0 6.0 0.8 
Mean <LOD 14.8 ± 6.2 2.6 ± 0.6 <LOD 1 ± 0.1 23.7 ± 1.4 <LOD <LOD 41.6 ± 16.6 17.8 ± 8.4 1.1 ± 0.3 
XCP <LOD 26.1 1.6 <LOD 1.9 30.9 <LOD <LOD 67.8 14.9 1.1 
 0.8 34.8 1.7 < LOD 1.8 46.2 2.2 4.1 70.3 15.8 1.5 
 0.9 42.1 1.6 <LOD 2.0 35.9 <LOD <LOD 51.5 15.7 1.1 
 0.8 45.7 1.6 <LOD 1.8 34.2 <LOD <LOD 46.0 15.2 1.8 
 <LOD 44.0 1.6 <LOD 1.7 42.5 <LOD <LOQ 99.5 14.1 1.1 
Mean 0.5 ± 0.5 24.5 ± 8.1 1.0 ± 0.1 <LOD 1.1 ± 0.1 37.9 ± 6.3 0.4 ± 1.0 1.6 ± 2.2 67.0 ± 20.9 15.2 ± 0.7 1.3 ± 0.3 
WN 0.8 24.0 4.3 <LOD 1.4 32.7 <LOD <LOD 37.0 9.0 0.7 
 0.8 19.6 5.8 < LOD 2.8 50.7 < LOD < LOQ 74.8 10.7 1.1 
 1.2 36.0 5.0 <LOD 1.9 69.9 <LOD 9.0 80.2 14.8 0.8 
 <LOD 12.1 2.0 <LOD <LOQ 13.1 <LOD <LOD 31.8 2.6 <LOD 
Mean 0.7 ± 0.5 22.9 ± 10.0 4.3 ± 1.6 <LOD 1.8 ± 0.8 41.6 ± 24.3 <LOD 3.3 ± 4.3 55.9 ± 25.1 9.3 ± 5.1 0.7 ± 0.5 
ZJB 0.8 16.7 2.8 <LOD 1.3 28.1 <LOD <LOQ 63.6 6.5 <LOD 
 0.9 19.3 2.4 < LOD 1.5 35.5 < LOD 6.9 72.2 9.7 0.6 
 1.3 15.3 2.5 <LOD 1.9 34.0 <LOD <LOD 75.7 8.6 0.9 
 1.3 23.5 2.5 <LOD 2.1 34.5 <LOQ <LOQ 93.4 16.9 1.7 
 1.1 24.0 3.9 <LOD 1.8 27.9 <LOD <LOD 56.5 10.5 0.5 
Mean 1.1 ± 0.2 19.8 ± 3.9 2.8 ± 0.6 <LOD 1.7 ± 0.3 32.0 ± 3.7 <LOD 3.0 ± 3.0 72.3 ± 14.0 10.4 ± 3.9 0.8 ± 0.6 
DDH <LOD 1.2 1.3 <LOD <LOQ 1.7 <LOD <LOD 8.8 <LOD <LOD 
 < LOD 4.4 1.6 < LOD < LOQ 3.4 < LOD < LOD 9.2 1.4 0.5 
 0.6 12.1 1.9 <LOD <LOQ 4.6 <LOD <LOD 28.2 1.6 <LOD 
 <LOD 16.9 2.3 <LOD 1.2 9.4 <LOD <LOD 38.6 3.8 0.6 
 <LOD 9.5 2.2 <LOD 1.5 3.5 <LOD <LOD 33.9 2.2 <LOD 
Mean 0.1 ± 0.3 8.8 ± 6.2 1.9 ± 0.4 <LOD 1.1 ± 0.2 4.5 ± 2.9 <LOD <LOD 23.7 ± 14.0 1.8 ± 1.4 0.2 ± 0.3 
KX <LOD 12.7 3.3 <LOD 1.1 20.8 <LOD <LOD 33.0 4.0 0.6 
 1.0 16.0 5.7 < LOD 1.6 41.7 < LOD < LOD 73.6 9.7 0.5 
 1.0 27.0 4.4 <LOD 1.6 45.2 <LOD <LOQ 58.7 11.4 1.2 
 0.6 9.8 7.5 <LOQ <LOQ 15.2 <LOD <LOD 50.2 3.6 0.5 
Mean 0.7 ± 0.5 16.4 ± 7.5 5.2 ± 1.8 <LOQ 1.3 ± 0.3 30.7 ± 14.9 <LOD <LOD 53.9 ± 16.9 7.2 ± 4.0 0.7 ± 0.3 
XS < LOD 18.6 7.4 < LOD < LOD 23.7 < LOD < LOD 24.2 4.0 0.7 
 <LOD 9.6 9.7 <LOQ <LOD 11.3 <LOD <LOD 20.0 3.2 <LOD 
 <LOD 11.6 7.9 <LOD <LOD 12.9 <LOD <LOD 12.9 4.2 0.7 
 <LOD 12.3 7.9 <LOD <LOD 14.8 <LOD <LOD 26.5 3.9 <LOQ 
 <LOD 9.1 6.2 <LOD <LOD 5.9 <LOD <LOD 16.3 2.3 <LOQ 
Mean <LOD 12.2 ± 3.8 7.8 ± 1.2 <LOQ <LOD 13.7 ± 6.5 <LOD <LOD 20.0 ± 5.6 3.5 ± 0.8 0.5 ± 0.3 
SJ <LOD 10.1 7.9 <LOD <LOQ 24.0 <LOD <LOD 12.8 4.9 0.6 
 <LOD 10.7 9.1 <LOQ <LOQ 29.8 <LOD <LOD 18.2 2.5 <LOQ 
 < LOD 17.8 10.7 < LOQ 1.5 24.9 < LOD < LOD 28.5 3.9 0.6 
 0.7 13.8 6.7 < LOD 1.1 26.5 < LOD < LOD 43.4 1.9 < LOQ 
 <LOD 11.4 8.0 <LOD <LOQ 16.0 <LOD <LOD 14.1 2.1 <LOQ 
 <LOD 14.1 9.4 <LOQ <LOQ 42.3 <LOD <LOD 45.4 3.2 <LOQ 
Mean 0.1 ± 0.3 13.0 ± 2.9 8.6 ± 1.4 <LOQ 1.1 ± 0.2 27.3 ± 8.7 <LOD <LOD 27.1 ± 14.5 3.1 ± 1.2 0.5 ± 0.1 
CJZ 0.6 21.0 8.7 2.9 <LOD 75.0 <LOD <LOD 65.3 7.9 1.5 
 0.6 25.5 8.6 2.3 < LOD 67.7 < LOD < LOD 45.1 8.9 0.8 
 <LOD 22.8 9.6 <LOQ <LOD 48.3 <LOD <LOD 40.5 6.9 1.4 
 <LOD 23.3 8.1 <LOQ <LOD 34.1 <LOD <LOD 46.2 5.9 1.0 
Mean 0.3 ± 0.3 23.2 ± 1.8 8.7 ± 0.6 2.3 ± 0.4 <LOD 56.3 ± 18.6 <LOD <LOD 49.3 ± 11.0 7.4 ± 1.3 1.2 ± 0.3 
SS <LOD 18.3 4.8 <LOD <LOQ 25.3 <LOD <LOD 78.2 6.1 <LOQ 
 < LOD 27.1 5.3 < LOD < LOQ 32.6 5.7 < LOD 60.5 7.0 1.2 
 <LOD 14.4 5.4 <LOD <LOQ 16.7 <LOD <LOD 60.0 3.6 0.9 
 <LOD 15.3 5.3 <LOD <LOQ 21.1 <LOD <LOD 55.9 3.8 0.9 
 0.7 18.9 7.4 <LOD 1.2 13.2 <LOD <LOD 65.8 4.0 0.8 
Mean 0.1 ± 0.3 18.8 ± 5.0 5.6 ± 1.0 <LOD 1.0 ± 0.1 21.8 ± 7.6 1.1 ± 2.5 <LOD 64.1 ± 8.6 4.9 ± 1.5 0.9 ± 0.3 
WWTPAlprazolamClozapineDiazepamNordiazepamEstazolamVenlafaxinePromethazineParoxetineSulpirideQuetiapine7-hydroxyquetiapine
ZHB <LOD <LOD 5.3 <LOD <LOD <LOQ <LOD <LOD 1.7 <LOQ <LOD 
 < LOD < LOD 5.7 < LOD < LOD 1.3 < LOD < LOD 3.0 < LOQ < LOD 
 <LOD <LOD 6.4 <LOD <LOD 1.1 <LOD <LOD 2.9 <LOQ <LOD 
 <LOD <LOD 17.1 <LOQ <LOD <LOD <LOD <LOD 2.7 <LOQ <LOD 
 0.6 8.1 1.4 <LOD <LOD 12.5 <LOD <LOD 34.9 3.8 <LOD 
Mean 0.1 ± 0.3 1.6 ± 3.6 7.2 ± 5.9 <LOQ <LOD 3.2 ± 5.2 <LOD <LOD 9 ± 14.5 1.2 ± 1.5 <LOD 
HC <LOD 10.8 3.0 <LOD 1.4 24.7 <LOD <LOD 77.3 7.4 1.3 
 < LOD 20.3 2.2 < LOD 1.3 31.2 < LOD < LOD 51.5 7.3 0.8 
 <LOD 26.1 2.1 <LOD 1.9 37.7 <LOQ <LOD 53.5 10.4 1.5 
 <LOD 24.8 1.8 <LOD 1.6 39.3 <LOD <LOQ 51.0 10.7 1.0 
 <LOD 15.8 2.5 <LOD 1.3 33.0 <LOD <LOD 35.9 5.7 0.7 
Mean <LOD 19.6 ± 5.7 2.4 ± 0.4 <LOD 1.5 ± 0.2 33.2 ± 5.8 <LOD <LOD 53.8 ± 14.9 8.3 ± 2.2 1.1 ± 0.3 
JXZ <LOD 7.6 1.9 <LOD <LOQ 22.3 <LOD <LOD 20.5 11.8 0.9 
 < LOD 9.4 3.0 < LOD 1.0 22.4 < LOD < LOD 38.0 24.4 1.3 
 <LOD 20.9 3.0 <LOD 1.2 24.7 <LOD <LOD 55.3 23.3 1.1 
 <LOD 20.5 3.1 <LOD <LOQ 25.4 <LOD <LOD 61.1 23.6 1.6 
 <LOD 15.7 1.9 <LOD <LOQ 23.6 <LOD <LOD 33.0 6.0 0.8 
Mean <LOD 14.8 ± 6.2 2.6 ± 0.6 <LOD 1 ± 0.1 23.7 ± 1.4 <LOD <LOD 41.6 ± 16.6 17.8 ± 8.4 1.1 ± 0.3 
XCP <LOD 26.1 1.6 <LOD 1.9 30.9 <LOD <LOD 67.8 14.9 1.1 
 0.8 34.8 1.7 < LOD 1.8 46.2 2.2 4.1 70.3 15.8 1.5 
 0.9 42.1 1.6 <LOD 2.0 35.9 <LOD <LOD 51.5 15.7 1.1 
 0.8 45.7 1.6 <LOD 1.8 34.2 <LOD <LOD 46.0 15.2 1.8 
 <LOD 44.0 1.6 <LOD 1.7 42.5 <LOD <LOQ 99.5 14.1 1.1 
Mean 0.5 ± 0.5 24.5 ± 8.1 1.0 ± 0.1 <LOD 1.1 ± 0.1 37.9 ± 6.3 0.4 ± 1.0 1.6 ± 2.2 67.0 ± 20.9 15.2 ± 0.7 1.3 ± 0.3 
WN 0.8 24.0 4.3 <LOD 1.4 32.7 <LOD <LOD 37.0 9.0 0.7 
 0.8 19.6 5.8 < LOD 2.8 50.7 < LOD < LOQ 74.8 10.7 1.1 
 1.2 36.0 5.0 <LOD 1.9 69.9 <LOD 9.0 80.2 14.8 0.8 
 <LOD 12.1 2.0 <LOD <LOQ 13.1 <LOD <LOD 31.8 2.6 <LOD 
Mean 0.7 ± 0.5 22.9 ± 10.0 4.3 ± 1.6 <LOD 1.8 ± 0.8 41.6 ± 24.3 <LOD 3.3 ± 4.3 55.9 ± 25.1 9.3 ± 5.1 0.7 ± 0.5 
ZJB 0.8 16.7 2.8 <LOD 1.3 28.1 <LOD <LOQ 63.6 6.5 <LOD 
 0.9 19.3 2.4 < LOD 1.5 35.5 < LOD 6.9 72.2 9.7 0.6 
 1.3 15.3 2.5 <LOD 1.9 34.0 <LOD <LOD 75.7 8.6 0.9 
 1.3 23.5 2.5 <LOD 2.1 34.5 <LOQ <LOQ 93.4 16.9 1.7 
 1.1 24.0 3.9 <LOD 1.8 27.9 <LOD <LOD 56.5 10.5 0.5 
Mean 1.1 ± 0.2 19.8 ± 3.9 2.8 ± 0.6 <LOD 1.7 ± 0.3 32.0 ± 3.7 <LOD 3.0 ± 3.0 72.3 ± 14.0 10.4 ± 3.9 0.8 ± 0.6 
DDH <LOD 1.2 1.3 <LOD <LOQ 1.7 <LOD <LOD 8.8 <LOD <LOD 
 < LOD 4.4 1.6 < LOD < LOQ 3.4 < LOD < LOD 9.2 1.4 0.5 
 0.6 12.1 1.9 <LOD <LOQ 4.6 <LOD <LOD 28.2 1.6 <LOD 
 <LOD 16.9 2.3 <LOD 1.2 9.4 <LOD <LOD 38.6 3.8 0.6 
 <LOD 9.5 2.2 <LOD 1.5 3.5 <LOD <LOD 33.9 2.2 <LOD 
Mean 0.1 ± 0.3 8.8 ± 6.2 1.9 ± 0.4 <LOD 1.1 ± 0.2 4.5 ± 2.9 <LOD <LOD 23.7 ± 14.0 1.8 ± 1.4 0.2 ± 0.3 
KX <LOD 12.7 3.3 <LOD 1.1 20.8 <LOD <LOD 33.0 4.0 0.6 
 1.0 16.0 5.7 < LOD 1.6 41.7 < LOD < LOD 73.6 9.7 0.5 
 1.0 27.0 4.4 <LOD 1.6 45.2 <LOD <LOQ 58.7 11.4 1.2 
 0.6 9.8 7.5 <LOQ <LOQ 15.2 <LOD <LOD 50.2 3.6 0.5 
Mean 0.7 ± 0.5 16.4 ± 7.5 5.2 ± 1.8 <LOQ 1.3 ± 0.3 30.7 ± 14.9 <LOD <LOD 53.9 ± 16.9 7.2 ± 4.0 0.7 ± 0.3 
XS < LOD 18.6 7.4 < LOD < LOD 23.7 < LOD < LOD 24.2 4.0 0.7 
 <LOD 9.6 9.7 <LOQ <LOD 11.3 <LOD <LOD 20.0 3.2 <LOD 
 <LOD 11.6 7.9 <LOD <LOD 12.9 <LOD <LOD 12.9 4.2 0.7 
 <LOD 12.3 7.9 <LOD <LOD 14.8 <LOD <LOD 26.5 3.9 <LOQ 
 <LOD 9.1 6.2 <LOD <LOD 5.9 <LOD <LOD 16.3 2.3 <LOQ 
Mean <LOD 12.2 ± 3.8 7.8 ± 1.2 <LOQ <LOD 13.7 ± 6.5 <LOD <LOD 20.0 ± 5.6 3.5 ± 0.8 0.5 ± 0.3 
SJ <LOD 10.1 7.9 <LOD <LOQ 24.0 <LOD <LOD 12.8 4.9 0.6 
 <LOD 10.7 9.1 <LOQ <LOQ 29.8 <LOD <LOD 18.2 2.5 <LOQ 
 < LOD 17.8 10.7 < LOQ 1.5 24.9 < LOD < LOD 28.5 3.9 0.6 
 0.7 13.8 6.7 < LOD 1.1 26.5 < LOD < LOD 43.4 1.9 < LOQ 
 <LOD 11.4 8.0 <LOD <LOQ 16.0 <LOD <LOD 14.1 2.1 <LOQ 
 <LOD 14.1 9.4 <LOQ <LOQ 42.3 <LOD <LOD 45.4 3.2 <LOQ 
Mean 0.1 ± 0.3 13.0 ± 2.9 8.6 ± 1.4 <LOQ 1.1 ± 0.2 27.3 ± 8.7 <LOD <LOD 27.1 ± 14.5 3.1 ± 1.2 0.5 ± 0.1 
CJZ 0.6 21.0 8.7 2.9 <LOD 75.0 <LOD <LOD 65.3 7.9 1.5 
 0.6 25.5 8.6 2.3 < LOD 67.7 < LOD < LOD 45.1 8.9 0.8 
 <LOD 22.8 9.6 <LOQ <LOD 48.3 <LOD <LOD 40.5 6.9 1.4 
 <LOD 23.3 8.1 <LOQ <LOD 34.1 <LOD <LOD 46.2 5.9 1.0 
Mean 0.3 ± 0.3 23.2 ± 1.8 8.7 ± 0.6 2.3 ± 0.4 <LOD 56.3 ± 18.6 <LOD <LOD 49.3 ± 11.0 7.4 ± 1.3 1.2 ± 0.3 
SS <LOD 18.3 4.8 <LOD <LOQ 25.3 <LOD <LOD 78.2 6.1 <LOQ 
 < LOD 27.1 5.3 < LOD < LOQ 32.6 5.7 < LOD 60.5 7.0 1.2 
 <LOD 14.4 5.4 <LOD <LOQ 16.7 <LOD <LOD 60.0 3.6 0.9 
 <LOD 15.3 5.3 <LOD <LOQ 21.1 <LOD <LOD 55.9 3.8 0.9 
 0.7 18.9 7.4 <LOD 1.2 13.2 <LOD <LOD 65.8 4.0 0.8 
Mean 0.1 ± 0.3 18.8 ± 5.0 5.6 ± 1.0 <LOD 1.0 ± 0.1 21.8 ± 7.6 1.1 ± 2.5 <LOD 64.1 ± 8.6 4.9 ± 1.5 0.9 ± 0.3 

Note. Sampling on weekends is in bold.

The mass loads of targets are shown in Table 2. Sulpiride had the highest average load (8.7 mg/1,000 inh/day), followed by venlafaxine (4.9 mg/1,000 inh/day), clozapine (3.6 mg/1,000 inh/day), quetiapine (1.5 mg/1,000 inh/day), and diazepam (0.8 mg/1,000 inh/day). In addition to being an antipsychotic, sulpiride has some effects on depression, and venlafaxine is also an antidepressant. These findings reflect the fact that depression has become one of the important components of mental disorders and has gradually become a major public health and social problem. The mean load of venlafaxine and quetiapine estimated in our study was lower than that reported in New York (Subedi & Kannan 2015). Diazepam loads reported previously in Beijing were similar to those of our study (Wang et al. 2017).

Table 2

Loads (mg/1,000 inh/day) of target compounds at the WWTPs

WWTPAlprazolam
Clozapine
Diazepam
Nordiazepam
Estazolam
Working daysWeekendsWorking daysWeekendsWorking daysWeekendsWorking daysWeekendsWorking daysWeekends
ZHB <0.1 <0.1 0.6(0.5a<0.2 1.1(1.0) 0.8 0.1(0.2) <0.3 <0.2 <0.2   
HC <0.1 <0.1 2.1(1.1) 1.6 0.2(0.1) 0.2 <0.2 <0.2 0.1(0.1) 0.1   
JXZ <0.1 <0.1 3.2(1.3) 2.2 0.5(0.1) 0.7 <0.5 <0.5 0.2(0.0) 0.2   
XCP 0.1(0.1) 0.1 5.5(1.6) 4.5 0.2(0.1) 0.2 <0.3 <0.3 0.3(0.1) 0.2   
WN 0.1(0.1) 0.1 2.4(1.0) 1.8 0.4(0.2) 0.5 <0.2 <0.2 0.1(0.0) 0.3   
ZJB 0.3(0.1) 0.2 5.4(1.0) 3.6 0.8(0.2) 0.4 <0.7 <0.7 0.5(0.1) 0.3   
DDH <0.1 <0.1 1.0(0.7) 0.6 0.2(0.1) 0.2 <0.2 <0.2 0.1(0.0) 0.1   
KX 0.1(0.1) 0.1 2.5(1.0) 1.6 0.8(0.3) 0.6 0.1(0.2) <0.2 0.2(0.0) 0.2   
XS <0.2 <0.2 3.2(0.6) 5.5 2.4(0.5) 2.2 0.2(0.3) <0.6 <0.3 <0.3   
SJ <0.1 0.1(0.1) 1.7(0.3) 2.4(0.4) 1.3(0.1) 1.3(0.4) 0.2(0.2) 0.2(0.2) 0.2(0.0) 0.2(0.0)   
CJZ <0.1 0.1 2.3(0.1) 3.0 0.9(0.1) 1.0 0.2(0.1) 0.3 <0.1 <0.1   
SS <0.1 <0.1 4.3(0.6) 7.0 1.5(0.3) 1.4 <0.5 <0.5 0.3(0.0) 0.3   
WWTPVenlafaxine
Promethazine
Paroxetine
Sulpiride
Quetiapine
7-hydroxyquetiapine
Working daysWeekendsWorking daysWeekendsWorking daysWeekendsWorking daysWeekendsWorking daysWeekendsWorking daysWeekends
ZHB 0.5(0.9) 0.2 <0.3 <0.3 <0.6 <0.6 1.6(2.4) 0.5 0.2(0.2) 0.1 <0.1 <0.1 
HC 3.5(1.2) 2.4 <0.2 <0.2 0.1(0.2) <0.3 5.6(2.0) 4.0 0.9(0.4) 0.6 0.1(0.1) 0.1 
JXZ 5.0(1.1) 5.2 <0.5 <0.5 <1.0 <1.0 8.6(4.1) 8.8 3.2(1.8) 5.6 0.2(0.1) 0.3 
XCP 4.9(1.0) 6.0 <0.3 0.3 0.1(0.2) 0.5 8.8(1.9) 9.1 2.1(0.6) 2.1 0.2(0.1) 0.2 
WN 3.8(2.9) 4.5 <0.2 <0.2 0.3(0.5) 0.4 4.9(2.6) 6.7 0.9(0.6) 1.0 0.1(0.0) 0.1 
ZJB 8.8(2.5) 6.6 0.1(0.3) <0.4 0.6(0.8) 1.3 20.6(6.9) 13.4 2.9(1.2) 1.8 0.2(0.2) 0.1 
DDH 0.5(0.3) 0.5 <0.2 <0.2 <0.5 <0.5 2.8(1.3) 1.2 0.2(0.2) 0.2 <0.1 0.1 
KX 4.0(1.7) 4.2 <0.4 <0.4 0.2(0.3) <0.4 7.2(0.6) 7.4 0.9(0.5) 1.0 0.1(0.0) 0.1 
XS 3.4(1.3) 7.0 <0.6 <0.6 <1.2 <1.2 5.6(1.7) 7.2 1.0(0.3) 1.2 0.1(0.1) 0.2 
SJ 4.2(1.7) 3.9(0.2) <0.3 <0.3 <0.6 <0.6 3.4(2.3) 5.4(1.6) 0.5(0.2) 0.4(0.2) <0.1 <0.1 
CJZ 5.6(2.8) 7.9 <0.2 <0.2 <0.5 <0.5 5.3(2.0) 5.2 0.7(0.2) 1.0 0.1(0.0) 0.1 
SS 4.9(1.4) 8.4 <0.5 1.5 <1.0 <1.0 16.7(2.5) 15.6 1.1(0.3) 1.8 0.2(0.1) 0.3 
WWTPAlprazolam
Clozapine
Diazepam
Nordiazepam
Estazolam
Working daysWeekendsWorking daysWeekendsWorking daysWeekendsWorking daysWeekendsWorking daysWeekends
ZHB <0.1 <0.1 0.6(0.5a<0.2 1.1(1.0) 0.8 0.1(0.2) <0.3 <0.2 <0.2   
HC <0.1 <0.1 2.1(1.1) 1.6 0.2(0.1) 0.2 <0.2 <0.2 0.1(0.1) 0.1   
JXZ <0.1 <0.1 3.2(1.3) 2.2 0.5(0.1) 0.7 <0.5 <0.5 0.2(0.0) 0.2   
XCP 0.1(0.1) 0.1 5.5(1.6) 4.5 0.2(0.1) 0.2 <0.3 <0.3 0.3(0.1) 0.2   
WN 0.1(0.1) 0.1 2.4(1.0) 1.8 0.4(0.2) 0.5 <0.2 <0.2 0.1(0.0) 0.3   
ZJB 0.3(0.1) 0.2 5.4(1.0) 3.6 0.8(0.2) 0.4 <0.7 <0.7 0.5(0.1) 0.3   
DDH <0.1 <0.1 1.0(0.7) 0.6 0.2(0.1) 0.2 <0.2 <0.2 0.1(0.0) 0.1   
KX 0.1(0.1) 0.1 2.5(1.0) 1.6 0.8(0.3) 0.6 0.1(0.2) <0.2 0.2(0.0) 0.2   
XS <0.2 <0.2 3.2(0.6) 5.5 2.4(0.5) 2.2 0.2(0.3) <0.6 <0.3 <0.3   
SJ <0.1 0.1(0.1) 1.7(0.3) 2.4(0.4) 1.3(0.1) 1.3(0.4) 0.2(0.2) 0.2(0.2) 0.2(0.0) 0.2(0.0)   
CJZ <0.1 0.1 2.3(0.1) 3.0 0.9(0.1) 1.0 0.2(0.1) 0.3 <0.1 <0.1   
SS <0.1 <0.1 4.3(0.6) 7.0 1.5(0.3) 1.4 <0.5 <0.5 0.3(0.0) 0.3   
WWTPVenlafaxine
Promethazine
Paroxetine
Sulpiride
Quetiapine
7-hydroxyquetiapine
Working daysWeekendsWorking daysWeekendsWorking daysWeekendsWorking daysWeekendsWorking daysWeekendsWorking daysWeekends
ZHB 0.5(0.9) 0.2 <0.3 <0.3 <0.6 <0.6 1.6(2.4) 0.5 0.2(0.2) 0.1 <0.1 <0.1 
HC 3.5(1.2) 2.4 <0.2 <0.2 0.1(0.2) <0.3 5.6(2.0) 4.0 0.9(0.4) 0.6 0.1(0.1) 0.1 
JXZ 5.0(1.1) 5.2 <0.5 <0.5 <1.0 <1.0 8.6(4.1) 8.8 3.2(1.8) 5.6 0.2(0.1) 0.3 
XCP 4.9(1.0) 6.0 <0.3 0.3 0.1(0.2) 0.5 8.8(1.9) 9.1 2.1(0.6) 2.1 0.2(0.1) 0.2 
WN 3.8(2.9) 4.5 <0.2 <0.2 0.3(0.5) 0.4 4.9(2.6) 6.7 0.9(0.6) 1.0 0.1(0.0) 0.1 
ZJB 8.8(2.5) 6.6 0.1(0.3) <0.4 0.6(0.8) 1.3 20.6(6.9) 13.4 2.9(1.2) 1.8 0.2(0.2) 0.1 
DDH 0.5(0.3) 0.5 <0.2 <0.2 <0.5 <0.5 2.8(1.3) 1.2 0.2(0.2) 0.2 <0.1 0.1 
KX 4.0(1.7) 4.2 <0.4 <0.4 0.2(0.3) <0.4 7.2(0.6) 7.4 0.9(0.5) 1.0 0.1(0.0) 0.1 
XS 3.4(1.3) 7.0 <0.6 <0.6 <1.2 <1.2 5.6(1.7) 7.2 1.0(0.3) 1.2 0.1(0.1) 0.2 
SJ 4.2(1.7) 3.9(0.2) <0.3 <0.3 <0.6 <0.6 3.4(2.3) 5.4(1.6) 0.5(0.2) 0.4(0.2) <0.1 <0.1 
CJZ 5.6(2.8) 7.9 <0.2 <0.2 <0.5 <0.5 5.3(2.0) 5.2 0.7(0.2) 1.0 0.1(0.0) 0.1 
SS 4.9(1.4) 8.4 <0.5 1.5 <1.0 <1.0 16.7(2.5) 15.6 1.1(0.3) 1.8 0.2(0.1) 0.3 

aStandard deviation.

In this study, a t-test was applied to compare drug loads between different regions (urban and suburban). The results showed that the load of most targets (alprazolam, diazepan, nordiazepam, estazolam, paroxetine, and quetiapine) had significant differences between urban and suburban (P < 0.05). As shown in Figure 2, diazepam showed a higher load in suburban WWTPs than in urban WWTPs. This may be related to the fact that large numbers of older people live in suburban areas in most regions of China. According to a recent survey, the use of hypnotic medications for insomnia and anxiety is higher in older individuals than in younger individuals (Smith & Tett 2009; Donoghue & Lader 2010). Additionally, SJ showed high loads of diazepam, indicating that diazepam may be abused among students, a finding confirmed by other studies (Li et al. 2005; Mancevska et al. 2008; Goreishi & Shajari 2013). Moreover, the loads of most target compounds with significant differences were lower in suburban WWTPs than in urban WWTPs, indicating that urban dwellers, rather than suburban dwellers, are more likely to suffer from mental disorders. This may be related to the fact that pressures from issues such as housing, education, and medical care are greater for urban people than for suburban people. No clear trend emerged with regard to average loads on the weekends relative to loads on the working days for any of the target compounds. ZHB and DDH showed lower loads for most targets than other WWTPs, which may be related to the surrounding environment of ZHB and DDH (more shopping malls, fewer residential areas, and more floating population).
Figure 2

Loads of diazepam (a), sulpiride (b), clozapine (c), venlafaxine (d), and quetiapine (e) in wastewater at the urban, suburban, and school WWTPs.

Figure 2

Loads of diazepam (a), sulpiride (b), clozapine (c), venlafaxine (d), and quetiapine (e) in wastewater at the urban, suburban, and school WWTPs.

Close modal

Correlations between targets

In this study, the correlations among targets were performed in Origin 2024 (OriginLab, USA) to explore the associations between patients with different mental illnesses. We found correlations among the venlafaxine, sulpiride, and clozapine loads (Figure 3). Venlafaxine is an antidepressant, whereas sulpiride and clozapine are important drugs for the treatment of schizophrenia. These results may therefore be associated with depression in patients with schizophrenia. The prevalence of depressive disorder in schizophrenia has been reported at around 40% (Upthegrove et al. 2017). In addition, the stage of illness (early vs. chronic) and state factors (acute or post-psychotic) influence these figures. In acute episodes, the rates of depression are up to 60%, while in post-psychotic schizophrenia, the rates of moderate to severe depression vary between 20% in chronic schizophrenia and 50% following treatment of the first episode (Upthegrove et al. 2010). In addition, some studies have found that the addition of sulpiride to clozapine resulted in clinical improvement in some patients with schizophrenia (Stubbs et al. 2000). Consequently, the observed correlations between venlafaxine, sulpiride, and clozapine loads are reasonable. According to the results of the study (Repo-Tiihonen et al. 2005), venlafaxine does not elevate clozapine plasma levels. If venlafaxine is combined with clozapine, special attention must be paid to the patient's clinical state because it is theoretically possible, due to the elevation of noradrenalin caused by both these agents, to produce blood pressure elevation. It is recommended that clozapine plasma levels and blood pressure be monitored.
Figure 3

Correlations between the loads of venlafaxine, sulpiride, and clozapine.

Figure 3

Correlations between the loads of venlafaxine, sulpiride, and clozapine.

Close modal

Monitoring the use of psychotropic drugs can reflect the mental disorders in a population. Currently, wastewater analysis is considered an effective method for monitoring the use of psychoactive pharmaceuticals. This study investigated the occurrence of 18 psychoactive drugs and their corresponding metabolites in 12 WWTPs in a southern city in China. For most compounds, lower loads of target compounds were found in the WWTPs in the suburbs than in the urban areas. On the contrary, diazepam showed higher loads in suburban WWTPs than in urban WWTPs, and school areas showed high loads of diazepam. In addition, a correlation analysis of the targets revealed a correlation between venlafaxine, sulpiride, and clozapine. These results may be associated with depression in patients with schizophrenia. Overall, wastewater analysis is expected to be a significant tool for monitoring the consumption of psychoactive pharmaceuticals.

The authors are grateful to the Shanghai Science and Technology Commission (24DZ3001800), the Opening Foundation of the Key Laboratory of Forensic Medicine of the Ministry of Justice (GY2025C-3), the Shanghai Forensic Service Platform (19DZ2292700), the Shanghai Key Laboratory of Forensic Medicine (21DZ2270800) and the Shanghai Association of Forensic Science (SHSFJD2023-009) for their financial support of this study.

All relevant data are included in the paper or its Supplementary Information.

The authors declare there is no conflict.

Calisto
V.
&
Esteves
V. I.
(
2009
)
Psychiatric pharmaceuticals in the environment
,
Chemosphere
,
77
(
10
),
1257
1274
.
Causanilles
A.
,
Ruepert
C.
,
Ibáñez
M.
,
Emke
E.
,
Hernández
F.
&
de Voogt
P.
(
2017
)
Occurrence and fate of illicit drugs and pharmaceuticals in wastewater from two wastewater treatment plants in Costa Rica
,
Science of the Total Environment
,
599
,
98
107
.
Donoghue
J.
&
Lader
M.
(
2010
)
Usage of benzodiazepines: a review
,
International Journal of Psychiatry in Clinical Practice
,
14
(
2
),
78
87
.
Du
P.
,
Thai
P. K.
,
Bai
Y.
,
Zhou
Z. L.
,
Xu
Z. Q.
,
Zhang
X.
,
Wang
J.
,
Zhang
C. B.
,
Hao
F. H.
&
Li
X. Q.
(
2019
)
Monitoring consumption of methadone and heroin in major Chinese cities by wastewater-based epidemiology
,
Drug and Alcohol Dependence
,
205
,
7
.
Goreishi
A.
&
Shajari
Z.
(
2013
)
Substance abuse among students of Zanjan's universities (Iran): a knot of today's society
,
Addiction & Health
,
5
(
1–2
),
66
72
.
Gurke
R.
,
Rössler
M.
,
Marx
C.
,
Diamond
S.
,
Schubert
S.
,
Oertel
R.
&
Fauler
J.
(
2015
)
Occurrence and removal of frequently prescribed pharmaceuticals and corresponding metabolites in wastewater of a sewage treatment plant
,
Science of The Total Environment
,
532
,
762
770
.
Hasselstrom
J.
&
Linnet
K.
(
2006
)
In vitro studies on quetiapine metabolism using the substrate depletion approach with focus on drug-drug interactions
,
Drug Metabolism and Drug Interactions
,
21
(
3–4
),
187
211
.
Kosjek
T.
,
Perko
S.
,
Zupanc
M.
,
Hren
M. Z.
,
Dragicevic
T. L.
,
Zigon
D.
,
Kompare
B.
&
Heath
E.
(
2012
)
Environmental occurrence, fate and transformation of benzodiazepines in water treatment
,
Water Research
,
46
(
2
),
355
368
.
Li
K.-Q.
,
Zhang
J.-J.
,
Lu
B.-l.
,
Liu
L.-Z.
&
Yin
P.
(
2005
)
Study on the use of benzodiazepine among middle school students in Wuhan City
,
Zhonghua liu xing bing xue za zhi=Zhonghua liuxingbingxue zazhi
,
26
(
8
),
596
599
.
Mancevska
S.
,
Bozinovska
L.
,
Tecce
J.
,
Pluncevik-Ghgoroska
J.
&
Sivevska-Smilevska
E.
(
2008
)
Depression, anxiety and substance use in medical students in the Republic of Macedonia
,
Bratislava Medical Journal-Bratislavske Lekarske Listy
,
109
(
12
),
568
572
.
Repo-Tiihonen
E.
,
Eloranta
A.
,
Hallikainen
T.
&
Tiihonen
J.
(
2005
)
Effects of venlafaxine treatment on clozapine plasma levels in schizophrenic patients
,
Neuropsychobiology
,
51
(
4
),
173
176
.
Steel
Z.
,
Marnane
C.
,
Iranpour
C.
,
Chey
T.
,
Jackson
J.
,
Patel
V.
&
Silove
D.
(
2014
)
The global prevalence of common mental disorders: a systematic review and meta-analysis 1980–2013
,
International Journal of Epidemiology
,
43
(
2
),
476
493
.
Stubbs
J. H.
,
Haw
C. M.
,
Staley
C. J.
&
Mountjoy
C. Q.
(
2000
)
Augmentation with sulpiride for a schizophrenic patient partially responsive to clozapine
,
Acta Psychiatrica Scandinavica
,
102
(
5
),
390
393
.
Sui
Q.
,
Wang
D.
,
Zhao
W. T.
,
Huang
J.
,
Yu
G.
,
Cao
X. Q.
,
Qiu
Z. F.
&
Lu
S. G.
(
2015
)
Pharmaceuticals and consumer products in four wastewater treatment plants in urban and suburb areas of Shanghai
,
Environmental Science and Pollution Research
,
22
(
8
),
6086
6094
.
Upthegrove
R.
,
Birchwood
M.
,
Ross
K.
,
Brunett
K.
,
McCollum
R.
&
Jones
L.
(
2010
)
The evolution of depression and suicidality in first episode psychosis
,
Acta Psychiatrica Scandinavica
,
122
(
3
),
211
218
.
Upthegrove
R.
,
Marwaha
S.
&
Birchwood
M.
(
2017
)
Depression and Schizophrenia: cause, consequence, or trans-diagnostic issue?
,
Schizophrenia Bulletin
,
43
(
2
),
240
244
.
Wang
C. C.
,
Hou
L. L.
,
Li
J.
,
Xu
Z. Q.
,
Gao
T. T.
,
Yang
J.
,
Zhang
H. F.
,
Li
X. Q.
&
Du
P.
(
2017
)
Occurrence of diazepam and its metabolites in wastewater and surface waters in Beijing
,
Environmental Science and Pollution Research
,
24
(
18
),
15379
15389
.
Wittchen
H. U.
,
Jacobi
F.
,
Rehm
J.
,
Gustavsson
A.
,
Svensson
M.
,
Jönsson
B.
,
Olesen
J.
,
Allgulander
C.
,
Alonso
J.
,
Faravelli
C.
,
Fratiglioni
L.
,
Jennum
P.
,
Lieb
R.
,
Maercker
A.
,
van Os
J.
,
Preisig
M.
,
Salvador-Carulla
L.
,
Simon
R.
&
Steinhausen
H. C.
(
2011
)
The size and burden of mental disorders and other disorders of the brain in Europe 2010
,
European Neuropsychopharmacology
,
21
(
9
),
655
679
.
Yuan
S.
,
Wang
X.
,
Wang
R. J.
,
Luo
R. X.
,
Shi
Y.
,
Shen
B. H.
,
Liu
W.
,
Yu
Z. G.
&
Xiang
P.
(
2020
)
Simultaneous determination of 11 illicit drugs and metabolites in wastewater by UPLC-MS/MS
,
Water Science and Technology
,
82
(
9
),
1771
1780
.
Zhou
H. D.
,
Wu
C. Y.
,
Huang
X.
,
Gao
M. J.
,
Wen
X. H.
,
Tsuno
H.
&
Tanaka
H.
(
2010
)
Occurrence of selected pharmaceuticals and caffeine in sewage treatment plants and receiving rivers in Beijing, China
,
Water Environment Research
,
82
(
11
),
2239
2248
.

Author notes

These two authors contributed equally to this work.

This is an Open Access article distributed under the terms of the Creative Commons Attribution Licence (CC BY-ND 4.0), which permits copying and redistribution with no derivatives, provided the original work is properly cited (http://creativecommons.org/licenses/by-nd/4.0/).

Supplementary data