Layout 6 ANNALS OF GEOPHYSICS, 60, SUPPLEMENT TO 6, SE668, 2017; doi: 10.4401/ag-7520 Attenuation characteristics, source parameters and site effects from inversion of S waves of the March 31, 2006 Silakhor aftershocks in Western Iran Somayeh Ahmadzadeh1, Stefano Parolai2,3, Gholam Javan-Doloei*,1, Adrien Oth4 1 International Institute of Earthquake Engineering and Seismology, Tehran, Iran 2 Helmholtz Centre Potsdam, GFZ German Research Centre for Geosciences, Potsdam, Germany 3 Now at: Istituto Nazionale di Oceanografia e di Geofisica, Sperimentale - OGS, Sgonico (TS), Italy 4 European Center for Geodynamics and Seismology, Luxemburg Article history Received August 15, 2017; accepted October 30, 2017. Subject classification: Generalized Inversion Technique; Source parameter; Site response; Bootstrap analysis; Zagros. SE668 ABSTRACT A nonparametric generalized inversion technique (GIT) has been used to derive attenuation characteristics, source pa- rameters and site amplification effects from S wave spectra of the 2006 Silakhor earthquake aftershocks recorded at 10 temporary stations in the epicentral area of the main shock along the main recent fault in Zagros region. We apply a total number of 971 spectral amplitudes from 213 after- shocks with ML 1.5-4.4 in the frequency range 0.5-35Hz. The obtained nonparametric attenuation curves decay monotonically with distance for the entire frequency range and the estimated S wave quality factor shows low Q values with fairly strong frequency dependence. We fit the omega- square source model to the inverted source spectra to derive earthquake source parameters. The estimated stress drops range from 0.04 to 5.3 MPa with an average of 0.9 MPa, within the range typically observed for crustal earthquakes in other regions of the World. The evaluated apparent stresses vary from 0.009 to 0.9 MPa and its relation with stress drops is generally consistent with the theoretically ex- pected relation assuming Brune model. The local magnitude ML shows a linear correlation with the moment magnitude MW that underestimates the MW scale by about 0.1-0.5 magnitude units for MW 1.5-3.0. We found a gradual in- crease in corner frequency values against seismic moments that may be considered as an evidence for a deviation from self-similarity scaling though the magnitude range of this study is fairly small to deduce a firm conclusion. The com- parison of estimated GIT site amplification functions with the evaluated horizontal to vertical (H/V) spectral ratios shows different levels of amplification in site responses due to amplification of the vertical component; however, the overall shape of site spectra provided by both methods are relatively similar. 1. Introduction The Iranian plateau among the Alpine-Hi- malayan seismic belt is limited between the conver- gent movements of the Arabian plate in the southwest and the Eurasian plate in the northeast. This conver- gence makes Iran one of the most seismically active areas of the World. Iran has been frequently shaken by great earthquakes accompanied with heavy casu- alties and economic losses. The Zagros fold - thrust belt, known as one of the youngest and most active continental collision zones on earth, is the result of the continued convergence between Arabia and Eura- sia plates [Snyder and Barazangi 1986]. Therefore, the Zagros region can be considered as one of the major active seismotectonic provinces with frequent and large earthquakes. The Main Recent fault (MRF) is an important element of tectonic units defined within the Zagros composed of a series of right-lateral strike- slip faults. Earthquakes occurring on or near the MRF have larger magnitudes than those in the other parts of the Zagros region. Also, the mechanism of earth- quakes and seismic potential of the MRF is quite dif- ferent from those earthquakes occurring within the rest of the Zagros range. The most devastating earth- quake of the region, the 1909 January 23 Silakhor earthquake (mb 7.4), occurred on the Main Recent fault and produced over 40 km surface rupture [Tchalenko and Braud 1974]. Based on the report of Ambraseys and Melville [1982], this historical event was strongly felt in a large area with more than 100 damaged villages and more than 6000 fatalities. The latest large event along the main recent fault (MRF), the Silakhor earthquake of 31 March 2006 with ML 6.1 (IIEES), occurred in the Silakhor plain in the south of Broujerd city in the Lorestan province in west Iran. A hypocentral depth ranging between 6 and 14 km was reported by several agencies [Peyret et al. 2008, Ghods et al. 2012]. According to field investiga- tions no significant surface rupture was observed and a maximum intensity of VIII was reported by Moghadam et al. [2006]. This event severely affected about 300 villages in the Silakhor plain with 68 fatali- ties and more than 1400 injuries. The 2006 Silakhor earthquake was preceded by two fairly large fore- shocks with magnitudes of 4.6 and 5.2, respectively, and followed by a large number of aftershocks. The occurrence of two relatively large foreshocks effec- tively alarmed many people to leave their homes lead- ing to a lower number of casualties than expected. Earthquake ground shaking at a specified loca- tion depends on three main region-specific parame- ters including seismic-source, wave propagation parameters and seismic-site effects. The study of these region-specific parameters is required to predict ac- curate ground-motion amplitudes in order to enable reliable seismic hazard assessment. The source, path and site effects can be separated from each other in the Fourier amplitude spectra (FAS) of recorded ground motion by implementing the common generalized inversion technique (GIT). This method was firstly proposed by Andrews [1986] recasting the spectral ratios method into a generalized inverse problem, then decomposing the body-wave spectra into source, site, and path contributions, si- multaneously. Several authors have used and devel- oped this approach to earthquake ground-motion data worldwide [Castro et al. 1990, Boatwright et al. 1991, Hartzell 1992, Parolai et al. 2000, Salazar et al. 2007 among others]. For instance, Castro et al. [1990] in- troduced a two-step nonparametric inversion method to isolate source, site and path effects along the Guer- rero, Mexico, subduction zone. In recent years, a few studies have been con- ducted to investigate seismic-source and wave propa- gation parameters in various regions of the Iranian plateau. Among these, however, only a limited num- ber are related to the Zagros region. Zafarani and Hassani [2013] applied the GIT to 148 S-wave ampli- tude spectra of 35 events recorded at 40 stations of the Iranian strong-motion network, which have been in- stalled by the Building and House Research Center (BHRC). Source parameters and site response of S- waves in two parts of the Zagros region were exam- ined for the events in the magnitude range from 4.2 to 6.2. Rahimi et al. [2008] estimated the coda Q, Qc, for the Zagros zone by analyzing the coda waves of 51 local earthquakes recorded on three stations of the Iranian National Seismic Network (INSN) with magni- tudes ranged from 3.1 to 4.9 during March and April 2006.Most of the analyzed events were foreshocks and aftershocks of the 2006 Silakhor earthquake. Using the SH-wave acceleration spectral data, Hamzehloo et al. [2010] estimated the source parameters and frequency dependent shear wave quality factor of the 2006 Silakhor earthquake and only one of its foreshocks and one aftershock that were recorded at 8 stations. In our study, we applied the nonparametric GIT to a large number of 2006 Silakhor earthquake after- shocks, with a dataset composed of 971 records from 213 events, to evaluate source, site and path effects, si- multaneously in the epicentral area of Silakhor earth- quake. The obtained nonparametric attenuation functions are parameterized to obtain an estimate of the S-wave frequency-dependent quality factor and the results are compared with other studies. The eval- uated source spectra are interpreted using the omega- square Brune [1970] model to retrieve seismic moments and stress drop values. Finally, we compare and discuss the site response functions estimated from GIT with those from horizontal-to-vertical (H/V) spectral ratios. 2. Database On 31 March 2006, the ML 6.1 (IIEES) Silakhor earthquake occurred along the Zagros main recent fault (MRF) in the Lorestan province, west Iran. The earthquake was followed by a large number of after- shocks. Most of them were recorded by a local tem- porary seismic network including 10 short-period stations that was deployed around the epicenter of the main shock, 5 days after the mainshock for about 2 months from 5 April to 10 June 2006. The network consisted of five three component short-period Gu- ralp CMG-6TD seismometers and five Guralp CMG- 5TD accelerometers with a sampling rate of 100 Hz. More than 1300 aftershocks were recorded during net- work installation. The local magnitudes of these events were determined using SEISAN software [Havskov and Ottemöller 1999]. We relocated events using the velocity model for the region [Sepahvand et al. 2012] and selected a subset of 213 events with well- constrained locations (location error <3 km and GAP <180°) with ML>1.5 (Figure 1). Each of the selected events is recorded by at least three stations and each of the used stations recorded at least three events. The final database is composed of 971 records from 213 earthquakes with a magnitude in range of ML 1.5-4.4 and hypocentral distances ranging from 3 km to 41 AHMADZADEH ET AL. 2 3 ATTENUATION, SOURCE AND SITE EFFECTS FROM INVERSION OF THE SILAkHOR AFTERSHOCkS km. Figure 1 shows the location of 2006 Silakhor mainshock and 213 selected aftershocks as well as temporary local stations around the Zagros Main Re- cent Fault with the source-station travel paths. The distribution of data in terms of magnitude versus hypocentral distance and focal depth is presented in Figure 2. The dataset spans the range of ML = 1.5-4.4 with focal depths ranging between 0.5-14 km. The majority of our events have magnitudes in the range 1.5 to 3 and the smallest hypocentral distance in the dataset is 3 km. Since near-field terms normally be- come most important at distances shorter than a few source dimensions, this effect is not expected to con- siderably affect the analysis. Figure 1. Epicenters of 2006 Silakhor earthquake (star) and aftershocks (circles). Inverse triangles indicate the location of 10 stations in- stalled around the epicenter of the main shock. The source-station ray paths are shown as gray lines. Figure 2. Distribution of earthquake magnitude (ML) versus hypocentral distance (left) and depth (right). For the spectral analysis, the raw recordings were corrected for instrument response and a time window of 5 seconds length was selected on the horizontal components starting 1 second before the S-wave ar- rival. Each selected window is cosine-tapered (5%) and the acceleration Fourier amplitude spectra (FAS) have been calculated. To compute the signal-to-noise ratio (SNR), a pre-event noise windows with the same AHMADZADEH ET AL. 4 Figure 3. Three instrument-corrected EW component velocity seismograms at stations ABSA, VANA and kABO (left) and calculated ac- celeration FAS of the seismograms (right). The smoothed FAS using the konno-Ohmachi window with b=20 are indicated with gray lines for both signal and noise spectra. 5 ATTENUATION, SOURCE AND SITE EFFECTS FROM INVERSION OF THE SILAkHOR AFTERSHOCkS length as the S-wave windows were considered (for some events the available pre-P noise were less than 5 seconds) and only records with a SNR larger than 3 at each frequency in the selected frequency range was taken into account for the present research. The selected frequency range of this study is lim- ited to 0.5-35 Hz, due to the signal-to-noise ratio (SNR) constraint. The obtained amplitude spectra were smoothed using the konno and Ohmachi [1998] windowing function with b=20 and sampled at 40 fre- quency points equidistant on log scale between 0.5 and 35 Hz. The root-mean-square average of the two horizontal components was considered in the analy- sis. The above data processing procedures are shown in Figure 3 for three examples. 3. Generalized Inversion Technique The generalized inversion technique (GIT) has been widely used to retrieve the source spectra, atten- uation characteristics and site effects in many regions of the world [e.g. Castro et al. 1990, 1996; Dutta et al. 2003; Oth et al. 2011]. Here we provide a brief de- scription of the GIT method and refer the reader to the above-mentioned papers for more in-depth details of this technique. The observed Fourier amplitude spectrum (FAS) of the ground motion can be expressed as multiplica- tion of the following three terms in frequency domain: Uij(f, Mi, Rij) = Si(f, Mi). A(f,Rij).Gj(f) (1) where Uij(f, Mi, Rij) is the observed FAS at site j from source i with magnitude Mi; Rij is the hypocentral dis- tance; Si (f, Mi) represents the source spectrum of the ith event; A (f, Rij) denotes the attenuation along the travel path, and Gj (f) is the site amplification function of the jth station. By taking the logarithm, the above Equation can be linearized as Equation 2: log10Uij(f, Mi, Rij) =log10 Si(f, Mi) +log10 A(f,Rij) + log10Gj(f) (2) Equation (2) provides a system of linear equa- tions of the form Ax=b, with data vector b including the observed spectral amplitudes, model vector x and system matrix A. This system of equations can be solved in a least squares sense for each frequency using appropriate inversion algorithms [e.g. SVD and LSQR] [Andrews 1986, Menke 1989]. According to Andrews [1986] one undetermined degree of freedom exists in Equation (2) that can be resolved by assum- ing a predefined source spectral model [Boatwright et al. 1991, Salazar et al. 2007] or setting the site response of one or more rock sites in the dataset equal to one [Andrews 1986, Castro et al. 1990, Parolai et al. 2000]. Following Castro et al. [1990] the inversion can be implemented using either a parametric or non- parametric approach. In parametric scheme [e.g. Salazar et al. 2007] a predefined function is considered for attenuation along the travel path while in the non- parametric approach [e.g. Oth et al. 2009] the attenuation is constrained to be a smooth function of distance including all attenuation effects like anelastic attenuation or geometrical spreading that takes the value of 1 at the reference distance R0, which is set to 9 km in present work. We solve the linear system of equations in one step implementing the LSQR algorithm [Paige and Saunders 1982] to determine the source spectra, site responses, and attenuation characteristics, simultane- ously. To eliminate the undetermined degree of freedom [Andrews 1986], one or several appropriate sites can be considered as a reference condition by set- ting either the amplification of one station or the average of several stations to be approximately one, irrespective of frequency. In order to select an appro- priate reference condition, the H/V spectral ratio and geological setting of all stations as well as the results of trial runs using various reference conditions have been considered. The H/V spectral ratio calculated from the amplitude spectra at two GOSH and VANA stations are almost flat and close to unity (Figure 4). On the other hand, these two stations are located on the rock outcrop based on field observations. For these reasons, in addition to the results of trial runs leading to reasonable spectral shapes of the source and site terms in the inversion, we conclude that selecting stations GOSH and VANA as reference stations and setting the average of site responses at these two stations to be equal to one is the most appropriate choice for this study. It is noteworthy that all site amplification and source spectra derived in the inversion are relative to the selected reference sites. The bootstrap method [Efron 1979)] is applied in order to evaluate the stability of inversion results fol- lowing Parolai et al. [2000, 2004]. In this method, a new dataset is obtained by means of randomly select- ing each row from the original system matrix with the same size. A certain row can either be selected several times or never. Repeated inversions of the bootstrap data set are performed at each frequency and then the means and standard deviations of these bootstrap samples assessed. We carried out 100 bootstrap inversions and calculated the mean and standard deviation of these 100 iterations for each model parameter. 4. Inversion results 4.1 Attenuation Characteristics Applying a one-step GIT inversion to the selected spectral amplitudes, the nonparametric attenuation functions are determined. The hypocentral distance ranging from 3km to 41km is divided into 13 bins of 3 km width. We acknowledge that with this limited dis- tance range available in the dataset, only limited conclusions can be drawn on the overall attenuation characteristics and in particular its frequency depen- dence. Nonetheless, as shown below, the attenuation curves are robust and the data coverage within the limited distance range is excellent. The obtained attenuation functions for two selected frequencies together with the observed amplitudes corrected for source and site terms are shown in Figure 5. The curves decay monotonically with distance and the observed data points are fairly described by the calculated attenuation curves. The attenuation functions for all 40 frequency points selected for the analysis between 0.5 Hz and 35 Hz are displayed in Figure 6. A 1/R decay function is also indicated for comparison. As we set the reference distance to 9 km, the attenuation function takes the value log A(f,R)=0 at 9 km. In general, the curves show a simple shape con- sistently decreasing with distance almost for the entire frequency range with a more rapid decay for the high frequency attenuation curves than the low frequency ones, as expected. The low frequency curves decay similarly to a 1/R function up to about 20km whereas at larger distances the curves attenuate more quickly than 1/R for the entire frequency band. AHMADZADEH ET AL. 6 Figure 4. The H/V spectral ratio calculated from the amplitude spectra at station GOSH (a) and VANA (b). The spectral ratios are almost flat and near to one. Figure 5. Nonparametric attenuation functions for two selected frequencies. The gray shaded area denotes one standard deviation cal- culated by bootstrap analysis around the mean. The observed spectral amplitudes corrected for source and site terms are indicated as gray crosses. 7 The obtained nonparametric attenuation term con- tains all attenuation characteristics including apparent quality factor and geometrical spreading. Assuming a simple geometrical spreading model, the apparent qual- ity factor can be estimated through fitting the following equation to the attenuation function: (3) where R is the hypocentral distance in km, C is the ref- erence distance for the fit, Q(f ) denotes the apparent S- wave quality factor, vS is the shear-wave velocity and b is the geometrical spreading exponent. Due to the existing trade off in estimation of b and Q(f ) together, we set the b=1 and restrict our calculation to the definition of Q(f ). The nonparametric attenuation function, A (f,R) is cor- rected for the geometrical spreading (G(r) = C/R) and then the frequency dependent quality factor is evaluated from the slope of a least squares regression to log A(f,r) - logG(r) versus distance. We fit a power function of the form Q(f) = Q0f N to the obtained Q(f ) values through a least square fit to determine Q0 and N. As shown in Fig- ure 7, the frequency dependent quality factor is approx- imated with the relation: Q(f)=35.99 f 0.77 (4) The defined apparent quality factor with N=0.77 demonstrates fairly strong frequency dependence in our studied area located in Zagros region. This evaluated ap- parent Q factor with low Q0 and value of N near to 1 in- dicates high attenuation characteristics in this region and is similar to the results have been found for other tec- tonically active regions of Iran and the world [Rovelli 1982, Hellweg et al. 1995, Yoshimoto et al. 1993, Cas- tro et al. 1996, Bindi et al. 2006, Motazedian 2006 amongst others]. It is important to point out that due to the trade- off between Q and the geometrical spreading function, we discuss Q in the framework of the simplified geo- metrical spreading model adopted here. In addition since the nonparametric attenuation function (A(f,r)) includes all effects leading to attenuation of seismic waves along the travel path, the apparent quality factor obtained by fitting Equation (3) to the attenuation curves include both intrinsic and scattering attenuation. However obtained apparent quality factor is frequency dependent, we cannot definitely conclude that the in- trinsic Q is also frequency dependent. 4.2 Source Spectra We fit the ω2 source model [Brune 1970, 1971] to the inverted nonparametric source spectra in order to retrieve earthquake source parameters including corner frequency, fC and seismic moment, M0. The acceleration source spectra in terms of the ω2 model can be expressed by: (5) and (6) where Rθφ is the average shear wave radiation pat- tern, V=1/√2 accounts for the division of S-wave en- ergy into two horizontal components, F=2 is free sur- ATTENUATION, SOURCE AND SITE EFFECTS FROM INVERSION OF THE SILAkHOR AFTERSHOCkS ! !!! ! ! !!! !"# !!"!! ! !!! ! !! ! Figure 6. Nonparametric attenuation functions versus hypocen- tral distance for all frequencies. For comparison, a 1/R curve is depicted as a dashed line. Note that logA(f,R)=0 at the reference distance R0=9 km. Figure 7. Solid black line indicates the obtained apparent Q(f ) model from the analysis of the attenuation functions shown in Figure 6. The individual Q values derived at each given frequency are denoted by dots and error bars. S( f )=(2! f )2 RθϕVF 4!ρυS 3R0 M( f ) M( f )= M0 1+( f / fC ) 2 face amplification, ρ = 2.8 g/cm3 and VS=3.4 km/s are density and shear wave velocity estimates in the vicin- ity of source region, respectively, R0 denotes the refer- ence distance, and M(f) represents the moment rate spectrum. Figure 8 shows four examples of inverted source spectra (solid black lines, with standard deviation as gray shading) along with an ω2 source model fitted to the spectra (dashed lines) using non-linear least squares. Estimated corner frequencies and seismic moments are defined by this fitting procedure. Taking into ac- count the standard deviations obtained from 100 boot- strap analyses, the evaluated source spectra are gener- ally very stable and the spectra are appropriately fitted by ω2 model. Under the assumption of the ω2 model the accel- eration spectra increase as f 2 at low frequencies with a plateau at high frequencies. Generally, the observed ac- celeration Fourier spectra show a decay at high fre- quencies above a certain frequency which is character- ized by an exponential of the form exp(-πκf ) [Anderson and Hough 1984]. The high-frequency spectral decay (κ) is often considered as a site effect originating from near surface attenuation [Hanks 1982, Anderson and Hough 1984] or a source related parameter [Papageorgiou and Aki 1983]. Since we did not consider this high-frequency κ decay term in the reference site constraint (stations GOSH and VANA), it can be reasonably expected to shift into the estimated source spectra as indicated by Oth et al. [2011]. In order to evaluate the influence of the exponential κ term on the source spectra we fitted equa- tion (7) to the high-frequency part of the spectra to es- timate the κ values Oth et al. [2011] as follows (7) where f E is the lower frequency limit. Taking into account the magnitude range of the events, we set f E = 12 Hz (we checked results from ap- plying both f E = 12 and f E = 15 Hz without remark- able differences) to be adequately above the corner fre- quency range of the events [Parolai and Bindi 2004] to obtain reliable κ estimates, however essentially f E would have to be properly selected for each individual spectrum. AHMADZADEH ET AL. 8 Figure 8. Examples of acceleration source spectra derived from GIT inversion (solid black lines) with related standard deviation (shaded area) of four events. The best fitted ω2 model is represented by dashed lines. The white squares trend indicates the high frequency κ cor- rection at frequencies larger than 12 Hz. 9 The distribution of κ values is depicted in Figure 9. The obtained values are normally distributed and very small, with the mean value of 0.011 s. These small κ values indicate the attenuation corrected spec- tra are subjected to relatively slight high frequency diminution effect. To investigate how these κ values affect the source spectra, the spectra have been cor- rected for κ decay for frequencies larger than 12 Hz (indicated with white squares at high frequencies in Figure 8). As can be seen from the Figure 8, the κ ef- fect slightly affects the high frequency part of source spectra in our dataset. We determine the stress drop values based on the estimated M0 and fC values using the following rela- tions [Brune 1970, 1971]: (8) (9) Moreover, we compute the apparent stress (σa) for each event as follows: (10) where μ is the rigidity modulus (μ=vs 2.ρ) and Es is the radiated S-wave energy that can be calculated from the inverted source spectra using the relation [e.g. Izutani and kanamori 2001]: (11) Figure 10a represents the evaluated corner fre- quency versus seismic moment with the correspond- ing standard deviations indicated as horizontal and vertical error bars, respectively. Dashed lines represent constant stress drop values from 0.01 to 10 MPa. The inferred corner frequencies range vary from 1.7 to 17 Hz whereas the seismic moments vary within 4 orders of magnitude, from 10^11 Nm to 10^15Nm. We found a less rapid decrease of corner frequency with increasing seismic moment than expected from self- similar scaling considerations, as can be seen in Fig- ure 10a. This behavior may be considered as an evidence for a break in self similarity scaling as re- ported in some studies [Malagnini et al. 2008, Mayeda et al. 2007, Mayeda 2009, 2010] though the magnitude range of our data (ML 1.5-4.4) with the vast majority of earthquakes (about 90 percent) with ML 1.5-3 is fairly limited to allow us to deduce a firm conclusion about non self-similar scaling. The evaluated stress drops range from about 0.04 to 5.2 MPa, with an average of 0.9 MPa (Figure 10b). These values are close to the stress drops reported by Hamzehloo et al. [2010] for three events of the Silakhor earthquake that are not included in our database. They obtained a stress drop value of 4.3 MPa for the mainshock (MW 6.1) and 3.2 MPa and 2.8 MPa for the selected foreshock and largest aftershock with MW = 5.1 and 4.9, respectively. The obtained apparent stresses range from about 0.016 to 2.5 MPa, with an average of 0.2 MPa. A com- parison of the evaluated apparent stress and stress drop values is represented in Figure 11. The solid line indicates the theoretical ratio of ∆σ⁄σa =4.3 [Singh and Ordaz 1994]. Considering the errors in Figure 11, the relation between estimated apparent stress and stress drop values is approximately consistent with the the- oretically expected relation. The comparison of estimated MW versus local magnitude ML is depicted in Figure 12. The moment magnitude MW and the local magnitude ML are linearly correlated. The local magnitude ML generally under- estimates MW by about 0.1 - 0.5 magnitude units for MW 1.5-3.0 and larger events have approximately the same MW and ML. The relation between the estimated moment magnitudes MW and local magnitude ML can be expressed using least squares regression as follows Mw = (0.68 ± 0.03) ML + (0.91 ± 0.07). (12) There are only a few studies that have investi- gated magnitude conversion relations for Iran [Shah- var et al. 2013, Zare et al. 2014, Shoja-Taheri et al. 2007]. Considering the magnitude range of events, to our knowledge this study is the first one looking at the relation between MW and ML within the magnitude ATTENUATION, SOURCE AND SITE EFFECTS FROM INVERSION OF THE SILAkHOR AFTERSHOCkS r = 2.34vS 2! fC !σ = 7M0 16r3 σa =µ ES M0 ES = 4! 5ρvS f  M( f ) 0 ! " 2  df Figure 9. Distribution of the estimated κ values evaluated from the high-frequency part of the inverted source spectra for f >12 Hz. The obtained κ values have approximately normal distribu- tion. range ML 1.5-4.4 in Iran and is in good agreement with previous studies in Iran that investigated different magnitude ranges [e.g. Shahvar et al. 2013 and Zare et al. 2014 relationships that are shown in Figure 12]. 4.3 Site Response The GIT site amplification terms have been de- rived for both horizontal (H) and vertical component (Z) data for all stations. We performed the GIT inver- sion for the vertical (Z) component applying the same S-wave windows and reference condition as the hori- zontal (H) component to investigate the Z component amplification effect. We provide a comparison between the GIT site response functions and the horizontal to vertical (H/V) spectral ratios calculated directly from the observed spectra. The H/V method was first ap- plied by Lermo and Chavez-Garcia (1993) to the S- waves of earthquake recordings to evaluate site effects. This technique assumes that the vertical component of the ground motion is not considerably affected by the local site conditions. Figure 13 depicts the evaluated GIT site amplifica- tion function ratio of the horizontal and vertical com- ponents (GIT H/Z) (solid line) and the H/V spectral ra- tios (dashed line) as well as the individual GIT amplification curves for horizontal (GIT H) and vertical (GIT Z) components for all stations. As previously men- tioned we considered stations GOSH and VANA as ref- erence constraint for which the H/V spectral ratios are reasonably flat and near to one over the entire fre- quency range of analysis (as can be seen from the Fig- ure 13). The amplification functions are generally well constrained with small standard deviations obtained from the bootstrap analysis indicated in the GIT H curves as gray shaded areas with the maximum ampli- fication of about 20 with respect to the reference sites. Overall the observed GIT H/Z and H/V curves are similar with small discrepancies. These small variations most probably originate from the differences in the ref- erence site's amplification function, as we set the aver- age response functions of stations GOSH and VANA to 1 for GIT method though it is not exactly one and slightly larger than unity for H/V case. Therefore H/V ratios underestimate the GIT H/Z results in most cases. AHMADZADEH ET AL. 10 Figure 11. Stress drop versus apparent stress with error bars. The theoretical relation is indicated by solid line. Figure 10. (a) Corner frequency versus seismic moment plot with constant stress drop values indicated by dashed lines. (b) Stress drop versus moment magnitude. Stress drop values range from 0.04 to 5.2 MPa. Figure 12. Local magnitude ML versus moment magnitude MW es- timated from source spectral fitting. The solid line shows the best least squares fit. The dashed line indicates a one-to-one scaling. 11 The largest amplification based on GIT H/Z is ob- served at station CHAG with a factor of 7 with peak near 1 Hz. The amplification peaks almost range from 1 to about 8 Hz for the majority of stations. The GIT H/Z and H/V site functions represent significant differences with GIT H curves in the level of amplification. These differences are due to amplifica- tion of vertical component as can be seen for instance at stations LENJ and CHAG at frequencies 10 and 25 Hz respectively. Based on geological maps the five sta- tions (i.e. VANA, DINA, GOSH, DEHT and kABO) are deployed on bedrock, for which the trends of the GIT H and H/V ratio are most consistent. The small differ- ences between these stations in the high frequency ATTENUATION, SOURCE AND SITE EFFECTS FROM INVERSION OF THE SILAkHOR AFTERSHOCkS Figure 13. Site amplification functions resulting from the GIT in comparison with H/V spectral ratios (dashed lines) directly calculated from the observed ground-motion spectra for all stations. The GIT H/Z (solid lines) denotes the GIT site response of H component di- vided by GIT site response of Z. The individual GIT site response functions for horizontal (solid lines with circles) and vertical (solid lines with squares) components are also indicated. The gray shaded areas mark the mean ± one standard deviation from 100 bootstrap analy- ses for the GIT H site responses. range are most probably due to the complex geological setting (e.g. sedimentary rocks; dolomite limestone and igneous rocks; granite and grano-diorite) around sta- tions. The remaining stations (i.e. ABSA, CHAL, LENJ, ZAGH and CHAG) are mostly placed on a thin layer of quaternary alluvium. The observed decrease in the ratio of H/V and GIT H/Z caused by amplification of vertical component in the high frequency range at sta- tions CHAG, LENJ and ZAGH is most probably due to a sharp contrast interface in this area and can be related to S-P conversions at the bottom of the soft layer in the case of high-impedance contrasts as Parolai and Rich- walski [2004] discussed. 5. Discussion and conclusions We applied the Generalized Inversion Technique (GIT) to the data recorded by 10 temporary stations installed around the epicenter of 2006 Silakhor earth- quake for 2 months. We investigated the source pa- rameters, attenuation characteristics and site response functions using the non-parametric one step GIT to the observed spectral amplitudes of 213 events within the magnitude range ML 1.5-4.4. The attenuation curves decay monotonically with distance for the entire frequency range with a more rapid decay for the high frequency attenuation curves than the low frequency ones. The attenuation is almost consistent with 1/R de- cay function for the low frequency curves up to about 20 km however beyond this distance the curves atten- uate more quickly than 1/R for all frequency range. A comparison of the obtained apparent quality factor in this study with other studies for Iran is de- scribed in Figure 14 (see also Table 1). Hamzehloo et al. [2010] applied the SH-wave acceleration spectral data of the 2006 Silakhor earthquake and only one of its foreshock and aftershock to estimate the shear wave quality factor, QS(f) = 121f 0.55 in the six selected frequency bands from 1 to 24 Hz. Hassani et al. (2011) estimated the S-wave quality factor as QS(f)=151f 0.75 using GIT with a dataset consisting of 40 earthquakes recorded in East-Central Iran for the frequency range 0.4-15 Hz. Farrokhi and Hamzehloo (2016) investi- gated the attenuation of P and S waves in Alborz and north central part of Iran. They applied the extended coda normalization method (CNM) for 380 local earthquakes in the seven frequency bands between 0.4 to 24 Hz to define QS(f)=83 f 0.99 and QS(f)= 68f 0.96 in Alborz region and North Central Iran, respectively. Another quality factor relation in North of Iran is re- ported by Motazedian [2006] as QS(f)=87f 1.46 based on the vertical component of 22 earthquakes. The es- timated S wave quality factor in our study shows fairly strong frequency dependence and high attenuation in seismically active region of Zagros which is compati- ble with results reported for other parts of Iran. How- ever, considering slight differences between these relations, it should be noted that in these studies dif- ferent methods have been applied to obtain quality factors implementing various geometrical spreading functions in different frequency ranges sampling dif- ferent crustal volumes. Furthermore, the study areas of these models are different (except Hamzehloo et al., 2010) and each relation is linked to various seis- motectonic provinces of Iran with distinct seismotec- tonic characteristics [e.g. Mirzaei et al. 1998]. The inverted source spectra were fitted by the omega square source model satisfactorily to retrieve earthquake source parameters. The inferred corner frequencies are ranging between about 1 and 17 Hz spanning within 4 orders of magnitude. We found a less rapid decrease of corner frequency with increas- ing seismic moment than expected from self-similar scaling considerations. This result can be considered as an evidence for a break in self-similar scaling, how- ever bearing in mind the small magnitude range of this study (ML 1.5-4.4) including the vast majority of earthquakes with ML 1.5-3 (about 90 percent), more work is needed to derive a definitive conclusion about non self-similarity in this sequence. The estimated stress drop values range between 0.04 and 5.2 MPa with an average of 0.9 MPa. The ob- tained rather low stress drops are near the lower limit of the typical worldwide observations [~0.1-100, e.g., kanamori 1994]. The values are however in the same range as observed for crustal earthquakes in Japan [Oth 2013] and are in agreement with the expected smaller stress drops for interplate earthquakes (such AHMADZADEH ET AL. 12 Figure 14. The comparison of the QS(f ) estimated in this study (solid black line) with other studies for Iran. 13 as Silakour events within tectonically active Zagros area) compared to higher stress drop values observed for intraplate events [e.g., Allmann and Shearer 2009]. Our estimated stress drop values were compared with the evaluated apparent stresses. The relation be- tween these two parameters is generally consistent with the theoretically expected relation assuming Brune model. Comparison of evaluated MW versus local mag- nitude ML indicates that the local magnitude ML gen- erally underestimates the moment magnitude MW in the order of 0.1-0.5 magnitude units. This result is in good agreement with previous studies in Iran [e.g. Shahvar et al. 2013 and Zare et al. 2014]. Considering the magnitude range of events, our study is the first one that investigated the relation between MW and ML within the magnitude range ML 1.5-4.4 in Iran. To investigate the influence of kappa parameter on the source spectra, this parameter was estimated from the high-frequency part of the spectra. The eval- uated κ values with small mean value of 0.011 sec in- dicate the attenuation corrected spectra are subjected to small high frequency diminution effect. The results of site amplification functions ob- tained from GIT, represent the amplification peaks, al- most range from 1 to about 8 Hz for the majority of stations with the maximum amplification of about 20 with respect to the reference sites. The different site responses for stations depend on the geological situa- tion beneath the seismic stations, individually. Some stations are not strongly affected by amplification ef- fects and are supposed to be located on rock. The pro- vided GIT Z site responses for the vertical components present considerable amplifications lead- ing to reduction in the H/V amplitude in comparison with the GIT H amplification functions. However, the overall shape of site spectra estimated by both meth- ods is almost similar. Finally, the results of our study including the re- gion-specific source, attenuation functions and site re- sponses provide important elements for further strong motion simulations and seismic hazard assessment. Acknowledgements. The authors acknowledge Inter- national Institute of Earthquake Engineering and Seismology (IIEES) for providing earthquake waveform database used in this study. This research was performed during the 7 months visit of S. Ahmadzadeh at the German Research Centre for Geosciences (GFZ) in Potsdam that was partially funded by the GFZ. We also appreciate the anonymous reviewers which their comments and suggestions improved our paper. References Allmann, B. P. and P. M. Shearer (2009). Global varia- tions of stress drop for moderate to large earth- quakes, J. Geophys. Res., 114, B01310, doi:10.1029/2008JB005821. Ambraseys, N. N. and C.P. Melville (1982). A History of Persian Earthquakes, Cambridge Univ. Press, Cambridge, 219 pp. Anderson, J. G. and S. E. Hough (1984). A model for the shape of the Fourier amplitude spectrum of ac- celeration at high frequencies, Bull. Seism. Soc. Am.,74(5), 1969-1993. Andrews, D. J. (1986). Objective determination of source parameters and similarity of earthquakes of different size, in Earthquake Source Mechanics, edited by S. Das, J. Boatwright, and C. H. Scholz, pp. 259-268, AGU, Washington, D. C. Bindi, D., S. Parolai, H. Grosser, C. Milkereit, and S. karakisa (2006). Crustal attenuation characteristics in northwestern Turkey in the range from 1 to 10 Hz, Bull. Seism. Soc. Am., 96, 200-214. Boatwright, J., J. Fletcher, and T. Fumal (1991). A gen- eral inversion scheme for source, site, and propa- gation characteristics using multiply recorded sets of moderate-sized earthquakes, Bull. Seism. Soc. Am., 81, 1764-1782. ATTENUATION, SOURCE AND SITE EFFECTS FROM INVERSION OF THE SILAkHOR AFTERSHOCkS Region Frequency range (Hz) Quality factor Geometrical spreading Reference Zagros region Six frequency bands between 1≤f ≤24 QS(f )=121f 0.55 R-1 Hamzehloo et al. [2010] East-Central of Iran 0.4≤f ≤15 QS(f )=151f 0.75 R -1 , R<60 (60 R)-0.5 , R≥60 Hassani et al. [2011] Alborz region North Central Iran Seven frequency bands between 0.4≤f ≤24 QS(f )=83f 0.99 QS(f )=68f 0.96 R-1 , R<90 (90 R)-0.5 , R≥90 Farrokhi and Hamzehloo. [2016] North of Iran 1150 Motazedian. [2006] Zagros region 0.5≤f ≤35 QS(f )=36f 0.77 R-1 This study Table 1. S-wave attenuation models reported for Iran. Brune, J. N. (1970). Tectonic stress and the spectra of seismic shear waves from earthquakes, J. Geophys. Res., 75, 4997-5009. Brune, J. N. (1971). Correction, J. Geophys. Res., 76, 5002. Castro, R. R., J. G. Anderson, and S. k. Singh (1990). Site response, attenuation and source spectra of S waves along the Guerrero, Me´xico, subduction zone, Bull. Seism. Soc. Am., 80, 1481 - 1503. Castro, R. R., F. Pacor, A. Sala, and C. Petrungaro (1996). S wave attenuation and site effects in the re- gion of Friuli, Italy, J. Geophys. Res., 101, 22355- 22369. Dutta, U., N. Biswas, A. Martirosyan, A. Papageor- giou, and S. kinoshita(2003). Estimation of earth- quake source parameters and site response in Anchorage, Alaska from strong-motion network data using generalized inversion method. Phys. Earth Planet. Inter., 137, 13-29. Efron, B. (1979). Bootstrap methods, another look at the jacknife, Ann. Stat., 7, 1-26. Farrokhi, M. and H. Hamzehloo (2016). Body wave at- tenuation characteristics in the crust of Alborz re- gion and North Central Iran, J. Seismol., doi:10.1007/s10950-016-9624-2. Ghods, A., M. Rezapour, E. Bergman, G. Mortezane- jad, and M. Talebian (2012). Relocation of the 2006 Mw 6.1 Silakhour, Iran, earthquake sequence: De- tails of fault segmentation on the main recent fault, Bull. Seism. Soc. Am., 102(1), 398-416. Hamzehloo, H., H. Rahimi, I. Sarkar, M. Mahood, H. MirzaeiAlavijeh, and E. Farzanegan (2010). Model- ing the strong ground motion and rupture charac- teristics of the March 31, 2006, Darb-e-Astane earthquake, Iran, using a hybrid of near-field SH- wave and empirical Green’s function method, J. Seismol., 14, 169-195. Hartzell, S. H. (1992). Site response estimation from earthquake data, Bull. Seism. Soc. Am., 82, 2308 - 2327. Hassani, B., H. Zafarani,J. Farjoodi, and A. Ansari (2011). Estimation of site amplification, attenua- tion and source spectra of S-waves in the East-Cen- tral Iran, Soil Dyn Earthq Eng 31, 1397-1413. Hellweg, M., P. Spandich, J. B. Fletcher, and L. M. Baker (1995). Stability of coda Q in the region of Parkfield, California: view from the U.S. Geologi- cal survey Parkfield dense seismograph array, J. Geopys. Res., 100, 2089-2102. Havskov, J. and L. Ottemoller (1999). SeisAn earth- quake analysis software, Seismological Research Letters, 70, 5, 532-534. kanamori, H. (1994). Mechanics of earthquakes, Annu Rev Earth Planet Sci., 22, 207-237. konno, k., and T. Ohmachi (1998). Ground-motion characteristics estimated from spectral ratio between horizontal and vertical components of mi- crotremors, Bull. Seism. Soc. Am., 88, 1228-1241. Lermo, J. and F. J. Cha´vez-Garcı´a (1993). Site effect evaluation using spectral ratios with only one sta- tion, Bull. Seism. Soc. Am., 83, 1574-1594. Malagnini, L., L. Scognamiglio, A. Mercuri, A. Akinci, and k. Mayeda (2008). Strong evidence for non-sim- ilar earthquake source scaling in central Italy, Geo- phys. Res. Lett., 35, L17303. Mayeda, k., L. Malagnini, and W. R. Walter (2007). A new spectral ratio method using narrow band coda envelopes: Evidence for non-self-similarity in the Hector Mine sequence, Geophys. Res. Lett., 34, L11303. Mayeda, k., and L. Malagnini (2009). Apparent stress and corner frequency variations in the 1999 Taiwan (Chi-Chi) sequence: evidence for a step-wise increase at Mw < 5.5, Geophys. Res. Lett., 36, L10308. Mayeda, k., and L. Malagnini (2010). Source radiation invariant property of local and near-regional shear- wave coda: Application to source scaling for th Mw 5.9 Wells, Nevada sequence. Geophys. Res. Lett. 37, L07306. Menke, W. (1989). Geophysical data analysis: discrete in- verse theory, in Int. Geophys. Series, R. Dmowska and J. R. Holton (Series Editors), Vol. 45, Academic Press, New York, 289 pp. Mirzaei, N., G. Mengtan and C, Yuntai (1998). Seismic source regionalization for seismic zoning of Iran: Major seismo tectonic provinces, J Earthquake Pred Res, 7, 465-95. Moghadam, A. S., k. Hesami,G, Javan Doloei, M. Mah- davifar and H. Hamzehloo (2006). The 2006 Darb-e- Astane earthquake field report. In: International in- stitute of earthquake engineering and seismology, 233 p. (In Farsi). Motazedian, D. (2006). Region-specific key seismic pa- rameters for earthquakes in northern Iran. Bull. Seism. Soc. Am., 96, 1383-1395. Oth, A., S. Parolai, D. Bindi and F. Wenzel (2009). Source spectra and site response from Swaves of in- termediate-depth Vrancea, Romania, earthquakes, Bull. Seism. Soc. Am., 99, 235 254. Oth, A., D. Bindi, S. Parolai and D. Di Giacomo (2011). Spectral analysis of k-NET and kik-net data in Japan, Part II: On attenuation characteristics, source spectra, and site response of borehole and surface stations. Bull. Seism. Soc. Am., 101, 667-687. AHMADZADEH ET AL. 14 15 Oth, A. (2013). On the characteristics of earthquake stress release variations in Japan, Earth Planet. Sci., 377-378, 132-141. Paige, C. C. and M. A. Saunders (1982). LSQR: An al- gorithm for sparse linear equations and sparse least squares, ACM Trans. Math. Software 8, 43-71. Parolai, S., D. Bindi and P. Augliera (2000). Application of the generalized inversion technique (GIT) to a microzonation study: numerical simulations and comparison with different site-estimation tech- niques. Bull. Seism. Soc. Am., 90, 286-297. Parolai, S., and D. Bindi (2004). Influence of soil-layer properties on k evaluation, Bull. Seism. Soc. Am., 94, 349-356. Parolai, S. and S. M. Richwalski (2004). The impor- tance of converted waves in comparing H/V and RSM site response estimates, Bull. Seism. Soc. Am., 94, 304-313. Parolai, S., D. Bindi, M. Baumbach, H. Grosser, C. Milkereit, S. karakisa and S. Zunbul (2004). Com- parison of different site response estimation tech- niques using aftershocks of the 1999 Izmit earth- quake, Bull. Seism. Soc. Am., 94, 1096. Peyret, M., F. Rolandrone, S. Dominguez, Y. Djamour and B. Meyer (2008). Source model for the Mw 6.1, 31 March 2006, Chalan-Chulan earthquake (Iran) from InSAR, Terra Nova20, 126-133. Rahimi, H. and H. Hamzehloo (2008). Lapse time and frequency-dependent attenuation of coda waves in the Zagros continental collision zone in South- western Iran, J. Geophys. Eng. 5, 173-185. Rovelli, A., O. Bonamassa, M. Cocco, M. Di Bona and S. Mazza (1988). Scaling laws and spectral parame- ters of the ground motions in active extensional ar- eas in Italy, Bull. Seism. Soc. Am., 78, 530-560. Salazar, W., V. Sardina and J. de Cortina (2007). A hy- brid inversion technique for the evaluation of source, path, and site effects employing S-wave spectra for subduction and upper-crustal earth- quakes in El Salvador, Bull. Seism. Soc. Am., 97, 208-221. Sepahvand, M. R., F. Yaminifard, M. Tatar and M. R. Abbassi (2012). Aftershocks study of the 2006 Silakhur earthquake (Zagros, Iran): seismological evidences for a pull-apart basin along the Main Re- cent Fault, Doroud segments, J. Seismol., 16, 233- 251. Shahvar, M. P., M. Zare and S. Castellaro (2013). A uni- fied seismic catalog for the Iranian plateau (1900 2011), Seismol. Res. Lett., 84, 233-249. Sharma, B., k. A. Gupta, k. D. Devi, D. kumar, S. S. Teotia and B. k. Rastogi (2008). Attenuation of high frequency seismic waves in kachchh region, Gujarat, India, Bull. Seism. Soc. Am., 98, 2325- 2340. Shoja-Taheri, J., S. Nasserrieh, and H. Ghofrani (2007). ML and MW scales in the Iranian Plateau based on the strong-motion records, Bull. Seism. Soc. Am., 97,661-669. Singh, S. k. and M. Ordaz (1994). Seismic energy re- lease in Mexican subduction zone earthquakes, Bull. Seism. Soc. Am., 84, 1533-1550. Tchalenko, J.S. and J. Braud (1974). Seismicity and structure of the Zagros (Iran): The Main Recent Fault between 33° and 35°N. Phil. Trans. R. Soc. London, 277, 1 -25. Yoshimoto, k., H. Sato and M. Ohtake (1993). Fre- quency dependent attenuation of P and S waves in the kanto area, Japan, based on the coda normal- ization method, Geophys. J. Int., 114, 165-174. Zafarani H, and B. Hassani (2013). Site response and source spectra of S-waves in the Zagros region, Iran, J. Seismol., 17, 645-666. Zare, M., H. Amini, P. Yazdi, k. Sesetyan, M.B. Demir- cioglu, D. kalafat, M. Erdik, D. Giardini, M. Asif khan and N. Tsereteli (2014). Recent developments of the Middle East catalog, J. Seismol., 18, 749-772. *Corresponding author: Gholam Javan-Doloei, International Institute of Earthquake Engineering and Seismol- ogy, Tehran, Iran; email: Javandoloei@iiees.ac.ir. © 2017 by the Istituto Nazionale di Geofisica e Vulcanologia. All rights reserved. ATTENUATION, SOURCE AND SITE EFFECTS FROM INVERSION OF THE SILAkHOR AFTERSHOCkS