Characteristics of PCDD/Fs in Flue Gas from MSWIs and HWIs: Emission Levels, Profiles and Environmental Influence

Received: December 27, 2019
Revised: March 25, 2020
Accepted: March 29, 2020

 Copyright The Author(s). This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are cited.

Cite this article:

Qiu, J., Tang, M., Peng, Y., Lu, S., Li, X. and Yan, J. (2020). Characteristics of PCDD/Fs in Flue Gas from MSWIs and HWIs: Emission Levels, Profiles and Environmental Influence. Aerosol Air Qual. Res. 20: 2085–2097. https://doi.org/10.4209/aaqr.2019.11.0610

HIGHLIGHTS

• Characteristics of PCDD/Fs in flue gas from MSWIs and HWIs was studied.
• Emission factors and annual emission amounts in the stack gas were estimated.
• PCDD/F levels and profiles in ambient air and soil around incinerators were discussed.

Keywords: PCDD/Fs; Flue gas; Emission factor; Air; Soil.

INTRODUCTION

Municipal solid waste (MSW) incineration in China has been boosted more than twelve fold in the past decade, in response to the rapid increase of MSW (Lu et al., 2017; Zhu et al., 2018). In addition, the amount of hazardous waste (HW) in China is increasing greatly from 10.79 million tons in 2007 to 69.37 million tons in 2017 as the development of economy (National Bureau of Statistics of China, 2008, 2018). Hazardous waste, including industrial hazardous waste, medical waste, and household hazardous waste (Duan et al., 2008), presents some particular characteristics of ignitability, corrosivity, reactivity, and toxicity. Trinh and Chang (2018) reported that municipal solid waste incinerators (MSWIs) were important emission sources of polychlorinated dibenzo-p-dioxins and dibenzofurans (PCDD/Fs, also referred to as dioxin) to ambient air (Ni et al., 2009; Hu et al., 2014; Wang et al., 2015; Bai et al., 2017; Ngo et al., 2018; Trinh and Chang, 2018). Dopico and Gomez (2015) and Zhan et al. (2019) also reported that the emission factor of PCDD/Fs from hazardous waste incinerators (HWIs) is significantly higher than those from MSWIs. Therefore, it is important to control the PCDD/F emissions from HWIs.

The emission of PCDD/Fs has recently been investigated extensively in the waste incineration process (Bie et al., 2007; Choi et al., 2008; Gao et al., 2009; Ni et al., 2009; Zhu et al., 2018). Choi et al. (2008) investigated dioxins in gases emitted from MSWIs and HWIs, finding the concentration of PCDD/Fs ranged from 0.05 to 609.27 ng TEQ Nm^–3. Bie et al. (2007) discovered the average dioxin emission in the flue gases from MSWIs by applying a twin internal fluidized bed incinerator in Harbin was 0.02 ng I-TEQ Nm^–3. Ni et al. (2009) analyzed the congener profiles from various MSWIs, estimating the emission factors of PCDD/Fs in the stack gas to be 0.169–10.7 µg TEQ t^–1 MSW with an average value of 1.73 µg TEQ t^–1MSW. Gao et al. (2009) revealed the stack gas emissions of PCDD/Fs from hospital waste incinerators (HWIs) exhibited a large variation (0.08–31.6 ng I-TEQ Nm^–3). Zhu et al. (2018) revealed that the concentrations of PCDD/F emissions in stack gas from 6 MSWIs were in the range of 0.007–0.059 ng TEQ Nm^–3. The emission factors of PCDD/Fs from 6 MSWIs varied from 0.027 to 0.225 mg I-TEQ t^–1, with a mean value of 0.17 mg I-TEQ t^–1 waste.

The PCDD/F concentrations and congeners distribution in the surrounding environment of different incinerators were investigated and studied in the past years (Oh et al., 2006; Yan et al., 2008; Li et al., 2018; Zhou et al., 2019). Oh et al. (2006) studied the influence of MSWI on the PCDD/Fs levels in ambient air and soil, and the concentration of PCDD/Fs in air samples and soil samples ranged from 0.22 to 1.16 ng I-TEQ Nm^–3 and 1.25 to 75.0 ng I-TEQ kg^–1, respectively. Zhou et al. (2019) studied the characteristics of PCDD/Fs in surface soil around an iron and steel plant and the PCDD/Fs content in the soil ranged from 0.16 to 4.5 ng I-TEQ kg^–1. However, few studies focused on the emission characteristics of PCDD/Fs and PCDD/Fs levels in the ambient air and soil samples from HWIs, and the similarities and differences between MSWIs and HWIs were still obscure in terms of PCDD/Fs.

In order to investigate the emission characteristics of PCDD/Fs the in stack gas and PCDD/Fs levels in the ambient air, soil samples from MSWIs and HWIs in China, stack gas samples from 5 MSWIs and 4 HWIs and ambient air and soil samples around the plants were collected and analyzed. The objectives of this study were: (1) to compare the PCDD/F concentration of stack gas and congener profiles of PCDD/Fs from 5 MSWIs and 4 HWIs; (2) to investigate the emission factors of PCDD/Fs and the annual emission of PCDD/Fs; (3) to analyze the possible grouping of similar emissions and demonstrate the PCDD/F congeners profiles using principal component analysis (PCA); (4) to study the PCDD/F concentration and patterns in the ambient air and soil samples from MSWIs and HWIs.

MATERIALS AND METHODS 

Sample Collection

Stack gas samples were collected from municipal solid waste incinerators (MW01–MW05) and hazardous waste incinerators (HW01–HW04) in China, and the basic information for these incinerators and their APCDs were given in Table 1. The MW01–MW05 treated municipal solid waste incinerators. The HW01 treated dye waste, coke, organic sludge, waste oil, paint waste water with rotary kiln. The HW02 dealt with waste gas and liquid produced during VCM production. The HW03 treated wastes composed of oily refuse, rock debris, mud, etc. The HW04 treated medical waste. 

Flue gas sampling was conducted using isokinetic samplers (Model KNJ23, KNJ corporation, Korea), complying with U.S. EPA Method 23a; and the specific sampling method was described in detail in previous study (Chen et al., 2008). All sampling points were located downstream to APCDs. Ambient air sampling was performed in winter and summer at the boundary of the plants over 2 meters in height, and down the wind. Ambient air samples were collected according to U.S. EPA Method TO-9A (U.S. EPA, 1999), as described in detail by Chen et al. (2015). The sampling time and air volumes were about 24 h and 900–1100 m³, respectively (Table 2). The soil samples (approximately 2 kg per sample) were taken from the upper 200 mm (MW01–MW05, HW02, HW04) and 200-400 mm (HW01, HW03) at the same sites. 

Sample Extraction and Analysis

The PCDD/Fs analysis of stack gas and ambient air samples were performed according to the U.S. EPA Method 23 (U.S. EPA, 1995). Detailed steps of analysis method of the PCDD/Fs could be found from our previous studies (Chen et al., 2014; Chen et al., 2016). The soil samples were dried until constant weight and then they were ground and screened by a 2-mm sieve. About 500 g of each soil sample was finally homogenized through a 60-mesh sieve, and was refrigerated until analysis. And the detailed steps of pretreatment and analysis could be found from our previous study (Yan et al., 2008).

Statistical Analysis

To better understand the congener profiles of dioxins in stack gases from various MSWIs and HWIs, principal component analysis (PCA) were used to evaluate the possible groupings of similar emissions and the dominant congeners in a defined grouping. Prior to analysis, the PCDD/F congeners levels were normalized with respect to the total PCDD/F concentration in the relevant sample. The statistical analyses were performed using the SPSS 20.0 software.

RESULTS AND DISCUSSION

TEQ Contribution of PCDD/Fs in Stack Gas

The gas emission levels of PCDD/Fs and the PCDF/PCDD weight ratio from 9 waste incineration facilities in China were investigated, and the results were presented in Fig. 1. The emissions of PCDD/Fs ranged from 0.016 to 0.29 ng I-TEQ m^–3 with an average value of 0.103 ng I-TEQ m^–3. These values except HW02 and HW03 were all below the emission limit of national standard in China (0.1 ng I-TEQ m^–3) (China MEP, 2014). However, the PCDD/F emissions from HW02 and HW03 were 0.29 and 0.25 ng I-TEQ m^–3 respectively, indicating more improvement space for PCDD/Fs control technologies in hazardous waste incinerators. The proportion of PCDF much higher than that of PCDD as in Fig. 1 indicating that PCDF was the main contributor to TEQ value in flue gas. Pham et al. (2019a) reported that concentrations of PCDD/Fs in flue gas were in the range of 13.4–17.6 pg TEQ m^–3 (mean: 15.5 pg TEQ m^–3) for MSWIs and 12.7–438 pg TEQ m^–3 (mean: 125 pg TEQ m^–3) for medical waste incinerators (MWIs), respectively. Zhu et al. (2018) reported that concentrations of PCDD/Fs in stack gas from MSWIs were in the range of 0.007-0.059 ng TEQ m^–3. 

Fig. 1. Emission levels of PCDD/Fs from the 9 waste incinerators.

As in Fig. 2, the 50% for HWIs was much higher than that for MSWIs, indicating that HWIs were generally characterized by higher emission levels than MSWIs in China.Three potential factors were identified to affect the concentration and the congener distributions of PCDD/Fs in incinerators: the types of furnace, the composition of APCDs, and waste composition (chlorine content in the waste) (Wang et al., 2003; Gao et al., 2009). Wastes incinerated at a low temperature yield a higher TEQ value than a high temperature. Besides, the absence or removal of chlorine in wastes effectively decreased the emission TEQ values (Wikstrom et al., 1996; Nzihou et al., 2012). Many studies focused on the relationship between chlorine and PCDD/Fs, and the effect of chlorine content in the feeding waste on PCDD/Fs formation was obviously. Wang et al. (2003) and Li et al. (2017) reported that MSW contains 0.2–0.8% chlorine, while the content of chlorine in hazardous waste is generally higher than municipal solid waste, which accelerates the formation of PCDD/F. The samples that did not contain chlorine or were not combusted with chlorides exhibited low emission TEQ values (Wikstrom et al., 1996; Nzihou et al., 2012). 

Fig. 2. PCDD/Fs statistical parameters for stack gas samples from MSWIs and HWIs: mean (▪), 25%, 50% and 75%.

Estimated Emission Factor and Emission of PCDD/Fs from MSWIs and HWIs

Much work so far focused on the gas emission of PCDD/Fs from MSWIs and HWIs in China by using the emission factor (EF) to calculate the annual emission (EA) of PCDD/Fs (Li et al., 2017; Zhu et al., 2018; Pham et al., 2019a, b).

EF (µg I-TEQ t^–1) = [concentration in flue gas (ng I-TEQ m^–3) × flue gas flow rate (m³·h^–1)]/[capacity level (t h^–1) × 10³] (1)

EA (mg I-TEQ y^–1) = [emission factor (µg I-TEQ t^–1) × capacity level (t h^–1) × average operating (h y^–1)]/10³        (2)

Estimated emission factors and annual emission amounts of PCDD/Fs (mean concentrations and the ranges) from flue gas of the 5 MSWIs, 4 HWIs and the cases from literatures were presented in Table 3. 

In this study, the emission factors of PCDD/Fs to the atmosphere from 5 MSWIs varied from 0.056 to 0.607 µg I-TEQ t^–1 MSW (mean: 0.235 µg I-TEQ t^–1). It was estimated that 10.2-88.0 mg I-TEQ of PCDD/Fs (mean: 35.0 mg I-TEQ of PCDD/Fs) were annually released to the atmosphere from MSWIs. In particular, the mean EF values and EA values of the MSW were similar, respectively. These values were in good consistent with the study of Zhu et al. (2018), in which the EF and EA values ranged in 0.027–0.225 µg I-TEQ t^–1 (mean value of 0.17 µg I-TEQ t^–1) and in 5.36–101 mg I-TEQ (mean: 35.0 mg I-TEQ), respectively.

The EF values ranged from 0.07 to 3.27 µg I-TEQ t^–1 for HWIs (mean: 1.03 µg I-TEQ t^–1), varying significantly for different HWIs. PCDD/Fs annual released to the atmosphere were estimated as 1.13–53.0 mg I-TEQ (mean: 17.1 mg I-TEQ) in HWI. The noticeable EF value of HW03 among HWIs could be attributed to the particularly treated wastes, which were composed of oily refuse, rock debris, mud, etc.

Until 2018, 352 MSWIs which had a yearly capacity of 89.4 million tons were running in China, then the estimated annual emission amounts of PCDD/Fs from flue gas was 22.3 g TEQ. And the scale of HWI treatment reached 6.06 million tons per year, the estimated annual emission amounts of PCDD/Fs from flue gas was 6.51 g TEQ.

Congener Profiles of PCDD/Fs in Stack Gas

Fig. 3 shows the scope plot of principal component analysis, for relative contributions of the 17 toxic dioxin congeners. The first principal component (PC1) took into account 48.0% of the variability of the data set, and the second principal component (PC2) accounted for 43.7% of the total variance. PC1 and PC2 accounted for 91.7% of the total variance. The score plot of PCA indicated that all investigated MSWIs and HWIs could be assigned to three groups. Group1 on the left were samples from HW01 and HW02, group 2 in the middle were samples from MSWIs and HW04, and group 3 on the right included samples from HW03. As in Fig. 3, the distribution of PCDD/F congeners was similar between the samples from MSWIs, but different between the samples from MSWIs and HWIs. Moreover, samples from HWIs were also distinct. 

Fig. 3. Score plot of principal component analysis of 17 toxic PCDD/F emission.

Fig. 4 presented the congener profiles of PCDD/Fs in each group. Congener profiles of 17 PCDD/F congeners in the samples were analyzed to understand the formation pathways and release behavior of dioxins in the investigated facilities. The mass ratios of PCDFs to PCDDs (R_F/D) was always used to estimate the formation of dioxins by the de novo synthesis (R_F/D > 1) or by precursor formation (R_F/D < 1) (Luijk et al., 1994; Huang and Buekens, 1995; Tuppurainen et al., 2000; Everaert and Baeyens, 2002; Pham et al., 2019b). The congener distributions were obtained by normalizing each congener to the total weight of all 2,3,7,8-substituted PCDD/F congeners. In Group 1, the proportion of PCDFs to the total amount of PCDD/Fs was 88.8% with a high proportion of 2,3,4,7,8-PeCDF (28.4%), 1,2,3,4,7,8-HxCDF (20.6%), 1,2,3,6,7,8-HxCDF (14.1%) and 2,3,4,6,7,8-HxCDF (10.3%). In Group 2, 2,3,4,7,8-PeCDF was the most abundant congener (35.3%) with higher fraction of 2,3,4,6,7,8-HxCDF (10.2%), and the total contribution of 2,3,7,8-substitued PCDFs was 70.8%. In Group 3, the congener 2,3,4,7,8-PeCDF contributed dominantly to the total TEQ with 31.4%. Other significant contributor were as follows: 2,3,7,8-TCDD (24.6%), 2,3,7,8-TCDF (17.5%) and 1,2,3,7,8-PeCDD (14.4%). PCDF isomer distributions indicated that low-chlorinated dibenzofurans were more than the high-chlorinated dibenzofurans. The ratios of PCDFs to PCDDs for all gas samples were over 1, implying the dominant role of de novo synthesis (Everaert and Baeyens, 2002; Ooi and Lu, 2011). 

Fig. 4. Representative congener patterns of PCDD/Fs in each group (Group 1: H01–H06, Group 2: M01–M39 and H10–H12, Group 3: H07–H09).

Concentration and TEQ Distribution of PCDD/F Congeners in Stack Gas from MSWIs and HWIs

Emission concentrations of 17 PCDD/Fs congeners and TEQ values in stack gas from 5 MSWIs and 4 HWIs were presented in Figs. 5 and 6. As shown in Fig. 5, the contribution of OCDD (37.6%) and 1,2,3,4,6,7,8-HpCDD (14.2%) in terms of total concentration were higher than other congeners. Notably, OCDD is the dominant congener in stack gas. However, for TEQ concentration, the 2,3,4,7,8-PeCDF, 2,3,4,6,7,8-HxCDF and 1,2,3,7,8-PeCDD were the dominant congeners of PCDD/Fs in stack gas, and accounted for 38.5%, 9.4% and 8.86%, respectively. They contribute over 50% of the total TEQ. Li et al. (2017) and Zhu et al. (2018) also found that 2,3,4,7,8-PeCDF was the largest contributor in MWIs, this conclusion was supported by the previous works (Bie et al., 2007; Choi et al., 2008; Gao et al., 2009; Yoon et al., 2017; Trinh and Chang, 2018; Zhu et al., 2018). 

Fig. 5. PCDD/F congeners distribution and TEQ in stack gas from MSWIs.

The congener profiles in the stack gas from HWIs were shown in Fig. 6. According to these TEQ results, the congener profiles of HWIs could be classified into three main groups: (1) PCDFs were more abundant than PCDDs with 2,3,4,7,8- PeCDF (28.4%), 1,2,3,4,7,8-HxCDF(20.6%) and 1,2,3,6,7,8-HxCDF(14.1%) as major congeners (e.g., HW01 and HW02); (2) the percentages of PCDFs and PCDDs were 59.17% and 40.83%, respectively, with 2,3,4,7,8-PeCDF (31.4%) and 2,3,7,8-TCDD (24.61%) as major congeners (e.g., HW03); (3) PCDFs were more abundant than PCDDs with high proportions of 2,3,4,7,8-PeCDF (37.34%) congener (e.g., HW04). We concluded that the de novo synthesis was always dominant. 

Fig. 6. PCDD/F congeners distribution and TEQ in stack gas from HWIs (a: HW01, b: HW02, c: HW03, d: HW04).

Influence of PCDD/F Emissions on the Surrounding Atmosphere

The concentrations and the PCDD/Fs patterns in the ambient air were studied to determine if the PCDD/F emissions from each plant affecting their surrounding environment.

As could be seen from Table 2, the PCDD/F concentrations in the ambient air varied between 0.088 and 13.9 pg I-TEQ m^–3 with an average value of 0.185 pg I-TEQ m^–3. Among the monitoring sites, MW01-A3 showed the highest PCDD/F concentration, which was roughly two order of magnitude higher than those in other sites. The possible factor was due to the location that was by the chimney. These values except MW01-A3 which was by the chimney were all below the emission limit of national standard in Japan (0.6 pg I-TEQ m^–3). Xu et al. (2009) founded that the PCDD/F concentrations in the ambient air varied from 0.059 to 3.03 pg I-TEQ m^–3 near the MSWIs in eastern China. Oh et al. (2006) reported that the levels of PCDD/Fs in ambient air near a MSWI in Korea were 0.22–1.16 pg I-TEQ m^–3 with an average of 0.66 pg I-TEQ m^–3. Li et al. (2018) reported that the concentration of PCDD/Fs in surrounding atmosphere of MSWIs in China was 0.05–0.12 pg I-TEQ m^–3.

It could be seen that PCDD/F concentrations was high in winter and low in summer in the same sites. The possible factors were the difference of waste composition, temperature conditions. It could be seen that PCDD/F concentrations of the samples collected by the chimney were higher than the concentrations of the samples collected in the boundary of the plants. we could observe a trend of declining PCDD/Fs levels as the distance increased. If PCDD/Fs levels tended to decrease as the distance from a certain site increase, then the site could be considered a source of PCDD/Fs (Oh et al., 2006). Comparing the sample MW01-A3 and HW03-A2, we can see that when the PCDD/F concentrations in the stack gas emitted from the plant were relatively low, the corresponding levels in the ambient air were also lower (Hu et al., 2014).

As could be seen in Fig. 7, the concentration of 2,3,4,7,8-PeCDF in the ambient air was higher than other congener of PCDD/Fs, and the PCDD/Fs congener profiles of them were similar with each other, the result was similar with the previous study by Chen et al. (2015) and Oh et al. (2006). 

Fig. 7. PCDD/F congener profiles in the ambient air of MSWI and HWI.

Fig. 8 shows comparison of PCDD/Fs patterns in ambient air and stack gas, the relative importance of the individual congener was normalized to the total I-TEQ value. Among all PCDD/F congeners, 2,3,4,7,8-PeCDF was the most predominant (21.1%–40.5% to the ∑PCDD/Fs), and the content of PCDFs was higher than PCDDs, Ratios of PCDF/PCDD varied from 0.688 to 0.831, and the less chlorinated PCDD/Fs were the main congeners in the ambient air, these were also true with peer investigation (Oh et al., 2006; Xu et al., 2009; Hu et al., 2014; Chen et al., 2015). Besides, PCDD/Fs patterns in ambient air were similar to the patterns in the stack gas samples. These results indicated that PCDD/F emission from the plants had an obvious impact on the environments surrounding the plants. 

Fig. 8. Comparison of PCDD/Fs patterns in ambient air and stack gas.

Congener Profiles and Levels of PCDD/Fs in Surrounding Soil

Table 4 summarized the content of PCDD/Fs in soils. The I-TEQ values of PCDD/F ranged from 0.68 to 3.76 ng I-TEQ kg⁻¹ in soils in the vicinity of MSWIs, with an average and a median value of 1.73 and 1.51 ng I-TEQ kg⁻¹. The I-TEQ values of PCDD/F ranged from 0.91 to 16.0 ng I-TEQ kg⁻¹ in soils in the vicinity of HWIs, with an average and a median value of 7.93 and 7.41 ng I-TEQ kg⁻¹. The values of TEQ in the soil samples in the vicinity of MSWIs and HWIs were comparable with those in surrounding soil of MSWI in China (0.39–5.04 pg I-TEQ g⁻¹) (Yan et al., 2008) and (7.62–15.5 ng I-TEQ kg⁻¹) (Li et al. 2018), the iron and steel plant (0.16–4.5 ng I-TEQ kg⁻¹) (Zhou et al., 2019), hospital waste incinerator (0.46–2.63 ng I-TEQ kg⁻¹) (Li et al., 2012) and lower than those in the vicinity of electronic wastes disposal sites in Indian (8–99 ng I-TEQ kg⁻¹) (Chakraborty et al., 2018). For the ratio of the PCDF/PCDD (F/D) ranged from 0.05–1.72 with a mean value of 0.52. Among 11 surrounding soil samples, there were 9 samples whose F/D ratio below 1, suggesting combustion was not the mainly source of PCDD/Fs in surrounding soil (Li et al., 2018). The 2 surrounding soil samples (F/D = 1.51, 1.72) from HW03 and HW04, indicating abundant enrichment of PCDFs in these soil sample, and the high TEQ of PCDD/Fs was formed because of the long-term deposition of pollutants in the vicinity of HWIs (Nganai et al., 2014). 

As shown in Fig. 9. For mass concentration of PCDD/Fs, among all soil samples, OCDD was the predominant congener, accounting for 28.3%–77.7% of the total concentration of 17 congeners for soil samples, followed by the high-chlorinated congeners 1,2,3,4,6,7,8-HpCDD, 1,2,3,4,6,7,8-HpCDF, and OCDF, which was consistent with those previously reported by other investigators (Yan et al., 2008; Meng et al., 2016; Zhou et al., 2019; Li et al., 2018). 

Fig. 9. Composition profiles of PCDD/Fs in the soil samples in the vicinity of MSWIs and HWIs.

CONCLUSIONS

In this study, all of the investigated MSWIs exhibited a small variation in PCDD/F emissions to the atmosphere (0.051–0.078 ng I-TEQ m^–3 with an average value of 0.058 ng I-TEQ m^–3), and the emission levels of 5 MSWIs were below the current standard in China (0.1 ng I-TEQ m^–3). However, the emission of PCDD/Fs from the investigated HWIs ranged from 0.016–0.29 ng I-TEQ m^–3 with an average value of 0.159 ng I-TEQ m^–3. Two HWIs had emission levels above the current standard in China, indicating that incineration technologies should be further improved. The emission factors of PCDD/Fs to the atmosphere were in the range of 0.056 to 0.607 µg I-TEQ t^–1 MSW (mean: 0.235 µg I-TEQ t^–1) and 0.07–3.27 µg I-TEQ t^–1 HW (mean: 1.03 µg I-TEQ t^–1), respectively. It was estimated that 10.2–88.0 mg (mean: 35.0 mg) and 1.13–53.0 mg (mean: 17.1 mg) I-TEQ of PCDD/Fs were annually released from MSWIs and HWIs to the atmosphere in China, respectively. Based on principal component analysis, most MSWIs belonged to Group 2, and 4 investigated HWIs could be divided into three groups because of different feeding wastes. Our results indicated that PCDD/F emission from the plants had an obvious impact on the PCDD/F concentration and patterns in ambient air. The I-TEQ values of PCDD/F in soil samples in the vicinity of MSWIs were lower than those around the HWIs.

ACKNOWLEDGEMENTS

This Project was supported by the National Key Research and Development Plan of China (No. 2018YFC1901300) and the Innovative Research Groups of the National Natural Science Foundation of China (No. 51621005) and the National Natural Science Foundation of China (No. 51976192).