Next Article in Journal
PROMANCOA Modular Technology for the Valorization of Mango (Mangifera indica L.) and Cocoa (Theobroma cacao L.) Agricultural Biowastes
Next Article in Special Issue
Mechanical Characteristics and Micro-Mechanism of Modified Dredged Sludge Based on Calcium-Containing Solid Waste Used as Landfill Cover Materials
Previous Article in Journal
Relationship between Total Phenolic Content, Antioxidant Capacity, Fe and Cu Content from Tea Plant Samples at Different Brewing Times
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Mechanistic Model to Assess the Fate of Naphthalene and Benzo(a)pyrene in a Chilean WWTP

1
Departamento de Ingeniería Civil, Universidad de Chile, Santiago 8370449, Chile
2
Centro Avanzado Para Tecnologías del Agua (CAPTA), Universidad de Chile, Santiago 8370449, Chile
*
Author to whom correspondence should be addressed.
Processes 2021, 9(8), 1313; https://doi.org/10.3390/pr9081313
Submission received: 1 June 2021 / Revised: 20 July 2021 / Accepted: 22 July 2021 / Published: 29 July 2021
(This article belongs to the Special Issue Activated Sludge Treatment Processes)

Abstract

:
Polycyclic aromatic hydrocarbons (PAHs) are a family of organic compounds of widespread presence in the environment. They are recalcitrant, ubiquitous, prone to bioaccumulation, and potentially carcinogenic. Effluent from wastewater treatment plants (WWTPs) constitutes a major source of PAHs into water bodies, and their presence should be closely monitored, especially considering the increasing applications of potable and non-potable reuse of treated wastewater worldwide. Modeling the fate and distribution of PAHs in WWTPs is a valuable tool to overcome the complexity and cost of monitoring and quantifying PAHs. A mechanistic model was built to evaluate the fate of PAHs in both water and sludge lines of a Chilean WWTP. Naphthalene and benzo(a)pyrene were used as models of low-MW and high-MW PAHs. As there were no reported experimental data available for the case study, the influent load was determined through a statistical approach based on reported values worldwide. For both naphthalene and benzo(a)pyrene, the predominant mechanism in the water line was sorption to sludge, while that in the sludge line was desorption. Compared to other studies in the literature, the model satisfactorily describes the mechanisms involved in the fate and distribution of PAHs in a conventional activated sludge WWTP. Even though there is evidence of the presence of PAHs in urban centers in Chile, local regulatory standards do not consider PAHs in the disposal of WWTP effluents. Monitoring of PAHs in both treated effluents and biosolids is imperative, especially when considering de facto reuse and soil amendment in agricultural activities are currently practiced downstream of the studied WWTP.

1. Introduction

Polycyclic aromatic hydrocarbons (PAHs) are widespread environmental pollutants formed during the incomplete combustion of organic materials [1]. They are hydrophobic, with low solubilities in water and low vapor pressure [2]. They can be classified as low molecular weight (LMW—two or three fused benzene rings) or as high molecular weight (HMW—four or more rings). As their molecular weight increases, their hydrophobicity increases, and their water solubility and vapor pressure decrease, which makes these compounds more recalcitrant [3]. Due to their persistence in the environment, and the carcinogenic and mutagenic effects of some PAHs, they are listed as US-EPA and EU priority pollutants to be monitored in the urban water cycle [4,5].
Even though their sources into the environment can be both natural and anthropogenic, the major source of PAHs into water bodies in urban areas are wastewater treatment plants (WWTPs) [6,7]. Several authors have studied and measured PAHs in WWTP [8,9,10,11,12], reporting a wide range of both influent and effluent concentrations (Table 1). In these studies, LMW PAHs were in most cases present in higher concentrations in both influent and effluent in the dissolved phase only, whereas HMW PAHs tended to sorb to particulate matter present in the sewers, thus entering the WWTP sorbed to particles [9,13]. Evidence suggests that the removal of PAHs depends on the treatment train and/or technology of the WWTP, as well as its operational parameters, influent wastewater composition, and seasonal/climatic conditions [7,14]. Thus, the influent concentrations are highly variable and inherent to each specific case study.
In conventional activated sludge plants, the main removal mechanisms of PAHs are adsorption to sludge, volatilization, and biotransformation [14,15], while abiotic degradation mechanisms such as photolysis and hydrolysis are not effective removal mechanisms [16]. The physicochemical properties of PAHs can also impact their removal. For instance, highly hydrophobic HMW PAHs are expected to be removed by sorption to sludge at higher rates than LMW PAHs [17], while it is expected that LMW PAHs are removed at higher rates by volatilization and biodegradation [14].
Understanding and quantifying PAH removal in WWTPs can provide a better awareness of their presence in the environment, especially when considering reuse applications for the treated wastewater and biosolids. However, testing and monitoring PAHs in WWTPs is not always feasible due to the costs and complexity associated with the experimental methods for their detection in the water and sludge matrices [18]. Whenever experimental data are unavailable, using modeling tools to assess the fate of PAHs in WWTP can be quite valuable [18,19]. To date, several authors have developed models to predict the fate of micropollutants in wastewater treatment. Although models can differ in their complexities (e.g., different number of state variables and treatment units), most of those applicable to PAHs consider two phases for the fate of micropollutants in the water matrix: as dissolved and as sorbed to particles. Models for the fate of PAHs in the sludge matrix consider either two or three phases: freely dissolved, sorbed to dissolved and colloidal matter (DCM), and sorbed to sludge particles. Table 2 summarizes some studies focusing on modeling the fate of different micropollutants in wastewater treatment. The available models for the water line are extensive, considering various treatment units such as conventional activated sludge, membrane bioreactors, and primary clarifiers. For the sludge treatment, however, available models are limited.
In this study, we constructed a mechanistic model to assess the presence and removal of PAHs in a Chilean WWTP. Naphthalene and benzo(a)pyrene were used as model contaminants to evaluate the fate of LMW and HMW PAHs. Results were compared to those obtained by different experiences available in literature, and the potential implications of the presence of PAHs for the reuse of the treated wastewater and biosolids are also discussed.
Table 1. Summary of influent and effluent concentration of PAHs in WWTPs worldwide. For WWTPs with industrial inputs, the percentage of domestic/municipal WW is provided.
Table 1. Summary of influent and effluent concentration of PAHs in WWTPs worldwide. For WWTPs with industrial inputs, the percentage of domestic/municipal WW is provided.
LocationInflow and Population ServedTreatment TechnologyNaphthalene [ng/L]Benzo(a)pyrene [ng/L]∑LMW PAHs [ng/L]∑HMW PAHs [ng/L]∑PAHs [ng/L]Reference
InfluentEffluentInfluentEffluentInfluentEffluentInfluentEffluentInfluentEffluent
Belo Horizonte, Brazil290,000 m3/d
1,600,000 PE
CAS114.346.125.827.6290.7183.950.354.4553.5483.7[20]
155,000 m3/d
1,100,000 PE
UASB with trickling filters217.5102.826.325.4483.5270.350.850699.6472.4
1600 m3/d
30,000 PE
WSP435.861.425.127.4996219.869.454.91210.8431.1
Konya, Turkey200,000 m3/d
1,300,000 population
Unknown content of industrial input
4 stage Bardenpho2282560724623865942291832615777[21]
16,500 m3/d
110,000 population
WSP1429370394315094361931811702617
210 m3/d
2000 people
CW3601034236444161197173641334
Guangzhou, China220,000 m3/d
Municipal ww
Inverted A2/O435.580.1NDND551.1143.899.079.9650.1223.8[17]
Harbin, China650,000 m3/d
80% domestic
A/O and BAF882471100ND1866.3859.61039.9114.12906.2973.7[12]
Hefei, China300,000 m3/d
60% municipal ww
A/O1437.3924.41383.5380.420731327.63685.8912.85758.82240.4[22]
Tai’an, China50,000 m3/d
60% municipal ww. In autumn
A2/O971.1207.9NDND1129.425827.53.11156.9261.1[23]
50,000 m3/day
60% municipal ww. In spring
205.656.7NDND659.5123.7488.95.51148.4129.2
Italy25,000 m3/d
30% municipal
CAS2505616ND73513321590950223[9]
15,000 m3/d
100% municipal
CAS9673NDND3071275612363139
19,500 m3/d
100% municipal
CAS11343NDND195438720350
118,000 m3/d
60% municipal
CAS6347414ND1315131356371671168
21,000 m3/d
90% municipal
CAS1033720ND32483189ND51383
Norway106,750 m3/d
Domestic ww
A/O3152613ND55830242180031[24]
35,450 m3/d
Domestic ww
Chemical treatment3254457ND72181519840925855
21,050 m3/d
Domestic ww
Chemical treatment8410516ND2451452255470150
23,125 m3/d
Domestic ww
Chemical treatment2501708ND38821515717545232
10,150 m3/d
Domestic ww
Mechanical treatment122105NDND1751457255250190
Venice, Italy100,000 m3/d
90% municipal
Denitrification, oxidation + nitrification18517290351214121176[8]
Thessaloniki, Greece40,000 m3/d
Municipal ww
CAS73005000224.810,49055471023.7137.811,513.75684.8[10]
Montreal, Canada1,300,000 m3/d
85% municipal
1,400,000 population
Coagulation-flocculation14788481582231603601185291[11]
CAS—Conventional Activated Sludge, UASB—Up-flow sludge blanket reactor, WSP—Waste stabilization pond, CW—Constructed wetland, A/O—anaerobic-oxic reactor, BAF—Biological Aerated Filter, A2/O—Anaerobic-anoxic-oxic, ∑LMW PAHs: naphthalene, acenaphthene, acenaphthylene, fluorene, phenanthrene, anthracene; ∑HMW PAHs: fluoranthene, pyrene, benzo(a)anthracene, chrysene, benzo(b)fluoranthene, benzo(k)fluoranthene, benzo(a)pyrene, benzo(g,h,i,)perylene, dibenzo(a,h)anthracene, indenol(1,2,3-c,d)pyrene; ∑PAHs: ∑LMW PAHs and ∑HMW PAHs. Ref. [20] does not consider benzo[k]fluoranthene. Ref. [21] only considers ∑PAHs: naphthalene, anthracene, fluoranthene, benzo(b)fluoranthene, benzo(k)fluoranthene, benzo(a)pyrene, benzo(g,h,i,)perylene, indenol(1,2,3-c,d)pyrene.
Table 2. Summary of studies focusing on modeling the fate of different micropollutants in wastewater treatment.
Table 2. Summary of studies focusing on modeling the fate of different micropollutants in wastewater treatment.
DescriptionModeled Treatment StagesModeled MicropollutantsProcessesReference
Studies using the 2-phase approach to model the fate of micropollutants
Fate of organic micropollutants in novel WWTPsCEPT, HRAS, RBF, PN-AMXPharmaceuticalsSorption, biodegradation (pseudo-first-order reaction)[25]
Modeling of micropollutant removal in full-scale MBRsMBR51 MPs: 24 metals, 10 beta blockers, 8 pharmaceutical products, 2 alkylphenols, 3 PAHs, and 4 pesticidesBiodegradation (dissolved and particulate phase) in oxic and anoxic conditions and sorption/desorption[26]
Mechanistic modeling approach to anaerobic systemsAD20 OMPsSorption, desorption, biotransformation in the liquid and solid phase[27]
Fate of micropollutants in aeration tank and secondary clarifierAS + SCMPs (androstenedione, ibuprofen, triclosan, and caffeine)Biological degradation (first-order reaction) and sorption[28]
Activated Sludge Model for Xenobiotics (ASM-X)ASPharmaceuticalsCo-metabolic biotransformation, retransformation in oxic and anoxic conditions and sorption/desorption[29,30]
Dynamic transport and fate of micropollutants in integrated urban wastewater and stormwater systemsSewers; WWTPs (PR, AS, SC, F, ST, AD, SD; river (water, sediment)MPsSettling, resuspension, volatilization, filtration/separation, sorption, desorption, hydrolysis, photolysis, biological biodegradation (aerobic, anoxic/anaerobic)[31]
Model for co-metabolic biotransformation of organic micropollutants in NAS reactorsNASPharmaceuticalsCo-metabolic biotransformation, sorption/desorption and volatilization[32]
Mechanistic model for fate of estrogens in inverted A/O reactorsAnoxic/aerobic reactorEstrogensCleavage, adsorption, desorption, aerobic biodegradation, anoxic biodegradation[33]
Micropollutant fate in SBRs during startup and steady stateSBRCadmium and leadAdsorption[34]
Steady-state model for the fate of hydrophobic and volatile compounds in activated sludgePR, ASVOCsBiodegradation, bubble and surface volatilization, sorption to particles[35]
Studies using the 3-phase approach to model the fate of micropollutants
Modelling PAHs removal in activated sludge process:
effect of disintegration
ASPAHsBiodegradation, volatilization, sorption/desorption to particulate matter and sorption/desorption to DCM.[15]
Dynamic model for bioavailability and co-metabolism of micropollutants during anaerobic digestionADPAHsCo-metabolic biodegradation, volatilization, sorption/desorption to particles and DCM[36]
A three-compartment model for micropollutants sorption in sludgeADPAHsBiodegradation, sorption to particles and DCM[37]
PR—primary, AS—activated sludge, SBR—sequencing batch reactor, MBR—membrane bioreactor, SC—secondary clarifier, F—filtration, CEPT—chemically enhanced primary treatment, HRAS—high-rate activated sludge, RBF—rotating belt filter, PN-AMX—partial nitritation unit, NAS—nitrifying activated sludge, ST—sludge thickener, SD—sludge dewatering, AD—anaerobic digestion. DCM—dissolved and colloidal matter.

2. Materials and Methods

2.1. Model Treatment Plant

The model WWTP is a conventional activated sludge plant located in Santiago, Chile. It has a influent flow of 8.8 m3/s, treating approximately 50% of the city’s domestic wastewater (equivalent to a population of 3.7 million people) [38]. The treated effluent is discharged to the Mapocho river, downstream from the WWTP, where it is further used for agricultural irrigation (de facto reuse). The plant’s biosolids (digested and stabilized sludge) are disposed either in agricultural lands or in designated landfills [39]. A simplified schematic of the treatment process and the plant’s data are presented in Figure 1.

2.2. PAH Influent Loads

Due to the lack of available experimental data, and to represent the variability in the influent wastewater characteristics, a bootstrap or resampling approach (n = 50,000) was applied to influent loads of PAHs from WWTPs around the world (Table 1). Naphthalene and benzo(a)pyrene were used as models of LMW and HMW PAHs. Naphthalene is usually found in higher concentrations in WWTPs and, due to its high vapor pressure, it can be used as a control compound to verify if volatilization occurs in the treatment process [9,12,22,23]. On the other hand, benzo(a)pyrene is considered among the IARC’s Group 1 carcinogens; its presence in WWTP has been reported in many studies and it can be used as a control compound to verify if sorption occurs [9,12,14,22]. Population served (>100,000 PE) and content of industrial waste (<50% of the influent flow) were the main criteria for data selection. Resulting loads of PAHs in the 10%, 50%, and 95% probability percentiles of occurrence were grouped, and the minimum, median, and maximum values of each group were used as inputs for the modeled WWTP. For the 95% percentile, these values correspond to 715.3, 760.1, and 795.4 g/day for naphthalene, and 82.9, 84.4, and 86.2 g/day for benzo(a)pyrene, respectively. Input loads resulting from the bootstrap are summarized in Figure 2. Additional data from the bootstrap calculations are provided in Tables S1 and S2.

2.3. Mechanistic Model for the Fate of PAHs in the WWTP

The constructed mechanistic model was based on a steady-state mass balance for each of the treatment stages. The mass balance equations for each of the phases in a treatment stage were solved simultaneously using a commercially available solver. The overall model was used to determine the mobility of the studied PAHs through the WWTP. The general mass balance equation can be expressed as
V · C i t = F in F out ± r i ,
where V is the reactor volume (L), Ci is the PAH concentration in the treatment stage (ng/L), Fin is the inlet PAH load (ng/day), Fout is the outlet PAH load (ng/day), and ri is reaction rate (ng/day). For each treatment stage, ri was modified to include the processes and/or partitioning mechanisms considered in the stage. Table 3 presents a summary of the selected models and transport/partitioning processes used in each of the treatment stages. Additional considerations and assumptions are itemized below:
  • According to different authors, the removal of micropollutants during preliminary treatment is close to null [40,41]. Moreover, experimental sampling by Manoli et al. (1999) reported negligible removal of PAHs by chlorination [10]. Thus, preliminary treatment and chlorination are not considered in the current model.
  • For the water line (Figure 1), PAHs were assumed to be distributed between the dissolved and sorbed-to-particle phase, and only the dissolved phase was bioavailable [42]. For the sludge line (Figure 1), three phases were considered: freely dissolved, sorbed to DCM, and sorbed to particles. In this case, only the PAHs freely dissolved and sorbed to the DCM were bioavailable [36].
  • Throughout all the treatment units, sorption was assumed to occur onto the total suspended solids fraction.
Table 3. Summary of the selected models and processes for each of the treatment stages.
Table 3. Summary of the selected models and processes for each of the treatment stages.
Treatment StageModeled ProcessesReference
Water linePrimary treatment:
Gravitational sedimentation
Advection, adsorption, surface volatilization,[35]
Secondary treatment: Conventional activated sludgeFor PAHs: advection, adsorption, biodegradation as co-metabolism (oxic and anoxic conditions), surface and bubble volatilization
For BOD: Activated Sludge Simple Model (ASSM)
[30,35,42,43]
Sludge linePrimary gravitational thickener
Secondary thickener (flotation)
Mixing chamber
Advection[37]
Anaerobic digestionAdvection, sorption/desorption to particles, sorption/desorption to DCM, volatilization, biodegradation as co-metabolism[36]
Dewatering with centrifugesAdvection[37]
The removal efficiency of each treatment stage was determined by
R % = ( F in F out ) F in · 100 ,
where R% is the removal efficiency, Fin is the inlet PAH load (ng/day), and Fout is the outlet PAH load (ng/day).
Details about the model equations are available in Tables S3–S7, as well as the summary of the parameters used in the model.

2.4. Sensitivity Analysis

The analysis was carried out as a function of parameters defining the governing partitioning mechanisms in each treatment stage. For the water line, the predominant mechanism was adsorption to sludge; thus, the solid–liquid partition coefficient (Kd) was used in the analysis for primary and secondary treatment. SRT was also considered in the latter. On the other hand, in the anaerobic digestion, the predominant mechanism is desorption from particles. Therefore, the first-order kinetic constant of sorption to particles (k1) and the equilibrium constant for PAH sorption to particle (KP) were evaluated.

3. Results

3.1. Fate of PAHs in a Chilean WWTP

Table 4 and Table 5 present the model results for the fate and distribution of naphthalene and benzo(a)pyrene throughout the treatment train. Results are presented for the 10%, 50%, and 95% probability percentiles of occurrence. For naphthalene in the water line, the influent load was 760.1 g/day, and the resulting effluent load was 238.9 g/day, predominantly in the dissolved phase. The overall removal efficiency of naphthalene from the water phase was 70.2%. In the sludge line, the resulting load to the biosolids (Figure 1, flow line 12) was 293.3 g/day, predominantly sorbed to particles. For benzo(a)pyrene, an influent load of 84.4 g/day resulted in an effluent load of 2.1 g/day, predominantly in the dissolved phase (97.7% overall removal). The resulting load to the biosolids was 80.9 g/day, predominantly sorbed to particles.

3.1.1. Distribution in the Water Line

Figure 3 and Figure 4 show the model results for the distribution of naphthalene and benzo(a)pyrene throughout the water line. Results are presented for the median input loads of the 95% percentile.
Considering each of the treatment units, the removal efficiencies were 22.9% for naphthalene and 56.5% for benzo(a)pyrene during primary treatment. For naphthalene, transport and sorption to sludge removed 15.8%, while volatilization removed 7.1%. In contrast, benzo(a)pyrene transport and sorption to sludge removed 55.8%, whereas volatilization was close to null (0.3%). The differences in the impact of sorption in the removal of naphthalene and benzo(a)pyrene can be explained by the solubility of PAHs. As their molecular weight increases, their solubility decreases; thus, they tend to adsorb at a higher rate [14]. For volatilization, the differences between the results obtained for naphthalene and benzo(a)pyrene can be explained by the differences in their Henry’s constant. Byrns et al., (2001) reported that compounds with log Kow between 3.3–5 present removal by sorption to sludge of 15–50%. They also reported that volatilization has a low impact on primary treatment removal rates, but is maximum for compounds with log Kow of 2 [44]. In the present study, results for naphthalene (log Kow = 3.36) and benzo(a)pyrene (log Kow = 5.97) agree with previous reports.
For naphthalene, 77% of the influent load went to secondary treatment units, most of which was in the dissolved phase, while 16% of the influent load went to the primary thickener, predominantly in the sorbed to particles phase. For benzo(a)pyrene, 43% of the influent load was transported to secondary treatment, while 56% was transported to primary thickener; in both cases, the sorbed to particle phase was predominant.
During secondary treatment, 64.3% of naphthalene and 95.7% of benzo(a)pyrene were removed. Once again, the predominant mechanism was adsorption to sludge for both PAHs. For naphthalene, volatilization removed 19.7% while biodegradation removed only 5.8%. In contrast, benzo(a)pyrene was hardly removed by biodegradation (0.4%) or volatilization (0.2%). Low biodegradation rates can be explained by two factors. First, naphthalene was more bioavailable than benzo(a)pyrene (61.3 g/day and 7.2 g/day in the dissolved phase, respectively). Second, as naphthalene is a shorter molecule (two benzene rings), its biodegradation rates are higher than benzo(a)pyrene (four-ringed molecule) [45].

3.1.2. Distribution in the Sludge Line

Figure 5 and Figure 6 show the model results for the distribution of naphthalene and benzo(a)pyrene along the sludge line. Results are presented for the median input loads of the 95% percentile.
Overall, naphthalene and benzo(a)pyrene removed in the sludge line (25.4% and 19.3%, respectively) were mostly returned to the water line. PAHs were not removed by the thickeners and dewatering units due to the model’s assumptions (See Table 3).
In the anaerobic digester, the removal efficiencies obtained for both PAHs were negligible (0.5% and 0.1%). Removal by volatilization for both PAHs was also negligible (less than 0.1% for both), in accordance with results obtained by other studies [36,37,46]. As for the biodegradation rates in the AD, naphthalene and benzo(a)pyrene were barely removed by this process (0.4% and 0.1%, respectively). These results can be explained because the influent load to AD (flow line 10, Figure 1) was predominantly sorbed to particles and, therefore, unavailable for biodegradation.

3.2. Sensitivity Analysis

Figure 7a,b present the results for the sensitivity analysis using the median of the 95th percentile influent load. The complete analysis is provided in Tables S8 and S9. Response variables include naphthalene loads in WWTP effluent and biosolids. Results are shown only for naphthalene because the model’s sensitivity to changes in benzo(a)pyrene parameters was similar to that to changes in naphthalene parameters (data not shown). Moreover, volatilization was more evident in the latter. Figure 7 shows that the parameter with greatest impact was the solid–liquid partitioning coefficient in the secondary treatment (Kd,sec). An increase in Kd,sec had the greatest impact on the load to the biosolids as the sorption to sludge increased. Furthermore, biodegradation and volatilization were negatively affected by less naphthalene present in the dissolved phase. Changes in SRT also affect the removal efficiencies described in the model. For a 50% increase in SRT, effluent loads of naphthalene decreased by approximately 12%. In contrast, a 50% decrease in SRT increasing the effluent load by approximately 19% (Table S8). Other studies also reported a direct correlation between SRT and PAH removal [9,46]; however, it is important to consider the feasibility of varying SRT during the operation. Increasing SRT decreases the frequency of sludge wasting, which potentially decreases the load of PAHs entering the AD unit. Longer SRTs could also promote changes in the mixed liquor (old sludge), thus promoting growth of undesired microorganisms and accumulation of extra polymeric substances. Changes in the chemical characteristics of the mixed liquor could potentially affect the partitioning coefficients (e.g., Kd,sec) and removal mechanisms of PAHs in the secondary treatment [47].
As for Kp and k1 (which determine sorption to particles during AD), results showed that these parameters impact naphthalene’s distribution in the sludge (e.g., Ffree, FDCM, and Fp), but do not alter the total load in the biosolids or in the water exiting the plant (Figure 7b). They could, however, impact the bioavailability and further biodegradation of naphthalene during anaerobic digestion.

4. Discussion

4.1. Fate of PAHs in a Chilean WWTP

Some clear differences could be observed in the behavior of the studied PAHs, mainly attributed to their physicochemical properties and associated partitioning, transport, and removal mechanisms in each treatment unit. Overall removal efficiencies for naphthalene and benzo(a)pyrene in the water line were 70.2% and 97.7% respectively, where volatilization removed 26.8% and 6.1% mostly in secondary treatment. Similar to other studies in WWTPs worldwide (Table 6), high sorption rates promoted their predominant transport to the sludge treatment units.
Biodegradation did not play an important role in the removal of PAHs in the modeled WWTP. When comparing the results obtained in the secondary treatment to other studies, Manoli et al., (2007) only obtained a 1–4% removal by biodegradation, while Byrns et al., (2001) obtained 80% removal for compounds with log Kow of 3.6 [1,44]. An important distinction, however, is that Byrns et al., (2001) considered both the dissolved and sorbed to particles phases as bioavailable. In this study, 36% of the naphthalene load in the activated sludge influent (Flow 2, Figure 1) was transported in the water line effluent (Flow 3, Figure 1), predominantly in the dissolved phase, while 39% was transported to secondary thickener (Flow 7, Figure 1). For benzo(a)pyrene, 4% of the influent load was transported in the water line effluent (Flow 2, Figure 1), while 95% was derived into the sludge line (Flow 7, Figure 1).
The predominant fate mechanism of PAHs in AD was desorption from sludge particles, while biodegradation was negligible. For example, even though the bioavailable fraction of naphthalene was 6.4% of the AD influent load, less than 1% was biodegraded. The remaining was sorbed to particles in this treatment stage and went back to the water line by desorption. Other available studies did not show a clear tendency in biodegradation of PAHs in AD. While some presented removal efficiencies on the order of 38% to 73% [36,48], other studies presented low or null removals in this stage [10,49,50]. Semblante et al. (2015) suggested that the inconsistency occurs because high removal rates are often obtained in pilot-scale studies [51].
In summary, for both naphthalene and benzo(a)pyrene, the predominant mechanism in the water line was sorption to sludge, while, in the sludge line, desorption was identified as the predominant mechanism. Higher removal efficiencies of benzo(a)pyrene in the water line by sorption to the sludge were expected due to its lower solubility in water. Compared to other studies in the literature, the model satisfactorily describes the mechanisms involved in the fate and distribution of PAHs in a WWTP [1,35,44].

4.2. Impacts on Reuse in Chile

Quality standards must be ensured in order to reuse both treated wastewater and biosolids; however, information is scarce regarding established regulatory standards worldwide. In Chile, PAHs are not included in the current regulation related to the disposal of treated effluents and/or biosolids [52]. For example, the DS 90/2000 regulates the discharge of treated effluent to marine and continental surface waters and does not include organic contaminants of any kind. Other related regulations, such as the NCh1333, which sets water quality requirements for different uses (including irrigation), only consider basic parameters for soil maintenance. Regarding reuse applications, the Ley 21.075/2018 regulates the collection and disposition of gray water in urban and rural areas, while it does not include emerging or persistent pollutants. In this case, international regulations/recommendations should be followed when considering water reuse or biosolids disposal projects in Chile.
Standards for PAH concentration in drinking water are variable across the world, but only a few organizations have set maximum concentrations for different micropollutants for direct potable reuse applications. For benzo(a)pyrene, WHO’s maximum concentration in water is 0.7 µg/L, which is based on their recommended drinking water standards [53]. The European Union, Council Directive 2008/105/EC on the quality of water intended for human consumption established a maximum limit for benzo(a)pyrene of 0.01 µg/L, and a maximum of 0.10 µg/L for the sum of benzo[b]fluoranthene, benzo(k)fluoranthene, benzo(g,h,i)perylene, and indenol(1,2,3-cd)pyrene [54]. The US federal government established regulatory standards and guidelines to protect people from potential health effects caused by eating, drinking, or breathing PAHs. Under the Safe Drinking Water Act, the US-EPA sets legal maximum limits on the level of benzo(a)pyrene in drinking water at 0.2 µg/L [55]. In Latin America, only a few countries have regulations for PAHs in water, specifically for benzo(a)pyrene: Brazil, with a maximum limit of 0.7 µg/L concentration in drinking water [56], and Argentina, with a maximum of 0.04 µg/L in water sources for human consumption [57].
In the case of Naphthalene, the US-EPA has defined a maximum concentration in surface runoff and in groundwater at a Health Reference Level (HRL) of 140 µg/L. Furthermore, the agency states that naphthalene is infrequently detected in public water supplies. When detected, naphthalene rarely exceeds the HRL or a value of one-half the HRL [55].
When considering other water reuse applications, such as irrigation of crops, no limit concentration of PAHs has been established, even though evidence suggests that PAH-contaminated water may lead to bioaccumulation of these compounds in crops [58,59].
A common practice worldwide is the disposal of stabilized biosolids for soil amendment in agriculture. Most international entities focus on establishing a limit concentration for heavy metals or pathogens and do not establish a maximum concentration of PAHs for soil amendment [60]. Only a few countries such as France and Italy have established maximum concentrations for benzo(a)pyrene and for the total concentration of PAHs as 2 μg/gDM and 6 μg/gDM, respectively [61].
Even though there is evidence of the presence of PAHs in Chilean urban centers [62,63,64,65,66,67], local regulations do not consider PAHs in the disposal of WWTP effluents, potential reuse scenarios, or biosolid use in agriculture. Monitoring of PAHs in both treated effluents and biosolids is imperative, especially when considering that de facto reuse in agricultural activities currently takes place downstream of the studied WWTP.
Table 6. Summary of PAHs presence in sludge/biosolids from literature.
Table 6. Summary of PAHs presence in sludge/biosolids from literature.
LocationSludge TypeNaphthalene [µg/kgDM]Benzo(a)pyrene [µg/kgDM]Reference
Thessaloniki, Greece2ary19066[10]
Mixed (1ary + 2ary)900140
Thickened340140
Stabilized800240
Dehydrated1.000240
Venice, Italy2ary3187[8]
Stabilized2886
Beijing, China2ary1.1005.000[68]
Wuxi, China50024
Dapu, China3107
Yuanlang, China120480
Shatian, China24033
GBD, China2ary98395[69]
BXH, China23327
Sibao, China2ary16.3204.540[70]
Ningbo, China1403.700
Spain2ary19734[71]
Spain27522
Lombardy, Italy1ary25956[50]
2ary27736
Thickened27944
Stabilized25462
Guadalete, Spain1ary776115[72]
2ary70787
Stabilized539170
This study1ary1.270271
2ary2.140378
Thickened (1ary + 2ary)1.646318
Stabilized2.692503
Dehydrated2.658498

4.3. Limitations

The results of this study present a first approach to quantifying the presence of PAHs in a Chilean WWTP. However, as experimental samples were unavailable, it is recommended to validate the model with in situ measurements. It is important to notice that the mechanistic model considers only influent PAHs loads in the dissolved phase. However, evidence suggests that, as urban wastewater takes hours to reach the WWTP, a sorption equilibrium condition is likely to be reached and PAHs could be found in the particulate phase [10]. Further studies should incorporate the particulate phase in the influent load.
The results also highlight the importance of estimating (or measuring) relevant parameters considered in the model (e.g., Kd), as they show how possible variations in them can impact the obtained concentrations. Calibration of these parameters is also recommended using current measured data.

5. Conclusions

The present study provides a modeling tool to estimate the presence and fate of PAHs in a conventional activated sludge WWTP in Chile, including both water and sludge lines. Despite the inherent variability of the treated wastewater, represented here by an influent concentration of 95% probability of occurrence, the model satisfactorily describes the main partitioning mechanism of PAHs. For instance, adsorption to sludge was the predominant mechanism in the water line, while desorption from particles was identified as the predominant mechanism in the sludge line. Clear differences could be established between naphthalene and benzo(a)pyrene mainly due to the differences in their physicochemical properties. Even though biodegradation and volatilization had negligible influence in the removal of benzo(a)pyrene during the treatment process, volatilization was an important removal mechanism of naphthalene in the secondary treatment. Further studies should include current in situ measurements to validate the model for the Chilean scenario.
Chilean regulatory standards for the disposal of WWTP’s effluents or biosolids do not consider maximum concentration for PAHs (or any other micropollutant). However, PAH’s presence in treated wastewater and biosolids should be monitored in order to ensure human and environmental safety, as evidence suggests that PAH accumulation in wastewater-irrigated soil and crops can occur. Considering the growing market for potable and non-potable reuse of treated effluents in Chile, this topic should be a matter of urgent concern.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/pr9081313/s1, Table S1: Influent naphthalene loads (g/d) obtained through the bootstrap for the 10%, 50% and 95% probability percentiles, Table S2: Influent benzo(a)pyrene loads (g/d) obtained through the bootstrap for the 10%, 50% and 95% probability percentiles, Table S3: Assumptions for each treatment stage in the PAH model, Table S4: Rate equations for the processes considered in the water line, Table S5: Parameters used in the water line, Table S6: Rate equations for the processes considered in the sludge line, Table S7: Parameters used in the sludge line, Table S8: Variations in the PAH load present in the effluent for variations of ±25%, ±50%, and ±75%, Table S9: Variations in the PAH load present in the biosolids for variations of ±25%, ±50%, and ±75%.

Author Contributions

A.A., methodology, investigation, formal analysis, and writing—original draft preparation; F.P., methodology, investigation, and formal analysis; M.R., methodology, investigation, and resources; A.L.P., conceptualization, data curation, writing—review and editing, supervision, project administration, and funding acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded in part by ANID PCI REDI170494 and by the Universidad de Chile, Vicerrectoría de Investigación y Desarrollo.

Acknowledgments

The authors acknowledge Nestor Rojas for providing background information on the plant’s original design and ASSM model considerations.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Manoli, E.; Samara, C. The removal of Polycyclic Aromatic Hydrocarbons in the wastewater treatment process: Experimental calculations and model predictions. Environ. Pollut. 2008, 151, 477–485. [Google Scholar] [CrossRef]
  2. Haritash, A.; Kaushik, C. Biodegradation aspects of Polycyclic Aromatic Hydrocarbons (PAHs): A review. J. Hazard. Mater. 2009, 169, 1–15. [Google Scholar] [CrossRef]
  3. Alegbeleye, O.O.; Opeolu, B.O.; Jackson, V.A. Polycyclic Aromatic Hydrocarbons: A Critical Review of Environmental Oc-currence and Bioremediation. Environ. Manag. 2017, 60, 758–783. [Google Scholar] [CrossRef] [PubMed]
  4. U.S. Environmental Protection Agency. Priority Pollutants. (EPA U.S.). 2011. Available online: http://water.epa.gov/scitech/methods/cwa/pollutants.cfm (accessed on 21 May 2017).
  5. Directive 2013/39/EU of the European Parliament and of the Council of 12 August 2013 Amending Directives 2000/60/EC and 2008/105/EC as Regards Priority Substances in the Field of Water Policy Text with EEA Relevance. Available online: https://eur-lex.europa.eu/LexUriServ/LexUriServ.do?uri=OJ:L:2013:226:0001:0017:EN:PDF (accessed on 6 June 2020).
  6. Qiao, M.; Qi, W.; Liu, H.; Qu, J. Occurrence, behavior and removal of typical substituted and parent polycyclic aromatic hydrocarbons in a biological wastewater treatment plant. Water Res. 2014, 52, 11–19. [Google Scholar] [CrossRef]
  7. Zhang, X.; Yu, T.; Li, X.; Yao, J.; Liu, W.; Chang, S.; Chen, Y. The fate and enhanced removal of polycyclic aromatic hydrocarbons in wastewater and sludge treatment system: A review. Crit. Rev. Environ. Sci. Technol. 2019, 49, 1425–1475. [Google Scholar] [CrossRef]
  8. Busetti, F.; Heitz, A.; Cuomo, M.; Badoer, S.; Traverso, P. Determination of sixteen polycyclic aromatic hydrocarbons in aqueous and solid samples from an Italian wastewater treatment plant. J. Chromatogr. A 2006, 1102, 104–115. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Fatone, F.; Di Fabio, S.; Bolzonella, D.; Cecchi, F. Fate of aromatic hydrocarbons in Italian municipal wastewater systems: An overview of wastewater treatment using conventional activated-sludge processes (CASP) and membrane bioreactors (MBRs). Water Res. 2011, 45, 93–104. [Google Scholar] [CrossRef]
  10. Manoli, E.; Samara, C. Occurrence and Mass Balance of Polycyclic Aromatic Hydrocarbons in the Thessaloniki Sewage Treatment Plant. J. Environ. Qual. 1999, 28, 176–187. [Google Scholar] [CrossRef]
  11. Pham, T.; Prolux, S. PCBs and PAHs in the Montreal urban community (Quebec, Canada) wastewater treatment plant and in the effluent plume in the St. Lawrence River. Water Res. 1997, 31, 1887–1896. [Google Scholar] [CrossRef]
  12. Sun, S.; Jia, L.; Li, B.; Yuan, A.; Kong, L.; Qi, H.; Ma, W.; Zhang, A.; Wu, Y. The occurrence and fate of PAHs over multiple years in a wastewater treatment plant of Harbin, Northeast China. Sci. Total Environ. 2018, 624, 491–498. [Google Scholar] [CrossRef]
  13. Dionisi, D.; Bertin, L.; Bornoroni, L.; Capodicasa, S.; Papini, M.P.; Fava, F. Removal of organic xenobiotics in activated sludges under aerobic conditions and anaerobic digestion of the adsorbed species. J. Chem. Technol. Biotechnol. 2006, 81, 1496–1505. [Google Scholar] [CrossRef]
  14. Jelic, A.; Katsou, E.; Malamis, S.; Bolzonella, D.; Fatone, F. Occurence, Removal, and Fate of PAHs and VOCs in Municipal Wastewater Treatment Plants: A Literature Review. In Wastewater Treatment: Occurrence and Fate of Polycylic Aromatic Hydro-carbons (PAHs); Forsgren, A., Ed.; CRC Press: Boca Raton, FL, USA, 2015; pp. 91–112. [Google Scholar]
  15. Mozo, I.; Bounouba, M.; Mengelle, E.; Lesage, N.; Sperandio, M.; Bessiere, Y. Modelling PAHs removal in activated sludge process: Effect of disintegration. Water Sci. Technol. 2019, 80, 794–805. [Google Scholar] [CrossRef] [PubMed]
  16. Margot, J.; Rossi, L.; Barry, D.; Holliger, C. A review of the fate of micropollutants in wastewater treatment plants. Wiley Interdiscip. Rev. Water 2015, 2, 457–487. [Google Scholar] [CrossRef] [Green Version]
  17. Liu, Z.; Li, Q.; Wu, Q.; Kuo, D.T.F.; Chen, S.; Hu, X.; Deng, M.; Zhang, H.; Luo, M. Removal Efficiency and Risk Assessment of Polycyclic Aromatic Hydrocarbons in a Typical Municipal Wastewater Treatment Facility in Guangzhou, China. Int. J. Environ. Res. Public Health 2017, 14, 861. [Google Scholar] [CrossRef] [Green Version]
  18. Clouzot, L.; Choubert, J.-M.; Cloutier, F.; Goel, R.; Love, N.G.; Melcer, H.; Ort, C.; Patureau, D.; Plosz, B.; Pomiès, M.; et al. Perspectives on modelling micropollutants in wastewater treatment plants. Water Sci. Technol. 2013, 68, 448–461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Plosz, B.; Benedetti, L.; Daigger, G.T.; Langford, K.H.; Larsen, H.F.; Monteith, H.; Ort, C.; Seth, R.; Steyer, J.-P.; Vanrolleghem, P. Modelling micropollutant fate in wastewater collection and treatment systems: Status and challenges. Water Sci. Technol. 2013, 67, 1–15. [Google Scholar] [CrossRef]
  20. Komolafe, O.; Mrozik, W.; Dolfing, J.; Acharya, K.; Vassalle, L.; Mota, C.; Davenport, R. Occurrence and removal of micropollu-tants in full-scale aerobic, anaerobic and facultative wastewater treatment plants in Brazil. J. Environ. Manag. 2021, 287, 112286. [Google Scholar] [CrossRef]
  21. Nas, B.; Argun, M.; Dolu, T.; Ateş, H.; Yel, E.; Koyuncu, S.; Dinç, S.; Kara, M. Occurrence, loadings and removal of EU-priority polycyclic aromatic hydrocarbons (PAHs) in wastewater and sludge by advanced biological treatment, stabilization pond and constructed wetland. J. Environ. Manag. 2020, 268, 110580. [Google Scholar] [CrossRef]
  22. Wang, X.; Xi, B.; Huo, S.; Sun, W.; Pan, H.; Zhang, J.; Ren, Y.; Liu, H. Characterization, treatment and releases of PBDEs and PAHs in a typical municipal sewage treatment plant situated beside an urban river, East China. J. Environ. Sci. 2013, 25, 1281–1290. [Google Scholar] [CrossRef]
  23. Tian, W.; Bai, J.; Liu, K.; Sun, H.; Zhao, Y. Occurrence and removal of polycyclic aromatic hydrocarbons in the wastewater treatment process. Ecotoxicol. Environ. Saf. 2012, 82, 1–7. [Google Scholar] [CrossRef]
  24. Vogelsang, C.; Grung, M.; Jantsch, T.G.; Tollefsen, K.E.; Liltved, H. Occurrence and removal of selected organic micropollu-tants at mechanical, chemical and advanced wastewater treatment plants in Norway. Water Res. 2006, 40, 3559–3570. [Google Scholar] [CrossRef]
  25. Taboada-Santos, A.; Behera, C.R.; Sin, G.; Gernaey, K.V.; Mauricio-Iglesias, M.; Carballa, M.; Lema, J.M. Assessment of the fate of organic micropollutants in novel wastewater treatment plant configurations through an empirical mechanistic model. Sci. Total Environ. 2020, 716, 137079. [Google Scholar] [CrossRef]
  26. Lejeune, A.; Choubert, J.-M. Modeling of micropollutant removal in full-scale membrane bioreactors: Calibration and operations to limit the emissions. Bioprocess Biosyst. Eng. 2019, 42, 1879–1892. [Google Scholar] [CrossRef] [PubMed]
  27. Gonzalez-Gil, L.; Mauricio-Iglesias, M.; Carballa, M.; Lema, J.M. Why are organic micropollutants not fully biotransformed? A mechanistic modelling approach to anaerobic systems. Water Res. 2018, 142, 115–128. [Google Scholar] [CrossRef]
  28. Baalbaki, Z.; Torfs, E.; Yargeau, V.; Vanrolleghem, P.A. Predicting the fate of micropollutants during wastewater treatment: Calibration and sensitivity analysis. Sci. Total Environ. 2017, 601, 874–885. [Google Scholar] [CrossRef] [PubMed]
  29. Polesel, F.; Andersen, H.R.; Trapp, S.; Plosz, B. Removal of Antibiotics in Biological Wastewater Treatment Systems—A Critical Assessment Using the Activated Sludge Modeling Framework for Xenobiotics (ASM-X). Environ. Sci. Technol. 2016, 50, 10316–10334. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Plósz, B.G.; Langford, K.H.; Thomas, K.V. An activated sludge modeling framework for xenobiotic trace chemicals (ASM-X): Assessment of diclofenac and carbamazepine. Biotechnol. Bioeng. 2012, 109, 2757–2769. [Google Scholar] [CrossRef] [PubMed]
  31. Vezzaro, L.; Benedetti, L.; Gevaert, V.; De Keyser, W.; Verdonck, F.; De Baets, B.; Nopens, I.; Cloutier, F.; Vanrolleghem, P.; Mikkelsen, P.S. A model library for dynamic transport and fate of micropollutants in integrated urban wastewater and stormwater systems. Environ. Model. Softw. 2014, 53, 98–111. [Google Scholar] [CrossRef]
  32. Fernandez-Fontaina, E.; Carballa, M.; Omil, F.; Lema, J. Modelling cometabolic biotransformation of organic micropollutants in nitrifying reactors. Water Res. 2014, 65, 371–383. [Google Scholar] [CrossRef]
  33. Lust, M.; Makinia, J.; Stensel, H.D. A mechanistic model for fate and removal of estrogens in biological nutrient removal activated sludge systems. Water Sci. Technol. 2012, 65, 1130–1136. [Google Scholar] [CrossRef]
  34. Dionisi, D.; Bornoroni, L.; Mainelli, S.; Majone, M.; Pagnanelli, F.; Papini, M.P. Theoretical and Experimental Analysis of the Role of Sludge Age on the Removal of Adsorbed Micropollutants in Activated Sludge Processes. Ind. Eng. Chem. Res. 2008, 47, 6775–6782. [Google Scholar] [CrossRef]
  35. Lee, K.-C.; Rittmann, B.E.; Shi, J.; McAvoy, D. Advanced steady-state model for the fate of hydrophobic and volatile com-pounds in activated sludge. Water Environ. Res. 1998, 70, 1118–1131. [Google Scholar] [CrossRef]
  36. Delgadillo-Mirquez, L.; Lardon, L.; Steyer, J.-P.; Patureau, D. A new dynamic model for bioavailability and cometabolism of micropollutants during anaerobic digestion. Water Res. 2011, 45, 4511–4521. [Google Scholar] [CrossRef]
  37. Barret, M.; Patureau, D.; Latrille, E.; Carrère, H. A three-compartment model for micropollutants sorption in sludge: Methodological approach and insights. Water Res. 2010, 44, 616–624. [Google Scholar] [CrossRef]
  38. Suez Degremont Water Handbook. Available online: https://www.suezwaterhandbook.com/case-studies/wastewater-treatment/La-Farfana-wastewater-treatment-plant-Chile (accessed on 25 May 2021).
  39. Superintendencia de Servicios Sanitarios (SISS). Estudio Grandes Plantas de Tratamiento de Aguas Servidas Aguas Andinas Proceso Tarifario 2015–2020. 2014. Available online: https://infota.siss.cl/concesiones/index.php?dir=empresas%2FAguasAndinas%2F09+Estudio+Tarifario+Intercambio+SISS%2F (accessed on 24 May 2020).
  40. Langford, K.H.; Lester, J.N. Fate and Behavior of Endocrine Disrupters in Wastewater Treatment Processes. In Endocrine Disrupters in Wastewater and Sludge Treatment Processes, 1st ed.; Birkett, J.W., Lester, J.N., Eds.; Lewis Publishers: London, UK, 2003; Volume 1, pp. 103–144. [Google Scholar]
  41. Meakins, N.; Bubb, J.; Lester, J. The fate and behaviour of organic micropollutants during wastewater treatment processes: A review. Int. J. Environ. Pollut. 1994, 4, 27–58. [Google Scholar]
  42. Plosz, B.; Leknes, H.; Thomas, K.V. Impacts of Competitive Inhibition, Parent Compound Formation and Partitioning Behavior on the Removal of Antibiotics in Municipal Wastewater Treatment. Environ. Sci. Technol. 2009, 44, 734–742. [Google Scholar] [CrossRef] [PubMed]
  43. Nanning, J.; Mena, M.P.; Rojas, N. ASSM (Activated Sludge Simple Model). In Proceedings of the XXXI Congreso Interamericano AIDIS, Santiago, Chile, 12–15 October 2015. [Google Scholar]
  44. Byrns, G. The fate of xenobiotic organic compounds in wastewater treatment plants. Water Res. 2001, 35, 2523–2533. [Google Scholar] [CrossRef]
  45. Bautista, L.F.; Sanz, R.; Molina, M.C.; González, N.; Sánchez, D. Effect of different non-ionic surfactants on the biodegradation of PAHs by diverse aerobic bacteria. Int. Biodeterior. Biodegrad. 2009, 63, 913–922. [Google Scholar] [CrossRef]
  46. Sponza, D.T.; Gok, O. Aerobic biodegradation and inhibition kinetics of poly-aromatic hydrocarbons (PAHs) in a petrochemical industry wastewater in the presence of biosurfactants. J. Chem. Technol. Biotechnol. 2012, 87, 658–672. [Google Scholar] [CrossRef]
  47. Semblante, G.U.; Hai, F.I.; Huang, X.; Ball, A.; Price, W.E.; Nghiem, L. Trace organic contaminants in biosolids: Impact of conventional wastewater and sludge processing technologies and emerging alternatives. J. Hazard. Mater. 2015, 300, 1–17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Geyik, A.G.; Kılıç, B.; Çeçen, F. Extracellular polymeric substances (EPS) and surface properties of activated sludges: Effect of organic carbon sources. Environ. Sci. Pollut. Res. 2015, 23, 1653–1663. [Google Scholar] [CrossRef] [PubMed]
  49. Torretta, V.; Katsoyiannis, A. Occurrence of polycyclic aromatic hydrocarbons in sludges from different stages of a wastewater treatment plant in Italy. Environ. Technol. 2013, 34, 937–943. [Google Scholar] [CrossRef] [PubMed]
  50. Wlodarczyk-Makula, M. The loads of PAHS in wastewater and sewage sludge of municipal treatment plant. Polycycl. Aromat. Compd. 2005, 25, 183–194. [Google Scholar] [CrossRef]
  51. Trably, E.; Patureau, D.; Delgenes, J. Enhancement of polycyclic aromatic hydrocarbons removal during anaerobic treatment of urban sludge. Water Sci. Technol. 2003, 48, 53–60. [Google Scholar] [CrossRef] [PubMed]
  52. Biblioteca del Congreso. Available online: https://www.bcn.cl/leychile/ (accessed on 30 May 2021).
  53. WHO (2017) World Health Organization. Guidance for Producing Safe Drinking-Water. Available online: http://apps.who.int/iris/bitstream/handle/10665/258715/9789241512770-eng.pdf?sequence=1&isAllowed=y (accessed on 17 July 2020).
  54. Consolidated Text: Directive 2008/105/EC of the European Parliament and of the Council of 16 December 2008 on Environmental Quality Standards in the Field of Water Policy, Amending and Subsequently Repealing Council Directives 82/176/EEC, 83/513/EEC, 84/156/EEC, 84/491/EEC, 86/280/EEC and Amending Directive 2000/60/EC of the European Parliament and of the Council. Available online: https://eur-lex.europa.eu/eli/dir/2008/105/oj (accessed on 17 July 2020).
  55. US-EPA. Safe Drinking Water Act. 2020. Available online: https://www.epa.gov/sdwa (accessed on 30 May 2021).
  56. Ministerio de Salud de Brasil. Ordenanza nº 518: Establece los Procedimientos y Responsabilidades Relativos al Control y la Vigilancia de la Calidad del Agua Para Consumo Humano y su Norma de Potabilidad, y Establece Otras Disposiciones. 2004. Available online: http://189.28.128.100/dab/docs/legislacao/portaria518_25_03_04.pdf (accessed on 18 August 2020).
  57. Subsecretaria de Recursos Hídricos de Argentina. Niveles Guía Nacionales de Calidad de Agua Ambiente: Ben-zo[a]pireno. 2003. Available online: https://www.argentina.gob.ar/sites/default/files/documento22.pdf (accessed on 18 August 2020).
  58. Khan, S.; Aijun, L.; Zhang, S.; Hu, Q.; Zhu, Y.-G. Accumulation of polycyclic aromatic hydrocarbons and heavy metals in lettuce grown in the soils contaminated with long-term wastewater irrigation. J. Hazard. Mater. 2008, 152, 506–515. [Google Scholar] [CrossRef] [PubMed]
  59. Samia, K.; Abdelkarim, C.; Yassin, H.; Teruo, H. The effect of long-term soil irrigation by wastewater on organic matter, polycyclic aromatic hydrocarbons, and heavy metals evolution: Case study of Zaouit Sousse (Tunisia). Arab. J. Geosci. 2012, 6, 4337–4346. [Google Scholar] [CrossRef]
  60. Forsgren, A.J. PAHs in Sewage Sludge, Soils and Sediments. In Wastewater Treatment: Occurrence and Fate of Polycylic Aromatic Hydrocarbons (PAHs); Forsgren, A., Ed.; CRC Press: Boca Raton, FL, USA, 2015; pp. 212–229. [Google Scholar]
  61. Collivignarelli, M.; Abbà, A.; Frattarola, A.; Miino, M.C.; Padovani, S.; Katsoyiannis, I.; Torretta, V. Legislation for the Reuse of Biosolids on Agricultural Land in Europe: Overview. Sustainability 2019, 11, 6015. [Google Scholar] [CrossRef] [Green Version]
  62. Bravo-Linares, C.; Ovando-Fuentealba, L.; Mudge, S.M.; Cerpa, J.; Loyola-Sepulveda, R. Source Allocation of Aliphatic and Polycyclic Aromatic Hydrocarbons in Particulate-Phase (PM10) in the City of Valdivia, Chile. Polycycl. Aromat. Compd. 2012, 32, 390–407. [Google Scholar] [CrossRef]
  63. Sienra, M.D.R.; Rosazza, N.G.; Préndez, M. Polycyclic aromatic hydrocarbons and their molecular diagnostic ratios in urban atmospheric respirable particulate matter. Atmos. Res. 2005, 75, 267–281. [Google Scholar] [CrossRef]
  64. Muñoz, C.; Droppelmann, A.; Erazo, M.; Aceituno, P.; Orellana, C.; Parro, J.; Mesias, S.; Marchetti, N.; Navas-Acien, A.; Iglesias, V. Occupational exposure to polycyclic aromatic hydrocarbons: A cross-sectional study in bars and restaurants in Santiago, Chile. Am. J. Ind. Med. 2016, 59, 887–896. [Google Scholar] [CrossRef] [Green Version]
  65. Pozo, K.; Estellano, V.H.; Harner, T.; Diaz-Robles, L.; Cereceda-Balic, F.; Etcharren, P.; Pozo, K.; Vidal, V.; Guerrero, F.; Vergara-Fernández, A. Assessing Polycyclic Aromatic Hydrocarbons (PAHs) using passive air sampling in the atmosphere of one of the most wood-smoke-polluted cities in Chile: The case study of Temuco. Chemosphere 2015, 134, 475–481. [Google Scholar] [CrossRef]
  66. Pozo, K.; Harner, T.; Rudolph, A.; Oyola, G.; Estellano, V.H.; Ahumada-Rudolph, R.; Garrido, M.; Pozo, K.; Mabilia, R.; Focardi, S. Survey of persistent organic pollutants (POPs) and polycyclic aromatic hydrocarbons (PAHs) in the atmosphere of rural, urban and industrial areas of Concepción, Chile, using passive air samplers. Atmos. Pollut. Res. 2012, 3, 426–434. [Google Scholar] [CrossRef] [Green Version]
  67. Bravo-Linares, C.; Ovando-Fuentealba, L.; Orellana-Donoso, S.; Villenas-Fernández, K.; Hernández-Millán, M.; Mudge, S.; Pinaud-Mendoza, J.P.; Loyola-Sepulveda, R. Source Apportionment of PAHs in Airborne Particulates (PM2.5) in Southern Chile. Polycycl. Aromat. Compd. 2016, 37, 189–202. [Google Scholar] [CrossRef]
  68. Cai, Q.-Y.; Mo, C.-H.; Wu, Q.-T.; Zeng, Q.-Y.; Katsoyiannis, A. Occurrence of organic contaminants in sewage sludges from eleven wastewater treatment plants, China. Chemosphere 2007, 68, 1751–1762. [Google Scholar] [CrossRef]
  69. Dai, J.; Xu, M.; Chen, J.; Yang, X.; Ke, Z. PCDD/F, PAH and heavy metals in the sewage sludge from six wastewater treatment plants in Beijing, China. Chemosphere 2007, 66, 353–361. [Google Scholar] [CrossRef]
  70. Hua, L.; Wu, W.-X.; Liu, Y.-X.; Tientchen, C.; Chen, Y.-X. Heavy Metals and PAHs in Sewage Sludge from Twelve Wastewater Treatment Plants in Zhejiang Province. Biomed. Environ. Sci. 2008, 21, 345–352. [Google Scholar] [CrossRef]
  71. Perez, S.; Guillamón, M.; Barceló, D. Quantitative analysis of polycyclic aromatic hydrocarbons in sewage sludge from wastewater treatment plants. J. Chomatography A 2001, 938, 57–65. [Google Scholar] [CrossRef]
  72. Villar, P.; Villar, M.; Callejón, M.; Pérez, J.; Jiménez, J.; Alonso, E.; Guiraúm, A. Evolution of polycyclic aromatic hydrocar-bons (PAHs) and heavy metals in sludge samples from conventional activated sludge wastewater treatment plants. Environmetrics 2009, 20, 561–574. [Google Scholar] [CrossRef]
Figure 1. Schematic of the model WWTP. Data correspond to flow rates (Q) and total suspended solids (TSS) concentration between treatment stages. The black arrows represent the water line (wastewater treatment units), the gray arrows represent the sludge line (sludge treatment units), and the dashed arrows represent return flows.
Figure 1. Schematic of the model WWTP. Data correspond to flow rates (Q) and total suspended solids (TSS) concentration between treatment stages. The black arrows represent the water line (wastewater treatment units), the gray arrows represent the sludge line (sludge treatment units), and the dashed arrows represent return flows.
Processes 09 01313 g001
Figure 2. Input loads of (a) naphthalene and (b) benzo(a)pyrene for 10%, 50%, and 95% probability percentiles.
Figure 2. Input loads of (a) naphthalene and (b) benzo(a)pyrene for 10%, 50%, and 95% probability percentiles.
Processes 09 01313 g002
Figure 3. Fate of naphthalene in the water line. Results are presented in g/day. Fdis is the resulting load in the dissolved phase, Fad is sorbed to particles in the water line, Ffree is freely dissolved, FDCM is sorbed to the dissolved and colloidal matter (DCM), and Fp is sorbed to particles going to the sludge line.
Figure 3. Fate of naphthalene in the water line. Results are presented in g/day. Fdis is the resulting load in the dissolved phase, Fad is sorbed to particles in the water line, Ffree is freely dissolved, FDCM is sorbed to the dissolved and colloidal matter (DCM), and Fp is sorbed to particles going to the sludge line.
Processes 09 01313 g003
Figure 4. Fate of benzo(a)pyrene in the water line. Results are presented in g/day. Fdis is the resulting load in the dissolved phase, Fad is sorbed to particles in the water line, Ffree is freely dissolved, FDCM is sorbed to the dissolved and colloidal matter (DCM), and Fp is sorbed to particles going to the sludge line.
Figure 4. Fate of benzo(a)pyrene in the water line. Results are presented in g/day. Fdis is the resulting load in the dissolved phase, Fad is sorbed to particles in the water line, Ffree is freely dissolved, FDCM is sorbed to the dissolved and colloidal matter (DCM), and Fp is sorbed to particles going to the sludge line.
Processes 09 01313 g004
Figure 5. Fate of naphthalene in the sludge line. Results are presented in g/day. Fdis is the resulting load in the dissolved phase, Ffree is the freely dissolved, FDCM is sorbed to the dissolved and colloidal matter (DCM), and Fp is sorbed to particles.
Figure 5. Fate of naphthalene in the sludge line. Results are presented in g/day. Fdis is the resulting load in the dissolved phase, Ffree is the freely dissolved, FDCM is sorbed to the dissolved and colloidal matter (DCM), and Fp is sorbed to particles.
Processes 09 01313 g005
Figure 6. Fate of benzo(a)pyrene in the sludge line. Results are presented in g/day. Fdis is the resulting load in the dissolved phase, Ffree is the freely dissolved, FDCM is sorbed to the dissolved and colloidal matter (DCM), and Fp is sorbed to particles.
Figure 6. Fate of benzo(a)pyrene in the sludge line. Results are presented in g/day. Fdis is the resulting load in the dissolved phase, Ffree is the freely dissolved, FDCM is sorbed to the dissolved and colloidal matter (DCM), and Fp is sorbed to particles.
Processes 09 01313 g006
Figure 7. Sensitivity analysis for a 50% increase in Kd,prim, Kd,sec, k1, Kp, and SRT. Response variables include naphthalene loads in WWTP (a) effluent and (b) biosolids.
Figure 7. Sensitivity analysis for a 50% increase in Kd,prim, Kd,sec, k1, Kp, and SRT. Response variables include naphthalene loads in WWTP (a) effluent and (b) biosolids.
Processes 09 01313 g007
Table 4. Fate of naphthalene in the modeled WWTP. Results are presented in g/day for 10%, 50%, and 95% percentiles. Combined values for flows 1, 2, 3, 6, 9, and 13 correspond to the aqueous phase (See Figure 1). Data are shown as min–max; median.
Table 4. Fate of naphthalene in the modeled WWTP. Results are presented in g/day for 10%, 50%, and 95% percentiles. Combined values for flows 1, 2, 3, 6, 9, and 13 correspond to the aqueous phase (See Figure 1). Data are shown as min–max; median.
Flow10%50%95%
FfreeFDCMFPFfreeFDCMFPFfreeFDCMFP
1170.7–221; 191.4---397.9–454; 426.4---454–795.4; 760.1---
2123.8–160.3; 138.914.7–19.1; 16.5288.7–329.4; 309.434.3–39.2; 36.8329.4–577.1; 551.539.2–68.6; 8.2
351.5–66.6; 57.72.2–2.9; 2.5119.9–136.8; 128.55.1–5.9; 5.5136.8–239.8; 229.15.9–10.3; 9.8
41.6–2.0; 1.74–5.2; 4.522.9–29.6; 25.63.7–4.2; 3.99.4–10.7; 10.153.3–60.8; 57.16.6–7.3; 716.9–18.8; 17.995.8–106.5; 101.8
50.1–0.1; 0.13.1–4.0; 3.521–27.2; 23.60.3–0.3; 0.37.2–8.2; 7.749–55.9; 52.50.5–0.5; 0.512.9–14.4; 13.788.1–97.9; 93.6
62.4–3.1; 2.71.8–2.4; 25.6–6.4; 64.3–4.9; 4.610.1–11.2; 10.77.7–8.6; 8.2
70.6–0.8; 0.72.3–3.0; 2.656.1–72.6; 62.91.4–1.5; 1.55.3–6.1; 5.7130.7–149.2; 140.12.4–2.7; 2.69.6–10.7; 10.2235–261.4; 249.8
8---1.4–1.9; 1.645.8–59.3; 51.4---3.4–3.8; 3.6106.8–121.9 114.4---6.1–6.7; 6.4192–213.5; 204
91.4–1.8; 1.610.3–13.3; 11.53.2–3.7; 3.523.9–27.3; 25.75.8–6.5; 6.243.1–47.9; 45.8
100.2–0.2; 0.24.4–5.6; 4.966.3–85.8; 74.90.4–0.4; 0.410.1–11.6; 10.9154.4–176.2; 165.60.7–0.7; 0.718.2–20.3; 19.4277.6–308.7; 295
111.2–1.6; 1.45.8–7.5; 6.563.4–82.1; 71.12.9–3.3; 3.113.6–15.5; 14.5147.7–168.6; 158.35.2–5.8; 5.524.4–27.1; 25.9265.6–295.3; 282.2
120.1–0.1; 0.15.3–6.9; 5.959.8–77.4; 670.2–0.3; 0.212.4–14.1; 13.3139.4–159.1; 149.40.4–0.5; 0.422.2–24.7; 23.6250.6–278.7; 266.3
131.6–2.1; 1.83.2–4.1; 3.63.8–4.3; 4.17.4–8.4; 7.96.8–7.6; 7.213.3–14.8; 14.1
Table 5. Fate of benzo(a)pyrene in the modeled WWTP. Results are presented in g/day for 10%, 50%, and 95% percentiles. Combined values for flows 1, 2, 3, 6, 9, and 13 correspond to the aqueous phase (See Figure 1). Data shown as min–max; median.
Table 5. Fate of benzo(a)pyrene in the modeled WWTP. Results are presented in g/day for 10%, 50%, and 95% percentiles. Combined values for flows 1, 2, 3, 6, 9, and 13 correspond to the aqueous phase (See Figure 1). Data shown as min–max; median.
Flow10%50%95%
FfreeFDCMFPFfreeFDCMFPFfreeFDCMFP
14.4–10.2; 6.2---24.9–33.9; 28.4---82.9–86.2; 84.4---
20.4–0.9; 0.50–4; 2.52.1–2.9; 2.49.9–13.4; 11.27.1–7.4; 7.232.8–34.1; 33.4
31.7–0.2; 0.1---0.5–0.7; 0.60.1–0.1; 0.11.7–1.8; 1.80.3–0.3; 0.3
4---2.7–6.3; 3.8---15.3–20.8; 17.4---50.9–53; 51.9------2.7–6.3; 4.2
5---2.5–5.8; 3.5---14.1–19.2; 16.10–0.1; 0.146.9–48.7; 47.7------2.5–5.8; 3.5
6---0.2–0.5; 0.30–0.1; 0.11.2–1.7; 1.40.1–0.1; 04.1–4.3; 4.2
7---2.4–5.6; 3.4---13.8–18.9; 12.8---45.9–47.7; 46.7------2.4–5.6; 3.4
8---2–4.6; 2.8---11.3–15.3; 12.8---37.5–39; 38.1------2–4.6; 2.8
9---0.4–1; 0.60–0.1; 0.12.5–3.4; 2.9---8.4–8.7; 8.6
10---4.4–10.3; 6.2---25–34.1; 28.6---83.4–86.7; 85.9------4.4–10.3; 6.2
11---4.1–9.5; 5.80.1–0.1; 023.2–31.6; 26.40.3–0.3; 0.377.2–80.3; 79.6---0.3–0.8; 0.54.1–9.5; 5.8
12---3.9–9; 5.5---21.9–29.8; 25---72.9–75.8; 75.1---0.3–0.7; 0.43.9–9; 5.5
130–0.1; 00.2–0.5; 0.30.3–0.3; 0.31.2–1.6; 1.30.8–0.8; 0.83.9–4; 4
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Arroyo, A.; Provoste, F.; Rodríguez, M.; Prieto, A.L. A Mechanistic Model to Assess the Fate of Naphthalene and Benzo(a)pyrene in a Chilean WWTP. Processes 2021, 9, 1313. https://doi.org/10.3390/pr9081313

AMA Style

Arroyo A, Provoste F, Rodríguez M, Prieto AL. A Mechanistic Model to Assess the Fate of Naphthalene and Benzo(a)pyrene in a Chilean WWTP. Processes. 2021; 9(8):1313. https://doi.org/10.3390/pr9081313

Chicago/Turabian Style

Arroyo, Antonia, Fabián Provoste, Montserrat Rodríguez, and Ana L. Prieto. 2021. "A Mechanistic Model to Assess the Fate of Naphthalene and Benzo(a)pyrene in a Chilean WWTP" Processes 9, no. 8: 1313. https://doi.org/10.3390/pr9081313

APA Style

Arroyo, A., Provoste, F., Rodríguez, M., & Prieto, A. L. (2021). A Mechanistic Model to Assess the Fate of Naphthalene and Benzo(a)pyrene in a Chilean WWTP. Processes, 9(8), 1313. https://doi.org/10.3390/pr9081313

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop