Study on Phosphorus Removal Pathway in Constructed Wetlands with Thermally Modified Sepiolite

: Constructed wetlands, as natural sewage treatment ecosystems, have been widely used in the fields of rural domestic sewage and sewage plant tailwater treatment. However, the synchro-nous removal of phosphorus in most constructed wetlands is not ideal. This study aimed to prepare thermally modified sepiolites with high phosphorus adsorption capacities and design a constructed wetland based on them. Thermal modification was adopted to enhance the adsorption capacity of sepiolite based on its high specific surface area and ion exchange capacity. The physicochemical properties and adsorption performance of thermally modified sepiolite were studied. The results showed that the specific surface area and adsorption capacity of thermally modified sepiolite were higher than those of natural sepiolite, reaching 19.494 mg·g −1 . The concentration of effluent and the removal of constructed wetlands based on thermally modified sepiolite was 0.07 mg·g −1 and 91.05%. An analysis of the phosphorus forms in constructed wetlands proved that the main phosphorus removal pathway is the adsorption of substrate, and the form of phosphorus was mainly Ca/Mg−P and Fe/Al−P.


Introduction
With the development of the social economy and large−scale agricultural and industrial production, a large amount of nitrogen and phosphorus pollutants are discharged into water bodies, which seriously endangers production and domestic water safety [1]. At present, the processes commonly used for phosphorus removal are coagulation sedimentation, ion exchange, electrodialysis, and biological methods [2][3][4][5]. However, most of these methods have high operational costs and are not suitable for the large−scale advanced treatment of phosphorus wastewater. Constructed wetlands (CW) have attracted much attention since the 1970s because of their low construction and operation costs, convenient maintenance, and environmental friendliness.
However, the application of constructed wetland technology worldwide has encountered many problems when dealing with nitrogen and phosphorus. Constructed wetlands have a good removal effect on COD (Chemical oxygen demand), BOD (Biochemical oxygen demand), SS (Suspended solids), and other pollutants, but the nitrogen and phosphorus removal effect is low [6][7][8]. According to reports, the nitrogen removal of most constructed wetlands is generally less than 30%, whereas the phosphorus removal of traditional constructed wetlands based on gravel is only approximately 15% [9,10]. The main mechanisms of phosphorus removal in constructed wetlands are natural sedimentation, absorption and the assimilation of plants, excessive absorption of phosphorus by denitrifying phosphate accumulating organisms (DPAO), and the adsorption of phosphorus by substrate, among which the most important part is the adsorption of phosphorus by the substrate [11]. At present, the main method to strengthen phosphorus removal in constructed wetlands is to study new absorbents for phosphate removal, such as fly ash, shale, dolomite, and kaoline [12][13][14][15][16][17]. Studies [14] have shown that slag, oyster shell, concrete and cinder have greater theoretical saturation adsorption capacity for TP (Total phosphorus). Zhen et al. [16] compared the removal effect and mechanism of phosphorus in the anaerobic digestion solution of swine wastewater by using zeolite, broken brick and oyster shell, etc. The experimental results showed that the removal effect of phosphorus by oyster shell was the best when the concentration was about 80 mg·L −1 and the removal was 90%. These adsorbents fixed P in sewage on the substrate surface through adsorption and chemical precipitation and achieved a good phosphorus removal effect, among which sepiolite is representative.
Sepiolite (SEP) is an inorganic magnesium−rich silicate clay mineral with a standard crystal chemical formula of Mg8[Si12O30][OH]4·8H2O [18]. According to the division of mineralogy, SEP is a kind of 2:1 layer of chain−like clay consisting of two layers of silicon oxygen tetrahedron sandwiched with a layer of magnesium oxygen octahedron, and the oxygen at the top angle of the silicon oxygen tetrahedron outside intersects with each other and forms a chain structure in the opposite direction, forming a series of 0.37 nm × 1.06 nm honeycomb pores [19,20]. SEP has a theoretical specific surface area of up to 900 m 2 /g, providing it with a strong adsorption capacity and numerous adsorption sites [21].
Previous studies on SEP have mainly focused on its performance in the laboratory instead of its practical application. Meanwhile, research on constructed wetlands also focuses on phosphorus removal performance. Many studies have shown that the phosphorus removal of CW is about 15-70%, which is greatly affected by the type of substrate material and temperature [22][23][24]. The mechanism of phosphorus removal by constructed wetland substrate was not clear and there were few studies in this field. It is speculated that the removal of phosphorus by constructed wetlands is a comprehensive process of physics, chemistry and biology, which can be roughly divided into three parts: (1) absorption and the assimilation of plants; (2) the normal absorption and excessive accumulation of phosphorus by phosphorus accumulating bacteria; (3) interception, filtration and the adsorption of substrates. However, the contribution of each part to phosphorus removal and how the phosphorus form changes in constructed wetlands are studied relatively infrequently.
Therefore, this paper intended to study the forms and distribution of phosphorus in constructed wetlands to reveal the pathway and mechanism of phosphorus removal in constructed wetlands. In this study, thermally modified sepiolite (T−SEP) was prepared and used to remove phosphorus from the effluent of a sewage treatment plant. Isotherm and kinetic models are discussed to explain the adsorption process. Characterization methods such as XRD (X−ray diffraction), FT−IR (Fourier transform infrared), and SEM (scanning electron microscope) were used to further understand the behavior and mechanism of adsorption. After thermal modification, SEP was used as a substrate material of constructed wetlands, and the removal efficiency of phosphorus, variation in pollutant concentration along the path, and form and distribution of phosphorous were studied to reveal the removal pathway of constructed wetlands based on T−SEP. Therefore, T−SEP was proved to have great potential for removal of phosphorus and constructed wetlands based on it had a strong ability to remove phosphorus.

Materials
Natural sepiolite (N−SEP) powder with a mesh number of 350 was provided by Chengyu Mining Co., Ltd. (Shijiazhuang City, China). All other chemicals in the study, including NH4Cl, HCl, NaOH, NaCl, KCl, CaCl2, Na2CO3, and Na2SO4, were purchased from Sinopharm Chemical Reagent Co., Ltd (Beijing, China)., were of AR grade, and were not subjected to further purification.

Modification of SEP
2.2.1. Purification of SEP N−SEP was added to distilled water at a solid−liquid ratio of 1:10 and stirred in a beaker for 10 min. The solution was poured into a separating funnel and left to stand until stratification was not obvious. The underlying SEP solid was added to distilled water. This process was repeated several times until no stratification occurred. After centrifugation, the supernatant liquid was removed and the SEP was washed several times with distilled water. The solid product was dried in an oven at 120 °C for 24 h. The purified SEP was then ground into powder and sieved through a 200−mesh screen for use.

Thermal Modification
10.0 g of purified N−SEP was placed in muffle furnace and calcined at 700 °C for 4 h (heating rate was 5 °C/min). After natural cooling to room temperature, the samples were washed with distilled water, dried in an oven at 120 °C for 24 h, screened through a 200− mesh sieve, and T−SEP was prepared.

Characterization of T−SEP
The XRD test was carried out using a German Brook D8 X−ray diffractometer under the following conditions: Cu target, step 0.02°, scanning angle range 5-80°, and scanning speed 12°·min − 1 . Fourier transform infrared (FT−IR) spectra were recorded on a Perki-nElmer 1725X FTIR spectrometer over the spectral range of 4000-400 cm − 1 using KBr pellets. The Brunauer-Emmett-Teller (BET) specific surface area was obtained using a Micromeritics ASAP 2020 apparatus, and the porous characteristics of the composites were obtained using theBarrett Joyner-Halenda (Barret−Joyner−Halenda) method through standard nitrogen adsorption/desorption measurements. An FEI Quanta 200 field−emission scanning electron microscope (SEM) was used for morphological analysis.

Adsorption Experiments
An adsorption experiment was conducted by 0.1 g T − SEP in 50 mL simulated wastewater containing phosphorus, and the pH of the solution was adjusted with 0.1 mol·L −1 HCl or NaOH. The mixture was then stirred at 200 rpm for the desired time at 25 °C. After centrifugation at 5000 rpm for 10 min, the supernatant was filtered through a 0.45 μm membrane, and its concentration of phosphorus was analyzed using the molybdate photometric method. The amount of phosphorus adsorbed by the adsorbent (qe, mg·g −1 ) was calculated using Equation (1): Additionally, the removal of phosphorus was estimated using Equation (2): C0 and Ce represent the initial and final concentrations (mg·L −1 ) of ammonium ions in the solution, respectively. V represents the volume of the solution (mL), and m represents the dose of the adsorbent (g).
The experimental data were fitted to the pseudo−first−order, pseudo−second−order, and Elovich models, as shown in Equations (3)- (5).
Pseudo−second−order model: q = (4) where t is the adsorption time (min), qe and qt are the amount of nitrogen adsorbed on the adsorbent (mg·g −1 ) at equilibrium and at time t (min), respectively, and k1 (1·min − 1 ) and k2 (1·min − 1 ) are the adsorption rate constants of the pseudo−first−order and pseudo− second−order models, respectively. Adsorption isotherm studies were conducted at different initial concentrations (0.5, 1, 2, 5, 10, 15, 20, 30, 40, 50, 80, 100, 150, and 200 mg·L −1 ). The other operating parameters were as follows: 50 mL of adsorption solution with 0.1 g adsorbent, the adsorption time was maintained at 10 h, the initial pH was fixed at 7.0 ± 0.1, and the magnetic stirring speed was 200 rpm.
The equilibrium data were fitted to the Langmuir, Freundlich, and Sips isotherm models, as shown in Equations (6)-(8): Freundlich model: q = K C Sips model: q = Ce (mg·L −1 ) and qe (mg·g −1 ) are the nitrogen and phosphorus concentrations, respectively, which are the corresponding adsorption capacities at equilibrium. qm (mg·g −1 ) is the maximum adsorption capacity, and KL is the Langmuir adsorption equilibrium constant, which is related to the adsorption capacity. KF and n are constants in the Freundlich isotherm model that are related to the adsorption capacity and intensity, respectively.

Design of Constructed Wetlands
As shown in Figure 1, the system consisted of two horizontal subsurface wetlands, each of which was rectangular polyvinyl chloride(PVC) tanks with the same dimensions(length: 1.80 m, width: 0.30 m, height: 0.50 m). The devices were divided into three areas: water distribution, treatment, and catchment. A peristaltic pump was used for inlet water. The length of the distribution area and catchment area was 0.15 m, and the gravel size was 30-50 mm. The length of the treatment area was 1.50 m, the thickness of substrate was 0.40 m, and the height of the running water level was 0.35 m.
The substrates and particle sizes in the two constructed wetlands were as follows: (1) W−G was 1-2 mm gravel and (2) W−S was 1-2 mm T−SEP. Iris pseudacorus L., which has developed a root system, strong resistance to pollution and clean pollution, and salt and alkali resistance [25], was selected as a wetland plant in this experiment.
The experimental device was built in the ecological greenhouse on the 4th floor of the Environment Building in Hohai University, Nanjing, China (32°03′44″ N, 118°46′01″ E). The device was built and tested in May, and the stable operation began from June to May, during which the temperature varied from 2 °C to 37 °C. The inlet water of the constructed wetland was obtained from the effluent of the Jiangxinzhou sewage treatment plant in Nanjing. The relevant water quality parameters are listed in Table 1. Four different hydraulic loads were adopted, and one month was taken as an inlet cycle. The inlet flow rates were designed as 10.98, 14.63, 21.95 and 43.90 L/d, respectively, and the corresponding hydraulic retention times were four, three, two and one day, respectively. During the operation of the device, the continuous flow mode was adopted, and a peristaltic pump was used to maintain the constant inflow rate of the inlet water.

Samples Analysis
Every three days, 50 mL of water was collected from each sampling point at 9:00 a.m. After passing through a 0.22 μm filter membrane, various water quality indicators were detected. The concentration of TP was analyzed using the molybdate photometric method.
In this study, the method of phosphorus form analysis in fillers suggested by Forbes [26] was adopted. The wet sample of the same weight as the dried sample was placed into the centrifuge tube, and an appropriate amount of leaching solution was added for leaching. At the same time, parallel samples were taken and dried to the balance weight, and the moisture content was measured. The remaining sample was added to the next leaching solution and was extracted in successive steps. After leaching, the sample was centrifuged at 3000 r/min for 10 min to analyze the organophosphorus and total phosphorus in the leaching solution, and the remaining sample was added to the next leaching solution and extracted in successive steps. In addition, the unused substrate was used as blank samples, and the background values of various forms of phosphorus in each filler were determined by the aforementioned method. The increment of each form of phosphorus in substrate was then obtained by subtracting the corresponding background values from the measured values of each form of phosphorus.
The plants in this experiment were harvested, washed, air dried, cut into small segments, and placed in the oven to dry to the balance weight. The biomass dry weight of the plants was then weighed. The dried samples were crushed and passed through a 100−mesh sieve. After digestion with H2SO4 and H2O2, the phosphorus content in the plant was determined by molybdate photometry.
After the operation of constructed wetlands, substrates with a depth of 10-15 cm were collected at each water intake. Collected samples were extracted using a bacterial genomic DNA extraction kit (Omega Bio−tek, Norcross, GA, USA). The bacterial 16S rRNA was then subjected to PCR amplification using the bacterial fused primers 338F (5′ − ACTCCTACGGGAGCCAGC − 3′) and 1492R (5′ − GGTTACCTTCTTACGACTT − 3′). The amplified products were detected by electrophoresis on 1% agarose gel. Purified DNA was sequenced using the Illumina Mi−Seq platform.

Statistical Analysis
Microsoft Excel 2016 and IBM SPSS Statistics 23.0 were used to analyze the experimental data, and three stable values were recorded for each measurement. The mean value and standard deviation were calculated.

Optimum Parameters for Preparation of T−SEP
The adsorption capacity of phosphorus by T−SEP calcined at different temperatures was shown in Figure 2a. According to the figure, the adsorption capacity of T−SEP increased with the increase of temperature and reached the maximum at 700 °C (18.634 mg/g). When the calcination temperature was below 300 °C, the morphology and structure of SEP hardly changed. Between 300 °C and 500 °C, the free water and impurities in SEP evaporated due to calcination, and the channels in SEP crystal were unblocked. As the temperature continued to rise to 700 °C, the crystal water and binding water in SEP crystal were removed, and the adsorption capacity was further improved. However, when the temperature exceeded 700 °C, the phosphorus adsorption capacity did not increase much, because the impurities and water in SEP have been almost completely removed at a lower temperature, and further increasing the temperature destroyed the chain layer structure of SEP instead of improving its adsorption performance, which is similar to that of previous studies [27].  Figure 2b showed the adsorption capacity of T−SEP at a calcination temperature of 700 °C and a different calcination time (one, two, three, four, and five hours). It could be observed form Figure 2b that prolonging the calcination time was beneficial to improving the adsorption performance of T−SEP before the calcination time reached 4 h. The adsorption capacity would be decreased in case the calcination time exceeded 4 h. The calcination temperature of 700 °C and calcination time of four hours were selected to be the optimum parameters for subsequent experimental research.

Characterization of T−SEP
X−ray diffraction diagrams of N−SEP and T−SEP are shown in Figure 3a. The diffraction pattern for SEP has peaks at 7.21°, 20.74°, and 26.48°, corresponding to the (100), (131), and (080) planes of SEP according to the PDF card [13−0595]. Compared with N− SEP, 2θangel at the characteristic diffraction peaks of T−SEP remains unchanged, which proves that modifications can generally preserve the structure of SEP [28]. However, the intensity of the peaks at T−SEP was reduced, which could be ascribed to the calcium ions entering the SEP ion channels, and calcination changes the structure of some channels. The peaks appearing at 28.89° can be indexed to CaCO3 (JCPDS [85−1108]). The characteristic diffraction peaks of CaCO3 at T−SEP were significantly enhanced, which suggests that Ca was successfully loaded onto the surface of the SEP. Active calcium carbonate covers the surface of SEP, forming active adsorption sites and carbonate hydroxyapatite precipitation with phosphate [29]. An FT−IR analysis was performed to verify the change in SEP during the modification process. As shown in Figure 3b, the FT−IR spectra of SEP before and after modification changed slightly, indicating that the modification process did not change the basic structure of the adsorbent. The peaks at 3694 cm − 1 , 3693 cm − 1 and 3695 cm − 1 were caused by the tensile vibration of hydroxy (V−OH) connected to octahedral Mg 2+ , which is the structural water of SEP [30,31]. The peaks at 3437 cm − 1 , 3420 cm − 1 , 3427 cm − 1 and 3438 cm The specific surface area and pore size distribution of the samples were calculated using the multi−point BET method, and the pore volumes of the samples were calculated using the adsorption capacity at a relative pressure P/P0 = 0.99. The pore size distribution and surface properties of the SEP are shown in Figure 3c,d and Table 2, respectively. Table  3 shows that the specific surface area of T−SEP increased from 12.84 m 2 ·g − 1 to 30.41 m 2 ·g − 1 , and the average pore size increased from 7.37 nm to 12.65 nm; the pore volume increased from 0.0083 cm 3 ·g − 1 to 0.0209 cm 3 ·g − 1 , proving that the specific surface area, pore size and pore volume of the T−SEP was improved, which gives it a stronger adsorption capacity.
As shown in Figure 3c,d, the pore size distributions of SEP and T−SEP are mainly mesopores between 2-50 nm, with a small number of micropores and macropores. Compared with N−SEP, the number of micropores in T−SEP decreased, whereas the number of macroporous and mesoporous pores increased, which could be ascribed to calcination. Owing to the removal of adsorbed, bound, and crystal water, the number of mesopores and macropores increased. As seen from Figure 3e,f, SEP and T−SEP was an aggregation of fiber bundles, needles, hairs, and sheets, and was covered with impurities, such as talc and calcite, which may affect its adsorption capacity. The surface of the T−SEP was relatively smooth, and the impurities were reduced. The fiber bundle, needle, hair, and flake aggregate structures were looser. The linear structure of T−SEP was broken, and the aggregate was significantly reduced because the surface structure of SEP partially collapses as a result of the thermal modification. Figure 4a shows the pseudo−first order, pseudo−second order and Elovich's kinetic model of the adsorption process of phosphorus by T−SEP. The relevant kinetic model parameters are shown in Table 3. It is obvious that the adsorption rate of phosphorus was very rapid in the first 100 min, yet decreased and gradually tended toward equilibrium after 200 min, which can be explained by the concentration difference in the solution. The adsorption process can be divided into three main steps: (1) external particle diffusion, (2) internal particle diffusion, and (3) adsorption reaction. The initial controlling step of the adsorption process was the external diffusion stage of the particles, the concentration of adsorbent in the solution was high, and the concentration of phosphorus on the surface of T−SEP was 0. The concentration difference between the solid and liquid phases formed the driving force, and phosphorus in the solution migrated to the liquid film on the surface of the adsorbent, hence the adsorption rate was rapid in this stage. As adsorption progressed, the concentration of phosphorus between the solid and liquid reached equilib-rium, and the controlling step of adsorption was the internal particle diffusion and adsorption reaction. The reaction rate of these two steps was slower than that of the first stage, thus the adsorption rate also decreased [34].

Adsorption Isotherm
As shown in Figure 4b, the Langmuir, Freundlich, and Sips isotherm models were used to describe the adsorption process of phosphorus by SEP, and the relevant parameters were listed in Table 4. The maximum adsorption capacities for phosphorus was 19.494 mg·g −1 . Compared with the Freundlich model (R 2 was 0.974), the adsorption processes of phosphorus were more consistent with the Langmuir and Sips models (R 2 was 0.989 and 0.993), which indicated that the adsorption of phosphorus by SEP may be dominated by monolayer and chemical adsorption. The 1/n parameter in the Freundlich model reflected the degree of favorable adsorption. When n > 1, adsorption was preferential, and, it was opposite when n < 1 [35]. The n values of SEP in this experiment were greater than 1, which indicated that the adsorption of phosphorus was preferential.

Effect of Hydraulic Retention Time on Total Phosphorus Concentration
The concentration of TP in inflow and outflow of the constructed wetlands based on gravel and SEP and its removal rate at different hydraulic retention times were shown in Figure 5a. The average effluent concentration of TP was 0.39~0.64 mg·L −1 in W−G and 0.07~0.24 mg·L −1 in W−S. It was apparent that removal of TP in W−S was superior to that in W−G, which is related to the excellent adsorption performance of T−SEP. It could be seen that removal of TP was positively correlated with HRT, because removal of phosphorus was related to the adsorption of substrate in the constructed wetland. Therefore, the reduction of HRT reduced the contact time of the substrate and microorganism to the pollutants. However, excessive HRT may also cause sewage retention and the consumption of DO in the constructed wetlands, resulting in the anaerobic state of the system Additionally, the adsorption capacity of phosphorus reached saturation, and the adsorption was difficult to improve further. Therefore, in the actual operation of the constructed wetlands, the appropriate HRT should be selected according to the specific situation. In this experiment, four days was the best.

Variation of Total Phosphorus along the Way
The variation of TP concentration along the way of W−G and W−S is shown in Figure  5b. The total phosphorus was rapidly adsorbed and degraded in the first quarter of the system, which decreased from the initial 0.48 mg·L −1 to 0.16 mg·L −1 , accounting for 72.73% of TP removal in the constructed wetlands. This indicated that in the front of the reactor, the absorption of substrate, the absorption and transformation of plants and the assimilation of microorganisms happened very violently. The TP removal in the second half of the constructed wetlands was mild, mainly to ensure the stability of effluent water quality.

Variation of Phosphorus Removal in Four Seasons
The variation of phosphorus removal in different seasons is shown in Figure 5c. It can be observed that the removal efficiency in W−S was obviously higher than that in W−G, which should be due to the adsorption of SEP. Another noteworthy phenomenon is that the phosphorus removal effect of W−G had an obvious seasonal effect, but the effluent of W−S was very stable. The phosphorus removal efficiency of W−G in winter was only 23.49%, which was only 1/3 of that in summer (63.28%). This phenomenon implied that the plant absorption and enrichment is one of the main pathways of phosphorus removal in W−G, and phosphorus removal is greatly affected by plant death in winter. Conversely, plants account for only a small part of the phosphorus removal by W−S, so the effect of phosphorus removal in winter is only slightly reduced

Spatial Distribution and Form of Phosphorus in Substrate
After one year of operation, the adsorption amount of phosphorus in the substrate basically reached saturation. After sampling, a continuous hierarchical extraction analysis was carried out. The experimental results are shown in Figure 6a. In W−G, the forms of phosphorus were mainly water−soluble phosphorus and organophosphorus, which indicated that phosphorus removal was mainly achieved through the biological phosphorus removal by DPAO and the surface adsorption of substrate in this system. However, the content of phosphorus in W−S was significantly higher than that in W−G, and the form of phosphorus was mainly calcium−magnesium bound phosphorus and iron−aluminum bound phosphorus, which further proved that the adsorption of phosphorus by substrate was the main way of phosphorus removal in constructed wetlands, and the co−precipitation of calcium and phosphorus played a key role in phosphorus removal.
In W−G and W−S, the contents of Org−P decreased along the way, and the contents of Org−P in the first quarter of the constructed wetlands were 5.46 and 3.28 mg·kg −1 , respectively, while the contents of Org−P in the last quarter were only 0.86 and 0.68 mg·kg −1 . When sewage was flowing through the front of the constructed wetland, the sewage gradually went through aerobic and anoxic processes, so the concentration of organophosphorus decreased rapidly. The contents of Ca/Mg−P in W−G and W−S along the way were nearly constant, which was about 13~17 and 38~43 mg·kg −1 . According to the spatial distribution characteristics of phosphorus, it can be seen that the phosphorus removal of the constructed wetland substrate was saturated.

Analysis of Phosphorus Content in Plants
According to Figure 6b, plant total phosphorus content in W−G and W−S had an obvious seasonal effect, and plant total phosphorus content in W−S was slightly higher than that in W−G. The plant total phosphorus content of W−G in summer, autumn, winter and spring was 1243, 312, 234 and 847 mg, respectively, and that of W−S was 1243, 1357, 357, 198 and 921 mg, respectively. The phosphorus content of plants varied greatly in different periods. This phenomenon may be due to the different ways of phosphorus removal in different seasons. In winter, large amounts of phosphorus was transferred to the roots to conserve nutrients; In the summer and autumn, phosphorus was transferred to the stems and leaves for growth.

Analysis of Phosphorus Removal Pathway in Constructed Wetland
There are two main methods of phosphorus removal in constructed wetlands, which are adsorption and precipitation of substrate and uptake and absorption of plants. Microbial synthesis also plays a role in the initial stage of initiation. Among them, the absorption and assimilation of plants synthesizes phosphorus into organic components such as ATP, DNA and RNA, which can be partially removed from the system by plant harvesting. Microbial synthesis includes the normal assimilation and absorption of phosphorus by microorganisms and the excessive uptake of phosphorus by phosphorus accumulating bacteria.
The system phosphorus removal pathway analysis is shown in Table 5. Other methods of phosphorus removal in the table include system interception of phosphorus removed by suspended matter and the excessive uptake of phosphorus by phosphorus accumulating bacteria normally assimilated and absorbed by phosphorus stored in the underground part of plants. According to Table 5, the main phosphorus removal pathway in W−G and W−S is adsorption of substrate, and the proportion is 67.2% and 78.2%, respectively. The phosphorus removal by substrate and total phosphorus removal are both significantly higher than those in W−G, which can be a ributed to the superior adsorption of T−SEP in W−S. The phosphorus removal by plants accounts for about 15% of the total phosphorus removal, which cannot be ignored, and has a significant seasonal effect. However, the CW based on T−SEP has a more stable operation effect throughout the year, and still has a better removal effect in winter by virtue of substrate adsorption.
As shown in Figure 7, The distribution of influent phosphorus mainly included effluent, plant absorption, microbial absorption and accumulation, substrate adsorption and sedimentation. Phosphorus sedimentation in CW referred to the process of phosphorus storage in wetland through the physical action caused by soluble phosphate in water. Soluble phosphate (HPO4 2− and H2PO4 − ) was absorbed and assimilated into the organic components of plants. The potential rate of phosphorus uptake by plants was limited by their net growth and the mass fraction of phosphorus in plant tissues. The storage capacity of phosphorus in plants depended on the mass fraction of phosphorus in plant tissues and the ultimate biomass accumulation potential. This study emphasized the effect of seasons, and obviously the plants were more vigorous in spring and summer than in autumn and winter, during which temperature is the main controlling factor.
The effects of microorganisms on phosphorus included normal absorption and excessive accumulation. The normal absorption of phosphorus by microorganisms can be divided into two types: autotrophic microorganisms and heterotrophic microorganisms, the former requiring high water temperature and sufficient sunlight, the latter requiring appropriate organic carbon sources. The excess accumulation of phosphorus depended on the over−uptake of phosphorus in aerobic environment and the release of phosphorus in anaerobic environment by DPAO, which usually absorbed much more phosphorus than it released.
From the material accounting of phosphorus, it could be seen that the adsorption of phosphorus by substrate was the main removal pathway of constructed wetland. Most of the phosphorus adsorbed by the substrate was stored in the form of Fe/ Al−P and Ca/ Mg−P. Compared with other material, the content of Ca/Mg−P in SEP was higher, which was mainly due to Mg 2+ in SEP crystals and Ca introduced by thermal modification. Numerous studies have shown that Proteobacteria contain most of the bacteria with a nitrogen removal function, and Proteobacteria are the dominant flora in the four constructed wetlands; therefore, the constructed wetlands in this experiment have a good nitrogen removal function. Proteobacteria, Bacteroidetes, Actinobacteria, and Acidobacteria were the dominant bacteria responsible for nitrogen removal. Bacteroidetes are a large family of natural chemotrophic bacteria that mainly play a role in COD degradation. Some denitrifying bacteria and most nitrogen−fixing bacteria also belong to this type of microorganism and function in nitrogen removal. These bacteria also play an important role in nitrogen removal in the constructed wetlands. In addition, it is worth noting that the relative abundance of Acidobacteria in A02 is significantly higher than that in other samples. Combined with the best decontamination effect of W−S in previous experiments, we can speculate that Acidobacteria have an important contribution to the process of pollutant removal in constructed wetlands. Figure 8b shows the top 10 species with relative abundances at the genus level in W−G and W−S. It can be seen that Bacteroidetes_vadinHA17, Anaeromyxobacter, Geobacter, Subgroup_7, MBNT15, UTCFX1, Subgroup_6, SC−I−84, Anaerolinea and KD4−96 were the dominant strains. However, the relative abundance of 17 genera of Bacteroides with the highest proportion was only 1.90-5.66%, while the 10 genera with the highest relative abundance only accounted for 18.6-25.01% of the microbial community, and there were no obvious dominant species, which indicated that the diversity of constructed wetlands at the genus level was high, and the differences between different constructed wetlands at the genus level were greater. This is similar to previous studies [36,37].

Conclusions
In this study, T−SEP was prepared, and its physicochemical and adsorption properties were studied. After modification, the structure of the SEP remained unchanged, the surface became smoother as the impurities were reduced, and the specific surface area slightly increased. The optimum calcination temperature was 700 °C and the optimum calcination time was 4 h. The maximum adsorption capacity of T−SEP for total phosphorus was 19.494 mg·g −1 . T−SEP was used as a substrate of the constructed wetlands and the removal of phosphorus from the tail water of the sewage treatment plant. The concentration of effluent and the removal of constructed wetlands based on T−SEP was 0.07 mg·g −1 and 91.05%. After analysis of the constructed phosphorus forms in constructed wetlands, it can be concluded that the main phosphorus removal pathway is the adsorption of substrate.
Funding: This study was financially supported by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).

Institutional Review Board Statement: Not applicable
Informed Consent Statement: Informed consent was obtained from all subjects involved in the study.
Data Availability Statement: Not applicable.

Conflicts of Interest:
The authors declare that they have no known competing financial interests or personal relationships that could have influenced the work reported in this study.