3.1. The impact of phytoremediation conditions on the CECs removal efficiency
The efficiency of phytoremediation depends to the greatest extent on the structure and chemical properties of the pollutants removed and the dynamics and nature of the plant's metabolism. The well-being of the plant used, and thus the purification results, can also be influenced by the conditions of the process. In the case of the system using a floating plant, the pH of the culture medium, the daily light exposure time, and the amount of the plant per volume of purified water were considered important for the course of the process. The experimental design chemometric approach was applied to assess the impact of selected parameters on the treatment efficiency. Such an approach reduces the number of experiments necessary for selecting optimal conditions. It saves both time and energy, which meets the assumptions of green chemistry. The characteristics of the course of experiments and the values of removal efficiency of TRC, DES, and E2 obtained for successive sets of conditions are listed in Supplementary Table S3. Based on this data, separate mathematical models were developed to characterize the relationships between the RE% and the pH, the mass of plants, and the time of exposure to the light. The statistical analysis (ANOVA) results of regression models generated for TRC, DES, and E2 are included in Supplementary Tables S4-S6.
The ANOVA results indicate that the obtained regression models are characterized by determination coefficients (R2) values of 0.88, 0.76, and 0.63 for DES, TRC, and E2, respectively. These values indicate that 88% of the variations in DES removal efficiency can be explained by changes in the independent variables selected for the construction of the model, and the model does not explain 12% of the observed changes. In the case of TRC and E2, the applied model explains 76% and 63% of the differences in RE% obtained during the experiments, respectively. The adjusted R2 values, which equal between 0.50 and 0.77, prove the high significance of the models. The optimization process was carried out for a mixture of pollutants, which causes more errors in the models but reflects the real conditions where pollutants do not occur individually.
The Pareto analysis shows the importance of independent variables and interactions in the developed models 32. Pareto charts obtained based on statistical analysis for analyzed EDCs are shown in Fig. 1.
On the graphs, only those parameters whose impact on the efficiency of removing is statistically significant (p < 0.05) were included. In the case of each EDC, all three considered parameters, i.e. lighting time, pH, and amount of plant, were significant. The obtained data indicate that the removal efficiency depends on linear values of the independent variables and quadratic values and interactions between them. In the case of DES and TRC, the greatest influence on the obtained RE% value is the weight of the plant per volume unit of the purified solution. The daily light exposure time was indicated as a factor whose significance is second in the case of TRC. pH of purified solution is the second most important factor in removing DES. The RE% of E2 depends the most on pH value, followed by light exposure time.
Based on CCD models that only considered statistically significant data, response surface plots (RSP) were generated. The RSPs show the influence of selecting two optimized parameters on the predicted removal efficiency. Figure 2 presents the predicted RE% values from the CCD models as a function of daily light exposure and pH.
The simultaneous impact of pH and plant density, as well as the plant density and light exposure time on CECs removal rates, are included in Supplementary Fig. S1 and Supplementary Fig. S2.
Effect of the mass of plant
Plant organisms, together with microorganisms coexisting with them, are, to the greatest extent, responsible for removing pollutants during phytoremediation. Organic micropollutants in water are adsorbed on the surface of floating plants and then collected inside the plant organism, where they undergo bioaccumulation or biodegradation. Therefore, it can be expected that the greater the mass of plants in the studied system, the greater the efficiency of the purification processes 33. The RE value increased with an increase in plant density from 0.5 g to 2.1 g per 100 mL of the purified solution was observed for examined CECs. Increasing plant density above these values does not increase or even decrease the observed degree of removal. This is related to the fact that density affects the growth and quality of plants 34. Too many neighboring plants disturb reproduction, biomass accumulation, and morphology due to competition between individuals 35. Plant density equal to 2 g per 100 mL of the purified solution was selected as the most optimal.
Effect of the light exposure time
The length of the light period in the diurnal cycle of 24 h is an important environmental signal for plants. They have developed mechanisms to measure the length of the photoperiod. This mechanism enables plants to synchronize developmental processes, such as the onset of flowering, with a specific time of the year. It is important in regulating responses to abiotic and biotic stresses as well as the redox state of plants 36. During the conducted research, the influence of light exposure time in the range of 10 to 16 hours a day on the results was examined. It turned out that extending the light exposure time improves the efficiency of CECs removal, but only to the level of 13 hours a day. Further extending the duration of light action, although it promotes the intensive course of life processes, not only does it not bring benefits but also reduces the effectiveness of processes. It can be assumed that excessive exposure to light causes abiotic stress, which disturbs metabolic processes and thus reduces the removal potential of CECs by floating plants. Finally, the photoperiod equal to 13 h per day should be considered optimal from the point of view of the system's efficiency.
Effect of the pH of purified solution
pH affects all chemical and biological processes that occur in water. This is one of the most important factors limiting the distribution and welfare of plants in aquatic habitats. Different species thrive in different pH ranges, with the optimum value for most aquatic plants being 6.5-8. Various environmental and anthropogenic factors can contribute to lowering or upgrading the pH of water outside the optimal range. The acidifying factors are acid-generating soils and rocks, industrial and agricultural wastewater, landfill leachate, and atmospheric acid precursors. High pH is less common than low pH in natural waters, as anthropogenic sources are more often acidic than alkaline. Alkalization of water may occur under the influence of alkaline rocks and soils and runoff from the production and use of asphalt, lime, cement, and soap. The influence of pH on the results of EDC removal was assessed in the range of 5 to 9. In the case of the considered substances, changes in pH impacted the RE values obtained, which turned out to be the highest when the pH of the culture solution was closest to neutral. It can be concluded that such conditions are optimal for small floating plants. Both too-acidic and too-alkaline living environments negatively affect aquatic plants, causing problems with osmoregulation, tissue damage, reduced growth and reproduction 37. Based on the tests, a pH equal to 7 was chosen as optimal.
3.2. Effectiveness and kinetics of CECs removal by floating plants
The course of W. arrhiza and L. minor in removing BPA, DEET, DES, TRC, E1, and E2 from the laboratory mineral culture medium was examined. For this purpose, the floating plant cultures were set up according to the procedure in the experimental section. The initial concentration of each CEC was 100 µg/L in the first research cycle. In addition, a second test cycle was carried out in the case of TRC, BPA, DES, and E1, in which the concentrations were 500 µg/L. Current concentrations of CECs were monitored using USAEME-GC/MS method sequentially after 6 hours and 1, 2, 3, 5, 7, and 14 days of W. arrhiza or L. minor cultivation in the first cycle of research and after 1, 3, 5, and 7 days of experiment in the second cycle of research. The obtained average reduction profiles are shown in Fig. 3. A gradual loss of each CEC's concentration was observed during the cultivation of both plants.
The results obtained in the first research cycle (see Fig. 3A, C0 = 100 µg/L) indicate that an intensive reduction of the pollutants concentrations occurred during the first hours of the experiment. After 6 hours of contact with the medium solution with plants, the concentration of the tested compounds decreased by 11–68% in the case of W. arrhiza and 5–32% in the case of L. minor. The reduction in concentration after one day ranged from 45 to 81% in the system with W. arrhiza and from 19 to 70% in the system with L. minor. Seven days of contact of the solution with the plants removes 88–98% of the initial amount of CECs in the case of W. arrhiza and 87–97% in the case of L. minor. After 14 days of the experiment, these values are 93–99.6% and 89–98%, respectively. The results obtained after the first day of the second experiment cycle (see Fig. 3B, C0 = 500 µg/L) indicate a similar or even higher reduction in concentrations compared to the first cycle in the case of W. arrhiza (62–76%) and slightly lower in the case of L. minor (12–57%). After seven days, the degree of concentration reduction is similar to the results obtained during the first experiment (94–96% W. arrhiza; 84–93% L. minor. By comparing the graphs ln C = f(t) and 1/C = f(t) it was determined that the rate of reduction meets the model of the pseudo-first-order kinetics. Therefore, the CECs removal rate constant (k) was determined based on the formula:
C t = C0e−kt, (4)
where Ct and C0 are the pollutants concentrations at time t and t = 0, respectively, (µg/L). The k values (day− 1) calculated for subsequent time intervals of the conducting experiments are summarized in Supplementary Tables S7 and S8. The relevant half-life values are reported in Supplementary Table S9. In the case of experiments with W. arrhiza, the highest CECs removal rate is observed during the first day of contact of the plant with the solution, and the observed k1 values range from 0.59 day− 1 to 1.67 day− 1. The k1 values for experiments with L. minor are in the range of 0.12 day− 1 − 1.20 day− 1, and the highest removal rate is observed for most analytes on the second day of contact with the plant. The removal rate of CECs by W. arrhiza on days 0–7 ranges from 0.30 day− 1 to 0.57 day− 1 in the case of a lower concentration of compounds (100 µg/L) and from 0.40 day− 1 to 0.46 day− 1 in the case of a higher concentration (500 µg/L). In the case of L. minor, the observed ranges are as follows: 0.29 day− 1 − 0.98 day− 1 and 0.26 day− 1 − 0.38 day− 1 for lower and higher concentrations, respectively. The results indicate that W. arrhiza immediately and without problems adapts to the presence of CECs in the medium at concentrations of 100 and 500 µg/L. In the case of L. minor, the increase in pollutant concentrations initially reduces the intensity of CECs removal. Still, the plant adapts well to the prevailing conditions in the following days. Similar behavior of both plants was observed during previous research 21,22.
3.3. Removal mechanism investigation
Removal of CECs during the contact of polluted waters with plants in the applied conditions of the experiment, as well as real conditions prevailing in nature, is the sum of biological and abiotic processes. Apart from uptake by the plant followed by bioconcentration and/or biodegradation, the removal of pollutants may also be caused by hydrolysis, photolysis, and sorption (evaporation can be omitted due to the low volatility of the tested compounds) 24,38. Additional experiments were performed under different conditions to determine the contribution of individual mechanisms to the removal of CECs. Placing the tested solutions in the dark without plants allows the evaluation of the kinetics of the hydrolysis process. When the solution is placed under lighting conditions but without plants, the removal of CECs is the sum of hydrolysis and photolysis processes. Determination of the total share of sorption, hydrolysis, and photolysis is possible when experiments are conducted with a dead plant with access to light. In studies conducted with live plants in the light, hydrolysis, photolysis, sorption, and uptake into the plant is responsible for removing CECs. The course of experiments to determine the share of individual processes in the total effect is presented in Table 2.
Table 2
Mechanisms responsible for the removal of CECs and the course of their determination
Experiment | T (°C) | Light | Plants | Removal mechanisms |
I | 22 ± 2˚C | Yes | Alive | Hydrolysis, photodegradation, adsorption, plant uptake |
II | 22 ± 2˚C | Yes | Dead (5days exposure for 2 mg/L NaN3) | Hydrolysis, photodegradation, adsorption |
III | 22 ± 2˚C | Yes | No | Hydrolysis, photodegradation |
IV | 22 ± 2˚C | No | No | Hydrolysis |
The rate constants of individual processes were determined based on the following relationships 39: kuptake = kI – kII; ksorption = kII – kIII; kphotodegradation = kIII – kIV; khydrolysis = kIV (where kI, kII, kIII, kIV are the rate constants determined for experiments I, II, III, and IV). The detailed data of the rate constants for all obtained removal mechanisms are gathered in Table 3.
Table 3
The values of the plant uptake, sorption, photodegradation, and hydrolysis rate constants
CECs | k (day− 1) |
uptake | sorption | uptake | sorption | photo-degradation | hydrolysis |
W. arrhiza | L. minor |
DEET | 0.385 | 0.038 | 0.023 | 0.131 | 0.029 | 0.113 |
TRC | 0.385 | 0.059 | 0.410 | 0.022 | 0.006 | 0.124 |
BPA | 0.109 | 0.269 | 0.251 | 0.083 | 0.061 | 0.123 |
DES | 0.753 | 0.060 | 0.559 | 0.254 | 0.005 | 0.169 |
E1 | 0.105 | 0.327 | 0.344 | 0.055 | 0.004 | 0.101 |
E2 | 0.055 | 0.229 | 0.074 | 0.192 | 0.050 | 0.155 |
Based on the averaging of the determined values, it can be concluded that the largest contribution to the removal of CECs is the plant uptake, with the mean values of kuptake equal to 0.299 day− 1 and 0.277 day− 1 for W. arrhiza and L. minor, respectively. Sorption is the dominant mechanism for removing BPA, E1, and E2 by W. arrhiza and DEET and E2 by L. minor. It should be emphasized that the processes of sorption and plant uptake are closely related because the pollutant must first be adsorbed on the surface of roots or leaves to be introduced into the plant organism 23. The average ksorption values are about half that of kuptake and are 0.164 day− 1 and 0.123 day− 1 for W. arrhiza and L. minor, respectively. Hydrolysis's contribution in removing CECs is similar to the contribution of sorption; the average value of khydrolysis for the tested compounds is 0.131 day− 1. Studies have shown that DEET, TRC, BPA, DES, E1, and E2 are very slightly susceptible to photodegradation; the rate constant of this process (kphotodegradation) equals only 0.026 day− 1.
As part of the mechanism assessment of CECs removal, the concentrations of the tested compounds in the plant material were determined after a 7-day contact with the enriched nutrient solution containing 100 µg of each compound per 1 g of W. arrhiza or L. minor, following the procedure described in subsection 2.5. The experiment conducted in this way allows for the assessment of what part of the tested compound that has been adsorbed and taken up by the plant remains accumulated in the plant. The results of the experiment are shown in Fig. 4.
As can be seen in Fig. 4, no DES was detected in the plants, while the content of the remaining ECs ranged from 0.15 µg/g (BPA) to 8.17 µg/g (E2) in the case of W. arrhiza and from 0.19 µg/g (TRI) to 12.23 µg/g (E2) in the case of L. minor. Based on the rate constant data listed in Table 3, it can be calculated that, under the conditions of the experiment, the total amount of individual CECs that were adsorbed and uptake by plants ranges from 58 µg/g (E2) to 82 µg/g (DES) in the case of W. arrhiza and from 52 µg/g (DEET) to 82 µg/g (DES) in the case of L. minor (because the sum of ksorption and kuptake constitutes from 58–82% and from 52–82% of the sum of the rate constants of all processes, for W. arrhiza and L. minor, respectively). The obtained results indicate that the plant takes up the vast majority of CECs after adsorption and that only a small part of the collected organic compounds is accumulated in their original form. Organic compounds introduced into a plant cell may be immobilized there unchanged and/or undergo enzymatic modifications, enzymatic degradation, or conjugation with other compounds, mainly glucose and glutathione 23. The uptake of organic compounds into plants and their further fate is affected by many factors, such as chemical hydrophobicity, molecular ionization, and the tendency of the compound to undergo sorption, which precedes the introduction of the compound into the plant 40,41. The values of the logarithm of the tested CECs partition coefficients in the n-octanol/water system (log Kow) are 2.02, 4.76, 3.32, 5.07, 3.43, and 3.94. The logarithm values of the soil organic carbon adsorption coefficient (log Koc) are 1.97, 4.28, 3.18, 4.68, 3.16, and 3.52 for DEET, TRC, BPA, DES, E1, and E2, respectively 42–45. In the studied group, DES is the compound with the highest log Kow and log Koc values. This compound is characterized by the highest level of sorption and uptake by plants (82% of total removal), and the most intense transformations in the plant (no presence in unchanged form in the tissues of W. arrhiza and L. minor). DEET, characterized by the lowest hydrophobicity and tendency to sorption, is much less subject to both uptake by the plant and changes inside it. However, there is no strict correlation between the log Kow and log Koc values and the fate of CECs because other factors, such as molecular volume, spatial structure, etc., also influence the processes. Ionization has also been shown to influence the plant's uptake and modifications of organic compounds because charged molecules have a reduced ability to permeate cell membranes 46. The values of the negative logarithm of the dissociation constants (pKa) of the tested compounds are 7.9, 9.6, 9.1, 10.8, and 10.7, respectively, for TRC, BPA, DES, E1, and E2 (no data for DEET) 42,43. Under the experimental conditions (pH 7), the tested compounds are in a non-ionized form, which is beneficial for phytoremediation.
3.4. Influence of the matrix on the course of phytoremediation
Experiments were carried out using raw municipal wastewater, treated municipal wastewater, and raw landfill leachate enriched with a mixture of TRC, BPA, DES, and E1 at 100 µg/L each as a culture medium, maintaining optimal conditions used in previous cultures. The study lasted seven days; the concentrations of the tested compounds were determined before the start of the experiment and after 1, 3, 5, and 7 days of contact of the plant with the solution. The CECs removal efficiency obtained using real matrices is shown in Fig. 5.
In the case of raw wastewater, the average removal efficiency of the tested compounds after seven days of treatment was 84% and 75%; in the case of treated wastewater, 93% and 89%, and in the case of landfill leachate, 59% and 56%, for W. arrhiza and L. minor, respectively. The observed effectiveness of W. arrhiza and L. minor is similar for each matrix used, with a slight advantage for W. arrhiza. A worse removal efficiency by L. minor is visible mainly in the first days of the experiment. This may indicate a more difficult acclimatization of this plant to the presence of a highly contaminated matrix. The high efficiency of CECs removal from raw wastewater indicates the resistance of the tested plants to the presence of pollutants. Raw municipal wastewater contains certain amounts of nutrients (N and P, see Table S1, Supplementary Material) that support the growth and development of plants. Additionally, there are high concentrations of easily digestible organic compounds, including carbohydrates, fats, and proteins, which may also positively affect the functioning of plants from the Lemnaceae family, given their well-documented ability to provide mixotrophic nutrition 47. The observed ability of floating plants to remove CECs from landfill leachates is significantly lower than in the case of municipal wastewater. Leachate is a liquid that is a toxic cocktail containing high concentrations of pollutants such as soluble organic matter, inorganic components, heavy metals, and xenobiotic organic compounds 48. Leachates are characterized by high phytotoxicity 49, so it can be assumed that the removal of the tested compounds during the experiment is carried out to a small extent by biological processes, mainly by sorption, hydrolysis, and photolysis.