## 6 Conclusions  
For the first time, a detailed field measured Bouguer gravity anomaly map has been obtained for the Eskisehir urban area and surrounding area. Linear features have been detected and mapped by the horizontal gradient magnitude (HGM) and Euler Deconvolution (ED) methods. It was observed that the locations of the subsurface structures obtained from both methods are quite compatible. It can clearly be seen in both methods that the general linearity characteristic of the Eskisehir Basin is in the northwest-southeast direction. This result was found to be compatible with the structure of the Eskisehir Fault Zone. It was inferred that the lineaments determined by edge detection, made with the methods used in the study, are the result of faulting. The most important indicator of this result is that the seismic activity in the region is compatible with these lineaments. It is suggested that the gravity data obtained from the edge detection and the lineaments in the northeast of the study area be investigated in detail. Lineaments within the Site A and Site B areas of Figure 8 were also found in the fault segments documented in the literature. Therefore, detailed research is recommended to reveal the potential dangers and risks arising from the faulting hazard in these regions. If the lineament maps obtained as a result of the HGM and ED are examined by superimposing them, a great deal of consistency is found (Figure 9). In particular, the HGM2 and ED5 lineaments and the HGM1 and ED3-4 lineaments show a high amount of consistency. Considering the differences between the HGM and ED, it has been observed that the ED solutions capture more linearity, since they also include vertical derivatives. However, both methods are extremely sensitive to noise in the data. In addition, deciding on the structural index and window size for the ED solutions is important for the reliability of the solutions. In other words, the HGM only offers more basic solutions; that is, it is crude compared to the ED and does not contain any information regarding the depth of the source. However, the choice of window size and structural index gains importance when applying the ED method. When the lineament solutions are examined together with the faults that fall within the study area, and are already known in the literature, a great deal of consistency is noticeable (Figure 10). While the edge detection solutions with the northern fault group show great consistency, the southern fault group with the HGM2 and ED5 lineaments seems to be somewhat offset towards the north, contrary to the faulting reported in the literature. The small-scale faults in the south are similar to the faulting in the literature. As a result of this study, most of the lineaments determined using the gravity method reveal buried faults revealed by previous studies. Only the ED1 and ED8 can be considered as newly revealed buried structures. Further investigation is needed to reveal the true nature of these lineaments, and to update the geological structure map of the area. However, it is clear that the edge detection using the land gravity survey and gravity data conducted for the first time in the region supports the known faults in the literature. The exact location of the southern group of faults proposed by other researchers was determined to be further north as a result of the gravity study. Most of the faults currently known in the literature have been identified through geological observations, and they are usually associated with surface ruptures. However, considering that the faults have different subduction angles, it is possible that they may have pushed the gravity anomaly further north than surface indications observed. It should also be mentioned that, based on our experience, it has been observed that the gravity method is much more advantageous in terms of both resource requirement and ease of implementation compared to other geological and geophysical methods previously applied in the field.