Acessibilidade / Reportar erro

Revisiting the derivation of bulk longshore sediment transport rates using meta-heuristic algorithms

ABSTRACT

During recent decades, different formulas have been developed to estimate longshore sediment transport rates through calibration using a wide variety of datasets, applicable for a range of particular wave and beach conditions. The equations that have shown the best capability to predict Bulk Longshore Sediment Transport Rate (BLSTR) are the formulas derived by CERC and by Kamphuis. In the present study, the five process parameters as used in the Kamphuis formula are accepted. The CERC formula includes only two of the five process parameters used in Kamphuis’ formula. A renewed optimization to derive the power of the five Kamphuis’ process parameters using an extensive dataset by Bayram was performed by Mil-Homens. In addition to this valuable effort, our contribution introduces two innovations. Firstly, the power coefficients of the five Kamphuis process parameters are optimized using a broad range of meta-heuristic algorithms. Secondly, the optimization is not based on the Bayram dataset as carefully collected and reviewed from published manuscripts but on a methodologically more homogeneous Iranian dataset acquired for port design and port management purposes. Independently from the results by Mil-Homens derived from the Bayram dataset, our study confirms these findings based on a totally different dataset. Specifically, the weaker impact of wave period and the stronger impact of the median grain diameter are in accordance with each other. The latter finding provides a stronger support for the mutual cancellation of the impact of slope and grain diameter in BLSTR, lending explanatory support to the CERC formula once beach slope and grain size are not known.

Descriptors:
Longshore sediment transport; Meta-heuristic optimization; Persian Gulf; Caspian Sea

INTRODUCTION

An important cause of coastline erosion and accretion is due to alongshore gradients in longshore sediment transport (CERC, 1984CERC (Coastal Engineering Research Center). 1984. Shore Protection Manual - Waterways Experiment Station. Washington, DC: Coastal Engineering Research Center., Fernandez et al., 2015FERNANDEZ, S., BAPITSA, P., MARTINS, V., SILVA, P., ABREU, T., PAIS-BARBOSA, J., BERNARDES, C., MIRANDA, P. V. L., ROSHA, M., SANTOS, F., BERNABEU, A. & REY, D. 2015. Longshore transport estimation on Ofir Beach in Northwest Portugal: sand-tracer experiment. Journal of Waterway, Port, Coastal, Ocean Engineering, 142(2), DOI: https://doi.org/10.1061/(ASCE)WW.1943-5460.0000319
https://doi.org/10.1061/(ASCE)WW.1943-54...
). Alongshore sediment transport gradients on coasts, not interrupted by headlands, outcrops or inlets, and caused by structural interventions, such as ports, are generally an order of magnitude larger than gradients caused by variations in natural coastline orientation (Bosboom and Stive, 2021BOSBOOM, J. & STIVE, M. J. F. 2021. Coastal dynamics. Netherlands: Delft University of Technology, DOI: https://doi.org/10.5074/T.2021.001
https://doi.org/10.5074/T.2021.001...
). Also, there exists an order-of-magnitude difference in longshore sediment transport rates between relatively higher or lower longshore transports in sandy or shingle environments, respectively. Hence, accurate estimates of the longshore sediment transport rates, especially for sandy coastal systems, are of high importance both for engineering and for managing these systems in a sustainable way.

During recent decades, different formulas have been developed to estimate the bulk longshore sediment transport rate (BLSTR). Formulas that have shown the best capability to predict BLSTR are the formulas of CERC (CERC, 1984CERC (Coastal Engineering Research Center). 1984. Shore Protection Manual - Waterways Experiment Station. Washington, DC: Coastal Engineering Research Center.), Kamphuis (Kamphuis, 1991KAMPHUIS, J. W. & READSHAW, J. S. 1978. A model study of alongshore sediment transport rates. Coastal Egineergin Proceedings, 1(16), 99, DOI: https://doi.org/10.9753/icce.v16.99
https://doi.org/10.9753/icce.v16.99...
) and Kamphuis-MilHomens (Mil-Homens, 2016MIL-HOMENS, J. 2016. Longshore sediment transport: bulk formulas and process based models. PhD. Netherlands: Delft University of Technology, DOI: https://doi.org/10.4233/uuid:f7703aba-2760-47b1-99c2-a076a84d9b0f
https://doi.org/10.4233/uuid:f7703aba-27...
). These formulas have been calibrated on broadly differing datasets applicable for a range of particular wave and beach conditions (Mil-Homens et al., 2013MIL-HOMENS, J., RANASINGHE, R., VAN THIEL DE VRIES, J. S. M. & STIVE, M. J. F. 2013. Influence of profile features on longshore sediment transport. In: Coastal Dynamics, 7th International Conference on Coastal Dynamics, 24-28 June, Arcachon, France , Bordeaux University.). Bayram et al. (2007)BAYRAM, A., LARSON, M. & HANSON, H. 2007. A new formula for the total longshore sediment transport rate. Coastal Engineering, 54(9), 700-710, DOI: https://doi.org/10.1016/j.coastaleng.2007.04.001
https://doi.org/10.1016/j.coastaleng.200...
shared the reviewed and extended Bayram dataset with Mil-Homens (2016)MIL-HOMENS, J. 2016. Longshore sediment transport: bulk formulas and process based models. PhD. Netherlands: Delft University of Technology, DOI: https://doi.org/10.4233/uuid:f7703aba-2760-47b1-99c2-a076a84d9b0f
https://doi.org/10.4233/uuid:f7703aba-27...
, who used the Bayram dataset to reevaluate the Kamphuis formula (here indicated as Kamphuis-MilHomens). The accuracy of the estimation of BLSTR is assumed to be suffering from the data being rather inhomogeneous due to methodological differences in data collection, believed to be one reason for an order-of-magnitude difference in estimates (Mil-Homens, 2016MIL-HOMENS, J. 2016. Longshore sediment transport: bulk formulas and process based models. PhD. Netherlands: Delft University of Technology, DOI: https://doi.org/10.4233/uuid:f7703aba-2760-47b1-99c2-a076a84d9b0f
https://doi.org/10.4233/uuid:f7703aba-27...
).

As observed by Fernandez et al. (2015)FERNANDEZ, S., BAPITSA, P., MARTINS, V., SILVA, P., ABREU, T., PAIS-BARBOSA, J., BERNARDES, C., MIRANDA, P. V. L., ROSHA, M., SANTOS, F., BERNABEU, A. & REY, D. 2015. Longshore transport estimation on Ofir Beach in Northwest Portugal: sand-tracer experiment. Journal of Waterway, Port, Coastal, Ocean Engineering, 142(2), DOI: https://doi.org/10.1061/(ASCE)WW.1943-5460.0000319
https://doi.org/10.1061/(ASCE)WW.1943-54...
, both the CERC and Kamphuis formulas can be expressed by a generic equation that includes all physical parameters that have been suggested to determine the longshore sediment transport on sandy beaches, but with different power values for the parameters. In the present study, the same generic equation is adopted.

Two innovations will be introduced here in comparison to the aforementioned studies. Firstly, the power coefficients of the most relevant BLSTR estimates are not optimized with traditional optimization algorithms but with twelve different meta-heuristic algorithms, with the expectation that this leads to improved accuracy or at least improved insights in the accuracy. Secondly, the optimization is not based on the extensive existing database from published manuscripts but rather on a new and methodologically more homogeneous Iranian dataset acquired for port design, engineering and management purposes. The results lead to converging insights regarding the accuracy of the earlier derivations of the data-fit for bulk longshore sediment transport rates.

METHODS

BLSTR FORMULATION REVIEWED

In this section, two major BLSTR formulations are reviewed based on their success as perceived from published manuscripts. All suggestions made in the past are still in place, and this confirms that this issue is fairly difficult to improve. This research is focused on increasing insight into the accuracy of BLSTR. Alternative approaches, detailing the longshore sediment transport as a function of the cross-shore location such as LITPACK (Deigaard et al., 1986DEIGAARD, R., FREDSOE, J. & HEDEGAARD, B. 1986. Mathematical model for littoral drift. Journal of Waterways, Ports, Coastal and Ocean Engineering, 112(3), 351-369. ) and UNIBEST (Bosboom and Stive, 2021BOSBOOM, J. & STIVE, M. J. F. 2021. Coastal dynamics. Netherlands: Delft University of Technology, DOI: https://doi.org/10.5074/T.2021.001
https://doi.org/10.5074/T.2021.001...
) are valuable once they incorporate additional data that are not considered in the BLSTR approach. In cases where these data do not exist, it does not necessarily lead to better accuracy to estimate BLSTR. Hence, there is practical justification for the use of BLSTR in engineering and management projects (Bosboom and Stive, 2021).

CERC FORMULA

The CERC formula (introduced in earlier editions of the Shore Protection Manual (CERC, 1984)) is the “mother” of all formulations derived to estimate the BLSTR. This formula was obtained based on the assumption that the BLSTR is proportional to the “longshore component of wave power”. The physical interpretation of this parameter remained unclear until 1970 when Longuet-Higgins (1970) presented a close relation of this parameter with the cross-shore gradient of the alongshore radiation stress integrated over the surfzone. Using the wave height and wave incidence angle at the breakerline, the CERC formula can be written as follows, where the k coefficient is dimensional and BLSTR (Q) is in (m ^ 3/s) (Mil-Homens et al., 2016):

(1) Q = k ρ g 1 / 2 16 γ b ρ s ρ 1 ρ H sb 5 / 2 sin 2 θ b

It is noted that this formula has only two physical process parameters lacking including sediment grain size.

Under the above conditions, the k coefficient was obtained through linear regression by Komar and Inman (1970)KOMAR, P. D. & INMAN, D. L. 1970. Longshore sand transport on beaches. Journal of Geophysical Research, 75(30), 5914-5927, DOI: https://doi.org/10.1029/JC075i030p05914
https://doi.org/10.1029/JC075i030p05914...
from field measurements on Silver Strand (USA) and at El Moreno (Mexico), resulting in k = 0.77 based on the root mean square of the wave height at the breaking point or k = 0.34 using the significant wave height at the breakpoint. CERC (1984) updated the k value in 1984 using an additional dataset. However, the remaining scatter (Rosati et al., 2002ROSATI, L. W. T. & KEVI, B. 2002. Coastal Engineering Manual. Engineering Manual (EM) 1102-1100. Coastal Sediment process. Washington, DC: U.S. Army Corps of Engineers. ) suggests that k should not be a constant but a function of other potentially relevant physical parameters. Below, several suggestions are reviewed but none of these have been accepted on a general basis.

Indeed, even before the earlier suggestion, several attempts were made to find the impact of potentially important parameters. The relationship of BLSTR with the surf similarity parameter also known as the Iribarren number (Battjes, 1974BATTJES, J. A. 1974. Surf similarity. In: Proceedings of the 14th International Conference on Coastal Engineering, 24-28 June, Copenhagen, Denmark, American Society of Civil Engineers, pp. 467-479, DOI: https:/doi.org/10.1061/9780872621138.029
https:/doi.org/10.1061/9780872621138.029...
) was suggested by Kamphuis and Readshaw (1987). A second attempt to include more physical parameters in the expression for k was introduced through an energy-based model (Bailard, 1981BAILARD, J. A. 1981. An energetics total load sediment transport model for a plane sloping beach. Journal of Geophysical Research: Oceans, 86(C11), 10938-10954, DOI: https://doi.org/10.1029/JC086iC11p10938
https://doi.org/10.1029/JC086iC11p10938...
). This author suggested the value of k should increase with larger incident wave angles and larger values of the Dean parameter.

A third attempt was due to Komar (1998)KOMAR, P. D. 1998. Longshore currents generated by obliquely incident seawaves. Journal of Geophysical Research, 75(33), 6778-6789., who critically reviewed the BLSTR data and introduced data from the Adra River Delta in Spain. The wide range of sediment grains ranging from 0.4 to 1.5 mm led to an experimental relation based on the mean diameter of grains.

Later evaluations by Schoones and Theron (1996) and Smith et al. (2009)SMITH, E. R., WANG, P., EBERSOLE, B. A. & ZHANG, J. 2009. Dependence of total longshore sediment transport rates on incident wave parameters and breaker type. Journal of Coastal Research, 25(3), 675-683.concluded that the first attempt to relate the value of k to the surf similarity parameter was the most successful. However, these suggestions, although interesting, are considered not relevant in the context of the present work.

KAMPHUIS FORMULA

Kamphuis (1991) developed a formula for BLSTR for an extensive range of data that included all earlier field data and a large new laboratory dataset acquired by that author, which may have led to a bias. Through dimensional analysis and using physics principles, this author obtained the following equation for the bulk volume of longshore sediment transport Im in kg/s and valid for regular waves:

(2) I m ρ H 3 T = K * H L 0 * P m b q H D 50 sin s 2 b

where H: wave height, T: wave period, d: depth of water, ρ: fluid density, m: beach slope, p and q and r and s are experimentally determined, k* is a calibration coefficient, L*0 is the wavelength in deep water (for regular waves), and D50 is the median sediment grain diameter. The beach slope in the surf zone region is calculated according to:

(3) m b = d b λ b

where λb is the distance from the coastline to the breaking point. The bulk longshore sediment transport rate is expressed by:

(4) Q 1 I m ρ s ρ 1 p

For irregular waves, where the subscript b refers to the wave parameters at the breakerline, Kamphuis (1991) also presented the following formula:

(5) Im ρ H sb 3 Tp = 1 . 3 * 10 3 H sb L 0 1 . 25 m b 0 . 75 H sb D 50 0 . 25 sin 0 . 6 2 b

where Hsb is the wave height at the breakerline, Tp is the wave period and L0=gTp22π is for deep water.

The above relation in a simplified form with k=k*ρgTp22π1.25=2.27 can be written as follows:

(6) I M = 2 . 27 H sb 2 T p 0 . 75 D 50 0 . 25 sin 0 . 6 2 b

However, using the same data, a new value for the calibration coefficient is equivalent to the values ks&T=2.27 for QI (m3/s)and K*=50,000 in terms of m3/year.

Kamphuis (1991), besides wave height and wave incidence at the breakpoint, considered three more process parameters, viz. coastal slope, wave period and median sediment grain size; in the CERC formula, these latter parameters are not considered. Interestingly, coastal slope and sediment grain size are related using the positive correlation between bed slope and the median grain size (Klein, 2003). As grain size decreases the coastal slope decreases, so these effects can cancel each other out, which is often explained as the justification for the absence of these parameters in the CERC formula.

RECENT OPTIMIZATIONS

Using the improved Bayram dataset, Mil-Homens (2013) derived best fits for CERC and Kamphuis with improved results in terms of bias and RMSE compared to all earlier studies mentioned.

The best fit for the CERC formula was achieved with a polynomial function,

(7) K CERC = 2232 . 7 H sb L 0 1 . 45 + 4 . 505 1

The best fit achieved for the Kamphuis formula reads:

(8) I m . new = 0 . 149 H sb 2 . 75 T p 0 . 89 D 50 069 sin 0 . 5 2 b

The above results were based on the extensive Bayram dataset; although accuracy improved, still only 53% to 56% of the data points fell within a factor of between 0.5 and 2 from the best fit. This can be attributed to the data being both inhomogeneous due to hydrodynamics, morphological and geological differences in the regions where the data are generated and the methods with which these are generated, leading to an order-of-magnitude difference in estimates.

CERC, Kamphuis and Kamphuis-Milhomens can all be expressed in the following form (Fernandez et al., 2015FERNANDEZ, S., BAPITSA, P., MARTINS, V., SILVA, P., ABREU, T., PAIS-BARBOSA, J., BERNARDES, C., MIRANDA, P. V. L., ROSHA, M., SANTOS, F., BERNABEU, A. & REY, D. 2015. Longshore transport estimation on Ofir Beach in Northwest Portugal: sand-tracer experiment. Journal of Waterway, Port, Coastal, Ocean Engineering, 142(2), DOI: https://doi.org/10.1061/(ASCE)WW.1943-5460.0000319
https://doi.org/10.1061/(ASCE)WW.1943-54...
):

(9) Q S = C H sb a 1 T p a m b a 3 D 50 a 4 sin a 5 2 θ b ρ s ρ 1 a

where BLSTR = Qs (m3/s) is the volume rate of longshore sediment transport, θb is the angle between crest and the coast normal at the initial breaking point, mb is the coast slope from the initial breakerline to the coastline, ρs sediment density is 2650 kg/m3, ρ water density is 1025 kg/m3 and the porosity index is 0.4.

Initially, coefficients for all sediment transport formulas, including all coefficients of α1, α2, α3, α4 and α5 and also C related to the CERC, Kamphuis and Kamphuis-MilHomens formulas were gathered totally in view of Equation 20. All these coefficients are presented in summary in Table 1.

Table 1
Coefficients and exponents of BLSTR formulas (modified and extended after Fernandez et al., 2015).

The experimental coefficient k for the CERC formula is 0.39 based on the initial breaking wave height parameter. Furthermore, g is acceleration due to gravity, yb is the breaker index, Hb/hb in which hb is the depth at the break point. Cf is a coefficient of friction that is 0.005 according Bayram formula (Mil-Homens, 2013) and ε is coefficient of sediment material transport that has been estimated in this formula.

In this study, the coefficient and exponents of Equation 9 are optimized with twelve different meta-heuristic algorithms, using a range of different optimization methods and a new dataset.

DATA AND DATA ANALYSIS METHODS

The extensive Bayram data set is very relevant and should not be ignored, but the data suffer from being inhomogeneous in terms of data-derivation techniques. Therefore, this study investigated whether our understanding of the accuracy improves by using a more homogeneous field dataset, i.e. data based on BLSTR field data derived according to similar methodologies. It is our hypothesis that a more homogeneous dataset could lead to improved accuracy and/or converging views.

The Iranian Ports and Maritime Organization (PMO, 2021PMO (Iranian Port and Maritime Organization). 2021. Homepage [online]. Iran: Iranian Port and Maritime Organization. Available at: https://irancoasts.pmo.ir/en/first [Accessed 17 April 2021].
https://irancoasts.pmo.ir/en/first...
) partly collected (through field observations) and partly constructed (through calibrated and validated numerical modelling) a more homogeneous dataset in the context of Integrated Coastal Zone Management Plans evaluating engineering interventions for 12 different ports in Iran (Figure 1, Figure 2 and Table 2). Comprehensive field measurements were conducted in four provinces both along uninterrupted stretches of the provincial coast and in the direct vicinity of 12 ports in these provinces, of which four are in the southern part of the Caspian Sea and eight in the Persian Gulf. The data include not only the relevant parameters in the BLSTR formulae such as grain size, mean beach slope, average wave period, mean wave height and dominant wave angle at the breaking point, but also detailed information on incident and cross-shore wave heights, periods, frequency of occurrence (see Figures 3, 4), longshore flows, cross-shore bed profiles and bottom roughness information. The collected detailed information justifies the use of a cross-shore resolving longshore wave, flow and sediment transport model for calibration. The state-of-the-art numerical model LITPACK (Deigaard et al., 1986DEIGAARD, R., FREDSOE, J. & HEDEGAARD, B. 1986. Mathematical model for littoral drift. Journal of Waterways, Ports, Coastal and Ocean Engineering, 112(3), 351-369. ) was applied. The calibration procedure was based on comprehensive modelling for the large coastal stretches in the four provinces where the ports are located. While for all variables default values were adopted, the only variables in the calibration were bed roughness and incident wave energy. (Mahmoodi et al., 2020MAHMOODI, A., LASHTEHNESHAEI, M. A., MANSOURI, A. & SHAFAEIBAJESTAN, M. 2020. Study of current- and wave-induced sediment transport in the Nowshahr port entrance channel by using numerical modeling and field measurements. Journal of Marine Science and Engineering, 8(4), 284, DOI: https://doi.org/10.3390/jmse8040284
https://doi.org/10.3390/jmse8040284...
). After calibration on the uninterrupted coastal stretches in the provinces, the LITPACK model was validated primarily based on satellite-based temporal areal changes on the updrift side of the ports (Jafari et all., 2016JAFARI, E., ALAEE, M. J., NAZARALI, M. & BALI, M. 2016. Sea water temperature observation and simulation in the Caspian Sea. In: Proceedings of the International Conference on Ports, Marine and Structures (ICOPMAS), 31 October - 2 November, Tehran, Iran, ICOPMAS.), without changing the model parameters.

Figure 1
Locations and dominant wave directions at the initial breaking point of 4 PMO managed Iranian ports in the Caspian Sea.

Table 2
Field data of different sites for longshore transport at different parts of the coast of Iran, after data analysis for sensitivity.

Figure 2
Locations and dominant wave directions at the initial breaking point of 8 PMO managed Iranian ports in the Persian Gulf.

Figure 3
Dominant wave directions at the initial breaking point of 4 Iranian ports in Caspian Sea.

Figure 4
Dominant wave directions at the initial breaking point of 8 Iranian ports in Persian Gulf.

The total number of data after analysis and normalization consists of more than 2000 items, describing detailed information on BLSTRs as a function of the relevant process parameters. These data are considered first to form an alternative, more homogenous basis for optimization other than the Bayram dataset and second an opportunity to explore meta-heuristic optimization algorithms as alternative, very efficient algorithms relative to the traditional optimizations. This allows using a wide range of such algorithms that should shed new light on the variability of the powers of the process parameters.

Meta-heuristic algorithms are a group of algorithms used to solve optimization challenges. In the field of computer sciences, artificial intelligence and optimization problems, these algorithms are a way to achieve optimization that classical solutions are slow to solve, or even more serious an approximate solution to problems that classical ways cannot find exact answers to. Most complex problems require evaluating a myriad of possible modes to determine an exact answer. Meta-heuristic algorithms play an effective role in resolving such issues by using methods that require less evaluation and provide answers within acceptable time limits. Meta-heuristic algorithms are general algorithmic frameworks that can provide specific solutions to the optimization problem. These algorithms belong to the class of various approximate optimization algorithms having solutions to exit from local optimal points and have the capability to be applied to an extensive spectrum of problems. Meta-heuristic algorithms use large-scale problems to provide satisfactory solutions within a reasonable amount of time. The common goal of all meta-heuristic algorithms is to solve difficult optimization problems (Dreo et al., 2006DREO, J., SIARRY, P., PETROWSKI, A. & TAILLARD, E. 2006. Metaheuristics for hard optimization: methods and case studies. Berlin: Springer.).

Meta-heuristic methods have the following common features. They are commonly applied to solve hybrid problems but can also be used in continuous problems. They are usually inspired by the concepts of biology, animal behavior and physics. One of the common drawbacks of these methods is the difficulty of setting and matching the parameters. Different criteria can be used to classify meta-heuristic algorithms, i.e. answer-based and population-based. Answer-based algorithms change a response during the search process, while population-based algorithms consider a population of responses during a search. Many meta-heuristic algorithms are inspired by nature. Some meta-heuristic algorithms lack memory, meaning that these types of algorithms do not use the information obtained during the search, while some search algorithms do. They use memory storing the information obtained during the search. Exact and possible: a definitive meta-heuristic algorithm solves the problem using definite decisions. However, in some meta-heuristic algorithms, a series of possible rules are used during the search.

In total, 12 meta-heuristic algorithms were explored that are among the most important ones inspired by the behavior of nature (Yang, X.S, 2010YANG, X. S. 2010. A new meta-heuristic bat-inspired algorithm. In: GONZÁLES, J. R., PELTA, D. A., CRUZ, C., TERRAZAS, G. & KRASNOGOR, N. (eds.). Nature inspired cooperative strategies for optimization (NICSO 2010). Berlin: Springer-Verlag, pp. 65-74, DOI: https://doi.org/10.1007/978-3-642-12538-6_6
https://doi.org/10.1007/978-3-642-12538-...
). Using the 12 meta-heuristic algorithms, an optimization of the BLSTR coefficients and exponents was conducted for all field data (Gholami et al., 2021GHOLAMI, Z., LARI, K., BIDOKHTI, A. A. & JAVID, A. 2021. Calculation of uniform longshore sediment transport rate in coastal zone using incomplete self-similarity theory. Hydrophysics [online], Article in Press. Available at: https://www.hydrophysics.ir/article_243538.html?lang=en [Accessed 17 April 2021].
https://www.hydrophysics.ir/article_2435...
). The coefficients and optimal capacities were extracted using meta-heuristic algorithms in MATLAB space for the data as summarized in Table 2.

From a more profound evaluation of the results, it is concluded that three algorithms, i.e. Ant Colony (ACO), Bath (BAT) and Whale (WOA), performed best. Although the results for all algorithms were close, these algorithms showed a fast convergence between the calculated and observed values, which is a good indicator of their optimization reliability. For reasons of conciseness only the results for the selected three algorithms are presented. These algorithms use swarm intelligence as a relatively new approach to problem solving that takes inspiration from the social behaviors of insects and of other animals.

Ant colony optimization (ACO) is a population-based meta-heuristic that can be used to find approximate solutions to difficult optimization problems (Dorigo et al, 2006DORIGO, M., BIRATTARI, M. & STUTZLE, E. 2006. Ant colony optimization. IEEE Computational Intelligence Magazine, 1(4), 28-39, DOI: https://doi.org/10.1109/MCI.2006.329691
https://doi.org/10.1109/MCI.2006.329691...
). Ant colony optimization takes inspiration from the foraging behavior of some ant species. These ants deposit pheromone on the ground in order to mark some favorable path that should be followed by other members of the colony. Ant colony optimization exploits a similar mechanism for solving optimization problems.

The bat optimization algorithm (BAT) is a novel meta-heuristic optimization algorithm introduced by Yang (2010)YANG, X. S. 2010. A new meta-heuristic bat-inspired algorithm. In: GONZÁLES, J. R., PELTA, D. A., CRUZ, C., TERRAZAS, G. & KRASNOGOR, N. (eds.). Nature inspired cooperative strategies for optimization (NICSO 2010). Berlin: Springer-Verlag, pp. 65-74, DOI: https://doi.org/10.1007/978-3-642-12538-6_6
https://doi.org/10.1007/978-3-642-12538-...
. The algorithm is inspired by the echolocation behavior of bats, which guides the bats’ movements during their flight and hunting even in complete darkness. This algorithm carries out the search process using artificial bats as search agents mimicking the natural pulse loudness and emission rate of real bats. (Fattahi et al., 2007FATTAHI, P., SAIDIMEHRABAD, M. & JOLAI, F. 2007. Mathematical modeling and heuristic approaches to flexible job shop scheduling problems. Journal of Intelligent Manufacturing, 18, 331-342, DOI: https://doi.org/10.1007/s10845-007-0026-8
https://doi.org/10.1007/s10845-007-0026-...
).

The whale optimization algorithm (WOA) is a biological heuristic algorithm presented by Mirjalili and Lewis (2016)MIRJALILI, S. & LEWIS, A. 2016. The whale optimization algorithm. Advances in Engineering Software, 95, 51-67, DOI: https://doi.org/10.1016/j.advengsoft.2016.01.008
https://doi.org/10.1016/j.advengsoft.201...
. WOA is a swarm intelligence optimization algorithm inspired by the unique humpback hunting method. Because of this unique optimization mechanism, WOA has a good global search capability. Therefore, the new algorithm has been widely proposed by the engineering community.

Table 3 lists some of the properties that characterize these three algorithms.

Table 3
Properties of the 3 selected meta-algorithms.

RESULTS

The meta-heuristic optimization results for the three selected algorithms are presented in Table 4 and compared to the Kamphuis and Kamphuis-MilHomens optimizations of Eq. 9. As a first observation it is promising to conclude that there is only moderate quantitative variation between these results mutually.

Table 4
Coefficients and exponents of Equation 12 for Kamphuis, Kamphuis-Milhomens using the Bayram data-set and for the 3 selected meta-heuristic algorithms) using the field data of table 2.

In order to check whether one of the algorithms performs better, the RMSE for these three algorithms was derived for the dataset of Manzadaran Province, considered to be one of the most homogeneous due to the monodirectionality of the alongshore sediment transport.

(10) RMSE = i = n log Qp . i log Qm . i 2 n

The RMSE value is a commonly used error measure. The sum of squares gives more weight to higher error values, and consequently higher error variances. The bias provides insight on any systematic offset of the data. Because logarithmic values (base 10) are considered in both statistical measures, the values indicate errors in terms of magnitude order, e.g., an RMSE value of 1 would mean that the predicted values are roughly, on average, 10 times larger or smaller than the measured ones. Logarithmic values were chosen because the data range extends through several orders of magnitude (Mil-Homens, 2013MIL-HOMENS, J., RANASINGHE, R., VAN THIEL DE VRIES, J. S. M. & STIVE, M. J. F. 2013. Re-evaluation and improvement of three commonly used bulk longshore sediment transport formulas. Coastal Engineering , 75, 29-39, DOI: https://doi.org/10.1016/j.coastaleng.2013.01.004
https://doi.org/10.1016/j.coastaleng.201...
).

The three algorithms, i.e., Ant Colony (ACO), Bath (BAT) and Whale (WOA), performed very similar with an RMSE of 0.012, 0.011 and 0.010 respectively, which does not imply a significant difference between the algorithms. These low values indicate that the dataset is very homogeneous in methodology and leads to two important conclusions. First, a more homogeneous dataset leads to more reliable, quantitative results for BLSTR formulations. Second, different combinations of values for the coefficients and the exponents in Equation 20 can lead to similar results for a homogeneous dataset.

DISCUSSION

During recent decades, different formulas have been developed to estimate BLSTR through calibration on regionally different and therefore often non-homogeneous datasets. Equations that have shown the best capability to predict BSLRT are the formulas of CERC, Kamphuis and Kamphuis-MilHomens. In the present study, the same process parameters as used in the Kamphuis and Kamphuis-MilHomens formulas are accepted. However, two innovations are introduced here. Firstly, the coefficients and exponents of the different formulations are optimized with twelve meta-heuristic algorithms, using a range of different optimization methods. Secondly, the optimization is not based on existing data from literature but on a unique, new and rather homogeneous Iranian dataset acquired for port design and port management purposes. Using meta-heuristic algorithms, a first optimization of the BLSTR coefficients and exponents was conducted for the whole Iranian field dataset. An evaluation of the results on the basis of two criteria concludes that three algorithms, i.e. Ant Colony (ACO), Bath (BAT) and Whale (WOA), performed best. These algorithms showed both a fast convergence between the calculated and observed values and are close to the quantitative values of the coefficients and exponents of the most recent formulation and optimization of Kamphuis-MilHomens. The three mentioned algorithms performed very similar with an RMSE of 0.012, 0.011 and 0.010 respectively. These low values indicate that the dataset is very homogeneous and leads to two important conclusions. First, a more homogeneous dataset leads to more reliable, quantitative results for BLSTR formulations. Second, different combinations of values for the coefficients and the exponents can lead to similar results for a homogeneous dataset.

The results lead to the conclusion that independently from the results by Mil-Homens (2016)MIL-HOMENS, J. 2016. Longshore sediment transport: bulk formulas and process based models. PhD. Netherlands: Delft University of Technology, DOI: https://doi.org/10.4233/uuid:f7703aba-2760-47b1-99c2-a076a84d9b0f
https://doi.org/10.4233/uuid:f7703aba-27...
derived from the Bayram dataset, the present study confirms these findings based on a totally different dataset and optimization methodology. Specifically, the weaker impact of wave period and the stronger impact of the median grain diameter are in accordance with each other. The latter finding provides a stronger support for the mutual cancellation of the impact of slope and grain diameter in BLSTR, lending explanatory and converging support to the CERC formula once beach slope and grain size are not known.

ACKNOWLEDGMENTS

The authors would like to acknowledge the Iran Port and Maritime Organization (PMO) for providing the data. The authors would also like to thank Dr Bayram and Dr Joao Mil-Homens for sharing their data. The authors are grateful for the thorough comments by two anonymous reviewers.

REFERENCES

  • BAILARD, J. A. 1981. An energetics total load sediment transport model for a plane sloping beach. Journal of Geophysical Research: Oceans, 86(C11), 10938-10954, DOI: https://doi.org/10.1029/JC086iC11p10938
    » https://doi.org/10.1029/JC086iC11p10938
  • BATTJES, J. A. 1974. Surf similarity. In: Proceedings of the 14th International Conference on Coastal Engineering, 24-28 June, Copenhagen, Denmark, American Society of Civil Engineers, pp. 467-479, DOI: https:/doi.org/10.1061/9780872621138.029
    » https:/doi.org/10.1061/9780872621138.029
  • BAYRAM, A., LARSON, M. & HANSON, H. 2007. A new formula for the total longshore sediment transport rate. Coastal Engineering, 54(9), 700-710, DOI: https://doi.org/10.1016/j.coastaleng.2007.04.001
    » https://doi.org/10.1016/j.coastaleng.2007.04.001
  • BOSBOOM, J. & STIVE, M. J. F. 2021. Coastal dynamics Netherlands: Delft University of Technology, DOI: https://doi.org/10.5074/T.2021.001
    » https://doi.org/10.5074/T.2021.001
  • CERC (Coastal Engineering Research Center). 1984. Shore Protection Manual - Waterways Experiment Station Washington, DC: Coastal Engineering Research Center.
  • DEIGAARD, R., FREDSOE, J. & HEDEGAARD, B. 1986. Mathematical model for littoral drift. Journal of Waterways, Ports, Coastal and Ocean Engineering, 112(3), 351-369.
  • DORIGO, M., BIRATTARI, M. & STUTZLE, E. 2006. Ant colony optimization. IEEE Computational Intelligence Magazine, 1(4), 28-39, DOI: https://doi.org/10.1109/MCI.2006.329691
    » https://doi.org/10.1109/MCI.2006.329691
  • DREO, J., SIARRY, P., PETROWSKI, A. & TAILLARD, E. 2006. Metaheuristics for hard optimization: methods and case studies Berlin: Springer.
  • FATTAHI, P., SAIDIMEHRABAD, M. & JOLAI, F. 2007. Mathematical modeling and heuristic approaches to flexible job shop scheduling problems. Journal of Intelligent Manufacturing, 18, 331-342, DOI: https://doi.org/10.1007/s10845-007-0026-8
    » https://doi.org/10.1007/s10845-007-0026-8
  • FERNANDEZ, S., BAPITSA, P., MARTINS, V., SILVA, P., ABREU, T., PAIS-BARBOSA, J., BERNARDES, C., MIRANDA, P. V. L., ROSHA, M., SANTOS, F., BERNABEU, A. & REY, D. 2015. Longshore transport estimation on Ofir Beach in Northwest Portugal: sand-tracer experiment. Journal of Waterway, Port, Coastal, Ocean Engineering, 142(2), DOI: https://doi.org/10.1061/(ASCE)WW.1943-5460.0000319
    » https://doi.org/10.1061/(ASCE)WW.1943-5460.0000319
  • GHOLAMI, Z., LARI, K., BIDOKHTI, A. A. & JAVID, A. 2021. Calculation of uniform longshore sediment transport rate in coastal zone using incomplete self-similarity theory. Hydrophysics [online], Article in Press. Available at: https://www.hydrophysics.ir/article_243538.html?lang=en [Accessed 17 April 2021].
    » https://www.hydrophysics.ir/article_243538.html?lang=en
  • JAFARI, E., ALAEE, M. J., NAZARALI, M. & BALI, M. 2016. Sea water temperature observation and simulation in the Caspian Sea. In: Proceedings of the International Conference on Ports, Marine and Structures (ICOPMAS), 31 October - 2 November, Tehran, Iran, ICOPMAS.
  • KAMPHUIS, J. W. & READSHAW, J. S. 1978. A model study of alongshore sediment transport rates. Coastal Egineergin Proceedings, 1(16), 99, DOI: https://doi.org/10.9753/icce.v16.99
    » https://doi.org/10.9753/icce.v16.99
  • KOMAR, P. D. 1998. Longshore currents generated by obliquely incident seawaves. Journal of Geophysical Research, 75(33), 6778-6789.
  • KOMAR, P. D. & INMAN, D. L. 1970. Longshore sand transport on beaches. Journal of Geophysical Research, 75(30), 5914-5927, DOI: https://doi.org/10.1029/JC075i030p05914
    » https://doi.org/10.1029/JC075i030p05914
  • MAHMOODI, A., LASHTEHNESHAEI, M. A., MANSOURI, A. & SHAFAEIBAJESTAN, M. 2020. Study of current- and wave-induced sediment transport in the Nowshahr port entrance channel by using numerical modeling and field measurements. Journal of Marine Science and Engineering, 8(4), 284, DOI: https://doi.org/10.3390/jmse8040284
    » https://doi.org/10.3390/jmse8040284
  • MIL-HOMENS, J. 2016. Longshore sediment transport: bulk formulas and process based models. PhD. Netherlands: Delft University of Technology, DOI: https://doi.org/10.4233/uuid:f7703aba-2760-47b1-99c2-a076a84d9b0f
    » https://doi.org/10.4233/uuid:f7703aba-2760-47b1-99c2-a076a84d9b0f
  • MIL-HOMENS, J., RANASINGHE, R., VAN THIEL DE VRIES, J. S. M. & STIVE, M. J. F. 2013. Re-evaluation and improvement of three commonly used bulk longshore sediment transport formulas. Coastal Engineering , 75, 29-39, DOI: https://doi.org/10.1016/j.coastaleng.2013.01.004
    » https://doi.org/10.1016/j.coastaleng.2013.01.004
  • MIL-HOMENS, J., RANASINGHE, R., VAN THIEL DE VRIES, J. S. M. & STIVE, M. J. F. 2013. Influence of profile features on longshore sediment transport. In: Coastal Dynamics, 7th International Conference on Coastal Dynamics, 24-28 June, Arcachon, France , Bordeaux University.
  • MIRJALILI, S. & LEWIS, A. 2016. The whale optimization algorithm. Advances in Engineering Software, 95, 51-67, DOI: https://doi.org/10.1016/j.advengsoft.2016.01.008
    » https://doi.org/10.1016/j.advengsoft.2016.01.008
  • PMO (Iranian Port and Maritime Organization). 2021. Homepage [online]. Iran: Iranian Port and Maritime Organization. Available at: https://irancoasts.pmo.ir/en/first [Accessed 17 April 2021].
    » https://irancoasts.pmo.ir/en/first
  • ROSATI, L. W. T. & KEVI, B. 2002. Coastal Engineering Manual. Engineering Manual (EM) 1102-1100. Coastal Sediment process Washington, DC: U.S. Army Corps of Engineers.
  • SMITH, E. R., WANG, P., EBERSOLE, B. A. & ZHANG, J. 2009. Dependence of total longshore sediment transport rates on incident wave parameters and breaker type. Journal of Coastal Research, 25(3), 675-683.
  • YANG, X. S. 2010. A new meta-heuristic bat-inspired algorithm. In: GONZÁLES, J. R., PELTA, D. A., CRUZ, C., TERRAZAS, G. & KRASNOGOR, N. (eds.). Nature inspired cooperative strategies for optimization (NICSO 2010) Berlin: Springer-Verlag, pp. 65-74, DOI: https://doi.org/10.1007/978-3-642-12538-6_6
    » https://doi.org/10.1007/978-3-642-12538-6_6

Edited by

Associate Editor: Eduardo Siegle

Publication Dates

  • Publication in this collection
    16 Aug 2021
  • Date of issue
    2021

History

  • Received
    25 Dec 2020
  • Accepted
    10 June 2021
Instituto Oceanográfico da Universidade de São Paulo Praça do Oceanográfico 191, CEP: 05508-120, São Paulo, SP - Brasil, Tel.: (11) 3091-6501 - São Paulo - SP - Brazil
E-mail: diretoria.io@usp.br