Stochastic ground-motion simulation of two Himalayan earthquakes: seismic hazard assessment perspective

The earthquakes in Uttarkashi (October 20, 1991, M _w 6.8) and Chamoli (March 8, 1999, M _w 6.4) are among the recent well-documented earthquakes that occurred in the Garhwal region of India and that caused extensive damage as well as loss of life. Using strong-motion data of these two earthquakes, we estimate their source, path, and site parameters. The quality factor (Q _β ) as a function of frequency is derived as Q _β (f) = 140f ^1.018. The site amplification functions are evaluated using the horizontal-to-vertical spectral ratio technique. The ground motions of the Uttarkashi and Chamoli earthquakes are simulated using the stochastic method of Boore (Bull Seismol Soc Am 73:1865–1894, 1983). The estimated source, path, and site parameters are used as input for the simulation. The simulated time histories are generated for a few stations and compared with the observed data. The simulated response spectra at 5% damping are in fair agreement with the observed response spectra for most of the stations over a wide range of frequencies. Residual trends closely match the observed and simulated response spectra. The synthetic data are in rough agreement with the ground-motion attenuation equation available for the Himalayas (Sharma, Bull Seismol Soc Am 98:1063–1069, 1998).     

Ground motion observations and simulation for local earthquakes in the Campi Flegrei volcanic area

Ground motion produced by low magnitude earthquakes can be used to predict peak values in high seismic risk areas where large earthquakes data are not available. In the present work 20 local earthquakes (M_D∈[?0.3, 2.2]) occurred in the Campi Flegrei caldera during the last decade were analyzed. We followed this strategy: empirical relations were used to calibrate synthetic modeling, accounting for the source features and wave propagation effects. Once the source and path parameters of ground motion simulation were obtained from the reference data set, we extrapolated scenarios for stronger earthquakes for which real data are not available. The procedure is structured in two steps: (1) evaluation of ground motion prediction equation for Campi Flegrei area and assessment of input parameters for the source, path and site effects in order to use the finite fault stochastic approach (EXSIM code); (2) simulation of two moderate-to-large earthquake scenarios for which only historical data or partial information are available (M_w4.2 and M_w5.4). The results show that the investigated area is characterized by high attenuation of peak amplitude and not negligible site effects. The stochastic approach has revealed a good tool to calibrate source, path and site parameters on small earthquakes and to generate large earthquake scenario. The investigated magnitude range represents a lower limit to apply the stochastic method as a calibration tool, due to the small size of involved faults (fault length around 200/300 m).     

Estimation of site amplification,attenuation and source spectra of S-waves in the East-Central Iran

Generalized inversion of the S-wave amplitude spectra from the strong-motion network data in the East-Central Iran has been used to estimate simultaneously source parameters, site response and the S-wave attenuation (Q_s). In this regard, 190 three-component records were used corresponded to 40 earthquakes with the magnitudes M3.5–M7.3. These earthquakes were recorded at 42 stations in the hypocentral distance range from 9 to 200 km. The inverse problem was solved in 20 logarithmically equally spaced points in the frequency band from 0.4 to 15 Hz. The frequency-dependent site amplification was imposed, as a constraint, on two reference site responses in order to remove the undetermined degree of freedom in the inversion and obtain a unique inverse solution. Also, a geometrical spreading factor was assumed for removing the trade-off between geometrical spreading and anelastic attenuation. Different source parameters, such as seismic moment (M₀), seismic energy (E_s), corner frequency (f_c) and Brune stress drop (Δσ), were estimated for each event by fitting an ω² model to the spectra obtained from the inversion. The stress drop values of earthquakes, obtained in this research, are in good agreement with those of other studies. Also average site response values were correlated to the average shear wave velocities in the uppermost 30 m, in high and low frequency bands. The peak frequencies of site amplifications, estimated by the generalized inversion method, where in good agreement with those of horizontal to vertical (H/V) spectral ratios for the S-wave portion of records. However, no perfect matching in amplitude was obtained due to the deficiencies of the H/V ratio technique. By supposing a free shape for Q factor, a frequency dependent function was found, the logarithm of which could be approximated by a linear function, Q(f)=151f^0.75. The uncertainties of model parameters have been evaluated by covariance matrix of least-square fit. The residuals were also analyzed in order to assess the validity of the model. The analysis of residuals with respect to magnitude and distance indicates that they are distributed normally with approximately zero mean. The robustness of the results has been studied concerning their sensitivities to the omission of different datasets, selected randomly from original database. The results obtained here can be used in predicting ground-motion parameters applying stochastic methods.     

A study of Guptkashi,Uttarakhand earthquake of 6 February 2017 (<Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript>5.3) in the Himalayan arc and implications for ground motion estimation

The 2017 Guptkashi earthquake occurred in a segment of the Himalayan arc with high potential for a strong earthquake in the near future. In this context, a careful analysis of the earthquake is important as it may shed light on source and ground motion characteristics during future earthquakes. Using the earthquake recording on a single broadband strong-motion seismograph installed at the epicenter, we estimate the earthquake’s location (30.546° N, 79.063° E), depth (H?=?19 km), the seismic moment (M₀?=?1.12×10¹⁷ Nm, M _w 5.3), the focal mechanism (φ?=?280°, δ?=?14°, λ?=?84°), the source radius (a?=?1.3 km), and the static stress drop (Δσ _s ~22 MPa). The event occurred just above the Main Himalayan Thrust. S-wave spectra of the earthquake at hard sites in the arc are well approximated (assuming ω^?2 source model) by attenuation parameters Q(f)?=?500f^0.9, κ?=?0.04 s, and f_max?=?infinite, and a stress drop of Δσ?=?70 MPa. Observed and computed peak ground motions, using stochastic method along with parameters inferred from spectral analysis, agree well with each other. These attenuation parameters are also reasonable for the observed spectra and/or peak ground motion parameters in the arc at distances ≤?200 km during five other earthquakes in the region (4.6?≤?M _w ?≤?6.9). The estimated stress drop of the six events ranges from 20 to 120 MPa. Our analysis suggests that attenuation parameters given above may be used for ground motion estimation at hard sites in the Himalayan arc via the stochastic method.     

QLg tomography in Gujarat,Western India

We propose a novel Lg attenuation tomography model (Q_Lg tomography) for the state of Gujarat, Western India, using earthquake data recorded by the Gujarat Seismic Network, operated by the Institute of Seismological Research in Gandhinagar. The waveform dataset consist of 400 3-component recordings, produced by 60 earthquakes with magnitude (M_L) spanning from 3.6 to 5.1, recorded at 60 seismic stations having epicentral distances spanning between 200 and 500 km. Spectral amplitude decays for Lg wave displacement were obtained by generalized inversion at 17 frequencies spanning between 0.9 and 9 Hz. Lg wave propagation efficiency was measured by Lg/Pn spectral ratio categorizing as efficient ratio ≥6 for 86%, intermediate ratio of 3–6 for 10% and inefficient ratio <3 for 4% paths of total 400 ray paths. The earthquake size and quality of waveform recorded at dense network found sufficient to resolve lateral variation of Q_Lg in Gujarat.Average power-law attenuation relationship obtained for Gujarat as Q_Lg(f) = 234f^0.64, which corresponds to high attenuation in comparison to peninsular India shield region and other several regions around the world. Q_Lg tomography resolves the highly attenuating crust of extremely fractured Saurashtra region and tectonically active Kachchh region. The Gujarat average attenuation is also lying in between them. The low attenuation in Cambay and Narmada rift basins and extremely low attenuation in patch of Surendranagar area is identified. This study is the first attempt and can be utilized as pivotal criteria for scenario hazard assessment, as maximum hazard has been reported in highly attenuating tectonically active Kachchh region and in low attenuating Cambay, Narmada and Surendranagar regions. The site and source terms are also obtained along with the Q_Lg inversion. The estimated site responses are comparable with observed local geological condition and agree with the previously reported site amplifications at the same sites. The source terms are comparable with local magnitude estimated from Network. The M_w (Lg) is nearly equivalent to M_L (GSN) and the slight differences are noted for larger magnitude events.     

Earthquake source parameters and scaling relationships in Hungary (central Pannonian basin)

Fifty earthquakes that occurred in Hungary (central part of the Pannonian basin) with local magnitude $M_\textrm{L}$ ranging from 0.8 to 4.5 have been analyzed. The digital seismograms used in this study were recorded by six permanent broadband stations and 20 short-period ones at hypocentral distances between 10 and 327 km. The displacement spectra for P- and SH-waves were analyzed according to Brune’s source model. Observed spectra were corrected for path-dependent attenuation effects using an independent regional estimate of the quality factor Q _S . To correct spectra for near-surface attenuation, the κ parameter was calculated, obtaining it from waveforms recorded at short epicentral distances. The values of the κ parameter vary between 0.01 and 0.06 s with a mean of 0.03 s for P-waves and between 0.01 and 0.09 s with a mean of 0.04 s for SH-waves. After correction for attenuation effects, spectral parameters (corner frequency and low-frequency spectral level) were estimated by a grid search algorithm. The obtained seismic moments range from 4.21×10¹¹ to 3.41×10¹⁵ Nm (1.7?≤?M _w ?≤?4.3). The source radii are between 125 and 1,343 m. Stress drop values vary between 0.14 and 32.4 bars with a logarithmic mean of 2.59 bars (1 bar = 10⁵ Pa). From the results, a linear relationship between local and moment magnitudes has been established. The obtained scaling relations show slight evidence of self-similarity violation. However, due to the high scatter of our data, the existence of self-similarity cannot be excluded.     

Estimation of site effects using strong motion data of BYTNet array in Turkey

Simultaneous estimation of effects of source, propagation path, and local site amplification was carried out using observed strong motion records in a frequency range from 0.8 to 20 Hz for the purpose of empirical evaluation of the local site effects in different geological conditions in the northwestern part of Turkey. The analyzed data are S-wave portions of 162 accelerograms from 39 shallow events observed at 14 sites of BYTNet array. A spectral separation method was applied to the observed S-wave spectra. The solutions for source spectra, inelasticity factor of propagation path for S-waves (Q _s-value), and factor of site amplification at each site were obtained in a least squares sense. In the analysis, we assumed that the factor of the site amplification at a reference site is the same as that of theoretical amplification of S-waves to the soil model whose bottom layer has an S-wave velocity around 2.15 km/s. The estimated Q _s-value of the propagation path is modeled as Q _s(f)?=?87.4f^0.78. The estimated site amplifications are characterized into three groups. The sites in the first group belong to rock site with no dominant peaks at a frequency range of 2 to 10 Hz. The second group of hard soil sites is characterized with moderately dominant peaks at a frequency of 5 Hz. The last group for soft soil sites has common peaks at a frequency of 4 Hz with larger amplitudes than those in the hard soil group. We, then, compare the amplifications with average S-wave velocity in top 30 m of the shallow S-wave profiles and proposed linear empirical formula between them at each frequency. We, furthermore, inverted the observed amplification factors into S-wave velocity and Q _s-value profiles of the deep soil over the basement.     

The study of the radiation characteristics and propagation of seismic waves in the North Caucasus by modeling the accelerograms of the recorded earthquakes

For evaluating the parameters of the vibrations of the Earth’s surface in the case of strong earthquakes, which are possible in the future, the regular patterns of the emission and propagation of seismic waves in the North Caucasus regions are investigated. The regional parameters of emission and propagation of seismic waves are evaluated by solution of the inverse problems of stochastic modeling of the accelerograms of the earthquakes, recorded by the seismic station in Sochi. The horizontal components of the strongest earthquakes (M _w ～ 3.9?5.6), that occurred in 2002–2006 within a radius of ～300 km from the seismic station, with source depths up to 60 km are modeled. For calculations of accelerograms, estimates of the quality are used, obtained earlier for this region in the form: Q(f) ～ 80 ～ f ^0.9. The parameter settings are carried out, which determine the shapes of the source spectra, the amplification of the seismic waves in the Earth’s crust, the weakening of the waves at high frequencies (κ), the parameters that determine the shape and duration of accelerograms, etc. Sufficiently good agreement of the calculated and recorded accelerograms is obtained, the regional characteristics of emission and propagation of seismic waves, which can be used for prediction of the parameters of strong motions in the North Caucasus, are evaluated; however, in the future these characteristics should be studied in more detail.     

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.     

Investigation on regional attenuation of Vrancea (Romania) intermediate-depth earthquakes

This paper aims at investigating possible regional attenuation patterns in the case of Vrancea(Romania) intermediate-depth earthquakes.Almost 500 pairs of horizontal components recorded during 13 intermediate-depth Vrancea earthquakes are employed in order to evaluate the regional attenuation patterns.The recordings are grouped according to the azimuth with regard to the Vrancea seismic source and subsequently,Q models are computed for each azimuthal zone assuming similar geometrical spreading.Moreover,the local soil amplification which was disregarded in a previous analysis performed for Vrancea intermediate-depth earthquakes is now clearly evaluated.The results show minor differences between the four regions situated in front of the Carpathian Mountains and considerable differences in attenuation of seismic waves between the forearc and backarc regions(with regard to the Carpathian Mountains).Consequently,an average Q model of the type Q(f) = 115×f~(1.25) is obtained for the four forearc regions,while a separate Q model of the type Q(f) = 70×f~(0.90) is computed for the backarc region.These results highlight the need to evaluate the seismic hazard of Romania by using ground motion models which take into account the different attenuation between the forearc/backarc regions.     

The 1986 Dharamsala earthquake of Himachal Himalaya – estimates of source parameters,average intrinsic attenuation and site amplification functions

The 26th April 1986 Dharamsala earthquake (m_b 5.5) occurred in the Kangra region of Himachal Himalaya, which lies in the rupture zone of great Kangra earthquake of 1905. This was the first moderate sized earthquake to be recorded at a few sites of the strong ground motion array in the NW Himalaya. The accelerograms of this earthquake have been used to estimate its source parameters, site amplification functions and to estimate the effective shear wave attenuation factor Q_β in the frontal region of Himachal Himalaya. A double couple fault plane solution for the earthquake has been obtained based on the spectra of the transverse component of the accelerograms. The estimated values of the source parameters are seismic moment: 2.1×10²⁴ dyne-cm, static stress drop (Δσ): 36 bars, source radius (r): 2.8 km and moment magnitude (M_w): 5.4. The estimated average values of effective shear wave attenuation factor Q_β for various sites are in the range of 125 to 300 with an overall spatial average of 239. The influence of local site effects on the observed PGA values have been examined on the basis of site amplification functions.     

A single-station coda solution for source, attenuation and site factors

Based on the empirical formulation formed from coda observations, the digital waveforms from 33 local events with magnitude M _L ranging between 2.1 and 3.5, recorded at Changli station of Beijing Telemetered Seismographic Network from 1989 to 1991, are used to compute coda Q for the Changli region and the source factors of all earthquakes by single-station coda method. Furthermore, assuming a certain source model, we have obtained the station site frequency response and source spectra, as well as source parameters such as corner frequencies, seismic moments and stress drops and so on. Their variations with time are monitored before and after larger earthquakes. Because the coda method can more effectively reduce the influence of source radiation pattern and a particular propagation path than direct wave method, more data can be used and more accurate results can be obtained, which provided a possible approach to study the source properties and reveal the variation of source parameters before larger earthquakes. Contribution No. 97A0107, Institute of Geophysics, SSB, China.     

Scattering and anelastic attenuation of seismic energy in northern Greece

The relative contribution of scattering (Q _s ^–1 ) and intrinsic (Q _i ^–1 ) attenuation to the totalS-wave attenuation for the frequencies of 1.5, 3.0, 6.0 and 12.0 Hz has been studied by applying the radiative energy transfer theory, Data of local earthquakes which occurred in northern Greece and were recorded by the permanent telementered network of the Geophysical Laboratory of the University of Thessaloniki have been used. The results show that in this area the scattering attenuation is dominant over all frequencies while intrinsic attenuation is significantly lower. The estimatedQ _s ^–1 andQ _i ^–1 values have frequency dependences off ^–0.72 andf ^–0.45, respectively. The frequency dependence ofQ _s ^–1 is the same as that of the codaQ _c ^–1 , obtained by applying the single scattering model, which probably implies that the frequency dependence of the coda wave attenuation is attributed to the frequency dependence of the scattering attenuation.Q _c ^–1 values are very close to scattering attenuation for short lapse times, (10–20 sec), and intermediate between scattering and intrinsic attenuation for the longer lapse times, (50–100 sec). This difference is explained as the result of the depth-dependent attenuation properties and the multiple scattering effects.     

Q value of anelastic S-wave attenuation in Yunnan region

The study of seismic attenuation property is a major subject in seismology. Seismic waves recorded by seismic stations （seismographs） contain source effect, seismic wave propagation effect, site response of seismic stations and instrumental response. The path effect of seismic wave propagation, site response of seismic stations and instrumental response must be taken out in the study of source property with seismic data. The path effect of seismic wave propagation （seismic attenuation） involves an important influential factor, the anelastic attenuation of medium, which is measured with quality factor Q, apart from geometric attenuation with the distance. As a basic physical parameter of the Earth medium, Q value is essential for quantitative study of earthquakes and source property （e.g. determination of source parameters）, which is widely used in earthquake source physics and engineering seismology.     

Low CodaQcin the Epicentral Region of the 2001 Bhuj Earthquake ofMw7.7

On 26 January, 2001 (03:46:55,UT) a devastating intraplate earthquake of M_w 7.7 occurred in a region about 5 km NW of Bhachau, Gujarat (23.42°N, 70.23°E). The epicentral distribution of aftershocks defines a marked concentration along an E-W trending and southerly dipping (45°) zone covering an area of (60 × 40) km². The presence of high seismicity including two earthquakes of magnitudes exceeding 7.7 in the 200 years is presumed to have caused a higher level of shallow crustal heterogeneity in the Kutch area; a site lying in the seismic zone V (zone of the highest seismicity for potentially M8 earthquakes) on the seismic zoning map of India. Attenuation property of the medium around the epicentral area of the Bhuj earthquake covering a circular area of 61,500 km² with a radius of 140 km is studied by estimating the coda-Q_c from 200 local earthquakes of magnitudes varying from 3.0–4.6. The estimated Q₀ values at locations in the aftershock zone (high seismicity) are found to be low in comparison to areas at a distance from it. This can be attributed to the fact that seismic waves are highly scattered for paths through the seismically active and fractured zone but they are well behaved outside the aftershock zone. Distribution of Q₀ values suggests that the local variation in Q₀ values is probably controlled by local geology. The estimated Q₀ values at different stations suggest a low value of Q=(102 ± 0.80)^*f^{(0.98 ± 0.02)} indicating an attenuative 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 attenuation associated with heterogeneities and/or intrinsic attenuation due to fluid movement in the fault zone and fluid-filled 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 heterogeneity is expected. Based on the attenuation curve estimated for Q₀=102, the ground acceleration at 240 km distance is 13% of 1 g i.e., 0.13 g agreeing well with the ground acceleration recorded by an accelerograph at Ahmedabad (0.11 g). Hence, it is inferred that the Q₀ value obtained from this study seems to be apt for prediction of ground motion for the region.     

The <Emphasis Type="Italic">Q</Emphasis>-factor estimates for the crust and upper mantle in the vicinity of Sochi and Anapa (North Caucasus)

The regularities in the radiation and propagation of seismic waves in the regions of the North Caucasus are analyzed for estimating the ground motion parameters during the probable future strong earthquakes. Based on the records of the regional earthquakes with magnitudes M_W ~ 3.9–5.6 within epicentral distances up to ~300 km obtained during the period of digital measurements at the Sochi and Anapa seismic stations, the Q-factors in the vicinities of these sites are estimated at ~55 f^0.9 and ~90f^0.7, respectively. The estimates were obtained by the coda normalization method developed by Aki, Rautian, and other authors. This method is based on the phenomenon of suppression of the earthquake (source) effects and local (site) responses by coda waves in the S-wave spectra. The obtained Q-factor estimates can be used for forecasting the ground shaking parameters for the future probable strong earthquakes in the North Caucasus in the vicinities of Sochi and Anapa.     

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.     

Estimation of frequency dependent coda wave attenuation structure at the vicinity of Cairo Metropolitan Area

Estimation of seismic wave attenuation in the shallow crust in terms of coda wave Q structure previously investigated in the vicinity of Cairo Metropolitan Area was improved using seismograms of local earthquakes recorded by the Egyptian National Seismic Network. The seismic wave attenuation was measured from the time decay of coda wave amplitudes on narrow bandpass filtered seismograms based on the single scattering theory. The frequency bands of interest are from 1.5 to 18 Hz. In general, the values obtained for various events recorded at El-Fayoum and Wadi Hagul stations are very similar for all frequency bands. A regional attenuation law Q _c = 85.66 f ^0.79 was obtained.     

Characteristics of radiation and propagation of seismic waves in the Baikal Rift Zone estimated by simulations of acceleration time histories of the recorded earthquakes

The regularities in the radiation and propagation of seismic waves within the Baikal Rift Zone in Buryatia are studied to estimate the ground motion parameters from the probable future strong earthquakes. The regional parameters of seismic radiation and propagation are estimated by the stochastic simulation (which provides the closest agreement between the calculations and observations) of the acceleration time histories of the earthquakes recorded by the Ulan-Ude seismic station. The acceleration time histories of the strongest earthquakes (M _W ~ 3.4–4.8) that occurred in 2006–2011 at the epicentral distances of ~96–125 km and had source depths of ~8–12 km have been modeled. The calculations are conducted with estimates of the Q-factor which were previously obtained for the region. The frequency-dependent attenuation and geometrical spreading are estimated from the data on the deep structure of the crust and upper mantle (velocity sections) in the Ulan-Ude region, and the parameters determining the wave forms and duration of acceleration time histories are found by fitting. These parameters fairly well describe all the considered earthquakes. The Ulan-Ude station can be considered as the reference bedrock station with minimum local effects. The obtained estimates for the parameters of seismic radiation and propagation can be used for forecasting the ground motion from the future strong earthquakes and for constructing the seismic zoning maps for Buryatia.