Modeling of novel processes for eliminating sidestreams impacts on full-scale sewage treatment plant using GPS-X7

The novel process consisted of two steps was established by combining all sidestreams lines (supernatant gravity thickener, underflow mechanical thickener, and centrate), treating them together away from the mainstream treatment plant, and returning treated sidestreams effluents to the plant outfall instead of plant head. The two steps novelty treatment combined degradation, nitrification, and dilution processes. To treat combined sidestreams, a novel pilot extended nutrient moving bed biofilm reactor was developed. The effects of sidestream elimination on a full-scale anaerobic/anoxic/oxic system were simulated using GPS-X7. The statistical results of R values greater than 0.8 and NMSE values near zero proved the calibrated model’s validation. The novel system successfully removed 98, 93, 100, 85, 98, 100, and 98% of BOD, COD, NH4, NO3, TSS, H2S, and PO4-P from sidestreams, respectively. Furthermore, the simulation results showed that eliminating sidestreams has reduced volumes of full-scale A2/O facilities, controlled hydraulic and pollutants shocks, and minimized cost and energy. The novel process proved successful in treating combined sidestreams and eliminating their impacts on the A/O2 system.

The activated sludge process (ASP) is the main biological system used in wastewater treatment plants (WWTPs) worldwide 1 . As part of ASP, the anaerobic/anoxic/oxic (A 2 /O) process is the most commonly utilized application 2 . One of the most critical challenges of the ASP system is sludge generation that requires treatment. Sludge treatment results in rejected water called sidestream (SS) 3 . The SSs in the A 2 /O system are three main lines: supernatant gravity thickener, underflow mechanical thickener, and centrate; all of these lines are returned to the head of the plant 4 . Although the quantity of the SSs ranges from 3 to 7% of the mainstream, the SSs contain high pollutants concentrations, especially nutrients, with about 30 and 10-80% nitrogen and phosphorous of the total load entering the plant, respectively. These high concentrations cause a tremendous burden on the plant's performance, especially for removing nutrients 5 . In the A 2 /O system, the return of SSs to the head of the plant creates hydraulic and mass shocks and operational problems. Most of the previous studies concerned the centrate line due to its high nutrient concentrations, where research was conducted to treat nutrients in this line before returning it to the head of the plant 6 .
Recently, several reactors have been developed to remove nitrogen using partial nitritation processes. Sharon's operations are the first of their kind in treating digested water 7 . After that, many new methods were developed under partial nitrification, such as anaerobic ammonium oxidation (ANAMMOX), DEamMONification (DEMON) process, completely autotrophic nitrogen-removal over nitrate (CANON) process, and others 8,9 . However, these reactors only reduced the nitrogen load without addressing other pollutants. No research has been done on treating combined SS lines away from the mainstream treatment system 10 .
Therefore, the activated sludge systems require an improvement in the treatment process to overcome its disadvantages and cope with restricted discharge standards 11 . However, the ASP is quite complex, considering the multivariable structures and multiple time scales in the internal process dynamics 12 . Designed wastewater treatment plants need effective optimization processes to obtain higher pollutants removal that is extremely time-consuming and has great uncertainty even for the most experienced engineers 13 . This uncertainty may be overcome by simulating process modification to get initial results for further scale-up process 14 .
The model simulation provides a valuable and effective tool to optimize complicated biological treatment processes. It has been recognized as an indispensable tool in redesigning and managing existing wastewater treatment plants 15 . Numerous commercial numerical simulation software (e.g., Simba, TOXCHEM, GPS-X, BioWin, and WEST) have incorporated activated sludge models (e.g., ASM1, ASM2, ASM2d, and ASM3) have been developed for engineering practices 16,17 . GPS-X Model is one of the best models for simulating wastewater treatment plants. This model performs many operations and functions and has high reliability in simulating treatment plants 18 . For process optimization, Andres, et al. 19 used GPS-X to simulate a complete mixed batch reactor and optimize the number of batch runs required to treat high COD concentration wastewater. The finding proved that the calibrated model gave an accurate prediction that simulated the actual results. El-Hoz Karbala WWTP mainstream and sidestreams quality characteristics. This study collected quality characteristics from the Karbala WWTP influent after the grit chamber and the effluent after the chlorination basin. The collected samples were transferred to the Karbala Sewer Directorate Laboratory. They immediately analyzed for pH, dissolved oxygen, Chemical Oxygen Demand (COD), Biochemical Oxygen Demand (BOD 5 ), NO 2 -N, NO 3 --N, NH 3 -N, PO 4 -P, H 2 S, and Total phosphorus (TP) according to procedures in "The Standard Methods for the Examination of Water and Wastewater" 21 . In this investigation, SSs from three lines were combined and treated separately from the mainstream treatment system: supernatant gravity thickener, underflow belt thickener, and reject water (centrate). Samples were collected for 12 months from the mainstream and SSs of Karbala WWTP. Table 2 shows the performance of the Karbala sewage treatment plant, while Table 3 shows the concentrations of pollutants in the sidestream lines.
Novel pilot EN-MBBR system configuration. A novel pilot has been designed and manufactured with an attached bacterial growth system, an upgraded system for a moving bed biofilm reactor (MBBR), as shown in Fig. 2. It is designed to treat a 30 m 3 /day SSs discharge with projected organic matter and nutrients removal efficiency up to 95%. The system contains three parts: the first part is designed to treat organic matter by heterotrophic bacteria, the second part is intended to treat nutrients by the growth of autotrophic bacteria, and the    www.nature.com/scientificreports/ third part is a sedimentation basin to remove suspended solids. The system's novelty was in the second part by increasing the hydraulic retention time (HRT) to 13 h to obtain complete nitrification and endogenous respiration, while HRT in the first part was 3.3 h. The aeration system has 2 air blowers and 16 diffusers to provide 65% of the dissolved oxygen (DO) required to oxidize pollutants and 35% of the amount of DO in the underflow mechanical thickener line. Based on the reactor design, the DO for the oxidation of pollutants is 3 mg/L. The first and second parts were filled with carriers at a rate of 40% of the reactor volume, and the surface area of the media was 500 m 2 /m 3 . Phosphates were then chemically treated by introducing 3 kg of calcium hydroxide (Ca(OH) 2 ) to the pilot system daily. This system used calcium hydroxide to efficiently remove phosphate and improve the nitrification process, resulting in a higher pH level, lower cost, and easier administration. It was injected into the pilot system's second stage, where the nitrification process occurs. The sedimentation basin was designed based on an HRT of 1.5 h, the number of plates 6, and an angle of inclination of 60°, and the distance between the plates is 5 cm. The system works electrically by a programmable logic controller (PLC).
Inoculum source for EN-MBBR system. Conventional activated sludge from nutrient biological removal units in the full-scale Karbala WWTP inoculated the EN-MBBR system. The seed used in the system's first stage was taken from the anaerobic phosphorus removal tank because it contains heterotrophic bacteria. The system's second phase was injected with seed from the internal recirculation from the oxic tank to the anoxic tank because it contains autotrophic bacteria. MLSS concentrations in heterotrophic and autotrophic seeds were 3500 mg/L and 3100 mg/L, respectively.
Pilot start-up and operational procedure. The pilot was operated on for 240 days, and the operating period is divided into two phases: the dynamic state (start-up) and the steady-state. Analyses carried out in the batch feeding phase were sampled from inside the pilot, and samples were taken from the outlet of the sedimentation basin in the case of continuous feeding. The start-up of the pilot was divided into two scenarios: the first scenario was batch feeding for 7 days, and the second scenario was continuous feeding for 21 days until reaching a steady state.
In the first scenario, the following procedure was applied 22 : (1) The pilot was fed by 3 m 3 /day of SSs on the first day, accounting for 10% of the total process flow. The air blower was turned on after SSs feeding until DO levels reached 4 mg/L. (2) An additional 10% of the SSs were fed into the EN-MBBR system on the second day, bringing the total SSs within 6 m 3 /day. The DO was kept at 4 mg/L, and media was introduced at a rate of 25% of the total design volume. (3) On the third day, 20% of the SSs were introduced, adding another 25% of media. (4) 20% of the SSs were fed into the pilot on the fourth day. The heterotrophic seed was added to the first phase of the pilot on this day, and the autotrophic seed was added to the second phase of the pilot on this day. The seed quantity per phase was 2 L. (5) On the fifth day, 20% of the SSs were introduced, along with a 25% increase in media. (6) On the sixth day, the SSs were fed to complete the total amount of media provided for the design, plus 2 L of seed per phase. DO concentrations were kept between 4 and 5 mg/L for the 6 days.
After that, the second scenario started after filling the pilot with the SSs until stability was reached. www.nature.com/scientificreports/ Karbala WWTP modeling in GPS-X. Hydromantis Environmental Software Solutions, Inc. 's GPS-X software version 8.1 (Educational license) was applied in this study. A popular standalone model includes integrated biological wastewater treatment procedures for ASP, anaerobic digestion system (ADS), and various physical and chemical interactions. In this study, a comprehensive ASP model was created by GPS-X software using the MANTIS2 model (carbon, nitrogen, and phosphorus processing), pH library, and simplified clarifier model. MANTIS2 model was readapted using ASM1 in the GPS-X software by integrating further changes relating to extra growth pathways in heterotrophic and autotrophic microorganisms. These pathways are (1) heterotrophic processes for anaerobic phosphorous removal, anoxic nitrate reduction by denitrification, and oxic organics degradation; (2) autotrophic processes for oxic ammonia removal by nitrification; (3) anaerobic ammonium oxidation (ANAMMOX) processes for ammonia and nitrite removal in anoxic condition; (4) endogenous respiration process of biomass oxidation 23 . More than 60 composite and state variables and various libraries of expressions characterizing the processes are included in the model, together with more than 30 stoichiometric and 24 kinetic input and output factors. This suggests that ASM1 was employed in the modeling process for carbon degradation, nitrification, denitrification, and phosphate removal 16 . With model limitations as per MANTIS2 24 , the assumptions for creating the model in this work are (1) At a content temperature, ASP runs; (2) There is sufficient mixing within the reactor and a consistent quantity of DO; (3) pH is stable and close to neutral; (4) there is simultaneous hydrolysis of organic and nitrogenous compounds; (5) There are sufficient amounts of inorganic nutrients to ensure adequate growth; (6) For any influent characteristics, the model's coefficients are assumed to be constants.
The following methods were used to model the data in this study using GPS-X: (1) the actual data of the Karbala sewage treatment plant required for GPS-X modeling was collected; (2) the effective fractionation of COD and nitrogen components has been modified and manually adjusted using the GPS-X influencer advisor to an acceptable condition and composite variables mass balance; (3) the model was run and calibrated by adjusting kinetics, stoichiometric, and other related matters typical parameters to get the best match between the predicted and the actual effluent quality data; (4) the calibrated model was validated using statistical analysis and a different set of wastewater quality data for the Karbala sewage treatment plant; (5) simulations were run under different scenarios to analyze the impact of related processes plant capacity and performance parameters in final effluent quality.

Calibration and validation of the simulated model by GPS-X. The model calibration simulation
was designed to estimate the best-fitted parameters for a specific collection of actual data obtained from Karbala WWTP. In this study, influent data for month 1 was used in model calibration, whereas the model validation input data was gathered from months 2 to 4. The Karbala WWTP's steady-state calibration and validation process steps are as follows: (1) Karbala WWTP biological treatment units and stages were simulated in GPS-X environment, as shown in Fig. 2; (2) the comprehensive carbon, nitrogen, phosphorus, and pH library from the GPS-X model (mantis2lib) was chosen; (3) to define the influent flow, COD, BOD, NO 2 , NO 3 , NH 4 , TSS, and PO 4 data from the plant were entered into the GPS-X influent advisor; (4) composite variables are passed to step 6 if the GPS-X influent suggests that the state's mass balance calculations requirements (i.e., organic, nitrogen, and MANTIS fractions) are met; (5) if there is an imbalance in the mass balance calculations in step 4 above, adjust manually the organic fractions, nitrogen, and MANTIS until a satisfactory state is reached and the composite variables are balanced may achieve; (6) calibration was performed using first-month data and GPS-X value's default kinematics and stoichiometric parameters. Calibration is completed when the model prediction fits all relevant effluent quality parameter data within acceptable limits; (7) if the GPS-X default model fails, the most sensitive parameters are initially screened and identified by optimizing the ASP and clarifier model parameters for 6 months of plant data by manually adjusting the relevant default parameters values one-by-one while visually observing the GPS-X output response collectively predicting the effluent quality parameters in terms of COD, BOD, NO 2 , NO 3 , NH 4 , TSS, and PO 4 ; (8) the first-month data calibration is restarted by changing values of the identified and screened parameters in step 7 to further improve the predictability; (9) the parameter values that contributed to the best-modeled prediction are determined in terms of standard effluent quality parameters, in this case, the final calibration is done; (10) the performance of the above-developed model is validated against a different set of actual Karbala WWTP data (i.e. months 2-4).
The schematic techniques for systematic model calibration and validation employed in this study are depicted in Fig. 3. To validate the simulated and measured data were quantitatively analyzed, and the goodness of fit was determined based on the performance measurement. To analyze the simulation results, the root mean square error (RMSE) and correlation coefficient (R) are utilized, as shown in the equations below: where Co is the actual data, C P is the modeled data, Co is the average of actual data, C P is the average of modeled data, and σ is the standard deviation over the dataset. These statistical criteria reasonable limits are 1 ≥ R > 0.8 and 0 ≤ RMSE < 1.5 11 . www.nature.com/scientificreports/

Modeling the impacts of novel process development by GPS-X. After reviewing previous studies
and research, no study has treated combined SSs (supernatant gravity thickener, underflow mechanical thickener, and centrate) away from the mainstream treatment system. Previous studies only referred to the treatment of the centrate line and the removal of ammonia from this line without concern for other pollutants. All studies are referred that after processing the centrate line, it is returned to the head of plant 25 . This study used a novel process to eliminate all the sidestream lines by treating them independently and mixing the treated sidestrams with the plant's outlet. Figure 4 depicts the proposed schematic diagram of the A 2 /O system for mainstream treatment and the novel process development of EN-MBBR for SSs treatment. By eliminating the SSs, the plant will be affected in terms of shocks, unit volumes of plant, bulk cost, and energy consumed. After successfully operating the pilot for more than 8 months without noticing any operational problems, a large-scale treatment plant similar to the pilot system will be designed according to data in Table 4. The proposed full-scale EN-MBBR has two parts: one to process organic substances and promote the growth of heterotrophic bacteria, and the other to treat nutrients because the active bacteria in this system are autotrophic. Calcium hydroxide is used to treat phosphorus in the EN-MBBR system chemically.

Results and discussion
Operating the pilot. Figure 5 depicts the start-up of organic material and nutrient fate until a steady-state is reached. Figure 5a shows that the mixed liquor suspended solids (MLSS) in the reactor increased in batch mode due to the conversion of the dissolved organic materials into suspended materials, which contributed to an increase in the suspended growth over the attach growth at the beginning of the process because the biofilm layer needs a period to the formation 26 . After completing the batch feeding period and passing MLSS to the sedimentation basin, a rapid decrease in suspended solids is observed due to the sedimentation process, decreasing from 610 to 11 mg/L. A slight decrease in COD and BOD concentrations was observed during the batch feeding phase, but during continuous feeding and after more than 10 days had passed, these concentrations decreased significantly and stabilized on day 17. The reason for the decrease in COD and BOD concentrations is due to the decomposition that occurred due to the presence and growth of heterotrophic bacteria present in the form of a suspension and attached, where all the appropriate conditions were provided to conduct this decomposition of the substrate through dissolved oxygen, temperature and nutrients 27 . Extended aeration was used in the second part of the EN-MBBR system for several reasons: the elimination of ammonia, the occurrence of endogenous respiration to reduce sludge 28 , and to remove a large part of the inert organic materials and no need for a primary sedimentation basin. The nitrification process in the EN-MBBR system eliminated the ammonium, as indicated in Fig. 5b. The nitrification process is carried out through two stages: the first is the oxidation of ammonia to nitrite by nitroso-bacteria, the second stage is the oxidation of nitrites to nitrates by the nitro-bacteria, both types of bacteria are autotrophic 29 . In the batch flow phase, a slight decrease in ammonia concentration was observed due to the consideration of the ammonia consumed in the first part of the pilot as a nutrient for heterotrophic bacteria to contribute to the removal of COD and BOD. A decrease in ammonia was observed after approximately the eighth day of operation due to the growth of nitrosobacteria to convert ammonia into nitrite. www.nature.com/scientificreports/ On the fourteenth day, ammonia stabilizes somewhat for two days due to the presence of a complex inversion between nitrite and nitrate, which somewhat contributed to the slowdown of the oxidation of ammonia. The slowdown in ammonia oxidation is due to the oxidation of nitrites and their conversion to nitrates. After the oxidation of nitrites by nitro-bacteria, ammonia begins to decline again. It was observed that the oxidation and stability of nitrite are faster than the oxidation and stability of ammonia. Nitrite concentrations increased significantly from the eighth day to reach their peak on the thirteenth day and then decreased rapidly to stabilize on the seventeenth day of starting operation. On the same day that nitrite concentrations decreased, nitrate concentrations reached 95 mg/L. The high concentrations of nitrates in the EN-MBBR system are due to several reasons, including the oxidation of high ammonia concentrations and their conversion to nitrates and proteins and amino acids present in its SSs, as well as due to endogenous respiration.
Furthermore, because hydrogen sulfide gas is impacted by numerous mechanisms, including stripping, biodegradation, and adsorption, hydrogen sulfide gas concentration declined faster than ammonium. After the stabilization of pilot biology, calcium hydroxide was injected into the nitrification zone, whereby adding this alkaline compound, phosphate concentrations decreased from 45 to 0.16 mg/L. This compound also increased  Karbala WWTP model calibration. Figure 6a observed that default predicted results (without calibration) were higher than predicted calibrated, especially for nutrients. The percentage difference between the default (uncalibrated) and the predicted (calibrated) COD, BOD, TSS, NO3 --N, NH4-N, and PO4-P were 30%, 60%, 60%, 28%, 97%, and 92%, respectively. This difference between default and predicted necessitated the modification of the most sensitive parameters affecting the substrate and nutrients together to produce results as close to the real as possible. This is done through steps 7-9 in "Calibration and validation of the simulated model by GPS-X" section. After adjusting the most sensitive parameters and achieving mass balance, the calibration results for the first month (Fig. 6a) were very close between the actual and the predicted.  www.nature.com/scientificreports/ The congruence was achieved by modifying composite variable stoichiometry, mixed tank volume, model stoichiometry parameters, kinetics parameters, and clarifier parameters, as shown in Tables 6 and 7. After checking that the simulated output was within the statistical confidence interval of the actual data values, the quality of the GPS-X calibration findings was determined (Table 6). If the actual and predicted results are within the acceptable framework stipulated in Eqs. (1) and (2), the difference is not significant. Furthermore, most of the calibrated model settings are substantially within the range of values published in the literature 16,30 . The lower or www.nature.com/scientificreports/ higher values of the parameters proposed by GPS-X, notably those related to nitrogenous substances (such as i XP, i XB , k A , and K NH4 ) in Table 7, contributed to the high GPS-X default calibration failure. These parameters were discovered to be extremely sensitive to small changes. Thus they had to be increased systematically by 3.2, 3.5, 3.75, and 3.8 times, respectively, to reach satisfactory calibration (  (1) and (2), model validation is defined as a high level of agreement between the model's predictions and a different collection of data that did not participate in the model's building. Months 2, 3, and 4 discrete average monthly effluent quality we're used to validating the model, and the Table 6. GPS-X input stoichiometry parameters (default and adjusted) based on GPS-X influent adviser.   Fig. 6b-d, respectively. Likewise, these validation results were in line with the current plant data and within the permitted range (Table 8). It's worth noting that all models' stoichiometry and kinetic parameters are similar and can be used for both calibration and validation (all month 1 to 4 simulated data). However, initial characterization of the influent stoichiometry fractions (Table 6) is required in each situation, dependent on influent wastewater quality factors, which vary from month to month (ranged values reported in Table 2). The model validation findings suggest that the calibrated models captured the biological processes at the Karbala WWTP for treating municipal wastewater well enough to be judged acceptable.

Effects of sidestreams elimination on mass loading shocks. The SSs contain high pollutants that
contribute significantly to the instability of the plant from time to time. Some studies have indicated that the SSs cause shocks to the plant while returning to the plant head. After adjusting the calibration and validation results of the plant, the plant was operated for 60 days according to two scenarios: the first (S1) was with the return of the SSs without treatment, and the second (S2) was not returned to the head of the plant, treated and mixed with the effluent. Figures 7 and 8 show the effect of SSs on pollutants' fate at dynamic and steady-state conditions. From Fig. 7, it is observed that the concentration BOD was the fastest pollutant to reach steady-state. This is due to the high biodegradability of BOD and the favorite growth of heterotrophic bacteria inside the system that contributed to the fast decrease of pollutant concentration 32 . Likewise, TSS quickly stabilized due to the formation of sufficient floc that contributes to an increase in the sludge mass and consequently the sedimentation increases, which reduces the concentrations of TSS at the outlet. The stability of TKN was observed before NO 3 because TKN contains organic nitrogen and ammonia, where organic nitrogen can be disposed of in conjunction with BOD due to the similarity of decomposing microorganisms 33 . Still, ammonia needs more than 7 days until the growth of the autotrophic bacteria responsible www.nature.com/scientificreports/ for the oxidation and transformation of ammonia to nitrites and nitrates 34 . It was observed that nitrites, nitrates, and phosphorus were stable after 12 days. When the MLSS concentrations reached 3500 mg/L, the WAS was started, and the necessary treatment for the sludge was carried out. After the necessary treatment of the sludge resulted in rejected water and supernatant from the sludge treatment units and processes, and this water after collecting it, which is called SSs, and returning it to the head of the plant caused a shock in the effluent of the plant lasted for more than 7 days. During this shock, high nutrient concentrations were observed, consistent with what previous studies reported on the SSs containing high pollutants of nutrients. Figure 8 shows the plant operation without returning the sidestream to the head of the plant, where a similar trend was observed in the pollutants kinetic. However, the shock was reduced to one day due to eliminating SSs and treating them in the EN-MBBR system. When the SSs were not returned to the head of the plant, this contributed to the absence of any shocks within the A 2 /O operations. The most critical issue addressed by the novelty of this study is the zeroing of shocks.
The feasibility of eliminating the sidestreams. Removing SSs from a full-scale A 2 /O process and their treatment separately by a full-scale novel EN-MBBR system resulted in economic and technical advantages.
The impact of removing sidestreams on the volume of the plant. After eliminating the quality and quantity of SSs and not returning it to the plant's head in Scenario 2, and based on the plant's new mass balance, each treatment unit was independently examined to calculate new volumes, as shown in Table 9. According to the results, the sizes of anaerobic and anoxic tanks were reduced by 4.61 and 4.04%, respectively. The size of these tanks is mainly determined by the raw wastewater entering the plant; thus, any decrease in inflow owing to the elimination of SSs will decrease tank volume. To clarify, the anaerobic process is responsible for treating orthophosphate  www.nature.com/scientificreports/ (PO4 -3 ), and the anaerobic tank's HRT should not exceed 3 h; otherwise, secondary phosphorus release may occur. As a result, reducing sewage discharge (S2) may increase the HRT more than necessary, resulting in secondary phosphorus release 35 . As a result, reducing tank volume is necessary to keep the process under balance. Meanwhile, nitrate reduction (NO 3 − ) is carried out using an anoxic procedure with an HRT of 1 to 3 h. On the other hand, low flow rates will increase the HRT beyond the limitations; thus, volume reduction maintains an appropriate HRT for NO 3 − reduction. The quantity and quality of effluent determine the aeration tank's design. A suitable reactor volume, aeration, and mixing are required to treat organic and nutrient loads. As a result of removing SSs and their loads, the volume of the aeration reactor can be reduced by up to 10%. The surface area of the secondary settling tank is related to the hydraulic load, and since more than 4000 m 3 /day of flowrate has been eliminated, the surface area has been reduced by more than 2.8%. The EN-MBBR system contributes to the decrease of sludge produced through endogenous respiration for sludge within the system, which is one of its advantages. As a result, the sludge load is lowered significantly, and because the drying bed design is based on the resulting sludge load, the drying bed area was reduced by 10%.
The impact of removing sidestreams on energy consumption. Because of the high concentrations of pollutants in the SSs, especially ammonia, and to treat these pollutants, sufficient quantities of oxygen and adequate mixing are required. The availability of these two factors consumes enormous amounts of energy, especially in the suspended growth of bacteria 36 . As indicated in Table 10, the EN-MBBR system has effectively contributed to the treatment of SSs and the economic feasibility of the WWTP by lowering energy consumption. During reactor design, all equipment is primarily dependent on the flow and volume of the reactor. After eliminating the SSs, this reduced the energy in the preliminary treatment phase by more than 16 kW/h due to the removal of more than 4000 m 3 /day. It was also noticed that during the primary treatment phase, the power decreased by more than 12 kW/h due to the removal of the SSs flow and not returning it to the plant head. The energy in the anaerobic basin decreased by 0.6 kW/h, while in the anoxic basin was reduced by 2.15 kW/h. Elimination SSs contributed to reducing more than 30% of the ammonia load in the plant; thus, this was reflected in the pumps for internal recirculation (IR) from the nitrification to denitrification process. Therefore, the power of the IR pumps was reduced to more than 20 kW/h. The highest decrease in the consumed energy was observed in the aeration basin due to not returning the substrate and nutrients to this basin, which contributed to the absence of the need to mix and large quantities of oxygen to oxidize the SSs pollutants 37 . This resulted in a reduction of energy by 145 kW/h. The energy consumed in the sedimentation basin with the RAS and WAS regulating process decreased by 8 kW/h. The energy consumed in the chlorination basin was not affected because the SSs was returned after treatment to this basin. It was observed that the power increased in the mechanical tanks by 2.85 kW/h due to the return of the sludge generated from the EN-MBBR system to this unit. By operating the EN-MBBR system, it needs 35 kW/h of power, but this system, and through the sludge generated from it and returning this sludge to the anaerobic digester, generated power of more than 270 kW/h. Considering the energy-saving and consumption in all treatment units, the total energy saved is 442.51 kW/h. The novel process reduced the energy consumed by 17% of the total plant power. This indicates that eliminating the SSs positively contributed to the energy-saving at the full-scale www.nature.com/scientificreports/ Karbala WWTP. Before SSs elimination, most of the energy was consumed to oxidize ammonia, organic matter, and hydrogen sulfide in the aeration tank, not to mention the technical and environmental issues associated with operating aeration tanks. The diversion of SSs treatment by the novel EN-MBBR system was supported by attached growth mechanisms that contributed to this energy reduction.
The impact of removing sidestreams on cost. The cost depends mainly on sewage treatment plants on the size of the units, energy consumed, maintenance, and land area. Table 11 shows the cost comparison of different units in S1 and S2. Construction, operating, and maintenance expenses are logically lowered after reducing sizes and energy consumption in S2. SSs are processed independently utilizing the EN-MBBR system, which uses attached growth processes and media with vast surface areas to minimize costs. A good reduction in the total cost of the plant was observed by more than 2 million USD $, due to the decrease in the proportion of the volumes of the units of the plant and the energy consumed. Most previous studies did not reach high-cost savings because they focused only on the separate centrate line treatment and ammonia removal without addressing other pollutants. Elimination of SSs, which amounts to more than 4% of the plant's total discharge, positively reflected all the factors that affected the cost. Eliminating SSs means removing large amounts of ammonia, organic matter, and hydrogen sulphide that need large amounts of oxygen to oxidize. This oxidation needs energy to operate the equipment. Eliminating the quantity and quality of the SSs and not returning it to the head of the plant contributed to saving sums of money in the primary lifting station and biological treatment by 13,500,066 and 9,693,000 $, respectively. The reduction in nitrogen concentrations resulting from eliminating the SSs contributes to reducing the demand for energy and alkalinity, which has been positively reflected on the internal circulation pumps, and sums of money are available estimated at 339,999 $. Elimination of the SSs also contributed to a significant reduction in the quality and quantity of water entering the secondary sedimentation basins due to reducing the quantitative and hydraulic loads of the sedimentation basins. The amounts available because of this elimination were 12,175,925 $. In addition to pollutants treatment and energy saving, the EN-MBBR system also reduces sludge production, which is positively reflected in sludge management and the spaces designated for drying beds.

Comparison of EN-MBBR system with other reactors.
After evaluating all past research, particularly those with patents, it becomes clear that they have all dealt with and concentrated on reject water (centre), implying that SSs have not been permanently eradicated. All SSs lines were deleted in this research, resulting in several benefits compared to other studies, as shown in Table 12 [38][39][40][41][42] .

Conclusion
To our knowledge, this is the first study of its sort that looks into the prospect of reducing the effects of SSs on sewage treatment plants by treating them separately and without returning them to the head of the plant. The majority of nutrients and organic contaminants are found in SSs. As a result, a novel prototype EN-MBBR system was successfully designed, manufactured, and operated to treat organic pollutants in the first half of the system, and ammonia, phosphorous, and H2S in the second portion of the system. The implications of eliminating SSs on a full-scale A 2 /O system were also effectively simulated using the GPS-X7 simulation. Eliminating SSs gives technical and economic feasibility, according to the modeling. These include the removal of loading shock and high cost, energy, and volume savings. Furthermore, measured data from full-scale A 2 /O and pilot EN-MBBR systems were used to validate the calibrated model.