Ten Years-Snapshot of the Occurrence of Emerging Contaminants in Drinking , Surface and Ground Waters and Wastewaters from São Paulo State , Brazil

Emerging contaminants have been considered one of the main concerns for ensuring the quality of water around the world. This work presents the results of 10 years of analyses carried out in the state of São Paulo (Brazil) that has the high population density and intense agricultural and industrial activities. In this work 58 compounds (9 hormones, 14 pharmaceuticals and personal care products, 8 industrial compounds, 17 pesticides and 10 illicit drugs) were determined from 2006 to 2015 in 708 samples including raw and treated sewage, surface and ground and drinking waters. A preliminary risk assessment for aquatic life protection identified potential risks for caffeine, paracetamol, diclofenac, 17α-ethynylestradiol, 17β-estradiol, estriol, estrone, testosterone, triclosan, 4-n-nonylphenol, bisphenol A, atrazine, azoxystrobin, carbendazim, fipronil, imidacloprid, malathion and tebuconazole. Drinking water criteria were available only for 22 compounds and for them no adverse effects were expected at the concentrations found, except for 17β-estradiol.


Introduction
A large number and quantity of substances are produced and consumed around the world to improve human life quality in a variety of activities, such as the protection of crops, the treatment of diseases, which has increased people's life expectance, or just to improve the comfort of everyday life. Residues of these substances reaches the environment through different pathways contaminating soils, waters and the atmosphere. Some of them are classified as emerging contaminants, i.e., natural or anthropogenic substances of emerging concern that are not commonly monitored but are found in the environment and may cause known or suspected adverse effects in the humans or wildlife. 1 Pesticides, pharmaceuticals and personal care products, hormones, sunscreen/UV filters, illicit drugs, perfluorinated compounds, disinfection by-products, nanomaterials, microplastics, are among the most investigated substances. 2 In fact, innumerous contaminants of emerging concern are being increasingly investigated around the world since pioneers works revealed negative effects on the reproduction of fish related to the presence of natural and synthetic hormones coming from wastewater treatment plants (WWTPs) discharges. 3,4 Although some emerging contaminants may demonstrate low acute toxicity, many of them are able to cause significant effects on the metabolism of living organisms at very low concentrations, i.e., around nanogram or microgram per liter, in a chronic exposure scenario and can affect the health of aquatic organisms as well as humans. 1 Over the past recent years, several studies have shown that these contaminants are not necessarily new molecules,

Ten Years-Snapshot of the Occurrence of Emerging Contaminants in Drinking, Surface and Ground Waters and Wastewaters from São Paulo State, Brazil
Cassiana C. Montagner but rather substances that were recently detected due to the improvement of analytical instrumentation that allowed the development of multi-residue methods that are capable of measure low concentrations of these contaminants in different water matrices such as surface, ground and drinking waters. 2 Occurrence data has been obtained over the years, which allowed the assessment of quality of waters by risk assessment approaches. Monitoring studies carried out in different countries [5][6][7][8][9][10] provide important information to the establishment of which compounds are the most relevant in each region, according to their social and economic characteristics. The quality and quantity of available data on the occurrence and effects of such contaminants have allowed the elucidation of contamination scenarios, as well as proposal of priority list of compounds considering their toxicity, concentration and/or frequency of detection in each region. This important step should be done before starting a regulation process, when it is necessary to choose few indicators in a list that includes thousands of non-regulated contaminants present in the environment. 11 This work was carried out in São Paulo state (Brazil), which has the highest density population in the country and 41 million inhabitants. It is also the most developed State, responsible for almost 30% of Brazilian gross domestic product (GDP). The most demanding activity for water is irrigation, followed by industrial, livestock and drinking water. 12 In a recent review about emerging contaminants in different Brazilian aquatic matrices, a list of 56 papers published between 1997 and 2017 reveal the occurrence of around two hundred compounds in waste, surface and drinking waters from 11 of the 26 states and Federal District. 7 The Southeastern region presented the largest number of studies 7 especially in the state of São Paulo, which has high patterns of consumption, similar to developed countries. However, this region also accumulates environmental problems that are typical to developing countries, especially related to the poor sanitation, where not all waste waters are properly treated before being released to surface waters. Depending on the level of contamination, several uses of water have been impaired. 13 In this paper occurrence data of 58 emerging contaminants of the 708 samples of drinking, surface, ground and waste waters (sampled in the state of São Paulo, Brazil, between 2006 and 2015) were collected to perform a preliminary risk assessment for two of the multiple uses of water, i.e., aquatic life protection and drinking water. The data were retrieved from different researches developed in the Environmental Chemistry Lab from the State University of Campinas, including published papers and unpublished data that were produced using the analytical methods developed in the research group.
Solid phase extraction (SPE) cartridges OASIS HLB 6cc, 500 mg (Waters, Milford, USA) were used for all compounds. Anionic-exchange cartridges 500 mg Strata SAX (Phenomenex, Torrance, USA) were used to remove humid substances before antimicrobials extraction. Glass fiber filters, used to preliminary filtration of the samples, were purchased from Sartorius Stedim Biotech (Goettingen, Germany).
Individual 400 mg L -1 stock solutions were prepared by dissolving appropriate amounts of each standard in methanol and kept in amber glass bottles at 4 ºC. Working solutions in the range from 0.005 to 10 mg L -1 were prepared by diluting appropriate aliquots of the stock solutions. Quantification was performed using external calibration according to Miller and Miller. 14

Study site and sample collection
Occurrence data considered a total of 708 samples, including raw and treated wastewaters, surface, ground and drinking waters, collected in the state of São Paulo at different locations and times between 2006 and 2015 ( Figure 1). Sampling locations were selected to represent urban areas with variable population densities and industrial activities. In addition, rural areas with different crop production were also included in this study. Details on sampling sites can be found in the Supplementary Information (SI) section.
Water samples were collected in amber glass bottles previously decontaminated with a detergent solution, tap and distilled water, and then rinsed with ethanol and acetone before baking at 400 ºC for 4 h. All containers were wrapped with aluminum foil. Before sampling, bottles were firstly rinsed and then, completed with the sampled water. Then, they were sealed and transported in a thermal cooler to the laboratory. Samples were kept refrigerated (4 ºC) until the extraction procedures within 48 h.

Analytical methods
The number of investigated compounds at Environmental Chemistry Lab has continuously increased over the years. Therefore, not all compounds were analyzed for all samples. Fifty eight emerging contaminants including 9 hormones, 13 pharmaceuticals and 1 personal care products, 8 industrial compounds, 17 pesticides and 10 illicit drugs, were quantified by liquid chromatography coupled to triple quadrupole mass spectrometry (LC(ESI) MS/MS) [15][16][17][18][19] (Agilent Technologies, USA) or high performance liquid chromatography coupled to ultravioletdiode array (HPLC-UV/DAD) or molecular fluorescence (HPLC-FLD) detectors (Shimadzu Corporation, Japan). 20 The methods used to quantify the emerging contaminants were already published. [15][16][17][18][19][20] Details about the methodology used for each sample are described in the SI section.

Preliminary risk assessment
Water quality criteria (WQC) were selected from peer-reviewed literature and/or from webpages of official environmental and health agencies around the world. Risk evaluation was performed in a worst-case scenario comparing the lowest water quality criteria found with the maximum environmental concentrations (MEC). The risk quotient method was employed dividing MEC by WQC. Values greater than one indicate that risk for the specific use of the water is expected. 21 Water quality criteria for protection of aquatic life were used to assess the risk of surface waters. Drinking water quality criteria were used to assess the risk of water for human consumption.

Surface water
A total of three hundred and twenty-nine (329) samples of surface waters were collected in the state of São Paulo between 2006 and 2015. The number of samples analyzed, the frequency of positive samples and the concentration of each compound are presented in Table 1 and in the SI section. WQC for aquatic life protection were found for 33 of the 58 substances investigated and the selected criteria and the risk quotient of the worst scenario were also presented in Table 1. Figure 2 shows a plot representation of the concentrations and the preliminary risk assessment.
Finally, other 6 substances commonly associated with their industrial origins, which are also notably found in consumer goods, were investigated, including bisphenol A that was quantified in 145 of the 217 samples in a wide range of concentrations, between 2 and 13016 ng L -1 . Among these substances, dibutylphthalate was investigated in 36 samples and was found in 34, corresponding to a frequency of detection 94%. Dioctylphthalate was found in the two samples, but other phthalates have never been found. The occurrence data is described below and classified by compound class.

Pesticides
Pesticide contamination of surface waters has been well documented worldwide as it represents a major issue when aspects related to water quality and uses are concerned. In this study, a wide range of concentrations was noticed for pesticides ( Figure 2) being the fungicide carbendazim, by far, the most frequent as well as the most abundant substance. Considering only positive data, the concentration of carbendazim varied from 0.8 to 4520 ng L -1 , reflecting different characteristics of the sampling points. This broadspectrum fungicide is widely used in different crops and forests, but also as a household biocide. 35,36 Palma et al. 37 verified surface water contamination with carbendazim at levels of up to 1200 ng L -1 in an agricultural-forestry basin in Southern Chile. However, the authors could not explain the origin of contamination since neither rain or pesticide application were registered during the sampling period. Lower concentrations of carbendazim, in comparison to our data, were also observed by Masiá et al. 38 in water samples of the Guadalquivir River Basin, an agriculture-based watershed in Spain. According to the authors, this fungicide presented maximum concentrations of 11 ng L -1 , with a frequency of detection of 17%. Besides the agriculture use, one can consider the contamination of urban waters due to carbendazim use as a household biocide. Burkhardt et al. 39 investigated the leaching of pesticides used in plants and building materials to sewer systems and concluded that, after the first flush, concentrations of the investigated substances, including carbendazim, exceeded the Swiss water quality standard. Chen et al. 40 investigated different biocides in aquatic environments of a highly urbanized region in China and also concluded that domestic sewage was the dominant pollution source for most biocides.
Triazoles, also used as broad-spectrum fungicides, were the second most abundant group of pesticides. In this group, tebuconazole was found in 31% of the samples and presented an average concentration of 39 ng L -1 . Other triazole fungicides presented average and maximum concentrations at the same order of magnitude. However, a maximum concentration of 1071 ng L -1 was noticed for tebuconazole in this work. Contamination of surface waters with tebuconazole by runoff from cropped areas is well documented, 41,42 but this fungicide is also used as a household biocide making it also commonly found in wastewaters and, consequently, in urban affected surface waters. 43,44 An average concentration up to 35 ng L -1 were estimated by Gerónimo et al. 45 in Argentinean agricultural river basins, whereas concentrations up to 15 ng L -1 were reported by da Silva et al. 46 investigating surface waters nearby rice crops in the southern region of Brazil.
Fungicides of the β-methoxyacrylates group were also investigated in this work and presented average concentrations varying from 3 to 13 ng L -1 . The most abundant fungicide of this group was azoxystrobin, which presented a maximum concentration of 431 ng L -1 . In the study carried out by Battaglin et al., 47 azoxystrobin was the most frequently detected fungicide (45%) in U.S. rivers with a maximum concentration of 1130 ng L -1 . According to the authors, its widespread occurrence was associated with its increasing use in a variety of agricultural and nonagricultural settings in the United States.
Atrazine, the second most frequently detected pesticide in this work, presented a maximum concentration of 611 ng L -1 . Although this triazine herbicide is included in several aquatic life standards, its presence in surface waters has an emerging concern due to possible endocrine disrupting effects at low concentrations. 48 Since the occurrence of atrazine is commonly associated to agricultural activities, 49,50 its concentrations are strongly dependent on the watershed characteristics. Steady concentrations around 500 ng L -1 has been noticed in the Mississipi River and its major tributaries during the past three decades suggesting that atrazine use in the area has not diminished. 51,52 Concentrations around 17 ng L -1 were detected by Moreno-González et al. 53 in a seasonal investigation carried out in a Mediterranean coastal lagoon. Similar levels were also found by Masiá et al. 54 in the Guadalquivir River Basin, in Spain. Lower concentrations in European samples may be a consequence of atrazine banishment in the European Union. 55 Other studies carried out in Brazil have found concentrations similar to those of the present work. 56,57 Other triazines were also investigated in this work, with emphasis to ametryn, which presented a maximum concentration of 148 ng L -1 in the investigated samples. The herbicides clomazone and bromacil were also quantified at similar concentrations, but the latter was detected only in 2% of the samples.
The insecticides imidacloprid, fipronil and malathion were also detected in the samples, being the former the most frequent in the samples. Average concentration for the insecticides varied between 10 and 26 ng L -1 .

Pharmaceuticals and personal care products
A wide range of concentrations was observed for pharmaceuticals ( Figure 2) with relatively high concentrations noticed for acetaminophen, acetylsalicylic acid and caffeine. Caffeine levels in this work varied from 19 to 127000 ng L -1 , reflecting the differences between sampling points according to their anthropic impact. Kolpin et al. 58 investigated more than 70 water samples collected upstream and downstream of selected urban areas from Iowa (USA) for caffeine analysis and other organic wastewater contaminants. Caffeine levels for low-flow streams were the highest amidst other flow conditions reaching a maximum of 294 ng L -1 , i.e., almost 430 times lower than the maximum concentration reported in this work. Ide et al. 59 investigated the presence of caffeine along the Iguaçu River, in Southern Brazil, and verified concentrations up to 5170 ng L -1 due to a large human influence in the upper part of the river.
Average concentrations around 6000 ng L -1 were noticed for acetylsalicylic acid and acetaminophen, two of the most consumed nonprescription pharmaceuticals in Brazil. Diclofenac was detected at an average concentration of 106 ng L -1 , which is a little higher than the annual average environmental quality standard (AA EQS = 100 ng L -1 ) proposed by European Union. 6 After a systematic literature review that included articles over the period from 1995 to 2015 about the occurrence of diclofenac in European surface waters, Tiedeken et al. 6 concluded that in general the concentrations were below AA EQS. Besides the high levels accounted for acetylsalicylic acid and acetaminophen in surface waters from the region of São Paulo, it is important to mention that these pharmaceuticals were not frequently detected in the samples, probably due to changes in the quantification limits of the methods over the years. Hu et al. 60 investigated the occurrence of 25 pharmaceuticals in a Chinese lake and observed average concentrations varying between 45 and 58 ng L -1 , depending on the season of the sampling. Tewari et al. 61 reported acetaminophen concentrations from 44 to 435 ng L -1 in six discharge canals of the Chao Phraya River, in Thailand as well as concentrations up to 70 ng L -1 in the receiving waters. The authors also showed that acetylsalicylic acid presented the highest concentration in the surface water samples ranging between 23 and 1100 ng L -1 . 61 An average concentration of 24 ng L -1 was depicted for triclosan in the positive samples, with minimum and maximum levels of 2 and 289 ng L -1 . Bedoux et al. 62 reviewed the occurrence of triclosan and its by-products in environmental samples revealing that concentrations worldwide are lower than those reported in this work. However, similar maximum concentrations were found in surface waters from Spain, 63 Costa Rica 64 and Brazil. 65 Higher maximum levels of triclosan were observed by dos Santos et al., 66 which investigated different rivers form the upper Iguaçu River Basin, in Brazil. Concentrations up to 415 ng L -1 were attributed to aspects such as direct raw sewage discharges into receiving waters, illegal sewer connections into pluvial drainage systems and inefficient removal at sewage treatment plants. 66

Illicit drugs
Considering the illicit drug class, average concentrations of cocaine and benzoylecgonine were 10 and 133 ng L -1 , respectively. The latter is commonly found at higher levels in environmental matrices due to its higher stability in aqueous samples. 67 Also, after the drug use, a higher proportion of benzoylecgonine is excreted in comparison with the unaltered drug. 68 Maximum concentrations of 62 and 1019 ng L -1 were also portrayed in Table 1 for cocaine and its metabolite, respectively. Concentrations in the investigated area were similar to those reported in countries such as Italy, 69 UK 70 and Spain. 71 Concentrations ranging from 1.1 to 115 ng L -1 for cocaine and from 2.3 to 520 ng L -1 for benzoylecgonine were reported in samples collected across Belgium. 72 In Brazil, Thomas et al. 73 investigated selected licit and illicit drugs in the urban streams of Manaus, in the Brazilian Amazon, and reported concentrations up to 5896 and 3582 ng L -1 for cocaine and benzoylecgonine, respectively. The authors verified that concentrations of both contaminants increased downstream of the urban region suggesting a continuous input along the entire length of the stream probably through raw sewage inputs.

Hormones
Hormones were detected in a wide range of concentrations as can be seen in Figure 2. Higher concentrations were noticed for the endogenous estrogen 17β-estradiol (average of 969 ng L -1 ), followed by the synthetic contraceptive 17α-ethynylestradiol (average of 777 ng L -1 ). Despite their relatively high concentrations in comparison with the other investigated hormones, they were not frequently detected in the samples probably due to limitations of the quantification methods. On the other hand, the concentrations exceeded AA EQS values proposed by EU for these compounds (0.04 ng L -1 for 17β-estradiol and 0.035 ng L -1 for 17α-ethynylestradiol). 6 According to Tiedeken et al. 6 the concentrations of 17β-estradiol and 17α-ethynylestradiol in European surface water samples are typically below 50 and 10 ng L -1 , respectively. Hormones such as estriol and estrone were more frequently detected, but at lower average concentrations (5 to 38 ng L -1 ). Levonorgestrel, other synthetic estrogen used as a contraceptive, was found at relatively high concentrations, i.e., from 1 to 663 ng L -1 . Quadra et al. 74 evidenced that hormone concentrations in Brazilian waters are significantly higher in comparison to data collected worldwide. In addition, the authors concluded that hormones were more frequently studied in the environment probably due to their potential impact on humans. However, it is important to point out that hormones are also potentially harmful to the aquatic biota, which justifies the large number of reports involving this class of emerging contaminants. Torres et al. 75 investigated the presence of hormones in the Piracicaba River, Brazil, and observed highest levels for 17α-ethynylestradiol and 17β-estradiol, in accordance with our results. Maximum hormone concentrations in Piracicaba River were assigned for 17α-ethinylestradiol (194 ng L -1 ), followed by 17β-estradiol (137 ng L -1 ), estriol (90 ng L -1 ), estrone (28 ng L -1 ), and progesterone (26 ng L -1 ). In water sources from the metropolitan area of Belo Horizonte, Brazil, natural and synthetic estrogens were detected in 15% of the samples and always at low concentrations, from 2 to 54 ng L -1 . 76 In the Jaboticabal Region, in Brazil, Lopes et al. 77 reported 17β-estradiol concentrations ranging from 16 to 30.6 ng L -1 in water samples from the Rico Creek only at the dry season.

Industrial compounds
Although six phthalates were considered in this work, only two were frequently investigated in surface water samples, i.e., diethylphthalate (DEP) and dibutyl phthalate (DBP), but only the latter was frequently detected. Dibutylphthalate levels ranged between 1.3 and 33.1 µg L -1 , being one of the most abundant substances together with the non-prescribed pharmaceuticals caffeine, acetylsalicylic acid and acetaminophen. The occurrence of three phthalates was investigated in German surface water samples by Fromme et al. 78 that shows DBP was found in minor concentrations in comparison with butylbenzylphthalate (BBP) and di(2-ethylhexyl)phthalate (DEHP), ranging from 0.12 to 8.80 µg L -1 . DEHP concentration in German waters varied between 0.33 to 97.8 µg L -1 and it was mainly influenced by water runoff from dumps. In the present work DEHP was investigated only in two samples being always present at an average concentration of 0.465 µg L -1 . Higher concentrations of DBP, in comparison with other phthalates, were also observed elsewhere. In rivers samples from Manchester, UK, Fatoki and Vernon 79 reported concentrations of DBP varying from 12.1 to 33.5 µg L -1 , while Vitali et al. 80 observed DBP concentrations up to 31.2 µg L -1 in surface water samples from the Rieti Basin in Italy.
Bisphenol A (BPA) levels in the present study varied between 2 to 13016 ng L -1 in the surface water samples, with an average concentration of 513 ng L -1 . A relatively high concentration of BPA was also reported by Jin et al. 81 in the Haihe River estuary, in China. The authors attributed an unusual concentration of 8.3 µg L -1 to a sewage outlet nearby the sampling point. 81 Kolpin et al. 58 observed an average BPA concentration of 740 ng L -1 for low flow U.S. streams, while at high flow conditions this contaminant was not detected. Santhi et al. 82 detected BPA in 93% of the investigated surface waters from the Langat River basin, in Malaysia. Concentrations up to 215 ng L -1 were reported in samples used as potable water, while six fold higher levels were detected in samples collected nearby industrial and sewage treatment plant outlets. 82 Huang et al. 83 reviewed the levels of BPA in water samples from China in comparison with other countries and concluded that concentrations are much higher in cities located at highly developed industrial and commercial regions. In Brazil, Moreira et al. 84 reported the presence of BPA in all surface water samples from das Velhas River with concentrations varying from 8.6 to 168 ng L -1 . In water samples from the upper Iguaçu River, in Brazil, Froehner et al. 85 reported BPA levels ranging between 620 to 12610 ng L -1 , i.e., very close to the results collected in this work. The authors concluded that the highest BPA levels were influenced by industrial and densely populated areas nearby the sampling points, which present a history of pollution by both domestic and industrial sewage inputs.

Drinking water
Two hundred and eighty-nine (289) samples of drinking water were analyzed by the group from 2007 to 2015. The number of samples analyzed, the frequency of positive samples and the concentration of each compound are presented in Table 2 and in the SI section. Figure 3 shows a plot representation of concentrations and the preliminary risk assessment.  (231), testosterone (215), 4-n-nonylphenol (189), triclosan (186) and atrazine (179).
The most prevalent substance among the most frequently investigated was caffeine (90%), atrazine (73%), BPA (19%) and triclosan (19%). All four substances, except for atrazine, appeared in the list of selected chemicals proposed by Glassmeyer et al. 89 that are useful as tracers of human wastewater. In fact, the presence of caffeine in Brazilian waters are commonly associated with sewage contamination, as previously reported elsewhere. 90 All hormones were investigated in more than 100 samples, being the most investigated class of emerging contaminants. Although their frequency of detection varied from zero to 4%.

Pesticides
Few data are available in the literature about the presence of pesticides in drinking waters. Average and maximum atrazine concentrations in this work were of 36 and 687 ng L -1 . In Brazil, lower concentrations, between 2.0 and 24 ng L -1 , were previously reported in drinking water samples from 16 capitals. 19 Sodré et al. 91 reported concentrations ranging from 2 to 3.3 ng L -1 in all samples collected from the Descoberto production system in the Brazilian Federal District. Median and maximum concentrations around 50 and 930 ng L -1 were observed in finished and distribution U.S. drinking waters by Benotti et al., 92 while Padhye et al. 93 noticed a median concentration of 15 ng L -1 in finished and drinking waters from a major U.S. urban region. Imidacloprid was detected in concentrations between 3 and 63 ng L -1 in this work and similar concentrations were detected in tap waters from Iowa City in the United States. 94 Otherwise, the average concentration of fipronil detected in this work was 5 ng L -1 , which is less than the average concentration (160 ng L -1 ) found in Vietnam. 95

Pharmaceuticals and personal care products
Recently, Machado et al. 19 investigated the presence of caffeine in drinking waters collected in different Brazilian capitals and showed higher average levels in samples collected in the city of Porto Alegre (1211 ng L -1 ), probably due to cultural peculiarities related to the consumption of yerba mate (Ilex paraguariensis) on a daily basis. The city of São Paulo was ranked in the top five capitals with an average concentration of 121 ng L -1 .
Triclosan was found at concentrations ranging from 2 to 37 ng L -1 in this work. Yavuz et al. 96 studied the levels of biocides in source and drinking water systems of Ankara, Turkey and noticed that the concentration of triclosan in all samples were below the limit of quantification of 0.87 ng L -1 . Benotti et al. 92 investigated finished and distribution waters from 19 U.S. water utilities and found triclosan in only one finished water sample at a concentration of 1.2 ng L -1 . Concentrations below 1 ng L -1 were reported by Servos et al. 97 investigating finished water produce from different types of source waters in the region of Ontario, Canada. In the nationwide survey carried out in Brazil, triclosan was detected in only one sample collected in the city of Porto Alegre at a concentration of 10 ng L -1 . 19 A triclosan concentration as high as those depicted in this work was found by Padhye et al. 93 investigating the removal of emerging contaminants in a drinking water supplier in the Southeast United States. The authors reported a frequency of detection of 63% for triclosan with median and maximum concentrations of 1.4 and 59.6 ng L -1 .

Illicit drugs
Illicit drugs were investigated only in few samples. Average concentration for benzoylecognine was higher (212 ng L -1 ) in comparison to cocaine (14 ng L -1 ). Huerta-Fontela et al. 98 showed that benzoylecgonine is persistent during conventional drinking water treatment, while cocaine is efficiently removed after post-chlorination and filtration in granulated activated carbon. Boleda et al. 99 investigated the presence of drugs of abuse in tap waters collected in Europe, Japan and Latin America and reported lower values for the two cocaine compounds in comparison to our data. According to the authors, an average concentration of 0.4 ng L -1 was accounted for both compounds in samples collected in Spanish tap waters, while concentrations up to 0.6 and 4.5 ng L -1 were observed for cocaine and benzoylecgonine, respectively, in Latin American countries. 99

Hormones
Concentrations of hormones up to 125,35,32,24 and 20 ng L -1 were observed for estriol, 17β-estradiol, 17α-ethinylestradiol, progesterone and estrone, respectively. These values are considered high in comparison with data reported elsewhere. Kuster et al. 100 verified that all drinking water samples investigated from the Llobregat River Basin (Barcelona, Spain) were free of estrogens with the exception of one sample where estriol was found at 11.6 ng L -1 . In German, it was reported maximum concentrations of estrone, 17α-estradiol, 17β-estradiol and 17α-ethinylestradiol varying from 0.3 to 2.1 ng L -1 . 101 Benotti et al. 92 investigated five hormones in U.S. drinking waters being all results presented as below the maximum residue limit, except for progesterone that was found in one finished water samples at a concentration of 0.57 ng L -1 . Similar result was reported by Cai et al. 102 investigating 17β-estradiol in drinking waters from Beijing, China. In Brazil, hormones were investigated in finished waters from different cities with wide range concentrations (0.3-2170 ng L -1 ). 92

Industrial compounds
The presence of BPA in Brazilian drinking waters has been also associated with the contamination of source waters by sewage discharges. Despite its industrial origins, BPA is part of our life style and it is still being used in a variety of manufactured consumer goods in Brazil. In this work, concentrations up to 178 ng L -1 were noticed in the investigated samples, with an average level of 23 ng L -1 . Kuch and Ballschmiter 101 found BPA in drinking water samples from Southern Germany in concentrations ranging from 0.3 to 2 ng L -1 . Concentrations ranging from 3.5 to 59.8 ng L -1 were reported by Santhi et al. 82 in tap water samples from the region of Kuala Lumpur, in Malaysia. The authors also suggested that higher concentrations were probably caused by polyvinyl chloride (PVC) fittings and filter devices connected to the taps. 82 Bisphenol A concentrations up to 128 ng L -1 , thus similar to our results, were also found in finished water samples from 62 drinking water treatment plants in 31 major cities across China. 103 Considering the alkylphenols 4-n-nonylphenol and 4-n-octylphenol, average concentrations of 114 and 225 ng L -1 , respectively, were observed, although a considerably high concentration of 4-n-octylphenol (1137 ng L -1 ) was also noticed. Kuch and Ballschmiter 101 also found alkylphenols in drinking waters from Germany, but the maximum concentrations were 4.9 and 16 ng L -1 for 4-tert-octylphenol and 4-n-nonylphenol, respectively. For the latter substance, Shao et al. 104 reported concentrations varying from 10 to 2700 ng L -1 in drinking water samples collected in the area of Chongqing (China).
Dioctylphthalate presented the highest levels among the investigated contaminants with an average concentration of around 7 µg L -1 . Casajuana and Lacorte 105 evaluated the release of phthalic esters in drinking water reporting dibutylphthalate concentrations varying between 0.016 and 0.032 µg L -1 in distribution waters collected from public fountains in the Vallbs area (Catalonia, Spain). In mineral waters bottled in polyethylene, the authors verified slightly higher concentrations when compared to the distribution of water samples. 105 Serôdio and Nogueira 106 showed that dibutylphthalate was the most abundant phthalate whether in tap (0.52 µg L -1 ) or in bottled mineral (0.35 µg L -1 ) waters from Lisbon, Portugal. An average concentration of 2.73 µg L -1 , thus similar to our results, were reported by Loraine and Pettigrove 107 investigating emerging contaminants in finished drinking water in San Diego (USA).

Groundwaters
The results collected for the analyses of the thirty-three groundwater samples are presented in Table 3. The samples were collected in Guarani Aquifer (n = 12), Bauru Aquifer (n = 14) and in wells used to public supply in the rural area of the city of Campinas (SP) (n = 7).
The most prevalent contaminant was the herbicide atrazine detected in 86% of samples, followed by caffeine (55%) and bisphenol A (50%). Triclosan, estrone, testosterone and progesterone were detected in less than 10% of the samples. This is of great concern since most of the samples were collected in deep wells (> 300 meter), in which anthropic contamination is expected to be more controlled.
Pesticides are the most investigated contaminants in groundwater samples around the world due to the potential of contamination by percolation through agricultural soils located nearby aquifer recharge areas. Dujaković et al. 108 investigated 14 pesticides, including atrazine, in 6 groundwater samples from Belgrade, Serbia, and reported results below the limit of detection (1.5 ng L -1 ). Schipper et al. 109 detected pesticides in 27% of the 771 groundwater samples from The Netherlands.
Loos et al. 110 investigated the occurrence of organic pollutants in 164 ground water samples collected from 23 European countries. Caffeine was quantified in 83% of the samples with a maximum concentration of 189 ng L -1 , thus being higher than the maximum concentrations reported in this work (53 ng L -1 ). In addition, although atrazine presented a high frequency detection, its concentrations varied between 2 and 5 ng L -1 , which were about 250 times lower that the levels found by Loos et al. 110 Bisphenol A presented the highest concentration in this work (643 ng L -1 ), but also lower than maximum levels found in Europe (2300 ng L -1 ). 110 Rabiet et al. 111 detected caffeine concentrations between 1.5 and 23 ng L -1 as well as other pharmaceuticals in wells used to supply drinking water to the Hérault watershed located in the south of France. Caffeine was also detected with other pharmaceuticals in groundwater samples from a drinking-water supply system located in the state of California (USA) in concentrations varying between 170 and 290 ng L -1 . 112 Despite the relatively high concentrations of caffeine in the samples, the authors found this contaminant only in 1% of the 1231 samples investigated.
It is known that the evaluation of the risks related to the presence of emerging contaminants in groundwater requires a more detailed study on the presence of metabolites 113 and the effects caused by the mixture of several compounds, which confers greater complexity in the choice of contaminants that must be monitored, but provides a more realistic scenario of the contamination stage.

Preliminary risk assessment
The preliminary risk assessment is shown in Figures 2  and 3 by colors in the plot representations. Green represents the concentrations below WQC (risk quocient smaller than 1), red represents the concentrations higher than the WQC (risk quocient greater than 1) and in grey are the compounds for which no WQC was found in the literature.
As no WQC were found, we used Binelli et al. 114 findings as reference for illicit drugs and their main active metabolites. They showed that zebra mussels exposed to 0.5 and 1.0 µg L -1 of benzoylecgonine induced changes in some proteins of gill cells that are crucial for the overall metabolism. The concentrations found in our study are greater than those values, therefore, representing a possible risk to the biota. 114 From the 39 substances analyzed in drinking water samples, drinking water criteria were available for 22 (Table 2). For pesticides, the quocient ranged from 10 -2 to 10 -7 , therefore, no adverse effects are expected at the concentrations found, except for atrazine, which presented a risk quotient of 0.3 indicating some level of concern. Only for 17β-estradiol the risk quotient was higher than 1, although the detection frequency was only 2% (Table 2, Figure 3).
Considering the number of groundwater samples analyzed in this study was limited to 33, the data are not representative for risk assessment, but the concentrations of testosterone, bisphenol A and triclosan detected in groundwater (Table 3) are in the range of risk for aquatic life (Table 1) and testosterone was detected in concentration also above the drinking water criteria ( Table 2).

Raw and treated sewage
The concentrations of each target compound in raw and treated effluents collected during 2010 and 2013 of the WWTPs located in the cities of Campinas and São José do Rio Preto in the state of São Paulo are presented in Figure 4 and in the SI section. Some emerging contaminants were still detected in relatively high concentrations after treatment process, but the average levels concentrations were below the levels suggested for secondary effluents by the Australian Guidelines for Water Recycling for caffeine (350 ng L -1 ), estrone (30 ng L -1 ), estriol (50 ng L -1 ), triclosan (350 ng L -1 ), bisphenol A (200 µg L -1 ) and atrazine (40 µg L -1 ). 115 The wastewater treatment technology applied in São Paulo State Brazil (based on 57 raw and treated wastewater samples) was not effective for the removal of most of the investigated compounds (Figure 4). Conventional secondary treatments, such as aerated lagoon, activated sludge, upflow anaerobic sludge blanket (UASB), moving bed biofilm reactor, physical-chemical treatments, and others, are usually non-effective to remove these contaminants. 115 Better removal rates are commonly achieved using tertiary treatments, such as photodegradation (UV/H 2 O 2 ), Fenton (Fe/H 2 O 2 ) or photo-Fenton (UV/Fe/H 2 O 2 ), ozone (O 3 ), electrochemical oxidation, granular activated carbon, or the use of membrane system (microfiltration, ultrafiltration, nanofiltration, reverse osmosis or bioreactor), [115][116][117] but these treatments are not applied in WWTPs from São Paulo.
Caffeine and benzoylecgonine were detected in 100% of raw and treated effluents. Caffeine remained in the effluent at relatively high concentrations (2707 ng L -1 ) probably due to its high solubility in water and molecular stability. 117,118 Furthermore, high levels of dissolved organic matter in water or wastewater decrease the efficiency of organic contaminants removal by conventional treatments.
The average of caffeine concentration found in treated wastewater (333 ng L -1 ) is higher than the levels reported by Blackbeard et al. 115 in a reuse water treatment plant in Australia (100 ng L -1 ), but similar to those reported by Arzate et al. 119 for a secondary wastewater treatment effluent in Spain (409 ± 52 ng L -1 ). To enhance the removal of caffeine and other emerging contaminants from secondary effluents, the latter authors employed a photo-Fenton approach to reach up to 93% of caffeine removal.
Other contaminants of concern such as estrone, BPA, and atrazine were not efficiently removed by the conventional wastewater treatment (Figure 4) operating in Campinas and São José do Rio Preto cities. Estrone was more frequently detected, so it could be more stable than the other endogenous hormones. Triclosan levels in the treated wastewater detected in this work were similar to those reported for secondary effluent from a reuse water treatment in Australia (20 ng L -1 ). 115 As observed for surface and drinking waters, the levels of cocaine in raw and treated wastewaters are lower in comparison to its major metabolite, benzoylecgonine. The low removal efficiency for benzoylecgonine during the wastewater treatment was also observed in India, where no significant differences between the concentrations in both influent and effluent were depicted. 120 However, concentrations of benzoylecgonine in the raw wastewater investigated in this work were almost 50 times higher than those reported in India. In Spain, Bijlsma et al. 121 reported levels of cocaine and benzoylecgonine in raw wastewater similar to those found in the present work. Despite the similar concentrations, the Spanish WWTP was able to remove almost 75% of these contaminants whereas in Brazil, a lower removal efficiency was observed to the treatment processes. It is possible to observe that the median concentration of caffeine in inland surface waters (approximately 1000 ng L -1 ) is almost ten times higher than the median concentration reported for treated wastewaters. Thus, this behavior is clear evidence of the role of untreated sewage discharges, since the caffeine content in the raw waste may be completely transferred to surface waters without any removal step in WWTP. This result shows that WWTP effluents are not the main sources of surface water contamination in the state of São Paulo. On the contrary, this role should be attributed to the direct disposal of raw wastewater into receiving waters. This same conclusion was made by Sodré et al., 122 who investigated the presence of caffeine and bisphenol A in drinking waters from the city of Campinas. The authors showed that high levels of caffeine in drinking waters were related to the low rates of sewage treatment in the region. In the present work, the data portrayed in Figure 5 not only corroborate the result raised by Sodré et al., 122 but also evidenced that the same route of contamination is prevalent in other cities of the state.

Removal of anthropogenic activity indicators by wastewater and drinking water plants
In the Brazilian scenario, the term reuse water holds an underlying meaning, which is the classical indirect and unplanned reuse water. Nowadays, the watersheds are used to solve the supply demands in one region at the expense of another. A water body is used as source water for one city, which also inputs its sewage downstream (mostly untreated sewage). Further, downstream another city uses the same river to supply their population with drinking water. As a result, some cities face the challenge of getting a drinking water treatment plant (DWTP) intended for treating surface water that actually has to treat water with the quality of sewage water.
Caffeine has been used as a marker of human activities since the work of Siegener and Chen, 123 which investigated how pharmaceuticals were entering Boston Harbor. The authors proposed the use of caffeine as a tracer of anthropogenic inputs into marine systems as highly caffeinated freshwaters were mixed with uncontaminated sea waters. Since then, a number of reports have consolidated the use of caffeine as an indicator of wastewater impacted waters, [124][125][126][127] including in Brazil. 90,122,128,129 Figure 6 shows concentrations of cocaine and benzoylecgonine determined in drinking water samples collected in 5 cities of the state of São Paulo (Limeira, Campinas, Espírito Santo do Pinhal, Santa Barbara d'Oeste and Piracicaba), in surface water samples from 16 different rivers, and in treated and raw wastewater samples collected in 2 WWTPs of the city of Campinas.
The concentrations of benzoylecgonine were higher than cocaine in all analyzed samples. This could be related to the excretion metabolism of this illicit drug since 25% of the consumed crack/cocaine is excreted in urine as benzoylecgonine. 130 In addition, considering that the main source of illicit drugs to the environment is the human excretion, a variety of studies show that benzoylecgonine is more stable than cocaine under different environmental conditions such as temperature and pH. 67,121,131 Different removal efficiencies for cocaine and benzoylecgonine were observed for each conventional WWTP evaluated. While a given WWTP could reduce cocaine in 55%, the metabolite concentration was only 8% reduced; the other WWTP was able to decrease more than 95% of the benzoylecgonine concentration and almost 100% of the parental drug. 15 These results may suggest that the concentration of illicit drugs, mainly benzoylecgonine, observed in surface and drinking water samples come from discharges of both treated and raw wastewaters into aquatic bodies. The differences observed during the treatment processes could also be observed in DWTPs. Benzoylecgonine concentrations in drinking water samples were similar to those observed in surface water samples, whereas cocaine concentration, when detected, was almost 10-fold less.
As for benzoylecgonine, many other emerging contaminants (e.g., psychoactive substances) pose a challenge to be removed by conventional treatment processes employed in WWTPs and DWTPs notably due to their high hydrophilicity. The use of water indicators based on different types of emerging contaminants have been considered to improve monitoring and to complement the current strategies related to watersheds management. 132,133 Considering the stability and the high frequency of detection of benzoylecgonine in environmental aquatic matrices, this substance may be a reliable biomarker for both surface and drinking waters, mainly as an indicator of the sanitation conditions, since this cocaine metabolite can persist in the aquatic environment even after two weeks from the discharge, 67 and it has been also detected in surface and drinking waters in countries with better sanitation index in comparison with Brazil. [134][135][136]

Conclusions
The set of data used in this work (58 target compounds and 708 samples collected from 2006 to 2015) represented the snapshot of the contamination in the state of São Paulo. Some compounds such as caffeine, estrone, 17β-estradiol, 17α-ethynylestradiol, bisphenol A, atrazine, carbendazim, fipronil, malathion and imidacloprid were highlighted as priority concern among the studied contaminants considering the frequency of detection and the number of positive samples above WQC.
The preliminary risk assessment for aquatic life protection were identified for caffeine, paracetamol, diclofenac, 17α-ethynylestradiol, 17β-estradiol, estriol, estrone, testosterone, triclosan, 4-n-nonylphenol, bisphenol A, atrazine, azoxystrobin, carbendazim, fipronil, imidacloprid, malathion and tebuconazole, but considering the drinking water criteria available in the literature no adverse effects to human health were expected to the studied scenario, except for 17β-estradiol.
Due to their stability in water and high frequency of detection in all studied aquatic matrices, caffeine and benzoylecgonine showed as good anthropogenic indicators, especially to the regions with high-density population and poor sanitation system.

Supplementary Information
Concentrations of the target substances in each sample (n = 708) collected from 2006 to 2015 as well as information about sampling location, year of sampling, instrumental tools and method used to determine the contaminants are given as supplementary information available free of charge at http://jbcs.sbq.org.br.