“…Gravity anomalies are frequently characterized by variations in rock density, and the inversion of this data is crucial for determining the mass properties and depths of the subsurface. Over the years, numerous geophysical methods for solving inverse problems have been developed, including robust simultaneous joint inversion [ 1 ], fair-function minimization [ 2 ], the s-curves method [ 3 ], non-linear least-squares [ 4 , 5 ], derivative-based approaches [ 6 , 7 ], moving average methods [ 8 , 9 ], multiple-linear regression [ 10 ], R-parameter imaging [ 11 ], and the Lanczos bidiagonalization method [ 12 ]. Metaheuristic optimization algorithms, such as genetic algorithms [ 13 ], particle swarm optimization [ [14] , [15] , [16] , [17] ], simulated annealing [ 18 ], cuckoo optimization algorithm [ 19 ], and ant colony algorithm [ 20 ], have also been used recently for gravity data inversion.…”