Modeling Compressive Strength of Eco-Friendly Volcanic Ash Mortar Using Artiﬁcial Neural Networking

: Forecasting the compressive strength of concrete is a complex task owing to the interactions among concrete ingredients. In addition, an important characteristic of the concrete failure surface is its six-fold symmetry. In this study, an artiﬁcial neural network (ANN) and adaptive neuro fuzzy interface system (ANFIS) were employed to model the compressive strength of natural volcanic ash mortar (VAM) by using the six-fold symmetry of concrete failure. The modeling was correlated with four parameters. To train and test the projected models, data for more than 150 samples were collected from the literature. Furthermore, mortar samples with varying proportions of volcanic ash were prepared in the laboratory and tested, and the results were used to validate the models. The performance of the developed models was assessed using numerous statistical measures. The results show that both the ANN and ANFIS models accurately predict the compressive strength of VAM with R-square above 0.9 and lower error statistics. The permutation feature analysis conﬁrmed that the age of specimens affects the strength of VAM the most, followed by the water-to-cement ratio, curing temperature, and percentage of volcanic ash.


Introduction
Construction industries are known for their hazardous effects on the environment [1,2]. Recent studies have shown that construction industries consume 40% of the world's energy [3,4]. Furthermore, they are responsible for 23% of air pollution, 45% of landfill waste, and 35% of CO 2 emissions [5,6]. Hence, efforts are being made worldwide to reduce the high carbon footprint of construction industries [7,8]. Concrete is an essential man-made material in the industry and the second-most utilized material in the world [9,10]. Cement, the binding material in concrete, causes 7% of the total CO 2 emissions and consumes around 15% of the total energy used by industries [11][12][13]. However, the demand for cement is at its peak and is estimated to reach 6000 million tons in the next four decades. To reduce the high carbon footprint and CO 2 emissions of construction industries, there is a need to reduce cement usage. The demand for cement can only be reduced if construction industries adopt alternative environmentally friendly binding materials [14,15].
Locally available natural materials have been successfully used as supplementary cementitious materials (SCMs) in several studies [16][17][18][19][20][21][22]. The main advantage of SCMs is that they reduce cement usage by as much as 30% [2,23]. Recent studies showed that replacing cement with SCMs enhances economical and sustainable construction [24][25][26]. CO 2 emissions per unit compressive strength of concrete are reduced when 20% cement is replaced with SCMs [27,28]. Besides the environmental benefits, replacing with SCMs increases the mechanical strength of concrete due to their pozzolanic nature [29][30][31]. In addition, it enhances the durability of concrete due to the improved pore structure [26,27,[32][33][34]. Hence, replacing cement with SCMs is a good choice to reduce the use of cement and enhance concrete strength and durability. However, using imported SCMs for concrete production in the Gulf area increases the overall cost of projects, which in turn affects the sustainability of the project [35,36]. Hence, in the last decade, engineers have focused on the potential and exploration of locally available SCMs, such as volcanic ash (VA), for use in construction projects. Recently, large reservoirs of volcanic pozzolan (containing natural basaltic VA) were found in several parts of western KSA [37,38]. These volcanic pozzolans can be used as SCMs in concrete without affecting the overall cost of the project as they are cheaper than imported SCMs.
Earlier research findings focused on the hydration of mortar and concrete in which cement is partially replaced with VA [39][40][41]. In addition, several experimental research studies have been performed on the effect of partial replacement of cement with VA on the strength of concrete [38,42]. However, to date, no model has been developed to accurately predict the compressive strength of concrete in which cement is partially replaced by VA. This is attributed to the complexity and composition of VA, as VA is a natural pozzolana and is found in various parts of the globe. The properties of VA vary with the origin, particle characteristics, mode of formation, and mineralogical composition [43,44]. The lack of accurate models for estimating the properties of VAM limits the practical use of VA. However, experimental evaluation of different mechanical and physical properties of VAM is tiresome, time-consuming, labor-intensive, and expensive. Moreover, the mechanical properties of mortars are affected by numerous factors including the skill of labor and accuracy of machines; thus, laboratory data are unreliable. In addition, an important characteristic of the concrete failure surface is its six-fold symmetry. This property can be applied in two steps. First by minimizing the number of experimental tests to one of the symmetric parts only and then evaluating the peers of each point to obtain the full surface. The second application of this property is to average the symmetric parts after fitting to lower the error and obtain a final surface consistent with the theoretical assumptions. Previously, solving the symmetry was time-consuming and required high-level mathematical skills; however, symmetrical problems can be solved easily with the help of artificial intelligence [45]. These reasons have increased the use of machine learning techniques in estimating the mechanical properties of mortars in recent years [46][47][48][49]. Different types of machine learning techniques have been successfully employed in different types of concrete and mortar [50][51][52][53]. Meanwhile, the application of machine learning in civil engineering is not restricted to the prediction of mechanical properties of mortar and concrete. Machine learning has been employed by practicing engineers and professionals to predict various important parameters [54][55][56][57][58]. Among different machine learning techniques, neural networks like ANN and ANFIS have gained more popularity due to their ability to predict responses in complex systems [59][60][61]. Both mentioned techniques are discussed in detail in a later section.
There is an indispensable requirement to develop accurate and fast methods for predicting the compressive strength of VAM and increase the use of such mortars. In this study, we investigated the effect of partially replacing cement in mortar with VA at different percentages as well as the effect of different curing temperatures at different ages of mortar. Several experimental tests were performed on mortar specimens to determine their maximum compressive strengths. Furthermore, we developed a computational model using a machine learning technique to accurately forecast the compressive strength of VAM. The results of this study will be vital in the development of green concrete to reduce global CO 2 emissions.

Experimental
Ordinary Portland cement (OPC), also known as Type-I cement, produced in a Saudi cement factory, was utilized to prepare mortar samples. Locally available VA, provided by the Super Burkani Blocks Company located in Jeddah, Saudi Arabia, was used. It was supplied in the form of a fine powder, with an average particle size of 30 µm. To obtain relatively very fine particles of VA, dry sieving was performed on the supplied VA with micro-sieves #400 and #635 to extract from it the fine VA (VF) and ultra-fine VA (VUF), respectively. The results of the physical and chemical analyses of the OPC and VA used in the current research are listed in Table 1.  [62]; ** LOI = loss on ignition.
The particle size analysis curves of the cement and VA are shown in Figure 1. Commercially available standard sand with a fineness modulus of 2.54 was used to prepare the mortar samples. Details of the sieve analysis of the sand are shown in Table 2.   Including the control mix (CM), a sum of four mixes were prepared by replacing cement with different percentages of VA (10%, 20%, and 30%). A consistent water-to-cement ratio (w/c) of 0.487 was exercised for all the mixes according to ASTM C192 [63]. The specific details of all prepared mixes are listed in Table 3. Including the control mix (CM), a sum of four mixes were prepared by replacing cement with different percentages of VA (10%, 20%, and 30%). A consistent water-tocement ratio (w/c) of 0.487 was exercised for all the mixes according to ASTM C192 [63]. The specific details of all prepared mixes are listed in Table 3. For all mixes, identical cubic specimens of size 50 mm were cast. After 24 h of casting, the specimens were demolded and moist-cured at different temperatures (20 • C, 40 • C, and 60 • C) until the testing age. The compressive strengths of the mortar specimens were evaluated at three ages (7, 28, and 91 days) according to ASTM C192. A universal compressive testing machine with a capacity of 300 kN and a loading rate of 1 mm/min was used to evaluate the compressive strength of the specimens. Three samples were tested for each age and the mean values were calculated. Further details of the experimental procedures are available in the literature [64].

Modeling
Recently, the use of ANN and ANFIS in predicting different engineering properties of concrete has improved exponentially in every field due to their numerous advantages [65,66]. Among all artificial intelligence techniques, many researchers in the field of civil engineering prefer ANN and ANFIS due to their outburst predictive capacity [67][68][69]. Models that can predict different properties of concrete based on experimental data can be developed using ANN and ANFIS. The detailed methodology adopted in this study for establishment of these models is presented herein.

Artificial Neural Network
ANNs exhibit high performance, require no detailed knowledge about the underlying mechanisms, and are more economical than other artificial intelligence methods. Previous studies employed various types of ANNs, including the spiking neural network, radial basis function network, Kohonen network, and Hopfield network. Each type of ANN has its own advantages and disadvantages for different types of data sets. The feed-forward neural network (FFNN) is widely used in predicting the compressive strength of concrete due to its simplicity and robustness.
FFNNs use unidirectional connections to create relationships between input parameters and target value [70,71]. The neurons present in various layers are also connected by unidirectional setups. Two types of FFNNs can be used to create predictive models. Single later perceptron is used for linear problems due to its simplicity, and multilayer perceptron is used for nonlinear complex problems due to its ability to efficiently solve complex problems. Here, we employed the multilayer perceptron. Multilayer perceptron consists of inputs, outputs, weights, neurons, mathematical functions, and hidden layers. Inputs are received from the user and transferred to the hidden layers. No computation is performed while transferring data from input to the hidden layers. All processing is performed in hidden layers using mathematical functions and neurons. After processing, the output is presented to the user for further analyses (if required). The whole process is summarized in Figure 2 for a couple of inputs and outputs, three multilayer perceptrons, and three neurons.       1. The summations of weighted inputs are calculated using weights and bias. 2. The output of the neurons is calculated using different mathematical functions, such as sigmoid and ramp. Here, the sigmoid function is employed due to its simplicity. 3. The final output is calculated using the following equations: 1.
The summations of weighted inputs are calculated using weights and bias.

2.
The output of the neurons is calculated using different mathematical functions, such as sigmoid and ramp. Here, the sigmoid function is employed due to its simplicity.

3.
The final output is calculated using the following equations: where w ik is the weight from the i − th node to the k − th target node and b k is the bias of the k − th target node.
Learning of ANNs starts after the finalization of biases and weights. Among all parameters, only the biases and weights affect the output values. Back-propagation is employed here to train the ANNs using multilayer perceptron, as recommended in previous studies.
The multi-layered FFNN method combined with back-propagation is used in this study. The nonlinear sigmoid function is used to connect two hidden layers. Four neurons are used as input layers, whereas a single neuron is obtained in the output layer. Two hidden layers are used to minimize the percentage error. Seven neurons are used in both hidden layers. Four parameters (percentage, cement, temperature, and day) are employed as inputs and the compressive strength is the output. Details of hyper-parameters utilized while performing an ANN modeling are provided in Table 4. Zadeh introduced the concept of fuzzy logic to represent complex systems [72]; it has gained popularity and also been effectively applied to a variety of civil engineering issues, particularly for forecasting mechanical properties of concrete composites made of waste ashes, like fly-ash, bagasse ash, and rice-husk-ash. Fuzzy logic has recently been strongly suggested for simulating properties of different types of concrete in order to address the intrinsic inaccuracy and ambiguous features [3]. The fundamental issue with fuzzy logic algorithms, however, is the unavailability of a consistent process for designing a fuzzy controlling algorithm. An ANFIS, on the other side, has the potential to understand from its surroundings (input to output combinations), modifying its architecture in a self-organized way, and adjusting it in an accessible condition. As a result, we suggest using the ANFIS technique to inherently re-organize neural network structure architecture and modify fuzzy system performance parameters for forecasting the compressive strength of VAM.
Unlike ANN, the ANFIS architecture combines ANN learning techniques and a fuzzy logical reasoning algorithm for a map showing an input space (representing an independent variable) to an output space (representing a dependent parameter). ANFIS has recently been used as an influential and powerful modeling technique for the estimation of numerous properties of concrete [3,73], due to its capacity to integrate the verbal and linguistic power of a fuzzy logic system with the numerical power (floating point data analysis) of a neural system. ANFIS is capable of learning, building, and categorizing. It has the benefit of extracting fuzzy logic rules from numerical information or professional expert information and constructing a rule framework flexibly. It can also fine-tune the complex transfer of human intellect to fuzzy networks. The major disadvantage of the ANFIS estimation method is the time required for structured methodology and hyper-parameter determination, which can take a long time. For the sake of clarity and simplicity, we can suppose that the fuzzy system in discussion contains two input variables (x 1 and x 2 ) and one outcome (z). A standard set of rules contains two IF-THEN rules for the interpretation of a first-order Takagi-Sugeno fuzzy model [74]. Equations (3) and (4) represent the two IF-THEN rules stated below.
The p i , q i , and r i ('i' is 1 or 2) dictate the linear or consequent parameters of the IF-THEN statements, also known as the consequent part. A concise explanation of the five layers as shown in Figure 4 is given below.

1.
First Layer (input nodes and fuzzification): This layer takes the inputs and recognizes the linked membership functions. Applying membership function parameters, all nodes in this layer produce membership scores indicating whether or not they correspond to each of the relevant fuzzy sets.

2.
Second Layer (rule nodes): In the rule layer, the firing strengths for the designated rules are generated. The AND operator is used to produce a single outcome that reflects the consequence of the predecessor (i.e., firing strength). Firing intensity refers to the degree to which the predecessor portion of a fuzzy set rule is fulfilled. It shapes and influence the rule's output function. interpretation of a first-order Takagi-Sugeno fuzzy model [74]. Equations (3) and (4) represent the two IF-THEN rules stated below.
The , , and ('i' is 1 or 2) dictate the linear or consequent parameters of the IF-THEN statements, also known as the consequent part. A concise explanation of the five layers as shown in Figure 4 is given below.  Like ANN, ANFIS was also executed in MATLAB. The same four inputs (percentage, cement, temperature, and day) and single outputs (compressive strength) were employed. After uploading the data arrays, the number of hyper-parameters provided in Table 5 were selected next to several initial trials. The Takagi-Sugeno fuzzy interface structure with triangular membership (trimf) in combination with a hybrid optimization procedure  Like ANN, ANFIS was also executed in MATLAB. The same four inputs (percentage, cement, temperature, and day) and single outputs (compressive strength) were employed. After uploading the data arrays, the number of hyper-parameters provided in Table 5 were selected next to several initial trials. The Takagi-Sugeno fuzzy interface structure with triangular membership (trimf) in combination with a hybrid optimization procedure (backpropagation and least-squares) gives the best performing model [74]. The ANFIS provides the leading performance in comparison with ANN because of its integrated algorithmic architecture of ANN and fuzzy logic rule set.  [32,[75][76][77][78][79][80][81][82][83] and analyzed to find the parameters affecting the properties of VAM the most. The literature review revealed that the properties of mortar with admixtures and SCMs are affected by the type and dosage of the admixtures [84]. Thus, we chose the proportion of VA (VA%), w/c, curing temperature, and age of the specimens as input parameters. The sand-to-cement ratio for the whole database was kept constant.

Data Representation Using Python
The database representation was performed using Python 3.7. The database consists of four input parameters named VA%, curing temp, age, and binder ratio, whereas compressive strength is the output. All the input parameters selected affect the compressive strength of the mortar directly. The complete dataset was divided into three sets; 70% was utilized to train the model, 15% to test the data, and 15% to validate the model. This distribution can be changed in different studies; however, it is recommended by numerous researchers. The maximum and minimum values for all parameters (input and output) and some statistical indicators are shown in Table 6. The correlation of input variables to the output variable, optimal dosage, and effect of all variables were calculated using built-in libraries of Python 3.7. The details are shown in Figure 5. It is well-known that the variables in a model have a significant impact on its performance. Deep learning is a method used in a neuron-based artificial intelligence approach to estimate material properties depending on the concentration of variables. Python was used to implement a machine learning technique, and the correlation plot was created with the help of the seaborn command. The performance evaluation and accurateness of the model are greatly affected by the selected input variables.

Python-Based Contour Maps
Using Python, five contour maps were drawn and are shown in Figure 6. As discussed in the previous section, the model performance is greatly reliant on the range of input variables. Thus, it is important to determine the optimal quantity of the variables rather than using experimental values. This offers a valuable representation to estimate the strength of VAM.

Performance Evaluation of Prediction Models
The accuracy of the developed model was evaluated using six performance indicators, including the mean absolute error (MAE), coefficient of determination (R-square), rootmean-square error (RMSE), root-squared-error (RSE), percentage of relative-root-meansquare-error (RRMSE%), and performance index (ρ). The greater the value of R-square and the lower the stated error statistics, the better the performance of the models will be. The performance index is the ultimate metric for performance evaluation, as it combines the effect of both error statistics and goodness of fit. It is recommended in the literature that the performance index and RRMSE% should be less than 0.2 and 10%, respectively, to classify the model as good [85,86]. The formulas for calculating the mentioned metrics of performance assessment are as follows.
where ex i and Po i are the actual measured and model predicted values of compressive strengths, respectively, and ex i and Po i are the mean values of the actual measured and model predicted compressive strengths, respectively.

Python-Based Contour Maps
Using Python, five contour maps were drawn and are shown in Figure 6. As discussed in the previous section, the model performance is greatly reliant on the range of input variables. Thus, it is important to determine the optimal quantity of the variables

Performance Evaluation of Prediction Models
The accuracy of the developed model was evaluated using six performance indicators, including the mean absolute error (MAE), coefficient of determination (R-square), root-mean-square error (RMSE), root-squared-error (RSE), percentage of relative-rootmean-square-error (RRMSE%), and performance index ( ). The greater the value of Rsquare and the lower the stated error statistics, the better the performance of the models will be. The performance index is the ultimate metric for performance evaluation, as it combines the effect of both error statistics and goodness of fit. It is recommended in the literature that the performance index and RRMSE% should be less than 0.2 and 10%, respectively, to classify the model as good [85,86]. The formulas for calculating the mentioned metrics of performance assessment are as follows.

Test Findings Compressive Strength in Respect of Curing Temperature and Age
As mentioned earlier, the main objective of this research was to explore the effect of VA replacement as cement up to 30% on the compressive strength development of mortar at different curing temperatures and aging. The compressive strengths of the various mortar samples at several aging and curing temperatures are listed in Table 7. In the specimens containing VA, cement was substituted up to 30% with increments of 10%.
Under standard curing conditions of 20 • C, the compressive strength of the specimens with VA was reduced as compared to CM at all ages, irrespective of the VA%. Moreover, the decrease in compressive strength increased as VA% increased up to 20%. Although the compressive strength decreased at early ages, most possibly because of the lesser or no pozzolanic activity, the mortar with 30% VA exhibited compressive strength comparable to that of CM at later ages (91 days). This could be credited to the delayed pozzolanic reaction of VA. These outcomes show that the partial substitution of cement with VA might not be beneficial, regarding the development of compressive strength. However, unlike standard curing conditions, the compressive strength of VA was significantly enhanced when subjected to moderate or high curing temperatures. This demonstrates an improved reactivity of the mortar with VA under high curing temperatures. For instance, the compressive strength of VAM samples with 10% and 20% VA (VF10 and VF20) was comparable to that of CM at a higher age (i.e., 91 days) and high curing-temperature (i.e., 60 • C). The satisfactory performance of VAM under higher curing temperature is a good indication and thus can be used under hot environmental conditions.  Figure 7 shows the regression analysis plot between the experimental and predicted training, testing, and validation datasets of the ANN predictive model. The linear leastsquares fit lines for each data set are also shown in the figure. The error between the experimental and predicted values is very small. R-square for the three datasets is greater than 0.5 and close to 1 i.e., 0.84 and 0.86 each, respectively, showing a close relationship between the predicted and experimental values. Table 8 lists the performance indicators for the three datasets. The MAE, RMSE, RRMSE, RSE, and ρ for the training set are 3.33, 4.45, 0.017, 0.16, and 0.1, respectively. The MAE is less than the MAE, showing the smaller percentage of larger error values [85]. In addition, the percentage of RRMSE is less than 10% along with ρ less than 0.2 [86]. Thus, the ANN model gives a good prediction of the compressive strength of VAM. The applicability of the developed ANN model was further evaluated by calculating the error distribution between the predicted and experimental values of training, testing, and validation datasets. Figure 8 shows the error distribution for the training, testing, and validation sets. The dark red shows the higher error values. The error value is low, which implies that the model accurately predicts the compressive strength of VAM. Table 9 presents the detailed statistical analysis and Figure 9 shows the regression analysis of training, testing, and validation set data of the ANFIS proposed model. The location of the points in the figures gives an idea to estimate the R-square. For the best model, the points will be nearer to the regression line. As shown in the figure, all the points in each stage are located nearer to the line. In all three stages, the R-square is greater than 0.9, i.e., 0.94, 0.95, and 0.96, for training, testing and validation sets. In addition, the error values are very low in comparison with the ANFIS predictive model. The MAE, RMSE, RRMSE, RSE, and ρ for the training set are 1.82, 2.51, 0.01, 0.058, and 0.05, respectively. Like ANN, for ANFIS the MAE values in each stage are lower than RMSE, depicting on the lower percentage of higher error values. Consequently, the RRMSEs for training, testing and validation sets are 0.01, 0.009, and 0.004, respectively, which are lower than 10% [85,86]. In addition, the performance index is nearer to 0, i.e., 0.05 and 0.2 for training/testing and validation sets, respectively. Comparatively, the ANFIS outperforms and gives an excellent prediction of the compressive strength of VAM. perature (i.e., 60 °C). The satisfactory performance of VAM under higher curing temperature is a good indication and thus can be used under hot environmental conditions. Figure 7 shows the regression analysis plot between the experimental and predicted training, testing, and validation datasets of the ANN predictive model. The linear leastsquares fit lines for each data set are also shown in the figure. The error between the experimental and predicted values is very small. R-square for the three datasets is greater than 0.5 and close to 1 i.e., 0.84 and 0.86 each, respectively, showing a close relationship between the predicted and experimental values. Table 8 lists the performance indicators for the three datasets. The MAE, RMSE, RRMSE, RSE, and ρ for the training set are 3.33, 4.45, 0.017, 0.16, and 0.1, respectively. The MAE is less than the MAE, showing the smaller percentage of larger error values [85]. In addition, the percentage of RRMSE is less than 10% along with ρ less than 0.2 [86]. Thus, the ANN model gives a good prediction of the compressive strength of VAM.  The applicability of the developed ANN model was further evaluated by calculating the error distribution between the predicted and experimental values of training, testing, and validation datasets. Figure 8 shows the error distribution for the training, testing, and validation sets. The dark red shows the higher error values. The error value is low, which implies that the model accurately predicts the compressive strength of VAM.  Table 9 presents the detailed statistical analysis and Figure 9 shows the regression analysis of training, testing, and validation set data of the ANFIS proposed model. The location of the points in the figures gives an idea to estimate the R-square. For the best model, the points will be nearer to the regression line. As shown in the figure, all the points  For further assessment of the performance of the ANFIS predictive model, the predicted and error values of every data point are presented on the heat map in Figure 10, which shows the range of errors between experimental values and ANFIS predicted results in the form of different colors. The lower error values are less in each data set, which shows the accurateness of the ANFIS predictive model.

Comparison between ANN and ANFIS Predictive Models
Each neural network (ANN and ANFIS) outperforms on the data provided for compressive strength results of VAM having four input variables. For comparing the performance of each predictive model, the targeted and predicted results along with the absolute error for every single data point are presented in Figures 11 and 12. It can be seen that the absolute error of both prediction techniques runs near the x-axis, showing the accurate prediction of compressive strength of VAM. However, the maximum absolute errors for ANN and ANFIS models are 13.94 MPa and 6.89 MPa, respectively. The maximum error in ANFIS is below 10 MPa. The ANFIS predictions nearly follow their targeted values, giving an outburst estimation of the compressive strength of VAM. Based on RRMSE and performance index, the ANFIS model gives 41% and 50%, respectively, providing a better prediction than the ANN predictive model. The leading performance of ANFIS is attributable to its integrated algorithmic structure of ANN and a fuzzy logic technique. For further assessment of the performance of the ANFIS predictive model, the predicted and error values of every data point are presented on the heat map in Figure 10, which shows the range of errors between experimental values and ANFIS predicted results in the form of different colors. The lower error values are less in each data set, which shows the accurateness of the ANFIS predictive model.

Comparison between ANN and ANFIS Predictive Models
Each neural network (ANN and ANFIS) outperforms on the data provided for compressive strength results of VAM having four input variables. For comparing the performance of each predictive model, the targeted and predicted results along with the absolute error for every single data point are presented in Figures 11 and 12

Permutation Feature Analysis
There are several applications of neural networks in the field of engineering and research. Engineers employ ANN and ANFIS in developing new rules for optimizing different mechanical properties of concrete. However, these are often provided to clients in the form of a black box containing different nonlinear internal algorithms to estimate outputs based on the provided inputs. Due to this complexity of neural-network-based models, it is impossible to check the mechanisms involved just on the basis of network bias and weights for every problem [74]. Hence, it is difficult to convert ANN models in the form of simple equations connecting input variables with output. In addition, the effect of

Permutation Feature Analysis
There are several applications of neural networks in the field of engineering and research. Engineers employ ANN and ANFIS in developing new rules for optimizing different mechanical properties of concrete. However, these are often provided to clients in the form of a black box containing different nonlinear internal algorithms to estimate outputs based on the provided inputs. Due to this complexity of neural-network-based models, it is impossible to check the mechanisms involved just on the basis of network bias and weights for every problem [74]. Hence, it is difficult to convert ANN models in the form of simple equations connecting input variables with output. In addition, the effect of a single input parameter on the output cannot be calculated using traditional statistical methods.
To further verify the performance of the developed models, permutation feature analysis (PFA) was conducted to check the most influential input variable affecting the compressive strength of VAM [1,3,5]. For performing PFA, the extension from the Python programming software was used. As shown in Figure 13, the effect of each input on the compressive strength of VAM is prominent. The most important factor affecting the compressive strength of VAM is the age of curing (43.59%), followed by water-to-cement ratio (22.79%), curing temperature (17.87%), and volcanic ash percentage (15.75%). analysis (PFA) was conducted to check the most influential input variable affecting the compressive strength of VAM [1,3,5]. For performing PFA, the extension from the Python programming software was used. As shown in Figure 13, the effect of each input on the compressive strength of VAM is prominent. The most important factor affecting the compressive strength of VAM is the age of curing (43.59%), followed by water-to-cement ratio (22.79%), curing temperature (17.87%), and volcanic ash percentage (15.75%).

Conclusions
This study had a two-fold objective. First, we modeled the compressive strength of volcanic ash mortar (VAM) using an artificial neural network (ANN) and an adaptive neuro-fuzzy interface system (ANFIS). The models were proposed and formulated on the basis of detail and a high range data set developed from the literature. Second, VA was characterized, and used in varying proportions (i.e., 10%, 20%, and 30%) partially replacing the cement in mortar. VAM specimens were prepared, and their mechanical properties were tested. The results show that the use of different natural pozzolans as a substitute for cement is indispensable. Regarding the sustainability, machine learning models perform a central part in its achievement. Based on the results, analysis, and discussion in the preceding sections, the subsequent conclusions are presented.
1. Compared to the standard curing temperature (20 °C), the compressive strength of VA-containing mortar improved significantly under high curing temperatures (40 °C and 60 °C). This is a good indication that such mortar could be applied under hot environmental conditions. 2. Among all mortar mixes, the mortar containing a high percentage (30%) of ultra-fine volcanic ash (VUF30) exhibited the highest compressive strength at both 40 °C and 60 °C. This indicates that VA can be used in high percentages as a partial substitute for cement.

Conclusions
This study had a two-fold objective. First, we modeled the compressive strength of volcanic ash mortar (VAM) using an artificial neural network (ANN) and an adaptive neuro-fuzzy interface system (ANFIS). The models were proposed and formulated on the basis of detail and a high range data set developed from the literature. Second, VA was characterized, and used in varying proportions (i.e., 10%, 20%, and 30%) partially replacing the cement in mortar. VAM specimens were prepared, and their mechanical properties were tested. The results show that the use of different natural pozzolans as a substitute for cement is indispensable. Regarding the sustainability, machine learning models perform a central part in its achievement. Based on the results, analysis, and discussion in the preceding sections, the subsequent conclusions are presented.

1.
Compared to the standard curing temperature (20 • C), the compressive strength of VA-containing mortar improved significantly under high curing temperatures (40 • C and 60 • C). This is a good indication that such mortar could be applied under hot environmental conditions. 2.
Among all mortar mixes, the mortar containing a high percentage (30%) of ultra-fine volcanic ash (VUF30) exhibited the highest compressive strength at both 40 • C and 60 • C. This indicates that VA can be used in high percentages as a partial substitute for cement. 3.
ANN and ANFIS algorithms were exercised to estimate the compressive strength of the VAM. The accuracy, reliability, and performance of the final proposed models were assessed using statistical indicators, including RMSE, RSE, RRMSE, MAE, ρ, and R-square. The predicted results displayed an exceptional correlation with the experimental data with R 2 value above 0.9 and RMSE and MAE below 5. Based on RRMSE and performance index, the ANFIS model gave 41% and 50%, respectively, providing a better prediction than the ANN predictive model.

4.
The permutation feature analysis (PFA) revealed that the developed model considers the effect of all the inputs. The most influential factor affecting the strength of VAM is the age of the specimen, followed by w/c, curing temperature, and VA%.