The Eastern Pontide orogenic belt is divided into three subzones ( northern, southern and axial zones) based on present lithologies and facies associations. NE-SW-, NW-SE- and E-W-trending fault systems, which play an important role in palaeotectonics and neotectonics of the Eastern Pontides, separate these zones. Three different methods were used to estimate the Moho depth from observed gravity value, namely (i) empirical relationship between Moho depth and Bouguer anomaly; (ii) spectral analysis of the radial wave number and; (iii) by the gravity inversion method. Power spectrum, a statistical approach, is a widely used technique to determine the depth of geological sources successfully. The empirical linear relations between the Bouguer anomalies and seismically determined crustal thicknesses have been used to compute depths to the Conrad and Moho discontinuities, which are consistent with the average depths obtained from the power spectrum method. The crustal structure was also determined from the inversion of a Bouguer anomaly profile along longitude 39 degrees. Gravity inversion results are consistent with those obtained from the empirical relations and power spectrum methods. We calculated the maximum crustal thickness of 43.8 km in the studied region by using the gravity inversion method, which also showed that crustal thickness increases from north to south.