Next Article in Journal
Uptake of Sb(V) by Nano Fe3O4-Decorated Iron Oxy-Hydroxides
Previous Article in Journal
Effect of Aluminium Salt Dosing on Activated Sludge Settleability Indicators: A New Settleability Model Development
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling of Pollutants Removal in Subsurface Vertical Flow and Horizontal Flow Constructed Wetlands

by
Wojciech Dąbrowski
1,*,
Beata Karolinczak
2,*,
Paweł Malinowski
1 and
Dariusz Boruszko
1
1
Faculty of Civil and Environmental Engineering, Bialystok University of Technology, 45E Wiejska St., Bialystok 15-351, Poland
2
Faculty of Building Services, Hydro and Environmental Engineering, Warsaw University of Technology, Nowowiejska 20, Warsaw 00-653, Poland
*
Authors to whom correspondence should be addressed.
Water 2019, 11(1), 180; https://doi.org/10.3390/w11010180
Submission received: 28 December 2018 / Revised: 16 January 2019 / Accepted: 19 January 2019 / Published: 21 January 2019
(This article belongs to the Section Wastewater Treatment and Reuse)

Abstract

:
Reject water is a by-product of every municipal and agro-industrial wastewater treatment plant (WWTP) applying sewage sludge stabilization. It is usually returned without pre-treatment to the biological part of WWTP, having a negative impact on the nitrogen removal process. The current models of pollutants removal in constructed wetlands concern municipal and industrial wastewater, whereas there is no such model for reject water. In the presented study, the results of treatment of reject water from dairy WWTP in subsurface vertical flow (SS VF) and subsurface horizontal flow (SS HF) beds were presented. During a one-year research period, SS VF bed reached 50.7% efficiency of TN removal and 73.8% of NH4+-N, while SS HF bed effectiveness was at 41.4% and 62.0%, respectively. In the case of BOD5 (biochemical oxygen demand), COD (chemical oxygen demand), NH4+-N, and TN (total nitrogen), the P-k-C* model was applied. Multi-model nonlinear segmented regression analysis was performed. Final mathematical models with estimates of parameters determining the treatment effectiveness were obtained. Treatment efficiency increased up to the specific temperature, then it was constant. The results obtained in this work suggest that it may be possible to describe pollutant removal behavior using simplified models. In the case of TP (total phosphorus) removal, distribution tests along with a t-test were performed. All models predict better treatment efficiency in SS VF bed, except for TP.

1. Introduction

In Europe, the biggest municipal and agro-industrial wastewater treatment plants (WWTPs) utilize anaerobic sewage sludge digestion with biogas production [1]. The possibility of heat and electric energy production from biogas is an advantage of such sewage sludge treatment, which decreases the cost of waste treatment significantly.
One of the major problems connected with anaerobic digestion chamber exploitation is reject water, which is generated in the process of dewatering stabilized sewage sludge. It is usually returned without pre-treatment to the biological part of a WWTP, having a negative impact on the treatment process and causing problems with reaching high effectiveness of nitrogen removal [2,3,4]. Reject water from anaerobic sewage sludge digestion is characterized by a high concentration of nitrogen in the form of ammonium nitrogen (NH4+-N) and irregular flow because of periodically working devices for sewage sludge dewatering [5,6,7,8]. The problem of reject water treatment in a conventional sludge activated system is caused by low biodegrability, indicated by proportions of easily biodegradable organic matter expressed by biochemical oxygen demand (BOD5) to total organic matter expressed by chemical oxygen demand (COD).
Constructed wetlands are considered environmentally friendly technology [9,10,11]. They have been used worldwide for the treatment of municipal and industrial wastewater, as well as reject water from anaerobic sewage sludge digestion in municipal WWTPs and landfill leachate, which is also characterized by high NH4+-N concentration [12,13,14,15]. Reject water generated during anaerobic digestion in dairy WWTPs differs from that from municipal WWTPs. There is a lack of experiments concerning constructed wetlands treatment of reject water generated during anaerobic digestion in dairy WWTPs [16].
The aim of the research was to define parameters, with their expected values, of the mathematical models describing the pollutants removal in subsurface vertical flow (SS-VF) and subsurface horizontal flow (SS-HF) constructed wetlands treating reject water from anaerobic sludge digestion in dairy WWTPs. The P-k-C* model was applied [12,17,18]. The issue of modeling the functioning of constructed wetlands is one of the most important challenges in that field [19,20]. The current models of pollutant removal concern the municipal and industrial wastewater treatment, whereas there is no such model for reject water treatment, so the necessity for such research goals was confirmed.

2. Materials and Methods

2.1. Study Sites

The research was conducted over a twelve-month period using a pilot-scale plant built in the biggest dairy plant in Poland, which Person Equivalent (PE) is up to 500,000, located in Wysokie Mazowieckie in the Podlaskie province. The average air temperature is 6.8 °C and precipitation is 562 mm year−1.
A problem with high concentration of nutrients in the reject water was observed after changing aerobic sewage sludge stabilization to anaerobic. After sewage sludge dewatering with centrifuge, reject water was returned to the main sewage line without separate treatment. In comparison with reject water obtained during aerobic sewage sludge stabilization, the concentration of ammonium nitrogen increased significantly [2].
The research installation (N 52.54.20-E 22.299.20) was designed by the authors using earlier experience with reject water treatment. Its main elements are SS VF and SS HF beds with 5 m2 surface area and 0.8 m height (Figure 1). The beds were used in parallel. In addition, the installation includes a sedimentation and retention tank, outflow, and sampling points (I, II, III). Figure 2 presents the cross section of SS VF and SS HF beds.
The SS VF bed was built of three layers of stones 0.2 m (fraction 16–60 mm), gravel 0.4 m (fraction 2–16 mm), and sand 0.2 m (fraction 0.5–2 mm). Two pipes (diameter 50 mm) with perforation 2 mm were used for passive aeration. The SS-HF bed was built of gravel (2–16 mm) and stones (16–60 mm). The gravel hydraulic conductivity was 4 × 10−2 ms−1. Both beds were planted with reeds (Phragmites australis). Reject water was taken directly from centrifuge which dewatered sludge after anaerobic digestion in a dairy WWTP.

2.2. Sampling and Analytical Procedures

The study was carried out between April 2015 and March 2016. Both beds were supplied from sedimentation and retention tanks with the same amount of reject water, which allowed comparison of these beds’ effectiveness. The daily hydraulic load was 0.1 m3 m−2 d−1 in both SS VF and SS HF beds. Hydraulic retention time (HRL) for SS HF bed was approximately 8 days. Samples were collected three times a month (influent to SS-VF and SS-HF and effluents from both beds). The air temperature during the research period varied from −11 to 26 °C, while reject water temperature varied from 4 °C to 20 °C.
The basic physical and chemical analyses were performed: BOD5, COD, total organic carbon (TOC), total suspended solids (TSS), total Kjeldahl nitrogen (TKN), ammonium nitrogen (NH4+-N), nitrate nitrogen (V) (NO3-N), nitrite nitrogen (III) (NO2-N), total phosphorus (TP), dissolved oxygen, and alkalinity. TN value was calculated as a sum of TKN, NO3–N, and NO2-N. To evaluate biodegrability of reject water, BOD5/COD and BOD5/TN ratios were determined. Determinations were conducted in a certified laboratory in accordance with the procedures set out in the Regulation of the Environmental Protection Minister [21] from November 18, 2014, and in accordance with the American Public Health Association (2005) [22].

2.3. Modeling of Pollutants Removal

Statistical analysis of NH4+-N and TN, as well as BOD5; COD removal was performed, based on the P-k-C* model [12]:
C o u t C C i n C = ( 1 + k ( T ) q P ) P
where Cout—output concentration (g/m3), Cin—input concentration (g/m3), C*—background concentration (g/m3), k(T)—chemical reaction coefficient (temperature dependent), q—hydraulic load (m/d), P—number of tanks in series.
Model (1) can be further simplified if some assumptions are made. The apparent removal efficiency η a p p can be defined as:
η a p p = 1 C o u t C C i n C
This apparent efficiency model’s real efficiency η in limit C * 0 . If limit P is taken, the right side of Equation (1) simplifies:
lim P ( 1 + k ( T ) q P ) P = exp ( k ( T ) / q )
Chemical reaction coefficient k is dependent on reject water temperature, using direct or modified first order Arrhenius dependency:
k d ( T ) = K 20 θ T 20 exp ( ln K 20 + [ T 20 ] ln θ )
k m ( T ) = K 20 θ T 20 θ m min ( T T k , 0 ) exp ( ln K 20 + [ T 20 ] ln θ + [ min ( T T k , 0 ) ] ln θ m )
where k d , k m —direct and modified chemical reaction coefficients; K 20 —direct or modified chemical reaction coefficients normalized for 20 °C; θ —temperature coefficient; θ m —modifier for temperature coefficient for temperatures less than specific temperature T k , where trend changes; min ( · , · ) —minimum function of 2 arguments.
Modified first order dependency allows for an additional change of the reaction coefficient below an estimated specific temperature. It was introduced to better fit modeled data, while not changing the basic rule of first order linear dependency on a logarithmic scale.
By using all available combinations of previously described equations, 16 models with different parameters subsets were estimated. The final parameter subset selection was performed using methodology given by Burnham & Anderson [23]. A detailed description of fitted models, their selection and fitting procedure are presented in Appendix A.
It is worth noting that a full analysis of a P-k-C* model is difficult and can lead to unexpected results [24]. In most of the cases, plug flow models can be too simple to describe the dynamics of particular pollutant removal. Sensitivity analysis of fitted models in literature [25] may suggest that some parameters, like P , may be insensitive and set to values obtained from earlier studies. There are some earlier results suggesting that temperature dependency can change [16].

3. Results and Discussion

3.1. Treatment Efficiency and Load Removal

Table 1 presents the characteristics of reject water before and after treatment with SS VF and SS HF beds during research period, while Table 2 shows removed pollutants load.
The BOD5/COD and BOD5/TN ratios give information about biodegradability [26]. Its value can be used to assess the reject water susceptibility for high efficiency of biological treatment in conventional systems (e.g., activated sludge method). An average BOD5/COD ratio was 0.59, while BOD5/TN ratio 0.43. Low BOD5/COD ratio pointed out low degradability of the organic compounds. Reject water from dairy WWTP with anaerobic sewage sludge digestion was discovered to have a higher BOD5/COD ratio than municipal WWTPs. In municipal, it ranges from 0.25 to 0.32 for reject water from the WWTP in Gdansk [27], and 0.2 for reject water from the WWTP in Minworth [8]. The WWTP in Gdansk records BOD5/TN ratios at the level of 0.37 up to 0.54. In this case, the conventional path of nitrogen removal cannot be applied as the content of easily biodegradable organics is insufficient [28,29]. A very high concentration of ammonium nitrogen should be viewed as the main reason for this.
The average calculated treatment efficiency in SS VF bed was: BOD5 82.12%, COD 79.79%, TOC 80.85%, TSS 85.30%, TN 51.46%, NH4+-N 74.06%, and TP 25.42%. The average calculated treatment efficiency in SS HF bed was: BOD5 74.13%, COD 75.34%, TOC 66.01%, TSS 89.90%, TN 42.58%, NH4+-N 62.52%, and TP 37.20%. After the treatment, the average ratios of BOD5/COD were 0.53 for SS VF bed and 0.63 for SS HF bed, and BOD5/TN were 0.16 and 0.20, respectively.
Through the analysis of the effect of organic matter removal (g m−2 d−1) expressed by BOD5, COD, and TOC (Table 2), similar values to household sewage treatment with constructed wetlands were achieved [10]. In the case of TN and NH4+-N, higher efficiency was observed mainly due to high concentration of nitrogen in reject water before treatment in SS VF and SS HF beds. The TP removal effect was similar to the observed in the case of household and municipal sewage treatment.

3.2. Modeling of Pollutants Removal

Statistical analysis of NH4+-N and TN, as well as BOD5, COD removal, presented in Appendix B, revealed existence of 3 group models, each with stable core parameters, which dominate the variability of logarithmic likelihood and residuals distribution. The most-preferred parameter set was T k , θ m , K 20 —a model without temperature dependency after T k , with Akaike weight greater than 0.94 in all but other cases except in one (TN removal in SS HF bed). Residual analysis suggests also T k , θ m , K 20 model selection in this case. Selected models are presented graphically in Figure 3, using efficiency scale.
In the cases of BOD5, COD, TN, and NH4+-N, the treatment efficiency increased until a specific temperature T k was reached, and then it was constant. Results obtained in this work suggest that it may be possible to effectively describe pollutant removal behavior in a consistent way using simplified models. Any additional parameter, with assumption of proper overall fit, will have a large confidence interval, suggesting overall insensitivity. Simplicity of obtained models does not interfere with their statistical validity.
Created models show that simple ratios obtained over the whole period of the research study are insufficient to properly describe behavior of beds, and a more dedicated approach, like P-k-C* modeling, is necessary. Only ratios in stable periods of bed work (after specific temperature T k ) are sufficient to describe it.
Values of specific temperature T k were marked in place of trend change and on the upper and lower section of figures along with 95% confidence intervals (CI). Parameter estimates of the fitted models, along with their 95% confidence intervals, and residual sum of squares were presented in Table 3.
Figure 4 presents standard a box-whisker plot of TP removal efficiency and reject water temperature during the whole research period. Mean values are described as red dots inside each box, and horizontal line represents median. Dashed lines inside efficiency vs. temperature plot represent mean values.
The calculated treatment efficiency of TP removal was on average 25.42% (SS VF) and 37.20% (SS HF). The efficiency of TP removal in SS VF bed was stable and higher than in SS HF bed. Phosphorous removal occurs mainly in the processes of chemical precipitation and sorption, which do not depend on temperature [12]. A distribution analysis of TP treatment efficiency revealed no significant deviation from normality (Shapiro-Wilk normality test, for SS HF: test statistics 0.952, p-value = 0.1, for SS VF: test statistics 0.96, p-value = 0.19). An F-test revealed no statistically significant difference between the variances (F test statistics 1.74, numerator and denominator degrees of freedom—37, p-value = 0.095). T-test revealed statistically significant differences between mean treatment efficiency of TP (t statistics 9.52, 74 degrees of freedom, p-value < 2×10−14). For all performed tests, significance level was set to α = 0.05.
The mathematical models found practical implementation. An application, which allows determining the effectiveness of reject water treatment in SS VF and SS HF constructed wetlands depending on the reject water temperature, was developed (Figure 5).

4. Conclusions

The study showed a high efficiency of SS VF and SS HF beds for removing main pollutants from the reject water generated during anaerobic sewage sludge digestion in dairy WWTP. Results for SS VF and SS HF, respectively, were: COD 79.79% and 75.34%, for BOD5 82.12% and 74.13%, for TN 51.46% and 42.58%, for NH4+-N 74.06% and 62.52%, for TP 25.42% and 37.20%. A higher efficiency of main organic pollutants removal was observed in the case SS VF bed during the whole research. The efficiency of TP removal was stable and higher for the SS HF bed.
The results have been used for modeling the work of SS VF and SS HF beds for reject water treatment. Models for determining the effectiveness of treatment depending on the temperature and type of beds were presented. P-k-C* model was sufficient to describe the observed variability of constructed wetland behavior. In the cases of BOD5, COD, TN, and NH4+-N, the treatment efficiency increased to the specific temperature, then it was constant. Results obtained in this work suggest that it may be possible to effectively describe pollutant removal behavior using simplified models.
All models predict better treatment efficiency in SS VF bed, except for TP. In the case of TP removal, distribution tests along with a t-test were performed. The created models show that simple ratios obtained over the whole period of study are insufficient to properly describe behavior of beds. A more dedicated approach, like P-k-C* modeling, is necessary. Only ratios in stable periods of bed work (after specific temperature) are sufficient to describe it.
An application in the form of software for modeling the efficiency of reject water treatment from sewage sludge digestion in dairy WWTPs was prepared on the basis of the results from the research.

Author Contributions

Conceptualization, W.D. and B.K.; Methodology, W.D., D.B., B.K.; Software, P.M.; Validation, P.M.; Formal Analysis, P.M., W.D.; Investigation, W.D., B.K.; Data Curation, W.D., P.M.; Writing—Original Draft Preparation, B.K., W.D.; Writing—Review and Editing, B.K., D.B.; Visualization, P.M.; Supervision, W.D.

Funding

This research received no external funding.

Acknowledgments

The study was conducted as a research project S/WBiIŚ/3/2014 in Faculty of Building and Environmental Engineering of Białystok University of Technology (Poland), and financed by Ministry of Science and Higher Education of Poland. Authors would like to thank dairy plant Mlekovita in Wysokie Mazowieckie for supporting research conducted due to scientific and technical cooperation between BUT and Mlekovita.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Statistical Modeling

The 16 models obtained from (1) were constructed by taking limit of certain variables simplifying the expression:
  • C * (allowed to vary or set to 0)
  • P (allowed to vary or taken limit resulting in (3) dependency)
  • θ (allowed to vary or set to 1 )
  • θ m and related T k (if θ m is present, then (5) form was used, otherwise (4)).
All models were estimated using non-linear least squares procedure with a variant of Levenberg-Marquardt’s algorithm for optimization [30]. It is worth noting that both sides of (1) are always constrained to a (0–1) interval. Residuals from modeling this dependency will be therefore non-homoscedastic. Final model estimation procedure reweighted them by term proportional to [31]:
s i ~ 1 r ^ i ( 1 r ^ i )
where r ^ i —estimated and right-hand side of (1), s i —residual correction for i-th observation.
Model parameter selection was performed using methodology from [23]. After all 16 models were built, they were weighted using Akaike weight s w i , which is a measure of a model fit to the data, normalized to (0–1) interval, relative to the best model from the set of analyzed ones. The Akaike weight was calculated from differences Δ i of second order [32] AIC [33]. Models with very similar logarithmic likelihood and residual distribution, containing only additional parameters, were removed from weighting procedure [23]. Top weighted models, forming 0.90–0.95 of summary Akaike weights, were retained for further analysis. Residuals structure of the remaining parameter sets was also checked, to further reject oversimplified ones using an approach based on [34]. Remaining parameter subsets were presented along with their basic statistics. Parameter T k , if present, was presented graphically in a similar way to [35].
Model calculation was performed using R version 3.5.1 [36]. The nlfb procedure from nlsr package [37] was used for optimization of models. Second order AIC was calculated using AICc procedure from AICcmodavg package [38].

Appendix B. Result of Statistical Modeling

Table A1, Table A2, Table A3 and Table A4 present fitted models grouped by almost identical values of logarithmic likelihood and residuals structure. Analysis revealed three preferred groups of models, each with stable core parameters, which dominate the variability of logarithmic likelihood and residuals distribution. For the simplest models from each group, containing the core parameters, Akaike weight was calculated. Any other model, with addition of one or more parameters, does not change its fit significantly.
Analysis results consistently supported choice of model with parameter subset T k , θ m , K 20 —a model with no temperature dependency after T k . In almost all cases, Akaike weight for this preferred model is greater than 0.95, with suggests strong support. The second one is θ , K 20 —a model with non-changing temperature dependency. Model θ , K 20 for total nitrogen for horizontal flow bed had approximately w i equal to 0.16, suggesting some support for further analysis. The third model contains only K 20 variable, and obviously does not depend on temperature. It has almost no support, and was not chosen for residual analysis.
Residual analysis, presented in Figure A1, Figure A2, Figure A3 and Figure A4, support selection for first preferred model in all cases; there is signal of non-linearity in models θ , K 20 , even a slight one for total nitrogen in horizontal flow bed. Residuals analysis suggests that models with T k , θ m , K 20 parameters also have properties allowing for further analysis of those models parameters—normal distribution, homoscedasticity and no visible autocorrelation—which were also performed.
Table A1. Fitted models for BOD5 with their statistics.
Table A1. Fitted models for BOD5 with their statistics.
Horizontal Flow Bed SS-HFVertical Flow Bed SS-VF
ModelVariablesklogLikAICc Δ i w i ModelVariablesklogLikAICc Δ i w i
1 T k , θ m , K 20 42.2724.6680.000>0.9991 T k , θ m , K 20 4−0.85510.9230.000>0.996
2 P , T k , θ m , K 20 52.4856.9052.237 2 θ , T k , θ m , K 20 5−0.26012.3941.471
3 θ , T k , θ m , K 20 52.2807.3152.6483 P , T k , θ m , K 20 5−0.60613.0882.164
4 C * , T k , θ m , K 20 52.2727.3302.6634 C * , T k , θ m , K 20 5−0.85513.5862.663
5 θ , P , T k , θ m , K 20 62.4859.7395.0725 θ , P , T k , θ m , K 20 6−0.03114.7723.849
6 C * , P , T k , θ m , K 20 62.4859.7405.0726 C * , θ , T k , θ m , K 20 6−0.26015.2294.306
7 C * , θ , T k , θ m , K 20 62.28010.1505.4837 C * , P , T k , θ m , K 20 6−0.60615.9224.999
8 C * , θ , P , T k , θ m , K 20 72.48512.7638.0968 C * , θ , P , T k , θ m , K 20 7−0.03117.7966.873
9 θ , K 20 3−8.79324.29319.6255.476 × 10−59 θ , K 20 3−7.78311.34911.3493.421 × 10−3
10 θ , P , K 20 4−7.89725.00720.339 10 θ , P , K 20 4−6.76322.73811.815
11 C * , θ , K 20 4−8.79326.79922.13111 C * , θ , K 20 4−7.78324.77813.855
12 C * , θ , P , K 20 5−7.89727.67023.00212 C * , θ , P , K 20 5−6.76325.40114.478
13 K 20 2−31.51067.36462.6965.430 × 10−1413 K 20 2−32.05268.44657.5233.218 × 10−13
14 C * , K 20 3−31.51069.72765.059 14 C * , K 20 3−32.05259.88670.809
15 P , K 20 3−31.51069.72765.05915 P , K 20 3−32.05259.88670.809
16 C * , P , K 20 4−31.51072.23367.56516 C * , P , K 20 4−32.05262.39273.316
Table A2. Fitted models for COD with their statistics.
Table A2. Fitted models for COD with their statistics.
Horizontal Flow Bed SS-HFVertical Flow Bed SS-VF
ModelVariablesklogLikAICc Δ i w i ModelVariablesklogLikAICc Δ i w i
1 T k , θ m , K 20 418.702−28.1930.000>0.9651 T k , θ m , K 20 417.528−25.8440.000>0.970
2 P , T k , θ m , K 20 518.823−25.7722.421 2 θ , T k , θ m , K 20 518.001−24.1271.716
3 θ , T k , θ m , K 20 518.711−25.5472.6453 P , T k , θ m , K 20 517.591−23.3072.537
4 C * , T k , θ m , K 20 518.702−25.5302.6634 C * , T k , θ m , K 20 517.528−23.1812.663
5 θ , P , T k , θ m , K 20 618.832−22.9555.2385 θ , P , T k , θ m , K 20 618.002−21.2954.549
6 C * , P , T k , θ m , K 20 618.823−22.9375.2566 C * , θ , T k , θ m , K 20 618.001−21.2934.551
7 C * , θ , T k , θ m , K 20 618.711−22.7135.4807 C * , P , T k , θ m , K 20 617.591−20.4725.372
8 C * , θ , P , T k , θ m , K 20 718.832−19.9318.2618 C * , θ , P , T k , θ m , K 20 718.002−18.2717.573
9 θ , K 20 314.117−21.5286.6643.449 × 10−29 θ , K 20 312.794−18.8836.9612.986 × 10−2
10 θ , P , K 20 414.821−20.4307.762 10 θ , P , K 20 413.497−17.7828.062
11 C * , θ , K 20 414.117−19.0229.17011 C * , θ , K 20 412.794−16.3769.468
12 C * , θ , P , K 20 514.821−17.76810.42512 C * , θ , P , K 20 513.497−15.11910.725
13 K 20 2−18.24440.83169.0239.922 × 10−1613 K 20 2−20.16544.67270.5164.726 × 10−16
14 C * , K 20 3−18.24443.19471.386 14 C * , K 20 3−20.16547.03572.879
15 P , K 20 3−18.24443.19471.38615 P , K 20 3−20.16547.03572.879
16 C * , P , K 20 4−18.24445.70073.89216 C * , P , K 20 4−20.16549.54275.385
Table A3. Fitted models for NH4+-N with their statistics.
Table A3. Fitted models for NH4+-N with their statistics.
Horizontal Flow Bed SS-HFVertical Flow Bed SS-VF
ModelVariablesklogLikAICc Δ i w i ModelVariablesklogLikAICc Δ i w i
1 T k , θ m , K 20 4−4.37117.9540.000>0.9801 T k , θ m , K 20 44.5830.0460.000>0.940
2 θ , T k , θ m , K 20 5−4.34220.5592.605 2 θ , T k , θ m , K 20 55.0601.7541.708
3 P , T k , θ m , K 20 5−4.37020.6152.6623 P , T k , θ m , K 20 54.8832.1092.063
4 C * , T k , θ m , K 20 5−4.37120.6172.6634 C * , T k , θ m , K 20 54.5832.7092.663
5 C * , θ , T k , θ m , K 20 6−4.34223.3945.4405 θ , P , T k , θ m , K 20 65.1174.4764.429
6 C * , P , T k , θ m , K 20 6−4.37023.4505.4976 C * , θ , T k , θ m , K 20 65.0604.5894.542
7 θ , P , T k , θ m , K 20 6−4.63823.9866.0327 C * , P , T k , θ m , K 20 64.8834.9444.897
8 C * , θ , P , T k , θ m , K 20 7−4.63827.0109.0568 C * , θ , P , T k , θ m , K 20 75.1177.4997.453
9 θ , K 20 3−9.50625.7197.7652.018 × 10−29 θ , 30.5665.5745.5275.932 × 10−2
10 θ , P , K 20 4−8.69226.5978.643 10 θ , P , K 20 41.7095.7935.747
11 C * , θ , K 20 4−9.50628.22510.27111 C * , θ , K 20 40.5668.0808.034
12 C * , θ , P , K 20 5−8.69229.26011.30612 C * , θ , P , K 20 51.7098.4568.410
13 K 20 2−37.86080.06262.1083.195 × 10−1413 K 20 2−33.10270.54770.5004.618 × 10−16
14 C * , K 20 3−37.86082.42564.471 14 C * , K 20 3−33.10272.91072.863
15 P , K 20 3−37.86082.42564.47115 P , K 20 3−33.10272.91072.863
16 C * , P , 4−37.86084.93166.97816 C * , P , K 20 4−33.10275.41675.370
Table A4. Fitted models for TN with their statistics.
Table A4. Fitted models for TN with their statistics.
Horizontal Flow Bed SS-HFVertical Flow Bed SS-VF
ModelVariablesklogLikAICc Δ i w i ModelVariablesklogLikAICc Δ i w i
1 T k , θ m , K 20 414.709−20.2060.000>0.8381 T k , θ m , K 20 45.926−2.6400.000>0.992
2 θ , T k , θ m , K 20 514.744−17.6132.593 2 θ , T k , θ m , K 20 56.384−0.8931.747
3 C * , T k , θ m , K 20 514.709−17.5432.6633 C * , T k , θ m , K 20 55.9260.0232.663
4 P , T k , θ m , K 20 514.709−17.5432.6634 P , T k , θ m , K 20 55.9260.0232.663
5 C * , θ , T k , θ m , K 20 614.744−14.7795.4285 C * , θ , T k , θ m , K 20 66.3841.9424.582
6 θ , T k , θ m , K 20 614.732−14.7555.4516 θ , P , T k , θ m , K 20 66.3841.9424.582
7 C * , P , T k , θ m , K 20 614.709−14.7095.4987 C * , P , T k , θ m , K 20 65.9262.8585.498
8 C * , θ , P , T k , θ m , K 20 714.732−11.7318.4758 C * , θ , P , T k , θ m , K 20 76.3844.9667.605
9 θ , K 20 311.810−16.9153.291>0.1619 θ , K 20 3−0.1827.0719.7107.729 × 10−3
10 θ , P , K 20 411.887−14.5625.644 10 θ , P , K 20 4−0.0849.38112.020
11 C * , θ , K 20 411.810−14.4095.79811 C * , θ , K 20 4−0.1829.57712.216
12 C * , θ , P , K 20 511.887−11.8998.30712 C * , θ , P , K 20 5−0.08412.04414.683
13 K 20 2−2.4059.15429.3603.532 × 10−713 K 20 2−9.12422.59025.2293.297 × 10−6
14 C * , K 20 3−2.40511.51731.723 14 C * , K 20 3−9.12424.95327.592
15 P , K 20 3−2.40511.51731.72315 P , K 20 3−9.12424.95327.592
16 C * , P , K 20 4−2.40514.02334.22916 C * , P , K 20 4−9.12427.45930.099
Figure A1. Diagnostic for BOD5 models. (A) Horizontal bed, model T k , θ m , K 20 ; (B) vertical bed, model T k , θ m , K 20 ; (C) horizontal bed, model θ , K 20 ; (D) vertical bed, model θ , K 20 .
Figure A1. Diagnostic for BOD5 models. (A) Horizontal bed, model T k , θ m , K 20 ; (B) vertical bed, model T k , θ m , K 20 ; (C) horizontal bed, model θ , K 20 ; (D) vertical bed, model θ , K 20 .
Water 11 00180 g0a1
Figure A2. Diagnostic for COD models. (A) Horizontal bed, model T k , θ m , K 20 ; (B) vertical bed, model T k , θ m , K 20 ; (C) horizontal bed, model θ , K 20 ; (D) vertical bed, model θ , K 20 .
Figure A2. Diagnostic for COD models. (A) Horizontal bed, model T k , θ m , K 20 ; (B) vertical bed, model T k , θ m , K 20 ; (C) horizontal bed, model θ , K 20 ; (D) vertical bed, model θ , K 20 .
Water 11 00180 g0a2
Figure A3. Diagnostic for NH4+-N models. (A) Horizontal bed, model T k , θ m , K 20 ; (B) vertical bed, model T k , θ m , K 20 ; (C) horizontal bed, model θ , K 20 ; (D) vertical bed, model θ , K 20 .
Figure A3. Diagnostic for NH4+-N models. (A) Horizontal bed, model T k , θ m , K 20 ; (B) vertical bed, model T k , θ m , K 20 ; (C) horizontal bed, model θ , K 20 ; (D) vertical bed, model θ , K 20 .
Water 11 00180 g0a3
Figure A4. Diagnostic for TN models. (A) Horizontal bed, model T k , θ m , K 20 ; (B) vertical bed, model T k , θ m , K 20 ; (C) horizontal bed, model θ , K 20 ; (D) vertical bed, model θ , K 20 .
Figure A4. Diagnostic for TN models. (A) Horizontal bed, model T k , θ m , K 20 ; (B) vertical bed, model T k , θ m , K 20 ; (C) horizontal bed, model θ , K 20 ; (D) vertical bed, model θ , K 20 .
Water 11 00180 g0a4

References

  1. Nowak, R.; Mechler, R.; Enderle, P. Evaluation of reject waters from co-digestion of solid wastes from agro-industries in municipal WWTP. Water Sci. Technol. 2007, 55, 37–44. [Google Scholar] [CrossRef] [PubMed]
  2. Dąbrowski, W.; Karolinczak, B. Application of subsurface vertical flow constructed wetlands to reject water treatment in dairy wastewater treatment plant. Environ. Technol. 2017, 38, 175–182. [Google Scholar] [CrossRef] [PubMed]
  3. Dąbrowski, W. The possibility of dairy WWTP load decreasing using constructed wetland for reject water treatment. Annu. Set Environ. Prot. 2013, 15, 901–913. [Google Scholar]
  4. Dąbrowski, W. A study of the digestion process of sewage sludge from dairy WWTP to determine the composition and load of reject water. Water Pract. Technol. 2014, 9, 71–78. [Google Scholar] [CrossRef]
  5. Janus, H.M.; van der Roest, H.F. Do not reject the idea of treating reject water. Water Sci. Technol. 1997, 35, 27–34. [Google Scholar] [CrossRef]
  6. Fux, C.; Boehler, M.; Huber, P.; Brunner, I.; Siegrist, H. Biological treatment of ammonium-rich wastewater by partial nitritation and subsequent anaerobic ammonium oxidation (anammox) in a pilot plant. J. Biotechnol. 2002, 99, 293–306. [Google Scholar] [CrossRef]
  7. Fux, C.; Lange, K.; Faesseler, A.; Huber, P.; Grueniger, B.; Siegrist, H. Nitrogen removal from digester supernatant via nitrite-SBR or SHARON? Water Sci. Technol. 2003, 48, 9–18. [Google Scholar] [CrossRef]
  8. Fux, C.; Valten, S.; Carozzi, V.; Solley, D.; Keller, J. Efficient and stable nitrification and denitrification of ammonium-rich sludge dewatering liquor using SBR with continuous loading. Water Res. 2006, 40, 2765–2775. [Google Scholar] [CrossRef]
  9. Kołecka, K.; Gajewska, M.; Obarska-Pempkowiak, H.; Rohde, D. Integrated dewatering and stabilization system as an environmentally friendly technology in sewage sludge management in Poland. Ecol. Eng. 2017, 98, 346–353. [Google Scholar] [CrossRef]
  10. Boruszko, D.; Butarewicz, A. Impact of effective microorganisms bacteria on low-input sewage sludge treatment. Environ. Prot. Eng. 2015, 41, 83–96. [Google Scholar]
  11. Kołecka, K.; Obarska-Pempkowiak, H. The quality of sewage sludge stabilized for a long time in reed basins. Environ. Prot. Eng. 2008, 34, 13–20. [Google Scholar]
  12. Kadlec, R.H.; Wallace, S.D. Treatment Wetlands, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2009. [Google Scholar]
  13. Karolinczak, B.; Dąbrowski, W. Effectiveness of septage pre-treatment in subsurface vertical flow constructed wetlands. Water Sci. Technol. 2017, 76, 2544–2553. [Google Scholar] [CrossRef] [PubMed]
  14. Kinsley, C.B.; Crolla, A.M.; Kuyucak, N.; Zimmer, M.; Lafléche, A. Nitrogen dynamics in a constructed wetland system treating landfill leachate. In Proceedings of the 10th International Conference on Wetland Systems for Water Pollution Control, Lisbon, Portugal, 23–29 September 2006; pp. 295–305. [Google Scholar]
  15. Obarska-Pempkowiak, H.; Kolecka, K.; Gajewska, M.; Wojciechowska, E.; Ostojski, A. A Sustainable Sewage Management in Rural Areas. Annu. Set Environ. Prot. 2015, 17, 585–602. [Google Scholar]
  16. Dąbrowski, W. Reject Water Treatment from Dairy WWTPs with Constructed Wetlands Systems; BUT Publishing House: Białystok, Poland, 2014. (In Polish) [Google Scholar]
  17. Kadlec, R.H. Large Constructed Wetlands for Phosphorus Control: A Review. Water 2016, 8, 243. [Google Scholar] [CrossRef]
  18. Cui, L.; Li, W.; Zhang, Y.; Wei, J.; Lei, Y.; Zhang, M.; Pan, X.; Zhao, X.; Li, K.; Ma, W. Nitrogen removal in a horizontal subsurface flow constructed wetland estimated using the first-order kinetic model. Water 2016, 8, 514. [Google Scholar] [CrossRef]
  19. Carvalho, P.N.; Arias, C.A.; Brix, H. Constructed Wetlands for Water Treatment: New Developments. Water 2017, 9, 397. [Google Scholar] [CrossRef]
  20. Langergraber, G. Applying Process-Based Models for Substrate Flow Treatment Wetlands: Recent Developments and Changes. Water 2017, 9, 5. [Google Scholar] [CrossRef]
  21. The Chancellery of the Sejm of the Republic of Poland. Regulations of the Minister of Environment from 18th of November 2014 on Conditions to be Met for Disposal of Treated Sewage into Water and Soil and Concerning Substances Harmful to the Environment (Dz.U. 2014. no. 1800); Chancellery of the Sejm—Poznaj Sejm: Warsaw, Poland, 2014. [Google Scholar]
  22. American Public Health Association (APHA). 2005 Standard Methods for Examination of Water and Wastewater, 21st ed.; American Public Health Association: Washington, DC, USA, 2005. [Google Scholar]
  23. Burnham, K.P.; Anderson, D.R. Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach, 2nd ed.; Springer: New York, NY, USA, 2002. [Google Scholar]
  24. Kadlec, R.H. The inadequacy of first-order treatment wetland models. Ecol. Eng. 2000, 15, 105–119. [Google Scholar] [CrossRef] [Green Version]
  25. Merriman, L.S.; Hathaway, J.M.; Burchell, M.R.; Hunt, W.F. Adapting the Relaxed Tanks-in-Series Model for Stormwater Wetland Water Quality Performance. Water 2017, 9, 691. [Google Scholar] [CrossRef]
  26. Struk-Sokołowska, J. Variability of dairy wastewater characteristics in Piatnica—One of the largest and most advanced milk processing plants in Poland. E3S Web Conf. 2018, 44, 00169. [Google Scholar] [CrossRef]
  27. Gajewska, M.; Obarska-Pempkowiak, H. Multistage treatment wetland for treatment of reject waters from digested sludge dewatering. Water Sci. Technol. 2013, 68, 1223–1232. [Google Scholar] [CrossRef] [PubMed]
  28. Struk-Sokołowska, J.; Rodziewicz, J.; Mielcarek, A. Effect of dairy wastewater on changes in COD fractions in technical-scale SBR type reactors. Water Sci. Technol. 2017, 156–169. [Google Scholar] [CrossRef] [PubMed]
  29. Struk-Sokołowska, J.; Mielcarek, A.; Wiater, J.; Rodziewicz, J. Impacts of dairy wastewater and pre-aeration on the performance of SBR treating municipal sewage. Desalin. Water Treat. 2018, 105, 41–50. [Google Scholar] [CrossRef]
  30. Nash, J.C. Compact Numerical Methods for Computers, Linear Algebra and Function Minimisation, 2nd ed.; Adam Hilger: Bristol, UK, 1990. [Google Scholar]
  31. Wedderburn, R.W.M. Quasi-likelihood Functions, Generalized Linear Models, and the Gauss-Newton Method. Biometrika 1974, 61, 439–447. [Google Scholar]
  32. Sugiura, N. Further analysis of the data by Akaike’s information criterion and the finite corrections. Commun. Stat. Theory Methods 1978, A7, 13–26. [Google Scholar] [CrossRef]
  33. Akaike, H. A new look at the statistical model identification. IEEE Trans. Autom. Control 1974, 19, 716–723. [Google Scholar] [CrossRef]
  34. Baty, F.; Ritz, C.; Charles, S.; Brutsche, M.; Flandrois, J.P.; Delignette-Muller, M.L. A Toolbox for Nonlinear Regression in R: The Package nlstools. J. Stat. Softw. 2015, 66, 1–21. [Google Scholar] [CrossRef]
  35. Muggeo, V.M.R. Segmented: An R Package to Fit Regression Models with Broken-Line Relationships. R News 2008, 8, 20–25. [Google Scholar]
  36. R Core Team. A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018; Available online: https://www.R-project.org/ (accessed on 1 December 2018).
  37. Nash, J.C.; Murdoch, D. nlsr: Functions for Nonlinear Least Squares Solutions, R Package Version 2018.1.28. Available online: https://CRAN,R-project,org/package=nlsr (accessed on 1 December 2018).
  38. Mazerolle, M.J. AICcmodavg: Model Selection and Multimodel Inference Based on (Q)AIC(c), R Package version 2.1-1. 2017. Available online: https://cran,r-project,org/package=AICcmodavg (accessed on 1 December 2018).
Figure 1. Scheme of research installation, with sampling points (I, II, III).
Figure 1. Scheme of research installation, with sampling points (I, II, III).
Water 11 00180 g001
Figure 2. Cross section of SS VF (subsurface vertical flow) and HF (subsurface horizontal flow beds) beds.
Figure 2. Cross section of SS VF (subsurface vertical flow) and HF (subsurface horizontal flow beds) beds.
Water 11 00180 g002
Figure 3. Selected model fit.
Figure 3. Selected model fit.
Water 11 00180 g003
Figure 4. TP removal efficiency and reject water temperature during the research period.
Figure 4. TP removal efficiency and reject water temperature during the research period.
Water 11 00180 g004
Figure 5. An application “Modeling of reject water treatment efficiency”.
Figure 5. An application “Modeling of reject water treatment efficiency”.
Water 11 00180 g005
Table 1. The characteristics of reject water before and after treatment in SS VF and SS HF beds (38 research series, 114 samples).
Table 1. The characteristics of reject water before and after treatment in SS VF and SS HF beds (38 research series, 114 samples).
Before Treatment (mg L−1)After Treatment (mg L−1)
Parameter SS-VFSS-HF
BOD5120.53 ± 16.9721.55 ± 8.2731.18 ± 10.04
COD201.92 ± 21.4840.82 ± 12.1449.79 ± 12.44
TOC53.18 ± 9.9110.18 ± 4.5218.08 ± 5.58
TSS141.01 ± 18.1220.60 ± 4.8426.88 ± 4.61
TN276.66 ± 42.71134.29 ± 19.44158.87 ± 20.29
NH4+-N195.53 ± 38.3250.71 ± 16.6073.29 ± 20.53
NO3-N0.69 ± 0.2812.68 ± 4.763.31 ± 1.33
NO2-N0.17 ± 0.070.13 ± 0.050.63 ± 0.20
TP22.17 ± 3.5416.43 ± 2.0013.95 ± 2.54
Mean ± standard deviation.
Table 2. Removed pollutant loads (mean values).
Table 2. Removed pollutant loads (mean values).
-Removed Load (g m−2 d−1)
ParameterSS-VFSS-HF
BOD59.908.93
COD16.1115.21
TOC4.303.51
TSS12.0011.41
TN14.2411.78
NH4+-N14.4812.22
TP0.570.82
Table 3. Parameters of selected models for BOD5, COD, NH4+-N and TN.
Table 3. Parameters of selected models for BOD5, COD, NH4+-N and TN.
Horizontal Flow Bed SS-HFVertical Flow Bed SS-VF
BOD5
ParameterEstimate/ValueLower 95% CIUpper 95% CIParameterEstimate/ValueLower 95% CIUpper 95% CI
T k 11.95910.30213.616 T k 12.58010.83614.323
θ m 1.1061.0731.140 θ m 1.0871.0611.113
K 20 0.1800.1700.190 K 20 0.2250.2130.238
RSS1.974-RSS2.327-
COD
ParameterEstimate/ValueLower 95% CIUpper 95% CIParameterEstimate/ValueLower 95% CIUpper 95% CI
T k 15.57813.29317.864 T k 17.00014.65619.344
θ m 1.0481.03501.06 θ m 1.0431.0331.052
K 20 0.1790.1720.187 K 20 0.2090.2000.218
RSS0.831-RSS0.884-
NH4+-N
ParameterEstimate/ValueLower 95% CIUpper 95% CIParameterEstimate/ValueLower 95% CIUpper 95% CI
T k 15,08712,8217,355 T k 16.00714.15317.862
θ m 1.1011.0731.131 θ m 1.0741.0581.089
K 20 0.1510.1390.164 K 20 0.1950.1840.207
RSS2.800-RSS1.748-
TN
ParameterEstimate/ValueLower 95% CIUpper 95% CIParameterEstimate/ValueLower 95% CIUpper 95% CI
T k 12.4019.13215.671 T k 10.3127.86012.763
θ m 1.0511.0241.078 θ m 1.0771.0331.122
K 20 0.0650.0600.069 K 20 0.0830.0770.089
RSS1.026-RSS1.629-

Share and Cite

MDPI and ACS Style

Dąbrowski, W.; Karolinczak, B.; Malinowski, P.; Boruszko, D. Modeling of Pollutants Removal in Subsurface Vertical Flow and Horizontal Flow Constructed Wetlands. Water 2019, 11, 180. https://doi.org/10.3390/w11010180

AMA Style

Dąbrowski W, Karolinczak B, Malinowski P, Boruszko D. Modeling of Pollutants Removal in Subsurface Vertical Flow and Horizontal Flow Constructed Wetlands. Water. 2019; 11(1):180. https://doi.org/10.3390/w11010180

Chicago/Turabian Style

Dąbrowski, Wojciech, Beata Karolinczak, Paweł Malinowski, and Dariusz Boruszko. 2019. "Modeling of Pollutants Removal in Subsurface Vertical Flow and Horizontal Flow Constructed Wetlands" Water 11, no. 1: 180. https://doi.org/10.3390/w11010180

APA Style

Dąbrowski, W., Karolinczak, B., Malinowski, P., & Boruszko, D. (2019). Modeling of Pollutants Removal in Subsurface Vertical Flow and Horizontal Flow Constructed Wetlands. Water, 11(1), 180. https://doi.org/10.3390/w11010180

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop