14
Archived Theses
0
DOIs Assigned
0%
DOI Rate
Discipline
A rapid estimation of moment magnitude mw for the October 23, 2011 Van earthquake using strong-motion records
It is well known that the Richter scale (local magnitude scale, ML) based on the measurement of the amplitude in a narrow- band time domain saturates for large crustal earthquakes, and the possibility of observing clipping on seismograms is high. Therefore it is insufficient to use the local magnitude for reliable magnitude estimation. A typical example of this situation is confronted in case of the October 23, 2011 Van earthquake ( ML 6.6 vs Mw 7.2) . Here Mw is the moment magnitude obtained from the far-field displacement spectrum of body waves. On the other hand, a better match is expected in the range of 3 < ML < 7 because in this range the empirical relation between seismic moment M0 and ML is almost the same with the relation between M0 and Mw (Thatcher and Hanks, 1973; Hanks and Kanamori, 1979). Hence, there is a need to find a way for the determination of Mw from near-field records too.The aim of this study is to determine moment magnitude of earthquakes very soon after occurrence, before any source or focal mechanism inversions have been performed. In this respect (Delouis et al.,2009) proposed the MWSYNTH method. In this method, observed displacement spectra are compared with synthetic spectra computed for a variety of finite dimension source models scaled with Mw. In this study we have tried to obtain a rapid estimation of Mw for the October 23, 2011 Van earthquake from the strong-motion records in the epicentral distances about 120km using MWSYNTH method. Using strong motion records of only two stations we obtained an agreeable estimation of moment magnitude Mw=6.9. The result may be improved (approach more seriously to the reference moment magnitude Mw ref=7.15) by doing computation for more stations. And it is foreseen that a correct recovery of Mw for the Van Earthquake within about 100 sec after origin time (recording+computing time) will be possible if the procedure is automated.
Analysis of the seismic activity in Tuz Gölü Fault Zone and surrounding areas following the 2023 East Anatolian earthquakes
The aim of this study is to examine whether the major earthquakes that occurred in the East Anatolia in 2023 had an impact on the Tuz Gölü Fault Zone and its surroundings and to carry out earthquake analysis of the region since the earthquakes. During the study, the seismic activity of the region in the 9-month period starting from the February 6 earthquakes was examined using Kandilli Observatory and AFAD data. Focal mechanism solutions and relocation studies were also carried out to obtain detailed results. Following the 2023 Eastern Anatolia earthquakes, significant seismic activity was detected in the south of the Tuz Gölü Fault Zone in Central Anatolia. In this study, earthquake relocations and focal mechanisms related to this activity are discussed. The relocation and focal mechanism solutions of the seismic activity of the region in the two-month period starting from the February 6 earthquakes shows that the activity occurred on two almost N-S structures that are to the west of Hasandağı Volcano. The mechanisms and epicenter locations are consistent with a right-lateral N-S fault for the eastern activity with some normal component. The activity might be related to the volcanism considering that they are on the proximity of fissures to the west of Hasandağı.
Bayesian modeling of co-seismic slip using GNSS data: 2002 Çay earthquake
In southwestern Turkey, in the Afyon-Sultandağı region, a mainshock of magnitude Mw 6.5 occurred at 07:11:28 (GMT) and was followed approximately two hours later at 09:26:49 (GMT) by an aftershock of magnitude Mw 5.8. The mainshock took place south of Eber Lake, causing surface discontinuities and significant surface deformation. Numerous researchers have studied the earthquake's source parameters using traditional teleseismic and geodetic data. However, to better visualize the co-seismic deformation and more accurately determine the geometry and slip distribution of the seismogenic fault, a Bayesian approach was adopted in this study. The Bayesian approach creates a posterior probability density function (PDF) that encompasses all possible values of model parameters by incorporating observations, prior physical knowledge of the source model, and the distribution of differences between forward model predictions and observations. In this study, geodetic datasets were utilized for both the rectangular source model and finite fault modeling. The results confirm that the 2002 Çay Earthquake predominantly involved a normal faulting mechanism with a significant left-lateral strike-slip component. The fault geometry, characterized by a strike angle of 277°, a dip angle of 43°, and a rake angle of -56°, aligns well with the regional tectonic framework. The uniform slip model, validated by the Bayesian method, effectively captured surface deformation patterns and achieved a variance reduction of 92%. This finding highlights the power of Bayesian methods in accurately modeling fault dynamics, particularly in regions with limited seismic data.
Determination of upper mantle heterogeneity beneath Aegean-Anatolian region from travel time tomography
The objective of this work is to determine the heterogeneities of the upper mantle in the Aegean-Anatolian domain using teleseismic tomography. A waveform dataset was prepared from 798 teleseismic earthquakes with magnitudes greater than 5.5 between January 2004 and December 2015. 417 stations from permanent and temporary networks with more than 64,000 direct P phases are used in the computations. The relative travel times of P waves with respect to the ak135 (Kennett et al. 1995) earth model are computed using waveform cross-correlations technique. The tomographic images are computed as perturbations with respect to ak135 earth model. An algorithm named as fast marching method (FMM) (Sethian, 1996a, 1996b) based on the solution of Eikonal equation is used in the forward computation of the travel times. The inversion is performed using subspace inversion scheme. Trade-off curves are plotted and several synthetic tests are performed in order to select optimum parameters (damping and smoothing) for tomography and the resolution and model roughness were investigated. The tomographic images obtained to a depth of 700 km. The computed tomographic images show a heterogeneous upper mantle structure in the Aegean-Anatolian domain. The results are similar to the previously published images mostly but provides higher resolution for the study area. Both Hellenic and Cyprus subductions are imaged to the depth of 700 km. The tear (Pliny-Strabo Tear) between two subduction zones is clearly observed reaching to 660 km discontinuity. A smaller scale tear (Antalya Bay Tear) is also observed on the Cyprus slab around Paphos Transform Fault. The Anatolian plate is underlined by low velocity mantle material with thickness increasing from west to east. The northern block of the North Anatolian Fault (NAF) is observed as high velocity body observable to a depth of 100-200 km. NAF has a sharp velocity contrast between the north and south.
Lithospheric structure of the Western Turkey and Aegean region
Aegean-Anatolia region undergoes an intense internal deformation as evidenced by the existence of major active faults, intense seismic activity and the marked thinning of the crust. It makes the region center of attraction to the study the interaction between the deep structure with the surface deformation. The aim of this study is to provide constraints on the crustal and uppermost mantle structure by using seismic data of permanent broad-band network of Kandilli Observatory and Earthquake Research Institute (KOERI-RETMC), and a temporary array of Seismic Imaging beneath Aegean-Anatolia Domain (SIMBAAD) experiment. Seismic stations of Republic of Turkey Prime Ministry Disaster and Emergency Management Presidency (AFAD), Incorporated Research Institutions for Seismology (IRIS) and previous experiment called Western Anatolia Seismic Recording Experiment (WASRE) were used to complement the network. In this regard we present two high resolution lithospheric images along a ~650 km transect crossing western Anatolia at 28°E longitude from the Black Sea to the Mediterranean and a ~550 km transect crossing central Anatolia at 30.5°E longitude. A total of 5250 receiver functions are computed from the records of teleseismic earthquakes at 40 broadband seismic stations for each of the profiles with an average spacing of ~ 15 km. Lateral variations of crustal thickness, Vp/Vs are inferred from both H-K, and common conversion point stacks (CCP). In order to have a better idea on the accuracy of the estimated crustal parameters we also performed a search scheme based on the Neigboorhood Algorithm. The receiver functions are inverted for a 1-D layered medium to determine the layer thicknesses, Vs and Vp/Vs. The CCP images reveals a long-wavelength variations of Moho depth from ~31 km in the Thrace basin to ~25 km beneath the Marmara Sea, ~25 km beneath the Menderes Massif and ~20 km on the coast of the Mediterranean on the western Anatolia transcent. On the eastern transect, a smooth Moho topography is observed with a sharp discontinuity at depths ranging from 34 km beneath the Black Sea coast, ~35 km beneath the Sakarya Zone with mafic composition to 43 km beneath the Antalya Bay on the central Anatolia profile. The Moho of the subducted African lithosphere is imprinted between ~40 and ~60 km depth at the southern end of the western Anatolia profile, dipping northward where the subducted Cyprus lithosphere is observed dipping northward with an angle of 40◦ between ~40 and ~100 km depths beneath the Antalya Bay on the central Anatolia transect.
Marmara Denizi'nde dinamik deprem kirilmasi simülasyonlari
The 1912 Mürefte and 1999 Izmit M_w7.4 earthquakes are the last devastating events of the western and eastern sections of the Marmara region, respectively. The center of the Sea of Marmara, the region between locations of these two earthquakes, is prone to creating another large earthquake. The main objective of our study is to determine 3D dynamic earthquake rupture scenarios, considering non-planar and heterogeneous stress distribution in the Sea of Marmara. Recent studies show that some segments of the North Anatolian Fault (NAF) beneath Marmara are partially creeping. In this study, it is the first time that we attempt to generate realistic earthquake scenarios by putting constrains on initial stress on the fault using regional stress from earthquake focal mechanisms, in addition to stress release during past earthquakes and strain accumulation during interseismic period using geodetical measurements on slip-rate and locking depth at various segments along the NAF beneath the Sea of Marmara. In order to constrain the regional stress in addition to our previous five cluster analysis a new earthquake cluster is analyzed in the Central Marmara Basin. We use 3D Finite Element Method (PyLith) for dynamic earthquake simulations and tetragonal mesh for better smoothing at the fault bends, which allows us to implement nonplanar fault geometry and initial stress heterogeneity using slip-weakening friction law. We place constraints on initial shear stress from geodetic and seismic studies of locking depth and interseismic strain accumulation. We consider 80 rupture scenarios and calculate slip distribution, rupture velocity and moment magnitude in addition to slip-rate and traction on the fault surface, and displacement and velocity on the ground surface. We find that for the most scenarios possible earthquake magnitude does not exceed Mw7.2. In addition, in none of the possible scenarios we obtain super-shear rupture velocity. We find that depending on the location of the initiation point, asperities in the partially creeping segments and loaded initial stress, the rupture may not extend into the Prince's Island Segment.
Studying seismotectonics of Eastern and Southern Anatolia using earthquake mechanisms
The Anatolia-Aegean domain provides a unique opportunity to explore plate interactions where oceanic subduction, continental collision and transform plate motions are observed simultaneously. High seismicity rates and diversity of the earthquake source mechanisms are the result of the accommodation of these relative plate motions. As the initial tectonic buildup involves the amalgamation of different tectonic units, it is natural that lithospheric segments with varying structural properties in this relatively small region also contributes to the complexities of the observations. Understanding interactions of these plates and related deformation requires an integrated analysis of various observations such as seismic tomography, earthquake slip models, geodetic observations and stress changes along with the seismicity and earthquake source mechanisms. In this thesis, 3 case studies in different tectonic settings are presented: the continental collision in the east, the extension due to roll back in the west and the transition between extension and compression. For these 3 case studies, the relation of earthquake source mechanisms to other seismological and geodetic data is used to better understand the present state of the seismotectonics of Easternmost Mediterranean including eastern Anatolia. The October 23, 2011 Mw7.1 Van, Eastern Anatolia earthquake which is on an EW trending thrust fault, in a region under N-S compression due to the convergence of the Arabian plate toward Eurasia. The three faults were activated during and after the coseismic rupture. The earthquake source mechanisms with consistent orientations are grouped in three clusters. An average fault mechanism is calculated for each cluster by the summation of moment tensors. The triggered faults have experienced Coulomb stress increase due to co-seismic rupture revealing a mechanism which accommodates NS shortening in the region. The June 20, 2017 Mw 6.6 Bodrum-Kos earthquake which occurred on an E-W trending normal fault is related to the roll back effect of Hellenic Subduction. The Bodrum-Kos event revealed that the extension in the western section of Gökova Bay is accommodated by a north dipping fault. Two different fault slip models, dipping to north and south, are used to compute the Coulomb stress changes at different depths. The coherency between the seismicity and the regions of increased stress is used to put a constraint on the dip of the ruptured fault. The gradual change of strikes of aftershock mechanisms from east to west is consistent with the rotation of the strain field region indicating that the observed earthquake pattern during the 2017 earthquake reflects the long term tectonic frame work in the region. In between these compressive tectonics of Eastern Anatolia and extension in the Aegean, Cyprus Arc region acts as a transitional zone which is tectonically less understood. Specifically how the convergence of Nubia toward Anatolia is accommodated remains unclear. By the analyses of novel earthquake source mechanisms, and other seismological and geodetic data, it is proposed that the segmentation of the subducting Nubian Plate has a significant contribution to the lithospheric deformation. The change in the orientations of the earthquake mechanisms around the Isparta Angle determines the eastern boundary of the N-S extension due to roll back of the Hellenic slab and is consistent with the counter clockwise rotation of AnatoliaAegean domain which is revealed by the recent GPS vector field. Thrust mechanism earthquakes along with Bouguer gravity, seismicity, and horizontal GPS velocities reveal the geometry of the subducting slab beneath Antalya Basin towards N-E. We suggest that the Antalya Slab deforms as an isolated block, responding in part to adjacent plates, including the Anatolian Plate that moves toward the west, overriding the remnant Antalya slab.
Interseismic behavior along the North Anatolian Fault in the Marmara Region using 3D structure
A series of earthquakes occurred along the North Anatolian Fault (NAF) during the 20th century, primarily migrating from east to west. The only part of the NAF that has not broken is under the Marmara Sea. The Main Marmara Fault (MMF), the NAF's northern branch, is the most active one, with the highest slip rate amongst the several branches of the NAF. Since the seismic gap of ~150 km is beneath the sea, the geodetic data is not sufficient to constrain the full fault coupling, particularly in the Central Marmara. Nevertheless, the current data does imply that the GNSS vectors along the northern coast of the Marmara Sea are smaller than expected. One interpretation is that the MMF has heterogeneous interseismic coupling with creeping and locked segments. Another explanation is that the fault is locked, but the strain is asymmetrically localized around the MMF as a result of the deep basins. In this study, the competing effects of weak interseismic locking of the MMF and deep basins around the fault are studied by developing a 3-D finite element model for the Marmara Region, which includes a realistic topography, the 3-D geometry of the main fault, and basins, and using the geodetic data as a constraint. Our findings show that the deep basins confine the interseismic strain in the fault vicinity, and using a homogeneous half-space model leads to a slight underestimation of the locking depth. Our 3-D model shows that while the basins have some effects on strain localization, the heterogeneity of interseismic coupling is necessary to explain the observed GNSS data. We infer a change in the locking depth at the Ganos Bend between the strongly coupled Ganos and the weakly coupled Western Marmara. Seismic studies also indicate that these two segments vary considerably in background seismicity. The 50 km creeping segment coincides well with repeating earthquakes and higher rates of diffuse seismicity. Variations in regional stresses and earthquake focal mechanisms, including the 2019 Silivri earthquake sequence, are compatible with the dilatational quadrants in the region due to the loading caused by the interseismic creep of the Western Marmara.
Simplektik geometri ve hamiltonian Monte Carlo metodu
Hamiltonian Monte Carlo (HMC) method is an application of a non-Euclidean geometry to an inverse problem. HMC is a probabilistic sampling method with the basis of Hamiltonian dynamics. One of the main advantages of HMC algorithm is to draw independent samples from the model space with a higher acceptance rate than other Markov Chain Monte Carlo (MCMC) methods. In order to understand how higher acceptance rate is achieved, I have studied HMC in the light of symplectic geometry. Hamiltonian dynamics is defined on the phase space (cotangent bundle), which has a natural symplectic structure, i.e. a differential two-form which is non-degenerate and closed. Hamiltonian function is defined on the phase space, which corresponds to the sum of misfit and the square of the generalized momentum. By using the non-degeneracy property of symplectic form, a vector field can be found in which Hamiltonian function is invariant along the integral curves of the vector field. The invariance of the Hamiltonian function results in high acceptance rate, where we apply accept-reject test to satisfy detailed-balance property. In this thesis, we define some basic concepts and theorems in symplectic geometry, then describe the relation between symplectic geometry and HMC, namely Hamiltonian dynamics. Lastly, we show an implementation for HMC algorithm to a 2D-tomography problem and analyze the tune parameters for application of HMC.
High resolution microseismicity and nearly-repeating events in the Marmara Sea
The Main Marmara Fault beneath the Marmara Sea has a prominent seismic gap that can produce a devastating earthquake and a serious risk for the surroundings. It is important to scrutinize the seismic activity in region and relate this activity to the deformation of the fault zone. In this study, a new micro-earthquake database is created for the Marmara Sea between 2014-2016 using the data mostly from ocean-bottom seismometers. The detected and located seismicity indicate that Tekirdağ Basin hosts a diffuse activity from ~7 km to about 18 km depth. A high micro-earthquake activity rate predominates beneath the Central Basin, at depths from 3 km to 15 km. The abundancy of earthquakes in the area can be attributed to a creeping zone, considering the conformity with the geodetic observations. On the other hand, Kumburgaz and the western part of Çınarcık Basins show sparse seismicity at depth ranges of 5-19 km and 3-18 km, respectively, signing to a locked fault compatible with the geodetic observations. In addition to micro-seismicity, the repeating events are detected using template matching method on the continuous waveforms from 2008-2021. The clusters of highly correlated detected earthquakes, which are closely spaced or partially overlapped, are attributed to the "near-repeating earthquakes". The nine nearly-repeating earthquake clusters beneath the Central Basin are observed at 8-13 km depths, suggesting seismic creep behavior together with a high seismicity rate. The fault mechanisms of the near-repeater clusters have strike-slip mechanism consistent with Main Marmara Fault zone. The nearly-repeating events have two different patterns of repeating intervals, as long-term and short-term type events. The amount of slip rates from the near-repeater clusters shows varying slip rates but comparable to geodetic rate. The number of near- repeating events decreased significantly after the 2018 and no repeating event is observed during 2019 which Mw 5.8 Silivri earthquake occurred.
Electrical resistivity structure beneath mt. Erciyes Central Anatolia, Turkey
The Central Anatolia (CA) is an intermediate belt between the Pontides and the Anatolide-Tauride blocks. Numerous basins, folds, faults, volcanic fields, and rich mineral resources in the region have attracted the attention of geoscientists for years. Mt. Erciyes is a Quaternary-aged stratovolcano located in CA, which has not been activated for a long time. In this thesis, I aimed to perform three-dimensional (3 - D) imaging of its electrical properties using magnetotellurics (MT) and interpret its structure. MT method has been extensively used in the exploration of volcanic areas. This study can be considered as the first systematic application of the MT performed in the vicinity of Mt. Erciyes. Between 2013 and 2018, deploying 48 wide-band MT soundings forming a grid surrounding the volcano and its location as part of the Continental Dynamics/Central Anatolian Tectonics project (CD/CAT). The obtained data set was used for determining the dimensionality of the survey area, and the realistic electrical structure of the region was achieved with 3 - D solution algorithms including 38 station data. The results showed that the regional geoelectric strike direction is ∼ N25◦E which coincides with the structural strike of the region. For determining the dimensionality of the MT data, phase tensor analyses were performed. The observed wide-band data were inverted using a 3 - D inversion algorithm known as the ModEM © in two distinct ways; with and without topography information to reveal the electrical structure. The 3 - D inversion results revealed three distinct conductivity features: a widespread conductive anomaly found beneath the volcano and two conductive areas in the vicinity of the volcano. When the final model is examined, it is possible to observe that the conductive regions are separated from each other. If the final model is examined together with the geological data, it is noticed that the boundaries of the conductor coincide with the faulting in the region.
Inversion for elasticity tensor of focal region using machine learning algorithms
The moment tensor is an essential tool in seismology to examine the structure of seismic sources. The deformation at the focal region, using combinations of force couples arranged in a 3 x 3 matrix, is represented by the moment tensor. The moment tensor can be expressed as a linear combination of the eigenvectors of the anisotropic focal region's elasticity tensor. The eigenvalues of a vertically transversely isotropic (VTI) elasticity tensor from an occurring moment tensor of a focal region can be obtained, and the precision of this determination depends on the degree of anisotropy. Machine learning optimization involves iteratively enhancing a machine learning model's precision by reducing the error level. Choosing an algorithm that can effectively sample the search space and identify optimal solutions is necessary to optimise a function. Many algorithms are available for function optimization, but it is crucial to set a baseline to determine the practicable solutions for a given problem. This thesis defines a new objective function (misfit function). The function is proposed for obtaining the elastic parameters of an anisotropic focal region, and these parameters are calculated by using machine learning algorithms such as Grid Search, Random Search, Simulated Annealing, and Nelder-Mead.
Potansiyel alanlarda otomatik değerlendirme teknikleri ve arkeolojik alanlara uygulanması
Bu çalışma, çoğunlukla otomatik veya yan otomatik değerlendirme yöntemleri olarak isimlendirilen jeofizik değerlendirme tekniklerinin (Normalize Edilmiş Tam Gradyan (NTG), Euler Dekonvolüsyon Yöntemi (EDY)) ve bu yöntemlere ek olarak Yatay Gradyantm maksimum noktalarının yerlerinin araştırılması yönteminin (Sınır Analiz Yöntemi (SAY)) arkeolojik alanlarda uygulanmasını içerir. İlk iki yöntem kaynak manyetik yönüne ait ön bilgi gerektirmemektedir. Sonuç olarak, bu yöntemler jeolojik bilginin az olduğu alanlarda başarılı şekilde uygulanabilirler. NTG yönteminde, manyetik veriler için Fourier katsayıları filon yöntemiyle hesaplanmış ve türevlerinden hareket ederek NTG değerleri hesaplanmıştır. Basit modellerin etkisi, NTG değerlerinin aşağı doğru analitik uzanımından elde edilen NTG kesitleri üzerinde araştırılmıştır. Aşağı doğru analitik uzanıma dayanan bu yöntemin temeli, uzanım değerinin doğru hesaplanması değil, tekil noktayla ilişkisini gösteren NTG alan karakteristiğinin elde edilmesi olduğuna değinmek önemlidir. Euler homojenite eşitliği, basit kaynaklar üzerinde manyetik anomalilerin değerlendirilmesi için, yeni bir teknik geliştirilmesinde kullanılmıştır. Euler homojenite ilişkisine dayanan EDY, manyetik verilerin otomatik olarak hızlı değerlendirilmesinde yaygın olarak kullanılır. Derinlik değerinin hesaplama kalitesi çoğunlukla manyetik kaynağın doğasıyla ilgili anomali azalma oranım tanımlayan ve yapısal indeks olarak bilinen parametreye bağlıdır.m Bu yöntemlere ek olarak, yapı sınırlarının belirlenmesi amacıyla, yatay gradyan genliklerinin maksimum noktaların yerlerinin harita düzlemi üzerindeki dağılımından yararlanılmıştır. Yukarıda değinilen bu yöntemler önce kuramsal modeller üzerinde sınanmış, daha sonra da Çorum-Ortaköy Şapinuva Hitit arkeolojik kentinden elde edilen manyetik verilere uygulanmış ve başarılı sonuçlar elde edilmiştir. Anahtar sözcükler : Arkeomanyetik Veri, Normalize Edilmiş Tam Gradyan, Tekil Nokta, Euler Dekonvolüsyon, Yapısal İndeks, Yatay Gradyan, Maksimum Nokta Yerinin Araştırılması
Heyelan duraylılığının jeolojik ve jeofizik yöntemlerle değerlendirilmesi (Gümüşhane Süleymaniye örneği)
Ülkemizde jeolojik, iklimsel, coğrafik özellikleri ve yanlış arazi kullanımı nedeniyle heyelan olayları sıkça yaşanmakta ve çoğu kez tekrarlanarak afete dönüşmektedir. Heyelanlar; akma, kayma, düşme ve devrilme gibi farklı hareket mekanizmalarda gerçekleşmekte olup depremler, ani-şiddetli yağışlar ve sıcaklık değişiklikleri ile tetiklenerek tehlikeli olabilecek afet olaylarına yol açmaktadır. Bu nedenle uzun süreli ve bol yağışın etkili olduğu, nemli iklime sahip Doğu Karadeniz Bölgemizde heyelanların afete dönüşerek büyük ölçüde can ve mal kaybına yol açmasını önleyici mühendislik çalışmalarına ihtiyaç duyulmaktadır. Kayan kütlenin ve altındaki yapının özelliklerinin yüksek doğruluk oranıyla belirlenmesi heyelan yapısının ortaya çıkarılmasında çok önem arz etmektedir. Bu çalışmada, Gümüşhane ili Süleymaniye Kayak tesis alanında jeolojik ve jeofizik yöntemler kullanılıp heyelan duraylılık analizleri yapılarak tehlike oluşturabilecek heyelan alanlarının ortaya çıkarılması amaçlanmıştır. Bu kapsamda, jeofizik yöntemlerle heyelanın sınırları, kayma yüzeyinin derinliği, heyelan kütlesinin farklı kayaç kesimleri, yeraltı suyu durumu, temel kayanın derinlik, yapı ve bileşim özellikleri tespit edilerek sonlu elemanlar ve limit denge analizleri ile düşey ve yanal deformasyonlar ortaya koyulacaktır. Bu değerlendirmeler sonucunda riskli alanlar tespit edilip en uygun iyileştirme yöntemleri önerilecektir.