Horizontal boundary identification of causative sources is an essential tool in potential field data interpretation due to the feasibility of automatically retrieving the boundary information of subsurface gravity or geomagnetic structures. Although many approaches have been proposed to address these issues, it is still a hot research topic for many researchers to derive novel methods or enhance existing techniques. We present two high-resolution edge detectors based on the Gudermannian function and the modifications of the second-order derivative of the field. The effectiveness of the newly proposed filters was initially tested on synthetic gravity anomalies and geomagnetic responses with different assumptions (2-D and 3-D; imposed and superimposed; noise-free and noise-contaminated). The obtained results verified that the two novel methods yield the capability of producing high-resolution, balanced amplitudes and accurate results for better imaging causative sources with different geometrical and geophysical properties, compared with the other nine representative edge enhancement techniques. Furthermore, the yielded results from the application of the two strategies to a real-world aeromagnetic data set measured from the Central Puget Lowland (C.P.L) of the United States and a gravity data set surveyed from the Jalal Abad area of Kerman province, Iran, with detailed comparative studies validated that the edges identified via the two methods are in good agreement with the major geological structures within the study areas and the determined lateral information using the tilt-depth, top-depth estimation method. These features make them valuable tools for solving edge detection problems.

Dealing with the ill-posed and non-unique nature of the non-linear geophysical inverse problem via local optimizers requires the use of some regularization methods, constraints, and prior information about the Earth's complex interior. Another difficulty is that the success of local search algorithms depends on a well-designed initial model located close to the parameter set providing the global minimum. On the other hand, global optimization and metaheuristic algorithms that have the ability to scan almost the entire model space do not need an assertive initial model. Thus, these approaches are increasingly incorporated into parameter estimation studies and are also gaining more popularity in the geophysical community. In this study we present the Barnacles Mating Optimizer (BMO), a recently proposed global optimizer motivated by the special mating behavior of barnacles, to interpret magnetic anomalies. This is the first example in the literature of BMO application to a geophysical inverse problem. After performing modal analyses and parameter tuning processes, BMO has been tested on simulated magnetic anomalies generated from hypothetical models and subsequently applied to three real anomalies that are chromite deposit, uranium deposit and Mesozoic dike. A second moving average (SMA) scheme to eliminate regional anomalies from observed anomalies has been examined and certified. Post-inversion uncertainty assessment analyses have been also implemented to understand the reliability of the solutions achieved. Moreover, BMO’s solutions for convergence rate, stability, robustness and accuracy have been compared with the solutions of the commonly used standard Particle Swarm Optimization (sPSO) algorithm. The results have shown that the BMO algorithm can scan the model parameter space more extensively without affecting its ability to consistently approach the unique global minimum in this presented inverse problem. We, therefore, recommend the use of competitive BMO in model parameter estimation studies performed with other geophysical methods.

SUMMARY: A gravity inversion procedure using the success-history-based adaptive differential evolution (SHADE) algorithm is presented to reconstruct the 3-D basement relief geometry in sedimentary basins. We introduced exponential population size (number) reduction (EPSR) to reduce the computational cost and used self-adaptive control parameters to solve this highly nonlinear inverse problem. Model parametrization was carried out by discretizing the sedimentary cover via juxtaposed right prisms, each placed below each observation point. Resolvability characteristics of the 3-D inverse problem were revealed through some cost function topography landscapes. The fine-tuned control parameter namely, population number allowed us to get best benefit from the algorithm. Additionally, a stabilizing function as a relative constraint was used to avoid undesired effects originated from the ill-posedness of the problem. In the synthetic data cases, the strategy we propose outperformed the linear population number reduction strategy which has won various IEEE–CEC competitions so far. Thorough uncertainty assessments via probability density function and principal component analysis demonstrated the solidity of the obtained inverse models. In the real data case, residual gravity anomalies of two well-known major grabens of Aegean Graben System (Türkiye), calculated thanks to the finite element method, were inverted. It was determined that the inverse solutions obtained for these basement reliefs, whose depths are still controversial, are statistically reliable. Moreover, these depths were found to be less than the depths reported in most previous studies. We conclude that the SHADE using EPSR strategy that we propose is a powerful alternative inversion tool for highly nonlinear geophysical problems.

Enhancing the detection accuracy of the edges of the geological features within the subsurface remains a significant objective in geophysical data interpretation. Despite numerous advancements, approaches stemming from the directional gradients of gravity and magnetic fields still grapple with challenges such as low-resolution outcomes and susceptibility to noise contamination. In this study, we introduce a novel filtering framework based on the total horizontal gradient and its derivatives, designed to yield more precise and coherent edges free from false boundaries or disruptive artifacts. Validation using synthetic Bishop complex magnetic and gravity datasets, alongside Tuangiao aeromagnetic data from Vietnam, substantiates the robustness and applicability of our modified approach. Furthermore, recognizing the inherent susceptibility of edge detection filters to noise contamination resulting from directional derivatives, we employ the recently developed modified non-local means (MNLM) algorithm to alleviate noise effects prior to the analysis of noisy synthetic and real datasets. Our findings confirm the efficacy of the proposed method in reducing false artifacts and identifying edges with heightened precision, positioning MHGA as a valuable alternative for processing potential field data.

Geomagnetic anomaly interpretation through inversion procedures often yields useful results in determining the key details of ore masses. However, the problem is complicated due to the known ambiguous phenomena of the inversion process. Thus, details such as location, depth and shape can only be estimated using an efficient algorithm. To this end, we presented here a novel global optimization algorithm called Hunger Games Search (HGS) for the inversion of geomagnetic anomalies caused by ore masses. HGS is a well‐organized metaheuristic inspired by the hunger‐driven instincts and social behavioral decisions of animals. This is the first work in the literature to introduce this optimizer for the inversion of geophysical anomalies. We revealed the model parameter dependencies and the optimum control parameter values of the algorithm by performing some modal analyses and parameter tuning studies, respectively. The capability of HGS was demonstrated on synthetically produced geomagnetic responses and on three real anomalies obtained from some exploration fields in India and the USA. Uncertainty appraisal studies showed the solidity of the outputs of the algorithm. Moreover, some comparative studies with the fine‐tuned standard Particle Swarm Optimization, a commonly used tool for the inversion of geophysical anomalies, indicated that HGS algorithm provides better results in terms of convergence characteristics and solution stabilities in the presented inverse problem. We therefore recommend the use of this metaheuristic in model parameter estimation studies when dealing with this kind of geophysical potential field problems.

The self-potential method (SP) has been used extensively to reveal some model parameters of various ore deposits. However, estimating these parameters can be challenging due to the mathematical nature of the inversion process. To address this issue, we propose here a novel global optimizer called the Modified Barnacles Mating Optimizer (MBMO). We improved upon the original approach by incorporating a variable genital length strategy, a novel barnacle offspring evolving method, and an out-of-bounds correction approach. The MBMO has not been previously applied to geophysical anomalies. Prior to inversion of real data sets, modal and sensitivity Analyzes were conducted using a theoretical model with multiple sources. The Analyzes revealed that the problem is modal in nature, model parameters have varying levels of sensitivity, and an algorithm that can well balance global exploration with local exploitation is required to solve this problem. The MBMO was tested on theoretical SP anomalies and four real datasets from Türkiye, Canada, India, and Germany. Its performance was compared to the original version under equal conditions. Uncertainty determination studies were carried out to comprehend the reliability of the solutions obtained via both algorithms. The findings indicated clearly that the MBMO outperformed its original version in estimating the model parameters from SP anomalies. The modifications presented here improved its ability to search for the global minimum effectively. In addition to geophysical datasets, experiments with 11 challenging benchmark functions demonstrated the advantages of MBMO in optimization problems. Theoretical and field data applications showed that the proposed algorithm can be used effectively in model parameter estimations from SP anomalies of ore deposits with the help of total gradient anomalies.