A differential search algorithm (DSA) application, which is a metaheuristic inspired by nature, for total field aeromagnetic data caused by volcanoes over a 2D dipping dyke is presented. Inversion of the total magnetic anomalies was performed by adding the background level in addition to the parameters of the dyke model (e.g., dip angle, the depth to the top, half-width, the distance from the origin to the reference point, and amplitude coefficient), which are often tried to be estimated in the literature studies. In synthetic dyke models, the efficiency of the DSA in parameter estimation of theoretically generated magnetic anomalies that do not contain noise and contain random noise at different levels has been demonstrated. Firstly, in the synthetic dyke model, the efficiency of the DSA in parameter estimation of theoretically generated noise-free magnetic anomaly is demonstrated. Additionally, different levels of random noise were added to the same synthetic model anomaly to test the performance of the algorithm in case the data contained noise. The results of the inversion show that the model parameters estimated from the DSA agree well with the correct ones. This fit was also statistically checked by calculating the probability density function. In the real case, the inversion approach was then used to interpret five prominent total aeromagnetic anomalies over the well-known Kula volcanic field located in western Turkiye. The depths and widths of these magmatic bodies lying underneath these volcanic cones are about 450 m and 470 m, respectively.