فهرست مطالب afsaneh nasrabadi

Crust velocity structure beneath two broadband seismic stations of Iran National Seismic Network (INSN), DAMV and THKV located in Central Alborz has been investigated by joint inversion of receiver function and Rayleigh wave phase and group velocity dispersion curves. The result suggests that Moho depth beneath the THKV and DAMV stations are 50-51 km and 52-54 km, respectively. Beneath the THKV station, there is a thin layer of very low-velocity materials at the surface and a sedimentary layer having a thickness of 10-12 km above a crystalline crust with a thickness of 34 km. Beneath the DAMV station, there is a thin sedimentary layer of low velocity with a thickness of 3-4 km, and also, a velocity change from 3.2 to 3.6 km/s at the depth of 14-16 km, indicating a discontinuity, which might be attributed to the border between the upper and lower crusts. The average Moho depth on the southern edge of Central Alborz is 52±2 km.

Crustal velocity structure in the Western Alborz have been investigated using joint inversion of receiver functions and Rayleigh wave group velocity dispersion curves. To determine the receiver functions, time domain iterative deconvolution and teleseismic events, which are recorded at five broadband seismic stations of the Iranian Seismological Center (IRSC), were used. The fundamental mode Rayleigh wave group velocity dispersion curves were provided by the study on the structure of crust and upper mantle of the Iranian Plateau. The results show that the average thickness of the crust in the southern margin of the Caspian Sea beneath of the CSN1 and RST1 stations is 38 km. Toward south, the depth of Moho increases up to 52, 54 and 50 km beneath the QALM, QCNT and QSDN stations. The low thickness of the crust in the southern shore of the Caspian Sea indicates that the Caspian crust is thin. Moreover, the moderate thickness of the crust in western Alborz, which is not in balance with its elevation, indicates the lack of root in Alborz.

Designing and launching rapid warning systems to reduce financial and human losses is one of the requirements of seismic areas. The main goal of these systems is to estimate the magnitude of the earthquake and the epicentral distance in a very short time and to announce a warning. A new method, called B-Δ, is used in these systems. It is based on information from a station and estimates the magnitude and epicentral distance of an earthquake using the initial seconds of the P wave. In this study, these two parameters have been estimated in east of Kerman Province after examining and processing of the vertical component of the accelerograms of the earthquakes in this area.

Crustal velocity structure and Moho discontinuity depth have investigated beneath 7 the broadband seismic stations, AFRZ, TKDS, TPRV, TNSJ, ANAR, KRSH of the Iranian Seismological Center (ISC) and YZKH of Iranian National Seismic Network (INSN) located in the center of Iran by joint inversion of receiver functions and Rayleigh waves group velocity dispersion. Three years (2012 to 2014) teleseismic waveforms (with epicentral distance 25o-90o) for computation receiver functions by iterative approach in time domain have been processed. The Rayleigh waves group velocity dispersion curves were incorporated into our joint inversion scheme from an independent surface wave tomography study. Receiver function is response of local structure of ground (located beneath the three–component broadband seismic station) to teleseismic P-wave, that is sensitive to seismic discontinuities. Since there is very little absolute-velocity information contained in the receiver function, its inversion for shear-wave velocity structure is non-unique (velocity-depth trade-off). On the other hand, dispersion curves are sensitive to the average velocity structure of the upper layers rather than to seismic discontinuities. So the non-uniqueness problem can be solved by combining receiver function inversion with surface-wave dispersion. Results from joint inversion in center of Iran indicates that Moho discontinuity depth depth beneath AFRZ, TKDS and TPRV stations is 40 Km, beneath TKDS 42 Km, beneath ANAR is 38 Km and beneath KRSH and YZKH stations are 44 Km. It was shown that the joint inversion method can cause ±2 kilometers of error. The average Moho depth is about 42±2 kilometers beneath center of Iran.

In this study, Seismic Hazard zoning was probed by studying seismogenic zones and earthquake data in a 150 km radius around the Fars province boundaries in the Zagros seismotectonic zone. Zagros is one of the most seismically active parts of Alpine-Himalayan seismic belt. Fars province located at latitude of 27.2° – 31.42° and longitude of 50.42° – 55.36°, respectively. In this study, seismotectonic zone was determined based on the location of interest, faults strikes, seismic records, and geological structure. In the first step of process, an earthquake catalogue of the instrumental and historical earthquakes was prepared by assimilating them to a uniform magnitude type (Mw) from 8th century until September 2016. Then, by using the Uhrhammer (1986) and Gardner and Knopoff (1974) methods in 1974, aftershocks and foreshocks were eliminated from the catalogue to have main earthquakes only. After the elimination of aftershocks and foreshocks, 2181 earthquakes remained as the main earthquakes in catalogue and their completeness magnitude was calculated by ZMap software. Afterwards, it was tried to recognize all the active faults around the sites; therefore, 23 area sources were introduced within the studied area from faults, historical and instrumental earthquakes. After introducing area sources, seismic parameters such as beta and lambda were calculated for each part using Kijko software. Similarly, space distribution function portion of each source was calculated. In addition, the results of six magnitude studies utilizing step method carried out by Mousvai et al. (2014) was used in this paper. In the deterministic approach, Pick Ground Acceleration (PGA) over bedrock was computed using two attenuation equations, Campbell and Bozorgnia (2008) and Boor and Atkinson (2008), using a MATLAB software for Fars province. Besides, in the probabilistic approach, PGA and pick ground velocity, PGV, were calculated using the same values in the attenuation equations that were used in the deterministic method using OpenQuake, one of the distinguished programs in seismic hazard analysis to determine the PGA on bedrock. The probabilistic approach was carried out on a grading area using a 10 km zonation map. PGA was calculated and covered 63% and 10% probability of exceedance in one life cycle of 50 and 475 years were presented, and hazard curves were calculated for 50 and 475 years using OpenQuake software. Probabilistic study resulted maximum PGA values in the periods of 50 and 475 years to be 0.23 g and 0.6 g respectively .These maximum PGA values centered on Kazeroon. After Kazeroon, Shiraz and Lamerd have the highest levels of hazard risk of cities in Fars province. Moreover, using the deterministic approach and the two attenuation equations, the obtained maximum PGA values fell between 0.68 g and 0.75 g. Finally, the results were compared with previous studies like Hamzeloo (2005), Ghafory-Ashtiany (1999) and Mousavi (2014).

Summary: Crustal velocity structure beneath four broadband seismic stations located in northeast of Iran, including Shahrood (SHRO) and Maraveh Tappeh (MRVT) stations set up by Iran National Seismic Network (INSN) and Sabzevar (SZBV) and Jarkhoshk (JRKH) stations set up by Iranian Seismology Center (IRSC), have been investigated by joint inversion of P receiver function and Rayleigh wave phase and group velocity dispersion curves. A three-year teleseismic data (2012 -2014) with epicentral distance of 25o-90o and magnitude more than 5.5 have been used to determine the receiver functions by iterative deconvolution in time domain proposed by Ligorria and Ammon (1999). Iterative deconvolution in time domain to determine the receiver functions are more stable with noisy data in comparison to frequency domain. The fundamental mode of Rayleigh wave group and phase velocity dispersion curves have been provided by the study on the structure of crust and upper mantle of the Iranian Plateau for the period interval of 10-100 seconds made by Rahimi (2010). A combined inversion of body wave receiver functions and Rayleigh wave velocities increases the uniqueness of the solution over the separate inversions, and also, facilitates explicit parameterization of the layer thickness in the model space. Moho discontinuity depth is one of the most important parameters for investigation of crustal structure. The results of this study indicate an average crustal thickness varying from 38 km beneath Maraveh Tappeh (MRVT) station in north of the study region up to 44 km beneath Shahrood (SHRO) station in west of the region. Moreover, the results of this study suggest that the average crust thickness beneath Sabzevar (SZBV) and Jarkhoshk (JRKH) stations located in center and east of the study region is 40 km. In general, northeastern Iran region has a thin crust compared to the crusts in the other areas investigated in this research work. It has also been shown that the joint inversion method can cause ±2 kilometers of error.
Introduction: Iran is situated in one of the world's seismic regions and the possibility of destructive earthquakes in most regions of the country has given great significance to recognition of Iranian seismic nature from a seismic and seismotectonic standpoint. The seismicity within Iran suggests that much of the deformation is concentrated in the Zagros, Alborz and Koppeh Dagh mountains, and in east of Iran, surrounding Central Iran and the Lut desert. The aim of this research is to study the crustal structure and Moho discontinuity of northeastern Iran region, Binalood mountains and Koppeh Dagh by the analysis of receiver function and surface waves dispersion.

Methodology and Approaches: Receivers functions are time series obtained from three-component seismometers, and are created by deconvolving the vertical component from the radial and transverse components of the seismogram to isolate the receiver site effects from the other information contained in a teleseismic P and S wave. The depth-velocity trade-off in receiver function causes nonuniqueness in the inverse problem. However, by incorporating information of absolute shear wave from dispersion estimates and joint inversion of these two datasets, this shortcoming can be compromised. To determine the receiver functions, we have used iterative deconvolution in time domain, proposed by Ligorria and Ammon (1999) that is more stable with noisy data in comparison to frequency domain. We have processed teleseismic events with epicentral distance of 25o-90o and magnitudes more than 5.5 that are recorded at a three-year time interval of 2012 to 2014. We have set the parameter a of the Gaussian filter to 1.00, which gives an effective high frequency limit of about 0.5 in the P wave. In order to eliminate the source, path and instrument effects, deconvolution of the vertical component from the horizontal components of the seismograms has been used. All receiver functions have been grouped by azimuth (<10◦) and distance (<15◦), and in order to improve the signal-to-noise ratio, the individual receiver functions within each group have been stacked. The fundamental mode of the Rayleigh wave group and phase velocities dispersion curves have been provided from the study carried out by Rahimi et al., (2014) on the structure of crust and upper mantle of the Iranian Plateau for the period interval of 10-100 seconds. Joint inversion of two independent data sets has been performed by considering appropriate weighting parameter obtained from Herrmann and Ammon program (2003). Minimizing standard error between real and predicted data is the criteria for getting the desired final and close to the earth real model. The inversion package requires that the real velocity structure is represented by a set of flat-lying, homogeneous, isotropic velocity layers. The starting model comprises of the layers having 1-km thick as the top 6 km of the model space, 2-km thick between the depths of 6 and 66 km, and 4 km thick between the depths of 66 and 78 km. The starting velocity for each layer in the model has been Vp=8.0 km/s, which equates to upper mantle velocity.

Results and Conclusions: The results of this study suggest that the average crust thickness beneath Shahrood (SHRO) station, located in west of the study region is 44 km and the average crust thickness beneath Sabzevar (SZBV) and Jarkhoshk (JRKH) stations located in center and east of the study region is 40 km. Furthermore, the crust thickness beneath Maraveh Tappeh (MRVT) station located in north of Koppeh Dagh region in is 38 km. In general, northeastern Iran region has a thin crust compared to the crust in other areas of northeast of Iran.

Today, due to the accumulation of resources in urban areas, urbanization has been increased and consequently, city limits have been increased. Thus, considering the population distribution in these areas as well as the existence of a large number of these areas in sedimentary region, the significance of the study of seismology and earthquake engineering for earthquake retrofitting and reduction of the risk of earthquakes has been increased. A phenomenon that enhances damage during the earthquake, even at great distances is the site effect. Mexico City earthquake, occurred in the epicentral distance of 300 km and caused damage, is a clear example of the enhanced damage caused by site effect. Thus, this study emphasizes on the site effect and underground structures affected by shear wave velocity in earthquakes.

Summary In this study, crustal velocity structure beneath two broadband seismic stations of Iran National Seismic Network (INSN), Ashtian-Arak (ASAO) and Naein (NASN) located in northwest of the Central Iran seimotectonic zone near the Ashtian and Nain cities have been investigated by joint inversion of P receiver function and of Rayleigh wave phase and group velocity dispersion curves. To determine the receiver functions, we have used iterative deconvolution in time domain proposed by Ligorria and Ammon (1999). which is more stable with noisy data in comparison to frequency domain. The fundamental mode Rayleigh wave group and phase velocities dispersion curves have been provided by the study of Rahimi et al. (2014) on the structure of crust and upper mantle of the Iranian Plateau for the period interval of 10-100 sec. The result of this study suggests that Moho discontinuity depth beneath Ashtian-Arak station (ASAO is 50 ± 2 km and beneath Naein station (NASN), it is 56 ± 2 km. Relative high crustal thickness beneath NASN station in comparison to other regions of central Iran can be attributed to abut the region to the Sanandaj–Sirjan zone (SSZ) and Urumieh– Dokhtar magmatic assemblage (UDMA). It can also attributed to the existence of thick Magma masses in Urumieh– Dokhtar magmatic assemblage and increase of the density and relative thickness of the area based on the isostasy theory. The average Moho depth in northwest edge of Central Iran is 53 ±2 km.
Introduction Iran is situated in one of the world seismic regions and the possibility of occurring destructive earthquakes in most regions of the country has given a great significance to recognition of Iranian seismic nature from a seismotectonic standpoint. The seismicity within Iran suggests that much of the deformation is concentrated in the Zagros, Alborz and Kopeh Dagh mountains, and also, in east Iran, surrounding Central Iran and the Lut desert, which are virtually aseismic and behave as relatively rigid . The aim of this research is the study of the crustal structure and Moho discontinuity of the northwest of the Central Iran from analysis of receiver function and surface waves dispersion.
Methodology and Approaches Receiver functions are the response of the local earth structure to the near-vertical arrival of p waves under a threecomponent seismogram and are susceptible to shear wave velocity contrasts. The depth-velocity trade-off in receiver function causes non-uniqueness in the inverse problem. However, by incorporating information of absolute shear wave from dispersion estimates and joint inversion of these two datasets, this shortcoming can be compromised. In this study, crustal velocity structure beneath two broadband seismic stations of INSN, i.e. ASAO and NASN located in northwest of the Central Iran seimotectonic zone near the Ashtian and Nain cities have been investigated by joint inversion of P receiver function and of Rayleigh wave phase and group velocity dispersion curves. To determine the receiver functions, we use iterative deconvolution in time domain proposed by Ligorria and Ammon (1999) which is more stable with noisy data in comparison with frequency domain and teleseismic events with source-receiver great circle paths larger than 30° and smaller than 90° with magnitudes more than 5.0 that are recorded at time period of 2009 to 2013. The 210 desired RFs have been recorded at two permanent stations. To remove high frequencies, Gaussian parameter 1.0 has been used. In order to eliminate the source, path and instrument effects, deconvolution of the vertical component from the horizontal components of the seismograms is used. For increasing signal to noise ratio, RFs have been clustered in 20˚ azimuthal and less than 15˚ epicentral distance ranges. Finally, the RFs are stacked. The fundamental mode Rayleigh wave group and phase velocities dispersion curves have been provided by the study of Rahimi et al. (2014) on the structure of crust and upper mantle of the Iranian Plateau for the period interval of 10-100 sec. In this way, more accurate information about the crustal structure can be obtained. Joint inversion of two independent data sets has been performed by considering combination of appropriate weighting parameter from Herrmann and Ammon program (2003). Minimizing standard error between real and predicted data is the criteria for getting to desired final and close to earth real model.
Results and Conclusions The result of this study suggests that Moho discontinuity depth beneath ASAO is 50 ± 2 km and beneath NASN is 56 ± 2 km. Relative high crustal thickness beneath NASN station in comparison to other regions of central Iran can be attributed to abut the region to the Sanandaj–Sirjan zone (SSZ) and Urumieh–Dokhtar magmatic assemblage (UDMA). It can also attributed to existence of thick Magma masses in Urumieh–Dokhtar magmatic assemblage and increase the density and relative thickness of the area based on the isostasy theory. The average Moho depth in northwest edge of Central Iran is 53 ± 2 km.

Iran is situated in one of the world's seismic regions and the possibility of destructive earthquakes in most regions of the country has given great significance to recognition of Iranian seismic nature from a seismic and seismotectonic standpoint. Study of the crust and upper mantle velocity structure in the Iranian plateau provides better understanding of its evolution and tectonic history of seismotectonic zones. Crustal velocity structure is used as initial information for various geological and geophysical studies, and therefore it is a basic and important issue in seismology. Receiver functions show Earth local structure response to P-wave vertical arrival approximately beneath of a three-component seismometer and are sensitive to shear-wave velocity impedance. Depth-velocity trade-off in RFs information is causing of inversion non-uniqueness problem, but one can overcome to this limitation by incorporating information from absolute velocity from dispersion estimations and joint inversion of this two data sets. By this, more exact constraints are provided about crustal structure. In this study, crustal velocity structure and Moho discontinuity depth beneath of four broadband stations of Kerman seismological network have been investigated from joint inversion of P-wave receiver functions (RFs) and Rayleigh wave group velocity dispersion. The teleseismic waveformes in time interval more than two years was used to compute RFs from the time domain iterative deconvolution procedure Ligorria and Ammon (1999) which has higher stability with noisy data compared to frequency-domain methods. The 165 desired RFs were computed from these waveforms that have magnitude bigger than 5.5 and have recorded at four permanent stations in epicentral distance 25˚-90˚. To delete high frequencies, Gaussian parameter 1.0 used. For increasing signal to noise ratio, RFs clustered in 10˚ azimuthal and less than 15˚ epicentral distance ranges. Finally, the RFs were stacked. This work performed under software SAC. Due to changes in group and phase velocity of surface waves with depth for different periods and dispersion in these waves and sensitivity of the waves dispersion curve to shear wave velocity, inversion of dispersion curve is an efficient method for determining the average shear wave velocity in a vast region of the depth between two seismic stations. Group velocity dispersion curves were incorporated into our joint-inversion scheme from an independent regional fundamental-mode Rayleigh waves tomography images for within the 20–80s period range in Iran by Rahimi et al. (2014). Joint inversion of two independent data sets was performed with considering combination weighting parameter appropriate performed from Herrmann and Ammon program (2003). Minimizing standard error between real and predicted data is the criteria for getting to desired final and close to earth real model.
The results from this study show that Moho discontinuity boundary is beneath of CHMN station at 52±2 km depth, beneath of KHGB station at 50±2 km depth, beneath of NGRK station at 54±2 km depth and beneath of TVBK station at 52±2 km depth. We used forward modeling test for error estimation and resulting models accuracy.
Relative high crustal thickness in this region compared to other regions of central Iran can be attributed to abut the region to the Sanandaj–Sirjan zone (SSZ) and Urumieh–Dokhtar magmatic assemblage (UDMA) that underthrusting of the Arabian plate beneath Central Iran along the main Zagros thrust fault is caused of thickening. It can also attributed to exist of thick Magma masses in Urumieh–Dokhtar magmatic assemblage and increase the density and relative thickness of the area based on the Isostasy theory.