An assessment of seismicity and near surface geophysical characteristics of potential solid waste landfill sites in the Eastern Black Sea Region of Türkiye

This study aimed to assess the suitability of the potential solid waste landfill sites in seven provinces (Samsun, Ordu, Giresun, Trabzon, Gümüşhane, Bayburt and Artvin) in the Eastern Blacksea Region of Türkiye. The earthquake hazard analysis for two major earthquakes which occurred in the region was first carried out. Then, the geophysical methods including seismic refraction tomography (SRT), electrical resistivity tomography (ERT) and Multichannel Analysis of Surface Waves (MASW) were conducted to find out the structural and physical properties of the subsurface which include the layering, soil classification based on VS30 and the groundwater content at 25 locations of 13 in target provinces. The integrated interpretation of whole data sets demonstrates that Işıktepe, Esence, Çamburnu and Kazantaş which are characterized by VP > 1200 m/s, VS30 ≥ 400 m/s, ρ > 70 Ohm-m, low earthquake hazard and seismicity are more suitable among others. Vezirköprü, Şebinkarahisar, Yenice, Bayburt-Center, Balkaynak and Murgul will be suitable after a geotechnical reclamation due to moderate seismic velocities and electrical resistivity which are 900 < VP ≤ 1200 m/s, 200 < VS30 < 400 m/s and 10 < ρ ≤ 70 Ohm-m representing stiff and wet soils. In addition, Bafra, Ağalık and Ovacık were considered to be unsuitable due to the presence of thick, water-saturated soft soil and extremely weathered rocks. Finally, this study shows that the joint interpretation of seismicity and geophysical data in potential waste landfill sites, extremely important for the planning and development of a city, can provide the valuable information which will enable to prevent possible deformations, environmental problems and economic losses after waste landfill.


Introduction
In many countries, due to the rapidly increasing population and industrial facilities, uncontrolled growth and expansion of new urban areas, the production of solid waste (domestic and/or industrial waste) brings problems with its landfill and disposal.In recent years, pollution of the underground water resources and the environment, especially by leachate from old and existing irregular and uncontrolled landfills, landslides and collapses in the dumping areas have posed a high risk to public health (EPA 1999;Özel et al. 2017).On the other hand, with the development of new-generation construction and manufacturing technologies, sanitary landfills and disposal sites are now built with appropriate procedures according to engineering standards and then they are constantly monitored and inspected by local governments and environmental agencies.However, inadequacies in the construction and operation of many landfills still remain relevant in small and medium-sized cities, particularly regarding existing but deactivated landfills (Soupios et al. 2007).The fact that most of the waste landfill areas are close to settlements, markets, farms, highways, industrial areas, etc. has revealed the necessity of carefully examining not only the environmental conditions but also the soil and rock characteristics of the selected places before landfilling.Therefore, the demand for examining the suitability of the sites considered for future landfill planning in terms of environmental and structural characteristics has increased significantly in the last few years due to the lack of space in dense settlements (Rao 1997;Dorhofer and Siebert 1998;Dorn and Tantiwanit 2001;Mondelli et al. 2012).For this reason, the problems of environmental pollution and waste management have become a major area of interest for geoscientists and researchers from other related sciences.
Currently, in international and national regulations, preliminary studies in terms of city planning, geological, hydrogeological, geophysical, environmental science, social and health sciences and geographic information systems have become mandatory for the selection, design and construction of suitable solid waste landfill areas (Choudhury and Savoikar 2009).In addition to knowing the local geological characteristics of potential waste landfill areas, it is necessary to know their geophysical characteristics in detail.On the other hand, alternative areas for solid waste landfill should be evaluated specially in terms of seismology, their distance to active faults, historical and instrumental seismicity and the earthquake characteristics (magnitudes and earthquake recurrence times) (Krinitzsky et al. 1997;Yeşilnaçar and Çetin 2005).Moreover, earthquake hazard and risk of the region should be explained according to scientific methods.
Although the mechanical methods (standard penetration, cone penetration, pressure metre tests, etc.) used for subsurface characterization of an area provide reliable information about soil thickness and type, these tests are costly (money, labour and time) and destructive (Kosugi et al. 2006).However, the main disadvantage of mechanical methods is that they provide information for subsurface soil characteristics from only one investigation point and therefore require many exploration drillings to determine complex geological structures or lateral facies changes.The main way to overcome this problem is to carry out geophysical investigations.Because geophysical methods help reduce the number of these traditional and more expensive mechanical investigations, they make a significant contribution to directly reducing project costs (Zhan et al. 2018;Boudreault et al. 2010;Fernández-Baniela et al. 2021).
Geophysical methods, which allow rapid exploration of large areas without changing their natural conditions at each site in depth and spatial dimensions, allow the non-invasive determination of physical properties and mechanical behaviour of the soils and rocks such as softness-stiffness characterization of the geological units, soil classification, possible structural elements (faults, fractures, voids, caves, cracks, etc.), geometry of the soil (soil thickness and depth to the bedrock), dynamic-elastic parameters and groundwater content.Accordingly, geophysical engineering studies have recently become an important area of use in the waste management workflow, and the use of geophysical methods has become widespread (Otto and Sass 2006;Schrott and Sass 2008;Kowalczyk et al. 2017).In this context, magnetic, electrical resistivity, electromagnetic, natural potential, ground radar penetrating and seismic (refraction, reflection and surface wave analysis) methods are frequently applied.Thus, the geophysical methods include providing the information about the subsurface distribution of certain physical parameters (density, seismic velocity, electrical resistivity and conductivity, magnetic susceptibility, dielectric permeability and electromagnetic velocity, etc.).Furthermore, the combined use of geophysical methods provides important opportunities for the evaluation and characterization of the pollution caused by solid wastes (Kayabalı 1996;Önal et al. 2013;Konstantaki et al. 2015;Çınar et al. 2016).For this reason, possible soil problems such as sliding, collapse, breaking and cracking in waste landfill areas, mixing of waste leachate with groundwater and polluting water resources, damage to the natural built environment and related health and socio-economic damages will be prevented (Tchobanoglous and Theisen 1993;Pomposiello et al. 2012).
Although geophysical methods have a wide range of applications, their use in detecting and monitoring of the waste and solving the problems arising from landfill areas is methodologically limited.Moreover, when the literature sources are examined, the reasons for using and applying geophysical methods for solid waste landfills can be grouped into two main groups: 1.In the past and up to the present, the first group of studies has mostly involved in monitoring the existing landfills, types of the old waste and their distribution, delimitating the waste dump area, detecting and solving the problems (leakage, collapse, outgassing, etc.) by using geophysical methods (Siracusa et al. 2005;Kaya et al. 2007;Balia and Littarru 2010;Genelle et al. 2012;Özel 2018;Dumont et al. 2018;Neyamadpour 2019).2. The second group of studies was the study that led to the final decision about the suitability of the potential and/ or alternative waste landfill area.They include detailed examination of seismicity of the region, flood, landslide hazards-risks and soil characteristics after preliminary studies which use the remote sensing methods for the selection of landfill area location (Reynolds and Taylor 1996;Silvestri and Omri 2008;Mondelli et al. 2012;Pomposiello et al. 2012;Soupios and Ntarlagiannis 2017;El-Kelani and Khader 2019).Thus, these studies, which include near surface geophysical methods, seismic hazard and risk assessments, have recently started to be preferred in the workflow where site selection decisions are made.
In this study, within the framework of previously mentioned second group of studies, we firstly overviewed the seismicity of the region and then we indicated an integration of geophysical (SRT, MASW, ERT) methods in 25 profiles to characterize subsoil materials and their potentially risks for waste landfill site assessments for seven provinces (Samsun, Ordu, Giresun, Trabzon, Gümüşhane, Bayburt and Artvin) in the Eastern Black Sea of Türkiye.This is demonstrated in terms of investigating the subsurface materials within the upper, approximately 30-m-thick layer, examining the near surface/ bedrock structure under landfills such as near-surface bedding, softness-stiffness of soils and low-high resistive zones and evaluating the suitability of the studied areas by geophysical data.

Study area and surveying locations
The study area includes provinces (Samsun, Ordu, Giresun, Trabzon, Gümüşhane, Bayburt and Artvin) in the area of responsibility of the Eastern Black Sea Project Regional Development Administration (DOKAP), which covers an area of approximately 49,000 km 2 in the Eastern Black Sea Region and operates to accelerate the development of the region (Fig. 1).The total population of these provinces is approximately 4,000,000 people.According to Yıldırım et al. (2022), in the DOKAP (2021) report, it is stated that, as of 2018, the largest amount of waste was dumped in Samsun (~ 402,000 tons), and the smallest amount of them was disposed in Bayburt (~ 16,000 tons).The province with the largest average amount of municipal waste was Rize (1.25 kg/day/person), and the province with the smallest average amount of municipal waste was Trabzon (0.81 kg/ day/person) (DOKAP 2021).
However, according to preliminary studies (environmental and sociological criteria in accordance with national and international legislation) which include the opinions of institutions, organizations and external stakeholders in Rize Municipality, no potential suitable landfill area could be determined.Therefore, seismicity and near-surface geophysical studies were not carried out for this municipal.

Regional geology
The Eastern Black Sea Region is located in the northeastern part of Türkiye which consists a mountain chain approximately 500 km long and 200 km wide parallel to the southeast coast of the Black Sea and is located on the Eastern Pontide Orogenic Belt tectonic unit (Ketin 1976).As shown in Fig. 2, this tectonic unit is divided into three different subunits depending on lithological differences, geological and geophysical features, tectonic structures and lithological changes from north to south: the north zone, south zone and axis zone (Bektaş et al. 1995;Eyüboğlu et al. 2011).The paleo-tectonic evolution of the Eastern Pontide Orogenic Belt was controlled Fig. 1 The provinces where landfill areas will be determined and the selected geophysical (SRT, MASW, ERT) measurement locations (yellow pins) by faults of three different directions.Because the small and large blocks (north, south and axis zone) formed by the faults move independently and relatively, the geological features of each block are different (Eyüboğlu et al. 2011).Folds that develop because of the horizontal and vertical movements of the blocks have the characteristics of drag or drape folds parallel or almost parallel to the block edges.While Mesozoic and Cenozoic aged volcanic rocks and granitic intrusions are generally dominant in the northern zone, sedimentary rock series constitute the dominant lithology in the Southern Zone.Güven (1993), who conducted one of the most detailed studies in the region, examined the succession in the northern zone from the bottom to the top and named all units as the formation and lithology.According to Güven (1993), volcano-sedimentary series and volcanic and intrusive rocks containing the products of magmatism that continued their development from the Early Jurassic until the end of the Eocene are common in the region where rock units developed in the Paleozoic-Quaternary period are exposed.Moreover, sedimentary successions accumulated during periods when magmatic activity stopped.
Türkiye is an active and complex region in which the Eurasian, Arabian and Anatolian plates meet (McKenzie 1972;Alptekin 1973;Ketin 1977).Consequently, the Anatolian plate moves westward between the dextral strike-slip North Anatolian Fault zone (NAFZ) and the left-lateral strike-slip Eastern Anatolian Fault zone (EAFZ).In this system, as shown in Fig. 3a, the region is under the influence of a thrust fault (inclined to the south) parallel to the coast in the north (Okay and Sahintürk 1997;Nikishin et al. 2003) and NE-SW extension strike-slip faults called Rize, Trabzon and Ordu in marine side, which have been determined based on geological, geophysical and GPS data collected for the last 30 years (Eyüboğlu et al. 2011;Meisner et al. 2009;Keskin et al. 2010).

Seismicity studies
Historical (Ambraseys and Jacson 1981;Gündoğdu and Altınok 1986;Ayhan et al. 1987) and instrumental data on earthquake occurred in the region that may affect potential waste landfill sites, national (Boğaziçi University Kandilli Observatory-KOERI) and international (International Seismological Center-ISC) were provided by the catalogues published by the US Geological Survey-USGS data centres.Different magnitudes (M b , M L , M d , M w ) in the catalogues of the instrumental period between 1900 and 2016 were converted to surface wave magnitude (M s ) by using the transformation equations for Türkiye and its surroundings given by Aydın (2016), and the catalogue was homogenized.In addition, the intensity scale describing historical earthquakes was transformed into a magnitude scale using the Eq.(1), developed by Sayıl (2014).
Here, I o , σ 1 = ± 0.02 and σ 2 = ± 0.24 indicate the earthquake intensity and errors in the least square linear regression approach, respectively.Thus, the (1) M s = 0.47 ± 1 I o + 2.05 ± 2 Fig. 3 Main faults and seismicity of Eastern Black Sea Region.a Faults in the land and offshore threatening the region (adapted from Eyüboğlu et al. 2011).Red thick line in (A) approximately indicates the location of the cross section in (B).b Epicentral distributions of earthquakes with M s ≥ 3.0  in the region.Faults are digitized by using Emre et al. (2013) magnitude-epicentre distribution for the lower limit M s ≥ 4.0 of the data set obtained by applying the completeness analysis of the earthquake data is shown in Fig. 3b.Most of the earthquake locations shown in Fig. 3b occurred in the NAFZ.However, the offshore earthquakes are thought to be related to the coastal parallel thrust fault shown in Fig. 3a and the strike-slip Rize, Trabzon and Ordu faults.Readers can find detailed information from Eyüboğlu et al. (2011).The total length of the North Anatolian Fault Zone (NAFZ), which is a righthand strike-slip fault system, is approximately 1500 km.While the eastern part of the fault was under compression, the western part extended.The focal mechanisms of earthquakes occurring so far along the NAFZ confirm these different stress regimes (Sayıl 2014).The major historical (M s ≥ 7.0) earthquakes along the NAFZ are listed in Table 1.
In the NAFZ, the 1939 (M s = 7.9) and 1992 (M s = 6.6) earthquakes and the 1942 (M s = 7.0) earthquake occurred in Erzincan causing significant loss of life and property by affecting seven municipalities covering the study area.
The NAFZ is a fault zone with a high seismicity.Therefore, the seismicity of the region was examined for the region and each province.The annual number of earthquakes with M s ≥ 4.0 in the study area (histogram) is shown in Fig. 4a and the cumulative number of earthquakes is shown in Fig. 4b.Accordingly, 265 earthquake data were used, and the magnitude range was taken as ∆m = 0.1 to evaluate the Gutenberg-Richter magnitude-frequency relationship of the region.
For an earthquake with a magnitude of M s ≥ 4.0 which may occur in the study area, the Poisson model was used to calculate the probability of an earthquake occurring in the past or in the future.The probability of the occurrence of  "n" earthquakes with magnitude M ≥ M 0 at time t in a region studied is given as: where ν is the average number of earthquakes per time unit (usually one year) in the studied region, and "N" is a variable showing the number of earthquakes in the region under study at time "t".Magnitude-frequency relations are the basis of earthquake statistics and provide valuable information about the earthquake activity.A relationship between occurrence frequencies and magnitudes of the earthquakes is generally formulated by Eq. ( 3) (Gutenberg and Richter 1954).
Here, N is the cumulative number of earthquakes, and M is the magnitude.Parameters a and b are the constants of Eq. ( 3), while parameter a depends on the observation period, the width of the study area and the level of earthquake activity the b value, which is the slope of the linear relationship, provides important information about earthquake formation physics and tectonic structure.Therefore, the b value is used in earthquake prediction processes and the preparation of seismotectonic zoning maps (Wiemer and Wyss 1997;Papazachos 1999;Öztürk 2012).
Thus, when Eq. ( 3) is obtained for a region, the probability that any earthquake of magnitude M will occur in ( 2) that region for an observation interval of T years is given by Eq. ( 4): and the recurrence period is calculated by Eq. ( 5): For this study, the regional magnitude-frequency relationship was obtained as follows.However, the magnitudefrequency relations for seven provinces are given in Table 2. Here, N(M) represents the cumulative number of earthquakes.The parameters, a = 5.3, b = 0.76 and R 2 are the correlation coefficient of the regression.
Assuming the earthquakes occurring between 1900 and 2016 for the study region fit the Poisson distribution, probability calculations for earthquakes of various magnitudes were made by using Eq. ( 6).Accordingly, the probability of exceeding the earthquakes in every 10-year in the next 100 years, R(M) and recurrence periods (Q) were calculated and are given in Table 3.In addition, the probability of occurrence of earthquakes (hazards) in 10-year periods is shown graphically in Fig. 5.
According to the results, the probability of mediumsized (M s ≥ 5.0) earthquakes within 10 years is 98%, and the probability of strong (M s ≥ 6.0) earthquakes within ( 4)  10 years is 54% for the region in general.In addition, the probability of a major earthquake (M s ≥ 7.0) was calculated as 14% for 10 years and 54% for 50 years.The recurrence periods for moderate (M s ≥ 5.0) and strong (M s ≥ 6.0) earthquakes were determined as 2.6 to 13 years, and 74.6 years for a large earthquake (M s ≥ 7.0).For the two major earthquakes (Erzincan, M s = 7.9 and Tokat, M s = 7.0) in the NAFZ that affected the region, the most and caused significant losses, a deterministic earthquake hazard assessment covering the study area was performed.Within this framework, possible earthquake intensity and acceleration values were calculated.To decrease the soil-independent intensity and horizontal acceleration representing the study areas, the following relations were used: where I means MSK (Medvedev-Sponheur-Karnig) intensity scale (macroseismic intensity scale), M s means surface wave magnitude, D is the closest (or perpendicular) (km) distance to the fault, and PA is the maximum horizontal acceleration (gal = cm/s 2 ).
The intensity and acceleration values calculated according to Eqs. ( 7) and ( 8) for both considered earthquakes are summarized in Table 4.
According to the data in Table 4, the locations where the intensity and acceleration are likely to be the highest for the Erzincan earthquake are Bayburt and Gümüşhane (I = VIII, PA = 0.21-0.56g), while the locations where they are likely to be the lowest are Samsun and Artvin (I = V, PA = 0.09-0.16g).On the other hand, for the Tokat earthquake, the locations where the intensity and acceleration are likely to be the highest and lowest were Samsun (I = VII, PA = 0.19-0.40g), and Trabzon, Gümüşhane,

Near surface geophysical methods, data acquisition and processing
The geophysical measurements included seismic refraction tomography (SRT), Multichannel analysis of surface wave (MASW) and electrical resistivity tomography (ERT) measurements.They were preferred and used because of their advantages of rapid data acquisition and low cost and providing the useful physical parameters and images to characterize the subsurface in detail.The surveying profiles were decided after the studies of geographic information systems and multi-criteria decision analyses including geomorphology, vegetation, distance to subsurface-surface water resources (drinking water and streams), odour effects, climatic conditions, agricultural usage, distance to settlement areas and ease of transportation.These methods allow obtaining P-wave, S-wave velocities (V P , V S ) and electrical resistivity (ρ) parameters with high accuracy, which are necessary for qualified subsurface characterization.For this purpose, 25 SRT, MASW and ERT datasets were acquired from the same profile.For seismic measurements, a 12-channel seismograph, 4.5-Hz vertical geophones (P-geophone) and 9.0-kg sledgehammer (20 × 20 × 5-cm steel plate) were used as seismic source.In each spread, the first receiver offset in the end-on and -off shots was 12 m, and the distance between the receivers was 3 m.For the SRT, five shots (source locations, 0, 19.5, 28.5, 37.5 and 57 m) were performed along the profile.For SRT, the time sampling interval was 0.25 ms, and the recording length was 0.5 s.For the MASW data, the sampling interval was 0.5 ms and the recording length was 1.0 s.In order to improve the signal/noise ratio of the data at each source point, three successive shots were stacked vertically (Table 5).
SRT provided a 2D P-wave velocity models of the subsurface in the sites, which are used to detect the bedrock depth   6 the refraction data of profile 1 in Bafra site.Briefly, the seismic data in Fig. 6b were acquisited, then first arrival times of the critically refracted waves (black line in Fig. 6a) and direct waves (blue line in Fig. 6a) are manually or automatically picked from each seismic data, and distancetime graph is drawn as in Fig. 6c.Lastly, 2D P-wave velocity sections (Fig. 6d) are obtained by tomographic inversion using the nonlinear least square algorithm.In Fig. 6d, the blue coloured area in the lower right corner of the SRT section does not actually show a high velocity compared to the overlying units, and this has been taken into account in the interpretation of the sections.Furthermore, the absence of critical refracted in some layers relative to the ray paths is entirely due to the lack of sufficient velocity difference.MASW method, which consists of the data acquisition, dispersion analysis and inversion steps, was used to estimate the depth profile of shear wave velocity (1D Vs) in order to determine the soil thickness, soil-rock definition, geotechnical properties (i.e.softness-stiffness) of soil and rocks in a site (Socco and Strobbia 2004;Zhan et al. 2008;Foti et al. 2011;Dumont et al. 2018).The most important advantage of this method over others methods is that the near-surface resolution (top 10 m) is much improved (Kondracka et al. 2021) and is not affected by velocity reversal problem as in seismic refraction (Glangeaud et al. 1999;Park et al. 1999).In MASW, the dispersion of surface wave field (Rayleigh type surface wave in this study) is concerned and generally the recording time, and sampling time can be chosen to be longer than the seismic refraction.The dispersion of the surface waves is extracted from the data presented in Fig. 7a by applying the frequency-velocity (f-v) transform (phase-shift technique in this study) (Fig. 7b), and the velocity value corresponding to the maximum energy at each frequency is selected as the phase velocity, named as fundamental mode dispersion curve (Fig. 7c).The 1D Vs profile in Fig. 7d also is obtained by inversion of the fundamental mode dispersion curve.This resulting V S -depth profile is assigned to the midpoint of the receiver array and is used to calculate V S30 (average S-wave velocity for 30 m depth) given by Eq. ( 9) (Borcherdt 1994) which is necessary for soil classification of a site.In this study, soil class assessment for potential landfill areas was accomplished according to TBEC (Türkiye Building Earthquake Code) (2018) adapted from NEHRP (National Earthquake Hazard Reduction Program) (BSSC 2003) (Table 6).According to the new classification of NEHRP (BSSC 2020), B, C and D soil classes are divided into subclasses such as BC, CD and DE.However, TBEC (2018) fully matches the old classification of NEHRP, but Table 6 is arranged in this study, taking into account the new classification of NEHRP.The expression letter "Z" in TBEC is the Turkish equivalent of "Soil" in English and is used as a prefix.Thus, by using V S30 , the parameters (amplification, dominant vibration frequency or period values) necessary to characterize the behaviour of soils in the potential landfill area under dynamic loads can be easily calculated using experimental formulas commonly used in the literature: In the equation, h i and V si are respectively layer thickness and shear wave velocity for each layer.N indicates the number of layers up to 30 m deep.
Electrical resistivity tomography method is a geophysical method considering the variation of the electrical resistivity of subsurface materials with the widest use in subsurface explorations and solution of geotechnical problems such as characterizing the subsurface structure of landfills (Vargemezis et al. 2015); detecting groundwater, clay content and seepage deposits (Campbell et al. 1999;Kaya et al. 2007); examining the potential sites for solid waste disposal facilities (Smith and Randazzo 2003;Bichet et al. 2016;El-Kelani and Khader 2019;Hu et al. 2019) and evaluating the locations of regions with different moisture levels in landfills, (9)  7 Steps for MASW method.a Data acquisition: shot record from the field.b data processing: f-v imaging by phase-shift transformation technique.c observed (or selected) and calculated (or inverted) fundamental mode curve.d 1D V S -depth profile.This profile may be expressed as a shear wave sounding and is generally placed in the middle of the geophone spreading Table 6 Site soil classification codes and definitions according to the V S30 value were adapted from TBEC (2018) and NEHRP (BSSC 2020).In the table, the soils include sand, clay and gravel units.The definition of the lithology is attributed to the mechanical properties, and borehole logs are not used For site class code E or ZE: Soil or any profile with more than 3 m of soft clay is defined as soil with PI > 20, w ≥ 40% and c u < 25 kPa.For site class code F or ZF: soils with the risk of collapse and potential collapse under the influence of earthquakes (liquefiable soils, highly sensitive clays, collapsible weak cemented soils, etc.), clays with a total thickness of more than 3 m of peat and/or high organic content, high plasticity (PI > 50) clays with a total thickness of more than 8 m, very thick (> 35 m) soft or medium solid clays PI Plasticity Index, w water content, c u undrained shear strength fractured-cracked structure, permeable-impermeable, mapping resistant and weak lithological units (Guerin et al. 2004;Loke 2004;Olona et al. 2010;Hossain et al. 2011).Therefore, in this study, ERT method was used specially to map the underground water and clay content and weathering degree of the rocks.The methodology of the ERT data acquisition and evaluation systematic is demonstrated in Fig. 8. Shortly, Fig. 8a shows the current and potential electrode positions (Wenner-Schlumberger in this study) within the multi-electrode systems which provide fast and less laborious (Loke and Barker 1996;Loke 1997) data acquisition.Figure 8 b shows the distribution of data points corresponding to the depth levels penetrated for different electrode intervals.To account for topographic changes over wide spans, the elevation values of each electrode are written to the data file.Obtained apparent resistivity values are converted to real resistivity values of the ground by tomographic inversion technique.By applying this process in 2D, ERT sections are created, and they are now ready for interpretation (Fig. 8c).When the resistivity changes are evaluated together with additional information (e.g.seismic velocity), it has the potential to make a significant contribution to the general description of both the physical and lithological nature of the geological material (soil: sand, clay, silt, gravel and rock: volcanic, sedimentary).The ERT data were acquisited with Wenner-Schlumberger electrode array using 23 electrodes of the Ambrogeo trademark, Mangusta model 48-electrode device.This array was used because it is highly sensitive to both horizontal and vertical structural-lithological changes and has a wider horizontal coverage and signal strength than other arrays.The electrode spacings were 3.0 m (Işıktepe-1, 2 in Ordu; Ağalık-1, 2 in Giresun; Ovacık-1, 2 in Trabzon), 4.0 m (Kazantaş-2 in Gümüşhane) and 5.0 m (Vezirköprü-1,2 and Bafra-1,2 in Samsun; Esence-1,2 in Ordu; Çamburnu-1,2 in Trabzon; Murgul-1,2 in Artvin; Yenice-1,2 and Kazantaş-1 in Gümüşhane; Center-1,2 and Balkaynak-1 in Bayburt), depending on the terrain conditions at the measurement sites.

Results
Within the scope of this study, quantitative and qualitative analyses of the seismic and electrical resistivity data acquisited in the seven provinces in 25 profiles were performed, and 13 of the obtained cross sections are presented in Figs. 9 and 10.While obtaining SRT, ERT sections and 1D-V S profiles, the RMS (root mean square) error amounts as a result of inversion were less than 3%, 7% and 5%, respectively.According to the V S30 value for all sites, the classification, dominant vibration period (T 0 ) and amplification (A 0 ) of the soil and physical description of the soil are presented in Table 7. Also, in Table 7, Ta and Tb values calculated by using T 0 indicate the corner frequencies required for the construction static project.According to these frequency values, ground-structure resonance period ranges are determined.In addition, the interpretations of the obtained SRT and ERT cross sections are presented in Table 8.It can be said that the results obtained from SRT and ERT are generally  According to the SRT section, the ground was very soft from the surface to a depth of 12 m (no bedrock!).According to ERT data, this low resistivity implies that the medium is highly saturated.This site is very risky in terms of ground settlement, collapse, fracture and liquefaction in case of both vertical and horizontal loads during the passage of future earthquake waves Moderately resistive soils According to the seismic parameters, the environment has a very solid soil near its surface.Especially after 3.2 m depth, the P-wave velocity implies a hard rock structure.The ERT section show local moistening at near the surface, but these were considered to be seasonal.Therefore, it is believed that no groundwater exists in the environment

66-173
Resistive soils There is massive bedrock starting at a depth of about 3.0 m.Moreover, the resistivity data show that the medium is dry  compatible.As seen in ERT sections, the range of resistivity values is quite wide (~ 2-2000 Ohm-m).Therefore, in the presentation of the ERT sections, it used individual colour scale to clearly visualize the anomalies.On the other hand, field observations and regional geology were taken into account while interpreting the ERT and SRT sections.
The subsurface geometry was defined as a three layers according to the P-wave velocity variation in the SRT sections.In particular, the bedrock depth and topographic changes were indicated in the sections, and the mechanical properties of the soil and rock layers (softness-stiffness of the soils and weathering degree of the rocks) were evaluated according to the P-wave velocity change utilized by the table of the seismic velocity-rippability-geotechnical definition given by Karslı et al. (2021).

Discussion
In order to determine the suitable solid waste landfill areas, seismological and near-surface geophysical measurements were carried out in alternative areas determined by preliminary studies in different provinces of the Eastern Black Sea Region.Since this study was carried out with the aim of making the suitable site selection of the waste landfill areas, it differs from the studies on the determination and solution of the problems which arise in the existing or abandoned landfills.These differences are not related to the geophysical methods used and the processing of the obtained data, but to the interpretation of the findings obtained from the data in order to prevent problems that may arise after landfill.
The statistical and deterministic evaluations (earthquake hazard analysis) based on seismological data were made by considering the geographical boundaries of the provinces where the landfill areas foreseen for solid waste are located.Therefore, rather than the locations of the landfill areas (or surveying locations), the distance of the geographical boundaries of the provinces to the earthquake epicentre (especially for NAFZ) was effective in these assessments.It should be noted that, although all other provinces (Samsun, Ordu, Giresun and Artvin) have a coast to the Black Sea, except for Gümüşhane and Bayburt provinces (Trabzon Province is bordered by Gümüşhane and Bayburt Provinces from the south), the southern parts of these provinces are quite close to the NAFZ.On the other hand, the foreseen landfill areas are mostly located in the northern parts of the provinces, in areas closer to the city centres.Therefore, according to Table 2,  The fact that the b-values calculated for Ordu, Gümüşhane and Trabzon are normally smaller than for the other provinces, especially those close to the NAFZ, is related to the limited number of earthquakes occurring within the municipal borders and near around of Ordu, Gümüşhane and Trabzon, especially with M s ≥ 6.0.Therefore, it may be misleading to make a comparative qualitative seismicity assessment of the b-values calculated for each province.Taking into account the values of parameter "b" determined locally, it is clear that the frequency of earthquakes causing the damage in the study areas is low and in good agreement with the probability calculations presented in Table 3.On the other hand, considering the probability of earthquake occurrence in 10-year period values in Fig. 5 and the b-value in Eq. ( 6), it is understood that the seismic hazard of the region is high, especially for a mediumsized earthquake in the near future.Furthermore, when evaluated independently of the soil, it is clear that the intensity and acceleration values in Table 4 were high in locations close to the earthquake epicentre and diminished as they moved away from it.However, the intensity and acceleration values may increase in local soils, far from the earthquake epicentre, which may cause the potential for soil amplification, increasing the vibration period and liquefaction and focusing due to the basement (or bedrock) topography.This situation requires site characterization and is explained by the nonlinear behaviour of the soils.
From the near-surface geophysical measurements, seismic measurements contributed to the soil characterization of the landfill areas, while electrical measurements provided information about water content, lithological changes and degree of rock weathering.In fact, this information directly explains the structure and behaviour of the subsurface against possible loads, in both vertical and horizontal direction.Accordingly, near-surface geophysical measurements performed in the site selection of solid waste landfill areas, compared with other determination methods (city planning, geological, hydrogeological, environmental science, social and health sciences and geographic information systems, etc.), have important advantages in terms of taking geotechnical measures for both the construction design of the landfill areas and the precaution of structural problems which may arise later.
Therefore, according to results shown in Table 7, the lowest V S30 value was measured in Bafra (V S30 = 223 m/s) and the highest in Işıktepe (V S30 = 1327 m/s).Although soil classes varied between D and B, it was mostly observed CD, C and BC soil classes for NEHRP (BSSC 2020) and ZC and ZB soil classes for TBEC (2018).These soil classes indicate the existence of geological environments containing moderately stiff to very stiff soil layers, which are not thick and cover less weathered rocks at the basement, which is consistent with the general geology of the study areas.The soft soils (D or ZD) are in Bafra and Ağalık, while the hard or medium hard rocks (B or ZB) are in Işıktepe and Çamburnu.Some soils (Şebinkarahisar-2, Kazantaş-2, Bayburt Merkez-2, Murgul) are ZB (less weathered or medium hard rocks) according to TBEC (2018) in the subclass BC (soft rocks) according to NEHRP (BSSC 2020).In such soils, it will be more appropriate to perform the necessary calculations considering the BC soil class in designed geotechnical projects in terms of construction safety.On the other hand, the minimum-maximum of them was calculated as T 0 = 0.09-0.47s and A 0 = 1.06-2.47,respectively.The soil amplification values were calculated according to the Midorikawa (1987) approach, and the lower and upper limits of the soil dominant vibration period value (T a = 0.5*T 0 , T b = 1.5*T 0 ), design period and frequency values, were obtained according to the quarter wavelength approach of Keçeli and Cevher (2018).For all profiles, the minimum and maximum T 0 were for Işıktepe (T 0 = 0.09 s) and Bafra (T 0 = 0.53 s), respectively, while the minimum and maximum A 0 were calculated at Çamburnu (A 0 = 1.06) and Bafra (A 0 = 2.65).A 0 < 1.0 means that there is no amplification in the soil, although numerically, the smallest soil amplification (A 0 = 0.9) in Table 7 is obtained for Işıktepe.As can be seen, the largest values for soil amplification and soil dominant vibration period were calculated at the Bafra location.
From the ERT sections, information was obtained regarding the water content (humidity and saturation) and possible clay content of the geological units according to the resistivity changes.The presence of water-saturated units in the ERT section corresponding to the fault locations determined in the SRT section was remarkable.This is associated with the fact that the faults are important transmission ways that allow the surface waters to infiltrate deep and the aquifer layer to be formed.Thus, the information provided by each geophysical data point in terms of soil and rock characteristics was interpreted to eliminate possible uncertainties.In this context: 1. Zones corresponding to V P = 1300-1550 m/s observed in SRT sections and low resistivity (ρ < 10 Ohm-m) in ERT sections were considered as saturated with water 2. Zones with low velocity (V P < 1000 m/s) but high resistivity (ρ > ~ 70 Ohm-m) are considered to be dry soil units 3. Zones with moderate velocity (V P = 1200-2100 m/s) and high resistivity are interpreted as very compact soil and soft and/or weathered rocks, which generally do not contain water 4. Zones with high velocity (V P > 2100 m/s) but low resistivity were evaluated considering the possibility that the rocks contain water or clay in their cracks or pores or conductive minerals 5.The shear wave velocity (V s ) is particularly sensitive to the mechanical properties of the soil, but very little sensitive to its water content.However, the V s value can be lower than 100 m/s, especially in soft clays, while it is more than 1000 m/s in stiff soil and rock materials.Thus, the 1D-V s depth profile was used to determine the softness-stiffness soils, soil-rock boundaries and soil classification 6.The ERT sections have been very helpful in evaluating the unconsolidated structure of the soil material at the measurement sites, the degree of weathering of the rocks, the presence of clayey intermediation and the water content, but they have not been successful enough to describe the soil-rock interfaces.However, these interfaces were identified from the SRT sections and marked on the ERT sections.On the other hand, the bedrock depth levels were precisely compatible in the SRT and ERT sections, particularly when the rocks were dry

Conclusions
This study includes seismicity analyses and subsurface characterization using near surface geophysical measurements (SRT, MASW, ERT) for planned waste landfill sites in seven provinces in the Eastern Black Sea Region of Türkiye.Seismicity analyses provided the probability of occurrence of earthquakes with the magnitude of medium, strong and large scale (M s ≥ 5.0, 6.0, 7.0) which can negatively affect the study areas as well as their recurrence periods.Furthermore, the acceleration and intensity predictions for two scenario earthquakes were performed.These analyses show that the seismic hazard of the region is quite high, especially for a medium-sized earthquake, at a high level for strong earthquakes and relatively low for major earthquakes, according to the 10-and 50-year recurrence probabilities.On the other hand, the calculated soil-independent acceleration and intensity values were quite high in provinces close to the scenario earthquake locations and low in the distant provinces.
While SRT showed soil geometry, bedrock depth and physical properties of geological materials, MASW effectively provided soil classification, soil-rock contrast and their mechanical behaviour (softness-stiffness profile, amplification and vibration dominant period) at 30-m depth.ERT is extremely effective in detecting moist, saturated and dry zones.The combined interpretation of SRT, MASW and ERT allowed to significantly reduce the inadequacy or uncertainty of interpretation possible with only one method.Thus, the subsurface structural condition of the possible landfill areas and the characterization of the soil and rocks were performed accurately and decisively.In particular, while the SRT and ERT results were highly compatible in determining the bedrock depth in dry environments, the ERT results were not successful in the presence of groundwater.According to the geophysical results, four sites (Işıktepe, Esence, Çamburnu and Kazantaş) were found to be suitable, and six sites (Vezirköprü, Şebinkarahisar, Yenice, Bayburt-Center, Balkaynak and Murgul) were found to be suitable after necessary geotechnical reclamation.Three sites (Bafra, Ağalık, Ovacık) were considered to be unsuitable because of the presence of thick, water-saturated soft soil and extremely weathered rocks.Soils at unsuitable sites are risky because they have the potential to amplify earthquake waves and increase the vibration period and are susceptible to possible deformations and landslides after landfilling.
As a result, this study has shown that the seismicity of possible waste landfill sites and the physical and mechanical properties of geological materials at these sites can be investigated in a non-invasive, rapid and easy way using multiple geophysical methods.This information will enable us to minimize and/or prevent possible deformations (fracturing, collapse, settlement, sliding, etc.), environmental problems (such as dirty water leaks) and economic losses after landfilling.Therefore, municipal governments will be able to make more accurate and less risky decisions in making final decisions about waste landfill sites, which are extremely important for the planning and development of a city.

Fig. 5
Fig.5The occurrence probability of earthquakes in 10-year periods for different magnitudes in the region

Fig. 6
Fig. 6 Methodology of the SRT method.a Data acquisition layout in the field.b The data according to source position (shots 3 and 4 are not shown and the time axis is truncated at 200 ms because of displaying purposes).c Observed and calculated (by inversion) first

Fig. 8
Fig. 8 Multi-electrode ERT data acquisition scheme.a Current (I; C 1 -C 2 )-potential electrodes (V; P 1 -P 2 ), connections and current flow lines.b Distribution of apparent resistivity data points read from the subsurface by Wenner-Schlumberger electrode array.c 2D ERT section from tomographic inversion.The black dots in b represent the

Fig. 9
Fig. 9 Sections of selected profiles from each survey location.a SRT and b ERT.The bedrock level (dotted line-dotted line) on the SRT sections is also indicated in the ERT sections, considering the approximately equivalent depth level.Saturated, moist and dry zones Figs. 9 and 10 and evaluation of suitability in terms of the waste landfill site According to the resistivity change, soft and compact soil units are mostly dry and partially moist at distances of 10-55 m, whereas the zone between 55 and 105 m distance and 9.0 m depth from the surface is saturated.The fact that V P = 1.4-1.6 km/s and ρ = 3-6 Ohm-m in this range indicates that the geological material is clay-rich.A clear discontinuity (probably a fault) is observed in the SRT section at 30-35 m.After 9 m depth, there are rock units suitable for foundation.Although the environment is considered suitable as a landfill area in terms of seismic parameters, resistivity data indicate the presence of groundwater up to 4.5 m deep in the environment soils The maximum-minimum resistivity variation was approximately 4-24 Ohm-m.
Medium-stiff soil and weathered rock were observed from the surface to a depth of 7 m.A discontinuity zone (probably lithological or facies change) was observed at intervals of 25-30 m in the SRT section and 65-70 m in the ERT section.However, the low resistivity zones (< 25 Ohm-m) within the first 10 m are thought to be caused by surface seepage water which moistens the weathered soils Units with a resistivity of ρ ≤ 20 Ohm-m observed from the surface to a depth of 13 m indicate the presence of clay-containing, water-saturated units.The low resistivity of the deeper rock material is associated with its highly weathered structure and therefore the presence of water seeping from its cracks.The risk of the site for any construction is are stiff and dry.The units below show very stiff and weathered rock characteristics.In addition, the decrease in the resistivity values of the rocks, which appear at high velocity from the surface to the depth, especially at distances beyond 55 m, was considered to be caused by water leaking from the thick and dry soil units up to a depth of 5 m in the environment.In the SRT section (Fig.10a(ii)), a fault, possibly a normal fault, was observed in the 30-35 m range, causing lateral lithological change.In addition, the 1.4-1.8km/s velocity unit placed at ~ 9 m from the surface after 30 m distance is seen as low resistivity (30-70 Ohm-m) in the ERT section (Fig.10b(ii)).This was clearly caused by the infiltration of seasonal water from the surface into soil materials or weathered require operations that include groundwater drainage, soil improvement (compaction, jet grouting, injection, pilling, etc.) and removal of soft top soil cover Table8 to a depth of 3.5 m, medium-stiff-very stiff dry soils and weathered rocks were observed, while deeper (> 3.5 m) dry and rather hard rocks The environment was very evenly bedded, and a medium-stiff-very stiff soil structure was observed from the surface to a depth of ~ 15 m.After this depth, moderately weathered dry rock units were placed.According to the ERT section, resistivity levels of 10-35 Ohm-m indicate moistened the SRT and ERT sections, the environment showed very dry soil characteristics up to a depth of 8 m.After this depth, the materials show rock structure and have massive character.However, the low resistivity of the massive rock unit at a depth of more than 8 m was associated with the fact that the rock contained conductive minerals rather than water, considering where the location of this site is an old mine (with sulfide mineral content) Especially, in the 3-9-m depth range, the average P wave velocity of 1500 m/s and the low resistivity value (5-13 Ohm) indicate the presence of water-saturated soil material here.In addition, a lateral discontinuity indication was observed in the 30-35-m distance range in the SRT section.Although the bedrock in this area is hard, geotechnical problems are likely to occur up to 9-m depth.Moreover, observed in the ERT section, the slope of both the surface and bedrock topography indicate a to a depth of 5 m, there are units with high moisture content, and deeper, resistant, moderately weathered and massive rock units are located.In addition, the slope of both the surface and bedrock topography and the water content of the materials on the bedrock indicate the presence of landslide the small b-values characterize the municipalities of Ordu (b = 0.47), Gümüşhane (b = 0.48) and Trabzon (b = 0.48), and the large b-values are characteristic of the provinces of Artvin (b = 0.87), Bayburt (b = 0.68), Samsun (b = 0.65), Giresun (b = 0.57) respectively.

Table 1
The large historical earthquakes (M s ≥ 7.0) in the

Table 3
The probability of an earthquake with a magnitude of M s ≥ 4.0 in the study area and its recurrence period values Magnitude (M s ) Seismic hazard, R(M)% period (year)

Table 4
Predicted intensity and acceleration values for each province (1 g = 1000 cm/s 2 )

Table 5
Equipment and parameters for field data acquisition in the study

Table 7
Soil classifications according to the V S30 value and calculated soil parameters.Although the distance between profiles was variable, it was less than 100 m

Table 8
Interpretation of SRT and ERT sections in