Analysis of the gravity data, with the help of the seismic refraction results (G.D.rnMackenzile, H .Thybo and P .K .H .Maguire 2005) and Euler deaconvolution results alsornproduced to simplify the modeling process of the study region and it was done to studyrnthe crurtal thickness variation under the subsurface layers of Metehara.rnWhich is located 39° 301 48.34011 to 40° 591 27.60011 E longitude and 8° 481 15.12011 torn9° 28141.26811 N latitude as shown in Fig. I. I.The Bouguer gravity anomaly which showsrnthe area of study relatively higher value with a strongly denser contour at the lowerrncomer of S-E of the anomaly map.rnThe Gravity modeling constrained with seismic refraction and Euler deconvelutivn result,rnindicates a crurtal thickness of about 34.8 and 36.6km under selected profiles AB andrnMN respectively. The model also shows two highest dense bodies as mafic intrusionsrnover profile AB.