Proposing a Wetland-Based Economic Approach for Wastewater Treatment in Arid Regions as an Alternative Irrigation Water Source

: Point and nonpoint wastewater sources have a detrimental, negative effect on agriculture, soil, surface


Introduction
Egypt is suffering from a scarcity of water supplies and increased rates of reused wastewater because of rapid urbanization, population growth, and the depletion of fresh-Additionally, subsurface vertical flow constructed wetlands improve the nitrification capacity compared to those with a horizontal structure by enhancing oxygen transfer from the air to the wastewater flow [29]. However, the granular media is not continuously inundated since the pulsed water circulation is discontinuous. Wastewater treatment in constructed wetlands is a multi-step process that involves sedimentation, sorption, filtration, microbial decomposition, and plant uptake [30], each of which affects the others.
In both tropical and arid countries, constructed wetlands (CWs) are one method that has demonstrated high removal rates of nutrients and organic [31][32][33]. Numerous HSSF-CWs proposed, according to the research that has been conducted across the globe, offer treatments which meet the regional water quality [34,35]. The pollutants in wastewater are viewed from an economic and environmental perspective as externalities that need to be handled in order to appropriately internalize all of the costs and benefits of the new wastewater treatment technologies that are recommended. There is a significant amount of existing research on wastewater treatment costs which focus on the cost functions for treating wastewater. Some studies consider quality factors (such as the contaminants eliminated or the quality of the influent and effluent), while others merely consider the amount of wastewater that has been treated. While some studies approach the operation and maintenance costs by solely predicting the energy expenses, the others estimate all of the cost factors.
Based on sample data from 22 Spanish WWTPs, Molinos-Senante et al. [36] presented information on the overall operating and maintenance expenses for each of the following five categories: energy, staff, reagents, waste management, and maintenance. The staff, accounting for one-third of the overall costs, was determined to be the most crucial component. The next two most significant expenses are maintenance and energy prices at 21% and 18%, respectively. The expenses for waste management and reagents have comparable percentage weights, making up 15% and 14% of the overall costs, respectively. According to their research, the average cost of plants with nutrient removal systems is EUR 0.21/m 3 , with this dropping to EUR 0.18/m 3 if plants do not them. The parametric approach has been utilized in numerous different investigations to create cost functions for wastewater treatment. These are predicated on the development of a functional relationship between the cost drivers (explanatory variables) and the expenses of wastewater treatment (the dependent variable).
The objective of this study was to propose, design, and cost estimate a wastewater treatment system comprised of a primary sedimentation tank, an HSSF-CW followed by a VF-CW, and a ground storage tank. This system will serve small populations in developing areas in arid and semi-arid climates. Additionally, it will protect agricultural lands and groundwater from pollution, and the treated wastewater can be utilized in irrigation purposes. The system is located in El-Moghra Oasis in the western desert of Egypt in new land reclamation and cultivation project of 1.5 million acres, where the irrigation water is based on groundwater resources [37]. The system services 6000 inhabitants and has a wastewater discharge of 780 m 3 /d. To achieve the study objective, (i) we proposed a feasible design approach, and a cost estimation for a wastewater treatment system, (ii) we collected the data regarding the climate, wastewater discharge, and influent pollutants concentration in the study area, (iii) the P-k-C* relaxed tanks in series were utilized to design the proposed HSSF-CW and VF-CW to remove the BOD, TP, and FC pollutants in light of the treated wastewater regulations in Egypt, and (iv) a mathematical model to compute the HSSF and VF-CW areas and removal efficiency was constructed and solved using an Excel sheet based on the influent discharge (Q), the inlet pollutant concentration (Ci), the first-order areal rate coefficient (Ka), the background concentration (C*), the number of RTS (P), and C was the effluent target for the BOD, TP, and FC inlet concentration removal, and (v) a design and investment cost curve for a population range of between 500 and 9000 was created for the treatment system application in similar regions in order to shorten the design phases and cost estimation analysis.

The Study Area
The study area is located in the El-Moghra Oasis western desert of Egypt. It is a part of the land reclamation project with the aim to reclaim 1.5 million acres, as shown in Figure 1. The projected cultivated area has 170 thousand acres. The winter crops are wheat, barley, and green beans, while the summer crops are sugar beet and maize. In addition, there are dates and olive trees. The project is based on groundwater wells [37]. The Ministry of Water Resources and Irrigation set sustainability guidelines for the well spacing, well depth, daily water withdrawal capacity from each well, submersible pump depth, and identifying the energy source used to power the well (solar energy). In addition, the agricultural and municipal runoff water should not be dumped outside of the project's specified beneficiary areas, instead, it should be mixed with well water after the treatment processes. The climate data from the Dabaa station (western desert) between 1990 and 2020 were collected. and 9000 was created for the treatment system application in similar regions in orde shorten the design phases and cost estimation analysis.

The Study Area
The study area is located in the El-Moghra Oasis western desert of Egypt. It is a of the land reclamation project with the aim to reclaim 1.5 million acres, as shown in ure 1. The projected cultivated area has 170 thousand acres. The winter crops are wh barley, and green beans, while the summer crops are sugar beet and maize. In addit there are dates and olive trees. The project is based on groundwater wells [37]. The M istry of Water Resources and Irrigation set sustainability guidelines for the well spac well depth, daily water withdrawal capacity from each well, submersible pump de and identifying the energy source used to power the well (solar energy). In addition agricultural and municipal runoff water should not be dumped outside of the proje specified beneficiary areas, instead, it should be mixed with well water after the treatm processes. The climate data from the Dabaa station (western desert) between 1990 2020 were collected.

Water Quality Data
The proposed HSSF-CW is designed to service a population of 6000 people (farm in the village. A single well irrigated an area of around 230 acres, which is shared by farmers, each of whom lives in the village. There are often six family members living gether in one home. The community has a population of roughly 6000 people, and

Water Quality Data
The proposed HSSF-CW is designed to service a population of 6000 people (farmers) in the village. A single well irrigated an area of around 230 acres, which is shared by ten farmers, each of whom lives in the village. There are often six family members living together in one home. The community has a population of roughly 6000 people, and its territory serves about 100 homes. The influent flow rate, 780 m 3 /d, is implied by the unit wastewater flow of 0.13 m 3 /capita/d. The BOD concentrations in the raw sewage water of 50 g/capita/d indicate a BOD concentration of 300 mg/L. In addition, the TN, TP, and FC influent concentrations are 35 mg/L, 10 mg/L, and 10 7 Counts/100 mL, respectively. The effluent is made to comply with Egypt's criteria [39] for wastewater reuse. The BOD and FC effluent concentrations are limited to 30 mg/L and 5000 Counts/100 mL, respectively. It is essential to perform the primary treatment prior to the HSSF-CW in order to reduce the amount of easily decomposable organic material that would otherwise accumulate in the wetland system's entry zone, leading to blockages, odors, and adverse impacts on the entry zone's plants. To accomplish this, similar technologies such as ponds, conventional primary treatment, Imhoff tanks, septic tanks, etc., can be used [11]. The main issue causing a mixed flow to occur in the HSSF-CW is media clogging [30]. Table 1 summarizes the influent and effluent pollutants concentrations during the primary sedimentation stage.

HSSF-and VF-CW's Model Descrpition
In this study, the Continuous Stirred Tank Reactor (CSTR) model and the k-C* model are combined in the relaxed tanks in series (RTS) model (Equation (1)) to effectively explain both the pollutant depletion and the flow hydrodynamics [40][41][42]. This model takes into account steady state circumstances, which include no evapotranspiration, no infiltration, and a constant flow of the water in the wetland.
where C = effluent concentration, mg/L, C * = background concentration, mg/L, C i = influent concentration, mg/L, P = apparent number of RTS, dimensionless, q = hydraulic loading rate, m/d, and Ka = first-order areal rate coefficient, which is given by: K 20 is the removal rate constant at 20 • C (m/year), θ is the dimensionless temperature coefficient, and T is wastewater temperature in • C. In addition, the hydraulic residence time (RT) in days is given by: where Q is the design discharge (m 3 /d); A is the surface area of the system (m 2 ); V is the system volume (m 3 ); y is the depth of flow (m); ϕ is the media porosity, which expresses the space available for the water to flow through the media, roots, and other solids in the HSSF-CW. In addition, the constructed wetland treatment efficiency (WE) is given as follows: In the subsurface engineered wetland systems, Reed et al. [42] suggested that the organic loading should not go over a limit value of 10 g/m 2 /d. This relationship can be expressed as follows:

Proposed HSSF and VF-CWs Construction Details
The HSSF-CW was selected as this kind of system does not have a clogging problem. The outlet flow comes from the primary sedimentation tank entering the HSSF-CW system (cell). The cell design comprises a rectangular bed that is surrounded by a 0.25 m masonry construction, the bed is lined by 10 cm of plain concrete or it is lined by geomembranes such as PVC sheets, which are durable, economical cost, and easily constructed, to prevent wastewater seepage. In this study, geomembrane lining was applied. Since the cell is 0.30 m above and 0.30 m below the earth, no groundwater from the natural ground surface can enter the cell. The cell in constructed wetland has a length-to-width ratio from 1:1 to 3:1. The wetland plants are Phragmites (Reeds) and Papyrus, the plant count is 9 per m 2 . Figure 2A,B shows the details of the reed bed media to avoid blocking in the HSSF and VF-CWs, respectively [30]. The thickness of the bed media (gravel and sand) for the HSSF-CW is 60 cm, and it is 75 cm for the V-CW. Therefore, the HSSF-CW media bed consists of four zones: an inlet zone with a gravel size of 4-8 cm, zone 1 with a gravel size of 1.6-3.2 cm, zone 2 (treatment zone) with a gravel size of 0.8-0.12 cm, and finally, an outlet zone with a gravel size of 4-8 cm. On the other hand, the VF-CW media bed comprises of a top zone (inlet) with a gravel size of 1.6-5.2 cm, which is followed by the treatment layer with a sand size of 0.5-4 mm, and this sand is supported by a gravel layer which has a size of 0.8-1.2 cm, and finally, the bed layer has a gravel of size 1.6-3.2 mm. In addition, the bed slope for the HSSF and VF-CWs is adjusted by 1-2% so that wastewater would flow towards the effluent end. During the preparation stage, the water was added to the HSSF and VF-CWs to help the soil settle and encourage the growth of macrophytes.

Modeling HSSF-CW and VF-CW Efficiencies Using the P-K-C*
The P-K-C* model, which is based on the relaxed tanks in series (RTS) flow model, is applied to design the HSSF-VF-CWs. Equation (1) is utilized to lower the inflow pollutants concentration to the desired level according to the wastewater treatment regulations in Egypt [39]. For the estimation of the area, the following variables are used: influent discharge (Q), the inlet concentration (C i ), the first-order areal rate coefficient (Ka), the background concentration (C*), and the number of RTS (P), and C is the effluent target for BOD, TP, and FC. The design discharge (Q) value is 780 m 3 /d, and the BOD, TP, and FC influent concentrations after primary sedimentation are 210 mg/L, 27 mg/L, and 10 6 Counts/100 mL, respectively. In this study, the background pollutant concentrations and the removal rate constant at 20 • C were taken from similar HSSF-CWs which were developed by Rashed [43]. Therefore, for the BOD, TP, and FC, the background pollutant concentrations (C * ) are 1 mg/L, 0.119 mg/L, and 4 Counts/100 mL, respectively. Additionally, the Ka values for BOD, TP, and FC are 0.662, 0.16, and 1.492 m/d, respectively. The water depth y = 0.6 m for the HSSF-CW and 0.75 m for the VF-CW, the media porosity ϕ = 0.3, the number of RTS, P = 3, and we apply Phragmites (Reeds) and Papyrus plants. A flow control device is used to divert the flow of the water. The CW sizing criteria are based on Egypt's regulations guideline of no more than 60 mg/L BOD effluent [39]. Table 2 summarizes the model input data for the proposed wastewater treatment system for the VF-CW and the HSSF-CW. In addition, the flow chart in Figure 3 shows the model steps to size the HSSF and VF-CWs and obtain the pollutants removal efficacies. The following steps are employed using an Excel spreadsheet:

1.
Use the input model data summarized in Table 2 and obtain the HSSF and VF-CWs areas for the BOD, TP, and FC by utilizing Equation (1).

3.
The design area is the area corresponding to the maximum RT in days.

5.
Obtain the length and width in (m) for the HSSF and VF-CWs based on an aspect length-to-width ratio of 1:1 6.
Divide the CW width into n cells based on a cell width of 8.0 m, and by doing so, the number of cells (n) = width/8.0.

Proposed HSSF and VF-CWs Construction Details
The HSSF-CW was selected as this kind of system does not have a clogging problem. The outlet flow comes from the primary sedimentation tank entering the HSSF-CW system (cell). The cell design comprises a rectangular bed that is surrounded by a 0.25 m masonry construction, the bed is lined by 10 cm of plain concrete or it is lined by geomembranes such as PVC sheets, which are durable, economical cost, and easily constructed, to prevent wastewater seepage. In this study, geomembrane lining was applied. Since the cell is 0.30 m above and 0.30 m below the earth, no groundwater from the natural ground surface can enter the cell. The cell in constructed wetland has a length-to-width ratio from 1:1 to 3:1. The wetland plants are Phragmites (Reeds) and Papyrus, the plant count is 9 per m 2 . Figure 2 A, B shows the details of the reed bed media to avoid blocking in the HSSF and VF-CWs, respectively [30]. The thickness of the bed media (gravel and sand) for the

Modeling HSSF-CW and VF-CW Efficiencies Using the P-K-C*
The P-K-C* model, which is based on the relaxed tanks in series (RTS) flow model, is applied to design the HSSF-VF-CWs. Equation (1) is utilized to lower the inflow pollutants concentration to the desired level according to the wastewater treatment regulations in  6. Divide the CW width into n cells based on a cell width of 8.0 m, and by doing so, the number of cells (n) = width/8.0.

The Cost Function for Treating Wastewater
Numerous studies have examined the cost function of treating wastewater. Andreia et al. [44] and Molinos-Senante et al. [45] evaluated the costs of small, decentralized, energy-efficient wastewater treatment systems. They concentrated on secondary treatment procedures for small agglomerations (a population equivalent (P.E.) of 2000). The suggested cost function (in EUR per equivalent population) is for the investment and operating expenses of secondary wastewater treatment technology. Therefore, the cost for investment for constructed wetlands (Equation (5)) is: Wetland operation and maintenance cost, where x is the P.E., and y is the total cost expressed as EUR/P.E. To ascertain the connection between the examined process's inputs and outputs, the methodological foundation of the shadow prices is built on econometrics and optimization procedures [46]. Gkika et al. [47] conducted a design database of nine constructed wetlands (CW) facilities of various capacities (population equivalent (P.E.)), which was used to estimate the construction and operation costs, and then, to derive empirical equations relating to the required facility land area and the construction cost per the P.E. The total construction cost C (EUR 10 3 ) is based on the contribution costs of the (i) inlet works, screening, (ii) Imhoff, settling, and sludge tank, (iii) first stage VF-CWs, (iv) second stage VF-CWs, (v) third stage HSSF or VF-CWs, (vi) sludge treatment bed, (vii) ancillary works (pipe network, siphon), (viii) disinfection (chlorination tank or maturation pond) (%), (ix) electrical, and (x) infrastructure and environmental restoration.
With a determination factor of R 2 = 0.97 and the total cost area (A t ) in ha: With a determination factor of R 2 = 0.97, the annual operating expense is EUR 17,853.74 or EUR 29.75 per P.E. or EUR 0.46 per m 3 . This operation cost consists of the following: salary for the facility operator (one person working part-time), electric energy costs for the pump operation for wastewater lifting and moving from one stage to another (e.g., all of the facilities), and the lighting for the facilities. For the computation, they used the values of EUR 127 for the electric energy per 1000 kWh and EUR 8 for operator hours.

HSSF and VF-CWs Removal Efficiencies
Using an Excel spreadsheet to solve Equations (1)-(3) utilizing the input data summarized in Table 2, the results show that the proposed HSSF-CW design treatment surface areas for the BOD, TP, and FC are 2198, 2375, and 5712 m 2 , respectively ( Table 3). The results show the area needed for the FC treatment is almost twice as large as the space needed for the TP therapy (2375 m 2 ) (Table 3). Therefore, for economic purposes, we selected the design area of the HSSF-CW is 2375 m 2 , with q = 0.33 m/d and the RT being 0.55 d. Consequently, the removal efficiencies in the HSSF-CW for the BOD, TP, and FC based on the design area (2375 m 2 ) are 78.2, 35.7, and 93.7%, respectively, as shown in Table 3. This indicates an increase in the removal efficiency from 76.2 to 78.2%, on the other hand, FC efficiency decreased from 99% to 93.7%. Therefore, in order to increase the FC removal efficiency, the outlet distribution channel is designed to be 4 m wide to utilize the solar hours in the disinfection of the FC [10,29,48]. Additionally, the proposed VF-CW design treatment surface area for the BOD, TP, and FC indicates areas of 2648, 2193, and 2700 m 2 , respectively (Table 3), whereas the water hydraulic residence time (RT) is 0.77, 0.63, and 0.78 d for the removal efficiencies of 80, 33.3, and 95% for the BOD, TP, and FC, respectively. Therefore, for economic purposes, as previously in the HSSF-CW, the selected design area of the VF-CW is 2193 m 2 , with q = 0.36 m/d and the RT being 0.63 d. The removal efficiencies in the VF-CW for the BOD, TP, and FC based on the design area (2193 m 2 ) are 75, 35.7, and 92.7%, respectively. In addition, Table 4 shows the overall cumulative efficiencies of the treatment system, which are 96.7, 70, and 100% for the BOD, TP, and FC, respectively.

HSSF and VF-CWs Sizing Design Curves
To facilitate design steps for the HSSF and VF-CWs, design curves for a population range from 500 to 9000 (discharge range 65-1170 m 3 /d) have been established utilizing an Excel spreadsheet. The design curves comprise (i) the HSSF-CW and VF-CW discharge and area function; (ii) the discharge, number of cells, and length of cell functions; (iii) the investment cost (LE) per capita function; (iv) the total investment cost function. To obtain the wetland design area for the HSSF and VF-CWs over a discharge range from 65 to 1170 m 3 /d, the influent and effluent pollutants concentrations from Table 2 were substituted into Equations (1)-(3) using an Excel sheet.
The design curves in Figure 4 are based on the BOD, TP, and FC removal efficiencies of 78.2%, 35.7%, and 93.7%, respectively, for the HSSF-CW. Additionally, the VF-CW design curves are based on the removal efficiencies for the BOD, TP, and FC of 75, 33.3, and 92.7%, respectively. The results of HSSF-CW design dimensions are listed in Table 5, and they include a discharge range of 65-1170 m 3 /d. In addition, the VF-CW design dimensions are listed in Table 6, and they include a P.E. range from 500 to 9000 and a discharge range of 65-1170 m 3 /d.
In addition, Figure 5 shows a detailed drawings of HSSF-CW and VF-CW, (A) cross-section elevation and (B) plan. Therefore, the system is intended for a 780 m 3 /d flow rate, the constructed wetland has a length-to-width ratio of 1:1, and it is designed to be a plug flow reactor. The HSSF-CW indicates a rectangular bed of area 2375 m 2 , which is divided into six cells, each cell is 48.7 m long and 8 m in width, and they house the common wetland plants Phragmites (Reeds) and Papyrus, which are the design element of this HSSF-CW.     Therefore, for a population equivalent of 6000 capita and a discharge value of 780 m 3 /d, the HSSF-CW design area is 2375 m 2 ( Figure 4A); this area is divided into six cells ( Figure 4B), and each cell is 8 m in width and 48.7 m in length, indicating a 1:1 rectangularity ratio. Consequently, the discharge (m 3 /d) and the HSSF-CW design area relationship function (linear equation) ( Figure 4A) based on Table 4 are obtained by using Equation (7). y = 3.0449 x + 0.2288 (7) where y is the design area (m 2 ), and x is the discharge (m 3 /d), with a determination coefficient of R 2 = 0.99. In addition, a relationship between the discharge and the number of cells based on a HSSF-CW aspect ratio of 1:1 is obtained by using Equation (8).

Width of Cell (B) (m)
where y is the number of cells in (m), and x is the discharge (m 3 /d), with a determination coefficient of R 2 = 0.99. Additionally, a relationship between the discharge in m 3 /d and the cell length based on a rectangularity of 1:1 is obtained by using a power function in Equation (9): y = 1.7514 x 0.4994 (9) where y is the HSSF-WC length of the cell, and x is the discharge (m 3 /d), with a determination coefficient of R 2 = 0.99. The results show that a design curve for the VF-CW for the discharge (m 3 /d) and the design area relationship function (linear equation) ( Figure 4C) indicates: y = 2.8002 x + 9.2157 (10) where y is the VF-WC design area (m 2 ), and x is the discharge (m 3 /d), with a determination coefficient of R 2 = 0.99. In addition, a relationship between the discharge (m 3 /d) and the number of cells based on a rectangularity ratio of 1:1 is obtained by using a power function ( Figure 4D) that is as the same that in Equation (9). Additionally, a relationship between the discharge and the cell length based on a rectangularity ratio of 1:1 is obtained by using a power function ( Figure 4D) such as: (11) where y is the length of cell, and x is the discharge (m 3 /d), with a determination coefficient of R 2 = 0.99. Figure 4 shows   In addition, the VF-CW has a rectangular bed of area 2193 m 2 , which is divided into six cells, and each cell 46.8 m long and 8 m in width, and they house the common wetland plants Phragmites (Reeds) and Papyrus.

HSSF and VF-CWs Water Balance
Tables 7 and 8 summarize the water balance for the HSSF and VF-CWs, respectively. The monthly ETo and rainfall values were computed by the FAO-CROPWAT 8 model and its attached CLIMWAT 2 model to obtain the daily average climatic parameters at the Dabaa station (western desert), Egypt, from 1991 to 2020.  The soil in the study region is sandy, and the reed beds have a PVC membrane covering them (zero infiltration). In terms of the water balance (Equation (7)) yields for the HSSF-CW, the inflow discharge is 780 m 3 /d, whereas the outflow discharge ranges from 746.7 to 765.2 m 3 /d, and there is a monthly water loss of 1.9 to 4.3%. In addition, in terms of the water balance yields for the VF-CW, the inflow discharge is 780 m 3 /d, whereas the outflow discharge ranges from 746.9 to 766.3 m 3 /d, and the monthly water loss ranges from 1.8 to 4.2%. Therefore, the total monthly water loss for the HF and VF-CWs indicates a range from 3.7 to 8.5%.

HSSF and VF-CWs Efficiencies
By comparing the results for the BOD, TP, and FC removal efficiencies in the HSSF-CW (78.2, 35.7, and 93.7%, respectively), the hydraulic residence time is 0.55 d, and the surface loading rate is 0.33 m/d according to previous work by the authors of [30,43,49,50], indicating that these are a good wetland operation condition. Table 9 summarizes this comparison for the proposed HSSF and VF-CWs and their predicted removal efficiency with the monitored case studies. Therefore, for the Agaa HSSF wastewater treatment in the delta of Egypt, the BOD removal efficiency is 76%, which is very close to those of the proposed HSSF-CW (76%). By comparing the removal efficiency results for the HSSF-CW with those of Diana et al. [51], they studied the efficiency of subsurface constructed wetland for domestic wastewater treatment in Jakarta, Indonesia, with an average BOD of 205.08 mg/L and a wastewater discharge rate of 0.79 m 3 /d. The determined that the wetland's organic loading rate was 283.53 kg/ha/d, q was 0.20 m 3 /m 2 /d, and Typha latifolia plants were used. The BOD pollutant degradation constant was 0.76/day, and the COD constant was 0.73/day. The retention time in the built wetlands (RT) was 1.6 days. It used 0.6 m 2 of created wetland per person. The removal efficiencies for the BOD and TP were 75.56 ± 7.06 and 85.26 ± 11.35, respectively. Compared to this study, the removal efficiencies for the BOD and TP are 75% and 33%, respectively, and the TP indicated low values in our study. The Ka values in our study were 0.662 and 0.16 for the BOD and TP, respectively. Therefore, the Ka value of the newly constructed wetland was calculated using a scale, the standards of which were determined by a shorter deteriorating period. In addition, the sun's rays, wastewater temperature, amount of precipitation, and humidity all would have an impact on how quickly the things degraded. In the tropical climate, degradation would occur more quickly during the fieldwork in these sites. The pace of the degradation would have an effect on the retention time. In addition, modeling the HSSF and VF using the P-k-C* model reported a compatible result for the BOD, TP, TN, and FC compared to those in [15][16][17]44].
How accurately with which the variables are characterized, particularly the first order removal rate coefficients, which are influenced by the temperature and the time of year in which inlet pollutant concentrations vary could, however, have an impact on the accuracy. The temperature of the wastewater is another factor that influences the Ka values. The primary force behind the biological reactions, thereby altering the reaction constants, is temperature. Therefore, the efficiency of the CW treatment depends on temperaturesensitive chemical and physical processes.
Climate, soil composition, topography, biota, and time are the main factors influencing the decomposition of organic matter. Organic matter degradation occurs 57% more quickly in tropical environments than it does in sub-tropical environments [52]. The plants' contributions to the constructed wetland's physical effects on the root system include the roots serving as a home for microorganisms, plant absorption, evapotranspiration, other plant activities in the CW, and plant production is just another of their roles [53]. For nutrients such as nitrogen and phosphate, as well as organics, the vegetation primarily has a positive impact, supporting greater treatment efficiency. Phragmites australis (Common reed) is unquestionably the most widely used plant in the world. The species belonging to the genus Typha (latifolia, angustifolia, domingensis, orientalis, and glauca) and Scirpus (such as lacustris, validus, californicus, and acutus) are additional often employed species. For wastewater treatment, the plants used in the HSSF and VF-CWs should have the following characteristics: (1) rich belowground organs (such as rhizomes and roots) to save attached bacteria and oxygen for the substrate in areas adjacent to the roots and rhizomes; (2) be accepting of high organic and nutrient loads; (3) have a high aboveground biomass for winter insulation in cold and temperate regions and for nutrient removal via harvesting [54].
Comparing the dimensions and the efficiency of the proposed HSSF-CW and VF-CW with those in the literature indicates a good wetland operation condition. Reduced pollution was mostly caused by the plants. Utilizing plants served four key purposes including the filtration of suspended particles, bacterial development, oxygen uptake by the plant roots, and substrate maintenance [55]. The soil in the constructed wetlands played a role in lowering pollution. Sand in the media has the impact of shorter retention times (from 0.75 to 1.5 days) [56]. Pollutant removal in sandy soils occurred more quickly than it did in fine soils, while fine soils have a higher efficiency when they are compared to that of sandy soils. Additionally, fertilizer delivery to the substrate prior to loading it on the created wetland served to achieve maximum plant growth in anticipation of the absorption of pollutants. For plants to grow gradually, it is stated that the topsoil should be at least 15 cm deep and contain 5% organic content [57].

Cost Estimation Analysis
Based on an examination of the CW costs and benefits, the feasibility of the HSSF-CW has been investigated. The costs have been determined by using current mark costs in Egypt in Egyptian pounds (EGP) (EUR 1 = 20 EGP ). A quantities list and the market cost of civil work included in the HSSF and VF-CWs structures were determined. Therefore, the total investment cost for the HSSF-CW is EGP 1,959,081.5, and the investment cost per capita is EGP 326.5. Additionally, for the VF-CW, the total investment cost is EGP 977,006.6, and the investment cost per capita is EGP 162.8. To facilitate estimation of the investment cost for the HSSF and VF-CWs, design curves have been established for a population range from 500 to 9000 utilizing an Excel spreadsheet. Figure 6 shows the cost estimation design curves for the HSSF-CW for the investment cost (EGP) per capita function ( Figure 6A) and the total investment cost (EGP) ( Figure 6B) as indicated in Equations (12) and (13), respectively.
where y is the investment cost (EGP) per capita, and x is the equivalent population (P.E.), with a determination coefficient of R 2 = 0.984. y = 3009.1 x 0.7447 (13) where y is total investment cost (EGP), and x is the equivalent population (P.E.), with a determination coefficient of R 2 = 1. In addition, for the VF-CW, the cost estimation functions are shown in Figure 6, which indicates the investment cost (EGP) per capita and the total investment cost indicated by Equation (14) and Equation (15), respectively: y = 384.82 x 0.8901 (14) where y is the investment cost (EGP) per capita, and x is the equivalent population (P.E.), with a determination coefficient of R 2 = 0.1.
where y is the total investment cost (EGP), and x is the equivalent population (P.E.), with a determination coefficient of R 2 = 0.

Conclusions
This study proposed, designed, and assessed the cost estimation of a wastewater treatment system composed of a sedimentation tank, a horizontal subsurface flow constructed wetland (HSSF-CW), a vertical subsurface flow constructed wetland (VF-CW), and a storage tank. The treatment system can be used as a cheap construction, maintenance, and operational option for wastewater treatment in developing countries with small populations that live in dry and semi-arid climates. The system services 6000 capita, and it has a wastewater discharge value of 780 m 3 /d in the developing arid region in El-Moghra Oasis western desert of Egypt, where the 1.5 million acres used for the land reclamation project are based on groundwater irrigation. The average winter and summer temperatures are 14.2 °C and 24.4 °C, respectively. The P-k-C* relaxed tanks in series are utilized in the design of the proposed HSSF and VF-CWs to remove the BOD, TP, and FC water pollutants. The relaxed P-k-C* tanks in series equations are modeled and solved using an Excel spreadsheet for the HSSF and VF-CWs sizing and pollutants removal efficacies.
The results show that: (i) For the HSSF-CW, the design treatment surface area for the BOD, TP, and FC is 2375 m 2 , which corresponds to a hydraulic retention time (RT) of 0.55 d, a surface loading rate (q) of 0.33 m/d, and the removal efficiencies for the BOD, TP, and FC were 78.2, 35.7, and 93.7%, respectively. (ii) For the VF-CW, the findings indicated that the design treatment surface area for the BOD, TP, and FC for the HSSF-CW is 2193 m 2 , which corresponds to an RT of 0.63 d, a q of 0.36 m/d, and the removal efficiencies for the BOD, TP, and FC were 75, 35.7, and 92.7%, respectively. (iii) The overall cumulative efficiencies of the treatment system were 96.7, 70, and 100 for the BOD, TP, and FC, respectively. (iv) To facilitate the design steps for the HSSF and VF-CWs, a design curves functions for a population range from 500 to 9000 demonstrated a discharge range of 65-1170 m 3 /d, and the wetland area, number of cells, and cell dimensions have been established.
(v) Detailed drawings for the HSSF and VF-CWs were produced, and the cost estimate function curves were carried out. (vi) The total investment cost analysis for the proposed system corresponding to a wastewater discharge of 780 m 3 /d indicates a total investment cost of EUR 146,804 and EUR 24.46/P.E. (vii) The total monthly water loss for the HF and VF-CWs indicates a range from 3.7 to 8.5%. (vii) The treated wastewater can be used for irrigation practices.

Conclusions
This study proposed, designed, and assessed the cost estimation of a wastewater treatment system composed of a sedimentation tank, a horizontal subsurface flow constructed wetland (HSSF-CW), a vertical subsurface flow constructed wetland (VF-CW), and a storage tank. The treatment system can be used as a cheap construction, maintenance, and operational option for wastewater treatment in developing countries with small populations that live in dry and semi-arid climates. The system services 6000 capita, and it has a wastewater discharge value of 780 m 3 /d in the developing arid region in El-Moghra Oasis western desert of Egypt, where the 1.5 million acres used for the land reclamation project are based on groundwater irrigation. The average winter and summer temperatures are 14.2 • C and 24.4 • C, respectively. The P-k-C* relaxed tanks in series are utilized in the design of the proposed HSSF and VF-CWs to remove the BOD, TP, and FC water pollutants. The relaxed P-k-C* tanks in series equations are modeled and solved using an Excel spreadsheet for the HSSF and VF-CWs sizing and pollutants removal efficacies.
The results show that: (i) For the HSSF-CW, the design treatment surface area for the BOD, TP, and FC is 2375 m 2 , which corresponds to a hydraulic retention time (RT) of 0.55 d, a surface loading rate (q) of 0.33 m/d, and the removal efficiencies for the BOD, TP, and FC were 78.2, 35.7, and 93.7%, respectively. (ii) For the VF-CW, the findings indicated that the design treatment surface area for the BOD, TP, and FC for the HSSF-CW is 2193 m 2 , which corresponds to an RT of 0.63 d, a q of 0.36 m/d, and the removal efficiencies for the BOD, TP, and FC were 75, 35.7, and 92.7%, respectively. (iii) The overall cumulative efficiencies of the treatment system were 96.7, 70, and 100 for the BOD, TP, and FC, respectively. (iv) To facilitate the design steps for the HSSF and VF-CWs, a design curves functions for a population range from 500 to 9000 demonstrated a discharge range of 65-1170 m 3 /d, and the wetland area, number of cells, and cell dimensions have been established. (v) Detailed drawings for the HSSF and VF-CWs were produced, and the cost estimate function curves were carried out. (vi) The total investment cost analysis for the proposed system corresponding to a wastewater discharge of 780 m 3 /d indicates a total investment cost of EUR 146,804 and EUR 24.46/P.E. (vii) The total monthly water loss for the HF and VF-CWs indicates a range from 3.7 to 8.5%. (vii) The treated wastewater can be used for irrigation practices.
The findings of this study as the proposed wastewater treatment system compared with those of the traditional wastewater treatment demonstrate that it is a low-cost operation and maintenance option. Consequently, this strategy can be employed by the decision makers for water resources management in the Middle East and the Mediterranean region to enhance water quality according to the social and economic criteria. This study has resulted in the successful implementation of ecological restoration projects as an affordable treatment system and the addition of an unconventional water source that can be used for irrigation.