Bearing capacity prediction of the concrete pile using tunned ANFIS system

The design process for pile foundations necessitates meticulous deliberation of the calculation pertaining to the bearing capacity of the piles. The primary objective of this work was to investigate the potential use of Coot bird optimization ( CBO ) techniques in predicting the load-bearing capacity of concrete-driven piles. Despite the availability of several suggested models, the investigation of Coot bird optimization ( CBO ) for estimating the pile-carrying capacity has been somewhat neglected in this research. This work presents and validates a unique approach that combines the Coot bird optimization ( CBO ) model with the Multi-layered perceptron ( MLP ) neural network and adaptive neuro-fuzzy inference system ( ANFIS ). The findings of 472 different driven pile static load tests were put in a database. The proposed framework’s building, validation, and testing stages were each accomplished utilizing the training set (70%), validation set (15%), and testing set (15%) of the dataset, respectively. According to the findings, MLP CBO and ANFIS CBO both offer remarkable possibilities for accurately predicting the pile-bearing capacity of a given structure. The R 2 values for ANFIS CBO during the training stage were 0.9874, while during the validating stage, they were 0.9785, and during the testing stage they were 0.987. After considering various kinds of performance studies and contrasting them with existing literature, it has been concluded that the ANFIS CBO model provides a more appropriate calculation of the bearing capacity of concrete-driven piles.


Introduction
To sustain superstructures with enormous loads or those situated on unstable ground, deep foundations are a typical and essential form of foundation [1].Driven piles composed of wood, steel, precast concrete, and composite are another cost-and quality-efficient alternative to drilled shafts.The axial pile-bearing capability is considered the most significant factor when designing a pile foundation.As a result, many theoretical and empirical geotechnical research has focused on calculating this value.
The pile-bearing capability may be assessed using five primary techniques: static analysis, dynamic analysis, dynamic testing, pile load test, and in-situ testing [2][3][4][5].It is often suggested to use the essential depth idea in design recommendations according to static analysis.The crucial depth, nevertheless, is an idealization that defies physical principles and lacks theoretical or trustworthy empirical evidence.
The hammer-pile-soil system is subject to dynamic analysis using wave mechanics as a foundation.Uncertainties in calculating bearing capability are brought on by the hammer impact effect's ambiguity and also by variations in soil power due to environmental factors at the time of pile drive and loading.The measurement of velocity and force close to the pile head during driving is the foundation of dynamic testing techniques.Nevertheless, a skilled individual is able to interpret the evaluations.The capability estimate is also not known until the pile is driven, which is a significant restriction [6].The full-scale pile settlement under a static load is measured in the field in a pile load test, which is thought to convey the most precise findings.This approach, nevertheless takes a lot of time and money [7].So, it is crucial to create a straightforward, cost-effective, and precise procedure.
Since the 1970s, in situ test techniques for measuring soil parameters have progressed quickly.This advancement is being accompanied by an increase in the estimation of pilebearing capability using data from in-situ tests.Standard penetration test ( SPT ), cone pen- etration test ( CPT ), flat dilatometer test ( DMT ), pressuremeter test ( PMT ), plate loading test ( PLT ), dynamic probing test ( DP ), press-in and screw-on probe test ( SS ), and field vane test ( FVT ) are a few examples of typical tests.In an effort to assess the material properties, every test employs several loading techniques to quantify the related soil reaction.The SPT is often utilized to forecast the piles' bearing capability between these in-situ test data [8,9].
In the research, many SPT data-based approaches for calculating pile-bearing capaci- ties have been developed.Direct and indirect techniques are the two categories into which they fall.Because of their simplicity in calculation, the direct techniques are preferred by field engineers.For instance, SPT direct approaches for sandy or clayed soil were presented in [10][11][12][13][14][15].The researchers used the limited component approach to assess the pile for a case investigation in Iran [15] and contrasted the results with four distinct SPT direct techniques to arrive at an accurate forecast of the pile's bearing capa- bility.Nevertheless, every one of these experimental formulations have certain shortcomings, based on [4].As a result, scholars have started searching for various techniques to forecast pile-bearing capability using SPT data.Utilizing machine learning techniques may be successful, according to earlier research [16].
Artificial neural networks ( ANN ), in particular, have been widely employed in ML-based models to forecast pile-bearing capability.Previous attempts in this manner may be seen in [42,43], which use ANN with error back propagation.Utilizing results from 50 dynamic load tests performed on prefabricated concrete piles, [44] combines ANN with genetic algorithm ( GA ), with the weights of ANN being adjusted by GA .Identical methods are sug- gested in [45], where the ANN connection weights are optimized using particle swarm opti- mization ( PSO ) along with GA .When using ANN to forecast the piles' bearing capability, GA is also utilized to identify the most crucial characteristics in the unprocessed dataset [3].Other ML-based approaches than ANN have also been taken into account; for example, Samui [46] employs support vector machine ( SVM ), Pham et al. [2] explore random forest, and Chen et al. [38] investigated the neuro-imperialism and neuro-genetic approaches.
Many studies employed a number of machine learning techniques to successfully measure the foundation's bearing capabilities.For the purpose of predicting the bearing capacity of the thin-walled foundation, ANN models were developed utilizing a total of 150 specimens, and an adaptive neuro-fuzzy inference system ( ANFIS ) was developed as well [47].The researchers obtained the coefficient of determination ( R 2 ) and root mean squared error ( RMSE ) values for the ANFIS and ANN models, which were determined to be 0.875 and 0.048, and 0.71 and 0.529, accordingly.In an alternative study, the determination of pile bearing capability included the use of 296 sets of data and the application of Gaussian process regression analysis, resulting in an R 2 value of 0.84 [48].In their study, Kulkarni et al. [49] used a dataset consisting of 132 data points to elucidate the behavior of rock-socketed piles.This was achieved by using a hybrid approach that combined the genetic algorithm ( GA ) technique with artificial neural networks ( ANN ).The researchers obtained a coefficient of determination ( R 2 ) value of 0.86 and a root mean square error ( RMSE ) value of 0.0093.Using this dataset and a hybrid model known as particle swarm optimization ( PSO)-ANN , another group of researchers were able to replicate the behavior of rock-socketed piles [50].Their results showed an R 2 value of 0.918 and an RMSE value of 0.063.The dataset used in previous empirical endeavors to construct models consisted of a total of 472 specimens [3,51].Models that were taken into account included a hybrid deep neural network ( DNN ) with genetic algorithm ( GA ), as well as a hybrid whale optimiza- tion ( WOA ) technique with extreme gradient boosting ( XGB ).The R 2 and RMSE values that were generated by the GA-DNN model came out to be 0.882 and 109.965,accordingly.The findings for WOA-XGB showed significant improvement, with R 2 , and RMSE values of 0.97 and 64 during the training step, and 0.94 and 87.03 at the testing step, respectively.
In conclusion, the present study has made some notable contributions, which may be outlined as follows: • In this study, a comprehensive dataset including 472 pile test findings is used to develop and validate machine learning models for the assessment of pile-bearing capacity.It is essential to acknowledge that prior research has often concentrated on limited datasets.• The investigation of coot bird optimization ( CBO ) for assessing the load-bearing capacity of piles has been somewhat neglected in scholarly research, despite the existence of numerous suggested models.• The selection of the appropriate model, which entails the selection of the appropriate hyper-parameters, is of the utmost importance.This work presents and validates a unique approach that combines the coot bird optimization ( CBO ) model the multi- layered perceptron ( MLP ) neural network and the adaptive neuro-fuzzy inference system ( ANFIS).

Data description and pre-processing
The dataset used for training and evaluating the machine learning approach described in this paper consists of findings of static load tests conducted on driven strengthened concrete piles.The 472 pieces of data that are now included in this collection were obtained for previous research [3].The adequacy of the data collection for the development and validation of advanced machine-learning models is attributed to its substantial size.In order to evaluate the capacity of the precast piles, hydraulic pile-driving machinery is used to drive the piles into the soil layers.The piles have closed tops.Figure 1 illustrates the testing environment, including the data gathering apparatus, as well as the specific testing site.The bearing capacity of piles was determined using a static load test experiment detailed in Fig. 1b. .Table 1 contains an overview of all ten of the conditioning variables that were considered while making a prediction about the dependent parameter (axial pile bearing capacity).These variables are denoted by the notation In 1 through In 10 .
Along with other pieces of information, this table also provides a statistical analysis of the independent parameters as well as the parameters that are being predicted.The dataset contains 472 pieces of data, and it was partitioned into three distinct subsets: the testing set, which comprises 15% of the dataset (or 71 data), the validation set, which also comprises 15% of the dataset (or 71 data), and the training set, which comprises 70% of the dataset (or 330 data) [3,51].The data points are extracted from the bigger database via a standard diffusion and are then picked at random.An investigation using statistics was performed so that the practicability of the input parameters could be evaluated.Additionally, the input and outcome factor distribution diagrams are shown in Fig. 3.
The degree to which the data points in a scatterplot conform to a linear pattern is measured by the Pearson correlation coefficient [52], which provides a numerical value.The primary objective of the Pearson correlation coefficient is to ascertain the degree of association between two variables.The range of the coefficient is from − 1 to + 1.A positive numerical value indicates a positive linear relationship, whereby there is a propensity for both parameters to increase together.A negative number indicates a negative linear relationship, whereby there is a tendency for the other variable to decrease as one parameter increases.Values closer to + 1 or − 1 indicate a powerful linear relationship, while values closer to 0 suggest a weak or negligible linear correlation.The Pearson correlation coefficient is shown here in Fig. 4. The study may have been not capable of effectively capturing the impact of substantial negative or positive Pearson correlation coefficient variables on the findings due to a potential deficiency in powerful technique.It is noteworthy to observe that there are not many deviations from the high association values that exist among the parameters.As a consequence of this, the development of models that make use of these inputs needs to be successful for the purpose of reaching the highest possible accuracy.It is evident that a significant proportion of the association values exhibit high magnitudes.Among the investigated variables, the highest positive value is seen for In 8 and In 3 , both having a value of 0.99.The most negative value is seen in In 5 and In 2 , measuring -0.96.

Coot bird optimization (CBO)
The COOT method, developed by Iraj Naruei et al. (2021), is a recently introduced meta-heuristic optimization technique that aims to replicate the collective behavior of a swarm of coot birds.In comparison to other optimization methods, including the particle swarm optimization method and the differential evolution method, the COOT  method has rapid convergence velocity and great precision in convergence.Furthermore, the COOT method has undergone validation in several engineering applications, including pressure vessel design, welded beam design, stepped cone pulley, and rolling bearing difficulties [53].Hence, this study employs the COOT method as a solution to address the HDEED issue.
The underlying concept of the COOT method may be described below (Figs. 5 and 6).
(1) The first placement of the coot population will be established, followed by the computation of the associated fitness levels for each of these populations.
The variable Q is formed by the process of randomly initializing people.The variable d reflects the dimension of the problem that has to be solved.Additionally, the variables Ub and Lb represent the highest and lowest bounds of the population positions.The matrix displaying the positions of the coot is shown below.
Here, Cpos l shows the entire coot position, and N shows the population's scale.The inclusion of the coot's position is included in the fitness function, resulting in the generation of a fitness matrix for the coot population.
Here, fit l i shows the fitness value of the i th individual at the time of the l th iteration.
(  Fig. 6 Readjust the leaders' locations relative to the proper location [53] A location is generated at random inside the spatial domain, allowing a common one to afterward update the location.The formula governing the accidental mobility of the individual is as outlined below: Here, the location of the i th the typical individual is denoted as Cpos(i) , where A repre- sents a decreasing integer ranging from 0 to 1, R2 indicates an accidental value between 0 and 1, and Q indicates an accidental location inside the search area.The following for- mula determines the values of A and Q in the following way: Here, l is the iterations' present number, and L is the iterations' highest number that may be performed, respectively.
A distance vector among two individuals updates the location throughout the typical population chain movement.The formula for the chain is outlined below: In the context of leader-based location modification, the typical individual is required to modify its location by means of the k th leader.
Where Lpos(k) is the k th leader's location, R1 represents an accidental integer within the range [0,1], and the variable k is selected in the manner outlined below: In this equation, m represents the total number of leaders, and i indicates the present member of the i th coot group.
(3) Position update of the leader population As part of the leader's location modification, it is not uncommon for the leader to abdicate the currently finest location in order to assist in the discovery of a location that is more advantageous.If a leader is looking for a finer location, the ideal location from the prior instant should be assigned to the leader; alternatively, the leader should return to the place they were in before.The following outlines the leader's location: The variables R3 and R4 denote accidental values within the range of 0 to 1, whereas gbest indicates the present ideal position.The modified position of variable B is outlined below: (4)

Multi-layered perceptron (MLP) neural network
According to [54], ANN models are strong non-linear modeling techniques that operate by mimicking the functioning of the human brain.They are computer models for information processing.When the input and outcome parameters are mapped into a complicated network, an ANN may assess nonlinear functions.This research utilized a popular ANN with one input level, one outcome level, and at minimum one concealed level that was first put out by [55].According to [56], the ANN might be represented in the following way: Here, x i stands for the i th nodal value, N for the nodes' number, in the prior level; f shows the j th nodal value, y i shows the j th node's bias, b i denotes the activation function, and w ji stands for a weight-linking x i and y i .
Numerous investigations have shown that the ANN model can predict complicated nonlinear functions for hydrological data with only one concealed level [57].According to our research's preliminary findings, the link between the groundwater level and the other hydrologic cycle elements may be roughly approximated by one concealed level.Although optimal practices for finding the volume of the concealed nodes have also been offered, calculating the volume of the concealed nodes is commonly carried out by examination and mistake and is a crucial component of the ANN .According to Huang's (2003) research, the following amount of concealed nodes is required to learn N i samples with a little tiny fault: Here, N H shows the concealed nodes' highest volume, and N i and N o show the input nodes and outcome nodes.In this work, the concealed nodes' size was changed from one to N H , and the ideal size was determined via a process of examination and mistake.For training, the Levenberg-Marquardt ( LM ) method, one of ANN s' most effective meth- ods, was employed [58].The back-propagation method is applied to learn the biases and weights, and the CBO algorithm determines the optimal neuron numbers in each hidden layer (Fig. 7).In this study, the maximum number of each hidden layer is assumed to be 25.According to the literature, two hidden layers are considered for MLP structure to predict the bearing capacity of the concrete pile.MLPs often have hyperparameters (e.g., learning rate, number of layers) that impact their performance.CBO was employed to perform hyperparameter tuning, searching through the hyperparameter space to find the best configuration for the MLP.MLPs may suffer from overfitting, especially with complex architectures.Regularization techniques were optimized using optimization algorithms to find the right balance between model complexity and generalization.

Adaptive neuro-fuzzy inference system (ANFIS)
The basic model in the present work uses ANFIS as a machine learning method.The Takagi- Sugeno fuzzy system approach, which underpins ANFIS , was first suggested by Jang et al. in 1993 [59].ANFIS is capable of using both fuzzy logic and neural networks in a single frame if one thinks of it as an integrated notion.The ANFIS method's proper structure is chosen based on the input data, membership level, and input and result membership functions.By (11) modifying the level of membership variables in accordance with the allowed rate of error throughout the training phase, input values are able to be closer to observed values.The ANFIS technique makes use of fuzzy logic and neural network learning methods to handle the nonlinear connection between the values of the input and the output.It has great potential for use in time series modeling and categorization problems [60].Based on the benefits of fuzzy rules, numeric data may be collected from a technically written rule, and the ANFIS is able to easily assess the intricate conversion of mankind's intelligence into fuzzy systems.The fuzzy if-then principles that may be learned to estimate nonlinear functions are used in its inference mechanism.ANFIS is hence regarded as a well-liked estimating tool in engi- neering domains.ANFIS is hence regarded as a well-liked estimating tool in engineering domains.In this work, fuzzy grouping according to the k-means grouping method was done using the fuzzy cluster means ( FCM ) technique, and the optimal optimization technique for the traditional ANFIS was found using the backpropagation technique [61].In the present work, the ANFIS model was implemented using fuzzy c-means ( FCM ) clustering by collect- ing the collection of rules required for producing the fuzzy inference system ( FIS ).Regarding ANFIS designs, a fundamental ANFIS design was created using the initial parameters.The CBO methodology was then used to optimize the developed ANFIS framework.The param- eters of the membership function for the proposed system were enhanced by the use of the CBO method in this research.The RMSE index was calculated as a fitness indicator in order to assess the precision of the optimization framework.Finally, a hybridized and optimized ANFIS CBO network was defined, whereby the parameters of the network, such as the num- ber of fuzzy words and the maximum number of iterations, were determined.

Indices to effectiveness evaluation
The performance comparison criteria were developed to fulfill the need for a standardized and quantifiable approach to assessing and comparing the overall efficiency of a Fig. 7 The schematic of CBO − MLP number of different models.The use of these metrics enables the stakeholders to perform an evaluation of the performance of a number of answers, determine regions in need of development, and make choices based on accurate information.In this study, the following metrics were calculated and included in the analysis: coefficient of determination ( R 2 ), root mean square error ( RMSE ), normalized root-mean-square ( NRMSE ), relative absolute error ( RAE ), root relative squared error ( RRSE ), and performance index ( PI ).Lower values of RMSE , RAE , RRSE , and PI are indicative of greater accuracy.Fur- thermore, a greater R 2 value signifies enhanced performance.
In these equations, m d , m , z d , and z are the observations, the average of the observa- tions, the simulations, and the average of the simulations, respectively.Also, D depicts the total number of datasets.

Model development
The determination of the bearing capacity of the concrete piles was conducted utilizing the hybridized CBO models as outlined in this study.The hybrid models denoted as MLP CBO and ANFIS CBO , were used to ascertain the maximum iterations and the count of fuzzy words via the employment of MLP and ANFIS techniques.The bearing capacity of the concrete piles was assessed and computed throughout the training, validation, and testing phases of the developed MLP CBO and ANFIS CBO techniques, as seen in Fig. 8.The error distribution is shown on the right side of Fig. 8.The metrics R 2 , RMSE , NRMSE , RAE , RRSE, and PI were used to conduct an analysis of how ( 13) effective MLP CBO and ANFIS CBO were in the research of procedure prediction (see Table 2 for more information).Furthermore, the reliability and efficacy of the developed models, namely XGB , WOA − XGB , and GA − DLNN , were assessed based on the results obtained from the present study, which primarily centered on the performance of these models.The findings were compared to those that were found in the previous research [3] and [51].
Based on the data, it seems that both the MLP CBO and ANFIS CBO models have sub- stantial potential for accurately predicting the pile-bearing capacity.During the training stage, the R 2 values for MLP CBO were determined to be 0.954, while during the validat- ing stage, they were 0.9667, and during the testing stage they were 0.9729.in a similar manner, the R 2 values for ANFIS CBO during the training, stage were 0.9874, while dur- ing the validating stage, they were 0.9785, and during the testing stage they were 0.987.Evaluating the dependability of a model based just on one metric is inadequate.In order to achieve this objective, it is important to conduct a comprehensive examination of a range of measures, including but not limited to RSME , NRMSE , RAE , RRSE , and PI .When one compares the percentages that are reported for MLP CBO and ANFIS CBO for a variety of metrics, it is easy to see that ANFIS CBO reports a significant drop in com- parison to MLP CBO .The observed drop in performance demonstrates the high level of accuracy shown by the ANFIS CBO in estimating the bearing capacity of concrete piles.
In order to assess the reliability of the models, a comprehensive evaluation is conducted by comparing them with existing literature, specifically considering the models XGB [51,51] WOA − XGB , and GA − DLNN [3].A comprehensive examination of Table 2 reveals that the ANFIS CBO model, as introduced in this study, had favorable results in comparison to previous investigations discussed in the existing literature.The assessment was conducted using consistent metrics for the training, validation, and testing datasets, namely R 2 (coefficient of determination) and RMSE (root mean square error).The best model, ANFIS CBO , has higher R 2 values and lower RMSE compared to XGB [51,51]WOA − XGB , and GA − DLNN [3], indicating more precision and strength in its conclusions.For example, according to the results obtained from the WOA − XGB model, the coefficient of determination ( R 2 ) increases from 0.97 to 0.9874 in the train- ing stage and from 0.94 to 0.987 in the testing stage.Additionally, the metrics based on RMSE error shown a drop from 62.74 to 39.6769 in the training stage and from 87.72 to 38.6543 in the testing stage.A comprehensive comparison between ANFIS CBO and GA − DLNN [3] may be undertaken, taking into account the improvements seen in the training, validating, and testing data sets.These improvements include an increase in the values of R 2 and a decrease in the values of RMSE.
The distribution of the predicted/observed ratio of the algorithms throughout both the training and testing phases can be viewed on the right side of Fig. 8. Demonstrating enhanced efficacy is shown by a limited margin of error, including both a minimum and maximum threshold.The findings indicate that ANFIS CBO has higher efficacy in two phases, as shown by a decrease in error variability.

Recommendations for future studies
The focus of the current study was on the installation of driven reinforced concrete piles and their interaction with certain soil conditions.In the next research, it may be investigated how effectively the previously constructed machine learning models function with various soil profiles and pile kinds, such as drilled heaps and wood piles.Investigating ensemble models (such as random forests, gradient boosting, or model stacking techniques) might potentially enhance prediction accuracy even more.This is accomplished by combining a large number of machine-learning models.In subsequent research, field validation on actual construction projects may play a role, with the purpose of determining how useful the suggested models are when applied to actual-life circumstances.There is also the possibility of doing research on the challenges and considerations that go into effective implementation.It is feasible that by performing case studies on a variety of construction projects and contrasting the findings to the standards of the industry, it will be possible to learn a great deal about the effectiveness of machine learning models as well as the potential cost savings they may provide.

Sensitivity analysis
Sensitivity analysis is a technique used in various fields, including finance, engineering, environmental science, and decision-making processes [18].The primary aim of sensitivity analysis is to assess how changes in the input parameters or assumptions of a model affect the output or outcomes.It helps in understanding the robustness and reliability of a model or system under different scenarios.Sensitivity analysis helps in identifying which variables or parameters have a significant impact on the model's output.By understanding the sensitivity of the model to different inputs, decisionmakers can focus on key factors that drive the results.In the present, the procedure mentioned by Khatti, and Grover was chosen as a sensitivity analysis ( S a )'s method [62].The results of the S a was provided in Fig. 9.As it is clear from this plot, most of the input variables have high impact on the target higher than 0.9444, with some exceptions.The highest value belonged to In3 at 0.9836, while the lowest value was for In4 at 0.7216 of S a (Table 3 in Appendix).

Conclusions
The hybridized CBO systems that were given in this study were used order to evalu- ate the bearing capacity of the concrete piles.The hybrid models denoted as MLP CBO and ANFIS CBO , in which the MLP and ANFIS techniques were specified for the pur- pose of model augmentation.Furthermore, the reliability and robustness of the developed models, namely XGB , WOA − XGB , and GA − DLNN , were assessed based on the results obtained from the present study, which mostly centered on the produced models.The results obtained in this study were compared with the findings reported in previous research studies [3,51].
1. Based on the data, it seems that both the MLP CBO and the ANFIS CBO models have a great lot of potential for accurately forecasting the pile-bearing capacity.During the training stage, the R 2 values for MLP CBO were determined to be 0.954, while during Fig. 9 The sensitivity analysis of variables on target the validating stage, they were 0.9667, and during the testing stage they were 0.9729.In a similar manner, the R 2 values for ANFIS CBO during the training stage were 0.9874, while during the validating stage, they were 0.9785, and during the testing stage they were 0.987.

Figure 2
is a graphic representation of the pile structure, which includes the geometrical factors associated with it as well as the location of the soil layers.The following are categories of input parameters: In 1 ) pile's diameter (mm) , In 2 ) the first soil layer's thickness that piles embedded (m) , In 3 ) the second soil layer's thickness that pile embedded (m) , In 4 ) the third soil layer's thickness that piles embedded (m) , In 5 ) top elevation of pile (m) , In 6 ) natural ground elevation (m) , In 7 ) the extra segment pile top's elevation (m) , In 8 ) the pile tip's depth (m) , In 9 ) the mean SPT blow count along the pile shaft, and In 10 ) the mean SPT blow count at the pile tip[3]

Fig. 1 a
Fig. 1 a Experimental position.b Static load test experiment to record bearing capacity of pile[3]

Fig. 8
Fig. 8 The performance of the models.(right side: error, left side: correlation)

Table 1
Statistical description of the input and output variables

Table 2
The CBO-based models' performance and literature comparison

Table 3
2. The comparison of the percentages of MLP CBO and ANFIS CBO across several error- based measurements reveals a notable decrease in ANFIS CBO in comparison to MLP CBO .The observed drop in performance demonstrates the high level of accuracy shown by the adaptive neuro-fuzzy inference system with coot bird optimization ( ANFIS CBO ) in estimating the bearing capacity of concrete piles.3. The results obtained from the ANFIS CBO model provides higher values of R 2 and lower RMSE compared to the XGB [51, 51]WOA − XGB , and GA − DLNN [3] models, indicating that the ANFIS CBO model exhibits more effectiveness and accu- racy.The outcomes obtained from the train, validate, and test datasets shown significant enhancements by increasing the R 2 values and decreasing the RMSE values.These outcomes may be considered when conducting a comprehensive comparison between ANFIS CBO and GA − DLNN [3]. 4. The findings of the distribution of the predicted/observed ratio indicated that ANFIS CBO has higher efficacy in two phases, as shown by a decrease in error variability.The partial dataset