新疆及邻区体波震级初步分析研究

依据《地震震级的规定（GB 17740—2017）》,分析了2009—2017年新疆地震台网所记录的新疆及邻区476次中深源地震事件,测定了13601个m_b（短周期体波震级）和12035个m_B（BB）（宽频带体波震级）的数据样本,回归分析m_b和m_B（BB）得到回归方程及量规函数,结果显示m_b和m_B（BB）相关系数为0.966,表明两者显著相关。因此,建议对于中强型中深源地震可以直接从原始速度型宽频带数字地震记录上测定长周期体波震级m_B（BB）,提高地震速报测定的速度和精度。通过震级偏差统计和台站场地响应计算,分析新疆地震台网中的XKR、HTA、ATS和KSZ等地震台站震级偏差较大的原因为砂岩、灰岩、砂土层等类型的台基放大了场地响应,说明台基类型对体波震级偏差的影响较大。与NEIC测定的体波震级对比时,发现新疆地震台网测定体波震级平均偏大0.42级,且偏差随着震源深度的增加有增大的趋势。     

Rediscovering Signal Complexity as a Teleseismic Discriminant

We re-examine the utility of teleseismic seismic complexity discriminants in a multivariate setting using United Kingdom array data. We measure a complexity discriminant taken on array beams by simply taking the logarithm of the ratio of the P-wave coda signal to that of the first arriving direct P wave (β_CF). The single station complexity discriminant shows marginal performance with shallow earthquakes having more complex signatures than those from explosions or deep earthquakes. Inclusion of secondary phases in the coda window can also degrade performance. However, performance improves markedly when two-station complexity discriminants are formed showing false alarm rates similar to those observed for network m_b − M_s. This suggests that multistation complexity discriminants may ameliorate some of the problems associated with m_b – M_s discrimination at lower magnitudes. Additionally, when complexity discriminants are combined with m_b – M_s there is a tendency for explosions, shallow earthquakes and deep earthquakes to form three distinct populations. Thus, complexity discriminants may follow a logic that is similar to m_b – M_s in terms of the separation of shallow earthquakes from nuclear explosions, although the underlying physics of the two discriminants is significantly different.     

Evaluation of magnitudes using USSR and US (NOS) techniques

Summary WWSSN seismograms for the Caucasus earthquake of May 14, 1970, and the Alma-Ata earthquake of June 5, 1970, were scaled form _b using techniques employed by USSR (P amplitude during first 20 seconds) and USA (P amplitude during first few cycles). USSR seismograms will be used in a second phase of this study. Using both techniques, the scatter at each station form _b (SPZ) was 0.0–0.3 with half of the stations having no scatter. The scatter among all of the stations was over 1 unit. Using both techniques, the scatter at each station form _b (LPZ) was zero except for NUR on May 14. The scatter among all of the stations was 1 unit on May 14 and 0.4 on June 5. This preliminary study indicates the scaling techniques have no effect in computingm _b (SPZ) andm _b (LPZ) but the differences are more likely related to source functions, crust and upper mantle response functions, and instrumental band pass.NOS Seismology personnel were transferred to NOAA-ESL in 1971.     

Teleseismic Lg of Semipalatinsk and Novaya Zemlya Nuclear Explosions Recorded by the GRF (Gräfenberg) Array: Comparison with Regional Lg (BRV) and their Potential for Accurate Yield Estimation

—?A comparison of regional and teleseismic log rms (root-mean-square) L g amplitude measurements have been made for 14 underground nuclear explosions from the East Kazakh test site recorded both by the BRV (Borovoye) station in Kazakhstan and the GRF (Gräfenberg) array in Germany. The log rms L g amplitudes observed at the BRV regional station at a distance of 690?km and at the teleseismic GRF array at a distance exceeding 4700?km show very similar relative values (standard deviation 0.048 magnitude units) for underground explosions of different sizes at the Shagan River test site. This result as well as the comparison of BRV rms L g magnitudes (which were calculated from the log rms amplitudes using an appropriate calibration) with magnitude determinations for P waves of global seismic networks (standard deviation 0.054 magnitude units) point to a high precision in estimating the relative source sizes of explosions from L g-based single station data. Similar results were also obtained by other investigators (Patton, 1988; Ringdal et?al., 1992) using L g data from different stations at different distances.¶Additionally, GRF log rms L g and P-coda amplitude measurements were made for a larger data set from Novaya Zemlya and East Kazakh explosions, which were supplemented with m _b (L g) amplitude measurements using a modified version of Nuttli's (1973, 1986a) method. From this test of the relative performance of the three different magnitude scales, it was found that the L g and P-coda based magnitudes performed equally well, whereas the modified Nuttli m _b (L g) magnitudes show greater scatter when compared to the worldwide m _b reference magnitudes. Whether this result indicates that the rms amplitude measurements are superior to the zero-to-peak amplitude measurement of a single cycle used for the modified Nuttli method, however, cannot be finally assessed, since the calculated m _b (L g) magnitudes are only preliminary until appropriate attenuation corrections are available for the specific path to GRF.     

<Emphasis Type="Italic">Lg</Emphasis> Body-Wave Magnitude Scaling for the Continental Margin around Korea and Japan

Regional body-wave magnitude scalings are essential for quantification of small and moderate-size earthquakes that are observed only up to regional distances. Crustally-guided shear waves, Lg, develop stably at regional distances in continental crusts and are minimally influenced by the source radiation patterns. Lg body-wave magnitude scalings, m_b(Lg),m_b(Lg), are widely used for assessment of sizes of regional crustal events. The m_b(Lg)m_b(Lg) scaling has rarely been tested in continental margins where Lg waves are significantly attenuated due to abrupt lateral variation of crustal structures. We test the applicability of m_b(Lg)m_b(Lg) scaling to the eastern margin of the Eurasian plate around the Korean Peninsula and Japanese islands. Both third-peak and root-mean-square (rms) amplitudes of Lg vary significantly according to the crustal structures along raypaths, causing apparent underestimation of m_b(Lg).m_b(Lg). Implementation of raypath-dependent quality factors (Q) allows accurate estimation of m_b(Lg),m_b(Lg), retaining the transportability of m_b(Lg)m_b(Lg) in the continental margin around Korea and Japan. The calibration constants for an rms-amplitude-based m_b(Lg)m_b(Lg) scaling are not determined to vary by region in the continental margin due to complicated crustal structures. The calibration constants are determined to be distance-dependent. Both the third-peak-amplitude-based and rms-amplitude-based m_b(Lg)m_b(Lg) scalings yield accurate magnitude estimates when raypath-dependent quality factors are implemented.     

The empirical relationships between M s, m b and M L for China and vicinity

Data from 753 earthquakes are used to determine a relationship between surface-wave magnitude (M _s) and bodywave magnitude (m _b), and from 541 earthquakes to determine a relationship between surface-wave magnitude (M _s) and local magnitude (M _L) for China and vicinity: M _s=0.9883 m _b-0.0420, M _s=0.9919 M _L-0.1773. The relationship of M _s versus m _b is obtained for 292 events occurred in the Chinese mainland in the time period from 1964 to 1996, 291 events occurred in Taiwan in the time period from 1964 to 1995 and 170 events occurred in the surrounding area. Standard deviation of the fitting is 0.445. Relationship of M _s versus M _L is obtained for 36 events occurred in the Chinese mainland, 293 events occurred in Taiwan, China and 212 events occurred in the surrounding area. The total amount is 541 events. Standard deviation of the fitting is 0.4673. The uncertainties of the converted M _s in different magnitude intervals can be estimated using complementary cumulative distribution function (CCDF). In the relationship of M _s versus m _b, taking ±0.25 as a range of uncertainties, in magnitude interval m _b 4.0–4.9, the probabilities for the converted M _s taken value less than (M _s-0.25) and more than (M _s+0.25) are 17% and 27% respectively. Similarly, we have probabilities for m _b 5.0–5.9 are 34% and 20% and that for m _b 6.0–6.9 are 11% and 47%. In the relationship of M _s versus M _L, if the range of uncertainties is still taken as ±0.25, the corresponding probabilities for magnitude interval M _L 4.0–4.9 are 22% and 38%, for M _L 5.0–5.9 are 20% and 15% and for magnitude interval M _L 6.0–6.9, are 15% and 29%, respectively. The relationships developed in this paper can be used for the conversion of one magnitude scale into another magnitude scales conveniently. The estimation of uncertainties described in this paper is more accurate and more objective than the usual estimation expressed by deviation. The estimations described in this paper indicate various dispersions in different magnitude intervals of original data. The estimations of uncertainties described by probabilities can be well connected with the total estimations of uncertainties in seismic hazard assessment.     

Source parameters of some presumed Semipalatinsk underground nuclear explosions

Summary Short-period vertical-componentP-wave spectra of seven presumed Semipalatinsk underground nuclear explosions, recorded by the Swedish seismic station network, are investigated. The events considered have closely spaced foci and cover the magnitude range fromm _b=5.5 tom _b=6.6. Spectra of six of these explosions show pronounced minima, varying from about 1.5 to 1.8 cps, which could be explained as principle minima due toP-pP interference. Supposing a nearsurfaceP-wave velocity at the test area of 4 km/sec, the shot depths are estimated to vary roughly from 750 to 1350 m. In order to obtain an estimate of the yield, the observed spectra are compared withHaskell's theoretical source spectra. For four events, relative yield estimates fit well the predicted values for explosions fired in a granitic medium. The behaviour of the remaining three explosions is discussed in detail.     

Comparison between different earthquake magnitudes determined by China Seismograph Network

By linear regression and orthogonal regression methods, comparisons are made between different magnitudes (lo-cal magnitude ML, surface wave magnitudes MS and MS7, long-period body wave magnitude mB and short-period body wave magnitude mb) determined by Institute of Geophysics, China Earthquake Administration, on the basis of observation data collected by China Seismograph Network between 1983 and 2004. Empirical relations between different magnitudes have been obtained. The result shows that: 1 As different magnitude scales reflect radiated energy by seismic waves within different periods, earthquake magnitudes can be described more objectively by using different scales for earthquakes of different magnitudes. When the epicentral distance is less than 1 000 km, local magnitude ML can be a preferable scale; In case M<4.5, there is little difference between the magnitude scales; In case 4.5MS, i.e., MS underestimates magnitudes of such events, therefore, mB can be a better choice; In case M>6.0, MS>mB>mb, both mB and mb underestimate the magnitudes, so MS is a preferable scale for deter-mining magnitudes of such events (6.08.5, a saturation phenomenon appears in MS, which cannot give an accurate reflection of the magnitudes of such large events; 2 In China, when the epicentral distance is less than 1 000 km, there is almost no difference between ML and MS, and thus there is no need to convert be-tween the two magnitudes in practice; 3 Although MS and MS7 are both surface wave magnitudes, MS is in general greater than MS7 by 0.2~0.3 magnitude, because different instruments and calculation formulae are used; 4 mB is almost equal to mb for earthquakes around mB4.0, but mB is larger than mb for those of mB≥4.5, because the periods of seismic waves used for measuring mB and mb are different though the calculation formulae are the same.     

P-wave amplitude spectra of Nevada underground nuclear explosions

Summary Short-period vertical-componentP-wave spectra have been determined for twelve Nevada underground nuclear explosions recorded by the Swedish seismic station network. Selected events cover the magnitude range fromm _b =5.9 tom _b =7.0 and the shot-depth range from 460 m to 1400 m. All explosion spectra show pronounced minima near 1 sec period. Within individual test areas the period at which the minimum occurs increases systematically with increasing shot-depth. This dependence makes it possible to interpret the observed minima in terms of interference betweenP- andpP-waves. One atmospheric explosion from Novaya Zemlya is analyzed to verify the inferred minima intrpretation.The minimum-period shows also a systematic linear dependence upon the magnitude. However, the physical explanation for this dependence is not evident and it may just be a consequence of a magnitude-depth relation.On leave from the Geophysical Institute, Czechosl. Acad. Sci., Prague 4-Spoilov, Boni II.     

Empirical Global Relations Converting M S and m b to Moment Magnitude

The existence of several magnitude scales used by seismological centers all over the world and the compilation of earthquake catalogs by many authors have rendered globally valid relations connecting magnitude scales a necessity. This would allow the creation of a homogeneous global earthquake catalog, a useful tool for earthquake research. Of special interest is the definition of global relations converting different magnitude scales to the most reliable and useful scale of magnitude, the moment magnitude, M _W. In order to accomplish this, a very large sample of data from international seismological sources (ISC, NEIC, HRVD, etc.) has been collected and processed. The magnitude scales tested against M _W are the surface wave magnitude, M _S, the body wave magnitude, m _b, and the local magnitude, M _L. The moment magnitudes adopted have been taken from the CMT solutions of HRVD and USGS. The data set used in this study contains 20,407 earthquakes, which occurred all over the world during the time period 1.1.1976–31.5.2003, for which moment magnitudes are available. It is shown that well-defined relations hold between M _W and m _b and M _S and that these relations can be reliably used for compiling homogeneous, with respect to magnitude, earthquake catalogs.     

Sources of Error and the Statistical Formulation of M S: m b Seismic Event Screening Analysis

The Comprehensive Nuclear-Test-Ban Treaty (CTBT), a global ban on nuclear explosions, is currently in a ratification phase. Under the CTBT, an International Monitoring System (IMS) of seismic, hydroacoustic, infrasonic and radionuclide sensors is operational, and the data from the IMS is analysed by the International Data Centre (IDC). The IDC provides CTBT signatories basic seismic event parameters and a screening analysis indicating whether an event exhibits explosion characteristics (for example, shallow depth). An important component of the screening analysis is a statistical test of the null hypothesis H ₀: explosion characteristics using empirical measurements of seismic energy (magnitudes). The established magnitude used for event size is the body-wave magnitude (denoted m _b) computed from the initial segment of a seismic waveform. IDC screening analysis is applied to events with m _b greater than 3.5. The Rayleigh wave magnitude (denoted M _S) is a measure of later arriving surface wave energy. Magnitudes are measurements of seismic energy that include adjustments (physical correction model) for path and distance effects between event and station. Relative to m _b, earthquakes generally have a larger M _S magnitude than explosions. This article proposes a hypothesis test (screening analysis) using M _S and m _b that expressly accounts for physical correction model inadequacy in the standard error of the test statistic. With this hypothesis test formulation, the 2009 Democratic Peoples Republic of Korea announced nuclear weapon test fails to reject the null hypothesis H ₀: explosion characteristics.     

Comparison between earthquake magnitudes determined by China seismograph network and US seismograph networks (I): Body wave magnitude

Introduction Gutenberg (1945a, b) introduced body wave magnitude based on P, PP and S waves (with a period of 0.5~12.0 s) of teleseismic events. Body wave magnitude includes mb determined with short-period seismograph and mB determined with middle- and long-period seismographs. Some-times it is written as m, which is referred to as unified earthquake magnitude. mb represents earth-quake magnitude measured with body wave amplitude around 1 s, while mB represents earthquake magnitude measured …     

Optimized Seismic Threshold Monitoring – Part 1: Regional Processing

v--v Continuous seismic threshold monitoring is a technique that has been developed over the past several years to assess the upper magnitude limit of possible seismic events that might have occurred in a geographical target area. The method provides continuous time monitoring at a given confidence level, and can be applied in a site-specific, regional or global context.¶In this paper (Part 1) and a companion paper (Part 2) we address the problem of optimizing the site-specific approach in order to achieve the highest possible automatic monitoring capability of particularly interesting areas. The present paper addresses the application of the method to cases where a regional monitoring network is available. We have in particular analyzed events from the region around the Novaya Zemlya nuclear test site to develop a set of optimized processing parameters for the arrays SPITS, ARCES, FINES, and NORES. From analysis of the calibration events we have derived values for beam-forming steering delays, filter bands, short-term average (STA) lengths, phase travel times (P and S waves), and amplitude-magnitude relationships for each array. By using these parameters for threshold monitoring of the Novaya Zemlya testing area, we obtain a monitoring capability varying between m_b 2.0 and 2.5 during normal noise conditions.¶The advantage of using a network, rather than a single station or array, for monitoring purposes becomes particularly evident during intervals with high global seismic activity (aftershock sequences), high seismic noise levels (wind, water waves, ice cracks) or station outages. For the time period November-December 1997, all time intervals with network magnitude thresholds exceeding m_b 2.5 were visually analyzed, and we found that all of these threshold peaks could be explained by teleseismic, regional, or local signals from events outside the Novaya Zemlya testing area. We could therefore conclude within the confidence level provided by the method, that no seismic event of magnitude exceeding 2.5 occurred at the Novaya Zemlya test site during this two-month time interval.¶As an example of particular interest in a monitoring context, we apply optimized threshold processing of the SPITS array for a time interval around 16 August 1997 m_b 3.5 event in the Kara Sea. We show that this processing enables us to detect a second, smaller event from the same site (m_b 2.6), occurring about 4 hours later. This second event was not defined automatically by standard processing.     

Spectral characteristics of natural and artificial seismic events in the Lop Nor test site, China

Seismic discriminants based on the spectral seismogram and spectral magnitude techniques have been tested to discriminate between three events; a nuclear explosion which took place in Lop Nor, China with m _b 6.1 and two earthquakes from the closest area with m _b 5.5 and 5.3, respectively. The spectral seismogram of the three events shows that the frequency content of the nuclear explosion differs from that of the earthquakes where the P wave is richier in high frequency content in the nuclear explosion than the corresponding earthquakes. It is also observed that the energy decays more rapidly for the nuclear explosion than for the earthquakes. Furthermore, the spectral magnitudes reveal significant differences in the spectra between the nuclear explosion and the two earthquakes. The observed differences appear to be quite enough to provide a reliable discriminant. The estimated stress drop from the magnitude spectra indicates a higher stress drop of the nuclear explosion relative to the earthquakes of the same tectonic region.     

Calibration of magnitude scales for earthquakes of the Mediterranean

In order to provide the tools for uniform size determination for Mediterranean earthquakes over the last 50-year period of instrumental seismology, we have regressed the magnitude determinations for 220 earthquakes of the European-Mediterranean region over the 1977–1991 period, reported by three international centres, 11 national and regional networks and 101 individual stations and observatories, using seismic moments from the Harvard CMTs. We calibrate M(M₀) regression curves for the magnitude scales commonly used for Mediterranean earthquakes (M_L, M_WA, m_b, M_S, M_LH, M_LV, M_D, M); we also calibrate static corrections or specific regressions for individual observatories and we verify the reliability of the reports of different organizations and observatories. Our analysis shows that the teleseismic magnitudes (m_b, M_S) computed by international centers (ISC, NEIC) provide good measures of earthquake size, with low standard deviations (0.17–0.23), allowing one to regress stable regional calibrations with respect to the seismic moment and to correct systematic biases such as the hypocentral depth for M_S and the radiation pattern for m_b; while m_b is commonly reputed to be an inadequate measure of earthquake size, we find that the ISC m_b is still today the most precise measure to use to regress M_W and M₀ for earthquakes of the European-Mediterranean region; few individual observatories report teleseismic magnitudes requiring specific dynamic calibrations (BJI, MOS). Regional surface-wave magnitudes (M_LV, M_LH) reported in Eastern Europe generally provide reliable measures of earthquake size, with standard deviations often in the 0.25–0.35 range; the introduction of a small (±0.1–0.2) static station correction is sometimes required. While the Richter magnitude M_L is the measure of earthquake size most commonly reported in the press whenever an earthquake strikes, we find that M_L has not been computed in the European-Mediterranean in the last 15 years; the reported local magnitudes M_WA and M_L do not conform to the Richter formula and are of poor quality and little use, with few exceptions requiring ad hoc calibrations similar to the M_S regression (EMSC, ATH). The duration magnitude M_D used by most seismic networks confirms that its use requires accurate station calibrations and should be restricted only to events with low seismic moments.     

m b(Pn) Scale for the Korean Peninsula and Site-Dependent Pn Amplification

The Korean Peninsula is located in the far-eastern Eurasian plate margin where crustal structures vary laterally, causing significant raypath-dependent modulations of crustal phases. The discriminative variations of crustal phases hinder application of conventional local magnitude scales in the continental margin. The mantle-lid phase is less affected by the crustal structures than the crustal phases, providing a better constraint to magnitude estimation. A regional body-wave magnitude scale based on the mantle-lid P wave (Pn), m _b(Pn),?is developed for regional events around the Korean Peninsula. The m _b(Pn) scale is determined to be m _b(Pn)?=?0.380 (±0.299)?+?log A?+?2.012 (±0.122) log d, where A is the peak-to-peak Pn amplitude in?μm and d is the epicentral distance in km. The m _b(Pn) estimates of regional events around the Korean Peninsula are determined. The m _b(Pn) estimates are compared with other available magnitude estimates (m _b(Lg),?M _L). The influence of structures beneath stations on Pn amplification is investigated from inter-station magnitude residuals. A characteristic spatial variation of inter-station magnitude residuals with strengths mostly between ?6 and 6 %, but with maximum strengths of?±10 %, is observed. The inter-station magnitude residuals appears to be correlated well with geological and seismic structures in the crust.     

Recent seismicity in Northeast India and its adjoining region

Recent seismicity in the northeast India and its adjoining region exhibits different earthquake mechanisms – predominantly thrust faulting on the eastern boundary, normal faulting in the upper Himalaya, and strike slip in the remaining areas. A homogenized catalogue in moment magnitude, M _W, covering a period from 1906 to 2006 is derived from International Seismological Center (ISC) catalogue, and Global Centroid Moment Tensor (GCMT) database. Owing to significant and stable earthquake recordings as seen from 1964 onwards, the seismicity in the region is analyzed for the period with spatial distribution of magnitude of completeness m _t, b value, a value, and correlation fractal dimension D _C. The estimated value of m _t is found to vary between 4.0 and 4.8. The a value is seen to vary from 4.47 to 8.59 while b value ranges from 0.61 to 1.36. Thrust zones are seen to exhibit predominantly lower b value distribution while strike-slip and normal faulting regimes are associated with moderate to higher b value distribution. D _C is found to vary from 0.70 to 1.66. Although the correlation between spatial distribution of b value and D _C is seen predominantly negative, positive correlations can also be observed in some parts of this territory. A major observation is the strikingly negative correlation with low b value in the eastern boundary thrust region implying a possible case of extending asperity. Incidentally, application of box counting method on fault segments of the study region indicates comparatively higher fractal dimension, D, suggesting an inclination towards a planar geometrical coverage in the 2D spatial extent. Finally, four broad seismic source zones are demarcated based on the estimated spatial seismicity patterns in collaboration with the underlying active fault networks. The present work appraises the seismicity scenario in fulfillment of a basic groundwork for seismic hazard assessment in this earthquake province of the country.     

Research on characteristics of magnitude structure of earthquake sequences

ＲｅｓｅａｒｃｈｏｎｃｈａｒａｃｔｅｒｉｓｔｉｃｓｏｆｍａｇｎｉｔｕｄｅｓｔｒｕｃｔｕｒｅｏｆｅａｒｔｈｑｕａｋｅｓｅｑｕｅｎｃｅｓＰＥＩ－ＱＩＮＧＳＵＮ（孙佩卿）ＱＩＮ－ＺＵＬＩ（李钦祖）ＹＩＮＧ－ＨＵＡＤＡＩ（戴英华）ＪＵＮＺＨＡＯ（赵军）Ｓｅ...     

传统与现代震级标度

首先阐述并比较传统震级(M_L、m_b、m_B、M_S及M_JMA)所用数据、计算方法及适用范围,系统梳理两种现代震级标度M_e和M_W的理论基础及计算方法。然后根据使用数据及适用场景不同,介绍M_WW、M_WC、M_Wb、M_wp及M_dt等矩震级测定方法;同时,随着越来越丰富的GPS观测数据,基于高频GPS观测也越来越多地被用于地震震级快速测定工作。最后对这些震级测定方法的时效性进行对比分析,讨论不同震级标度的适用场景及其稳定性。     

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.