The increased pumping of freshwater from coastal aquifers, to meet growing demands, causes an environmental problem called saltwater intrusion. Consequently, proper management schemes are necessary to tackle such a situation and permit the optimal development of coastal groundwater basins. In this research, a probabilistic search algorithm, namely Probabilistic Global Search Lausanne (PGSL), is used to calculate optimal pumping rates of unconfined coastal aquifer. The results of using PGSL are compared with a stochastic search optimization technique, Shuffled Frog Leaping Algorithm (SFLA). The finite element method is applied to simulate the hydraulic response of the steady state homogenous aquifer. The lower and upper (LU) decomposition method is adapted to invert the conductance matrix, which noticeably decreases the computation time. The results of both the PGSL and the SFLA are verified through the application on the aquifer system underlain the City of Miami Beach in the north of Spain. Multiple independent optimization runs are carried out to provide more insightful comparison outcomes. Consequently, a statistical analysis is performed to assess the performance of each algorithm. The two optimization algorithms are then applied on the Quaternary Aquifer of El-Arish Rafah area, Egypt. The results show that both algorithms can effectively be used to obtain nearly global solutions compared with other previous published results.