Removal of Hexavalent Chromium in Aqueous Solution by Cellulose Filter Paper Loaded with Nano-Zero-Valent Iron: Performance Investigation and Numerical Modeling

Cr(VI) pollution in water bodies is very harmful to human health and the environment. Therefore, it is necessary to remove Cr(VI) from water. In this study, the composite (FP-nZVI) was prepared by loading nano-zero-valent iron (nZVI) onto cellulose filter paper (FP) using a liquid-phase reduction method to improve the dispersibility and oxidation resistance of nZVI. In batch experiments, the effects of iron loading of FP-nZVI, initial concentration of Cr(VI), temperature, and pH on Cr(VI) removal were particularly investigated. The maximum removal rate of 98.6% was achieved at 25 °C, pH = 5, initial concentration of Cr(VI) of 20 mg/L, and FeCl3·6H2O solution concentration of 0.8 mol/L. The removal of Cr(VI) by FP-nZVI conformed to a pseudo-second-order kinetic model and Langmuir isotherm model. The mechanism of Cr(VI) removal was a multi-step removal mechanism, involving adsorption, reduction, and coprecipitation. Column experiments investigated the effect of flow rate (1 mL/min, 3 mL/min, and 5 mL/min) on Cr(VI) removal. We found that increasing flow rate slightly decreased the removal rate of Cr(VI). The transport of Cr(VI) in composite porous media was simulated using HYDRUS-1D, and the results show that the two-site model can well simulate the reactive transport of Cr(VI). This study may provide a useful reference for the remediation of groundwater contaminated with Cr(VI) or other similar heavy metals using FP-nZVI.


Introduction
Chromium contamination in water is extremely dangerous to human health and ecosystems. Cr(VI) has greater solubility, mobility, and toxicity than Cr(III). The treatment methods for Cr(VI) in water mainly include the biological method, ion-exchange method, membrane separation technology, adsorption method, and adsorption-reduction method. The biological method [1] mainly includes the enrichment of Cr(VI) by plants and the removal of microorganisms. Guo [2] et al. used Pseudomonas to explore its adsorption capacity for Cr(VI) and found that under the optimal conditions, the removal rate of Cr(VI) by Pseudomonas could reach 72.05%. The ion-exchange method [3] is that in which Cr(VI), which mainly exists in water in the form of oxygen anion, is separated from water by replacing it with functional groups available for reaction in anion-exchange resin. Rivero [4] et al. found that it was feasible to use ion-exchange resin Lewatit MP-64 to treat Cr(VI)-polluted groundwater. Membrane separation technology [5] relies on some kind of mass transfer power (such as pressure difference, concentration difference, or potential difference) to move each component of water selectively through the membrane, so as to achieve Cr(VI) interception. Two new positively charged ultrafiltration membranes prepared by Yao [6] et al. showed high Cr(VI) adsorption capacity at low pH. Adsorption method [7] refers to the accumulation of pollutants on the adsorbent surface to achieve the purpose of Cr(VI) removal. Using activated carbon, Chen et al. found that the highest X-ray diffraction (XRD), X-ray photoelectron spectroscopy (XPS), and Brunauer-Emmett-Teller (BET). The removal effect of FP-nZVI on Cr(VI) in water was investigated by batch and column experiments. The transport of Cr(VI) in a porous medium composed of FP-nZVI and sand was simulated by HYDRUS-1D. The aims of this study are to design a filter paper composite system loaded with nZVI nanoparticles with good dispersion and antioxidant stability through a convenient preparation process and explore its performance and mechanism of Cr(VI) removal.

Preparation of FP-nZVI
In the experiment, 1.0812, 2.1624, 5.4059, 8.6495, and 10.8119 g FeCl 3 ·6H 2 O were dissolved in five parts of 40 mL ethanol solvent to prepare iron solutions with different concentrations (0.1, 0.2, 0.5, 0.8, and 1.0 mol/L). Five parts of 0.02 g NaOH were dissolved in 100 mL of deionized water, and twenty pieces of filter papers were divided into five equal portions and immersed in the solution simultaneously. After 30 min, the filter papers were removed and washed with deionized water until the deionized water was approximately neutral; the washed filter papers were immersed in the iron solution for 40 min and then taken out and air-dried to obtain yellow filter papers. Finally, five parts of 0.95 g NaBH 4 were weighed and dissolved in five parts of 100 mL deionized water to prepare five parts of the reducing agent. The five parts of yellow filter papers were immersed directly into the NaBH 4 solution for 45 min to produce black filter papers (FP-nZVI). The FP-nZVI composites were thoroughly cleaned with deionized water and anhydrous ethanol and finally dried under vacuum for 600 min in a freeze-drying oven (Shanghai Bilan Instrument Manufacturing Co., Ltd., FD-1A-50, Shanghai, China).

Characterization
The nanoscale morphology of FP-nZVI was obtained by SEM (JEOL, JEM-7001M, Tokyo, Japan) at an operating voltage of 10 kV. The surface functional groups of FP-nZVI were detected using an FTIR analyzer (Shimadzu, IRAffinity-1S, Tokyo, Japan). The crystal structures of FP-nZVI were analyzed by Bruker D8 Advance XRD instrument. The surface element composition and valence states of FP-nZVI were obtained by XPS (Thermo Fisher, ESCALAB-250Xi, Waltham, MA, USA). The BET surface areas of FP and FP-nZVI were determined by adsorption-desorption instrument (Quantachrome, Autosorb-iQ, Boynton Beach, FL, USA).

Batch Experiments
A total of four series of batch analyses were conducted to investigate the effects of different Fe 3+ concentrations used to prepare FP-nZVI, initial Cr(VI) concentration, temperature, and pH on the removal of Cr(VI). The reaction conditions are shown in Table 1. 0.8 20 25 3, 5, 7, 9, 11 Batch experiment 1 was to investigate the effect of FP-nZVI prepared with different Fe 3+ concentrations on the removal effect of Cr(VI). Five 250 mL triangular flasks were filled with 200 mL of Cr(VI) solution at a concentration of 20 mg/L, and the Cr(VI) solution was deoxygenated with N 2 for 20 min. Then, one piece of FP-nZVI was put into each triangular flask filled with Cr(VI) solution. (Five pieces of FP-nZVI used had been prepared with FeCl 3 ·6H 2 O solutions of 0.1 mol/L, 0.2 mol/L, 0.5 mol/L, 0.8 mol/L, and 1.0 mol/L, respectively.) With initial pH = 5, sealed with sealing film, and keeping the temperature at 25 • C, the triangular flasks were placed in a constant-temperature shaker (Shanghai Bilang Company, COS-100B, Shanghai, China) and shaken at 150 r/min. Samples were taken by syringe at 2,5,10,20,30,60,90,120,150,180, and 210 min. The concentrations of Cr(VI) were determined by diphenylcarbonyl dihydrazide spectrophotometric method (GB7467-87) after passing through 0.22 µm filter membrane.
Batch experiment 2 was to investigate the effect of initial concentration of Cr(VI) on the removal effect of Cr(VI). In total, 200 mL of the Cr(VI) solutions at concentrations of 5 mg/L, 10 mg/L, 20 mg/L, 30 mg/L, and 40 mg/L was added in five triangular flasks, respectively, and the Cr(VI) solution was deoxygenated with N 2 for 20 min. Then, one piece of FP-nZVI prepared with 0.8 mol/L FeCl 3 ·6H 2 O solution was put into each triangular flask. The concentration of Cr(VI) was determined after shaking in a shaker and taking samples. Initial pH, temperature, sampling scheme, and testing method were the same as batch experiment 1.
Batch experiment 3 was to investigate the effect of temperature on the removal effect of Cr(VI). Five 250 mL triangular flasks were filled with 200 mL of Cr(VI) solution at a concentration of 20 mg/L, and the Cr(VI) solution was deoxygenated with N 2 for 20 min. Then, one piece of FP-nZVI prepared with 0.8 mol/L FeCl 3 ·6H 2 O solution was put into each triangular flask. The initial pH = 5 was adjusted, and the control temperatures were kept at 15 • C, 20 • C, 25 • C, 30 • C, and 35 • C. The concentration of Cr(VI) was determined after shaking in a shaker and taking samples. The sampling scheme and testing method were the same as batch experiment 1.
Batch experiment 4 was to investigate the effect of initial pH on the removal effect of Cr(VI). Five 250 mL triangular flasks were filled with 200 mL of Cr(VI) solution at a concentration of 20 mg/L, and the Cr(VI) solution was deoxygenated with N 2 for 20 min. Then, one piece of FP-nZVI prepared with 0.8 mol/L FeCl 3 ·6H 2 O solution was put into each triangular flask. The initial pH was adjusted to 3, 5, 7, 9, and 11, and the temperature was controlled at 25 • C. The concentration of Cr(VI) was determined after shaking in a shaker and taking samples. The sampling scheme and testing method were the same as batch experiment 1.
Parallel experiments were set up for all experimental procedures, and the experimental results were averaged to reduce accidental errors.
The removal rate is calculated by Equation (1): where C 0 is the initial Cr(VI) concentration at t = 0 and C t is the remaining Cr(VI) concentration in the solution at time t. At 210 min of reaction, the removal rates of Cr(VI) by FP-nZVI prepared with 0.1, 0.2, and 0.5 mol/L iron concentration were 66.97%, 80.04%, and 98.22%, respectively. When the iron concentration was 0.8 and 1.0 mol/L, the removal rate of Cr(VI) had reached 98% at 30 min of reaction.
Since the removal effect of FP-nZVI on Cr(VI) was essentially the same for FP-nZVI prepared with iron solution concentrations of 0.8 mol/L and 1.0 mol/L, FP-nZVI prepared with iron solution concentration of 0.8 mol/L was used for subsequent experiments. These experiments included FTIR, SEM, XRD, and XPS analysis; the study of effects of initial concentration, reaction temperature, reaction pH, and aging time on the removal of Cr(VI); and the isothermal adsorption analysis, removal mechanism analysis, and column experiments.

Fe Content Measurement
First, 60 mL of 1:2 hydrochloric acid solution was measured in a triangular flask. Then, one piece of FP-nZVI was put into the triangular flask. The triangular flask was sealed with sealing film, and the temperature was kept at 25 • C. The flask was placed in a constant temperature shaker and shaken at 150 r/min for one hour to make the iron powder in FP-nZVI come off completely. Then, 1 mL of the solution was taken to be measured in a 100 mL colorimetric tube, and the iron loading of FP-nZVI was determined according to GB/T3049-2006 1,10-phenanthroline spectrophotometric method. Repeating the above procedure, the iron loading of FP-nZVI prepared by FeCl 3 ·6H 2 O solution with concentrations of 0.1 mol/L, 0.2 mol/L, 0.5 mol/L, 0.8 mol/L, and 1.0 mol/L could be determined.

Column Experiments
The column experimental device consisted of a Plexiglas cylinder (internal diameter = 3 cm, thickness = 0.3 cm, height = 23 cm), a storage container for Cr(VI) solution, a peristaltic pump, and a hose for solution delivery. The Plexiglas cylinder was filled with acid-washed quartz sand, FP-nZVI, and glass beads, with a total filling height of 18 cm, and the filling position of each material is shown in Figure 1. The quartz sand was purchased from Shanghai Macklin Biochemical Technology Co., Ltd. (Shanghai, China), with a particle size of 0.55-1 mm. The glass beads were purchased from Beijing Mreda Technology Co., Ltd. (Beijing, China), with a particle size of 3-4 mm. The peristaltic pump used in the experiment was BT100-1L-type (Baoding Lange Constant Flow Pump Co., Ltd., Baoding, China); the storage container for Cr(VI) solution was a 1L volumetric flask. The material of the hose for solution delivery was silica gel.
Each time before formal injection of Cr(VI) solution, deionized water was injected The peristaltic pump used in the experiment was BT100-1L-type (Baoding Lange Constant Flow Pump Co., Ltd., Baoding, China); the storage container for Cr(VI) solution was a 1L volumetric flask. The material of the hose for solution delivery was silica gel.
Each time before formal injection of Cr(VI) solution, deionized water was injected through the bottom of the column with a peristaltic pump to saturate the whole medium in the column and achieve a stable flow state. Then, deionized water was replaced with Cr(VI) solution. Cr(VI) solution with initial pH = 5 and initial concentration of 20 mg/L was continuously injected at flow rates of 1 mL/min, 3 mL/min, and 5 mL/min. Samples were taken at certain time intervals at the sampling port, and the concentration of Cr(VI) was determined.
In addition, the porosity and permeability coefficient of experimental medium were measured. Using the conservative substance Cl − as a tracer, transport experiments were performed to obtain the dispersion coefficient of the medium. See Text S1-S3 and Table S1-S3 for details. These parameters were used for subsequent numerical simulations of Cr(VI) in the intracolumn medium.

FTIR Analysis
FTIR analysis was performed for FP and FP-nZVI, and the results are shown in Figure 2. In Figure 2, 3394 cm −1 is the stretching vibration peak of -OH; 2914 cm −1 is the stretching vibration peak of -CH 2 ; 1645 cm −1 can be attributed to H-O-H stretching [35]; and 1044 cm −1 is the C-O stretching vibration peak [36]. When nZVI was loaded onto FP, a new peak of FP-nZVI appeared at 509 cm −1 , which can be attributed to Fe-O stretching vibration [37]. Thus, FTIR spectra indicated that nZVI was successfully loaded on the FP.

SEM Analysis
The SEM images of FP, nZVI, and FP-nZVI are shown in Figure 3. Figure 3a shows the SEM image of FP. It can be seen that FP presented a porous microstructure, and its fiber strip surface was slightly rough, which was conducive to the loading of nZVI. Moreover, the entry of nZVI into the inner space of the fiber channel can effectively avoid oxidation. Figure 3b shows a clear agglomeration of bare nZVI, which was caused by van der Waals forces and magnetic effects between the nZVI particles [38]. Figure 3c,d is the SEM images of FP-nZVI. It can be seen that the iron nanoparticles of FP-nZVI existed as separately dispersed and the particles were regular spheres with a particle size of about 100 nm. The reason was that the hydroxyl groups rich in cellulose filter paper fibers had a strong interaction with nZVI [31], which prevented the large-scale agglomeration of nZVI particles. Wang et al. also found that the presence of hydroxyl groups on the membrane surface can significantly improve iron loading [39].

SEM Analysis
The SEM images of FP, nZVI, and FP-nZVI are shown in Figure 3. Figure 3a shows the SEM image of FP. It can be seen that FP presented a porous microstructure, and its fiber strip surface was slightly rough, which was conducive to the loading of nZVI. Moreover, the entry of nZVI into the inner space of the fiber channel can effectively avoid oxidation. Figure 3b shows a clear agglomeration of bare nZVI, which was caused by van der Waals forces and magnetic effects between the nZVI particles [38]. Figure 3c,d is the SEM images of FP-nZVI. It can be seen that the iron nanoparticles of FP-nZVI existed as separately dispersed and the particles were regular spheres with a particle size of about 100 nm. The reason was that the hydroxyl groups rich in cellulose filter paper fibers had a strong interaction with nZVI [31], which prevented the large-scale agglomeration of nZVI particles. Wang et al. also found that the presence of hydroxyl groups on the membrane surface can significantly improve iron loading [39].  Figure 4 shows the XRD diffraction patterns of FP and FP-nZVI with 2θ of 15.0°, 16.6°, and 23.1°, which stand for the (110), (110), and (200) diffraction peaks of cellulose, respectively [40]. Diffraction patterns of FP-nZVI showed a 110 crystal plane diffraction peak of α-Fe at 2θ of 44.68°, which was a typical peak of nZVI [41], indicating that the prepared iron existed in the form of zero-valent iron. The diffraction peak of Fe 0 was observed at 44.8° for nZVI-BC prepared by Huang et al., which was similar to the results of this study [42]. The crystallite size of FP-nZVI was calculated to be 9.06 nm by Scherrer equation, which was smaller than the SEM result (about 100 nm). The finer sizes of the primary nZVI particles can lead to the amorphous structure of the stabilized nZVI [43,44]. The Scherrer equation obtained the crystallite size, i.e., primary particle size, while SEM obtained the particle size, i.e., secondary particle size. It is generally believed that the reason why the secondary particle size was much larger than the primary particle size was due to the agglomeration and growth of crystallites [45]. Therefore, the particle size observed by SEM was larger than the crystallite size calculated by Scherrer equation. Similar results have been found in other literature [46,47].  [40]. Diffraction patterns of FP-nZVI showed a 110 crystal plane diffraction peak of α-Fe at 2θ of 44.68 • , which was a typical peak of nZVI [41], indicating that the prepared iron existed in the form of zero-valent iron. The diffraction peak of Fe 0 was observed at 44.8 • for nZVI-BC prepared by Huang et al., which was similar to the results of this study [42]. The crystallite size of FP-nZVI was calculated to be 9.06 nm by Scherrer equation, which was smaller than the SEM result (about 100 nm). The finer sizes of the primary nZVI particles can lead to the amorphous structure of the stabilized nZVI [43,44]. The Scherrer equation obtained the crystallite size, i.e., primary particle size, while SEM obtained the particle size, i.e., secondary particle size. It is generally believed that the reason why the secondary particle size was much larger than the primary particle size was due to the agglomeration and growth of crystallites [45]. Therefore, the particle size observed by SEM was larger than the crystallite size calculated by Scherrer equation. Similar results have been found in other literature [46,47].   Figure 5a shows the XPS full-spectrum scans of FP and FP-nZVI. Binding ener near 286.08 eV and 531.2 eV were observed in both FP and FP-nZVI, which were attribu to C1s and O1s, respectively. When the filter paper was loaded with iron, one can obse the electron binding energy Fe2p at 711.08 eV and 725.1 eV in the resolution scann image, which was the photoelectron energy measured when the 2p orbital electron iron atoms were excited. In order to better explore the existing form of iron in FP-nZ narrow-spectrum scans of Fe2p were performed using XPS (as shown in Figure 5b). lines with different colors are the fitting curves of the intensity of the corresponding lence states of the iron element. The star curve is the fitting curve of the total stren curve of iron element. At the binding energies of 710.2 eV and 723.8 eV, there were 2 characteristic peaks and 2p1/2 characteristic peaks corresponding to Fe 2+ , respectively  Figure 5a shows the XPS full-spectrum scans of FP and FP-nZVI. Binding energies near 286.08 eV and 531.2 eV were observed in both FP and FP-nZVI, which were attributed to C1s and O1s, respectively. When the filter paper was loaded with iron, one can observe the electron binding energy Fe2p at 711.08 eV and 725.1 eV in the resolution scanning image, which was the photoelectron energy measured when the 2p orbital electrons of iron atoms were excited. In order to better explore the existing form of iron in FP-nZVI, narrowspectrum scans of Fe2p were performed using XPS (as shown in Figure 5b). The lines with different colors are the fitting curves of the intensity of the corresponding valence states of the iron element. The star curve is the fitting curve of the total strength curve of iron element. At the binding energies of 710.2 eV and 723.8 eV, there were 2p 3/2 characteristic peaks and 2p 1/2 characteristic peaks corresponding to Fe 2+ , respectively. At the binding energies of 711.7 eV and 725.3 eV, there were 2p 3/2 characteristic peaks and 2p 1/2 characteristic peaks corresponding to Fe 3+ , respectively. These iron oxidations occurred mainly during the synthesis process and slight oxidation of the nZVI surface [48]. The characteristic peaks corresponding to Fe 0 were found at binding energies of 707.0 eV and 720.1 eV, which proved the successful loading of nZVI particles on FP.

Int. J. Environ. Res. Public Health 2023, 20, x FOR PEER REVIEW 10 o
The characteristic peaks corresponding to Fe 0 were found at binding energies of 707.0 and 720.1 eV, which proved the successful loading of nZVI particles on FP.

BET Analysis
N2 adsorption-desorption curves of FP and FP-nZVI are shown in Figure 6. The s cific surface areas of FP and FP-nZVI were 1.517 m 2 /g and 3.963 m 2 /g, respectively. enlarged specific surface area of FP-nZVI can provide more active sites and contact a for the degradation of pollutants, thus improving the removal efficiency [49]. The aver pore size of FP-nZVI was 5.770 nm, and the pore volume was 0.010 cc/g, indicating

BET Analysis
N 2 adsorption-desorption curves of FP and FP-nZVI are shown in Figure 6. The specific surface areas of FP and FP-nZVI were 1.517 m 2 /g and 3.963 m 2 /g, respectively. The enlarged specific surface area of FP-nZVI can provide more active sites and contact area for the degradation of pollutants, thus improving the removal efficiency [49]. The average pore size of FP-nZVI was 5.770 nm, and the pore volume was 0.010 cc/g, indicating that the material was structurally porous and could also provide many reduction and adsorption sites for pollutants.

Effect of Iron Loading on the Removal of Cr(VI)
The effects of iron loading of FP-nZVI prepared with different concentratio chloride hexahydrate solution (0.1, 0.2, 0.5, 0.8, and 1 mol/L) on the removal of C investigated at 25 °C , pH = 5, and an initial concentration of Cr(VI) of 20 mg/L loadings of FP-nZVI prepared with five solutions were 55.62, 61.85, 65.28, 85.35 mg, respectively. As shown in Figure 7, the reaction can be divided into two s first 20 min was a relatively fast reaction stage. From after 20 min to the end ti basically a low-speed reaction process, and it eventually reached equilibrium. T the initial concentration of iron solution, the higher the removal rate of Cr(VI) rium. When the concentration of iron solution was 0.1, 0.2, and 0.5 mol/L, th rates of Cr(VI) reached 66.97%, 80.04%, and 98.22%, respectively, at 210 min o Furthermore, when the concentration of iron solution was increased to 0.8 m mol/L, the removal rate of Cr(VI) could reach 98% at 30 min of reaction. This w as the concentration of iron solution increased, the iron loading of FP-nZVI beca and therefore more reaction active sites were available. At the same time, it cou more contact reaction opportunities for Cr(VI), thus consuming more Cr(VI). A the reaction rate was faster and the removal rate of Cr(VI) was improved [50].  Figure 7, the reaction can be divided into two stages. The first 20 min was a relatively fast reaction stage. From after 20 min to the end time, it was basically a low-speed reaction process, and it eventually reached equilibrium. The higher the initial concentration of iron solution, the higher the removal rate of Cr(VI) at equilibrium. When the concentration of iron solution was 0.1, 0.2, and 0.5 mol/L, the removal rates of Cr(VI) reached 66.97%, 80.04%, and 98.22%, respectively, at 210 min of reaction. Furthermore, when the concentration of iron solution was increased to 0.8 mol/L and 1 mol/L, the removal rate of Cr(VI) could reach 98% at 30 min of reaction. This was because as the concentration of iron solution increased, the iron loading of FP-nZVI became larger, and therefore more reaction active sites were available. At the same time, it could provide more contact reaction opportunities for Cr(VI), thus consuming more Cr(VI). As a result, the reaction rate was faster and the removal rate of Cr(VI) was improved [50].

Effect of Initial Concentration on the Removal of Cr(VI)
Under the condition of 25 °C and pH = 5, the effects of different initial concentrations of Cr(VI) (5, 10, 20, 30, 40 mg/L) on the removal effect of Cr(VI) by FP-nZVI prepared by 0.8 mol/L ferric chloride hexahydrate solution were investigated. As shown in Figure 8, the larger the initial concentration of Cr(VI), the lower the removal rate. Fazlzadeh et al. also found that the removal rates of TV-Fe, UD-Fe, and RD-Fe prepared with three green plant extracts decreased with the increase in initial concentration of Cr(VI) [51]. This is consistent with the results of this study. When the initial concentration of Cr(VI) was 5 mg/L and 10 mg/L, the Cr(VI) removal rate reached more than 90% for the first 30 min of the reaction. When the initial concentration of Cr(VI) increased to 20 mg/L, the Cr(VI) removal rate was about 80% for the first 30 min of the reaction. When the initial concentration of Cr(VI) was further increased to 30 mg/L and 40 mg/L, the removal rates were only 59% and 51% for the first 30 min of the reaction. This was because the reactive active

Effect of Initial Concentration on the Removal of Cr(VI)
Under the condition of 25 • C and pH = 5, the effects of different initial concentrations of Cr(VI) (5, 10, 20, 30, 40 mg/L) on the removal effect of Cr(VI) by FP-nZVI prepared by 0.8 mol/L ferric chloride hexahydrate solution were investigated. As shown in Figure 8, the larger the initial concentration of Cr(VI), the lower the removal rate. Fazlzadeh et al. also found that the removal rates of TV-Fe, UD-Fe, and RD-Fe prepared with three green plant extracts decreased with the increase in initial concentration of Cr(VI) [51]. This is consistent with the results of this study. When the initial concentration of Cr(VI) was 5 mg/L and 10 mg/L, the Cr(VI) removal rate reached more than 90% for the first 30 min of the reaction. When the initial concentration of Cr(VI) increased to 20 mg/L, the Cr(VI) removal rate was about 80% for the first 30 min of the reaction. When the initial concentration of Cr(VI) was further increased to 30 mg/L and 40 mg/L, the removal rates were only 59% and 51% for the first 30 min of the reaction. This was because the reactive active sites that FP-nZVI can provide were relatively sufficient when the Cr(VI) concentration was low. With the increase in Cr(VI) concentration, more Cr(VI) would contact and react with FP-nZVI. nZVI would be rapidly oxidized, and the active sites tended to saturate, thus reducing the reaction activity. Furthermore, the high strength of Cr(VI) would reduce the electrochemical corrosion of nZVI, limit the formation of Fe(II), and reduce the reaction rate. At the same time, the reduction of Cr(VI) was incomplete, reducing the removal rate [52]. sites that FP-nZVI can provide were relatively sufficient when the Cr(VI) concentration was low. With the increase in Cr(VI) concentration, more Cr(VI) would contact and react with FP-nZVI. nZVI would be rapidly oxidized, and the active sites tended to saturate, thus reducing the reaction activity. Furthermore, the high strength of Cr(VI) would reduce the electrochemical corrosion of nZVI, limit the formation of Fe(II), and reduce the reaction rate. At the same time, the reduction of Cr(VI) was incomplete, reducing the removal rate [52].  Figure 9, it can be seen that both the reaction rate and the Cr(VI) removal  Figure 9, it can be seen that both the reaction rate and the Cr(VI) removal rate was accelerated or increased with the rising of reaction temperature. After 90 min, the removal rates of Cr(VI) at reaction temperatures of 15 • C, 20 • C, and 25 • C were 60.15%, 80.42%, and 87.23%, while the removal rates at reaction temperatures of 30 • C and 35 • C were greater than 95%. It can be seen that the increase in reaction temperature would promote the removal of Cr(VI). This was due to the reaction between nZVI and Cr(VI), which was endothermic [53]. Increasing the temperature was conducive to the positive chemical reaction, thus increasing the removal rate of Cr(VI). The increase in temperature led to a more intense movement of Cr(VI) ions in the solution, so that Cr(VI) would have a greater chance to contact and react with nZVI, thus accelerating the reaction rate. rate was accelerated or increased with the rising of reaction temperature. After 90 min, the removal rates of Cr(VI) at reaction temperatures of 15 °C , 20 °C , and 25 °C were 60.15%, 80.42%, and 87.23%, while the removal rates at reaction temperatures of 30 °C and 35 °C were greater than 95%. It can be seen that the increase in reaction temperature would promote the removal of Cr(VI). This was due to the reaction between nZVI and Cr(VI), which was endothermic [53]. Increasing the temperature was conducive to the positive chemical reaction, thus increasing the removal rate of Cr(VI). The increase in temperature led to a more intense movement of Cr(VI) ions in the solution, so that Cr(VI) would have a greater chance to contact and react with nZVI, thus accelerating the reaction rate.

Effect of Reaction pH on the Removal of Cr(VI)
Under the condition of 25 °C and initial concentration of Cr(VI) of 20 mg/L, the effects of different pH levels (3,5,7,9,11) on the removal effect of Cr(VI) by FP-nZVI prepared with 0.8 mol/L ferric chloride hexahydrate solution were investigated. As shown in Figure  10, the reaction pH had a significant impact on the removal rate of Cr(VI). The lower the pH, the faster the reaction rate and the better the removal effect. The removal mainly occurred within 2 h. When the reaction time was 60 min, the removal rate of Cr(VI) was as high as 96% when pH = 3, but the rates were only 42% and 43% when pH = 9 and 11. This showed that the pH of the solution affected the redox reaction between nZVI and Cr(VI) [54] and that acidic conditions were more favorable for the removal of Cr(VI). This was because under acidic conditions, HCrO4 − was the main form of Cr(VI) [55], and there was more H + . On the one hand, the increase in H + can accelerate the dissolution of nZVI, make nZVI release more reductive electrons, and accelerate the conversion of more Cr(VI) to the low-valence state. On the other hand, with the decrease in pH, the positive net charges of the adsorbent increase, and the electrostatic attraction between the adsorbent surface and the chromium oxyanions increases, which is conducive to the adsorption of Cr(VI) [56]. Therefore, the lower the pH was, the faster the reaction rate was and the more Cr(VI) was removed. In contrast, in alkaline environments, the nZVI surface made it easy to produce iron oxide/hydroxide precipitates, which would produce a passivating effect [10] and hinder the electron transfer of nZVI, resulting in a decrease in the reaction rate and thus inhibiting the reaction process. In addition, CrO4 2− was the main form of Cr(VI) in the alkaline environment, and the competition between high concentration of OH − and chromium oxyanions for sorption sites was intensified [57], which further

Effect of Reaction pH on the Removal of Cr(VI)
Under the condition of 25 • C and initial concentration of Cr(VI) of 20 mg/L, the effects of different pH levels (3,5,7,9,11) on the removal effect of Cr(VI) by FP-nZVI prepared with 0.8 mol/L ferric chloride hexahydrate solution were investigated. As shown in Figure 10, the reaction pH had a significant impact on the removal rate of Cr(VI). The lower the pH, the faster the reaction rate and the better the removal effect. The removal mainly occurred within 2 h. When the reaction time was 60 min, the removal rate of Cr(VI) was as high as 96% when pH = 3, but the rates were only 42% and 43% when pH = 9 and 11. This showed that the pH of the solution affected the redox reaction between nZVI and Cr(VI) [54] and that acidic conditions were more favorable for the removal of Cr(VI). This was because under acidic conditions, HCrO 4 − was the main form of Cr(VI) [55], and there was more H + . On the one hand, the increase in H + can accelerate the dissolution of nZVI, make nZVI release more reductive electrons, and accelerate the conversion of more Cr(VI) to the low-valence state. On the other hand, with the decrease in pH, the positive net charges of the adsorbent increase, and the electrostatic attraction between the adsorbent surface and the chromium oxyanions increases, which is conducive to the adsorption of Cr(VI) [56]. Therefore, the lower the pH was, the faster the reaction rate was and the more Cr(VI) was removed. In contrast, in alkaline environments, the nZVI surface made it easy to produce iron oxide/hydroxide precipitates, which would produce a passivating effect [10] and hinder the electron transfer of nZVI, resulting in a decrease in the reaction rate and thus inhibiting the reaction process. In addition, CrO 4 2− was the main form of Cr(VI) in the alkaline environment, and the competition between high concentration of OH − and chromium oxyanions for sorption sites was intensified [57], which further reduced the removal of Cr(VI). This was consistent with most literature reporting that the adsorption and reduction of Cr(VI) removal decreased with increasing pH [58,59]. reduced the removal of Cr(VI). This was consistent with most literature report adsorption and reduction of Cr(VI) removal decreased with increasing pH[58

Effect of Aging on Cr(VI) Removal
The effect of natural oxidation of FP-nZVI on the removal of Cr(VI) was in Six pieces of FP-nZVI prepared with 0.8 mol/L ferric chloride hexahydrate so placed in the same clean and open environment for different time periods (0, and 35 days) to oxidize naturally in air. Six pieces of FP-nZVI were immersed gular flasks, each filled with 200 mL of a Cr(VI) solution of 20 mg/L. The p adjusted, the flasks were sealed with sealing film, and the temperature was c 25 °C . The concentration of Cr(VI) was determined after shaking in a shaker samples. The sampling scheme and testing method were the same as batch ex As shown in Figure 11, the removal rate of Cr(VI) decreased with incre time. The main reason was that the formation of the oxide layer on the surfa hindered the electron transfer. However, the removal rate still remained abov cating that the FP-nZVI prepared in this study had good oxidation resistance in the condition of Cr(VI) concentration of 10 mg/L and pH = 5, the removal rat by CaCO3-coated nZVI prepared by Cheng et al. were 83.2% and 82.6%, r when it was aged in air for 2 days and 30 days [60], which were lower than t rate (86.3%) of Cr(VI) by FP-nZVI aged for 30 days in this study under the c Cr(VI) concentration of 20 mg/L and pH = 5. This further showed that the FP pared in this study had good oxidation resistance.

Effect of Aging on Cr(VI) Removal
The effect of natural oxidation of FP-nZVI on the removal of Cr(VI) was investigated. Six pieces of FP-nZVI prepared with 0.8 mol/L ferric chloride hexahydrate solution were placed in the same clean and open environment for different time periods (0, 7, 14, 21, 28, and 35 days) to oxidize naturally in air. Six pieces of FP-nZVI were immersed in six triangular flasks, each filled with 200 mL of a Cr(VI) solution of 20 mg/L. The pH = 5 was adjusted, the flasks were sealed with sealing film, and the temperature was controlled at 25 • C. The concentration of Cr(VI) was determined after shaking in a shaker and taking samples. The sampling scheme and testing method were the same as batch experiment 1.
As shown in Figure 11, the removal rate of Cr(VI) decreased with increasing aging time. The main reason was that the formation of the oxide layer on the surface of nZVI hindered the electron transfer. However, the removal rate still remained above 85%, indicating that the FP-nZVI prepared in this study had good oxidation resistance in air. Under the condition of Cr(VI) concentration of 10 mg/L and pH = 5, the removal rates of Cr(VI) by CaCO 3 -coated nZVI prepared by Cheng et al. were 83.2% and 82.6%, respectively, when it was aged in air for 2 days and 30 days [60], which were lower than the removal rate (86.3%) of Cr(VI) by FP-nZVI aged for 30 days in this study under the condition of Cr(VI) concentration of 20 mg/L and pH = 5. This further showed that the FP-nZVI prepared in this study had good oxidation resistance.

Kinetic Analysis
In order to understand the removal process of Cr(VI) by FP-nZVI, under the condition of 25 • C, pH=5, and Cr(VI) concentration of 20 mg/L, the kinetic experiment on Cr(VI) was carried out with FP-nZVI prepared from 0.8 mol/L ferric chloride hexahydrate solution. Pseudo-first-order and pseudo-second-order kinetic models and an intraparticle diffusion model were used to fit the experimental data. The pseudo-first-order and pseudo-second-order kinetic models and intraparticle diffusion model could be described by Equations (2)-(4): where q t (mg/g) and q e (mg/g) are the adsorption capacities at time t and equilibrium, respectively; k 1 (min −1 ) is the pseudo-first-order kinetic rate constant; k 2 (g·mg −1 ·min −1 ) is the pseudo-second-order kinetic rate constant; k 3 (g·mg −1 ·min −0.5 ) is the rate constant of the intraparticle diffusion; and c is a constant. The fitted results and the parameters are shown in Figure 12 and Table 2. The R 2 value of the pseudo-second-order kinetic model was 0.9386, which was greater than that of the pseudo-first-order kinetic model (0.8636) and the intraparticle diffusion model (0.7470). The pseudo-second-order kinetic model fitted the experimental data more accurately than the pseudo-first-order kinetic model and the intraparticle diffusion model, indicating that the adsorption process was mainly chemisorption. This was consistent with the results of Cr(VI) removal by amino-modified biochar loaded with nZVI [61].

Kinetic Analysis
In order to understand the removal process of Cr(VI) by FP-nZVI, under tion of 25 °C , pH=5, and Cr(VI) concentration of 20 mg/L, the kinetic experimen was carried out with FP-nZVI prepared from 0.8 mol/L ferric chloride hexahyd tion. Pseudo-first-order and pseudo-second-order kinetic models and an intrap fusion model were used to fit the experimental data. The pseudo-first-order an second-order kinetic models and intraparticle diffusion model could be des Equations (2)-(4): Figure 11. Removal effect of aging FP-nZVI on Cr(VI). Figure 13 shows the fitting curves of the isothermal model of FP-nZVI and Cr(VI), and the correlation parameters are listed in Table 3. The correlation coefficient R 2 of the Langmuir model was 0.9886, which was better than the coefficients 0.9635 of the Temkin model and 0.8944 of the Freundlich model. Therefore, the adsorption of Cr(VI) by FP-nZVI was more consistent with the Langmuir model, indicating that the adsorption of Cr(VI) by FP-nZVI was mainly based on monolayer chemisorption and that adsorption occurred on specific and homogeneous adsorption sites [62]. The maximum adsorption capacity obtained from the Langmuir model was estimated to be 19.7619 mg Cr g −1 FP-nZVI. Shu et al. used biocharloaded nZVI prepared from almond shell for Cr(VI) removal and found that the isothermal data were consistent with the Langmuir model [63]. The comparison of adsorption capacity with other nZVI-based adsorbents is shown in Table S4 [17,19,42,58,[64][65][66][67]. In addition, this study also investigated the adsorption of Cr(VI) by quartz sand and performed isothermal adsorption analysis to obtain parameters as the basis for the numerical simulations below. The results show that the adsorption of Cr(VI) by quartz sand was weak. The specific process of Cr(VI) adsorption study by quartz sand is shown in Text S4. An adsorption fitting diagram and related parameters are shown in Figure S1 and Table S5.   Figure 13 shows the fitting curves of the isothermal model of FP-nZVI and Cr(VI), and the correlation parameters are listed in Table 3. The correlation coefficient R 2 of the Langmuir model was 0.9886, which was better than the coefficients 0.9635 of the Temkin model and 0.8944 of the Freundlich model. Therefore, the adsorption of Cr(VI) by FP-nZVI was more consistent with the Langmuir model, indicating that the adsorption of Cr(VI) by FP-nZVI was mainly based on monolayer chemisorption and that adsorption occurred on specific and homogeneous adsorption sites [62]. The maximum adsorption capacity obtained from the Langmuir model was estimated to be 19.7619 mg Cr g −1 FP-nZVI. Shu et al. used biochar-loaded nZVI prepared from almond shell for Cr(VI) removal and found that the isothermal data were consistent with the Langmuir model [63]. The comparison of adsorption capacity with other nZVI-based adsorbents is shown in Table S4 [17,19,42,58,[64][65][66][67]. In addition, this study also investigated the adsorption of Cr(VI) by quartz sand and performed isothermal adsorption analysis to obtain parameters as the basis for the numerical simulations below. The results show that the adsorption of Cr(VI) by quartz sand was weak. The specific process of Cr(VI) adsorption study by quartz sand is shown in Text S4. An adsorption fitting diagram and related parameters are shown in Figure S1 and Table S5.

Thermodynamic Analysis
The thermodynamic parameters were calculated by Equations (5)- (7): where Q e (mg/g) is the adsorption capacity at equilibrium, C e (mg/L) is the concentration of Cr(VI) at equilibrium, K d is the adsorption equilibrium constant, T (K) is temperature, R (8.314 J/mol/K) is the gas constant, enthalpies are represented by ∆H 0 (KJ/mol), entropy change by ∆S 0 (J/mol/K), and Gibbs free energy by ∆G 0 (KJ/mol). The linear curve was fitted using 1/T versus lnK d , as shown in Figure 14. The corresponding thermodynamic parameters are shown in Table 4. The ∆H 0 of Cr(VI) adsorption by FP-nZVI was positive, demonstrating an endothermic adsorption mechanism. ∆S 0 was positive, indicating that during the adsorption process, the level of freedom in the procedures as well as the instability of the solid-liquid boundary increase [68]. ∆G 0 below zero indicated that there was a spontaneous adsorption mechanism. In addition, when the adsorption temperature went up, the absolute value of ∆G 0 grew, suggesting an increase in the degree of spontaneous Cr(VI) adsorption by FP-nZVI.

Removal Mechanism
FP-nZVI prepared with 0.8 mol/L ferric chloride hexahydrate solution was used to remove Cr(VI) at 25 °C , pH = 5, and an initial concentration of Cr(VI) of 20 mg/L. The FTIR analysis of the post-reaction material (FP-nZVI@Cr) is shown in Figure 15. It can be seen from Figure 15 that a new absorption peak appeared at 616 cm −1 , which may be attributed to the complexation effect of Fe and Cr during the coprecipitation process [37,69]. The

Removal Mechanism
FP-nZVI prepared with 0.8 mol/L ferric chloride hexahydrate solution was used to remove Cr(VI) at 25 • C, pH = 5, and an initial concentration of Cr(VI) of 20 mg/L. The FTIR analysis of the post-reaction material (FP-nZVI@Cr) is shown in Figure 15. It can be seen from Figure 15 that a new absorption peak appeared at 616 cm −1 , which may be attributed to the complexation effect of Fe and Cr during the coprecipitation process [37,69]. The surface element composition and chemical valence state of FP-nZVI@Cr were studied by XPS. The lines with different colors are the fitting curves of the intensity of the corresponding valence states of the iron element or chromium element. The star curve is the fitting curve of the total strength curve of iron element or chromium element. It can be seen from Figure 16a that after the removal reaction of Cr(VI) by FP-nZVI, the peak area of Fe 0 decreased significantly from 12.98% to 4.12% and Fe(III) increased from 33.35% to 45.31%, indicating that Fe 0 and Fe(II) were oxidized in the Cr(VI) removal process. Furthermore, as shown in Figure 16b, the binding energies of 578.59 eV and 586.74 eV are the characteristic peaks of Cr2p 3/2 and Cr2p 1/2 of Cr(VI), and the binding energies of 576.56 eV and 586.71 eV are the characteristic peaks of Cr2p 3/2 and Cr2p 1/2 of Cr(III). The peak areas of Cr(VI) and Cr(III) were estimated semi-quantitatively, indicating that 40.79% was in the Cr(VI) state, while about 59.21% was in the Cr(III) oxidation state. The existence of Cr(VI) indicated adsorption between FP-nZVI and Cr(VI), which may be due to the electrostatic interaction between FP-nZVI and Cr(VI) [57,70]. The presence of Cr(III) indicated that part of the adsorbed Cr(VI) was reduced to Cr(III). Similar results have been found in other research [71].  Therefore, the possible mechanisms for Cr(VI) removal can be inferred from the analysis results: (1) mass transfer of Cr(VI) occurred from bulk solution to the surface of FP-nZVI, and it was adsorbed; (2) then, Cr(VI) was reduced to Cr(III) by Fe 0 or Fe(II), resulting in Fe 0 being oxidized to Fe(II), and then Fe(II) continued to participate in the reduction reaction and eventually oxidized to Fe(III); (3) Fe(III) and Cr(III) were coprecipitated in the form of Fe x Cr 1−x (OH) 3 [72,73]. The removal mechanism is shown in Figure 17. The main reaction equations are shown in Equations (8)-(13) [41]: Figure 15. FTIR spectra of FP-nZVI and FP-nZVI@Cr.

Column Experiments
Using ae breakthrough curve is a typical method to assess the ability of an adsorbent to remove contaminants [74]. The experimental column was a sand column with FP-nZVI interlayer, which represented the reaction barrier. A Cr(VI) solution with an initial concentration of 20 mg/L was pumped into the sand column at three flow rates (1, 3, and 5 mL/min). The outlet Cr(VI) concentration was measured according to the time interval. The variation curves of Cr(VI) concentration in the effluent with time and PV (pore volume) were plotted, as shown in Figure 18. When the flow rates were 1, 3, and 5 mL/min, the breakthrough times were 960, 150, and 120 min, respectively. With the increase in flow rate, the breakthrough time was shortened and the curve became steeper. This was because the larger the injection velocity, the stronger the convective dispersion effect and the larger the mass of Cr(VI) outflow per unit time. At the same time, the interaction time between Cr(VI) and FP-nZVI became shorter, which led to a faster increase in Cr(VI) concentration in the effluent and a steeper curve. This was consistent with the results of previous studies [75]. However, from the perspective of PV, the start of Cr(VI) breakthrough was related to the volume of Cr(VI) solution passed through. At the three flow rates, although the start of breakthrough time was different, all curves started to breakthrough at

Column Experiments
Using ae breakthrough curve is a typical method to assess the ability of an adsorbent to remove contaminants [74]. The experimental column was a sand column with FP-nZVI interlayer, which represented the reaction barrier. A Cr(VI) solution with an initial concentration of 20 mg/L was pumped into the sand column at three flow rates (1, 3, and 5 mL/min). The outlet Cr(VI) concentration was measured according to the time interval. The variation curves of Cr(VI) concentration in the effluent with time and PV (pore volume) were plotted, as shown in Figure 18. When the flow rates were 1, 3, and 5 mL/min, the breakthrough times were 960, 150, and 120 min, respectively. With the increase in flow rate, the breakthrough time was shortened and the curve became steeper. This was because the larger the injection velocity, the stronger the convective dispersion effect and the larger the mass of Cr(VI) outflow per unit time. At the same time, the interaction time between Cr(VI) and FP-nZVI became shorter, which led to a faster increase in Cr(VI) concentration in the effluent and a steeper curve. This was consistent with the results of previous studies [75]. However, from the perspective of PV, the start of Cr(VI) breakthrough was related to the volume of Cr(VI) solution passed through. At the three flow rates, although the start of breakthrough time was different, all curves started to breakthrough at about 16 PV and reached the dynamic equilibrium state at approximately 60 PV, and the removal effect was basically the same. This indicated that the flow rate had an effect on the breakthrough time, which had a slight effect on the overall treatment effect. PV was the main factor affecting the equilibrium state and removal effect.
With an integral of breakthrough curve of C/C 0 -T or C/C 0 -PV, the capture capacity caused by reduction and adsorption could be calculated. The capture capacities at equilibrium time were 13.50, 16.46, and 14.45 mg Cr(VI) g −1 FP-nZVI at the flow rates of 1, 3, and 5 mL/min, respectively. In the column experiment, the maximum retention rate of Cr(VI) by the sand column during the period from the time of injection at the inlet to the time when the Cr(VI) concentration in the effluent was basically stable was calculated by Equations (14) and (15): Rate capture = M capture /M total in (15) where Rate capture is the maximum retention rate of Cr(VI) by the sand column; M capture (mg) is the total mass of Cr(VI) retained by the sand column; M total in (mg) is the total mass of Cr(VI) flowing in from the inlet of the sand column; and M total out (mg) is the total mass of Cr(VI) flowing out from the outlet of the sand column. about 16 PV and reached the dynamic equilibrium state at approximately 60 PV, and the removal effect was basically the same. This indicated that the flow rate had an effect on the breakthrough time, which had a slight effect on the overall treatment effect. PV was the main factor affecting the equilibrium state and removal effect. With an integral of breakthrough curve of C/C0-T or C/C0-PV, the capture capacity caused by reduction and adsorption could be calculated. The capture capacities at equilibrium time were 13.50, 16.46, and 14.45 mg Cr(VI) g −1 FP-nZVI at the flow rates of 1, 3 and 5 mL/min, respectively. In the column experiment, the maximum retention rate of Cr(VI) by the sand column during the period from the time of injection at the inlet to the time when the Cr(VI) concentration in the effluent was basically stable was calculated by Equations (14) and (15): The overall maximum retention rates of Cr(VI) in the sand column were 57.91%, 53.15%, and 50.39% at the flow rates of 1, 3, and 5 mL/min, respectively. The retention rate of Cr(VI) by the column decreased slightly with the increase in flow rate. Digiacomo et al. used sulphate green rust and S-nZVI to remove Cr(VI) in the column and found that the flow rate had a lower impact on the removal rate, but the higher the flow rate, the lower the removal rate [76]. This was consistent with the results of this study. The main reason for the decrease in the retention rate was that the larger the flow rate, the shorter the contact reaction time between Cr(VI) and FP-nZVI, and the reaction was not sufficient. With the injection of Cr(VI), the remaining active sites of FP-nZVI continued to decrease, and the Cr(VI) consumed by the reaction gradually decreased. Therefore, the Cr(VI) concentration in the effluent gradually increased.

Reactive Transport Model of Cr(VI) in Column
To further verify the reactive transport of Cr(VI), the breakthrough procedure was modeled using the CDE model in HYDRUS-1D. The water flow model was described by the modified Richards equations [77], given as follows in Equations (16) and (17): where h (cm) is the pressure head, θ is the volume water content, t (min) is the time, and S (min −1 ) is the source sink term. α is the angle between the flow direction and the vertical direction, and in this study, α = 0 • . x (cm) is the vertical coordinate, and upward direction is defined as positive. K (cm/min) is the unsaturated hydraulic conductivity, K r (cm/min) is the relative hydraulic conductivity, and K s (cm/min) is the saturated hydraulic conductivity. The water flow was stable, and the composite porous media was saturated. The lower boundary was assumed as the constant flow boundary, and the upper boundary was the given head boundary. The water quality transport model considered two-site adsorption.
The two-site model is a chemical non-equilibrium model, in which exchange sites are divided into two types. Type 1 sites (S 1 ) are in equilibrium with the solution phase, and on type 2 sites (S 2 ), sorption is considered to be time-dependent. In this study, sorption is described by the Langmuir equation, and the two-site model is given as Equation (18) [78]: where f is the equilibrium site fraction, α (min −1 ) is the first-order kinetic rate coefficient, K L (mL/g) is the Langmuir adsorption coefficient, and b (mg/kg) is the maximum adsorption capacity of the medium to solute. S 2 is the solid phase concentration at the type 2 site, and θ is the volumetric water content. c (mg/mL) is the solute concentration in the liquid phase, D (cm 2 /min) is the dispersion coefficient, v (cm/min) is pore water velocity, x (cm) is the spatial position, ρ (g/cm −3 ) is the bulk density of porous media, and t is time.
HYDRUS-1D requires input parameters such as q q = Q A , θ, λ λ = D v , b, K L , f, and α to build a nonlinear non-equilibrium model (NLNE). In the water quality transport model, the solute concentration at each position of the column was assumed as 0 at the initial moment. The lower boundary was the given concentration boundary, and the upper boundary was assumed as the flux boundary.
It had been found that sorption parameters obtained from static batch experiments may not be suitable for determining solute transport in column experiments [79]. Because the interaction between FP-nZVI and Cr(VI) under static conditions was confined to a relatively stable experimental environment, however, the soil solution system composed of soil particles and pore solutions under dynamic conditions would have inhomogeneity, which could affect the chemical reaction process of heavy metal ions in soil. Therefore, the reaction state and parameters expressed in static experiments cannot fully reflect the mass transfer processes such as sorption, reaction conditions, and transport state in the soil [80]. K L and b obtained from static batch experiments can be substituted into the simulation model as initial parameters. The calculated values were fitted with the observed breakthrough data and curves of the column experiments for iterative adjustment, and thus, the parameters K L and b were trimmed [81]. The kinetic parameters (f and α) were obtained automatically by the inverse solution module of HYDRUS-1D, and the model parameters are shown in Table 5. The decrease in f indicated that more equilibrium adsorption was transformed into non-equilibrium adsorption. α described the adsorption/desorption ratio of Cr(VI) in time. The comparison between the observed breakthrough curves and calculated values of Cr(VI) at three flow rates is shown in Figure 19. The overall changing trend of the calculated and observed values was consistent. The root mean square error (RMSE) ranged from 0.489% to 1.204%, and R 2 values were all greater than 0.99. This indicated that the two-site model of HYDRUS-1D can simulate well the transport of Cr(VI) in composite porous media. Montalvo et al. used the two-site model to fit the reactive transport of arsenate in porous media, and its kinetic parameters were also obtained by inversion. The fitted breakthrough curve could match the observed data [78]. The model was helpful to understand the transport mechanism of Cr(VI) in porous media. In this study, both K L and b after adjustment were smaller than those of batch experiments. This was because, in the static experimental environment, FP-nZVI can provide more reaction sites, and the reaction time of FP-nZVI and Cr(VI) was more sufficient, but column experimental conditions did not have these characteristics. At the end of the column experiment, C/C 0 = 0.8, and there still existed a retaining effect on nZVI in the column, which had not reacted completely and not reached the adsorption capacity of FP-nZVI. This was one of the reasons why the adsorption parameters of the column experiment were smaller than those of the static batch experiment. In addition, the K L and b obtained from the inversion can provide some basis for the prediction of the transport of Cr(VI) in groundwater.

Conclusions
Batch experiments showed that the removal rate of Cr(VI) increased with the increase in iron loading and reaction temperature, and acidic conditions were more conducive to

Conclusions
Batch experiments showed that the removal rate of Cr(VI) increased with the increase in iron loading and reaction temperature, and acidic conditions were more conducive to the removal of Cr(VI) by FP-nZVI. The removal rate was higher than 86% at 35 days of aging, indicating that FP-nZVI had very good oxidation resistance. The removal of Cr(VI) by FP-nZVI conformed to the pseudo-second-order kinetic model, and FP-nZVI had a good adsorption capacity for Cr(VI), with the maximum adsorption capacity of b = 19.7619 mg Cr g −1 by FP-nZVI. The removal process of Cr(VI) by FP-nZVI included adsorption, reduction, and coprecipitation. Column experiments discovered that increasing flow rate slightly decreased the removal rate of Cr(VI) and pore volume ratio of fluid passing through was the main factor affecting the removal rate of Cr(VI). HYDRUS-1D simulation results show that the two-site model can simulate well the transport of Cr(VI) in composite porous media. The inversion parameters provided a basis for the transport of Cr(VI) in groundwater. This study demonstrated that FP-nZVI had a strong ability for removal of Cr(VI) and has potential practical application value in the remediation of heavy-metal-contaminated groundwater.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/ijerph20031867/s1. Text S1: Porosity experiment; Text S2: Permeability coefficient experiment; Text S3: Dispersion coefficient experiment; Text S4: Adsorption Experiment of Cr(VI) by quartz sand; Figure S1: Adsorption Fitting Curve of Cr(VI) by quartz sand; Table S1: Porosity-related parameters and results; Table S2: Permeability-coefficient-related parameters and results; Table S3: Dispersion-coefficient-related parameters and results; Table S4: Comparison of removal capacity for Cr(VI) by FP-nZVI with other adsorbents; Table S5. Isothermal model parameters of quartz sand to Cr(VI), standard error was in parentheses.