P. aeruginosa interactions with other microbes in biofilms during co-infection

1.   Introduction

The role of water quality in streams, lakes, and seas can be stated as organic, synthetic, and environmental condition of waterbody ^[1,2]. The items for water quality include the diverse features such as dissolved oxygen (DO), chemical oxygen demand (COD), biochemical oxygen demand (BOD₅), total organic carbon (TOC), total phosphorus (T-P), total nitrogen (T-N), suspended solids (SS), turbidity (TU), potential of Hydrogen (pH), electrical conductivity (EC), water temperature (WT), and chlorophyll-a (Chl-a) and so on. The quantitative evaluation of water quality items is significant for the transaction of integrated water resources ^[3].

The water quality items can be identified in three ways: situ-measurement class (e.g., TU, pH, EC, DO, and WT), lab-measurement class (e.g., T-P, T-N, TOC, SS, COD, and Chl-a), and incubated-measurement class (e.g., BOD₅). Among water quality items, BOD₅, the representative incubated-measurement indicator, was considered as a reference to appraise the organic pollution of waterbody by the American Public Health Association Standard Methods Committee (APHASMC) ^[4]. Also, the concentration of BOD₅ can be recommended as the necessity of DO to cut down the organic matter of fluid at specific temperature ^[5]. It can, therefore, be estimated using the quantity of oxygen used up per liter of inspected data dependent on the 5-day period at 20 Celsius (℃) ^[6], and was assessed as one of essential river water quality items for the preservation and management of eco-environmental systems ^[7].

Although different machine learning and deep learning paradigms have been implemented for estimating the incubated-measurement indicator in rivers, this article proposes the unique technique for the accurate prediction of BOD₅ concentration. Hybrid neuroscience approaches involving the diverse data preprocessing coupled with the neuroscience techniques promote the evolution of more complex models based on the higher precision of estimated problems in natural behavior ^[2,8]. The double-stage synthesis models, one of hybrid neuroscience approaches, combining the wavelet transformation (WT) and different neuroscience models were developed and implemented to boost the predictive accuracy of BOD₅ concentration in Hwangji and Toilchun stations, South Korea. The standalone models such as extreme learning machine (ELM), support vector regression (SVR), and deep echo state network (Deep ESN) were also employed for integrating and evaluating novel double-stage synthesis model's scheme clearly. The novel double-stage synthesis models (i.e., Wavelet-ELM, Wavelet-SVR, and Wavelet-Deep ESN), therefore, demonstrates the efficient and accurate estimation of highly complex and nonstationary problem in rivers. The distinguished attraction of double-stage synthesis models motivates to explore the accurate prediction of BOD₅ concentration.

To the best of our knowledge and recognition from the previous information such as published articles, documents and reports, the double-stage synthesis models in the addressed article have not been frequently implemented for predicting BOD₅ concentration among the various water quality items. This article discusses the performance of implemented models (ELM, SVR, Deep ESN, Wavelet-ELM, Wavelet-SVR, and Wavelet-Deep ESN) for predicting BOD₅ concentration. They are evaluated by utilizing three mathematical formulae (R², NS, and RMSE) and four graphical aids (Scatter diagram, boxplot, violin plot, and Taylor diagram), respectively.

The rest of addressed research is arranged as follows. A brief review of BOD₅ concentration estimation and prediction is presented in section 2. The detailed description of machine learning and deep learning paradigms are provided in section 3. Also, the wavelet transformation is discussed. In section 4, report for data available and the criteria of model assessment are provided in detail. In section 5, a case study is presented by using the standalone and double-stage synthesis models based on water quantity and quality items collected in Hwangji and Toilchun stations, South Korea. In section 6, the advantages of standalone and double-stage synthesis models using mathematical formulae and graphical aids are discussed. In the end, the conclusions are drawn up.

2.   Literature review on BOD₅ concentration estimation and prediction

Various machine learning and deep learning paradigms for the estimation and prediction issues of water quality have been extensively reported in numerous articles and documents. ^[9] developed the hybrid model utilizing SVR and firefly algorithm (FFA) for predicting water quality indicator in the Euphrates River, Iraq. They found that the SVR-FFA model could predict the water quality indicator accurately. ^[10] implemented four standalone and twelve hybrid models to predict the Iran water quality indicator. The BA-RT, one of hybrid models, provided the best performance to predict the Iran water quality indicator. ^[11] employed seven standalone and three hybrid models to predict the diverse water quality indicators in China. Results showed that the decision tree (DT), random forest (RF), and deep cascade forest (DCF) models produced the outstanding achievements to predict water quality indicators in major rivers and lakes. ^[12] reviewed the recent advances in water quality remote sensing system using 200 datasets of water quality indicators. They demonstrated that the deep learning model outperformed the other proposed models to predict water quality indicators in Midwestern United States. ^[13] investigated the ELM, RF, group method of data handling (GMDH), classification and regression tree (CART), and Bat-ELM models to predict the chlorophyll-a concentration in river and lake systems, USA. They concluded that the Bat-ELM model predicted the chlorophyll-a concentration precisely compared to other models. ^[14] proposed the deep learning models including the recurrent neural network (RNN), long-short term memory (LSTM), and gated recurrent unit (GRU) to predict the drainage water quality indicator in Southern China. They showed that the deep learning models produced better prediction compared to the multiple linear regression (MLR) and multilayer perception (MLP) models.

However, limited techniques and methods have been implemented to estimate and predict BOD₅ concentration ^{[15,16,17,18]}. ^[19] employed the regression tree (RT) and SVR models to estimate total suspended solids (TSS), total dissolved solid (TDS), COD, and BOD₅ concentration using the datasets from National Stormwater Quality Database (NSQD), USA. Results showed that the applied models could estimate BOD₅ concentration accurately. ^[20] developed the adaptive neuro-fuzzy inference system (ANFIS) and wavelet SVR (WSVR) models to predict BOD₅ concentration in Karun River, Iran. They demonstrated that the WSVR model provided better prediction compared to the ANFIS model. ^[1] estimated BOD₅ concentration employing the RF, gradient boosting regression tree (GBRT), ELM, and Deep ESN in the Han River, South Korea. It can be found from ^[1]'s article that the Deep ESN5 model supplied the most accurate predictions of BOD₅ concentration among the developed models. Also, ^[2] developed two-stage and standalone neuroscience models to predict BOD₅ concentration in the Nakdong River, South Korea. Considering the developed models, the DWT-RF5 and DWT-GRNN4 models were the best model for predicting BOD₅ concentration. ^[21] utilized the SVR, RF, artificial neural networks (ANNs), long short-term memory (LSTM), convolutional neural networks (CNN)-LSTM, and Bi-LSTM models for forecasting COD and BOD₅ concentrations in the Yamuna River, India. This investigation provided that the Bi-LSTM model supplied the best performance for forecasting COD and BOD₅ concentrations. ^[22] implemented four standalone (ANN, RF, support vector machines (SVMs), and gradient boosting machines (GBM)) and six hybrid (RF-SVMs, ANN-SVMs, GBM-SVMs, RF-ANNs, GBM-ANNs, and RF-GBM) models to predict BOD₅ concentration in the Buriganga River, Bangladesh. They found that the RF-SVMs model provided the best predictive accuracy among the developed models. In contrast with the above-mentioned machine learning and deep learning paradigms, the novel double-stage synthesis models were introduced to find the optimal models between BOD₅ concentration and well-known water quality items based on five input associations. The addressed research can highlight how the novel double-stage synthesis models enhance the predictive results of BOD₅ concentration.

3.   Implemented models and supplementary method

The implemented models in the addressed article were machine learning (ELM and SVR) and deep learning (Deep ESN) paradigms, and the supplementary method was classified as the wavelet transformation, which is one of data preprocessing techniques used in various research fields. It can be seen from Figure 1 that the comprehensive mechanism of research process is underlined. Successive sub-phases explain the implemented models and supplementary method.

3.1.   Extreme Learning Machine (ELM)

^[23] initially recommended the ELM model as a rapid and effective category of feedforward neural networks (FFNN) (refer to Figure 2). It involves a single-middle-layer, which receives a particular scheme for training the parameters of networks compared to the conventional multilayer perceptron (MLP) model. The ELM model can map using a single-middle-layer with M independent input indicators and be written as follows:

where g(.) is the activation function, which supplies the output in the middle layer; ${\beta _i}$ is the weight of output for connecting the middle neurons to the output neuron; and L is the neuron number of the middle layer. The output indicator can be given by the following formula (2):

where $\varepsilon$ is the error. The Gaussian and sigmoid functions are the most employed mapping ones in the ELM model's category. The underlying formula (3) expresses the Gaussian function:

where a and c refer to the activation functions. During training phase, the connection weight is fixed in the ELM model's category. That is, random values are allowed directly to neurons' activation functions instead of requesting an iterative process to update them. The connection weights for output neuron can be achieved continuously utilizing the least squares method. In other words, the fitting error ought to be reduced by computing ${\left\| {{\bf{H \pmb{\mathsf{ β}} }} - {\bf{T}}} \right\|^2}$ for the connection weight (${\bf{ \pmb{\mathsf{ β}} }}$), where T is the matrix for target and H is the randomized matrix corresponding to the middle layer:

The connection weight for output is resolved, based on the linear equation system such as ${\bf{ \pmb{\mathsf{ β}} }}{\text{ = }}{{\bf{{\rm H}}}^{\text{ + }}}{\bf{T}}$, where ${{\bf{{\rm H}}}^{\text{ + }}}$ is the generalized inverse function of Moore-Penrose ^[1,24].

3.2.   Support Vector Regression (SVR)

The SVR model (refer to Figure 3), which is a special type of SVMs, has been applied in various fields, including stock index prediction, bioinformatics engineering, chemical synthesis, and production process control and so on ^[25,26]. The generalization of conventional ANNs models may reduce to a local optimized generalization, while a universal optimization is insured for the SVR model ^[27,28,29].

The fundamental concepts of SVR model are as follows. Recognizing the training sample $\left({{{\text{x}}_{\text{i}}}{\text{, }}{{\text{y}}_{\text{i}}}} \right)$, where ${{\text{x}}_{\text{i}}} \in {\Re ^n}$ is a specific value of input indicator x, and ${{\text{y}}_{\text{i}}} \in {\Re ^n}$ is the matching value of surveyed model output. Also, a nonlinear transfer function ($\Phi (\cdot)$) and a linear function (${\text{f(}} \cdot {\text{)}}$) can be defined between input and output indicators. The actual output, therefore, is expressed by formula (5):

where $\overline {\text{y}}$ is the actual output; and w and B are the adjustable parameters of the model. In the SVR model, the empirical risk can be written as the following formula (6):

where R_emp is the empirical risk; and ${\left| {{{\rm{y_i}}} - {{\overline {\rm{y}} }_{\rm{i}}}} \right|_{\rm{ \mathsf{ ε} }} }$ is Vapnik's ε-insensitive loss function. The adjustable parameters (i.e., w and B in formula (5)) of the model can be calculated by obtaining the minimum cost function ^[27]. In the addressed article, the following cost function was used:

where ${{\rm{ \mathsf{ψ} }}_{\rm{ \mathsf{ ε} }}}{\text{(w}}, {\rm{ \mathsf{ ξ} }}, {{\rm{ \mathsf{ ξ} }}^ * }{\text{)}}$ is the cost function; ${{\rm{ \mathsf{ ξ} }}_{\text{i}}}{{, {\rm{ \mathsf{ ξ} }}}}_{\text{i}}^ *$ are the positive slack variables; and C is the cost constant. In addition, the constraints for formula (7) can be classified as: (1) ${{\text{y}}_{\text{i}}} - {\overline {\text{y}} _{\text{i}}} \leqslant {\rm{ \mathsf{ ε} }} + {{\rm{ \mathsf{ ξ} }}_i}$ i = 1, 2, ..., N; (2) $- {{\text{y}}_{\text{i}}} + {\overline {\text{y}} _{\text{i}}} \leqslant {\rm{ \mathsf{ ε} }} + {\rm{ \mathsf{ ξ} }}_{\text{i}}^ *$ i = 1, 2, ..., N; and (3) ${{\rm{ \mathsf{ ξ} }}_{\text{i}}}{{, {\rm{ \mathsf{ ξ} }}}}_{\text{i}}^ * \geqslant {\text{0}}$ i = 1, 2, ..., N.

3.3.   Deep Echo State Network (Deep ESN)

The recurrent neural networks (RNN), including the echo state network (ESN), has the most broadly utilized reservoir computing (RC) method ^[1,31,32]. Since the RNN model is powerful and accurate for computing the complicated and nonlinear problems, the Deep ESN model is effective for historical data. The Deep ESN model contains an order of deformed recurrent layers named as reservoir, where every layer output performs as the following layer input. Figure 4 explains a conceptual diagram of Deep ESN model with the recurrent structure of reservoir from the viewpoint of a discrete-time dynamic system. The reservoir dynamic state can be renewed by recognizing the leaky integration ESN (i.e., LI-ESN) as below:

where α is the leaky coefficient; u(t) is the outside input based on time t; x(t) is the reservoir state in the corresponding layer based on time t; W_in is the matrix of input connection weight for the reservoir; and W_R is the matrix of recurrent connection weight. Because the input indicator to the following reservoir can be supplied by the output indicator of its prior reservoir, an ordinary equation for the Deep ESN model is organized for expanding the function of state transition as below.

where $\ell$ is a layer (reservoir) in the structure of RC; ${{\text{W}}^\ell }_{\text{R}}$ are the connection weights between the layer $\ell$ and the prior one $\ell - 1$; and L is the quantity for the layers of reservoir ^[1,32,33]. Here, five layers were employed in the reservoir for the Deep ESN model to estimate BOD₅ concentration.

3.4.   Wavelet Transformation (WT)

The WT method, which is one of multi-resolution signal procedure methods, is employed to build the double-stage synthesis models. The original data can be isolated into various frequency components involving an approximation and numerous details handling the WT method. In the addressed article, discrete-based wavelet transformation, which has been utilized for the data preprocessing in diverse fields, was selected. In fact, the discrete-based WT method can be accomplished by implementing the Mallat method ^[34]. The bottom line of Mallat method is two-route filters comprising two filters as low-pass and high-pass ^[2,35]. ^[36] defined that the coefficients for the wavelet (high-pass) and scaling (low-pass) in the j^th level of decomposition is outlined as

where ${W_{j, \; t}}$ and ${V_{j, \; t}}$ are the elements for corresponding ${{\bf{W}}_j}$ and ${{\bf{V}}_j}$. The WT decomposes the complex and original input time series into the components (approximation and details) which show relatively simpler patterns than the original input time series. The different components obtained from WT were implemented as input association of corresponding double-stage synthesis model. Evolving double-stage synthesis models for the components separately and summing their predicted values can improve the predictive accuracy of double-stage synthesis models compared to performance of standalone models for the original input time series with high complexity. A flowchart for dual-step discrete-based WT is shown in figure 5. Here, two details (D₁ and D₂) and an approximation (A₂) are achieved from the original input time series. Also, Figure 6 illustrates the sequential diagram for evolving the double-stage synthesis models.

4.   Report for data and assessment criteria

4.1.   Preparation of utilized data

The original data can be isolated into various frequency In the addressed article, Hwangji (longitude 129°05′07″E; latitude 37°06′74″N) and Toilchun (longitude 128°44′46″E; latitude 36°47′09″N) stations were employed to predict BOD₅ concentration using diverse physical and chemical variables such as T-N, T-P, TOC, DO, WT, SS, COD, pH, EC, and station discharge (DIS) in South Korea. Figure 7 shows the illustrative map of Hwangji and Toilchun stations.

The surveyed data (2008/02-2020/12 for Hwangji and 2011/07-2020/12 for Toilchun stations) for water quantity and quality items can be directly accessed and collected from official website (http://water.nier.go.kr) of National Institute of Environmental Research (NIER), South Korea. The full data file consisted of training and validation samples. The training sample involved 80% (data = 398 from Hwangji and data = 294 from Toilchun stations) and the validation sample applied the last 20% (data = 99 from Hwangji and data = 74 from Toilchun stations) of full data file.

Recognizing the source code and software of machine learning and deep learning paradigms, the ELM model was evolved by employing the R (https://www.r-project.org, a free software environment for statistical computing and graphics) package and the elmNNRcpp (https://cran.r-project.org/web/packages/elmNNRcpp/index.html). In case of SVR model, it was implemented by the DTREG predictive modeling software (www.dtreg.com). In addition, the Deep ESN model was developed by utilizing the MATLAB programing language (https://www.mathworks.com), which is a freely available MATLAB toolbox for the Deep ESN (https://it.mathworks.com/matlabcentral/fileexchange/69402-deepesn).

The optimal number of hidden nodes for the ELM and Wavelet-ELM models was determined using a trial and error process. As the number of hidden nodes was changed from 1 to 5m (where, m is the number of input indicators), the number of hidden nodes with the minimum RMSE value was decided as the optimal value. The logistic sigmoid function and linear function were used for activating hidden and output nodes, respectively. In addition, epsilon type of SVR model with radial basis function (RBF) kernel was employed for predicting BOD₅ concentration using the SVR and Wavelet-SVR models. Also, the V-fold cross-validation were applied to validate the SVR and Wavelet-SVR models, and the grid search algorithm found the optimal parameters by minimizing total errors. Finally, the optimal number of layers and reservoirs units were decided based on the trial and error process for the Deep ESN and Wavelet-Deep ESN models.

Table 1 explains the computed results for the correlation coefficients and P values between individual input indicators and BOD₅ concentration. It can be judged from Table 1 that TOC (e.g., CC = 0.721, P value = 0.0001 at Hwangji and CC = 0.563, P value = 0.0001 at Toilchun stations) and COD (e.g., CC = 0.721, P value = 0.0001 at Hwangji and CC = 0.626, P value = 0.0001 at Toilchun stations) items exhibited high correlation and statistically significant with BOD₅ concentration among various input indicators. In the addressed article, all indicators can be categorized as class 1 (i.e., in situ-measurement items (pH, EC, DO, and WT), class 2 (i.e., lab-measurement items (SS, COD, T-N, T-P, and TOC) and incubated-measurement item (BOD₅)), and class 3 (i.e., water quantity item (river discharge)), respectively.

4.2.   Mathematical assessment criteria of standalone and double-stage synthesis models

To assess the performance of standalone (ELM, SVR, and Deep ESN) and double-stage synthesis (Wavelet-ELM, Wavelet-SVR, and Wavelet-Deep ESN) models, three mathematical formulae, which have been recognized and utilized worldwide, were employed. The coefficient of determination (R²) criterion ^[37,38] is clarified as the square of correlation between surveyed and estimated BOD₅ concentrations (see formula (11)). The Nash-Sutcliffe (NS) efficiency criterion ^[39] can resolve the models' effectiveness between surveyed and estimated BOD₅ concentrations (see formula (12)). Also, the disparity between surveyed and estimated BOD₅ concentrations can be referred by handling the root mean square error (RMSE) criterion ^[40]. The RMSE criterion can be computed by employing formula (13).

where $BO{D_{sur}}$ and $BO{D_{est}}$ = surveyed and estimated BOD₅ concentrations; ${\overline {BOD} _{sur}}$ and $\overline B O{D_{est}}$ = surveyed and estimated mean BOD₅ concentrations; and n = the number of full data available.

5.   Case study

The addressed article utilized the miscellaneous water quantity and quality items to predict BOD₅ concentration in Hwangji and Toilchun stations, South Korea. As defined formerly, the assessment of standalone and double-stage synthesis models to estimate BOD₅ concentration was the essential view of this article.

Among water quality items, some indicators, including pH, EC, DO, and WT, were directly surveyed by utilizing a commercial mechanical tool. Also, the indicators including SS, COD, T-N, T-P, and TOC were surveyed through the laboratory assistant system indirectly. BOD₅ concentration, however, can be indirectly surveyed via the incubation system, based on the 20 ℃ for the 5-day period ^[41]. Since the plan of addressed article was scheduled to predict BOD₅ concentration utilizing the standalone and double-stage synthesis models, this behavior could save and protect the time and effort to estimate and incubate BOD₅ concentration.

From the correlation coefficients and P values of water quantity and quality items (refer to Table 1), divergent organizations were provided to choose the best input association for given stations. To employ the same input indicators on both stations among them, some input indicators with positive (e.g., COD, TOC, T-P, and SS) and negative (e.g., WT, pH, and DIS) correlations were selected for diverse input associations in Hwangji station because input indicators with negative correlations can also contribute for predicting BOD₅ concentration. However, many input indicators based on positive (e.g., COD, TOC, T-P, SS, WT, pH, and DIS) correlation were implemented for different input associations in Toilchun station.

Hence, the standalone and double-stage synthesis models were evolved for predicting BOD₅ concentration, based on five input associations (so called, 1^st–5^th distributions). Because TOC and COD items were picked out as the underlying water quality items for given stations, the addressed article determined the consolidation of TOC and COD items as the 1^st distribution. Table 2 presents the diverse input associations of water quantity and quality items to predict BOD₅ concentration. All developed models in Table 2 can be categorized into five distributions.

5.1.   Predicting BOD₅ concentration at Hwangji station

5.1.1.   Application of standalone models

The results of three mathematical formulae (R², NS, and RMSE) for the standalone models are summed up in Table 3 for Hwangji station. Table 3 explains that the outcomes of SVR1 (R² = 0.905, NS = 0.891, and RMSE = 0.299 mg/L) are more excellent than the ELM1 and Deep ESN1 in the validation phase dependent on the 1^st distribution. In the 2^nd distribution, the SVR2 (R² = 0.908, NS = 0.905, and RMSE = 0.279 mg/L) performs more excellent than the ELM2 and Deep ESN2. And, the SVR3 (R² = 0.925, NS = 0.915, and RMSE = 0.264 mg/L) surpasses the ELM3 and Deep ESN3 clearly in the validation phase for the 3^rd distribution. The contrast of standalone models in the 4^th distribution, furthermore, indicates that the ELM4 (R² = 0.902, NS = 0.893, and RMSE = 0.295 mg/L) dominates the SVR4 and Deep ESN4 in the validation phase. In the end, the SVR5 (R² = 0.905, NS = 0.884, and RMSE = 0.309 mg/L) is more accurate than ELM5 and Deep ESN5 in the validation phase for the 5^th distribution.

Recognizing the impressive models from the 1^st–5^th distributions, the best accomplishment of standalone models can be found from the ELM (the 4^th distribution), SVR (the 3^rd distribution), and Deep ESN (the 2^nd distribution) among diverse input associations in the validation phase. Table 3 tells us that the desirable performance of SVR3 gives more accurate than ELM4 and Deep ESN2 in the validation phase. Therefore, it can be said that the SVR3 is the most accurate for predicting BOD₅ concentration among the desirable standalone models at Hwangji station.

5.1.2.   Application of double-stage synthesis models

The results of three mathematical criteria for the double-stage synthesis models are also arranged in Table 4 at Hwangji station. From Table 4, it is clear that the outcomes of Wavelet-SVR1 (R² = 0.904, NS = 0.895, and RMSE = 0.293 mg/L) are more dominant compared to the Wavelet-ELM1 and Wavelet-Deep ESN1 in the validation phase dependent on the 1^st distribution. Based on the 2^nd distribution, the Wavelet-SVR2 (R² = 0.911, NS = 0.911, and RMSE = 0.271 mg/L) is more excellent than the Wavelet-ELM2 and Wavelet-Deep ESN2. Also, the Wavelet-SVR3 (R² = 0.920, NS = 0.912, and RMSE = 0.269 mg/L) outperforms the Wavelet-ELM3 and Wavelet-Deep ESN3 regarding the 3^rd distribution in the validation phase. Moreover, the contrast of double-stage synthesis models in the 4^th distribution demonstrates that the Wavelet-SVR4 (R² = 0.926, NS = 0.915, and RMSE = 0.264 mg/L) performs superior to the Wavelet-ELM4 and Wavelet-Deep ESN4 in the validation phase. In the end, Wavelet-SVR4 (R² = 0.919, NS = 0.914, and RMSE = 0.266 mg/L) is more efficient than the Wavelet-ELM5 and Wavelet-Deep ESN5 in the validation phase for the 5^th distribution.

Contemplating the outstanding models from the 1^st-5^th distributions, the admirable performance of double-stage synthesis models can be judged from the Wavelet-ELM (the 2^nd distribution), Wavelet-SVR (the 4^th distribution), and Wavelet-Deep ESN (the 3^rd distribution) among diverse input associations in the validation phase. It can be noticed from Table 4 that the Wavelet-SVR4 provides more effective outcomes than the Wavelet-ELM2 and Wavelet-Deep ESN3 in the validation phase. For that reason, the Wavelet-SVR4 is more trustworthy than the Wavelet-ELM2 and Wavelet-Deep ESN3 for predicting BOD₅ concentration among the desirable double-stage synthesis models at Hwangji station.

5.1.3.   Graphical aids of model accomplishment

To verify the accuracy of desirable standalone and double-stage synthesis models using graphical aids, Figures 8(a)–(f) present the scatterplots for the surveyed and estimated BOD₅ concentration values at Hwangji station. The linear functions and values of NS efficiency criterion are presented for the corresponding standalone and double-stage synthesis models, respectively. It can be concluded from NS values and the slopes of linear functions that an apparent discrepancy can be followed among the desirable standalone and double-stage synthesis models (ELM4, SVR3, Deep ESN2, Wavelet-ELM2, Wavelet-SVR4, and Wavelet-Deep ESN3). Therefore, the SVR3 and Wavelt-SVR4 performs the most reliable accuracy for predicting BOD₅ concentration values clearly, whereas the Wavelet-ELM2 was the worst among the desirable models at Hwangji station.

Additional portraits can evaluate the performance of standalone and double-stage synthesis models using the boxplot, violin plot ^[42], and Taylor diagram ^[43]. Figures 9(a)–(c) show the diverse graphical aids for the desirable standalone and double-stage synthesis models at Hwangji station. It can be found from Figure 9(a) that the estimated BOD₅ concentrations of SVR3 and Wavelet-SVR4 yield more analogous configuration to the surveyed values for median, interquartile ranges and dispersion, adjacent values, and sign of skewness compared to other desirable models. Another graphical aid for the distribution of surveyed and estimated BOD₅ concentration values utilizing the desirable models can be provided with the violin plots (Figure 9(b)). The violin plot can be defined as one of approaches to discern the distribution of assigned numerical values. Figure 9(b) supplies a close shape pattern for the SVR3 and Wavelet-SVR4 concerning the median, interquartile, and distribution of assigned values. In addition, the Taylor diagram (Figure 9(c)) utilizes three statistical indices, including correlation coefficient, normalized standard deviation, and root mean square error. The principal approach of Taylor diagram can be explained as to find the closest estimated model with the corresponding surveyed BOD₅ concentration based on standard deviation (polar axis) and correlation coefficient (radial axis). The Taylor diagram, therefore, demonstrates the accuracy and efficiency of SVR3 and Wavelet-SVR4 over the other desirable models (ELM4, Deep ESN2, Wavelet-ELM2, and Wavelet-Deep ESN3).

.

5.2.   Predicting BOD₅ concentration at Toilchun station

5.2.1.   Utilization of standalone models

The outputs of three mathematical formulae (R², NS, and RMSE) for the standalone models are summed up in Table 5 for Toilchun station. Table 5 shows that the estimates of ELM1 (R² = 0.571, NS = 0.472, and RMSE = 0.472 mg/L) are preferable to the SVM1 and Deep ESN1 in the validation phase considering the 1^st distribution. During the 2^nd distribution, the SVR2 model (R² = 0.722, NS = 0.691, and RMSE = 0.361 mg/L) are more remarkable compared to the ELM2 and Deep ESN2. Also, the performance of SVR3 (R² = 0.701, NS = 0.661, and RMSE = 0.378 mg/L) exceeds the ELM3 and Deep ESN3 obviously regarding the 3^rd distribution in the validation phase. The contradiction of standalone models subjected to the 4^th distribution, besides, approves that the SVR4 (R² = 0.868, NS = 0.854, and RMSE = 0.248 mg/L) outperforms the ELM4 and Deep ESN4 in the validation phase. Eventually, the SVR5 (R² = 0.876, NS = 0.870, and RMSE = 0.234 mg/L) is more reliable than ELM5 and Deep ESN5 in the validation phase with the 5^th distribution.

Contemplating the magnificent models among the 1^st–5^th distributions, the best capability of standalone models can be discovered from the ELM (the 5^th distribution), SVR (the 5^th distribution), and Deep ESN (the 5^th distribution) among diverse input associations in the validation phase. As seen from Table 5, the improved performance of SVR5 contributes better prediction compared to the ELM5 and Deep ESN5 in the validation phase. As a result, the SVR5 is most reliable for predicting BOD₅ concentration among the improved standalone models at Toilchun station.

5.2.2.   Utilization of double-stage synthesis models

The outputs of three mathematical criteria for the double-stage synthesis models are still organized as in Table 6 at Toilchun station. As observed from Table 6, the estimates of Wavelet-SVR1 (R² = 0.662, NS = 0.647, and RMSE = 0.386 mg/L) are more prevalent than the Wavelet-ELM1 and Wavelet-Deep ESN1 in the validation phase utilizing the 1^st distribution. Favoring the 2^nd distribution, the Wavelet-SVR2 (R² = 0.866, NS = 0.845, and RMSE = 0.255 mg/L) is more exquisite than the Wavelet-ELM2 and Wavelet-Deep ESN2. Likewise, the performance of Wavelet-SVR3 (R² = 0.688, NS = 0.646, and RMSE = 0.386 mg/L) exceeds the Wavelet-ELM3 and Wavelet-Deep ESN3 viewing the 3^rd distribution during validation phase. Likewise, the contradiction of double-stage synthesis models in the 4^th distribution demonstrates that the Wavelet-SVR4 (R² = 0.922, NS = 0.917, and RMSE = 0.187 mg/L) surpasses the Wavelet-ELM4 and Wavelet-Deep ESN4 definitely in the validation phase. Eventually, the Wavelet-SVR5 (R² = 0.780, NS = 0.775, and RMSE = 0.308 mg/L) is more effective than the Wavelet-ELM5 and Wavelet-Deep ESN5 in the validation phase with the 5^th distribution.

Envisaging the eminent models from the 1^st–5^th distributions, the attractive performance of double-stage synthesis models can be evaluated from the Wavelet-ELM (the 5^th distribution), Wavelet-SVR (the 4^th distribution), and Wavelet-Deep ESN (the 5^th distribution) among various input associations in the validation phase. It can be seen from Table 6 that the Wavelet-SVR4 yields more reliable outcomes compared to the Wavelet-ELM5 and Wavelet-Deep ESN5 in the validation phase. As a result, the Wavelet-SVR4 performs superior to the Wavelet-ELM5 and Wavelet-Deep ESN5 for predicting BOD₅ concentration among the enhanced double-stage synthesis models at Toilchun station.

5.2.3.   Visual aids of model accomplishment

To validate the precision of desirable standalone and double-stage synthesis models using visual aids, Figures 10(a)–(f) provide the scatterplots for the surveyed and estimated BOD₅ concentration values employing the desirable standalone and double-stage synthesis models at Toilchun station. The linear formulae and values of NS efficiency criterion are inserted for the corresponding standalone and double-stage synthesis models, respectively. It can be inferred from NS values and the slopes of linear formulae that a definite inconsistency can be traced among the desirable standalone and double-stage synthesis models (ELM5, SVR5, Deep ESN5, Wavelet-ELM5, Wavelet-SVR4, and Wavelet-Deep ESN5). Therefore, the Wavelet-SVR4 accomplishes the most reliable precision for predicting BOD₅ concentration values obviously, while the Deep ESN5 yields the least precise among the desirable models at Toilchun station.

Additional pictures can anticipate the performance of standalone and double-stage synthesis models using the boxplot, violin plot, and Taylor diagram. Figures 11(a)–(c) support the various visual aids for the desirable standalone and double-stage synthesis models at Toilchun station. It can be seen from Figure 11(a) that the boxplots of estimated BOD₅ concentration utilizing the Wavelet-SVR4 can resemble that of surveyed BOD₅ concentration intimately. Another visual aid for the distribution of surveyed and estimated BOD₅ concentration values employing the desirable models can be displayed with the violin plots (Figure 11(b)). The violin plot can be described as one of schematic techniques to reveal the distribution of mandated numerical values. Figure 11(b) demonstrates a similar contour pattern for the Wavelet-SVR4 employing the median, interquartile, and distribution of mandated values. Figure 11(c) supplies the Taylor diagram employing the desirable standalone and double-stage synthesis models for Toilchun station. It can be seen from Figure 11(c) that the point of Wavelet-SVR4 which has the minimal RMSE value visualizes the straight distance from the surveyed one, while the point of Deep ESN5 displays the longest distance from the surveyed point.

6.   Discussion

The addressed research explored the nonlinear behavior (e.g., hard to predict) of BOD₅ concentration employing standalone and double-stage synthesis models in Hwangji and Toilchun stations, South Korea. Since both stations (i.e., Hwangji and Toilchun) yielded the different high-quality accuracies for the desirable standalone models, it was hard to judge which model predicted BOD₅ concentration with accuracy. Also, the outputs of three mathematical formulae explained that the SVM models with diverse input associations could predict BOD₅ concentration precisely compared to the ELM and Deep ESN models based on the corresponding distribution on both stations. Because all standalone models enforced the various theoretical structures and inference, the accurate prediction was changed for diverse input associations of standalone models.

The main aim for developing the double-stage synthesis models was to enhance the accurate prediction of BOD₅ concentration compared to corresponding standalone models. Unfortunately, the Wavelet-ELM models could not boost the accurate prediction for corresponding ELM models from the perspective of double-stage synthesis models' performance, based on NS values at Hwangji station. Among the Wavelet-SVR models, the Wavelet-SVR1 (0.4% for SVR1), Wavelet-SVR2 (0.7% for SVR2), Wavelet-SVR4 (3.7% for SVR4), and Wavelet-SVR5 (3.4% for SVR5) models slightly enhanced the accurate prediction. Also, all the Wavelet-Deep ESN models increased the predictive accuracy on a small scale, including the Wavelet-Deep ESN1 (2.5% for Deep ESN1), Wavelet-Deep ESN2 (0.1% for Deep ESN2), Wavelet-Deep ESN3 (3.5% for Deep ESN3), Wavelet-Deep ESN4 (6.0% for Deep ESN4), and Wavelet-Deep ESN5 (1.0% for Deep ESN5). Noticing the desirable models' categorization for the standalone and double-stage synthesis models, the Wavelet-SVR4 model, which yielded the best accuracy, improved the accurate prediction by 12.7% (Wavelet-ELM2), 9.8% (Wavelet-Deep ESN3), 2.5% (ELM4), and 10.1% (Deep ESN2), respectively.

Regarding the performance evaluation of double-stage synthesis models by NS values at Toilchun station, only the Wavelet-ELM1 (19.7% for ELM1) model could boost the estimated efficiency obviously among the Wavelet-ELM models. Also, Wavelet-SVR1 (35.6% for SVR1), Wavelet-SVR2 (22.3% for SVR2), and Wavelet-SVR4 (7.4% for SVR4) models increased the accurate prediction clearly among the Wavelet-SVR models. In addition, all the Wavelet-Deep ESN models, including the Wavelet-Deep ESN1 (31.2% for Deep ESN1), Wavelet-Deep ESN2 (0.5% for Deep ESN2), Wavelet-Deep ESN3 (36.7% for Deep ESN3), Wavelet-Deep ESN4 (3.5% for Deep ESN4), and Wavelet-Deep ESN5 (8.9% for Deep ESN5) models, boosted the precise efficiency, respectively. Considering the desirable models' categorization for the standalone and double-stage synthesis models, the Wavelet-SVR4 model which produced the best accuracy, reinforced the accurate prediction by 32.8% (Wavelet-ELM5), 38.6% (Wavelet-Deep ESN5), 13.4% (ELM5), 5.2% (SVR5), and 51.0% (Deep ESN5), respectively. The double-stage synthesis models, therefore, could not always reinforce the accurate prediction of corresponding standalone models on both stations. This experience pursued the previous works of ^[2] and ^[44]. ^[44] predicted DO concentration employing the single and hybrid machine learning models in Florida, USA. Results demonstrated that the hybrid machine learning models could not regularly improve the predicted accuracy of single machine learning models. Also, ^[2] implemented the single and combinational paradigm to predict BOD₅ concentration in South Korea. They found that the combinational paradigm could not always increase the predictive accuracy of single models clearly.

Therefore, the process which embeds the different data preprocessing algorithms ^{[45,46,47,48,49]} in the diverse standalone (i.e., machine learning and deep learning) models, is required to increase the accurate prediction and efficiency of BOD₅ concentration for the continuous research.

7.   Conclusions

The addressed research explored the precision and efficiency of the standalone and double-stage synthesis models for predicting BOD₅ concentration in Hwangji and Toilchun stations, South Korea. Five input associations (1^st–5^th distributions) were resolved for developing the standalone and double-stage synthesis models based on seven water quantity and quality items. For the modeling and prediction of standalone and double-stage synthesis models, the assembled data were divided into training and validation samples, respectively. Three mathematical formulae (R², NS, and RMSE) and four graphical aids (scatter diagram, boxplot, violin plot, and Taylor diagram) were used to evaluate the accurate prediction of addressed models.

Considering the best models from the 1^st-5^th distributions, the SVR3 (R² = 0.925, NS = 0.915, and RMSE = 0.264 mg/L) and Wavelet-SVR4 (R² = 0.926, NS = 0.915, and RMSE = 0.264 mg/L) models were the most precise compared to other desirable models (ELM4, Deep ESN2, Wavelet-ELM2, and Wavelet-Deep ESN3) based on standalone and double-stage synthesis models in the validation phase at Hwangji station. Also, the Wavelet-SVR4 model (R² = 0.922, NS = 0.917, and RMSE = 0.162 mg/L) provided more precise results than other desirable models (ELM5, SVR5, Deep ESN5, Wavelet-ELM5, and Wavelet-Deep ESN5) for predicting BOD₅ concentration in the validation phase at Toilchun station. However, it was found the addressed research that explained that the precision and efficiency of BOD₅ concentration estimated by the standalone models could not be reinforced by the double-stage synthesis models on both stations. Therefore, using the credible water quantity and quality items from the available data groups can confirm the outputs of the addressed research, and perform the best prediction of BOD₅ concentration employing the different standalone and double-stage synthesis models in river.

Conflict of interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.