Next Article in Journal
Deep Learning Based Electric Pylon Detection in Remote Sensing Images
Next Article in Special Issue
Hybrid Spatial–Temporal Graph Convolutional Networks for On-Street Parking Availability Prediction
Previous Article in Journal
Deriving Stand Structural Complexity from Airborne Laser Scanning Data—What Does It Tell Us about a Forest?
Previous Article in Special Issue
U-Net-Id, an Instance Segmentation Model for Building Extraction from Satellite Images—Case Study in the Joanópolis City, Brazil
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Permeable Breakwaters Performance Modeling: A Comparative Study of Machine Learning Techniques

1
Department of Civil Engineering, College of Engineering, University of Tehran, Tehran 1417613131, Iran
2
Department of Civil and Environmental Engineering, College of Engineering, Shiraz University, Shiraz 7134851156, Iran
*
Author to whom correspondence should be addressed.
Remote Sens. 2020, 12(11), 1856; https://doi.org/10.3390/rs12111856
Submission received: 27 March 2020 / Revised: 28 April 2020 / Accepted: 4 May 2020 / Published: 8 June 2020
(This article belongs to the Special Issue Data Mining and Machine Learning in Urban Applications)

Abstract

:
The advantage of permeable breakwaters over more traditional types has attracted great interest in the behavior of these structures in the field of engineering. The main objective of this study is to apply 19 well-known machine learning regressors to derive the best model of innovative breakwater hydrodynamic behavior with reflection and transmission coefficients as the target parameters. A database of 360 laboratory tests on the low-scale breakwater is used to establish the model. The proposed models link the reflection and transmission coefficients to seven dimensionless parameters, including relative chamber width, relative rockfill height, relative chamber width in terms of wavelength, wave steepness, wave number multiplied by water depth, and relative wave height in terms of rockfill height. For the validation of the models, the cross-validation method was used for all models except the multilayer perceptron neural network (MLP) and genetic programming (GP) models. To validate the MLP and GP, the database is divided into three categories: training, validation, and testing. Furthermore, two explicit functional relationships are developed by utilizing the GP for each target. The exponential Gaussian process regression (GPR) model in predicting the reflection coefficient (R2 = 0.95, OBJ function = 0.0273), and similarly, the exponential GPR model in predicting the transmission coefficient (R2 = 0.98, OBJ function = 0.0267) showed the best performance and the highest correlation with the actual records and can further be used as a reference for engineers in practical work. Also, the sensitivity analysis of the proposed models determined that the relative height parameter of the rockfill material has the greatest contribution to the introduced breakwater behavior.

Graphical Abstract

1. Introduction

Originally, ports were located naturally in a safe place away from natural hazards. However, with the increasing need for port space, engineers are always looking to build artificial structures named breakwaters in a cost-efficient manner with the least impact on the environment to eliminate wave energy and protect ports [1]. To build traditional types of breakwaters, such as rubble mound and the gravity type, a very high volume of material is required, which leads to an increase in the width of these breakwaters [2]. This volume of material reduces the naval area while changing the environment, ecology system, and landscaping [3]. Also, due to the disruption of the lateral drift process as well as water circulation, the construction of these structures causes severe erosion, heavy loading on nearby beaches, and a deterioration in the quality of water [2]. In addition, these breakwaters cause problems such as scour due to the high reflection of wave energy [1].
To solve the problems mentioned above, permeable breakwaters were introduced [4]. Although these types of breakwaters have a lower degree of protection compared to traditional breakwaters, their main advantage is to reduce the volume of materials and the structural problems of previous breakwaters. Also, these breakwaters cause the reduction of environmental problems and further increase the quality of water in ports [5]. The hydrodynamic performance of these breakwaters is measured based on wave reflection and transmission, which are expressed as coefficients of reflection and wave transmission [4]. In the following study, we have attempted to introduce an innovative hydrodynamic permeable breakwater behavior as a solution to improve the performance of the generation of this breakwater through lab experience. The introduced permeable breakwater is a two-layer vertical permeable type filled with high porosity rockfill materials (Figure 1). These materials change the flow patterns, and if used in the right way, could play a significant role in energy absorption and reduction of the transmission coefficient, which leads to the improvement of the level of protection and reliability of this type of breakwater, in addition to the advantages of the common permeable breakwater. This new type of permeable breakwater has many applications and can be used as a main or secondary breakwater and loading or unloading jetty. Also, it could be applied in different seabed conditions such as coastal slopes as well as various types of soil. In this study, the front permeable panel has a constant value of 50% permeability, and the rear panel permeability is variable.
Machine learning methods have been widely used for the modeling of sophisticated engineering and sensing problems [6,7,8]. This study aims to model the hydrodynamic behavior of the breakwater with machine learning regressors (MLR). For comparison purposes, several widely used regression methods are employed, such as linear regression models, regression trees, support vector machine, Gaussian process regression models, ensembles of trees, and genetic programming. Based on the performance of the models, the best models will be selected, and knowledge will be extracted from optimal models.
Numerous studies have been conducted by the authors to investigate the hydraulic performance of these types of breakwaters, such as the total and semi-immersed, single, and double rows of vertical slotted pipes, and different kinds of wave screen breakwaters that consist of horizontal slots. In general, most of these studies have been conducted through hands-on experience in the lab. In addition, several analytical models have been used with the eigenfunction expansion of different wave theories [9,10,11,12,13]. These models alone are not capable of predicting the hydraulic performance of these breakwaters. In all models, the friction coefficient and added mass were obtained through the model calibration of the experimental results [9]. In 2013, Koraim et al. tabulated the specifications of the most important theoretical and laboratory studies in this field [4].
Due to the limitations of mathematical models, soft computing is widely used to solve a variety of classifications and prediction problems in science, medicine, and engineering [10]. In recent years, there has been increasing interest in using machine learning algorithms in marine engineering problems, especially those related to the modeling of breakwater behavior [11,12]. Elbisy (2015) used multiple additive regression trees (MART) and multi-layer perceptron neural networks (MLP) to quantify the regular wave runup on smooth slopes of perforated coastal structures constructed on sloping beaches. The results showed that the accuracy of the MART method is higher than the MLP methods [13]. In another study in 2016, Negm and Nassar presented the data for a new formula for wave reflection on smoothed, roughed, and perforated sloped seawalls using nonlinear regression. The results are consistent with previous studies; however, they remain limited in their application [14].
Yagci et al. (2005) modeled the damage rate of different breakwaters with three methods, such as artificial neural network, multilinear regression, and fuzzy logic model [15]. They argue that, unlike multilinear regression, which does not have satisfactory results, the other two methods would be able to interpolate data in the absence of abundant data. Recently, in another study, Pourzangbar et al. (2017), used two methods, namely support vector regression and model tree algorithm, to predict the maximum depth of scoring at breakwater toe from an average of 95 different experimental results. They recognized that the model obtained by the model tree algorithm was more accurate and showed better results than the experimental relationships [16]. Similarly, with the experimental data of Van Der Meer [17], Koç et al. (2017) studied the stability of rubble-mound breakwaters and proposed an explicit model by utilizing genetic programming. They showed GP models have better predictive performances than Van der Meer’s stability equations.
Previous studies have shown that each machine learning regressor has its strengths and weaknesses, and their performance depends on the characteristics of the subject under study [18]. Therefore, it is desirable to compare their performance in the modeling of introduced breakwater hydraulic behavior.
In this study, the hydrodynamic behavior of permeable breakwater under irregular waves has been studied. To fulfill this purpose, several machine learning regressors are employed to relate wave variables and breakwater parameters to the reflection and transmission properties of the structure by using 360 lab results.

2. Experimental Set-Up

The data used in this study are obtained from the experimental results on the physical model of the mentioned breakwater. The experiments are conducted in the laboratory of the Soil Conservation and Watershed Management Research Institute (SCWMRI) of the Ministry of Agriculture. Six wave probes are utilized to measure water elevation, which results in measuring reflection and transmission coefficients.

2.1. Physical Model Scale

To adjust the model scale, the equality of the Froude number that allows the correct reproduction of the gravitational and inertia forces of the flow was used, with the simulation scale of 1:25 which means, the length and time scales are 1:25 and 1:5 respectively. The range of model parameters and spectral wave variables have been listed in Table 1.

2.2. Test Facility

The flume has 32.5 m long, 5.5 m wide, and 1-m depth that was equipped with a DHI wavemaker system. A pedal wavemaker is located at one end of the flume, and the other end of the flume contains a passive wave absorber in the shape of a porous beach, consisting of a 1:12 slope rock ramp to absorb energy and prevent reflection. The piston-generating pedal allows the generation of irregular waves perpendicular to the structure with proper steepness and different periods. Irregular waves have been used to reach the sea conditions, and the impact wave spectrum in experiments is the modified JONSWAP wave spectrum. Figure 2 shows the details of the experimental setup.

2.3. Physical Model Details

The introduced permeable breakwater is a two-layer vertical permeable breakwater where the chamber between the two panels is filled with high porosity rockfill materials (Figure 1). These materials modify the flow pattern and, if used correctly, can play a significant role in reducing energy absorption and the transmission coefficient. The front panel permeability is constant and set at 50%, while the rear panel permeability is set at 30% and 50% respectively.

3. Methodology

As mentioned, the study of wave information necessary to understand their hydrodynamic behavior [4]. These fundamental elements transmission and reflection of these structures provides key can be obtained through different methods, such as physical modeling and on-site investigation. However, implementing these methods is often cost-intensive and, in most cases, impossible [19]. Another approach is the use of mathematical models, however in practice, due to the complexity and nonlinearity of the problems, it remains challenging to obtain mathematical models. On the other hand, if such relationships exist, they should be validated with real data. Soft computing has been introduced to avoid the problems mentioned. The methods of machine learning used in this study are described below.

3.1. Linear Regression Model

Linear regression models have predictors that are linear in the model parameters, easily interpretable, and fast to perform. These properties make linear regression models the most popular models for initial use. However, the very limited shape of these models indicates that they often have poor prediction accuracy. After the application of a linear regression model, it is better to apply more flexible models, and thereafter compare the results.

3.2. Regression Tree

The classification and regression tree (CART) introduced by Breiman [20] is a kind of decision-tree algorithm that is made by recursive binary partitioning. CART is a method of creating a set of rules that ends in one class or value. CART can be used for classification and regression problems. If the target value is a categorical type, CART generates a classification tree, and when the response variable is continuous, it makes a regression tree. The result of using CART is a hierarchical binary tree that is created by dividing the subsets of the data, using all the predictor variables to generate two branches per node. This tree begins with the entire dataset. In each division, each predictor is evaluated to obtain the best cut-off point based on the degree of improvement or impurity. In the regression tree, the least-squares deviation impurity measure is used to divide the rules and integrity of fit criteria. This measurement is expressed by R(t) as follows [21]:
R ( t ) = 1 N w ( t ) i t w i f i ( y i ȳ ( t ) ) 2
ȳ ( t ) = 1 N w ( t ) i t w i f i y i
N w ( t ) = i t w i f i
where Nw (t) is the weighted number of records in node t, Wi is the weighting field value for the record i (if any), fi is the repeat field value (if any), Уi is the target field value, and ȳ ( t ) is the mean value of the responding variable at node t.
The least-squares deviation (LDS) function for the division of s at node t is calculated as follows:
Q ( s , t ) R ( t ) R ( t L ) R ( t R )
where R(tR) is the sum of the squares of the right child node, R(tL) is the sum of the squares of the left child node, and the division of s is also chosen to maximize the value of Q(s,t).

3.3. Multilayer Perceptron Artificial Neural Network

McCulloch et al. established the neural network method in the early 1940s [22]. In general, a neural network is a predictive tool for constructing a mathematical model of an unknown system. One of the well-known or perhaps the most well-known of these models is the multilayer perceptron neural network, or MLP, which usually has a feed-forward architecture. The multilayer perceptron neural networks is often trained via the back-propagation algorithm.
The MLP network consists of an input and output layer and at least one hidden layer. Each of these layers has several neurons and contains the processing unit (s), and each unit is fully connected to the next layer units by a weighted connection (Wij) [23]. The Y output is obtained by transferring the sum of the previous outputs and mapping by an activation function. Figure 3 basically illustrates the multilayer perceptron neural network. In this figure, Xi represents the inputs, and b/B represents the bias between the different layers. For complex and nonlinear problems, the hyperbolic tangent function or the sigmoid (or log-sigmoid) function can be used.
A neural network is trained by a set of data called training data. During the training process, the weights of the network are optimized to reach the stopping criteria. The training process has two basic steps, whereby the first step is to give the initial values to the system (initialization), and the second step is optimization. In the initialization step, the initial values of weights and bias are allocated to neuron networks. Initial values of weights and bias can be obtained randomly or by using a global optimization method such as simulated annealing.

3.4. Support Vector Machine

The support vector machine, which was originally introduced by Boser et al. in 1992 [24], is inspired by the concept of statistical learning theory. To model the hydrodynamic behavior of the breakwater that is studied, the back-vector machine will be employed as a regression method by introducing an ε-insensitive loss of function. The following is a brief description of the use of a support vector machine in regression problems. For more information, refer to [25,26]. Consider a set of data for model training; {(x1,y1),…(xn,yn), x ϵ Rn, y ϵ r, where x is input, y is output, Rn is the n-dimensional vector space, and r is the one-dimensional vector space. The ε-insensitive loss of function can be expressed as follows:
L ε ( y ) = 0   For   | f ( x ) y | <   ε   or   L ε ( y ) = | f ( x ) y | ε
This expression defines a ε-domain such that if the predicted value is in the domain, the loss is zero, while if the predicted value is outside the domain, the loss will be equal to the absolute value of standard deviation subtract ε. The primary purpose of a vector machine is to find a function f (x) that gives the deviation of ε from the actual output value and is uniform.
The final Equation in the support vector machine can be written as follows [25]:
f ( x ) = i = 1 nsv ( α i α i * ) K ( x i , x j ) + b .  
In the above relationship, αi and αi * are Lagrangian coefficients, nsv is the number of supporting vectors, and K(xi,xj) is a kernel function. Some common kernels such as homogeneous and non-homogeneous polynomials, radial basis function, and Gaussian function. In the present study, Gaussian, linear, quadratic, and cubic kernel functions are used to investigate the model.

3.5. Gaussian Process Regression

Gaussian process regression (GPR) models are probabilistic models based on non-parametric kernels. Consider the training dataset {(xi, yi); i = 1, 2, …, n} such that xi ∈ ℝd and yi ∈ ℝ are derived from an unknown distribution. A GPR model with respect to the new input vector xnew and the training data supports the prediction of the variable response value of ynew. A linear regression model is as follows:
y = xTβ + ε
where ε ~ N ( 0 , σ 2 ) . The error variance σ2 and the coefficient β are calculated from the data. A GPR model describes the response by introducing latent variables from a Gaussian process, f(xi), i = 1,2, …, n, and the explicit basic functions h. The covariance function of the hidden variables introduces smoothness of the response, and the base functions transfer the inputs of X to a p-dimensional feature space. A Gaussian process is a set of random variables such that any finite number of them, have a Gaussian distribution. If {f(x),x∈ℝd} is a Gaussian process and n has observations x1, x2, …, xn, the distribution of random variables f (x1), f (x2), …, f (xn) is also Gaussian. A Gaussian process is defined by the mean function m(x) and the covariance function k (x, x′). That is, if {f(x),x∈ℝd} is a Gaussian process, then E(f(x)) = m(x) and Cov[f(x),f(x′)] = E[{f(x)−m(x)}{f(x′)−m(x′)}] = k(x,x′). Now consider the following model:
h(x)Tβ + F(x)
Therefore, f(x)~ Gaussian process (0, k(x,x′)), which means that f (x) is obtained from a Gaussian process with a mean of zero and the covariance function of c(x, x′). h (x) is a set of basic functions that transfers the original feature vector x from Rd to the new feature vector h (x) in Rd. β is a p-by-1 vector of the coefficients of the base functions. This model represents a GPR model. An example of a response y can be modeled as follows:
P ( y i | f ( x i ) , x i ) ~ N ( y i | h ( x i ) T β + f ( x i ) , σ 2 )
Hence, a GPR model is a probabilistic model. The hidden variable f (xi) is introduced for each xi observation, which makes the GPR model non-parametric [27].

3.6. Genetic Programming

GP was originally inspired by the evolution process [27] An individual program refers to a parse tree that consists of functions (nodes) and terminals (leaves). A function set may include mathematical functions (such as sin, cos, sqrt, exp), mathematical operators (+, , ÷, ×), Boolean operators, conditional functions, iterative functions, and any user-defined functions. Also, a terminal set, contains independent variables, dependent variables, random coefficients, and constant values.
Producing the final result, called solution in GP, was achieved by performing the following step-by-step procedure. As indicated in Figure 4, initial models (populations) are randomly created from the terminal set and function set. Each model (individual) is evaluated based on defined statistical criteria such as RMSE to select next-generation parents. Then, GP uses selection methodologies, such as ranking, to randomly pick up a certain number of the fittest models at several times to select the next-generation parents. Next-generation offspring are produced using GA operators such as reproduction, crossover, and mutation. This will continue until the termination condition means the number of generations or the desired accuracy.
According to the above mentioned, the following steps should be taken to create a model in GP [28]: (1) Specify the terminal set. (2) Choose the appropriate functional set. This step can be challenging and misleading as choosing the wrong functions can lead to models that are not physically sound. In this paper, the functional set is selected based on trial and error. (3) Determine a fitness criterion to evaluate the accuracy of models. This criterion determines which models survive to produce the next generation. In this study, root mean square error (RMSE) was used as a fitness criterion. (4) Specify the model controlling parameters. This parameter includes gene linking function, GP operators, and chromosomes. These parameters are capable of controlling the size and accuracy of models and, if not considered, can produce nested models that are hard to interpret and computationally expensive. (5) Criterion for termination conditions: This criterion is specified in either the number of generations or the specific accuracy achieved.

4. Modeling of Reflection and Transmission Coefficients

According to the above, reflection and transmission coefficients play a key role in determining the behavior of permeable-type breakwaters. Reflection and transmission coefficients have been expressed in many studies as a function of wave variables such as wavenumber and wave steepness and main breakwater characteristics such as chamber width and wall permeability. It is also shown in these studies that some parameters like panel thickness, the shape of vertical panel piles, and their roughness have a secondary influence or no specific influence on the hydrodynamic behavior of permeable-type breakwaters [2,29,30].
In this study, due to the use of rockfill materials in the absorber chamber, the height of rockfill materials has been considered as a parameter for the first time. The other most important parameters affecting the behavior of reflection and transmission coefficients are selected by reviewing the technical literature related to the permeable breakwater. As a result, the relationship between the reflection and transmission coefficients can be expressed as follows:
C t ,   C r = f ( B h , d h , B L p , H s L p ,   Kh , H s d , p )
where Ct and Cr, respectively, are transmission and reflection coefficients, B is the width of breakwater chamber (distance between front and back panel), h is the depth of water, d is the height of rockfill materials, Lp is wavelength associated with peak period (Tp) wave spectrum, Hs significant wave height of wave spectrum, K is wave number or angular frequency (is equal to 2 π L p ), and p is the permeability of the back wall, defined as the ratio of permeable wall openings to the entire vertical surface of the wall from the water surface to the bed of the flume. Previous studies on Jarlan’s multi-layer breakwaters without the use of rockfill materials have suggested that the presence of a back wall with less porosity compared to the front wall will improve the performance of the breakwater [31,32]. In some other studies on multi-layer permeable breakwater, it has been shown that if the back wall has more than 40% porosity, it is ineffective. Given the fact that the introduced breakwater behaves differently from the previous breakwaters due to the presence of rockfill material in its core, this study seeks to investigate the effect of back wall porosity on suggested breakwater [33,34].
The water depth is assumed to be constant and 40 cm in all tests. Table 1 gives definitions and ranges of input parameters that are summarized for the main model used in experiments.

4.1. Data and Preprocessing

The experimental data were used to develop the model. The database covers a wide range of parameters influencing the reflection and transmission coefficients. The database consisting of 360 records and predictors in this study are relative chamber width (B/h), relative rockfill height (d/h), relative chamber width in terms of wavelength (B/Lp), wave steepness (Hs/Lp), wave number multiplied by water depth (kh), and relative wave height in terms of rockfill height (Hs/d). The primary purpose of using dimensionless parameters is to apply the results of a smaller-scale physical model to real work. Selected parameters have been considered and applied previously in several studies related to permeable breakwater [3,4,9] and have specific concepts in coastal engineering. For example, investigating the effect of Hs/Lp and kh on Ct and Cr is essential to understand the hydrodynamic characteristics of the present breakwater for coastal and deepwater regions. Also, investigating the effect of B/Lp, d / h ,   B / Lp , Hs/d, and P on Ct and Cr is required to select the appropriate and optimized structure configuration.
The porosity of the rockfill material was initially considered as the parameter candidate. However, in later studies, it was found that there are two limitations to the nominal size of the rocks. First, the size of the rocks must be large enough to prevent them from coming out of the front and back walls. Second, the limitations related to the implementation and supply of material issues, which limits the size of the rocks. With these two limitations, the range of porosity of the rocks is practically low and, therefore, it is not considered as a parameter. The effect of the shape and porosity of the materials may be studied by using concrete block armor with different shapes and porosities and can be the subject of future studies.
Table 2 summarizes the definitions and ranges of input parameters of the original model used in the experiments and shows the domain of model parameters and the standard deviation.
A cross-validation statistical method was used to evaluate the performance of machine learning models. This method is commonly used in machine learning approaches to compare and select a model in a predictive modeling problem, due to its ease in understanding, implementation, and in estimating modeling skills. In general, this method is considered to be less biased over other methods. This approach involves randomly dividing the dataset into Kc groups or folds of approximately equal size. The first fold is treated as a validation set (testing), the modeling is fitted to the remaining (Kc – 1) folds, and this partition is applied until each fold is tested one time and will continue (Kc − 1) times to fit the model [35].
Poor selection of Kc may lead to a false idea of the model skill level, such as having a high variance score (which may vary greatly based on the data used to fit the model) or high bias (which may be due to overestimating the model’s skills). Kc is usually 5 or 10, but there is no formal rule. As Kc increases, the difference between training datasets and subsequent iterations subsets becomes smaller. As the difference decreases, the bias becomes smaller. In this study, considering the number of data, k is equal to 8 (28). As a result, from a total of 360 sets of lab results, 315 datasets in 7-fold were used for the training process, and 45 were used in each step for model testing.
In the case of multilayer perceptron neural network, due to its nature instead of using the cross-validation method, the data are divided into three sections: training, validation, and testing with ratios of 70%, 15%, and 15%, respectively, and after constructing the model, its features are obtained. The number of datasets for training, validation, and testing data are 252, 54, and 54, respectively.
Frank and Todeschini (1994) proposed a minimum ratio of data number to the number of input variables to accept the model, which is 3, and also recommended that this ratio be greater than 3 [36]. In the present study, this ratio for the training and test datasets is 315/3 = 105 and 45/3 = 15, respectively, both of which are much higher than the proposed ratio.
For the neural network model, this ratio for training, validation, and test data are 252/3 = 83, 54/3 = 18, and 54/3 = 18, all of which are much higher than the proposed ratio.

4.2. Performance Measurement

Correlation between model outputs and measured values that are expressed by the correlation coefficient is one of the important tools in evaluating the predictive performance of a machine learning model.
Whenever the correlation coefficient is closer to one, the model is more correlated to the actual results. According to a rational principle, Smith stated that, if a model had a correlation coefficient of more than 0.8 (R > 0.8), there was a strong correlation between the model output and the actual results [37].
Another indicator of model performance is the prediction error rate of the models, which are mean absolute error (MAE) and root mean square error. Finally, the optimal model is selected based on two objectives:
(1)
Provide the best correlation expressed by the correlation coefficient
(2)
Provide the least error represented by the RMSE index
For this purpose, the following objective function is used to satisfy both goals simultaneously. Then, choosing the optimal model is achieved by minimizing the following objective function.
ρ = RMSE h ¯ i 1 R + 1
where h ¯ i is the mean of the predicted values for the reflection and transmission coefficients. Due to the lack of a definitive evaluation criterion [38], in addition to the correlation coefficient, RMSE, and objective function, several statistical indicators, including the relative error (RE) or index of agreement, the scatter index (SI), and the Nash–Sutcliffe efficiency coefficient (E). have been used to evaluate the performance of the models. In general, when the values of RE and E are close to 1, the model has a better predictive performance. Conversely, the lower the value of SI, the better agreement between predicted and measured values [39,40].
The R, RMSE, MAE, RE, SI, and E functions can be found in the Appendix A.

5. Results and Comparison

5.1. Comparison of Models and Model Selection

In this study, 20 linear and nonlinear regression models for modeling hydrodynamic behavior of breakwater using the parameters listed in Table 2. The results for the RMSE, R, and ρ models are presented in Table 3 and Table 4. The results show that, for the Cr output for the whole exponential Gaussian process model dataset, the least error and highest correlation are obtained simultaneously, and the objective function value has the lowest value, 0.0273. It is also evident that linear, linear SVM, and robust linear models are less capable of predicting the pass-through coefficient, although other GPR models perform well near the exponential Gaussian process model.
In the case of the Ct output, similarly, the exponential Gaussian process model has the least error and the highest correlation, and its objective function value is 0.0267, which is the least value. Concerning the target parameter, Ct, linear, linear SVM, and robust linear models are also weaker than other models, and other GPR models include squared exponential, matern5/2, and rational Quadratic models demonstrate the best performance after the exponential GPR model.
Therefore, according to the introduced objective function and the numbers obtained in Table 3 and Table 4, the exponential Gaussian process model is generally the best model to predict the reflection and transmission coefficients. The prediction results of the top four models of each output are also visualized in Figure 5.
The introduced innovative breakwater has not been studied before, and therefore there are no empirical equations in predicting reflection and transmission coefficients as a reference for engineers in practical work. As indicated in Section 3, in contrary to most machine learning approaches, genetic programming represents an explicit functional relationship between input and output variables. Using this capability, the following formula has been developed respectively for reflection and transmission coefficients:
Cr = d 4 h + ( d h ) 0.2 + ( H s L p ) 2 × ( Kd ) + ( Tan h ( B h ) ( B L p ) + 1 ) 16.85   ×   Kd 1.3   ×   [ ( d h ) 2   ×   kd   ×   ( Arctg ( d h ) 1 )   ×   ( H s d + H s L p ) ] ( 1 15 ) 2.07
Ct = Arctg [ ( d h ) 3 ( 2.07   ×   kd   ×   Ln ( d h ) ) ( 1 3 ) ) ] 0.2 + [ ( 0.31   ×   B L p ) ( 3.06   ×   H s L p ) + ( B h ) 0.2 ( 2 × d h ) ] 9 ( B h ) 0.25   ×   ( H s L p ) 0.2   ×   ( d h )   ×   Tanh ( L p 11.41   ×   H s )
Although GP models are not among the best-performing models, they still demonstrate a robust performance, according to Smith [38] (see Table 3 and Table 4), and can be a guide for engineering purposes for both targets. By looking at Equations (10) and (11), it could be found that the parameter P is not present in the Equations, and it probably implies that the reflection and transmission coefficients are not sensitive to this parameter. In the sensitivity analysis section, the relationship between model parameters and outputs is examined in detail.

5.2. Parameter Settings of Selected Models

To obtain effective parameter settings, several performances with different parameters have been performed, and the best parameters have been selected based on the trial and error approach and the values recommended in previous studies. The parameters of the exponential Gaussian process for the modeling of both targets are shown in Table 5.

5.3. Selected Models

After evaluating the models through cross-validation, the exponential Gaussian process models are proposed to predict the behavior of reflection and transmission coefficients. A comparison of the predicted values of reflection and transmission coefficients with the experimental data for the preferred models is shown in Figure 6. It should be noted that the results shown in Figure 6 correspond to one of the models extracted from the cross-validation method.
The error distribution of the predicted and actual values is shown in Figure 7. As can be seen, the noise is almost Gaussian (a normal distribution with mean zero), has a maximum of about zero, and has not deviated in any particular way. Also, the median (cumulative distribution of 0.5) is in the range of zero.

5.4. Sensitivity Study

Sensitivity analysis studies the contribution of input parameters in predicting outputs. A simple process was used to perform sensitivity analysis. The sensitivity percentage of output to each input parameter is obtained using the following Equations:
N i = f max ( x i ) f min ( x i )
S i = N i j = 1 i N j   ×   100
In the above Equations, f max ( x i ) and f min ( x i ) are the highest and the lowest predicted output values on the i-th input domain, respectively, so that the other parameters are equal to their mean values. The results of the sensitivity analysis of the model presented are shown in Figure 8.
For both reflection and transmission coefficients, the highest contribution is related to the d/h parameter, which is close to 40. From the outset, due to the importance of rockfill materials in eliminating wave energy and increasing wave resistance, this parameter was expected to have the greatest influence on the hydrodynamic behavior of the breakwater. Also according to previous studies, the probability of wave resonance at specified ratios (B/Lp = 0.5 + 0.5n) is another effective parameter in the reflection coefficient which is reflected in the B/Lp parameter.
However, the relative chamber width parameter (B/h) has little contribution to the reflection coefficient parameter, which means that the reflection coefficient is not sensitive to the breakwater width. Conversely, in the transmission coefficient case, the B/h parameter is the second most influential parameter. Changing the values of this parameter means increasing/decreasing the interaction of the incident wave with the rockfill materials in the absorber chamber that significantly leads to altering the wave energy, and consequently, the transmission coefficient. Another influential parameter is the wave steepness coefficient that is mentioned in previous studies. The least contributing parameter in both target coefficients is the permeability of the back wall, which may be related to slight variations in this parameter due to practical limitations (30%–50%).

6. Conclusions

In this study, 20 regression modeling methods, including linear and nonlinear models, have been used to model the hydraulic behavior of the innovative breakwater. The models presented are based on a database of 360 laboratory tests. The proposed models relate the reflection and transmission coefficients to seven dimensionless parameters ( B h ,   d h ,   B L p ,   H s L p ,   Kh ,   H s d ,   p ). The results show that linear models are not capable of predicting reflection and transmission coefficients, but GPR models, especially exponential GPR, present highly correlated models. The values of the objective function for the coefficients of reflection and transmission are 0.0273 and 0.0267, respectively. Also, several statistical indices have been calculated for each model.
The error study also showed that the error distributions are approximately Gaussian, and the selected models are not biased. The sensitivity analysis of the models showed that the reflection coefficient model has the most sensitivity to the relative rockfill height (d/h) and the relative chamber width in terms of wavelength (B/Lp) parameters, and the back-wall porosity (P) parameter is the least important in this model. In the case of transmission coefficient modeling, the results show the sensitivity of the proposed model to the relative rockfill height (d/h) and the relative chamber width (B/h). Similarly, the model has the least sensitivity to the P parameter. Therefore, it has been concluded the wave hydrodynamic parameters (reflection and transmission coefficients) depend strongly on the relative rockfill height. In general, the reflection coefficient increases with the relative rockfill height (d/h), and conversely, the transmission coefficient decreases with an increasing relative rockfill height (d/h).
Providing the optimal model and examining the parameters influencing the model constitute an important tool for engineers in future studies. Furthermore, two explicit functional relationships are developed by utilizing the GP that can be used for practical purposes.

Author Contributions

Conceptualization, M.G. and M.D.P.; methodology, M.G., M.R.N.; software, M.G.; validation, M.D.P., I.V. and M.R.N.; formal analysis, M.G.; investigation, M.G.; resources, M.G.; data curation, M.G., I.V.; writing—original draft preparation, M.G.; writing—review and editing, M.D, I.V.; visualization, M.G.; supervision, M.D.; project administration, M.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

I wish to thank my sister, Zeinab Gandomi, for her valuable and constructive suggestions and encouragement throughout my study. I would also like to extend my thanks to Karami, Hosein poor, and other lab members for their support while I used equipment in Soil Conservation and Watershed Management Research Institute lab during my experiments.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Error and Statistical Indicators Functions

MAE = i = 1 n | h i t i | n
RMSE = i = 1 n ( h i t i ) 2 n
R = i = 1 n ( h i h ¯ ) ( t i t ¯ ) i = 1 n ( h i - h ¯ ) 2 i = 1 n ( t i - t ¯ ) 2
RE = 1 - ( i = 1 N ( h i - t i ) 2 / i = 1 N ( | t i - h ¯ | + | h i - h ¯ | ) 2 )
SI = 1 N i = 1 N ( h i - t i ) 2 / h ¯
E = 1 - ( i = 1 N ( t i - h i ) 2 / i = 1 N ( h i - h ¯ ) 2 )
where hi and ti are, respectively, measured and predicted Cr or Ct values for the ith output,   h   ¯ and   t   ¯ are the average of the measured and predicted outputs, and n is the total number of tests.

References

  1. Nikoo, M.R.; Varjavand, I.; Kerachian, R.; Pirooz, M.D.; Karimi, A. Multi-objective optimumA design of double-layer perforated-wall breakwaters: Application of NSGA-II and bargaining models. Appl. Ocean Res. 2014, 47, 47–52. [Google Scholar] [CrossRef]
  2. Somervell, L.T.; Thampi, S.G.; Shashikala, A.P. Hydrodynamic characteristics of vertical cellular breakwater. J. Waterw. Port. Coast. Ocean Eng. 2017, 143, 04017015. [Google Scholar] [CrossRef]
  3. Shih, R.-S.; Weng, W.-K.; Chou, C.-R. The performance characteristics of inclined highly pervious pipe breakwaters. Ocean Eng. 2015, 100, 54–66. [Google Scholar] [CrossRef]
  4. Koraim, A.S.; Iskander, M.M.; Elsayed, W.R. Hydrodynamic performance of double rows of piles suspending horizontal c shaped bars. Coast. Eng. 2014, 84, 81–96. [Google Scholar] [CrossRef]
  5. Rageh, O.S.; Koraim, A.S. Hydraulic performance of vertical walls with horizontal slots used as breakwater. Coast. Eng. 2010, 57, 745–756. [Google Scholar] [CrossRef]
  6. Lary, D.J.; Alavi, A.H.; Gandomi, A.H.; Walker, A.L. Machine learning in geosciences and remote sensing. Geosci. Front. 2016, 7, 3–10. [Google Scholar] [CrossRef] [Green Version]
  7. Peres, D.J.; Iuppa, C.; Cavallaro, L.; Cancelliere, A.; Foti, E. Significant wave height record extension by neural networks and reanalysis wind data. Ocean Model. 2015, 94, 128–140. [Google Scholar] [CrossRef]
  8. Browne, M.; Castelle, B.; Strauss, D.; Tomlinson, R.; Blumenstein, M.; Lane, C. Near-shore swell estimation from a global wind-wave model: Spectral process, linear, and artificial neural network models. Coast. Eng. 2007, 54, 445–460. [Google Scholar] [CrossRef] [Green Version]
  9. Somervell, L.T.; Thampi, S.G.; Shashikala, A.P. Estimation of friction coefficient for double walled permeable vertical breakwater. Ocean Eng. 2018, 156, 25–37. [Google Scholar] [CrossRef]
  10. Yaghouby, F.; Ayatollahi, A. An arrhythmia classification method based on selected features of heart rate variability signal and support vector machine-based classifier. In Proceedings of the World Congress on Medical Physics and Biomedical Engineering, Munich, Germany, 7–12 September 2009; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  11. Kundapura, S.; Arkal, V.H.; Pinho, J.L. Below the data range prediction of soft computing wave reflection of semicircular breakwater. J. Mar. Sci. Appl. 2019, 18, 167–175. [Google Scholar] [CrossRef]
  12. Kuntoji, G.; Manu, R.; Subba, R. Prediction of wave transmission over submerged reef of tandem breakwater using PSO-SVM and PSO-ANN techniques. Ish J. Hydraul. Eng. 2018, 1–8. [Google Scholar] [CrossRef]
  13. Elbisy, M.S. Estimation of regular wave run-up on slopes of perforated coastal structures constructed on sloping beaches. Ocean Eng. 2015, 109, 60–71. [Google Scholar] [CrossRef]
  14. Negm, A.; Karim, N. Determination of wave reflection formulae for vertical and sloped seawalls via experimental modelling. Procedia Eng. 2016, 154, 919–927. [Google Scholar] [CrossRef] [Green Version]
  15. Yagci, O.; Mercan, D.E.; Cigizoglu, H.K.; Kabdasli, M.S. Artificial intelligence methods in breakwater damage ratio estimation. Ocean Eng. 2005, 18, 2088–2106. [Google Scholar] [CrossRef]
  16. Pourzangbar, A.; Brocchini, M.; Saber, A.; Mahjoobi, J.; Mirzaaghasi, M.; Barzegar, M. Prediction of scour depth at breakwaters due to non-breaking waves using machine learning approaches. Appl. Ocean Res. 2017, 63, 120–128. [Google Scholar] [CrossRef]
  17. Van der Meer, J.W. Rock Slopes and Gravel Beaches under Wave Attack. Ph.D. Thesis, Delft University of Technology, Delft, The Netherlands, 1988. [Google Scholar]
  18. Chen, W.; Xie, X.; Wang, J.; Pradhan, B.; Hong, H.; Bui, D.T.; Duan, Z.; Ma, J. A comparative study of logistic model tree, random forest, and classification and regression tree models for spatial prediction of landslide susceptibility. Catena 2017, 151, 147–160. [Google Scholar] [CrossRef] [Green Version]
  19. Alavi, A.H.; Gandomi, A.H. Prediction of principal ground-motion parameters using a hybrid method coupling artificial neural networks and simulated annealing. Comput. Struct. 2011, 24, 2176–2194. [Google Scholar] [CrossRef]
  20. Breiman, L.; Friedman, J.; Stone, C.J.; Olshen, R.A. Classification and regression trees, wadsworth statistics. In Probability Series, Belmont, California: Wadsworth; CRC: Boca Raton, FL, USA, 1984. [Google Scholar]
  21. Samadi, M.; Jabbari, E.; Azamathulla, H.M. Assessment of M5′ model tree and classification and regression trees for prediction of scour depth below free overfall spillways. Neural Comput. Appl. 2014, 24, 357–366. [Google Scholar] [CrossRef]
  22. Perlovsky, L.I. Neural Networks and Intellect: Using Model.-Based Concepts; Oxford University Press: New York, NY, USA, 2001. [Google Scholar]
  23. Hossein Alavi, A.; Hossein Gandomi, A.; Mollahassani, A.; Akbar Heshmati, A.; Rashed, A. Modeling of maximum dry density and optimum moisture content of stabilized soil using artificial neural networks. J. Plant Nutr. Soil Sci. 2010, 173, 368–379. [Google Scholar] [CrossRef]
  24. Boser, B.E.; Guyon, I.M.; Vapnik, V.N. A training algorithm for optimal margin classifiers. In Proceedings of the Fifth Annual Workshop on Computational Learning Theory, Pittsburgh, PA, USA, 27–29 July 1992; ACM: New York, NY, USA, 1992. [Google Scholar]
  25. Smola, A.J.; Schölkopf, B. A tutorial on support vector regression. Stat. Comput. 2004, 14, 199–222. [Google Scholar] [CrossRef] [Green Version]
  26. Williams, C.K.; Carl, E.R. Gaussian Processes for Machine Learning; MIT Press: Cambridge, MA, USA, 2006; Volume 2, p. 4. [Google Scholar]
  27. Koza, J.R. Genetic Programming: On the Programming of Computers by Means of Natural Selection; MIT Press: Cambridge, MA, USA, 1992; Volume 1. [Google Scholar]
  28. Pourzangbar, A.; Losada, M.A.; Saber, A.; Ahari, L.R.; Larroudé, P.; Vaezi, M.; Brocchini, M. Prediction of non-breaking wave induced scour depth at the trunk section of breakwaters using genetic programming and artificial neural networks. Coast. Eng. 2017, 121, 107–118. [Google Scholar] [CrossRef]
  29. Ahmed, H.; Schlenkhoff, A.; Oertel, M. Stokes second-order wave interaction with vertical slotted wall breakwater. Coast. Struct. 2013, 2, 691–703. [Google Scholar]
  30. Elbisy, M.S.; Mlybari, E.M.; Helal, M.M. Hydrodynamic performance of multiple-row slotted breakwaters. J. Mar. Sci. Appl. 2016, 15, 123–135. [Google Scholar] [CrossRef]
  31. Kondo, H. Analysis of breakwaters having two porous walls. Proc. Coast. Struct. 1979, 2, 962–977. [Google Scholar]
  32. Kono, T.; Tsukayama, S. Wave deformation on a barrier reef. Jpn. Soc. Civ. Eng. 1981, 1981, 27–38. [Google Scholar] [CrossRef] [Green Version]
  33. Huang, Z. Wave interaction with one or two rows of closely spaced rectangular cylinders. Ocean Eng. 2007, 12, 1584–1591. [Google Scholar] [CrossRef]
  34. Hagiwara, K. Analysis of upright structure for wave dissipation using integral equation. Coast. Eng. 1985, 1984, 2810–2826. [Google Scholar]
  35. Hastie, T.; James, G.; Witten, D.; Tibshirani, R. An Introduction to Statistical Learning; Springer: New York, NY, USA, 2013; Volume 112. [Google Scholar]
  36. Frank, I.E.; Todeschini, R. The Data Analysis Handbook; Elsevier: Amsterdam, The Netherlands, 1994; Volume 14. [Google Scholar]
  37. Smith, G.N. Probability and statistics in civil engineering. In Collins Professional and Technical Books; Smith, G.N., Ed.; Nichols Publishing Company: New York, NY, USA, 1986; Volume 244. [Google Scholar]
  38. Koç, M.L.; Özdemir, Ü.; İmren, D. Prediction of the pH and the temperature-dependent swelling behavior of Ca2+-alginate hydrogels by artificial neural networks. Chem. Eng. Sci. 2008, 63, 2913–2919. [Google Scholar] [CrossRef]
  39. Koç, M.L.; Balas, C.E.; Koç, D.İ. Stability assessment of rubble-mound breakwaters using genetic programming. Ocean. Eng. 2016, 111, 8–12. [Google Scholar] [CrossRef]
  40. Coulibaly, P.; Baldwin, C.K. Nonstationary hydrological time series forecasting using nonlinear dynamic methods. J. Hydrol. 2005, 4, 164–174. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of the introduced permeable breakwater.
Figure 1. Schematic diagram of the introduced permeable breakwater.
Remotesensing 12 01856 g001
Figure 2. Details of wave flume, model position, and wave probes locations; (a) Passive wave absorber; (b) Model setup; (c) Wave probes; (d) Wave generating system
Figure 2. Details of wave flume, model position, and wave probes locations; (a) Passive wave absorber; (b) Model setup; (c) Wave probes; (d) Wave generating system
Remotesensing 12 01856 g002
Figure 3. Basic representation of a Multilayer Perceptron Neural Network.
Figure 3. Basic representation of a Multilayer Perceptron Neural Network.
Remotesensing 12 01856 g003
Figure 4. Typical genetic programming flowchart [27].
Figure 4. Typical genetic programming flowchart [27].
Remotesensing 12 01856 g004
Figure 5. Measured versus predicted values for the top four models; (AD) reflection prediction model and, (E,F) transmission prediction model.
Figure 5. Measured versus predicted values for the top four models; (AD) reflection prediction model and, (E,F) transmission prediction model.
Remotesensing 12 01856 g005
Figure 6. Measured versus predicted (A) Cr using exponential GPR (B) Ct using exponential GPR.
Figure 6. Measured versus predicted (A) Cr using exponential GPR (B) Ct using exponential GPR.
Remotesensing 12 01856 g006
Figure 7. The error distribution of predicted and actual values for (A) Reflection Coefficient parameter (B) Transmission Coefficient parameter.
Figure 7. The error distribution of predicted and actual values for (A) Reflection Coefficient parameter (B) Transmission Coefficient parameter.
Remotesensing 12 01856 g007aRemotesensing 12 01856 g007b
Figure 8. Parameter contribution ratio for (A) Cr coefficient (B) Ct Coefficient.
Figure 8. Parameter contribution ratio for (A) Cr coefficient (B) Ct Coefficient.
Remotesensing 12 01856 g008
Table 1. Ranges of principal experimental parameters.
Table 1. Ranges of principal experimental parameters.
SymbolDefinitionRangeStandard Deviation.
LpSpectrum wave length (cm)0.85–5.151.16
HsSpectrum significant wave (cm)2.91–14.203.21
dRockfill Height (cm)30–569.63
BChamber width (cm)16–8022.63
Table 2. Definitions and ranges of scaled model parameters.
Table 2. Definitions and ranges of scaled model parameters.
ParameterDefinitionRangeAverageStandard Deviation
PBack wall porosity (%)30–504010
B/hRelative chamber width0.4–21.20.57
d/hRelative rockfill height 0.75–1.41.090.24
B/LpRelative chamber width in terms of wavelength0.03–0.860.190.15
Hs/dRelative wave height in terms of rockfill height 0.05–0.450.180.09
KhThe wavenumber multiplies the water depth0.49–2.961.010.57
Hs/LpWave steepness0.01–0.060.030.01
CrReflection coefficient0.16–0.540.370.09
CtTransmission coefficient0.11–0.770.440.1
Table 3. Proposed ML regressors and their result for the reflection coefficient.
Table 3. Proposed ML regressors and their result for the reflection coefficient.
ML Regressor Reflection Coefficient
ERelative Error (RE)Scatter Index (SI (%))RMSER SquaredOBJ Function (ρ)
Linear Regression0.83630.95369.70.0367940.830.0520
Interactions Linear0.93480.99706.10.0251970.920.0348
Robust Linear0.83600.98219.70.0368590.830.0521
Stepwise Linear0.92590.99576.50.0254460.920.0351
Fine Tree0.97510.99793.80.0259850.920.0358
Medium Tree0.94940.99675.40.0294670.890.0410
Coarse Tree0.88250.99098.20.0351100.850.0494
Linear SVM0.83550.98869.70.0367580.830.0520
Quadratic SVM0.94760.99605.50.0230150.930.0317
Cubic SVM0.96790.99774.30.0221480.940.0304
Medium Gaussian SVM0.97120.99794.10.0213940.940.0294
Coarse Gaussian SVM0.86470.98998.80.0339730.860.0476
Boosted Trees0.93580.99546.10.0267070.910.0369
Bagged Trees0.95820.99694.90.0249710.920.0344
Squared exponential GPR0.97300.99813.90.0203570.950.0279
Matern 5/2 GPR0.97970.99853.40.0202220.950.0277
Exponential GPR0.99990.99990.0140.0199490.950.0273
Rational Quadratic GPR0.97530.99823.80.0203660.950.0279
GP0.92070.99406.70.0251170.920.0346
ANN0.91000.99357.20.0267910.910.0370
Table 4. Proposed ML regressors and their result for the transmission Coefficient.
Table 4. Proposed ML regressors and their result for the transmission Coefficient.
ML RegressorTransmission Coefficient
ERelative Error (RE)Scatter Index (SI (%))RMSER SquaredOBJ Function (ρ)
Linear Regression0.79600.940415.100.0691520.780.0835
Interactions Linear0.94480.98567.850.0396210.930.0458
Robust Linear0.78260.935615.580.0710680.770.0860
Stepwise Linear0.94230.98498.020.0389900.930.0451
Fine Tree0.98100.99514.600.0380050.930.0440
Medium Tree0.94890.98677.550.0437230.910.0509
Coarse Tree0.87400.965411.860.0656570.810.0785
Linear SVM0.78660.935715.440.0705990.780.0852
Quadratic SVM0.95840.98936.820.0344230.950.0396
Cubic SVM0.98060.99504.650.0279940.960.0321
Medium Gaussian SVM0.98040.99484.670.0333010.950.0383
Coarse Gaussian SVM0.83510.948713.570.0644710.810.0771
Boosted Trees0.96080.98946.610.0349540.950.0402
Bagged Trees0.97050.99195.740.0375620.940.0433
Squared exponential GPR0.98970.99743.380.0246120.970.0282
Matern 5/2 GPR0.99450.99862.470.0245260.970.0281
Exponential GPR0.99990.99990.020.0234020.980.0267
Rational Quadratic GPR0.99310.99822.760.0247670.970.0284
GP0.91950.97849.480.0421690.940.0487
ANN0.93830.98388.30.0369000.940.0426
Table 5. Parameters setting of the exponential GPR in the modeling of reflection and transmission coefficients.
Table 5. Parameters setting of the exponential GPR in the modeling of reflection and transmission coefficients.
TargetParameterValue
Reflection coefficient
Kernel FunctionExponential
Kernel Scale1.664473
Basic FunctionConstant
Signal Standard Deviation0.063185
sigma0.063185
Transmission coefficient
Kernel FunctionExponential
Kernel Scale1.664473
Basic FunctionConstant
Signal Standard Deviation0.105257
sigma0.105257

Share and Cite

MDPI and ACS Style

Gandomi, M.; Dolatshahi Pirooz, M.; Varjavand, I.; Nikoo, M.R. Permeable Breakwaters Performance Modeling: A Comparative Study of Machine Learning Techniques. Remote Sens. 2020, 12, 1856. https://doi.org/10.3390/rs12111856

AMA Style

Gandomi M, Dolatshahi Pirooz M, Varjavand I, Nikoo MR. Permeable Breakwaters Performance Modeling: A Comparative Study of Machine Learning Techniques. Remote Sensing. 2020; 12(11):1856. https://doi.org/10.3390/rs12111856

Chicago/Turabian Style

Gandomi, Mostafa, Moharram Dolatshahi Pirooz, Iman Varjavand, and Mohammad Reza Nikoo. 2020. "Permeable Breakwaters Performance Modeling: A Comparative Study of Machine Learning Techniques" Remote Sensing 12, no. 11: 1856. https://doi.org/10.3390/rs12111856

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop