A New Polyvinylidene Fluoride Membrane Synthesized by Integrating of Powdered Activated Carbon for Treatment of Stabilized Leachate

Stabilized landfill leachate contains a wide variety of highly concentrated non-biodegradable organics, which are extremely toxic to the environment. Though numerous techniques have been developed for leachate treatment, advanced membrane filtration is one of the most environmentally friendly methods to purify wastewater effectively. In the current study, a novel polymeric membrane was produced by integrating powdered activated carbon (PAC) on polyvinylidene fluoride (PVDF) to synthesize a thin membrane using the phase inversion method. The membrane design was optimized using response surface methodology (RSM). The fabricated membrane was effectively applied for the filtration of stabilized leachate using a cross-flow ring (CFR) test. The findings suggested that the filtration properties of fabricated membrane were effectively enhanced through the incorporation of PAC. The optimum removal efficiencies by the fabricated membrane (14.9 wt.% PVDF, 1.0 wt.% PAC) were 35.34, 48.71, and 22.00% for COD, colour and NH3-N, respectively. Water flux and transmembrane pressure were also enhanced by the incorporated PAC and recorded 61.0 L/m2·h and 0.67 bar, respectively, under the conditions of the optimum removal efficiency. Moreover, the performance of fabricated membranes in terms of pollutant removal, pure water permeation, and different morphological characteristics were systematically analyzed. Despite the limited achievement, which might be improved by the addition of a hydrophilic additive, the study offers an efficient way to fabricate PVDF-PAC membrane and to optimize its treatability through the RSM tool.


Introduction
Sanitary landfills are the widely applied technique to tackle municipal solid waste (MSW). Inappropriately, the majority of these landfills do not fulfill the normal discharged limits [1]. In developing countries such as Malaysia, more than 80% of the MSW produced was received by open duping and landfill sites [2]. This resulted in the generation of highly contaminated leachate, which is the liquid generated due to the precipitation above these solid litters and could be toxic to the surrounding environment. This leachate could contaminate the sources of fresh water if not carefully treated before discharging to the environment [3]. Stabilized leachate, which is more than ten years old, has lower BOD 5 /COD ratio. Thus, it is almost impossible to treat this kind of leachate using some

Collection of Leachate
Leachate sample was taken form Sahom landfill site located in Perak, Malaysia, which is an operative landfill site with a daily production of 100 tonnes of MSW in average [31]. After collection of leachate sample, it was stored in a refrigerator at 4 • C. Initial leachate characterization was performed using standardized methods of water and wastewater [32]. All measurements, including dissolved oxygen (DO), colour, chemical oxygen demand (COD), 5-day biochemical oxygen demand (BOD 5 ), and ammoniacal nitrogen (NH 3 -N), were undertaken in triplicate.

Materials
The PVDF polymer (Kynar ® 740) was purchased from Afza Maju trading (Terengganu, Malaysia), and utilized after drying for 24 h at 70 • C. 1-Methyl-2-pyrrolidone (NMP, 99.5%) was purchased from Sigma-Aldrich. Methanol, (99.8%) was supplied by Chem Soln. Ultrapure distilled water (DI) was utilized throughout the experiments. PAC was purchased from R&M Chemicals. The AC was charcoal-based, and consists of sulfide, chloride, calcium, sulphate, iron, lead, zinc, and copper. This PAC density was 1.8-2.1 kg/m 3 with pH (4)(5)(6)(7). Particle size analysis (PSA) and field emission scanning electron microscopy (FESEM) tests were used to investigate the distribution and the size of PAC particles, respectively. All these chemical materials were of analytical grade, and used without additional treatment.

Experiment's Design and Optimization Process
Central Composite Design (CCD) is the design method used in response surface methodology (RSM) for the membrane fabrication's experimental design [33]. Both CCD and RSM were run by version 8 from the Design Expert. For membrane dope solution design, two factors, the polymer (PVDF) weightage and the additive (PAC) weightage, were set into the CCD. Based on preliminary experiments and the extensive literature [34,35], the total mass of fabricated membrane dope was fixed at 100 g, which represents 100% of the dope weight, thus each 1 g of the dope element is equivalent to 1% weightage. The dosage of the PVDF was set within the range of 10 g to 18 g, and the amount of PAC was set within the range of 0 g to 2 g. Regarding the CCD, the alpha value was selected to be 1.0, and thus the centre points were 14.0 and 1.0 wt.% for the polymer content and additive content, respectively. The rest of the dope weight (to complete 100 g) is the NMP solvent. The total concentration of PVDF/PAC was kept at 20% (as maximum) and 10% (as minimum), as concentration higher than 20% resulted in solutions of extremely high viscosity, and was difficult to be casted on the glass plate, while clumsy, non-thick membrane was the result of using concentration less than 10%. Five responses, which are the removing efficiencies of COD, colour, and NH 3 -N, as well as maximum transmembrane pressure (max. TMP), and pure water flux, were also set into the CCD to have the full design of experiments. The influence of various parameters was optimized by RSM using a combination of statistical and numerical techniques. In the current work, nine experiments were reinforced with four replications to assess the pure error [36]. The 13 different membranes were applied in double repetition and have their effluent collected. The quadratic model for every response was investigated by analysis of variances (ANOVA) to identify the results significancy, and to find the represented quadratic model after eliminating irrelevant terms. The frontal sign of each model term signifies to either antagonistic or synergistic effect on the response when it is positive or negative, respectively [4]. In RSM, it does mention that Prob > F less than 0.050 indicates model terms are significant, and Prob > F with the values greater than 0.10 indicates model term is not significant. "Not significant", in the description of lack of fit, is regarded a decent model, as it means the experimental reading is fitting the model [37]. Additionally, a good experimentally fitted data will have a higher coefficient (R 2 ) value. The higher the R 2 value, the closer the experimental data towards the predicted graph model by the RSM [38,39]. Selection of the best membrane takes into consideration the membrane purification performance. Desirability value closer to 1.0 used to be selected as the ideal design for the data.

PVDF-PAC Membrane Fabrication 2.4.1. Dope Preparation
To produce the polymeric membrane, PVDF and NMP were applied as polymer and solvent, respectively. Figure 1 presents the process used for the dope preparation. Initially, the polymeric PVDF was entirely dissolved in the NMP solvent at a temperature ranged between 60 and 70 • C using a heating mantle (Figure 1a). In order to achieve a better permeate flux of the synthesized membrane, the heating mantle temperature should always be maintained within the above stated range [21]. The dope solution containing dissolved PVDF polymer in the NMP solvent was then infused into a clean Schott bottle. After that, the required amount of PAC was inserted into the dope solution to generate the dope for hybrid membrane. Lastly, the Schott bottle containing the dope solution was placed into a sonicator bath (Cole-Parmer, Vernon Hills, IL, USA) for eight hours to confirm the homogeneous mixing of the additives without any air bubbles raised in the prepared dope [40].
significant", in the description of lack of fit, is regarded a decent model, as it means the experimental reading is fitting the model [37]. Additionally, a good experimentally fitted data will have a higher coefficient (R 2 ) value. The higher the R 2 value, the closer the experimental data towards the predicted graph model by the RSM [38,39]. Selection of the best membrane takes into consideration the membrane purification performance. Desirability value closer to 1.0 used to be selected as the ideal design for the data.

Dope Preparation
To produce the polymeric membrane, PVDF and NMP were applied as polymer and solvent, respectively. Figure 1 presents the process used for the dope preparation. Initially, the polymeric PVDF was entirely dissolved in the NMP solvent at a temperature ranged between 60 and 70 °C using a heating mantle (Figure 1a). In order to achieve a better permeate flux of the synthesized membrane, the heating mantle temperature should always be maintained within the above stated range [21]. The dope solution containing dissolved PVDF polymer in the NMP solvent was then infused into a clean Schott bottle. After that, the required amount of PAC was inserted into the dope solution to generate the dope for hybrid membrane. Lastly, the Schott bottle containing the dope solution was placed into a sonicator bath (Cole-Parmer, Vernon Hills, IL, USA) for eight hours to confirm the homogeneous mixing of the additives without any air bubbles raised in the prepared dope [40].

Membrane Casting
A semi-automated membrane casting machine (TECH INC, Chennai, India) was applied to synthesize a flat sheet membrane using the dry-wet phase process, as illustrated in Figure 2a. The membrane was produced at temperature 27 °C to 30 °C with an approximate thickness of 60 µm based on literature reports [41,42]. After 60 s of membrane casting above the glass board, it was submerged into a distilled water (DW) basin for 180 s (Figure 2b). As a result, a thin layered polymeric film was generated, which separated from the glass plate. Later, the newly produced membrane was transferred into a DW coagulation bath and remained there for 24 h. Afterwards, a methanol bath was used for 8 h, as shown in Figure 2c, to perform a post-treatment to ensure the excess solvent in the membrane can be removed completely [43]. Finally, the membrane was dried 24 h at the ambient temperature with 60% humidity, as shown in Figure 2d, to be ready to use in the filtration process [13].

Membrane Casting
A semi-automated membrane casting machine (TECH INC, Chennai, India) was applied to synthesize a flat sheet membrane using the dry-wet phase process, as illustrated in Figure 2a. The membrane was produced at temperature 27 • C to 30 • C with an approximate thickness of 60 µm based on literature reports [41,42]. After 60 s of membrane casting above the glass board, it was submerged into a distilled water (DW) basin for 180 s (Figure 2b). As a result, a thin layered polymeric film was generated, which separated from the glass plate. Later, the newly produced membrane was transferred into a DW coagulation bath and remained there for 24 h. Afterwards, a methanol bath was used for 8 h, as shown in Figure 2c, to perform a post-treatment to ensure the excess solvent in the membrane can be removed completely [43]. Finally, the membrane was dried 24 h at the ambient temperature with 60% humidity, as shown in Figure 2d, to be ready to use in the filtration process [13].

Membrane Performance and Characterization
The produced membranes have been characterized to investigate their treatmen ficiencies, fouling, and permeability properties and surface morphologies. To ensure accuracy of the findings, all of the tests have been duplicated. Each time, a fresh membr has been utilized to investigate their characteristics and performance.

Treatment Efficiency
The membrane filtration performance was investigated using laboratory scale cr flow filtration setup with a 3.34 cm disc diameter, as exhibited in Figure 3. The membr rejection capabilities were studied against the treatment of landfill leachate. Before e experiment, initial characterization of leachate was measured to eliminate the small er which occurred due to the minor changes in organics concentration with time. The ste flux for all individual membranes was acquired by a constant (200 mL/min) flow for min. The volume of permeate, along with the recorded transmembrane pressure, w noted down under the flow of 200 mL/min for different intervals of time (0.5, 1, 2, 5 20, 40, 90, and 120 min).

Membrane Performance and Characterization
The produced membranes have been characterized to investigate their treatment efficiencies, fouling, and permeability properties and surface morphologies. To ensure the accuracy of the findings, all of the tests have been duplicated. Each time, a fresh membrane has been utilized to investigate their characteristics and performance.

Treatment Efficiency
The membrane filtration performance was investigated using laboratory scale crossflow filtration setup with a 3.34 cm disc diameter, as exhibited in Figure 3. The membrane rejection capabilities were studied against the treatment of landfill leachate. Before each experiment, initial characterization of leachate was measured to eliminate the small errors which occurred due to the minor changes in organics concentration with time. The steady flux for all individual membranes was acquired by a constant (200 mL/min) flow for 120 min. The volume of permeate, along with the recorded transmembrane pressure, were noted down under the flow of 200 mL/min for different intervals of time (0.5, 1, 2, 5, 10, 20, 40, 90, and 120 min).  Final leachate characterizations were evaluated in terms of removing efficie the COD, colour, and NH3-N pollutants using Equation (1): where CF is the contaminant concentration at the feed (mg/L) and CP is the cont concentrations in the permeated solution (mg/L). All contaminants' concentrati checked using the UV-V spectrophotometer (Hach DR6000, Loveland, Colorado prior and post of filtration practice.

Productivity of Membrane
Pure flux plays a dynamic role in the membrane productivity evaluation. Pe ity of membrane was investigated through the pure water flux, which was meas a dead-end filtration apparatus, as illustrated in Figure 4. A metallic ring havi average pore size and 8.76 cm 2 effective permeate area was applied to support t brane. Initially, the impurities present in the membrane were removed by subme membrane in DW for 30 min. Then, a stable flux was achieved by pre-compa membrane with N2 gas at a pressure of 30 KPa for 2 min. After 30 min, the pe water volume was noted at a similar pressure of 30 KPa. The pure water flux ca culated using the Equation (2): where V is the permeated pure water volume (L), A is the membrane effective sur (m 2 ), and t is the time of permeation (h). Final leachate characterizations were evaluated in terms of removing efficiencies for the COD, colour, and NH 3 -N pollutants using Equation (1): where C F is the contaminant concentration at the feed (mg/L) and C P is the contaminants concentrations in the permeated solution (mg/L). All contaminants' concentrations were checked using the UV-V spectrophotometer (Hach DR6000, Loveland, CO, USA) in prior and post of filtration practice.

Productivity of Membrane
Pure flux plays a dynamic role in the membrane productivity evaluation. Permeability of membrane was investigated through the pure water flux, which was measured via a deadend filtration apparatus, as illustrated in Figure 4. A metallic ring having 5 mm average pore size and 8.76 cm 2 effective permeate area was applied to support the membrane. Initially, the impurities present in the membrane were removed by submerging the membrane in DW for 30 min. Then, a stable flux was achieved by pre-compacting the membrane with N 2 gas at a pressure of 30 KPa for 2 min. After 30 min, the permeated water volume was noted at a similar pressure of 30 KPa. The pure water flux can be calculated using the Equation (2): where V is the permeated pure water volume (L), A is the membrane effective surface area (m 2 ), and t is the time of permeation (h).

Antifouling Valuation
Throughout the membrane filtration process, the overall decrease in flux, alongside the improvement of transmembrane pressure, were mainly caused by either membrane fouling, concentration polarization, or a combination of both [44]. Both of these components can be attained from the experimental data using both of the leachate permeate flux and maximum transmembrane pressure (Max. TMP) values which are measured by the cross-flow ring test. Max. TMP was applied to indicate the antifouling ability of fabricated membranes [45].

Morphological Characteristics
It is a well-known fact that the membrane properties and performance are highly dependent on its morphology (pore size, surface texture, and microstructure). Therefore, investigation of membrane morphologies is considered a significant factor in the effectiveness evaluation of the produced membranes.
EDX is a chemical microanalysis method used for quantitative, qualitative, and elemental mapping examination. Octane Silicon Drift Detector (SDD, EADX Inc., Mahwah, NJ, USA) was used at high voltage of 15 kV, using Mn Kα as source of energy. The fabricated PVDF-PAC membranes with different compositions were measured by INCA Energy 400 software (Firmware INCA, Version V1.09R13), along with the image taken by the Quanta FEG 450 instrument.
FESEM (Quanta FEG 450, FEI, Hillsboro, USA) was applied to record the cross-sectional and surface morphologies of the fabricated membrane. The cross-sectional morphologies were investigated by fracturing the membranes in liquid nitrogen and immediately cutting them after air drying. FESEM measurement starts by placing the sample on carbon tape, which was attached with the sample stub. The sample was also coated with the platinum nanoparticles in auto fine coater (JFC-1600, SUTD-MIT International Design Centre, Singapore) before performing the analysis.

Antifouling Valuation
Throughout the membrane filtration process, the overall decrease in flux, alongside the improvement of transmembrane pressure, were mainly caused by either membrane fouling, concentration polarization, or a combination of both [44]. Both of these components can be attained from the experimental data using both of the leachate permeate flux and maximum transmembrane pressure (Max. TMP) values which are measured by the cross-flow ring test. Max. TMP was applied to indicate the antifouling ability of fabricated membranes [45].

Morphological Characteristics
It is a well-known fact that the membrane properties and performance are highly dependent on its morphology (pore size, surface texture, and microstructure). Therefore, investigation of membrane morphologies is considered a significant factor in the effectiveness evaluation of the produced membranes.
EDX is a chemical microanalysis method used for quantitative, qualitative, and elemental mapping examination. Octane Silicon Drift Detector (SDD, EADX Inc., Mahwah, NJ, USA) was used at high voltage of 15 kV, using Mn Kα as source of energy. The fabricated PVDF-PAC membranes with different compositions were measured by INCA Energy 400 software (Firmware INCA, Version V1.09R13), along with the image taken by the Quanta FEG 450 instrument.
FESEM (Quanta FEG 450, FEI, Hillsboro, OR, USA) was applied to record the crosssectional and surface morphologies of the fabricated membrane. The cross-sectional morphologies were investigated by fracturing the membranes in liquid nitrogen and immediately cutting them after air drying. FESEM measurement starts by placing the sample on carbon tape, which was attached with the sample stub. The sample was also coated with the platinum nanoparticles in auto fine coater (JFC-1600, SUTD-MIT International Design Centre, Singapore) before performing the analysis.
An atomic force microscopy (AFM, Dimension 5000, Bruker AXS, Santa Barbara, CA, USA) was also applied to study the surface morphologies and roughness of the synthesized membranes. Herein, membranes were cut into small square pieces (1 × 1 cm) and pasted on a glass slide. Sample scanning were performed using a probe-optical microscope on tapping mode and images of 10 µm × 10 µm were taken by AFM. The root-mean-square roughness (R q ) and average roughness (R a ) was applied to measure the surface roughness for each membrane.
Porosity of membrane could be easily defined as the pore's volume divided by the membrane total volume. Wet membranes were weighed after carefully wiping the surface (Ww). Afterwards, these membranes were dried in an oven at 50 • C for 24 h and weighed again (Wd). The porosity of membrane ε (%), was measured by gravimetric method using where, Ww is the weight of wet membrane (kg), Wd represents the weight of dry membrane (kg), ρw is the density of water (1000 kg/m 3 ), and ρp, the polymer density (1770 kg/m 3 for PVDF). Based on the measured distilled water flux, the average pore size (d) of the membrane was calculated by the Guerout-Elford-Ferry equation, Equation (4) [46].
Herein, ε is membrane porosity (%), δ, the water viscosity (8.9 × 10 −4 Pa s), l represents membrane thickness (60 × 10 −6 m), V is the volume of the distilled water penetrating through the membrane (m 3 ), t is the experimental time interval (s), A, the effective membrane surface area (m 2 ), and ∆P is the working pressure (30 kPa). Table 1 displays the key characteristics of the raw leachate sample of more than 10 years in age. The lower BOD 5 to COD ratio (0.074) was another strong indication of highly stabilized leachate sample [3]. The other quality parameters of leachate, such as COD, BOD 5 , NH 3 -N, colour, and pH values, were around 1188 mg/L, 89 mg/L, 313 mg/L, 1360 PtCo/L, and 8.33, respectively. These obtained values were also compared with the standard discharged limits set by the Malaysian Environmental Quality was conducted (Table 1) [47]. As shown in Table 1, the COD, colour, and NH 3 -N concentrations were found to be far greater than the standard discharged limits.

PAC Characterization
Analysis test of the particle size was conducted to investigate the particle size distribution of fine samples in terms of volume. The particle size distribution of PAC sample is shown in Figure 5a. It can be seen from Figure 5 that PAC has small particle sizes which varied between (0.02-50 µm) in diameter. The average particle diameter of the PAC is 25 µm. It is evident from Figure 5a that the distribution curve of PAC particles could be counted a uniform-distribution curve. The percentage of adsorption is higher for those adsorbents have smaller particle size due to the availability of more surface area [48]. The surface morphology of PAC was visualized via FESEM, with a magnification of 10,000×, as shown in Figure 5b. FESEM micrographs of PAC, shows uniform size particles, which confirmed the results obtained from the particle size analysis. To some extent, the PAC surface having small cavities, pores, and more rough surfaces indicates the presence of an interconnected porous network. Increasing the particles' number of an adsorbent material by decreasing its particles size resulted in increasing the adsorption surface area, and thus the material adsorption characteristics [49].
Water 2021, 13, x FOR PEER REVIEW 9 of 26

PAC Characterization
Analysis test of the particle size was conducted to investigate the particle size distribution of fine samples in terms of volume. The particle size distribution of PAC sample is shown in Figure 5a. It can be seen from Figure 5 that PAC has small particle sizes which varied between (0.02-50 µm) in diameter. The average particle diameter of the PAC is 25 µm. It is evident from Figure 5a that the distribution curve of PAC particles could be counted a uniform-distribution curve. The percentage of adsorption is higher for those adsorbents have smaller particle size due to the availability of more surface area [48]. The surface morphology of PAC was visualized via FESEM, with a magnification of 10,000×, as shown in Figure 5b. FESEM micrographs of PAC, shows uniform size particles, which confirmed the results obtained from the particle size analysis. To some extent, the PAC surface having small cavities, pores, and more rough surfaces indicates the presence of an interconnected porous network. Increasing the particles' number of an adsorbent material by decreasing its particles size resulted in increasing the adsorption surface area, and thus the material adsorption characteristics [49].

Membrane Filtration and Experimental Results
Herein, the relationship among the independent factors (PVFD and PAC dosage in membrane) and responses (COD, NH3-N, colour removal, max. TEM, and pure water flux) were thoroughly investigated. There were 13 different experiments performed on the PVDF and PAC composition based on the central RSM composite design, as shown in Table 2. CFR test was performed to investigate the pollutants removal efficiency together with the max. TEM, while dead-end test was executed to measure the pure water flux.
The COD, colour, and NH3-N removal efficiencies were found to be around 14.8-37.2, 14.6-56.3, and 7.5-23.8%, respectively, while the pure flux and max. TMP were ranged between 26.2-127.7 L/m 2 ·h and 0.42-1.00 bar, respectively. ANOVA analysis was performed for the further investigation on the obtained experimental results.
It is observed from Table 2 that an increase in both PVDF and PAC concentrations on the membrane leads, to some extent, to an increase in the contaminants removal. When PVDF and PAC concentration are higher than 14 wt.% and 1.0 wt.%, respectively, the removal efficiency starts to decrease with increasing the amount of PVDF and PAC. This behaviour was attributed to the combination effect between polymer and additive in dope. This leads to the creation of large volume voids with increasing polymer dosage, and allows the small particles of contaminants to pass through the membrane [50].

Membrane Filtration and Experimental Results
Herein, the relationship among the independent factors (PVFD and PAC dosage in membrane) and responses (COD, NH 3 -N, colour removal, max. TEM, and pure water flux) were thoroughly investigated. There were 13 different experiments performed on the PVDF and PAC composition based on the central RSM composite design, as shown in Table 2. CFR test was performed to investigate the pollutants removal efficiency together with the max. TEM, while dead-end test was executed to measure the pure water flux. The COD, colour, and NH 3 -N removal efficiencies were found to be around 14.8-37.2, 14.6-56.3, and 7.5-23.8%, respectively, while the pure flux and max. TMP were ranged between 26.2-127.7 L/m 2 ·h and 0.42-1.00 bar, respectively. ANOVA analysis was performed for the further investigation on the obtained experimental results.
It is observed from Table 2 that an increase in both PVDF and PAC concentrations on the membrane leads, to some extent, to an increase in the contaminants removal. When PVDF and PAC concentration are higher than 14 wt.% and 1.0 wt.%, respectively, the removal efficiency starts to decrease with increasing the amount of PVDF and PAC. This behaviour was attributed to the combination effect between polymer and additive in dope. This leads to the creation of large volume voids with increasing polymer dosage, and allows the small particles of contaminants to pass through the membrane [50]. Table 3 depicts the empirical model using the data obtained from COD, colour, and NH 3 -N removals. F-values of the model, together with the low probability values (P > F > 0.05), clearly suggest that the models were significant for all responses. The significant model terms for COD removals in the ANOVA analysis were sorted in descending order depending upon the influential terms (AB, B 2 , A, B, and A 2 ). It was clearly seen that the PVDF and PAC (AB) had the highest impact on the COD removal with an F-value of around 32.25, followed by the quadratic term of PAC concentration (B 2 ), PVDF concentration (A), PAC concentration (B), and finally the quadratic term of PVDF concentration (A 2 ) with an F-value of 12.18, 4.34, 4.19, and 0.42, respectively. The quadratic terms of PVDF concentration together with the linear terms of PAC and PVDF contents caused a positive effect on the COD removal. Nonetheless, interaction and quadratic terms of PAC exhibited negative effects. In fact, an increase in the COD removal was recorded upon the change in the liner terms of PVDF and PAC concentrations, and PVDF concentration with quadratic term from lower to higher level. Hence, this change is complemented by the outstanding COD removal using PVDF-PAC membrane. On the other hand, a decline in COD removal was recorded when the interaction term and quadratic term of PAC was in the higher level.

Removal Efficiency of Contaminants
The quadratic term of PVDF contents (A 2 ) has the most significant effect towards the colour removal rate. This is due to the highest F-value (97.03), where other terms had the values of 21.10, 5.69, 3.89, respectively. The PAC content (B) had a progressive influence on the colour removing. However, the quadratic term of PVDF, PVDF concentration, and interaction among the PVDF and PAC displayed a negative effect. Thus, the removal of colour was enhanced with the enhancement of the PVDF contents in membrane fabrication until the optimum amount (>14 wt.% PVDF).
Additionally, in case of NH 3 -N removal, the A, B, B 2 , A 2 , and AB were sorted in descending order of their effecting strength. The highest F-value of 55.81 was recorded for the linear term of PVDF concentration (A), and thus it had the huge effect in NH 3 -N removal. On the other hand, the lowest F-value of 0.032 was recorded for interaction term, which regarded to have a negligible effect on the model. The PVDF linear term only offered a strong influence on removing of NH 3 -N, while the remaining terms were found to be the negligible influencers. Hence, the NH 3 -N removal was increased upon enhancing the PVDF contents in membranes. However, for PAC concentration after the point (PAC = 1.0 wt.%); when either the quadratic term of PVDF or PAC, or the interaction term is in the significant level, the NH 3 -N removal starts to decrease.
The lack of fit F-statistic was statistically not significant, as the values of (P) were higher than 0.05. A significant lack of fit suggests that there may be some systematic variation unaccounted for the proposed models. This may be due to the exact replicate values of the independent variables in the models that provide an estimate of pure error [15]. The correlation coefficient value (R 2 ) resulted in the present study for COD removal (0.9482), colour removal (0.9411), and NH 3 -N removal (0.9456), indicating that only 5.18, 21.09, 5.89, and 5.44% of the total dissimilarity might not be explained by the empirical models. Zielinska et al. [10] stated that the correlation coefficient should be more than 0.80 for a good fit of a model. Moreover, the C.V.% of the obtained models for COD, colour, and NH 3 -N removals were 6.42%, 10.26%, and 7.52%, respectively, which designates an adequate model [51].
In the current study, all insignificant model terms which have limited effects were eliminated from the study to improve the model. Based on the findings, the response surface models for COD, colour, and NH 3 -N removal efficiency were constructed to predict responses, which were considered reasonable. The final regression models, in terms of their coded factors, are expressed by the second-order polynomial equations, and are presented in Table 3.
Typically, it is vital to study the effect of the operational factors on the different responses. The effect of PVDF and PAC concentration on the responses of COD, colour, and NH 3 -N removals over PVDF-PAC membranes could be evaluated using perturbation and three-dimensional (3D) response surface plots ( Figure 6). Perturbation plots show the comparative effects of independent variables on the responses. For instance, in Figure 6, the different sharp curvatures in PVDF concentration (A) and PAC concentration (B) show that the three responses (COD, colour, and NH 3 -N removal efficiency) were very sensitive to the fabrication variables, but with different behaviours. In other words, PVDF and PAC contents have a major function in the treatment process under the experimental conditions. This is another confirmation of the important effects of the independent variables (PVDF and PAC concentrations) on the treatment removal efficiency. Therefore, the 3D surface response and contour plots of the quadratic models were utilized to assess the interactive relationships between independent variables and responses. The 3D response surface was introduced as a function of PVDF and PAC concentrations. Figure 6a,c shows a symmetrical 3D surface response for both COD and NH 3 -N removals. In the meantime, the removal of colour presents a different 3D surface (Figure 6b), which indicates that colour removal was influenced differently by experimental factors than the other responses.

Pure Flux and Transmembrane Pressure Studies
By applying the factorial regression analysis on the experimental data related to PVDF-PAC membranes, both max. TMP and pure water flux responses were well agreed to a linear model of the second degree, as shown in the ANOVA analysis presented in Table 4. In a general linear model or a multiple regression model: Y = ß0 + ∑ ß + ,  Figure 6a,c indicated that the responses for COD and NH 3 -N removal rate was sufficiently enhanced upon the increase in PVDF contents in applied membranes. On the other hand, the increase in PAC contents in membrane fabrication led, to the removal of COD and NH 3 -N to some extent. It was seen that, when the PAC contents in membrane were higher than 1.0 wt.%, the removal rate for COD and NH 3 -N began to decline. According to Figure 6c, for the removal of NH 3 -N, the effect of interaction between PVDF and PAC concentrations have a noteworthy influence on removal percent. The NH 3 -N removal were gradually increased with the increasing of PAC concentration to some extent, which means that the incorporated PAC has enhanced the membrane performance in terms of NH 3 -N removal, in addition to the main separation action gained by the membrane texture itself. This good result might be ascribed to the high adsorption characteristics of the used PAC, which significantly improved the fabricated membrane efficiency [5,52].
However, PVDF concentration has limited effect on COD removal efficiency compared with the PAC content. Where 35.5 and 38.5% of COD were removed at minimum and maximum PVDF concentration (10.0 and 18.0 wt.%), respectively, 18.5 and 35.5% of COD removal were removed at minimum and medium PAC concentration (0.0 and 1.0 wt.%), respectively. Likewise, the minimum NH 3 -N removal was found to be 7.5% at membrane concentration of 10.0 wt.% PVDF and 2.0 wt.% PAC, while the maximum NH 3 -N removal (24.5%) was observed at the PVDF and PAC concentration of 18.0 and 1.0 wt.%, respectively.
On the other hand, the 3D response surface in Figure 6b displays a different effect of interaction between the experimental factors on the colour removal rates. It was observed that an increase in the concentration of PVDF in the membrane leads to an improvement in the colour removal to some degree. When the concentration of PVDF was higher than (14 wt.%), the colour removal performance starts to decrease. This behaviour was credited to the combined effects of additive and polymer in the dope. This leads to creating large volume voids with increasing polymer dosage, and lets the fine particles from contaminants to permeate through the membrane [50]. Meanwhile, the enhancement of PAC concentration in a membrane drove a steady increase in the colour removal efficiency. As witnessed in Figure 6b, the predicted minimum and maximum efficiencies of colour removal were 15.0 and 56.5% present at fabrication concentrations of (18.0 wt.% PVDF, 0.0 wt.% PAC), and (14.0 wt.% PVDF, 1.0 wt.% PAC), respectively. This also confirms the effectiveness of PAC content in enhancing the removal performance of the filtration process using PVDF fabricated membrane.
Despite the incorporation of PAC into membrane enhancing the COD, colour, and NH 3 -N removal, the filtered leachate still did not meet the Malaysian Discharge Standard (Table 1). This is due to the highly concentrated pollutants of leachate that resulted in a reduction in membrane efficiency owing to the clogging caused by influent SS component. Therefore, a pre-treatment process such as PAC adsorption is suggested to be used before the membrane treatment [33].

Pure Flux and Transmembrane Pressure Studies
By applying the factorial regression analysis on the experimental data related to PVDF-PAC membranes, both max. TMP and pure water flux responses were well agreed to a linear model of the second degree, as shown in the ANOVA analysis presented in Table 4.
In a general linear model or a multiple regression model: where: Y is the response, X i is the independent factor, k is the number of variables, ß 0 is the constant term, ß i represents the coefficient of the linear, and ε is the random error or noise [53]. The final linear models obtained for each response has been expressed by the first order polynomial equation, as presented in the last raw of Table 4.
The fitted model for the pure water flux suggests a large F-value (53.56), suggesting that the model is significant. As the value of Prob > F of all terms is less than 0.050, this suggests that all the model terms are significant. Based on their F-values, the PVDF concentration term (A) has the highest influence on the model, followed by PAC concentration term, and lastly the combination term. The term of PAC concentration presents a positive effect on pure flux, while the other two terms have been found to be negative influencers. Hence, the pure water flux was raised only with enhancing PAC contents in the membrane while, in contrast, it is decreased with the increasing of the PVDF content of a membrane. On the other hand, the suggested model of max. TMP was significant with a high F-value (49.62), as can be seen from Table 4. Based on its effect on the model from the highest to the lowest, the model terms can be arranged as follows: PVDF content, PAC content, and the combination of both, with F-values of 131.07, 13.01, and 4.78, respectively. However, the PVDF concentration is the only factor which showed a positive influence on the max. TMP, due to the positive sign of its term; this indicates a worse impact on the max. TMP, as it could be increased with the increasing of PVDF content on the fabricated membrane. On the other hand, PAC concentration exhibited a better effect on the max. TMP, which showed a reduction in max. TMP occurred due to the increasing of the PAC content.
Additionally, both of the models display a non-significant lack of fit F-value, which indicates that well fitted models have been selected to present the experimental results with minor pure errors [15].
The R 2 values obtained in the present study for pure flux and max. TMP were 0.9470 and 0.9430, respectively. The high value of R 2 represents good agreement between the observed and the calculated results within the experimental ranges [37]. Moreover, C.V. % for the water flux and TMP were 10.50% and 6.56%, respectively. Where these small values indicate good fitness of the models [51].
Based on these findings, the resulted response surface models in the current work for predicting the two responses (pure flux and max. TMP) were considered reasonable.
The influence of integrated PAC and the interaction of content's concentrations on the max. TMP can be explored by the plots of perturbation and 3D response surface, as shown in Figure 7. From perturbation plots at Figure 7, it is easy to notice that pure flux and max. TMP responses are very sensitive to the experimental factors, and to conclude that both have a different (inversed) behaviour regarding the PVDF and PAC concentration values. As can be seen from Figure 7a, increasing of PVDF concentration (A) resulted in a linear decrease in pure water flux and increase in max. TMP, which attributed to the reduction in membrane porosity due to the increase in polymer concentration, which is well recognized for the system of a single polymer casting solution [50]. However, PAC concentration (B) showed a different effect, as any increase in its value causes a linear increment on the pure water flux, but a decrease in max. TMP.

Fabricated Membrane Characterization
The morphology of produced membrane can explain the effect of dope composition on membrane performance. A collection of membranes composed from different concentrations of PVDF and PAC (wt.%) were chosen from the fabricated membranes to represent the different membrane compositions, and consequently to be investigated by the morphological studies. These membranes were: FM1 with the content of (10.

Fabricated Membrane Characterization
The morphology of produced membrane can explain the effect of dope composition on membrane performance. A collection of membranes composed from different concentrations of PVDF and PAC (wt.%) were chosen from the fabricated membranes to represent the different membrane compositions, and consequently to be investigated by the morphological studies. These membranes were: FM1 with the content of (10. The FTIR spectrum of PVDF-PAC fabricated membranes with the various compositions is illustrated in Figure 8. It is clearly observed from Figure 8 that membranes displayed semi-typical distinctive spectra along the range of 4000 and 400 cm −1 . Characteristic chemical groups are witnessed in the band of all membranes at waves with lengths 3020, 2990, 2370, 1400, 1070, 875, 590, and 490 cm −1 with altered vibrations of strength depends on the different membrane compositions. The spectrum shows bands at 2990 and 3020 cm −1 which are attributed to the symmetric and asymmetric stretching vibrations of C-H coming from ketones and carboxylic acids [54], where vibrations at 1070 and 1400 cm −1 presented the deformation peaks of C-F related to PVDF.

R PEER REVIEW
17 of 26 coming from ketones and carboxylic acids [54], where vibrations at 1070 and 1400 cm −1 presented the deformation peaks of C-F related to PVDF. The notable peaks of the various membranes at 2370, 875, and 590 cm −1 , assigned to CO2, CO3 −2 , and C-O-groups, respectively, were the features distinctive of neutralization methanol, used after membrane casting [55]. Moreover, the OH group detected at 490 cm −1 is attributed to the DW used for membrane solidification during the casting process [56]. Figure 8 also confirmed that the recorded wave numbers in the spectrum of both membranes without PAC (FM1 and FM4) have higher frequencies in comparison with the spectrums of the other two membranes with PAC content (FM2 and FM3).
Furthermore, it could be observed that the peaks of the membrane with higher content of PAC (FM2) have lower vibrations compared with the membranes with lower PAC content (FM3). Evidently, the peaks become narrow with less strength at the increasing of PAC weight, indicating that the hydrogen bonds were constructed well between PVDF polymer chains and the hydroxyl groups from PAC, which reduces the PVDF hydrophobic tendency [57]. These outcomes confirmed that PAC was well integrated to PVDF membranes, and partially relocated on the membrane surface, which leads to membrane treatment efficiency enhancement.
To investigate the elemental composition present in the fabricated PVDF-PAC membranes with different compositions, EDX analysis was recorded in the binding energy region from 0 to 15 keV as exhibited in Figure 9. The PVDF characterized elements C and F The notable peaks of the various membranes at 2370, 875, and 590 cm −1 , assigned to CO 2 , CO 3 −2 , and C-O-groups, respectively, were the features distinctive of neutralization methanol, used after membrane casting [55]. Moreover, the OH group detected at 490 cm −1 is attributed to the DW used for membrane solidification during the casting process [56]. Figure 8 also confirmed that the recorded wave numbers in the spectrum of both membranes without PAC (FM1 and FM4) have higher frequencies in comparison with the spectrums of the other two membranes with PAC content (FM2 and FM3).
Furthermore, it could be observed that the peaks of the membrane with higher content of PAC (FM2) have lower vibrations compared with the membranes with lower PAC content (FM3). Evidently, the peaks become narrow with less strength at the increasing of PAC weight, indicating that the hydrogen bonds were constructed well between PVDF polymer chains and the hydroxyl groups from PAC, which reduces the PVDF hydrophobic tendency [57]. These outcomes confirmed that PAC was well integrated to PVDF membranes, and partially relocated on the membrane surface, which leads to membrane treatment efficiency enhancement.
To investigate the elemental composition present in the fabricated PVDF-PAC membranes with different compositions, EDX analysis was recorded in the binding energy region from 0 to 15 keV as exhibited in Figure 9. The PVDF characterized elements C and F were clearly observed in the spectra of the pure PVDF membranes (without PAC), while the AL element, which characterizes the presence of PAC, appeared only at the PVDF membranes incorporated with PAC [33]. Figure 9b,d shows the EDX analysis of 2.0 and 1.0 wt.% PAC, respectively. It is clearly witnessed that the presence of PAC was presented well.  Table 5 shows the atomic percentages of the different elemental compositions of the selected membranes (FM1-FM4). From the EDX findings, the weight percentages of elemental AL on the FM2 and FM3 were determined as 1.04 and 0.79, respectively, which confirmed the presence PAC with representative weights on the integrated membranes.  Figure 10 presents the FESEM images for produced membranes with different compositions, which show the top surface morphology of membranes, along with its cross section. As can be seen from Figure 9a-d, there were many small pores available on the surface of FM1 membrane which contains the lowest PVDF polymer content (10.0 wt.%).  Table 5 shows the atomic percentages of the different elemental compositions of the selected membranes (FM1-FM4). From the EDX findings, the weight percentages of elemental AL on the FM2 and FM3 were determined as 1.04 and 0.79, respectively, which confirmed the presence PAC with representative weights on the integrated membranes.  Figure 10 presents the FESEM images for produced membranes with different compositions, which show the top surface morphology of membranes, along with its cross section. As can be seen from Figure 9a-d, there were many small pores available on the surface of FM1 membrane which contains the lowest PVDF polymer content (10.0 wt.%). Furthermore, the number and size of these pores start to be decreased, first on membrane FM3, with PVDF content 14.0 wt.% and PAC content 1.0 wt.%, followed by FM2 membrane with the highest PAC content (2.0 wt.%), while the membrane FM4 has a semi-impermeable surface due to its high PVDF polymer content (18.0 wt.%) with no PAC content. This was in agreement with the findings earlier discovered by Kunst and Sourirajan [58]. Furthermore, an AFM test was carried out to investigate the membrane top surface, along with its roughness, as shown in Figure 11. The FM2 membrane might contain some extra PAC particles which made its top surface rougher compared to others (Figure 11b). Having less depth of facial peaks and valleys, the FM4 membrane surface (Figure 11d) is relatively smooth due to containing only PVDF polymer which received a homogeneous Referring to membrane cross sections on Figure 10a-d, all membranes display the formation of macrovoid with loosely packed structures. Typically, the membrane consists of two layers, which are a spongy porous support layer and a dense top finger-like layer. The establishment of these configurations can be attributed to the instantaneous demixing of polymer and solvent during the process of phase inversion.
FM1 membrane, with only PVDF and the weight of 10.0 wt.%, displayed an unimproved finger-like formation and a sponge-like support layer containing large, unconnected pores, delimited by polymer walls (see Figure 10a). The finger-like voids turn become flat, bigger, and even strained to the bottommost of the fabricated membranes with an increase in PAC concentration (i.e., in FM2 and FM3), and the spherical voids of the spongelike structures connect more closely with themselves (Figure 10b,c). However, the FM4 membrane, containing the highest concentration of PVDF, gives thin, smaller, and nonstretched figure-like pores, with less connection to the little sponge-like pores located on the cross-section's bottom. This produces low membrane flux due to the greater amount of polymer contributing a higher membrane viscosity, which lead to a decease in the membrane porosity and pore size. The overall FESEM micrographs have proved the significant effect of the PAC presence in improving the fabricated membrane characteristic in terms of membrane rejection, and therefore removal rate of contaminants.
Furthermore, an AFM test was carried out to investigate the membrane top surface, along with its roughness, as shown in Figure 11. The FM2 membrane might contain some extra PAC particles which made its top surface rougher compared to others (Figure 11b). Having less depth of facial peaks and valleys, the FM4 membrane surface (Figure 11d) is relatively smooth due to containing only PVDF polymer which received a homogeneous mixing at the preparation phase of dope solution [59]. However, the peaks and valleys of FM1 and FM3 membranes reduced gradually compared to FM2, where FM3 has the smoothest surface compared with other membranes (see Figure 11a-d). To confirm all above observations, the values of membrane surface roughness (R q and R a ) given in Figure 11 can be considered.
Water 2021, 13, x FOR PEER REVIEW 21 of 26 above observations, the values of membrane surface roughness (Rq and Ra) given in Figure  11 can be considered. For membrane permeability analysis, the impact of PAC addition to membrane permeability, in terms of porosity and average pore size, were evaluated for the produced PVDF membranes. As presented in Table 6, the porosity and average pore size of the fabricated PVDF membranes incorporated with PAC were higher than the other membranes without PAC. Based on Table 6, the resulted fabricated membranes were "micro-filtration", and the highest mean values of porosity and average pore size were achieved by For membrane permeability analysis, the impact of PAC addition to membrane permeability, in terms of porosity and average pore size, were evaluated for the produced PVDF membranes. As presented in Table 6, the porosity and average pore size of the fabricated PVDF membranes incorporated with PAC were higher than the other membranes without PAC. Based on Table 6, the resulted fabricated membranes were "micro-filtration", and the highest mean values of porosity and average pore size were achieved by FM2 membrane at the values 77.48% and 24.43 µm, respectively. On the other hand, the lowest values of the same corresponding permeability parameters were found using membrane FM4 at 48.38% and 12.15 µm, respectively. These findings agreed with the above morphological results.

Membrane Treatment Optimization
The best synthesized membrane has been selected using the RSM tool, where the membrane efficiencies of COD, colour, and NH 3 -N removal were optimized during this study.
Based on the DoE software, the operational conditions (PVDF weight and PAC weight) were targeted to be within the range. While the dependents of treatment performance (COD, colour, and NH 3 -N removal) were chosen as ''maximum" to achieve the ultimate filtration treatment. The other responses were remained "within the range". Accordingly, the optimization tool assimilates the singular desirability into a particular number, and then aims to optimize the function.
Consequently, the composition of the optimum membrane, together with respective rates of removal efficiency, were obtained. The optimum removals and their corresponding water flux and max. TMP are presented in Table 7. The membrane with 14.9 wt.% of PVDF and 1.0 wt.% of PAC was found to be the optimum, and thus selected as the best membrane design, having optimum removal efficiency according to its highest desirability (0.870) [60].
As shown in Table 7, 35.34, 48.71, and 22.00% removal of COD, colour and NH 3 -N, respectively, was predicted by the software under optimized operational conditions. The corresponding (non-optimized) water flux and max. TMP were found at the val-ues 61.00 L/m 2 ·h and 0.67 bar, respectively. An additional experimentation was then performed to confirm the optimum findings.
As illustrated in Table 7, the error column indicates the differences between the predicted and laboratory values, which shows that the lab experiments agree well with the response values estimated by the software. However, less agreement between the predicted and the laboratory result was obtained in case of NH 3 -N removal (8.36% error).

Membrane Performance Comparison with Other Reported Studies
The performance of the optimum fabricated membrane with other reported PVDF produced membranes is shown in Table 8. It can be noticed from Table 8 that the current study offered the smoother surface among the existing works based on the average roughness (R a = 36.39 nm), which accordingly improves the removing performance and antifouling properties of the created membrane [61]. There exists few values of pure flux that are higher than the achieved in the current work, such as the flux of 143.24 L/m 2 ·h produced by Penboon et al. [62]. The low value of pure water flux of the current work (61.00 L/m 2 ·h) could be ascribed to the differences in the experimental characteristics such as the type of wastewater or the rates of feed flow. In addition, the rejection efficiency in the current work is lower than previously reported, which could be solved through further enhancement of the produced membranes using hydrophilic additives such as PVA or PVP [25,57]. Based on previous studies, after saturation, membrane corroborated PAC can be washed back and reused [13,63].

Conclusions
The adsorbent material PAC was used to fabricate a novel PVDF membrane for the treatment of stabilized landfill leachate. The fabricated PVDF flat sheet membranes integrated with PAC showed better performance when compared with PVDF membrane (without PAC). The addition of PAC effectively enhanced the removal rate and the fouling control parameters of produced membranes. Increasing PAC content to a certain value has a positive influence on the removal efficiency of COD, colour, and NH 3 -N, as well as on membrane characteristics. Operational optimization was performed using RSM to select the optimum membrane design in terms of the removal efficiency. The best membrane composition was found at (14.9 wt.%) PVDF and (1.0 wt.%) PAC, which removed 36.63% of COD, 49.50% of colour, and 23.84% of NH 3 -N. This was in agreement with the predicted removals. The corresponding experimental values of water flux and max. TMP also agreed with the prediction, with the values of 61.10 L/m 2 ·h and 0.64 bar, respectively. The performance and structure of fabricated membranes were investigated by filtration tests, FTIR, FESEM, and AFM spectroscopy. In general, this work shows the potential of treatment and hydrophilic improvement of hydrophobic PVDF polymer membranes using PAC. For further removal efficiency, membrane properties or practice could be improved by either adding a hydrophilic material, or applying pre-treatment process such as adsorption via PAC.