The Effects of Attenuation and Site on the Spectra of Microearthquakes in the Shillong Region of Northeast India

Microearthquake spectra from the Shillong region are analyzed to observe the effect of attenuation and site on these spectra. The spectral ratio method is utilized to estimate the Q values for both P- and S-waves in the subsurface layer, wherein the ratio of spectral amplitudes at lower and higher frequencies are taken into consideration for three stations at varying epicentral distances. Average estimates of Q _P and Q _S are 178 and 195. The ratio of Q _S to Q _P is estimated to be greater than 1 in major parts of the Shillong area, which can be related to the dry crust prevailing in the Shillong region. Typically, the variation in corner frequencies for these spectra is inferred to be characteristic of the site. Simultaneously, observations from spectral content of local earthquakes recorded at two different stations with respect to the reference site yield greater amplification of incoming seismic signals in the frequency range of 2–5 Hz, which is found to be well supported by the existing local lithology pertinent to that region.     

The seismic wave absorption in the crust and upper mantle in the vicinity of the Kislovodsk seismic station

The Q-factor estimates of the Earth’s crust and upper mantle as the functions of frequency (Q(f)) are obtained for the seismic S-waves at frequencies up to ~35 Hz. The estimates are based on the data for ~40 earthquakes recorded by the Kislovodsk seismic station since 2000. The magnitudes of these events are M_W > 3.8, the sources are located in the depth interval from 1 to 165 km, and the epicentral distances range from ~100 to 300 km. The Q-factor estimates are obtained by the methods developed by Aki and Rautian et al., which employ the suppression of the effects of the source radiation spectrum and local site responses in the S-wave spectra by the coda waves measured at a fixed lapse time (time from the first arrival). The radiation pattern effects are cancelled by averaging over many events whose sources are distributed in a wide azimuthal sector centered at the receiving site. The geometrical spreading was specified in the form of a piecewise-continuous function of distance which behaves as 1/R at the distances from 1 to 50 km from the source, has a plateau at 1/50 in the interval from 50–70 km to 130–150 km, and decays as \({\raise0.7ex\hbox{$1$} \!\mathord{\left/ {\vphantom {1 {\sqrt R }}}\right.\kern-\nulldelimiterspace} \!\lower0.7ex\hbox{${\sqrt R }$}}\) beyond 130–150 km. For this geometrical spreading model and some of its modifications, the following Q-factor estimates are obtained: Q(f) ~ 85f^0.9 at the frequencies ranging from ~1 to 20 Hz and Q(f) ~ 75f^1.0 at the frequencies ranging from ~1 to 35 Hz.     

f max and fault zone property of Lushan earthquake of 20 April 2013, Sichuan,China

In this study, we determined f _max from near-field accelerograms of the Lushan earthquake of April 20, 2013 through spectra analysis. The result shows that the values of f _max derived from five different seismography stations are very close though these stations roughly span about 100 km along the strike. This implies that the cause of f _max is mainly the seismic source process rather than the site effect. Moreover, according to the source–cause model of Papageorgiou and Aki (Bull Seism Soc Am 73:693–722, 1983), we infer that the cohesive zone width of the rupture of the Lushan earthquake is about 204 with an uncertainty of 13 m. We also find that there is a significant bulge between 30 and 45 Hz in the amplitude spectra of accelerograms of stations 51YAL and 51QLY, and we confirm that it is due to seismic waves’ reverberation of the sedimentary soil layer beneath these stations.     

Aftershock Activity and Frequency-dependent Low Coda Qc in the Epicentral Region of the 1999 Chamoli Earthquake of Mw 6.4

— On 28 March, 1999 (19:05:10.09, UT) a significant earthquake of M _w 6.4 occurred in the Garhwal Himalaya (30.555°N, 79.424°E). One hundred and ten well-recorded aftershocks show a WNW-ESE trending northeasterly dipping seismic zone extending from a depth of 2 to 20?km. As the main shock hypocenter occurred at the northern end of this seismic zone and aftershocks extended updip, it is inferred that the main-shock rupture nucleated on the detachment plane at a depth of 15?km and then propagated updip along a NE-dipping thrust plane. Further, the epicentral distribution of aftershocks defines a marked concentration near a zone where main central thrust (MCT) takes a significant turn towards the north, which might be acting as an asperity in response to the NNE compression due to the underthrusting of Himalayan orogenic process prevalent in the entire region. Presence of high seismicity including five earthquakes of magnitude exceeding 6 and twelve earthquakes of magnitude exceeding 5 in the 20th century is presumed to have caused a higher level of shallow crustal heterogeneity in the Garhwal Himalaya, a site lying in the central gap zone of the Himalayan frontal arc. Attenuation property of the medium around the epicentral area of the 1999 Chamoli earthquake, covering a circular area of 61,500?km² with a radius of 140?km, is studied by estimating the coda Q _c from 48 local earthquakes of magnitudes varying from 2.5–4.8. These earthquakes were recorded at nine 24-bit REFTEK digital stations; two of which were equipped with three-component CMG40T broadband seismometers and others with three-component L4-3D short-period seismometers. The estimated Q _o values at different stations suggest on average a low value of the order of (30?±?0.8), indicating an attenuating crust beneath the entire region. The frequency-dependent relation indicates a relatively low Q _c at lower frequencies (1–3?Hz) that can be attributed to the loss of energy due to scattering on heterogeneities and/or the presence of faults and cracks. The large Q _c at higher frequencies may be related to the propagation of backscattered body waves through deeper parts of the lithosphere where less heterogeneities are expected. An important observation is that the region north of MCT (more rigid highly metamorphosed crystalline rocks) is less attenuative in comparison to the region south of MCT (less rigid slightly metamorphosed rocks (sedimentary wedge)). The acceleration decays to 50% at 20?km distance and to 7% at 100?km. Hence, even 1g acceleration at the source may not cause significant damage beyond 100?km in this region.     

Characteristics of seismic wave attenuation in the crust and upper mantle of the Northern Caucasus

The quality factor Q as a function of frequency in an S wave range of 1–8 Hz is estimated from records of ～60 earthquakes (M _w > 3.9 and source depths of 1–60 km) obtained at the Sochi seismic station at epicentral distances of less than ～300 km. Methods of Q estimation used in the paper were developed in works by Aki, Rautian, and others; they are based on the suppression of source-related and local effects in S wave spectra with the help of coda waves measured at a fixed time from the first arrival. To compensate for directivity effects, averaging was performed over the set of events whose sources were located in a wide range of back azimuths. The geometric divergence is represented as a three-segment function: 1/R, 1, and 1/√R at epicentral distances of 1/50–1/70 to 50–70 km, 50–70 to 130–150 km, and greater than 130–150 km, respectively. The geometric divergences in this model yielded the following estimates of the quality factor: Q(f) ～ 80f ^0.9 with a base of 35–280 km and Q(f) ～ 110f ^0.8 with a base of 60–280 km. The resulting combinations of the propagation path effects (Q and the geometric divergence) can be used for predicting strong motion parameters in the Northern Caucasus.     

Stochastic finite-fault simulation of ground motion from the August 11, 2012, <Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript> 6.4 Ahar earthquake,northwestern Iran

In this study, the 11 August 2012 M _w 6.4 Ahar earthquake is investigated using the ground motion simulation based on the stochastic finite-fault model. The earthquake occurred in northwestern Iran and causing extensive damage in the city of Ahar and surrounding areas. A network consisting of 58 acceleration stations recorded the earthquake within 8–217 km of the epicenter. Strong ground motion records from six significant well-recorded stations close to the epicenter have been simulated. These stations are installed in areas which experienced significant structural damage and humanity loss during the earthquake. The simulation is carried out using the dynamic corner frequency model of rupture propagation by extended fault simulation program (EXSIM). For this purpose, the propagation features of shear-wave including \( {Q}_s \) value, kappa value \( {k}_0 \), and soil amplification coefficients at each site are required. The kappa values are obtained from the slope of smoothed amplitude of Fourier spectra of acceleration at higher frequencies. The determined kappa values for vertical and horizontal components are 0.02 and 0.05 s, respectively. Furthermore, an anelastic attenuation parameter is derived from energy decay of a seismic wave by using continuous wavelet transform (CWT) for each station. The average frequency-dependent relation estimated for the region is \( Q=\left(122\pm 38\right){f}^{\left(1.40\pm 0.16\right)}. \) Moreover, the horizontal to vertical spectral ratio \( H/V \) is applied to estimate the site effects at stations. Spectral analysis of the data indicates that the best match between the observed and simulated spectra occurs for an average stress drop of 70 bars. Finally, the simulated and observed results are compared with pseudo acceleration spectra and peak ground motions. The comparison of time series spectra shows good agreement between the observed and the simulated waveforms at frequencies of engineering interest.     

Influence of attenuation and site on microearthquakes’ spectra in Shillong region, of Northeast, India: A case study

We examine the influence of attenuation and site on the spectra of microearthquakes having origin within the Shillong region. The ratios of spectral amplitudes at lower and higher frequencies are measured for three different stations at varying epicentral distances to estimate Q value for both P- and S-wave in near and sub-surface layer. The average estimates of Q _P and Q _S are found to be 178 and 195. The ratio of Q _S to Q _P emerges to be greater than unity in major parts of the Shillong area, suggesting dominance of dry crust prevailing in Shillong region. The variation in corner frequencies for these spectra is inferred to be characteristics of the site. Besides, the disparity in spectral content with reference to hard rock site yields the inference that the incoming seismic signals get amplified considerably while traversing from southern part to northeastern part of Shillong, best outlined at 2 to 5 Hz, which is well corroborated by the existing lithology.     

Separation of Source, Propagation and Site Effects from S Waves¶of Local Earthquakes in Bursa Region, Northwestern Turkey

—?We have used micro-earthquake recordings (M= 1.8–4.1) of local events in the distance range of 5–60?km in order to quantify the attenuation and site effects in the vicinity of the Bursa city, Marmara region, Turkey. The data set consists of 120 three-component recorded accelograms from 69 earthquakes, recorded at six stations. Each station is deployed on different geologic units, such as massive limestone, slope deposit and Quaternary young sediments, in the framework of the Marmara Poly-Project.¶In this study a nonparametric inversion method was applied to acceleration records from the Bursa region to estimate source, site and path effects using a two-step inversion. At the first step, we determined attenuation functions by analyzing the distance dependence of the spectral amplitudes and retrieved values of Q _s (f) = 46.59f ^0.67. At the second step, the corrected S-waves spectral records for the attenuation function, including the geometrical spreading effect, were inverted to separate source and site response for 21 different frequencies selected between 0.5 and ～25?Hz. The near-surface attenuation, κ value, was also estimated by using the model proposed by Anderson and Hough (1984) at each site. We observed that κ₀ is smaller for stations located on rock site (I?dιr, SIGD, κ₀～0.004) compared to the one that is located on Neogene sediment (Çukurca, SCKR, κ₀～0.018).¶Site amplifications from inversion showed that the station located within the Bursa basin, Çukurca (SCKR), is the most important site with about 4.0 amplification value at 1.8?Hz. Demirta? (SDEM) amplifies the spectral amplitudes about 3.0 times at 2.0?Hz, SHMK about 3.5 times between 2.5 and 3.5?Hz and SHMT nearly reaching 3.5 times between 1.5 and 4.0?Hz. However, stations located on the Uluda? Mountain Massif (SKAY and SIGD), which correspond to a deep limestone geological unit, have the smallest amplification, that values between 0.6 and 1.4.     

Coda Q Estimates in the Koyna Region, India

—The coda Q, Q _c ?, have been estimated for the Koyna region of India. The coda waves of 76 seismograms from thirteen local earthquakes, recorded digitally in the region during July–August, 1996, have been analyzed for this purpose at nine central frequencies viz., 1.5, 2.0, 3.0, 4.0, 6.0, 8.0, 12.0, 16.0 and 24.0 Hz using a single backscattering model. All events with magnitude less than 3 fall in the epicentral distances less than 60 km and have focal depths which range from 0.86 to 9.43 km. For the 30 sec coda window length the estimated Q _c values vary from 81 to 261 at 1.5 Hz and 2088 to 3234 at 24 Hz, whereas the mean values of Q _c with the standard error vary from 148 ± 13.5 at 1.5 Hz to 2703 ± 38.8 at 24 Hz. Both the estimated Q _c values and their mean values exhibit the clear dependence on frequency in the region and a frequency dependence average attenuation relationship, Q _c = 96f ^1.09, has been obtained for the region, covering an approximate area of 11500 km² with the surfacial extent of about 120 km and depth of 60 km.¶Lapse time dependence of Q _c has also been studied for the region, with the coda waves analyzed at five lapse time windows from 20 to 60 sec duration with the difference of 10 sec. The frequency dependence average Q _c relationships obtained at these window lengths Q _c = 66f ^1.16 (20 sec), Q _c = 96f ^1.09 (30 sec), Q _c =131f ^1.04 (40 sec), Q _c = 148f ^1.04 (50 sec), Q _c = 182f ^1.02 (60 sec) show that the frequency dependence (exponentn) remains mostly stationary at all the lapse time window lengths, while the change in Q ₀ value is significant. Lapse time dependence of Q _c in the region is also interpreted as the function of depth.     

大容量气枪震源长江定点激发信号检测

地学长江计划“安徽段实验”是大容量气枪震源在长江的首次激发。本文针对布设在气枪固定激发点附近的流动台和周边固定台接收到的气枪信号进行线性叠加分析近场和远场信号的时频特性,利用叠加结果检测气枪信号的传播特性,分析不同环境因素对信号传播距离的影响。结果表明：①近岸首台可以接收到清晰的压力脉冲、气泡脉冲的体波和面波信号;②气枪信号主频为5Hz左右,随震中距的增加,压力脉冲信号衰减很快,信号主频频带变窄;③对信号传播距离进行初步检测,最近的传播距离为180km,最远共有3个激发点传播达到260km,夜晚激发信号传播距离较远。     

Influence of wind turbines on seismic stations in the upper rhine graben,SW Germany

By analysing long- and short-term seismological measurements at wind farms close to the town of Landau, SW Germany, we present new insights into ground motion signals from wind turbines (WTs) at local seismic stations. Because of their need to be located in similar regions with sparsely anthropogenic activities, wind turbines impact seismic stations and their recordings in a way that is not yet fully understood by researchers. To ensure the undisturbed recording tasks of a regional seismic array or a single station by a protected area around those endangered stations, it is very important to investigate the behavior of WTs as a seismic source. For that reason, we calculate averaged one-hour long spectra of the power spectral density (PSD) before and after the installation of a new wind farm within the investigated area. These PSD are ordered according to the rotation speed. We observe a clear increase of the PSD level after the WT installation in a frequency range of 0.5 to 10 Hz up to a distance of 5.5 km away from the WT. By analysing seismic borehole data, we also observe a decrease of the PSD of wind dependent signals with depth. The impact of wind-dependent signals is found to be much more pronounced for the shallower station (150 m depth) than for the deeper one (305 m depth). Using short-term profile measurements, we fit a power-law decay proportional to 1/r ^b to the main WT-induced PSD peaks and differentiate between near-field and far-field effects of ground motions. For low frequencies in the range from 1 to 4 Hz, we determine a b value of 0.78 to 0.85 for the far field, which is consistent with surface waves. The b value increases (up to 1.59) with increasing frequencies (up to 5.5 Hz), which is obviously due to attenuating effects like scattering or anelasticity. These results give a better understanding of the seismic wavefield interactions between wind turbines (or wind farms) with nearby seismic stations, including borehole installations, in a sedimentary setting.     

High-frequency P- and S-wave Attenuation in the Earth

—Investigations of the spectral characteristics of teleseismic body waves revealed that the spectral falloff rate between 1 Hz and 10 Hz is primarily controlled by anelastic attenuation along the path. In addition, the amount of high-frequency energy in teleseismic body waves is far above the level expected on the basis of Q estimates at low frequencies, thus leading to the idea of frequency dependence in Q. Q variations in the earth’s mantle can be investigated by mapping out the variations of high frequency (4 - 10 Hz) energy relative to the low frequency (1 - 3 Hz) energy in teleseismic P waves, and similar ratios at lower frequencies in teleseismic S waves. Because of the extreme sensitivity of spectral content of short-period body waves to Q variations, large uncertainties in other factors affecting spectral content can be tolerated in such studies. With the increasing number and density of broadband seismic stations recording at high sampling rates, tomographic studies of Q at high frequencies become possible.     

Research on the regionalQ values of coda waves for the aftershocks in Luquan,Yunnan of 1985

Based on the scattering coda model by which local and regional earthquakes are interpreted (K. Aki, 1969), and using observational coda data of 68 aftershocks of the 1985 Luquan, Yunnan earthquake registered by the VGK seismographs installed at 12 stations in the Yunnan regional short-period network, theQ-values of coda waves are calculated respectively for 6 time intervals. It is observed that within the frequency range of 0.40–1.65 Hz of the observed data, theQ-values are closely related with the frequencies and the calculated codaQ ranges between 80–240 with the coefficient of frequency dependence η=0.45. The calculated source factorsB(f> _p) of the coda waves which indicate the scattering strength are mostly within the order 10^?23–10^?24. Areas with lowQ-values present high scattering. It should be noted that by comparing data obtained before and after the Luquan earthquake, clear changes can be detected in theQ-values measured at stations close to the epicentral region, and that theQ-values of the aftershock coda are less than about one half of the pre-shock values. It may be mentioned that the time-dependent regional variations of theQ-values might possibly bring about practical significance in earthquake prediction. Moreover, aftershock focal parameters are determined. Through discussions on the quantitative relations between the focal parameters, we get: 1gE=1.59M _L+ 11.335;E=(2.10 × 10^?5)M ₀; length of focal rupturea=0.40?0.80 km for 3.0≤M _L<5.0 events; stress drop Δσ=(6.0–130) ×10⁵ Pa. Through interpretation of the data, we have also learned the important characteristics that there is no linear relation between the stress drops and the earthquake magnitudes.     

Frequency relationship for seismic Qβ of central southern Italy from accelerograms for the Irpinia earthquake (1980)

The frequency dependence of Q_β for seismic waves in a distance range with a maximum of 150 km from the epicentre of the Irpinia earthquake of November 23, 1980 has been sought using displacement spectral ratios computed from strong-motion accelerograms recorded in the region. The method has been applied to calculate the behaviour of Q_β as a function of frequency in the band 0.1–25 Hz, and to investigate whether azimuthal variations appear in seismic Q_β for the lithosphere in central southern Italy. The same result is obtained using data from stations in western south Italy as using data from eastern south Italy, namely,

$\text{Q}{}_{\mathrm{\beta }}\text{(f) = 40f (}\text{Hz}\text{)}$

The linear relationship suggest that apparent Q_β depends more on the scale of heterogeneity of the lithosphere, affecting reflection and scattering mechanisms, than on intrinsic energy losses related to the anelasticity of the materials through which the seismic waves propagate.The existence of a peak in Q_β^?1 has been investigated in the low-frequency band (0.1–2.5 Hz) using a higher resolution power. A stable result in this low-Q_β zone is not possible on the basis of the available data: only in six Q_β(f) profiles does an evident minimum exist, between 0.2 and 1 Hz, while in nine cases the curves are monotonically increasing from the lowest observable frequencies; a further nine cases appear of uncertain interpretation.     

Near-source attenuation of high-frequency body waves beneath the New Madrid Seismic Zone

Attenuation characteristics in the New Madrid Seismic Zone (NMSZ) are estimated from 157 local seismograph recordings out of 46 earthquakes of 2.6?≤?M?≤?4.1 with hypocentral distances up to 60 km and focal depths down to 25 km. Digital waveform seismograms were obtained from local earthquakes in the NMSZ recorded by the Center for Earthquake Research and Information (CERI) at the University of Memphis. Using the coda normalization method, we tried to determine Q values and geometrical spreading exponents at 13 center frequencies. The scatter of the data and trade-off between the geometrical spreading and the quality factor did not allow us to simultaneously derive both these parameters from inversion. Assuming 1/R ^1.0 as the geometrical spreading function in the NMSZ, the Q _P and Q _S estimates increase with increasing frequency from 354 and 426 at 4 Hz to 729 and 1091 at 24 Hz, respectively. Fitting a power law equation to the Q estimates, we found the attenuation models for the P waves and S waves in the frequency range of 4 to 24 Hz as Q _P?=?(115.80?±?1.36) f ^{(0.495?±?0.129)} and Q _S?=?(161.34?±?1.73) f ^{(0.613?±?0.067)}, respectively. We did not consider Q estimates from the coda normalization method for frequencies less than 4 Hz in the regression analysis since the decay of coda amplitude was not observed at most bandpass filtered seismograms for these frequencies. Q _S/Q _P?>?1, for 4?≤?f?≤?24 Hz as well as strong intrinsic attenuation, suggest that the crust beneath the NMSZ is partially fluid-saturated. Further, high scattering attenuation indicates the presence of a high level of small-scale heterogeneities inside the crust in this region.     

Attenuation of P- and S-waves in the Chamoli Region, Himalaya, India

The attenuation properties of the crust in the Chamoli region of Himalaya have been examined by estimating the frequency-dependent relationships of quality factors for P waves (Q_α) and for S waves (Q_β) in the frequency range 1.5–24 Hz. The extended coda normalization method has been applied on the waveforms of 25 aftershocks of the 1999 Chamoli earthquake (M 6.4) recorded at five stations. The average value of Q_α is found to be varied from 68 at 1.5 Hz to 588 at 24 Hz while it varies from 126 at 1.5 Hz to 868 at 24 Hz for Q_β. The estimated frequency-dependent relations for quality factors are Q_α = (44 ± 1)f^(0.82±.04) and Q_β = (87 ± 3)f^(0.71±.03). The rate of increase of Q(f) for P and S waves in the Chamoli region is comparable with the other regions of the world. The ratio Q_β/Q_α is greater than one in the region which along with the frequency dependence of quality factors indicates that scattering is an important factor contributing to the attenuation of body waves in the region. A comparison of attenuation relation for S wave estimated here (Q_β = 87f^0.71) with that of coda waves (Q_c = 30f^1.21) obtained by Mandal et al. (2001) for the same region shows that Q_c > Q_β for higher frequencies (>8 Hz) in the region. This indicates a possible high frequency coda enrichment which suggests that the scattering attenuation significantly influences the attenuation of S waves at frequencies >8 Hz. This observation may be further investigated using multiple scattering models. The attenuation relations for quality factors obtained here may be used for the estimation of source parameters and near-source simulation of earthquake ground motion of the earthquakes, which in turn are required for the assessment of seismic hazard in the region.     

Determination and application of path duration of seismic ground motions based on the K-NET data in Sagami Bay,Japan

Duration models are one of the important parameters in ground-motion simulations. This model varies in different study areas, and plays a critical role in nonlinear structural response analysis. Currently, available empirical models are being globally used in ground-motion simulations, with limited research focusing on path duration in specific regions. In this study, we collected 6,486 sets of three-component strong-motion records from 29 K-NET stations in the Sagami Bay, Japan, and its surrounding areas between January 2000 to October 2018. We extracted the effective duration of 386 pieces of ground-motion records by manually picking up the S-wave arrival time and calculating the significant duration. We then obtained the path duration model of the study area based on the empirical equation of dynamic corner frequency and source duration of [7]. Compared with the results of the available empirical models, the Fourier spectrum of the simulated ground motion from our effective duration model showed higher accuracy in the long-term range, with less fitting residuals. This path duration model was then applied to simulate two earthquakes of M_W5.4 and M_W6.2, respectively, in the region using the stochastic finite-fault method with a set of reliable source, path, and site parameters determined for the study area. The simulation results of most stations fit well with observation records in the 0–30 Hz frequency band. For the M_W5.4 earthquake, the simulated ground motions at KNG005/KNG010/SZO008 stations were relatively weak in the mid to high frequency band (1–30 Hz) because the quality factor and geometric diffusion model used in the simulation were the averages of the entire Sagami Bay region, causing a bias in the results of a few stations owing to local crustal velocity anomalies and topographic effects. For the M_W6.2 earthquake, the simulated ground motions were relatively weak at all SZO and TKY stations, mainly because of the close distance from these stations to the epicenter and the complex seismic-wave propagation paths. The analysis suggests that the differences between the simulation results of the two earthquakes were mainly related to complex geological conditions and seismic-wave propagation paths.     

Rupture process and aftershock focal mechanisms of the 2022 M6.8 Luding earthquake in Sichuan

According to the China Earthquake Networks Center, a strong earthquake of M6.8 occurred in Luding County, Ganzi Tibetan Autonomous Prefecture, Sichuan Province, China (102.08°E, 29.59°N), on September 5, 2022, with a focal depth of 16 km. Rapid determination of the source parameters of the earthquake sequence is vital for post-earthquake rescue, disaster assessment, and scientific research. Near-field seismic observations play a key role in the fast and reliable determination of earthquake source parameters. The numerous broadband seismic stations and strong-motion stations recently deployed by the National Earthquake Intensity Rapid Report and Early Warning project have provided valuable real-time near-field observation data. Using these near-field observations and conventional mid- and far-field seismic waveform records, we obtained the focal mechanism solutions of the mainshock and M ≥ 3.0 aftershocks through the waveform fitting method. We were further able to rapidly invert the rupture process of the mainshock. Based on the evaluation of the focal mechanism solution of the mainshock and the regional tectonic setting, we speculate that the Xianshuihe fault formed the seismogenic structure of the M6.8 strong earthquake. The aftershocks formed three spatially separated clusters with distinctly different focal mechanisms, reflecting the segmented nature of the Xianshuihe fault. As more high-frequency information has been applied in this study, the absolute location of the fault rupture is better constrained by the near-field strong-motion data. The rupture process of the mainshock correlates well with the spatial distribution of aftershocks, i.e., aftershock activities were relatively weak in the maximum slip area, and strong aftershock activities were distributed in the peripheral regions.     

龙门山地区的kappa（κ0）模型及汶川MS8.0地震的强地震动模拟

在地震动模拟中,高频衰减参数（κ₀）是一个重要的参数,它控制了傅里叶谱高频部分的衰减特性.本文利用汶川M_S8.0地震和芦山M_S7.0地震主余震的1671组强震动观测记录,计算了龙门山地区50个断层距小于150 km的强震动台站的κ₀,基于随机有限断层法模拟了汶川地震中这些台站的加速度时程、傅里叶振幅谱和反应谱,并与前人的研究进行比较.结果表明,合理计算的κ₀可以有效地改善加速度时程振幅和高频谱（>1 Hz）的模拟结果.另一方面,本文基于κ₀与地形高程的相关性,建立了龙门山地区的κ₀经验模型.分别采用该经验模型和κ₀= 0.04 s模拟了汶川地震的峰值加速度分布,并与观测记录进行比较.结果表明,采用本文提出的κ₀模型可以更好地重现汶川地震的峰值加速度分布,特别是在断层破裂方向的反方向区域和山区.综上,汶川地震中山区的峰值加速度明显大于盆地地区的现象,不仅与断层破裂产状有关,还与山区和盆地地区的κ₀之间显著的差异有关.     

Features of Radiation and Propagation of Seismic Waves in the Northern Caucasus: Manifestations in Regional Coda

Based on the Anapa (ANN) seismic station records of ~40 earthquakes (M_W > 3.9) that occurred within ~300 km of the station since 2002 up to the present time, the source parameters and quality factor of the Earth’s crust (Q(f)) and upper mantle are estimated for the S-waves in the 1–8 Hz frequency band. The regional coda analysis techniques which allow separating the effects associated with seismic source (source effects) and with the propagation path of seismic waves (path effects) are employed. The Q-factor estimates are obtained in the form Q(f) = 90 × f 0.7 for the epicentral distances r < 120 km and in the form Q(f) = 90 × f1.0 for r > 120 km. The established Q(f) and source parameters are close to the estimates for Central Japan, which is probably due to the similar tectonic structure of the regions. The shapes of the source parameters are found to be independent of the magnitude of the earthquakes in the magnitude range 3.9–5.6; however, the radiation of the high-frequency components (f > 4–5 Hz) is enhanced with the depth of the source (down to h ~ 60 km). The estimates Q(f) of the quality factor determined from the records by the Sochi, Anapa, and Kislovodsk seismic stations allowed a more accurate determination of the seismic moments and magnitudes of the Caucasian earthquakes. The studies will be continued for obtaining the Q(f) estimates, geometrical spreading functions, and frequency-dependent amplification of seismic waves in the Earth’s crust in the other regions of the Northern Caucasus.