Determining Discharge Characteristics and Limits of Heavy Metals and Metalloids for Wastewater Treatment Plants ( WWTPs ) in China Based on Statistical Methods

Industrial wastewater and sewage are both important sources of heavy metals and metalloids in urban wastewater treatment plants (WWTPs). China has made great efforts to control heavy metal and metalloid pollution by setting discharge limits for WWTPs. There is, however, limited discharge data and no systematic methodology for the derivation of discharge limits. In this study, 14 heavy metals and metalloids (Hg, alkyl mercury, As, Cd, Cr, Cr6+, Pb, Ni, Be, Ag, Cu, Zn, Mn, Se) that are listed in the Discharge Standard of Pollutants for Municipal Wastewater Treatment Plant (GB 18918-2002) were selected for the analysis of discharge characteristics while using the supervised monitoring data from more than 800 WWTPs located in nine provinces in China. Of the 14 heavy metals and metalloids, all but alkyl mercury were detected in the discharge water. There was a high rate of detection of As, Cu, Mn, Se, and there were some samples that exceeded the standard concentrations of Cr, Cr6+, Pb, and Ni. Removal rates of Hg, As, Cd, Cr, Cr6+, Pb, Ni, Cu, Zn, Mn, and Se were higher than 40%, comparable to values from other countries. Hg and As were selected to analyze the influencing factors of effluent and derive discharge limits of WWTPs using a statistical method, because these two metals had more detected data than other metals. The study used supervised monitoring data from Zhejiang WWTPs with 99 for Hg and 112 for As. Based on the delta-lognormal distribution, the results showed that geographic location was significantly closely correlated with Hg (P = 0.027 < 0.05) and As (P ≈ 0 < 0.05) discharge concentrations, while size (for Hg P = 0.695 > 0.05, for As P = 0.088 > 0.05) and influent concentration (R2 < 0.5) were not. Derived Hg and As discharge limits suggest that it is necessary to establish stricter discharge limits for WWTPs, which is more consistent with the real-world situation in China. The study here comprehensively researches the discharge characteristics of heavy metals and metalloids in effluent of WWTPs in China, and developed for the first time in China heavy metals and metalloids discharge limits based on statistical methods. The results may inform special discharge limit settings for WWTPs in China.


Introduction
Industrial wastewater and sewage are major sources of heavy metals and metalloids at urban wastewater treatment plants (WWTPs) Waste Water and Sewage Sludge [1] indicated that washing, drug use, and daily care were the main sources of heavy metal emissions.Oliver et al. [2] studied the heavy metals concentration in influent of WWTPs in Canada, which showed that the plants that received more industrial wastewater had higher concentrations of heavy metals than the plants that received only sewage, and Li et al. [3] had the similar conclusion when studied the influent of WWTPs received more industrial wastewater in China.Strengthening the monitoring and control of heavy metals and metalloids in WWTPs can play a substantial role in improving environmental quality and reducing environmental risks.Like those of most countries, China's discharge standards have been designed to regulate the end-of-pipe wastewater discharges.The national standard "Integrated Wastewater Discharge Standard (GB 8978-1996)" [4] first established separate discharge limits for WWTPs.With development of environmental standard systems in the country, in 2002, the national discharge standard of WWTPs (GB 18918-2002) [5] was established, and it has played an important role in controlling water pollutants discharged from WWTPs.Because some of the WWTPs in China process industrial wastewater, GB 18918-2002 has set the limits for 14 heavy metals and metalloids to prevent environmental risks from both industrial wastewater and sewage, including total mercury (Hg), alkyl mercury, total arsenic (As), total cadmium (Cd), total chromium (Cr), hexavalent chromium (Cr 6+ ), total nickel (Ni), total lead (Pb), total silver (Ag), total copper (Cu), total zinc (Zn), total beryllium (Be), total manganese (Mn), and total selenium (Se).However, the setting of limits for heavy metals and metalloids was based primarily on comparison to the limits for industrial wastewater control [6], without enough discharge data to analyze the difference to industrial wastewater.The limits setting also took some account of public health using times of surface water quality standards as the discharge limits considering the dilution [6] with a shortage of measured data to evaluate the real-world discharge level.
In recent years, China has made great efforts to collect discharge data, and these data, especially from WWTPs, have been used to reveal statistical distributions for pollutants discharge concentrations and evaluate the rationality of discharge standards [7][8][9].However, these studies mainly concerned the conventional pollutants as chemical oxygen demand (COD) and ammonia nitrogen rather than heavy metals and metalloids or other toxic pollutants.Furthermore, there has been a lack of an overall perspective regarding heavy metals and metalloids discharge situations of WWTPs and a systematic method for deriving the discharge limits based on statistical approaches in China.
This study aims to comprehensively research the discharge characteristics of heavy metals and metalloids in effluent of WWTPs in China, and for the first time to derive discharge limits of heavy metals and metalloids based on measured data analysis.In the present study, 14 heavy metals and metalloids controlled in GB 18918-2002 were selected for analysis of the discharge characteristics, and the Hg and As with more detected data were selected to derive discharge limits of WWTPs using a statistical method.We also addressed factors that may affect Hg and As discharge concentrations from different WWTPs.The result of this study may provide guidance for revisions to standards for heavy metals and metalloids discharge in China.

Data Source
In our study, supervised monitoring data carried out by local environmental protection authorities from Oct. 2016 to Mar. 2017 of concentrations of 14 heavy metals and metalloids discharged from more than 800 WWTPs located in nine provinces were used.These nine provinces were Jilin, Shanxi, Shandong, Zhejiang, Hubei, Hunan, Guangxi, Gansu, and Sichuan, which covered the northeast, northwest, north, east, south, and the southwest of China (Figure 1), including the higher heavy metal pollution level regions as Hunan, Guangxi, and Gansu [10].The WWTPs from data sources of different sizes covered various ratios of industrial wastewater treated, classified as <50%, 50-70% (including 50%) and ≥70%.In this way, the selected WWTPs were representative of Chinese WWTPs for locations, sizes, and ratios of the received industrial wastewater.data sources of different sizes covered various ratios of industrial wastewater treated, classified as <50%, 50-70% (including 50%) and ≥70%.In this way, the selected WWTPs were representative of Chinese WWTPs for locations, sizes, and ratios of the received industrial wastewater.Hg and As were selected to derive the discharge limits using statistical method because of higher ratios of detected data.When using the statistical method, there should be minimum data requirements of at least seven daily averages, and among these averages there should be at least three daily averages higher than the method detection limits [11].Based on this principle, for the sake of more valid data, we adopted the supervised monitoring data from 2015 to 2017 of Hg and As discharged by WWTPs in Zhejiang Province as the data source, which were also representative of Chinese WWTPs.There were 99 WWTPs for Hg and 112 WWTPs for As conforming to the data principles.

Data Analysis
The data collected were daily average concentrations of 14 heavy metals and metalloids sampled once every 2 h for 24 h and mixed for analysis, while using the Environment Monitoring Analytical Method Standards in China (Table 1).We obtained the averages, maximums, minimums, mediums for each pollutant to study the discharge characteristics as rate of detection, rate of exceeding standards, and similar outcomes.In order to investigate the statistical distribution of heavy metals and metalloids discharged by WWTPs in China, the Hg and As with more available data and higher detected rate were selected for analysis while using Origin Pro 7.0 (OriginLab Corporation, Northampton, MA, USA).Based on the fitted statistical distribution, we used One-way analysis of variance to study the factors influencing discharge level by PASW Statistics 18.0 (IBM, New York, NY, USA), including locations, sizes, and influent concentrations level of WWTPs.Hg and As were selected to derive the discharge limits using statistical method because of higher ratios of detected data.When using the statistical method, there should be minimum data requirements of at least seven daily averages, and among these averages there should be at least three daily averages higher than the method detection limits [11].Based on this principle, for the sake of more valid data, we adopted the supervised monitoring data from 2015 to 2017 of Hg and As discharged by WWTPs in Zhejiang Province as the data source, which were also representative of Chinese WWTPs.There were 99 WWTPs for Hg and 112 WWTPs for As conforming to the data principles.

Data Analysis
The data collected were daily average concentrations of 14 heavy metals and metalloids sampled once every 2 h for 24 h and mixed for analysis, while using the Environment Monitoring Analytical Method Standards in China (Table 1).We obtained the averages, maximums, minimums, mediums for each pollutant to study the discharge characteristics as rate of detection, rate of exceeding standards, and similar outcomes.In order to investigate the statistical distribution of heavy metals and metalloids discharged by WWTPs in China, the Hg and As with more available data and higher detected rate were selected for analysis while using Origin Pro 7.0 (OriginLab Corporation, Northampton, MA, USA).Based on the fitted statistical distribution, we used One-way analysis of variance to study the factors influencing discharge level by PASW Statistics 18.0 (IBM, New York, NY, USA), including locations, sizes, and influent concentrations level of WWTPs.Based on the fitted statistical distribution of heavy metals and metalloids discharge concentrations, discharge limits were derived by Equation (1) while using the long-term average (LTA) multiplied by the variability factor (VF).This equation has been used by the U.S. EPA in developing limits for many industries, such as the Organic Chemicals, Plastics and Synthetic Fibers (OCPSF) industry [25].In this equation, the average of concentrations and the fluctuations in the treatment system are both considered, which reflects a state that most treatment systems are capable of achieving.
LTA is the target value that a plant's treatment system should achieve on an average basis.In the present study, we used the average concentrations to derive the Hg and As LTAs for each WWTP when considering both the undetected and detected data [11,25].VF is the ratio of strong effluent to the average level, which expresses the relationship between large values and average treatment performance levels that a well-designed and operated treatment system should be capable of achieving all the time [11].The U.S. EPA used the 99th and 95th percentiles to express the daily maximum and monthly average VFs, respectively [11,25].Because the concept of monthly average was not used in China, we adopted the daily maximum 99th percentile VF in our study.

The Heavy Metal and Metalloid Discharge Concentrations Level of WWTPs
As shown (Table 2), there were 13 pollutants that were detected in the discharged water of WWTPs, except alkyl mercury.In relation to the detection rate, four pollutants were higher than 40%, including As, Cu, Mn, and Se, which indicated the higher exposure of these four pollutants in WWTPs of China.The detection rate of Cr, Hg and Zn were higher than 20%, which indicated the relative high exposure of the three pollutants also.In the point view of discharge concentrations, among the 13 detected pollutants, there were cases of over the standard for four pollutants, including Cr, Cr 6+ , Pb, and Ni, with the standard exceeding ratio being 0.9%, 0.4%, 0.1%, and 17.0%, respectively, which most occurred in the WWTPs received wastewater from electroplate plants.The ratios of average concentrations to discharge limits in GB 18918-2002 for other nine pollutants were between 0.02 and 0.4, which were 1-2 orders of magnitude lower than the limits and indicated the low level of discharge concentrations.The ratios of average concentrations to mediums were all ≥1.0 for 13 pollutants, which indicated that most of the concentrations were lower than the average.When comparing the discharge concentrations of WWTPs in China to other countries, they were equivalent (Table 3).

Removal Rate of Heavy Metals and Metalloids in WWTPs in China
When researching the removal rate of heavy metals and metalloids in WWTPs, the discharge concentrations of undetected data were calculated with the Equation ( 2).The removal rates of 14 heavy metals and metalloids are shown in Table 4.

Discharge concentration of undetected data
where: MDL-Method Detection Limit.As noted in Table 4, the removal rates of Hg, Cd, Cr, Cr 6+ , Pb, Ni, Cu, Zn, and Se in the WWTPs are higher than 50% ranging from 50.8% to 69.8%, while the removal rate of As and Mn are 47.3% and 42.8%, respectively.The rate of Be is lower than 30%.
Removal rates of heavy metals and metalloids in China's WWTPs were comparable to those reported in research findings in other countries.Obarska-Pempkowiak [32] showed that the removal rate of zinc in WWTPs was 50%, copper was 60%, lead was 79%, and cadmium was 80%.Olive et al. [2] found that the removal rate of aluminum, cadmium, chromium, copper, iron, lead, mercury, and zinc were all higher than 70% in Canadian WWTPs, but they reported lower values for manganese, nickel and strontium.Kulbat et al. [29] found that the rate of zinc removal was 78.1%-86.2%,copper was 78.1%-93.2%,chromium was 45.3%-66.7%,silver was 28.3%-66.7%,lead was 48%-83% in Polish WWTPs, and these research results indicated higher removal rate of heavy metals.
It is difficult to discuss the influence of treatment process to the removal rate of heavy metals and metalloids under any but the most controlled conditions.Most metals in wastewater are soluble, but if mixed with sewage, the changes in pH and the act of mixing with other organics can render many types of metals insoluble and convert them to sludge, leading to high removal rates in wastewater [2].The removal rate is affected by many factors, including type of metal, concentration in the influent, interactions with microbes in the sewage treatment system, and the treatment processes that are employed by the WWTPs [29].

Statistical Distribution of Heavy Metal and Metalloid Discharge Concentrations
Among the 1080 Hg discharge figures that were collected from 872 WWTPs, 792 indicated that no mercury was detected and 288 indicated mercury had been detected.Among the 925 As discharge data from 858 WWTPs, 387 were not detected, while 537 were detected.While using the detected daily average data, as shown in Figures 2 and 3, it the lognormal distribution was visible for both Hg and As detected according to the fitting lines with the logarithm of Hg or As concentrations on the horizontal x-axis against the normal probabilities on the vertical y-axis.This distribution was also verified while using the P-P graph test (Figures 4 and 5) with the fitting lines of cumulative probability of monitoring data after logarithmic transformation against the expected cumulative probability.The lognormal distribution of effluent concentrations is also consistent with other research [7,33].When considering some of the undetected data, the delta-lognormal distribution provides a reasonable and practical basis for further analysis of the heavy metal and metalloid data and determining the discharge limits, which is in accordance with the EPA regulations [11,25,34].The delta-lognormal distribution consists of two parts, undetected numbers and detected numbers conforming to the lognormal distribution.The distribution function of delta-lognormal distribution can be expressed as [11,25,34]: where: δ -the ratio of non-detected, -the function of lognormal distribution.

Factors Influencing Heavy Metal and Metalloid Discharge Concentrations
There are some data requirements for One-way analysis of variance, conforming to normal distribution, independence, and homogeneity variance.The logarithm of detected data have been shown to fit normal distribution, and because the data were all from independent samples, they met the independence requirement.The homogeneity of variance test was taken with the One-way When considering some of the undetected data, the delta-lognormal distribution provides a reasonable and practical basis for further analysis of the heavy metal and metalloid data and determining the discharge limits, which is in accordance with the EPA regulations [11,25,34].The delta-lognormal distribution consists of two parts, undetected numbers and detected numbers conforming to the lognormal distribution.The distribution function of delta-lognormal distribution can be expressed as [11,25,34]: where: δ-the ratio of non-detected, 0 ≤ δ ≤ 1; when x 0 = ND, I(x 0 )= 1, or else I(x 0 )= 0; g(x)-the function of lognormal distribution.

Factors Influencing Heavy Metal and Metalloid Discharge Concentrations
There are some data requirements for One-way analysis of variance, conforming to normal distribution, independence, and homogeneity variance.The logarithm of detected data have been shown to fit normal distribution, and because the data were all from independent samples, they met the independence requirement.The homogeneity of variance test was taken with the One-way analysis of variance, if any sample did not pass the test, the Kruskal-Wallis test was used for analysis.
3.4.1.Regional Differences Discharge concentrations of heavy metals and metalloids in WWTPs from nine provinces were analyzed.Figure 6 indicates that the Hg concentrations from Shandong, Zhejiang, Guangxi, Gansu, and Shanxi WWTPs were relatively high with about 50% WWTPs in these regions being higher than 0.0001 mg/L.Figure 7 indicates that the As concentrations from Zhejiang, Hunan, Guangxi, Gansu, and Shanxi WWTPs were relatively high with average concentrations of As in these regions higher than 0.002 mg/L.According to China's "Environment Statistical Annual Report 2015" [10] of China, Guangxi, and Gansu are top two regions in which more Hg was discharged in wastewater, because of large number of facilities that engage in nonferrous metallurgy, chemical raw materials, and chemical products producing.This was consistent with the results of data analysis.
The effect of geographical location on Hg and As discharge concentration levels in WWTPs was analyzed while using One-way analysis of variance with the independent variable being the logarithm of concentration data, and the variables being different regions.Because the data did not pass the homogeneity of variance test, Kruskal-Wallis test was used and the results showed there to be significant differences in Hg (P = 0.027 < 0.05) and As (P ≈ 0 < 0.05) discharge concentrations for different regions, with the Hunan Province showing significant differences from all other regions.

Scale Differences
The numbers of samples from WWTPs in nine provinces that consisted of the scale information were 95 for Hg and 187 for As.Among these data, the average concentrations of Hg and As for largescale (treatment capacity ≥ 0.1 million m 3 /day), mid-scale (10 thousand m 3

Scale Differences
The numbers of samples from WWTPs in nine provinces that consisted of the scale information were 95 for Hg and 187 for As.Among these data, the average concentrations of Hg and As for largescale (treatment capacity ≥ 0.1 million m 3 /day), mid-scale (10 thousand m 3 /day ≤ treatment capacity Water 2018, 10, 1248 9 of 17

Scale Differences
The numbers of samples from WWTPs in nine provinces that consisted of the scale information were 95 for Hg and 187 for As.Among these data, the average concentrations of Hg and As for large-scale (treatment capacity ≥ 0.1 million m 3 /day), mid-scale (10 thousand m 3 /day ≤ treatment capacity < 0.1 million m 3 /day), and small-scale (treatment capacity < 10 thousand m 3 /day) of WWTPs are shown in Table 5.We can find the Hg concentrations from different scale WWTPs are close, but As concentrations from large scale are higher.The latter may be attributable to 6 of the 18 large-scale WWTPs with high influent concentrations.The scale factor is not significant to Hg (P = 0.695 > 0.05) according to the One-way analysis of variance and not significant for As (P = 0.088 > 0.05) either, according to Kruskal-Wallis test.

Relationship between the Concentrations of Influent and Effluent
The relationship between the concentrations of influent and effluent was analyzed (Figures 8 and 9).The results show that the linear relationship between the concentrations of influent and effluent is weak (R 2 < 0.5).Gbondo-Tugbawa et al. [26] also reported that there is no significant relationship between the mercury influent concentration and effluent.
Water 2018, 10, x FOR PEER REVIEW 10 of 18

Relationship between the Concentrations of Influent and Effluent
The relationship between the concentrations of influent and effluent was analyzed (Figures 8  and 9).The results show that the linear relationship between the concentrations of influent and effluent is weak (R 2 < 0.5).Gbondo-Tugbawa et al. [26] also reported that there is no significant relationship between the mercury influent concentration and effluent.Research has been performed on influent and effluent concentration data of Hg and As from WWTPs in Zhejiang Province.From the Table 6 we can conclude that the standard deviation of Hg and As influent concentration is higher than the effluent concentration, which suggests that the biochemical treatment system could play a role in buffering the fluctuation of heavy metals and metalloids, such as Hg and As, and there was the same conclusion in Balogh et al. [35] and Li et al.The relationship between the concentrations of influent and effluent was analyzed (Figures 8  and 9).The results show the linear relationship between the concentrations of influent and effluent is weak (R 2 < 0.5).Gbondo-Tugbawa et al. [26] also reported that there is no significant relationship between the mercury influent concentration and effluent.Research has been performed on influent and effluent concentration data of Hg and As from WWTPs in Zhejiang Province.From the Table 6 we can conclude that the standard deviation of Hg and As influent concentration is higher than the effluent concentration, which suggests that the biochemical treatment system could play a role in buffering the fluctuation of heavy metals and metalloids, such as Hg and As, and there was the same conclusion in Balogh et al. [35] and Li et al. [36] research, which could explain the weak linear relationship between influent and effluent Research has been performed on influent and effluent concentration data of Hg and As from WWTPs in Zhejiang Province.From the Table 6 we can conclude that the standard deviation of Hg and As influent concentration is higher than the effluent concentration, which suggests that the biochemical treatment system could play a role in buffering the fluctuation of heavy metals and metalloids, such as Hg and As, and there was the same conclusion in Balogh et al. [35] and Li et al. [36] research, which could explain the weak linear relationship between influent and effluent concentrations.

Long-Time Averages
Based on the delta-lognormal distribution, the long-time averages were calculated, as follows [11,25,34]: where: LTA j -the long-time average of heavy metal or metalloid of WWTP j ; δ-the rate of below method detection limit; D-method detection limit; X i -the concentration higher than method detection limit; n-the number of concentration data up the method detection limit.
According to the equation above, the long-time averages of Hg for 99 WWTPs and As for 112 WWTPs are shown in Tables 7 and 8.The LTAs for Hg are from 0.00004 to 0.00081 mg/L with an average of 0.00015 mg/L and medium of 0.00009 mg/L.The LTAs for As are from 0.0003 to 0.0250 mg/L with an average of 0.0021 mg/L and medium of 0.0012 mg/L.
We can conclude that the LTAs of Hg and As are close to the level of surface water quality limits Grade I~III in "Surface water environmental quality standards" (GB 3838-2002) [37], and also comparable to the results of other works.For example, research shows the influent concentration of Hg was about 0.001 mg/L, varying with time and location [38].The average influent concentration of Hg of 3 WWTPs in Canada was 0.000061 mg/L, and the discharge concentrations were from 0.000003-0.000014mg/L [39].Some studies in China also showed that the Hg concentrations from sewage were 0.0002-0.0024mg/L, and the discharge concentrations were 0.00001-0.000029mg/L, which was much lower than the limits in GB 18918-2002 [5].Lu et al. [40] studied on the As concentration of WWTPs in Chongqing Province, showing that the influent concentration was about 0.00105 mg/L, and the discharge concentration was about 0.0007 mg/L.Based on the delta-lognormal distribution of Hg and As discharge concentrations, the daily maximum variability factors were calculated while using the following equations: where: When δ ≥ 0.99, P99 = D; When δ < 0.99, P99 = max(D, e µ+zσ ) and, where: VF(1)-daily maximum variable factor; P99 -estimated 99th percentile of Hg or As discharge concentration at a given plant; Ê(x)-estimated expected value of Hg or As discharge concentration at a given plant; y i -natural logarithm of Hg or As discharge concentration other than ND; φ-normal cumulated distribution function value.
Based on the method described above, the Hg daily maximum variable factors for 99 WWTPs and As daily maximum variable factors for 112 WWTPs are shown in Tables 9 and 10.As shown in the tables, it can find that the Hg daily maximum variable factors are 1.34-11.07with an average of 3.89 and medium of 3.33.The As daily maximum variable factors are 1.44-16.05with an average of 3.77 and medium of 3.35.We obtained the Hg and As discharge limits for each WWTP (Tables 9 and 10) while using Equation (1).From these results, it can find that the Hg discharge limits are 0.00007-0.00902mg/L with an average of 0.00067 mg/L and medium of 0.00035 mg/L.These values are lower than the limit (0.001 mg/L) in GB 18918-2002.Derived discharge limits for As are 0.0005-0.1482mg/L with an average of 0.0092 mg/L and medium of 0.0041 mg/L, which are also lower than the limit (0.1 mg/L) in GB 18918-2002.When considering that the national discharge limits should be suitable for the entire country, we used the medium value of LTAs and average value of VF(1)s as the suggested Hg and As discharge limits [11] (Table 11).

Conclusions
This study analyzed the discharge characteristics of 14 heavy metals and metalloids at WWTPs in China.The study used Hg and As as an example to derive the discharge limits based on statistical methods.In this study, 13 of the 14 heavy metals and metalloids, with the exception of alkyl mercury, were detected in the discharge water with high rates of detection of As, Cu, Mn, and Se.In many cases, Cr, Cr 6+ , Pb, and Ni exceeded standards.Removal rates of Hg, As, Cd, Cr, Cr 6+ , Pb, Ni, Cu, Zn, Mn, and Se exceeded 40%, which is consistent with values in other countries.
Results from this study suggest that the delta-lognormal distribution was suitable for the heavy metal and metalloid discharge concentrations of WWTPs.Based on the statistical distribution, the study found that geographic location was strongly correlated with Hg and As discharge concentrations, while size and influent concentration were not.The derived LTAs for Hg of 99 WWTPs and As of 112 WWTPs showed that they were much lower than the limits in GB 18918-2002.The derived VF(1)s of WWTPs revealed that the Hg and As fluctuations were not large, suggesting that WWTPs in China are relatively well-managed.Lastly, the derived Hg and As discharge limits suggest that it is necessary to establish 1-2 orders of magnitude lower discharge limits for WWTPs in China.

Figure 1 .
Figure 1.Location of nine investigated provinces of China.

Figure 1 .
Figure 1.Location of nine investigated provinces of China.
Water 2018, 10, x FOR PEER REVIEW 7 of 18

Figure 2 .
Figure 2. Statistical distribution of detected daily average concentrations of Hg of the WWTPs.

Figure 2 .
Figure 2. Statistical distribution of detected daily average concentrations of Hg of the WWTPs.

Figure 2 .
Statistical distribution of detected daily average concentrations of Hg of the WWTPs.

Figure 3 .
Figure 3. Statistical distribution of detected daily average concentrations of As of the WWTPs.

Figure 4 .
Figure 4. P-P graph of logarithm of the detected daily average concentrations of Hg.

Figure 3 .
Figure 3. Statistical distribution of detected daily average concentrations of As of the WWTPs.

Figure 2 .
Figure 2. Statistical distribution of detected daily average concentrations of Hg of the WWTPs.

Figure 3 .
Figure 3. Statistical distribution of detected daily average concentrations of As of the WWTPs.

Figure 4 .
Figure 4. P-P graph of logarithm of the detected daily average concentrations of Hg.

Figure 4 . 18 Figure 5 .
Figure 4. P-P graph of logarithm of the detected daily average concentrations of Hg.Water 2018, 10, x FOR PEER REVIEW 8 of 18

Figure 5 .
Figure 5. P-P graph of logarithm of the detected daily average concentrations of As.

Water 2018 ,
10, x FOR PEER REVIEW 9 of 18

Figure 6 .
Figure 6.Boxplot of discharge concentrations of Hg from different regions WWTPs.Note: Because there was not enough data from Jilin WWTPs, so there was no boxplot for Jilin.

Figure 7 .
Figure 7. Boxplot of discharge concentrations of As from different regions WWTPs.

Figure 6 .
Figure 6.Boxplot of discharge concentrations of Hg from different regions WWTPs.Note: Because there was not enough data from Jilin WWTPs, so there was no boxplot for Jilin.

Figure 7 .
Figure 7. Boxplot of discharge concentrations of As from different regions WWTPs.

Figure 7 .
Figure 7. Boxplot of discharge concentrations of As from different regions WWTPs." " stands for average value, and "×" stands for 95% percentiles.

Figure 8 .
Figure 8. Relationship between the concentration of influent and effluent for Hg from WWTPs.

Figure 9 .
Figure 9. Relationship between the concentration of influent and effluent for As from WWTPs.

Figure 8 .
Figure 8. Relationship between the concentration of influent and effluent for Hg from WWTPs.

Figure 8 .
Figure 8. Relationship between the concentration of influent and effluent for Hg from WWTPs.

Figure 9 .
Figure 9. Relationship between the concentration of influent and effluent for As from WWTPs.

Figure 9 .
Figure 9. Relationship between the concentration of influent and effluent for As from WWTPs.

Table 1 .
Environment monitoring analytical methods used for analysis of 14 heavy metals and metalloids.

Table 2 .
Discharge concentrations of 14 heavy metals and metalloids in wastewater treatment plants (WWTPs) of China.

Table 3 .
Comparison of discharge concentrations of 14 heavy metals and metalloids in WWTPs of China and other countries.

Table 4 .
Removal rates of heavy metals in WWTPs of China.
1 "-" indicates that the data is not enough to get the results.

ShandongZhejiang Sichuan Hubei Hunan Guangxi Gansu Shanxi 0.0000 0.0001 0.0002 0.0003
Boxplot of discharge concentrations of Hg from different regions WWTPs.Note: Because there was not enough data from Jilin WWTPs, so there was no boxplot for Jilin." " stands for average value, and "×" stands for 95% percentiles.
/day ≤ treatment capacity Water 2018, 10, x FOR PEER REVIEW 9 of 18

Table 5 .
Discharge concentrations of Hg and As for different scales of WWTPs.

Table 6 .
Average of standard deviations of influent and effluent for Hg and As.

Table 9 .
VF(1)s and derived discharge limits of Hg for 99 WWTPs.

Table 10 .
VF(1)s and derived discharge limits of As for 112 WWTPs.

Table 11 .
Suggested Hg and As discharge limits of WWTPs in China.