In this work, gravity anomalies have been analyzed using gradient analytic signal (GAS) obtained from the square root of the sum of the squares of the second complex and vertical gradients. The gravity anomalies have been decomposed at 1, 2 and 3 levels with Haar mother wavelet. The DWT leads to a decomposition of the approximation coefficients in four distinct components: the approximation, horizontal, vertical and diagonal. I have tested the maxima of the magnitude computed from the square root of the sum of the squares of the horizontal, vertical and diagonal components (HVDM), and maxima of GAS in imaging the source edges in theoretical examples, with and without random Gaussian noise. The basins and uplifts of western Anatolian have received more attention in geothermal potential research in terms of tectonic position. I utilize the GAS and DWT of Bouguer gravity data to image the geological discontinuities. These methods detected many faults, compared with the mapped faults in the surface geology. In addition, a 3D model of the basement undulation from residual gravity field is designed by means of Parker-Oldenbures algorithm to find the relationship between the lineaments from the maxima of GAS and HVDM. The algorithm reveals the presence of basement depths between 0.1 and 3.5 km in the region bounded by NW-SE and E-W trending lineaments. The lineaments from all methods and 3D basement model are well correlated. (C) 2014 Elsevier B.V. All rights reserved.