Geology of the ignimbrites and the associated volcano–plutonic complex of the Ezine area, northwestern Anatolia

The Ezine region is located in the northwestern part of Anatolia where young granitic and volcanic rocks are widespread and show close spatial and temporal association. In this region magmatism began with the Kestanbol granite, which intruded into metamorphic basement rocks, and formed contact metamorphic aureole. To the east and southeast the pluton is surrounded by hypabyssal rocks, which in turn, are surrounded by volcanic associations. The volcanic rocks may be divided into two main groups on the basis of their lithological properties. Lavas and lahar deposits dominate the northern sector while ignimbrites dominate the southern sector. The ignimbrite eruptions were formed partly coevally with the plutonic and the associated volcanic rocks during the early Miocene. They appear to have been associated in a caldera collapse environment. Geochemical properties of the plutonic and the associated volcanic assemblages indicate that the magmas are hybrid and co-genetic and, were formed from a similar mantle source, under a compressional regime prior to the opening of the present E–W-trending graben of the Aegean western Anatolian region.     

Relationships between volcanic patterns and neotectonics in Eastern Anatolia from analysis of satellite images and DEM

The late Neogene to Quaternary volcanism in Eastern Anatolia is related to the Arabia–Eurasia convergence but a clear deformation pattern has not yet been established in this region. We have used the distribution and shape of volcanoes and fault geometry as indicators of the tectonic regime. Volcanic edifices and related faults were analyzed in vertical view using SAR–ERS, Spot images and a Digital Elevation Model (DEM). In several places, adjacent volcanoes that form linear clusters or elongated volcanoes are clearly rooted on vertical tension fractures. These are compatible with horizontal σ₃ striking 90°N, associated with σ₁ horizontal (strike-slip regime) or vertical (extensional regime). We mapped the recent faults that are directly associated to volcanoes. Volcanic vents are related to tail-crack, horsetail or releasing bend structures. In this work, it has been possible to define the ESE-striking, 270-km-long Tutak–Hamur–Çaldiran fault that forms a releasing bend testifying to right-lateral motion. Extension is well documented for few places but no recent fold has been observed. Since 8 Ma, the tectonic system is principally strike-slip. Most of the tension fractures being 2 to 10 km in length, so we infer that they affect only part of the crust. Most strike-slip fault zones are of several tens to a few hundred kilometers long and thus not of lithospheric scale. Therefore, the channels used by the magma to reach the surface are crustal structures.     

Magma evolution of Quaternary minor volcanic centres in southern Peru,Central Andes

Minor centres in the Central Volcanic Zone (CVZ) of the Andes occur in different places and are essential indicators of magmatic processes leading to formation of composite volcano. The Andahua–Orcopampa and Huambo monogenetic fields are located in a unique tectonic setting, in and along the margins of a deep valley. This valley, oblique to the NW–SE-trend of the CVZ, is located between two composite volcanoes (Nevado Coropuna to the east and Nevado Sabancaya to the west). Structural analysis of these volcanic fields, based on SPOT satellite images, indicates four main groups of faults. These faults may have controlled magma ascent and the distribution of most centres in this deep valley shaped by en-echelon faulting. Morphometric criteria and ¹⁴C age dating attest to four main periods of activity: Late Pleistocene, Early to Middle Holocene, Late Holocene and Historic. The two most interesting features of the cones are the wide compositional range of their lavas (52.1 to 68.1 wt.% SiO₂) and the unusual occurrence of mafic lavas (olivine-rich basaltic andesites and basaltic andesites). Occurrence of such minor volcanic centres and mafic magmas in the CVZ may provide clues about the magma source in southern Peru. Such information is otherwise difficult to obtain because lavas produced by composite volcanoes are affected by shallow processes that strongly mask source signatures. Major, trace, and rare earth elements, as well as Sr-, Nd-, Pb- and O-isotope data obtained on high-K calc-alkaline lavas of the Andahua–Orcopampa and Huambo volcanic province characterise their source and their evolution. These lavas display a range comparable to those of the CVZ composite volcanoes for radiogenic and stable isotopes (⁸⁷Sr/⁸⁶Sr: 0.70591–0.70694, ¹⁴³Nd/¹⁴⁴Nd: 0.512317–0.512509, ²⁰⁶Pb/²⁰⁴Pb: 18.30–18.63, ²⁰⁷Pb/²⁰⁴Pb: 15.57–15.60, ²⁰⁸Pb/²⁰⁴Pb: 38.49–38.64, and δ ¹⁸O: 7.1–10.0‰ SMOW), attesting to involvement of a crustal component. Sediment is absent from the Peru–Chile trench, and hence cannot be the source of such enrichment. Partial melts of the lowermost part of the thick Andean continental crust with a granulitic garnet-bearing residue added to mantle-derived arc magmas in a high-pressure MASH [melting, assimilation, storage and homogenisation] zone may play a major role in magma genesis. This may also explain the chemical characteristics of the Andahua–Orcopampa and Huambo magmas. Fractional crystallisation processes are the main governors of magma evolution for the Andahua–Orcopampa and Huambo volcanic province. An open-system evolution is, however, required to explain some O-isotopes and some major and trace elements values. Modelling of AFC processes suggests the Charcani gneisses and the local Andahua–Orcopampa and Huambo basement may be plausible contaminants.     

petrological and geochemical characteristics of quarternary volcanic rocks in haixing area,eastern north china

Field investigation and lab analysis on samples were carried out for Quaternary volcanoes, including Xiaoshan volcano, Dashan volcano and Bianzhuang hidden volcano, in Haixing area, east of North China. Results show that Xiaoshan volcano with the eruptive material of volcanic scoria, crystal fragments and volcanic ash is a maar volcano, the eruptive pattern is pheatomagmatic eruption, and the influence scope is near the crater. Dashan volcano exploded in the early stage, and then the magma intruded, forming the volcanic neck. The eruption strength and scale are limited, and the eruptive materials are scoria, volcanic agglomerate and dense lava neck. The volcanic rocks in Bianzhuang are porosity and dense volcanic rocks and volcanic breccia, reflecting the pattern of weak explosive eruption and lava flow, and the K-Ar age dating on volcanic rocks indicates that the eruption happened in early Pleistocene. Xiaoshan volcanic scoria and Bianzhuang hidden volcanic rocks are mainly basaltic, Dashan volcanic rocks with lower SiO₂ content are nephelinite in composition. Their oxide contents have no linear relationship, indicating that there is no magma evolution relationship between these magmas from the three places. Three volcanic rocks all have enrichment of light rare earth. The Bianzhuang volcanic rocks are rich in large ion lithophile elements, and have no high field strength elements Zr and Hf, Ti losses. The volcanic materials from Xiaoshan and Dashan are intensively rich in Th, U, Nb and Ta, and significantly poor in K and Ti. Although the magmas from these three places in Haixing area may all come from asthenosphere, the volcanic materials have different petrological and geochemical features, and relatively independent volcanic structures, therefore, they experienced different magma processes.     

Evolution of the Bayramiç magmatic complex, northwestern Anatolia

During the Oligocene–Middle Miocene period widespread magmatic activity developed in Western Anatolia, following the continental collision between the Sakarya continent and the Tauride–Anatolide platform. This produced both intrusive and extrusive rocks, which appear to be associated in space and time, as exemplified from the Bayramiç area. In the Bayramiç area, the magmatic activity started with the intrusion of the Evciler granite, and the coeval lower volcanic association. This was followed by the development of the upper volcanic association. These rock groups form collectively the Bayramiç magmatic complex, which was generated under an on-going north–south compressional regime. The Bayramiç magmatic complex has a subalkaline composition, displaying a calcalkaline trend. Trace elements and REE contents resemble to island-arc and collision-related magmas. According to the isotope values the Bayramiç magmatic complex was derived from the magmas of lithospheric mantle origin, which were later contaminated, while passing through the thick continental crust, in a post-collisional tectonic setting, during the Oligocene–Early Miocene period. The latest product of the magmatism is the Late Miocene–Pliocene basalt lavas. Their geochemical properties are clearly different from the Oligocene–Early Miocene magmatic rocks. The basalts were generated when the north–south compression gave way to the north–south extensional regime.     

Curie Point Depth Investigation of Central Anatolia,Turkey

The residual aeromagnetic total field intensity anomalies in central Anatolia were calculated from the regional aeromagnetic anomalies surveyed by the Mineral Research and Exploration (MTA) of Turkey. The residual aeromagnetic data were analyzed to produce Curie point estimates by the method of OKUBO et al. (1985). The Curie point depth of central Anatolia varies from 7.9 km and 22.6 km. The shallowest Curie point depths were observed around the Cappadocia and Erciyes Volcanic complexes in central Anatolia. A good correlation was deduced between the Curie point depths and the heat-flow data measured previously, which is most certainly important for the geothermal resources of the region. The shallow Curie point depths also correlate well with the hot spring locations in central Anatolia.     

The Effect of Regional Borders when Using the Gutenberg-Richter Model, Case Study: Western Anatolia

A probabilistic seismic hazard analysis (PSHA) for Western Anatolia is presented using the Gutenberg-Richter (G-R) frequency-magnitude relation. Since the modeling is sensitive to the location of seismotectonics boundaries, to use the information content of the observed earthquake data, as a general rule the borders of the affected area are extended. In this study, the effect of the region’s definition on the G-R model is debated on the Western Anatolian region, which is one of the most seismically active and rapidly deforming regions of the world throughout the ages. Calculations are carried out for two subregions and one combined region as a whole using the seismic catalog from 1900 to 2005. The data sets are determined by the region’s borders with the parameters computed according to these data sets by the least-squares and maximum likelihood methods, and then future predictions are estimated via these parameters. Comparing the results with historical earthquake records, most appropriate regional borders for Western Anatolia are defined, and for this new region G-R model parameters are obtained.     

The effects of secular disequilibrium on (U-Th)/He systematics and dating of Quaternary volcanic zircon and apatite

The (U-Th)/He dating method applied to U-rich phases such as zircon and apatite has sufficient sensitivity and precision to be of potential use for dating relatively recent geologic events such as volcanic eruptions. However, in phases with crystallization ages less than ∼1 Ma, chemical fractionation within the ²³⁸U decay series may modify the He ingrowth rate, causing He ages computed from the secular equilibrium age equation to be incorrect. The resulting systematic error depends on the [²³⁰Th/²³⁸U] activity ratio of the dated phase when it is erupted, and on the eruption age. Zircons, which exclude Th relative to U, will likely have secular equilibrium He ‘ages’ that underestimate the eruption age by up to a few tens of %, decreasing with increasing eruption age. Apatites tend to accommodate U and Th with little fractionation, so apatite secular equilibrium He ages will be nearly concordant with eruption age. If minerals are erupted immediately after crystallization, the disequilibrium effect can be reasonably accounted for based on Th/U systematics. However, crystals are likely to reside for unknown but potentially long periods in a magma chamber, such that the degree of secular disequilibrium will be reduced prior to the onset of He accumulation. (U-Th)/He analyses of co-genetic phases that fractionate the U/Th ratio differently, like apatite and zircon, can be used to better constrain eruption age, as well as to provide insights into magma chamber residence time. We illustrate this approach with (U-Th)/He analyses of zircons and apatites of the Pleistocene-age Rangitawa Tephra, New Zealand.     

Geochemical characteristics of Yamadağı volcanics in central east Anatolia: an example from collision-zone volcanism

Neogene Yamadağı volcanic rocks consist of basaltic trachyandesite, trachyandesite, andesite, and dacite. The major- and trace-element chemistry indicates that the lavas are dominantly calc-alkaline and mildly alkaline in character, sodic in nature, and intermediate to acidic in composition. REE and trace-element patterns of volcanic rocks are similar to those typical of within plate magmatics. Volcanic rocks have low ⁸⁷Sr/⁸⁶Sr (0.70389–0.70633) and high ¹⁴³Nd/¹⁴⁴Nd ratios (0.51267–0.51276) and mostly plot within the mantle array of the isotope ratio diagram. The linear correlations among ⁸⁷Sr/⁸⁶Sr−¹⁴³Nd/¹⁴⁴Nd, SiO₂−⁸⁷Sr/⁸⁶Sr and SiO₂−¹⁴³Nd/¹⁴⁴Nd isotope ratios in the volcanics suggest that fractional crystallization combined with minor assimilation was an important process within the collision zone.     

An application of the time- and magnitude-predictable model to long-term earthquake prediction in eastern Anatolia

In order to estimate the recurrence intervals for large earthquakes occurring in eastern Anatolia, this region enclosed within the coordinates of 36^∘–42^∘N, 35^∘–45^∘E has been separated into nine seismogenic sources on the basis of certain seismological and geomorphological criteria, and a regional time- and magnitude-predictable model has been applied for these sources. This model implies that the magnitude of the preceding main shock which is the largest earthquake during a seismic excitation in a seismogenic source governs the time of occurrence and the magnitude of the expected main shock in this source. The data belonging to both the instrumental period (M_S≥ 5.5) until 2003 and the historical period (I₀≥ 9.0 corresponding to M_S≥ 7.0) before 1900 have been used in the analysis. The interevent time between successive main shocks with magnitude equal to or larger than a certain minimum magnitude threshold were considered in each of the nine source regions within the study area. These interevent times as well as the magnitudes of the main shocks have been used to determine the following relations:

Alkali mapping of the japanese quaternary volcanic rocks

Of ca. 200 Quaternary volcanic centers in the Japanese islands, ca. 120 are characterized by K₂O and Na₂O concentrations of their ejecta normalized to SiO₂ = 60% K₂O values give a geographically well-defined trend of regularly increasing away from the volcanic front. Na₂O values show no geographical trend. Relationships between the SiO₂-normalized K₂O and the depth of the deep seismic zone (K-h relation) are different for different arcs. The Southwest Japan arc shows distinctly higher K₂O and less regular geographic trend than the others probably due to crustal contamination. In the normal segments of the Northeastern Japan belt, K₂O increases ca. 2.5 to 3.0% per 100 km depth of the seismic zone. Groups of volcanic centers near the junction between the Northeast Honshu and Izu-Mariana arcs give characteristically lower K₂O while Rishiri volcano far behind the arc junction between the Kurile and Honshu arcs also gives a low K₂O value.     

Geochemical approach to magmatic evolution of Mt. Erciyes stratovolcano Central Anatolia, Turkey

Erciyes stratovolcano, culminating at 3917 m, is located in the Cappadocian region of central Anatolia. During its evolution, this Quaternary volcano produced pyroclastic deposits and lava flows. The great majority of these products are calc-alkaline in character and they constitute Kocdag and Erciyes sequences by repeated activities. Alkaline activity is mainly observed in the first stages of Kocdag and approximately first-middle stages of Erciyes sequences. Generally, Kocdag and Erciyes stages terminate by pyroclastic activities. The composition of lavas ranges from basalt to rhyolite (48.4–70.5 wt.% SiO₂). Calc-alkaline rocks are represented mostly by andesites and dacites. Some compositional differences between alkaline basaltic, basaltic and andesitic rocks were found; while the composition of dacites remain unchanged. All these volcanics are generally enriched in LIL and HFS elements relative to the orogenic values except Rb, Ba, Nb depleted alkaline basalt. ⁸⁷Sr/⁸⁶Sr and ¹⁴³Nd/¹⁴⁴Nd isotopic composition of the volcanics range between 0.703344–0.703964, 0.512920–0.512780 for alkaline basalts and change between 0.704322–0.705088, 0.512731–0.512630 for alkaline basaltic rocks whereas calc-alkaline rocks have relatively high Sr and Nd isotopic ratios (0.703434–0.705468, 0.512942–0.512600). Low Rb, Ba, Nb content with high Zr/Nb, low Ba/Nb, La/Yb ratio and low Sr isotopic composition suggest an depleted source component, while high Ba, Rb, Nb content with high La/Yb, Ba/Nb, low Zr/Nb and low ⁸⁷Sr/⁸⁶Sr ratios indicate an OIB-like mantle source for the generation of Erciyes alkaline magma. These elemental and ratio variations also indicate that the different mantle sources have undergone different degree of partial melting episodes. The depletion in Ba, Rb, Nb content may be explained by the removal of these elements from the source by slab-derived fluids which were released from pre-collisional subduction, modified the asthenospheric mantle. The chemically different mantle sources interacted with crustal materials to produce calc-alkaline magma. The Ba/Nb increase of calc-alkaline samples indicates the increasing input of crustal components to Erciyes volcanics. Sr and Nd isotopic compositions and elevated LIL and HFS element content imply that calc-alkaline magma may be derived from mixing of an OIB-like mantle melts with a subduction-modified asthenospheric mantle and involvement of crustal materials in intraplate environments.     

Investigation into the regional wrench tectonics of inner East Anatolia (Turkey) using potential field data

The residual aeromagnetic and gravity anomalies of inner East Anatolia, surveyed by the Mineral Research and Exploration (MTA) of Turkey, display complexities. Some faults, which are known and new lineaments, are drawn from maxspot map derived from the location of the horizontal gradient of gravity anomalies. Tectonic lineaments of inner East Anatolia exhibit similarities to the direction of East Anatolian Fault Zone. Anticlockwise rotation, approximately −30°, defined from disorientations of aeromagnetic anomalies. The lineaments obtained from maxspots map produced from the gravity anomalies and disoriented aeromagnetic anomalies are in-line with the mobilistic system revealed by the palaeomagnetic data. These Alpine age continental rotations caused westward wrenching of the global lithosphere and led to significant tectonic reactivation and deformations. GPS measurements, current tectonic knowledge and the results of the evaluation of potential field data were combined in a base map to demonstrate similarities.     

Long-term forecasting of large volcanic eruptions

It is shown that the subareal volcanic system of the Earth evolves like individual volcanoes, preserving a constant average power. A simple model that allows the forecasting of extreme eruptions is proposed and discussed.     

Geochemistry and Ar/Ar geochronology of Miocene volcanic rocks from the Karaburun Peninsula: Implications for amphibole-bearing lithospheric mantle source,Western Anatolia

New geochemical and ⁴⁰Ar/³⁹Ar age data are presented from the Neogene volcanic units of the Karaburun Peninsula, the westernmost part of Western Anatolia. The volcanic rocks in the region are associated with Neogene lacustrine deposition and are characterized by (1) olivine-bearing basaltic-andesites to shoshonites (Karaburun volcanics), high-K calc-alkaline andesites, dacites and latites (Yaylaköy, Arma?anda? and Kocada? volcanics) of ~ 16–18 Ma, and (2) mildly-alkaline basalts (Ovac?k basalt) and rhyolites (Urla volcanics) of ~ 11–12 Ma. The first group of rocks is enriched in LILE and LREE with respect to the HREE and HFSE on N-MORB-normalised REE and multi-element spider diagrams. They are comparable geochemically with volcanic rocks in the surrounding regions such as Chios Island and other localities in Western Anatolia. The Ovac?k basalt is geochemically similar to the first stage early–middle Miocene volcanic rocks but differs from NW Anatolian late Miocene alkali basalts.     

K–Ar ages of the Akan-Shiretoko volcanic chain lying oblique to the Kurile trench: Implications for tectonic control of volcanism

The Akan‐Shiretoko volcanic chain, situated in the Southwestern Kurile arc, consists mainly of nine subaerial andesitic stratovolcanoes and three calderas. The chain extends in a SW–NE direction for 200 km, situated oblique to the Kurile trench at an angle of 25 degrees. Thirty‐seven new K–Ar ages, plus previous data, suggest that volcanic activity along the Akan‐Shiretoko volcanic chain began at ca 4 Ma at Akan, at the southwestern end of the chain, and systematically progressed northeastward, resulting in the southwest‐northeast‐trending volcanic chain. This spatial and temporal distribution of volcanoes can be explained by anticline development advancing northeastward from the Akan area, accompanied by magma rising through northeast‐trending fractures that developed along the anticlinal axis. The northeastward development of the anticline caused uplifting of the Akan‐Shiretoko area and changed the area from submarine to subaerial conditions. Anticline formation was likely due to deformation of the southwestern Kurile arc, with southwestward migration of the Kurile forearc sliver caused by oblique subduction of the Pacific plate. The echelon topographic arrangement of the Shiretoko, Kunashiri, Etorofu and Urup was formed at ca 1 Ma.     

Eocene volcanism during the incipient stage of Izu–Ogasawara Arc: Geology and petrology of the Mukojima Island Group,the Ogasawara Islands

The Ogasawara Islands mainly comprise Eocene volcanic strata formed when the Izu–Ogasawara–Mariana Arc began. We present the first detailed volcanic geology, petrography and geochemistry of the Mukojima Island Group, northernmost of the Ogasawara Islands, and show that the volcanic stratigraphy consists of arc tholeiitic rocks, ultra‐depleted boninite‐series rocks, and less‐depleted boninitic andesites, which are correlatable to the Maruberiwan, Asahiyama and Mikazukiyama Formations on the Chichijima Island Group to the south. On Chichijima, a short hiatus is identified between the Maruberiwan (boninite, bronzite andesite, and dacite) and Asahiyama Formation (quartz dacite and rhyolite). In contrast, these lithologies are interbedded on Nakodojima of the Mukojima Island Group. The stratigraphically lower portion of Mukojima is mainly composed of pillow lava, which is overlain by reworked volcaniclastic rocks in the middle, whereas the upper portion is dominated by pyroclastic rocks. This suggests that volcanic activity now preserved in the Mukojima Island Group records growth of one or more volcanoes, beginning with quiet extrusion of lava under relatively deep water followed by volcaniclastic deposition. These then changed into moderately explosive eruptions that took place in shallow water or above sea level. This is consistent with the uplift of the entire Ogasawara Ridge during the Eocene. Boninites from the Mukojima Island Group are divided into three types on the basis of geochemistry. Type 1 boninites have high SiO₂ (>57.0 wt.%) and Zr/Ti (>0.022) and are the most abundant type in both Mukojima and Chichijima Island Groups. Type 2 boninites have low SiO₂ (<57.1 wt.%) and Zr/Ti (<0.014). Type 3 boninites have 57.6–60.7 wt.% SiO₂ and are characterized by high CaO/Al₂O₃ (0.9–1.1). Both type 2 and 3 boninites are common on Mukojima but are rare in the Chichijima Island Group.     

Q c of three component seismograms of volcanic microearthquakes at Campi Flegrei volcanic area — Southern Italy

Digital recordings of three component microearthquake codas from shallow seismic events in the volcanic region of Campi Flegrei — Southern Italy — were used with an automatic technique to calculate the attenuation factorQ _c (codaQ) in the hypothesis of singleS toS backscattering.Results show the same value ofQ for each of the three components. This result is interpreted as due to isotropicS wave radiation pattern.A check of the coda method was performed using a single station method based on simple assumptions on the direct SH wave spectrum. Single stationQ was averaged over the stations and over the earthquakes. Results show that the two methods lead to comparable results.A frequency dependence quite different from that evaluated in active tectonic regions was found for coda attenuation, comparable to other volcanic areas throughout the world. This is interpreted as due to the presence of magma that affects anelasticity and scattering.     

Hazard assessment at volcanic fields: the Campi Flegrei case history

Volcanological analysis of the 10 000 yr –1538 explosive activity at Campi Flegrei shows that the most common explosive eruptions are characterized by the emplacement of flow or surge deposits, originating from the interaction between magma and shallow and/or sea water. The minimum volumes of pyroclastic products range between 0.04 and 0.7 km³; the proximal areas covered by these products range from 3–4 to 40–50 km². The pyroclastic flow and surge deposits occurring inside the caldera have been strongly controlled by pre-existent morphology; because of this, the area of present Napoli city was blanketed by approximately 5 m of pyroclastic deposits, during the last 5000 yr.Previous analysis suggests that the presence of even very low topographic obstacles may influence pyroclastic density current run out such that future eruptive deposits would mainly be confined inside the caldera rim. We suggest that a future eruption at Campi Flegrei would not seriously involve the urbanized area of Napoli city located on the hills. On the contrary, the plains located on the eastern side of the caldera (Fuorigrotta, Bagnoli) would be the most damaged area.     

Investigation of Crustal Thickness in Eastern Anatolia Using Gravity,Magnetic and Topographic Data

The tectonic regime of Eastern Anatolia is determined by the Arabia-Eurasia continent-continent collision. Several dynamic models have been proposed to characterize the collision zone and its geodynamic structure. In this study, change in crustal thickness has been investigated using gravity, magnetic and topographic data of the region. In the first stage, two-dimensional low-pass filter and upward analytical continuation techniques were applied to the Bouguer gravity data of the region to investigate the behavior of the regional gravity anomalies. Next the moving window power spectrum method was used, and changes in the probable structural depths from 38 to 52 km were determined. The changes in crustal thickness where free air gravity and magnetic data have inversely correlated and the type of the anomaly resources were investigated applying the Euler deconvolution method to Bouguer gravity data. The obtained depth values are consistent with the results obtained using the power spectrum method. It was determined that the types of anomaly resources are different in the west and east of the 40° E longitude. Finally, using the obtained findings from this study and seismic velocity models proposed for this region by previous studies, a probable two-dimensional crust model was constituted.