Evaluation and optimization of ASM1 parameters using large-scale WWTP monitoring data from a subtropical climate region in Brazil

This study aimed at providing a set of optimal kinetic and stoichiometric parameters of ASM1 representative of wastewater from a subtropical climate region in Brazil. ASM1 was applied on STOAT program, and the model parameters were evaluated and optimized with sensitivity analysis and Response Surface Methodology (RSM) to reach minimum prediction errors of effluent TSS, COD, and NH3. Six sensitive parameters were identified: YH, YA, μA, KNH, bA, and kOA. Predictions of RSM regression models were strongly correlated to the STOAT predictions. YH mainly affected TSS and COD, and the other parameters affected NH3. ASM1 calibration with estimated optimal values of sensitive parameters resulted in approximately null prediction errors for modeling state variables. NH3 presented similar results in the ASM1 validation; meanwhile, TSS and COD presented high errors related to the increase in YH due to the RSM optimization. The optimal parameters, mainly YA, μA, KNH, bA, and kOA, constitute references for other studies on ASM1 modeling using wastewater data from a subtropical climate region. YH optimal value should be evaluated as well as the effect of sludge wastage methods and the simulation periods.

to model an activated sludge (AS) process. STOAT and ASM1 may be used with scarce data sets available from standard monitoring of large-scale Wastewater Treatment Plants (WWTP) for setting up a facility and modeling the performance of AS under steady-state conditions (Andraka et al. 2018).
The calibration is one of the most critical steps in the modeling procedure to provide reliable predictions accordingly with the specific conditions of AS process (Rieger et al. 2013). ASM1 modeling and calibration depend on wastewater composition fractionation and model parameters determination, respectively. Such fractionation and determination may need extra laboratory analysis, which is not commonly viable for some plants (Borzooei et al. 2019). In some cases, typical ratios of municipal wastewater may be used for fractionation (Henze & Comeau 2008, p. 36), and kinetic and stoichiometric ASM1 parameters calibration on trial and error procedure may be adopted. However, this calibration procedure is not appropriate due to the importance of prioritizing methods that guarantee as much information as possible to form a suitable parameter combination (Petersen et al. 2003).
Alternatively, some authors have proposed a systematic calibration approach to find optimal parameters for ASMs. The strategy comprises modeling in association with sensitivity analysis for parameter selection, and design of experiments and Response Surface Methodology (RSM) for parameter evaluation and optimization (Kim et al. 2009;Lim et al. 2012;Ahn et al. 2014).
Most studies on the evaluation and optimization of ASM1 parameters were developed in Europe, North America, and Asia (Hauduc et al. 2011). Furthermore, ASM1 presents default values of parameters for European wastewater characteristics (Ahn et al. 2014). Thus, applications in other regions are relevant due to the variation of model parameters correspondingly to wastewater composition along space and time (Von Sperling et al. 2020).
The geographic expansion may also contribute to evaluating ASM1 modeling using STOAT in other climate regions. That program is broadly used in the temperate climate of the United Kingdom and was validated in tropical climate conditions in India, where the microorganisms' growth rate is higher than the temperate climate conditions (Sarkar et al. 2010).
Particularly in Brazil, AS modeling is essential to improve wastewater treatment. The country presents 354 WWTPs with such treatment process (10% of the total) (ANA 2020), and more than half of the population lack wastewater treatment (Brasil 2019). Some authors have used ASM1 and STOAT to model a WWTP using monitoring data of domestic wastewater treatment from a subtropical climate region in Brazil (Pistorello 2018;Baptista 2020). However, to the best of our knowledge, no studies have performed an estimation of optimal ASM1 parameters associating the use of monitoring data of a large-scale WWTP within a subtropical climate region with sensitivity analysis and RSM.
Therefore, this study aimed at providing a set of optimal kinetic and stoichiometric parameters of ASM1 representative of wastewater from a subtropical climate region in Brazil. Thus, these optimal values may compose the references of ASM1 parameters and be effective on AS modeling using the STOAT program and monitoring data of large-scale WWTP in the same climate conditions.

WWTP of study
The research was conducted on domestic wastewater generated from São João Navegantes WWTP in Porto Alegre, Brazil, located in a subtropical climate region (29°59 0 29″S; 51°11 0 43.5″W). The average raw-water flow is 0.44 m 3 /s, corresponding to 150 thousand inhabitants. The wastewater system comprises extended aeration activated sludge (Figure 1).

Data collection and reconciliation
It was determined the wastewater characteristics: flow, temperature, total suspended solids (TSS), chemical oxygen demand (COD), ammonia (NH 3 ), and pH. Influent and effluent means of these variables were considered after treating missing and censored data, as well as outliers of the WWTP monitoring data set, as described in von Sperling et al. (2020). Annual means of those variables for 2018 and 2019 were used for ASM1 calibration and validation, respectively (Table 1).

Activated sludge modeling
The ASM1 and Version 3 models were used, respectively, for the aeration and secondary sedimentation tanks. Both models are available on STOAT © simulation program and it was considered the guidelines of Rieger et al. (2013) and WRC PLC (1994) for modeling conduction. Modeling was developed for one of the four parallels flows of the treatment plant. Information of tank sizes and operational data were provided by the WWTP (Table 2).
On Version 3, the wastage method selected was fixed-rate (0.00153 m 3 /s) over variable time (0-8 h) to maintain a specific Mixed-Liquor Suspended Solids (MLSS) set-point (3,000 mg/l). By assuming extended aeration correlation, sewage calibration data for the same model considered Stirred Sludge Volume Index (SSVI) (3.5 g/l) equal to 100 ml/g (average for fair sludge sedimentation) (von Sperling 1994).
AS Solids Retention Time (SRT) was estimated using Equation (1) (Tchobanoglous et al. 2014), according to the WWTP monitoring data (Table 2). V is the volume of the aeration tank, X is the MLSS concentration, Q W is the sludge wastage flow (0.00153 m 3 /s considering a pump run time of 8 h), and X R is the sludge TSS. The first simulation period was assumed as 3 SRT (60 d) (Rieger et al. 2013).
The influent profile assumed on STOAT was the sinusoidal pattern. Thus, the following conditions were selected for most wastewater variables (Table 1): 0 h for phase, 30% for amplitude, and frequency equal to 0.261799 (representing daily fluctuations). Specifically for temperature, the frequency was 7.27E-4 (annual fluctuations), and amplitude was 0% for the same variable and pH. Typical ratios of municipal wastewater were applied at the fractionation of TSS, COD, and N contents ( Table 3).
The modeling state variables were the mean effluent concentrations of TSS, COD, and NH 3 ( Table 1). The following runs started from the end of the first simulation.
ASM1 parameters and respective default values are presented in Table 4. Some kinetic parameters display different units on STOAT, i.e., h À1 for 15°C (P 15°C ). By using the temperature coefficients (θ) provided on

Sensitivity analysis
A one-way sensitivity analysis was performed to indicate the critical parameters of ASM1 (Table 4). This analysis consisted of sequentially varying each parameter while keeping the others constant. A 10% increase was applied for each parameter to calculate its sensitivity coefficient (S i j ) (Equation (3)) related to each state variable; where Y  is the state variable output, P is the ASM1 parameter value, and the subscripts 0 e 1 stand for the default and the changed values, respectively. A parameter is considered sensitive if S i j ! 0:25 (Liwarska-Bizukojc et al. 2011;Andraka et al. 2018;Chen et al. 2020). Herein, the period of each simulation was equivalent to the estimated SRT (Table 2). (3)

Parameters optimization
Sensitive kinetic and stoichiometric parameters of ASM1 were optimized aiming at minimum prediction errors for state variables on simulations, according to the WWTP monitoring data (Lim et al. 2012). Response Surface Methodology (RSM) and Central Composite Design (CCD) were conducted using Minitab Statistical Software (Version 20.2), whose input consisted of the difference between the observed data recorded by the WWTP (Table 1) and the result of simulation.
Oxygen hsc for autotrophs Correction factor for anoxic growth of heterotrophs Hsc for hydrolysis of slowly biodegradable substrate Regression models were created for each modeling state variable, and ANOVA and standardized effects were used to evaluate the influence of changes in sensitive parameters on predictions errors of effluent TSS, COD, and NH 3 . Regression models' adequacy was verified using the Lack-of-fit test. Also, it was analyzed the correlation between results of simulations on STOAT and predictions of the regression models, as well as the RMSE (Equation (4)) of each regression model (Kim Rao & Yoo 2009;Lim et al. 2012;Ahn et al. 2014).
X is the simulation prediction using STOAT, Y is the regression prediction via RSM, and n is the number of simulations designed by CCD.
In the multiple optimization, the desirability function based on a target value was used to estimate optimal values for sensitive parameters, which would lead to minimum prediction errors of effluent TSS, COD, and NH 3 (Kim Rao & Yoo 2009;Lim et al. 2012). According to the modeling objective, it was assumed a range of minimum prediction error for each state variable: +5 mg/l for TSS and COD, and +1 mg/l for NH 3 (Rieger et al. 2013).
Finally, optimized parameters of ASM1 were applied on simulations on STOAT for model calibration and validation. It was verified whether or not the results complied with the previous range of prediction errors determined for each state variable. The period of each run for parameter optimization was also the estimated SRT (Table 2).

Sensitive parameters of ASM1
The sensitivity analysis required 20 runs on STOAT which is considered a fast way to evaluate sensitive model parameters. The present study identified six sensitive parameters of ASM1 on modeling using domestic wastewater data from the São João Navegantes WWTP (Figure 2). COD presented no sensitive parameters; meanwhile, NH 3 was the most affected state variable due to changes in the ASM1 kinetic and stoichiometric parameters. Y H was sensitive for the TSS and NH 3 outputs. Although the same parameter was not sensitive for COD, it presented a coefficient of 0.22, indicating an effect on that state variable. Lim et al. (2012) observed similar results and commented that the relationship between Y H and the state variables TSS and COD are related to sludge production. The heterotrophic yield also depends upon the nature of the substrate and the population of microorganisms performing degradation (Henze et al. 1987).
Other sensitive parameters only affected NH 3 predictions. In particular, the sensitivity coefficients of μ A , K NH , and b A presented high values (!1), characterizing a substantial effect on the NH 3 outputs of the WWTP AS process modeling. Contrarily, μ A also affected effluent TSS in another study (Lim et al. 2012). Nevertheless, the high sensitivity of μ A may be related to elevated ammonia loads and to a lack of complete nitrification (Levy 2007). This parameter is the most critical one for characterizing the growth of the autotrophic biomass and is influenced by many environmental factors, such as pH and temperature (Henze et al. 1987). On the other hand, the value of b A is considered difficult to measure and may be assumed (Levy 2007), as well as the values of Y A and K OA that also affected effluent NH 3 (Figure 2(c)). K NH and Y H are recommended to be evaluated for each type of wastewater (Henze et al. 1987).
It is noteworthy that 13 parameters of ASM1 were not sensitive for the predictions of effluent TSS, COD, and NH 3 of the São João Navegantes WWTP, regarding the specific modeling conditions. It meant that fewer parameters need to be evaluated, which sped up the ASM1 calibration. However, even though the parameters I XB , μ H , k S , and b H were not considered sensitive, they deserve attention as long as they affected COD and NH 3 predictions (Figure 2(b) and 2(c)).
Half of the six sensitive parameters identified (Y H , μ A , and b A ) were also subject to change in another study that used ASM1 modeling results of 18 WWTPs in Europe, three in Asia, and one in North America (Hauduc et al. 2011). Other authors confirmed the relevance of Y H , μ A , K NH , and K OA for either steady-state or dynamic AS simulation conditions (Petersen et al. 2003). Furthermore, sensitive parameters determined for steady-state conditions can significantly support the dynamic calibration (Liwarska-Bizukojc et al. 2011). Therefore, the results of the present study indicate a set of sensitive parameters of ASM1, which besides being representative of a subtropical climate region, may be used in dynamic modeling practices.

Central composite design
CCD resulted in 90 cases randomly designed to find optimal values for sensitive parameters that minimize the prediction errors of effluent TSS, COD, and NH 3 . The levels of the test were coded and the ranges are shown in Table 5.
The experimental ranges of Y A , K NH , and K OA were close to the default values of these parameters (Table 4). It corroborates other findings, as changes in the Y A aren't usual, and half-saturation coefficients (K NH and K OA ) depend on environmental conditions and don't require significant changes (Hauduc et al. 2011). However, Y H , μ A , and b A were tested in different ranges (Table 5), which may be related to the variability of such parameters (Hauduc et al. 2011).
ASM1 is a deterministic model and repeated runs result in the same outputs. Since RSM is a statistical model, it requires variability. Thus, random changes were performed on the other parameters which affected the state variables but were not considered sensitive by the 0.25 limit (Figure 2) to provide such variability. These changes were made in 13 out of 14 repeated cases designed by CCD, and the parameters changes were i XB , b H , μ H , and K S (Table 6).

Regression models
The prediction errors of effluent TSS, COD, and NH 3 were calculated from the results of runs designed by CCD, which means the difference between the observed data recorded by the WWTP and the STOAT simulation's prediction. Such prediction errors were input onto RSM to generate the regression models for each state variable. The regressions coefficients estimated by RSM are presented in Table 7. The respective non-   Uncorrected Proof significant quadratic relationships and interactions between parameters weren't considered to simplify the regression models. ANOVA showed similar results for TSS and COD. Changes in most sensitive parameters affected TSS and COD predictions (p , 0.001). Only K NH e k OA didn't affect those state variables (TSS: F (1,89) ¼ 0.22; p ¼ 0.639 and F (1,89) ¼ 0.28; p ¼ 0.598, respectively; COD: F (1,89) ¼ 0.36; p ¼ 0.552 and F (1,89) ¼ 0.45; p ¼ 0.505, respectively). All quadratic relationships also influenced TSS and COD, especially Y H *Y H for both variables (F (1,89) ¼ 197.60 and 109.12, respectively; p , 0.001). However, interactions between sensitive parameters weren't statistically significant for TSS and COD and were disregarded in their regression models (Table 7).
The ANOVA of the NH 3 regression model showed that Y H was the only sensitive parameter that did not affect the predictions of such state variable (F (1,89) ¼ 1.69; p¼ 0.198). Although only quadratic relationships of μ A and b A were statistically significant in that regression model (F (1,89) ¼ 303.27 and 35.70, respectively; p , 0.001), four interactions between parameters affected effluent NH 3 prediction errors: Y A *μ A ; μ A *b A ; μ A *K OA and b A *K OA (Table 7). Response surfaces of these significant interactions are shown in Figure 3.
The relevant impact of μ A *b A interaction (F (1,89) ¼ 189.17; p , 0.001) indicates the many combinations that may be assumed for such parameters to estimate minimum prediction errors for effluent NH 3 (Figure 3(a)). Hauduc et al. (2011) highlighted the variability of values for those parameters. Such relationship is usually inversely proportional, as reduction of μ A and b A increases and decreases, respectively, ammonia effluent concentration. Also, an increase in b A may inhibit cell growth and replace nutrients in the system, elevating NH 3 concentration. However, the effect of b A is not the most important for effluent NH 3 control (Levy 2007).
In the interactions of Y A *μ A (F (1,89) ¼ 11.43; p ¼ 0.001) and μ A *K OA (F (1,89) ¼ 10.08; p ¼ 0.002), the prediction errors of effluent NH 3 generated by the regression models were null for higher tested values of μ A (Figure 3(b) and 3(c)). Changes in Y A and K OA on such interactions didn't show a distinct effect on predictions errors, as observed by other authors (Petersen et al. 2002).
Finally, the interaction between b A *K OA (F (1,89) ¼7.83; p¼0.006) was the less significant in the NH 3 regression model. In that relationship, the whole experimental range tested for K OA was suitable to obtain minimum effluent NH 3 prediction errors (Figure 3(d)), though only lower values adopted for b A appropriately worked.

Uncorrected Proof
The standardized effects of the regression models pointed out the influence of Y H on TSS (87%) and COD (109%) prediction errors (Figure 4(a) and 4(b)). Linear and quadratic relationships of other parameters presented effects lower than 15%, even though they were significant in the ANOVA results.
For NH 3 predictions errors, μ A and b A stood out and indicated contrary effects on responses: À40 and 31%, respectively (Figure 4(c)). This situation is also shown in Figure 3(a).
The lack-of-fit tests showed that p values were not statistically significant for regression models of TSS and COD, which means that the alternative hypothesis that indicates the lack of fit of those regression models wasn't selected. On the other hand, the NH 3 regression model presented opposite results, revealing that such model didn't show the goodness of fit (Table 8).
For that specific case, to correct the lack of fit of the regression model, the original variable predicted on simulation using STOAT was used (prediction of NH 3 effluent concentration instead of the prediction error). Then, it  Uncorrected Proof was expected to obtain only positive values, which could be transformed to improve the fit of the regression model. However, it didn't show different results, and the first regression model was maintained (Table 7). The estimated prediction errors of all state variables provided by RSM were compared with those predicted using STOAT. Such relationships presented a strong correlation (R² .0.97), as shown in Figure 5. Thus, the regression models explain at least 97% of the variations on prediction errors of effluent TSS, COD, and NH 3 . The statistical error via RMSE presented lower values, i.e., close to 0, indicating that there wasn't a discrepancy between prediction errors estimated by RSM and predicted via STOAT. Those results are similar to the findings of other authors (Kim et al. 2009;Ahn et al. 2014).

Parameter optimization
The generated regression models of RSM were used to estimate optimal values for the sensitive parameters to be applied on the calibration and validation of ASM1. The target value onto the desirability function in the optimization procedure was zero to estimate parameters to reach minimum prediction errors of effluent TSS, COD, and NH 3 . RSM optimization generated the output shown in Figure 6. The optimal values for the sensitive parameters are displayed in red within a specific range. The optimal prediction errors for each state variable are in blue. The vertical lines in red show the optimal parameters within the experimental range of CCD (Table 5). Figure 6 also presents the relationships between each parameter and the modeling state variables. Changes in Y H mainly affected prediction errors of effluent TSS and COD. Variations in the values of the other sensitive parameters affected NH 3 prediction errors. It's important to emphasize the effect of μ A and b A , whose interaction was the most significant in the ANOVA results of the NH 3 regression model. Nonetheless, μ A showed two optimal points: on the central and upper regions of the parameter experimental range. The first one was more prominent and consisted of the optimal value determined for the parameter. That situation may be related to the great importance of μ A for the NH 3 predictions (Petersen et al. 2002), as shown in Figure 3(a)-3(c).
The single value of D, which lies between zero and one, is the overall desirability of the combined response levels. D close to one means that the balance of the properties becomes more feasible. On the other hand, D close to zero indicates that one of the response variables is unacceptable (Lim et al. 2012). The composite desirability showed D¼0.9181, pointing out a balance among the responses regarding the optimal values estimated for Uncorrected Proof the sensitive parameters (Kim et al. 2009;Ahn et al. 2014). Each state variable also presented an acceptable value of d (.0.87). Those ds may also be related to the lower and upper limits of prediction error selected for the state variables in the optimization procedure (+5 mg/l for TSS and COD, and +1 mg/l for NH 3 ). Then, the higher the limit (i.e., for TSS and COD) the higher was the estimated prediction error, even though all of the three ones were within +1 mg/l (in blue in Figure 6).
Moreover, the NH 3 desirability function (d¼0.99331) was the highest among all state variables, besides presenting the lower RMSE (Table 8). Despite the lack of fit of the NH 3 regression model, d and RMSE values are strong reasons to use this regression model in the parameter optimization (Kim et al. 2009;Ahn et al. 2014).
The optimal value estimated for Y H was 0.9903, placed between the minimum and central levels of the experimental design (Table 5). That value is relatively higher than the default one (0.67) (Henze et al. 1987). Also, other authors optimized the same parameter and had a maximum value of 0.67 (Hauduc et al. 2011). Baptista (2020) used Y H 0.45 on AS modeling of the same WWTP of study. However, that author considered a period of monitoring data between 2006 and 2010 which presented lower values of effluent COD (43 mg/l) and TSS (15 mg/l), compared to the means obtained in the present study for 2018 (COD: 52 mg/l; TSS: 30 mg/l) (Table 1). Such difference may explain the lower value of heterotrophic yield used on ASM1 calibration by that author.
Furthermore, the optimal value calculated in the present study for Y H is close to 0.9380, defined in the optimization performed by Lim et al. (2012). The authors also worked on AS modeling using ASM1, even though the physical conditions (flow, sludge wastage, and RAS) were lower than the records of the São João Navegantes WWTP. This could have accounted for the difference in the Y H value. In addition, the experimental range tested by those authors was also lower, and the optimal value estimated in that study comprised the higher level of test, i.e., close to the one found in the present study. Those results highlight the variability of Y H and its relationship with the influent wastewater composition (Henze et al. 1987;Hauduc et al. 2011).
The other sensitive parameters of ASM1 evaluated via RSM are related to NH 3 predictions, mainly μ A and b A , whose changes had a considerable effect on that state variable ( Figure 6). The optimal values estimated for those parameters were 0.012 h À1 (15°C) and 0.0054 h À1 (15°C), respectively. These values were converted to the default unit of ASM1 using Equation (2) and resulted in 0.5 d À1 (20°C) and 0.22 d À1 (20°C). The optimal Uncorrected Proof value estimated for μ A remained in the usual range presented in the literature (0.2-1) (Ahn et al. 2014). For b A , the optimal estimative was higher than the maximum value observed (0.05-0.15) (Henze et al. 1987).
The optimal values of Y A 0.2575 and K OA 0.4536 ( Figure 6) were similar to the default values of such parameters: 0.24 and 0.4, respectively (Ahn et al. 2014). On the other hand, K NH 1.4 also showed a higher value than the usual range of 0.75-1 (Hauduc et al. 2011). Table 9 shows the comparison of observed data, default results, and verification of parameter optimization results. Prediction error ranges and prediction errors estimated by RSM and via simulations on STOAT are also presented.
The predictions of the simulation using default parameters resulted in high prediction errors for all state variables (.45%) compared to the observed data of 2018. Effluent NH 3 prediction of that simulation presented the highest variation (90%). Such difference was also found by Baptista (2020); meanwhile, for TSS and COD, the author obtained higher prediction errors compared to the observed data of that study. Therefore, the simulation results using default parameters reinforce the importance of ASM1 calibration for every type of wastewater. Most parameters of ASM1, especially kinetic ones, vary accordingly to the wastewater characteristics throughout space and time (Von Sperling et al. 2020).
The ASM1 calibration used the optimal parameters resulting in approximately null prediction errors. These errors were different from the ones estimated via RSM, though they were within the range selected for each state variable. Thus, ASM1 was calibrated by applying the optimal parameters estimated via RSM.
The ASM1 validation was performed to verify the applicability of the estimated optimal parameters in simulation using data of a different period, i.e., 2019 (Table 1). The results of validation presented around null prediction error for effluent NH 3 (Table 9). However, high differences were obtained for TSS and COD.
Such differences are related to the increase in Y H compared to the parameter default value (Table 4), i.e., 0.67 to 0.99. Even though the regression models of TSS and COD presented goodness of fit, the values tested for Y H (Table 5) may not have considered the whole experimental domain of the parameter. Thus, a new experimental range may be evaluated in association with other ranges determined for the other sensitive parameters identified in the present study.
Another explanation for the discrepancy observed consists in the differences in effluent concentrations of TSS and COD between 2018 and 2019, used on ASM1 calibration and validation, respectively. The calibration data presented higher values of TSS and COD compared to 2019 (9 and 7 mg/l higher, respectively), requiring an increase in the heterotrophic yield to reduce the prediction errors for both state variables. The effluent concentration of NH 3 was also higher in 2018, though it presented a slight difference related to 2019 (0.7 mg/l higher). The increase in Y H raises the COD transformation efficiency, determining the biomass concentration in the reactor. Hence, it has a high effect on the microbial growth rates and on the sludge production (Petersen et al. 2002;Levy 2007;Chen et al. 2020). Therefore, adjustments in the sludge wastage method may lead to obtaining a suitable Y H value for ASM1 calibration and validation, aiming at minimum prediction errors of TSS and COD for the specific modeling conditions of the present study.
After that, the sludge wastage method adopted in the present study, i.e., fixed-rate over variable time, may also be the reason for the TSS and COD oscillation peaks observed in the ASM1 calibration and validation simulations ( Figure 7).
As estimated values of the AS operation informed by the WWTP (Table 2) were used to determine the sludge wastage method, it may have affected the process performance. Furthermore, the SRT estimative was also based on WWTP data and is linked to the period determined for simulations (20 d¼480 h). Baptista (2020) modeled the same WWTP in steady-state conditions and applied 10 days as the period of simulation. In the ASM1 calibration Figure 7 | Effluent concentration (mg/l) of the modeling state variables predicted on STOAT using estimated optimal parameters. Measured means are from the WWTP's monitoring dataset (Table 1).
Water Practice & Technology Vol 00 No 0, 14 Uncorrected Proof of the present study, the oscillation peaks were at the end of the period hence after 10 days of simulation (Figure 7(a)). Although in the ASM1 validation the peaks started after 100 hours (Figure 7(b)), i.e., less than 10 days of simulation, the SRT estimative may have been overrated. Therefore, accurately AS operation data would aid in improving process performance on modeling. Also, other methods of sludge wastage (constant or fixed time over variable rate), and changes in operational parameters of ASM1 and Version 3 may contribute to control predictions of effluent TSS and COD.
Moreover, such peaks always happened when the Y H upper level of test was chosen in the simulations performed via CCD (Table 5) (see Results of CCD simulations as Supplementary Material). Then, those peaks are also related to the increase of Y H and its influence on the biomass concentration in the treatment process (Levy 2007).
The heterotrophic yield sensitivity is highlighted in optimizing ASM1 parameters (Hauduc et al. 2011). It's important to determine the value of such parameter for every type of wastewater, as Y H is commonly changed in ASM1 modeling (Henze et al. 1987) and each WWTP should be modeled assuming specific conditions related to its reality and challenges (Borzooei et al. 2019). Then, the wastewater physical and chemical composition, especially temperature, may have influenced Y H hence the predictions of effluent TSS and COD on STOAT using the optimal value estimated for the parameter. Also, the results of the present study support that attention should be dedicated to the calibration of Y H .
Otherwise, in addition to the high variation of NH 3 prediction in the simulation using default parameters (Table 9), all six sensitive parameters affected this state variable (Figure 2(c)). Five of these parameters are related to the autotrophic microorganisms (Y A ), N content, and specifically to ammonia (μ A , K NH , b A , and k OA ). The optimal values estimated for all of them were appropriate for both calibration and validation of ASM1, and the NH 3 prediction errors were lower than 1 mg/l. Therefore, such parameters are related to wastewater concerning subtropical climate conditions being a reference for other applications.

CONCLUSIONS
AS process of a large-scale WWTP within a subtropical climate region in Brazil was modeled using standard monitoring data set, the ASM1 model, and the STOAT program.
With sensitive analysis, six parameters of ASM1 were identified as sensitive to effluent TSS, COD, and NH 3 predictions: two stoichiometric (Y H and Y A ) and four kinetic (μ A , K NH , b A , and k OA ). Sensitive parameters were evaluated and optimized via RSM aiming at minimum prediction errors for modeling state variables. Y H mainly affected TSS and COD predictions; meanwhile, NH 3 predictions were influenced by changes in the other sensitive parameters.
Estimated optimal parameters were used on ASM1 calibration resulting in approximately null prediction errors for effluent TSS, COD, and NH 3 . In the ASM1 validation, the same was obtained for NH 3 , even though high prediction errors were observed for TSS and COD. Such errors are related to an increase in Y H through parameter optimization that highlights the sensitivity of this parameter on ASM1 modeling.
The parameters optimized in the present study, mainly Y A , μ A , K NH , b A , and k OA , whose optimal values were appropriate for effluent NH 3 predictions in the ASM1 calibration and validation, constitute references for other studies on modeling using wastewater data from a subtropical climate region. In particular, those optimal parameters may be effective on either steady-state or dynamic AS modeling using the STOAT program and monitoring data of large-scale WWTP in the mentioned climate conditions.
The quantity and quality of the WWTP monitoring data set used to characterize the AS process and wastewater, and the typical ratios adopted at influent fractionation restrict the application of the optimized parameters. Also, factors related to the sinusoidal simulation pattern and the sludge wastage method consist of specific modeling conditions and deserve attention.
The optimal value estimated for Y H needs to be evaluated considering tests on a new experimental range and/ or adjustments in the sludge wastage method.
The optimal parameters may be tested on a ASM1 cross-calibration/validation, i.e., data of validation (2019) for calibration and calibration (2018) for validation; on dynamic modeling studies; and applied at the modeling of similar WWTPs. It is recommended to use the optimal parameters for evaluating the simulation period and different sludge wastage methods on ASM1 and STOAT modeling. Furthermore, it's important to assess and optimize