Average crustal structure of Sweden

Summary Records obtained at the permanent stations of the Swedish seismograph network from explosions carried out in Scandinavian waters in June 1969 are evaluated. The study includes determination of velocities for all crustal phases observed, furthermore of layer thicknesses, Poisson ratios and amplitude ratios. The purpose of the study is partly to provide a first approximation to the crustal structure in Sweden, partly to provide regional data for location of earthquakes and explosions in the area in the future. Average velocities (km/sec) are forPn 7.88±0.05,Pg1 6.25±0.08,Pg2 5.70,Sn 4.58±0.04,S ^* 3.70±0.04,Sg1 (Lg1) 3.58±0.03,Sg2 (Sg) 3.40±0.03,Rg 3.02±0.07. The average thickness is 12 km for the granitic layer, and 23 km for the basaltic layer, thus making the average crustal thickness equal to 35 km. Relative amplitudes plotted versus distance complete the dynamical side of the study and they are useful for identification of waves. A regional travel-time table is presented for the distance range 0°–10° with entries for each 0.1° and including all crustal phases read.     

A crust-mantle model for the NORSAR area

Summary Elastic waves from explosions were recorded at NORSAR and at a number of field stations, and the data were used for determining a crust-mantle model under the array. The number of explosions was eleven distributed on seven shot points. The total number of recording points was fifty-one, and the interpretation was based on 350 individual records.The velocities obtained for the crustal phases were 6.2, 6.6 and 8.2 km/sec for theP _g ,P _g andP _n waves respectively. A deep crustal phase with a velocity of about 7.4 km/sec was observed. The mean depths to the discontinuities within the crust were determined to be 17 and 26 km. The depth to Moho varied greatly across the array from 31.5 km in the central part to 38 km under the C-ring. The maximum dip observed for the Moho was 12^o.Contribution No. 57 to Norwegian Geotraverse Project.     

On the seismic phasesP n andS n as recorded in Fennoscandia

Summary P _n velocities determined from seismic refraction measurements, show significant differences between Southern Finland (7.96 km/sec) and Southeastern Norway (8.20 km/sec). TheP _n/S_n velocity ratios (k) were determined from earthquake and explosion data, and the observed variation ofk indicates lateral variations in theP _n and/orS _n velocities in Fennoscandia.     

Ray parameters of diffracted long period P and S waves and the velocities at the base of the mantle

The derivation of P and S velocities at the core-mantle boundary (CMB) from long-period diffracted waves by the use of the simple ray-theoretical formulav _CMB=r _c /p (v _CMB=velocity at the CMB;r _c =core radius;p=ray parameter) yields apparent velocity values which differ from the true velocities. Using a dominant period of about 20 sec for calculating theoretical seismograms, we found a linear relation between the apparent velocity and the average velocity in a transition zone at the base of the mantle with fixed velocity on top.The ray parameters determined from long-period earthquake data are found to be 4.540±0.035 and 8.427±0.072 sec/deg for P_diff and S_diff, respectively. These values yield apparent velocities of 13.378±0.103 for P and 7.207±0.062 km/sec for S waves. By means of the theoretical relation between apparent and average velocity and under the assumption of linear variation of velocity with depth, one can invert the apparent velocities into true CMB velocities of 13.736±0.170 and 7.320±0.124 km/sec. These results imply positive velocity gradients at the base of the mantle and hence no significant departures from adiabaticity and homogeneity.Contribution No. 211 of the Geophysical Institute, University of Karlsruhe.     

A continental channel wave,guided by the intermediate layer in the crust

Summary A study has been made of a new channel wave, denotedLi, using a total of 83 observations from the seismic records of Swedish stations, mainly from earthquakes at normal depth.Li resembles theLg waves in several respects: it propagates only through continental structures, it has a similar particle motion, i.e. mainly transverse horizontal, and only slightly larger period. ButLi has a higher velocity, 3.79±0.07 km/sec, and it is believed to propagate in the intermediate layer in the crust in a way similar to the propagation of theLg waves in the granitic layer.Li is identical withS ^* in records of near-by earthquakes in the same way asLg2 is identical withSg. Li usually exhibits no clear dispersion.

---

Zusammenfassung Es wurden Untersuchungen angestellt über eine neue Kanalwelle, welcher die BezeichnungLi gegeben wurde, wobei insgesamt 83 Erdbebenregistrierungen von schwedischen Stationen Verwendung fanden. Hauptsächlich waren es Erdbeben mit normaler Herdtiefe. DieLi-Wellen haben in verschiedener Hinsicht Ähnlichkeit mit denLg-Wellen: Sie pflanzen sich nur im Bereich kontinentaler Struktur fort und sie haben eine ähnliche Partikelbewegung, d.h. hauptsächlich horizontaltransversal.Li hat eine nur unwesentlich höhere Periode als dieLg-Wellen. AberLi hat eine wesentlich höhere Geschwindigkeit, 3.79±0.07 km/sec, und es ist anzunehmen, dass sie sich in der Basaltschicht der Kruste in ähnlicher Weise fortpflanzt wieLg in der Granitschicht.Li ist identisch mitS ^* in Aufzeichnungen von Nahbeben, so wieLg2 identisch mitSg ist.Li weist gewöhnlich keine deutliche Dispersion auf.

     

A neural-tree-based system for automatic location of earthquakes in Northeastern Italy

A neural network system for P and S-picking and location of earthquakes in Northeastern Italy is described. It is applied to 7108 seismograms corresponding to 1147 earthquakes occurring in Northeastern Italy and surrounding area in the period 2000–2003. Its results are compared with two sets of manual picks and with the picks performed by the existing seismic alert system. The new system recognizes 89% and 67% of P and S arrival times, respectively, which allows locating 92% of the earthquakes. P and S-picks differ from the best available manual picks by 0.00 ± 0.07 s and 0.00 ± 0.18 s, respectively. The corresponding earthquake locations differ by −0.18± 0.77 km in longitude, 0.10± 0.62 km in latitude and 0.1± 2.0 km in depth. These results suggest its use for alert purposes and rapid production of preliminary bulletins.Considering a subset of picks that are common to all the available data sets, the absolute accuracy (i.e., the inverse of the standard deviation of differences between the estimated and the true, unknown arrival times) of each picking method is estimated. The best available manual data set has standard deviation 0.03 s for P waves and 0.07 s for S waves, while for the new system it is 0.06 s and 0.18 s for P and S waves, respectively.This paper has not been submitted elsewhere in identical or similar form, nor will it be during the first three months after its submission to Journal of Seismology.     

Body-wave Attenuation in the Region of Garda, Italy

We analyzed the spectral amplitude decay with hypocentral distance of P and S waves generated by 76 small magnitude earthquakes (M_L 0.9–3.8) located in the Garda region, Central-Eastern Alps, Italy. These events were recorded by 18 stations with velocity sensors, in a distance range between 8 and 120 km. We calculated nonparametric attenuation functions (NAF) and estimated the quality factor Q of both body waves at 17 different frequencies between 2 and 25 Hz. Assuming a homogeneous model we found that the Q frequency dependence of P and S can be approximated with the functions Q _P = 65 f ^0.9 and Q _S = 160 f ^0.6 , respectively. At 2 Hz the Q _S /Q _P ratio reaches the highest value of 2.8. At higher frequencies Q _S /Q _P varies between 0.7 and 1.7, suggesting that for this frequency band scattering may be an important attenuation mechanism in the region of Garda. To explore the variation of Q in depth, we estimated Q at short (r ≤ 30 km) and intermediate (35–90 km) distance paths. We found that in the shallow crust P waves attenuate more than S (1.3 < Q _S /Q _P < 2.5). Moreover, P waves traveling along paths in the lower crust (depths approximately greater than 30 km) attenuate more than S waves. To quantify the observed variability of Q in depth we considered a three-layer model and inverted the NAF to estimate Q in each layer. We found that in the crust Q increases with depth. However, in the upper mantle (~40–50 km depth) Q decreases and in particular the high frequency Q _S (f > 9 Hz) has values similar to those estimated for the shallow layer of the crust.     

A study of the crustal structure beneath the Himalayas from body waves

Summary The crustal structure beneath the Himalayas has been investigated using body wave data from near earthquakes having epicentres over the Himalayas and recorded by the observatories situated over, or very near, the foothills of the mountains. A three-layered crustal model, without the top sedimentary layer, with velocities for theP wave group in Granite I, Granite II and the Basaltic layer as 5.48, 6.00 and 6.45 and for theS wave group as 3.33, 3.56 and 3.90 km/sec respectively, has been interpreted. The upper mantle velocity for theP wave has been observed to be 8.07 km/sec and for theS wave as 4.57 km/sec. Average thickness for the Granite I layer has been computed as 22.7 km, for the Granite II layer as 16.3 km and for the Basaltic layer as 18.7 km. Crustal and sub-crustal velocities indicate a lower trend under the mountain. A thicker crust has been obtained beneath the Himalayas.     

Attenuation of High Frequency P and S Waves in the Gujarat Region,India

The local earthquake waveforms recorded on broadband seismograph network of Institute of Seismological Research in Gujarat, India have been analyzed to understand the attenuation of high frequency (2–25 Hz) P and S waves in the region. The frequency dependent relationships for quality factors for P (Q _P) and S (Q _S) waves have been obtained using the spectral ratio method for three regions namely, Kachchh, Saurashtra and Mainland Gujarat. The earthquakes recorded at nine stations of Kachchh, five stations of Saurashtra and one station in mainland Gujarat have been used for this analysis. The estimated relations for average Q _P and Q _S are: Q _P = (105 ± 2) f ^{0.82 ± 0.01}, Q _S = (74 ± 2) f ^{1.06 ± 0.01} for Kachchh region; Q _P = (148 ± 2) f ^{0.92 ± 0.01}, Q _S = (149 ± 14) f ^{1.43 ± 0.05} for Saurashtra region and Q _P = (163 ± 7) f ^{0.77 ± 0.03}, Q _S = (118 ± 34) f ^{0.65 ± 0.14} for mainland Gujarat region. The low Q (<200) and high exponent of f (>0.5) as obtained from present analysis indicate the predominant seismic activities in the region. The lowest Q values obtained for the Kachchh region implies that the area is relatively more attenuative and heterogeneous than other two regions. A comparison between Q _S estimated in this study and coda Q (Qc) previously reported by others for Kachchh region shows that Q _C > Q _S for the frequency range of interest showing the enrichment of coda waves and the importance of scattering attenuation to the attenuation of S waves in the Kachchh region infested with faults and fractures. The Q _S/Q _P ratio is found to be less than 1 for Kachchh and Mainland Gujarat regions and close to unity for Saurashtra region. This reflects the difference in the geological composition of rocks in the regions. The frequency dependent relations developed in this study could be used for the estimation of earthquake source parameters as well as for simulating the strong earthquake ground motions in the region.     

Travel-times and Attenuation Relations for Regional Phases in the Barents Sea Region

A database containing 45 events in the Barents Sea region has been compiled and analyzed with the aim of evaluating crustal models, travel-times and attenuation relations in the context of performing regional detection threshold monitoring of this region. The 45 events are mostly located around the circumference of the study area due to the virtually aseismic nature of the Barents Sea itself. Regional P _n and S _n phases were observable for most events in the database, while P _g and L _g phases were only observable for events with raypaths that do not cross the tectonic structures in the Barents Sea. This corroborates a number of previous observations of L _g -wave blockage within the Barents Sea. Three existing velocity models were evaluated, with a model having slightly lower S velocities than earlier assumed in the upper mantle giving the overall best fit to the observed arrivals. In order to estimate magnitudes, short-term average (STA) and spectral amplitude values were calculated in several frequency bands for all phase arrivals in the database. There were no significant differences between spectral and STA amplitudes, so the latter were used as this parameter is more efficient to calculate in real-time processing. An inversion was performed in order to determine an attenuation relation specific for this region. The resulting magnitudes based on P _n , P _g , S _n and L _g phases gave an internally consistent, reasonably stable set of values, which can be calibrated towards any existing global or regional scale.     

Fe-Mg interdiffusion in magnesiowüstite up to 35 GPa

Fe-Mg interdiffusivities in (Fe,Mg)O magnesiowüstite have been measured in experiments conducted at pressures of 7-35 GPa and temperatures of 1573-1973 K using a Kawai-type high-pressure apparatus. The diffusion profiles were measured across the interface between MgO and (Fe_0.5,Mg_0.5)O samples by electron microprobe analysis, and the Fe-Mg interdiffusivities were determined as D_Fe-Mg=D₀exp{−E*(1+PV*_Mg/E*_Mg)/RT}, where D₀=4.1(^+16.1_−3.3)×10⁻⁷ m²/s, E*=(1−C_Mg)E*_Fe+C_MgE*_Mg (activation energy for the concentration of Mg, where E*_Fe=113(±74) kJ/mol and E*_Mg=226(±32) kJ/mol), the activation volume V*_Mg=1.8(±1.2)×10⁻⁶ m³/mol. By extrapolating these results to the P-T conditions of the core-mantle boundary, we conclude that the interdiffusivity of Fe-Mg in magnesiowüstite at the bottom of the lower mantle is at least one order of magnitude larger than that at the top of the lower mantle.     

Attenuation of coda waves in northern Greece

The single scattering model has been applied for the estimation of codaQ values for local earthquakes that occurred in northern Greece during the period 1983–1989 and recorded by the telemetered network of the Geophysical Laboratory of the University of Thessaloniki. CodaQ estimations were made for four frequency bands centered at 1.5 Hz, 3.0 Hz, 6.0 Hz and 12.0 Hz and for the lapse time windows 10–20 sec, 15–30 sec, 20–45 sec, 30–60 sec and 50–100 sec. The codaQ values obtained show a clear frequency dependence of the formQ _c =Q ₀ f ⁿ , whileQ ₀ andn depend on the lapse time window.Q ₀ was found equal to 33 andn equal to 1.01 for the time window of 10 to 20 sec, while for the other windowsQ ₀ increased from 60 to 129, withn being stable, close to 0.75. This lapse time dependence is interpreted as due to a depth dependent attenuation. The high attenuation and the strong frequency dependence found are characteristic of an area with high seismicity, in agreement with studies in other seismic regions.     

Receiver function study in northern Sumatra and the Malaysian peninsula

In this receiver function study, we investigate the structure of the crust beneath six seismic broadband stations close to the Sunda Arc formed by subduction of the Indo-Australian under the Sunda plate. We apply three different methods to analyse receiver functions at single stations. A recently developed algorithm determines absolute shear-wave velocities from observed frequency-dependent apparent incidence angles of P waves. Using waveform inversion of receiver functions and a modified Zhu and Kanamori algorithm, properties of discontinuities such as depth, velocity contrast, and sharpness are determined. The combination of the methods leads to robust results. The approach is validated by synthetic tests. Stations located on Malaysia show high-shear-wave velocities (V _S) near the surface in the range of 3.4–3.6 km s^− 1 attributed to crystalline rocks and 3.6–4.0 km s^− 1 in the lower crust. Upper and lower crust are clearly separated, the Moho is found at normal depths of 30–34 km where it forms a sharp discontinuity at station KUM or a gradient at stations IPM and KOM. For stations close to the subduction zone (BSI, GSI and PSI) complexity within the crust is high. Near the surface low V _S of 2.6–2.9 km s^− 1 indicate sediment layers. High V _S of 4.2 km s^− 1 are found at depth greater than 6 and 2 km at BSI and PSI, respectively. There, the Moho is located at 37 and 40 km depth. At station GSI, situated closest to the trench, the subducting slab is imaged as a north-east dipping structure separated from the sediment layer by a 10 km wide gradient in V _S between 10 and 20 km depth. Within the subducting slab V _S ≈ 4.7 km s^− 1. At station BSI, the subducting slab is found at depth between 90 and 110 km dipping 20° ± 8° in approximately N 60° E. A velocity increase in similar depth is indicated at station PSI, however no evidence for a dipping layer is found.     

Introductory notes on shock remanent magnetization and shock demagnetization of igneous rocks

Summary Effects of mechanical shocks of about 0.5 msec in duration on the remanent magnetization of igneous rocks are experimentally studied. The remanent magnetization acquired by applying a shock (S) in the presence of a magnetic field (H), which is symbolically expressed asJ _R (H₊S H_o), is very large compared with the ordinary isothermal remanent magnetization (IRM) acquired in the same magnetic field.J _R (H₊S H_o) is proportional to the piezo-remanent magnetization,J _R (H₊P₊P_o H_o).The effect of applyingS in advance of an acquisition of IRM is represented symbolically byJ _R (S H₊ H_o).J _R (S H₊ H_o) can become much larger than the ordinary IRM, and is proportional to the advance effect of pressure on IRM,J _R(P₊ P₀ H₊ H₀).The effect of shockS applied on IRM in non-magnetic space is represented by the shock-demagnetization effect,J _R(H₊ H₀ S), which also is proportional toJ _R(H₊ H₀ P₊ P₀).Because, the duration of a shock is very short, a single shock effect cannot achieve the final steady state. The effect ofn-time repeated shocks, is represented byJ ₀+J ^*(n), whereJ ₀ means the immediate effect and J ^*(n) represent the resultant effect of repeating, which is of mathematical expression proportional to [1–exp {–(n–1)}].

---

Zusammenfassung Die Effekte des mechanischen Stosses mit der Dauer von etwa 0.5 ms auf der remanenten Magnetisierung wurden experimentell nachgesucht. Das erworbene Remanenz der Magnetisierung nach dem Stoss (S) unter dem magnetischen Feld (H), das hier symbolisch alsJ _R(H₊ SH₀) bezechnet wird, ist sehr stark im Vergleich mit der normalen isothermischen remanenten Magnetisierung (IRM) unter demselben magnetischen Feld.J _R(H₊ S H₀) ist im Verhältnis zur piezoremanenten Magnetisierung,J _R(H₊ P₊ P₀ H₀).Der Effekt vom Stoss vor der Erwerbung von IRM wird symbolisch alsJ _R(S H₊ H₀) bezeichnet.J _R(S H₊ H₀) kann viel stärker als die normale IRM werden, im verhältnis zum Effekt des vorausgegebenen Drucks auf IRMJ _R(P₊ P₀ H₊ H₀).Der Effekt des Stosses auf IRM im Raum ohne magnetisches Feld wird mit dem Stossentmagnetisierungseffekt dargestellt,J _R(H₊ H₀ S), der auch proportional zuJ _R(H₊ H₀ P₊ P₀) ist.Da die Dauer einzelnen Stosses sehr kurz ist, kann der Effekt des einmaligen Stosses den endgültigen stabilen Zustand nicht erreichen. Der Effekt nachn-maligen wiederholten Stossen wird alsJ ₀+J ^*(n) bezeichnet, wobeiJ ₀ den unverzüglichen Effekt bedeutet, und J ^*(n) beschreibt den resultanten Effekt der Stosswiederholung, dessen mathematische Darstellung proporational zu [1–exp {–(n–1)}] ist.

     

A scattering model for seismic attenuation and its global applications

The attenuation characteristics of Indian lithosphere and its comparison with different tectonic settings in the world are determined from the observations of the Q for L_g(Q_Lg)-, and S(Q_S)-waves in the 1-30 Hz frequency range. The scattering is approximated with a Gaussian distribution of spherical scatterers. To approximate single scattering, we use Dainty's [Geophy. Res. Lett. 8 (11) (1981) 1126] model that attenuation is given by 1/Q(ω) = 1/Q_i + g(ω)v/ω, where Q_i is intrinsic Q due to anelastic attenuation, v is shear wave velocity, ω is angular frequency, g = ∫n(a)σ da is the total scattering coefficient for S-to-S scattering, n(a) da is the number of scattering spheres of radius a per unit volume, and σ is the scattering cross-section for the sphere. We find that if n(a) is described by a simple two parameter (a₀ and c) Gaussian of amplitude c and standard deviation and mean a₀, the attenuation data for different regions of the world are well approximated over the frequency band of seismic observations. Our major findings are: (1) the maximum effect of scattering on attenuation occurs at 0.84 Hz or a wavelength of 4.16 km; (2) the values of g are frequency dependent. Values of g are of the order of 10⁻³ km⁻¹ at 1-30 Hz, varying from 0.0031 to 0.01 and 0.001 to 0.0083 km⁻¹ for tectonically active and stable regions, respectively; (3) regions of active tectonics and seismicity generally have lower Q_i values (1000) than that in stable regions (2000); and (4) regions of high Q_i value exhibit low intensity of scattering.     

Coda Q in the Kachchh Basin, Western India Using Aftershocks of the Bhuj Earthquake of January 26, 2001

Q_C-estimates of Kachchh Basin in western India have been obtained in a high frequency range from 1.5 to 24.0 Hz using the aftershock data of Bhuj earthquake of January 26, 2001 recorded within an epicentral distance of 80 km. The decay of coda waves of 30 sec window from 186 seismograms has been analysed in four lapse time windows, adopting the single backscattering model. The study shows that Q_c is a function of frequency and increases as frequency increases. The frequency dependent Q_c relations obtained for four lapse-time windows are: Q_c=82 f^1.17 (20–50 sec), Q_c=106 f^1.11 (30–60 sec), Q_c=126f^1.03 (40–70 sec) and Q_c=122f^1.02 (50–80 sec). These empirical relations represent the average attenuation properties of a zone covering the surface area of about 11,000, 20,000, 28,000 and 38,000 square km and a depth extent of about 60, 80, 95, 110 km, respectively. With increasing window length, the degree of frequency dependence, n, decreases marginally from 1.17 to 1.02, whereas Q₀ increases significantly from 82 to 122. At lower frequencies up to 6 Hz, Q_c⁻¹ of Kachchh Basin is in agreement with other regions of the world, whereas at higher frequencies from 12 to 24 Hz it is found to be low.     

The preliminary interpretation of deep seismic sounding in western Yunnan

The preliminary interpretation of Project western Yunnan 86–87 is presented here. It shows that there obviously exists lateral velocity heterogeneity from south to north in western Yunnan. The depth of Moho increases from 38 km in the southern end of the profile to 58 km in its northern end. The mean crustal velocity is low in the south, and high in the north, about 6.17–6.45 km/s. The consolidated crust is a 3-layer structure respectively, the upper, middle and lower layer. P ₁ ⁰ is a weak interface the upper crust, P ₂ ⁰ and P ₃ ⁰ are the interfaces of middle-upper crust and middle-lower crust respectively. Another weak interface P ₃ ^0′ can be locally traced in the interior of the lower crust. Interface Pg is 0–6 km deep, interface P ₁ ⁰ 9.2–16.5 km deep, and interfaces P ₂ ⁰ and P ₃ ⁰ respectively 17.0–26.5 km, 25.0–38.0 km deep. The velocity of the upper crust gradually increases from the south to the north, and reaches its maxmium between Nangaozhai and Zhiti, where the velocity of basement plane reaches 6.25–6.35 km/s, then it becomes small northward. The velocity of the middle crust varies little, the middle crust is a low velocity layer with the velocity of 6.30 km/s from Jinhe-Erhai fault to the north. The lower crust is a strong gradient layer. There exists respectively a low velocity layer in the upper mantle between Jinggu and Jingyunqiao, and between Wuliangshan and Lancangjiang fault, the velocity of Pn is only 7.70–7.80 km/s, it is also low to the north of Honghe fault, about 7.80 km/s. Interface P_6/0 can be traced on the top of the upper mantle, its depth is 65 km in the southern end of the profile, and 85 km in the northern end. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,15, 427–440, 1993.     

FUNDAMENTAL FUNCTIONS FOR HORIZONTALLY STRATIFIED EARTH*

The potential distribution and the wave propagation in a horizontally stratified earth is considered and the analogy of the mathematical expression for seismic transfer function, electromagnetic and electric kernel functions, and magnetotelluric input impedance is discussed. Although these specific functions are conveniently treated by a separate expression in each method, it is indicated that the function for seismic and electromagnetic methods is mathematically the same with a change in the physical meaning of the variables from one method to the other. Similarly, the identity of the mathematical expressions of the resistivity kernel function and magnetotelluric input impedance is noticed. In each method a specific geophysical function depends on the thickness and the physical properties of the various layers. Every specific function involves two interdependent fundamental functions, that is P_n and Q_n, or P_n and P*_n, having different physical meaning for different methods. Specific functions are expressible as a ratio P_n/Q_n or P*_n/P_n. Fundamental functions may be reduced to polynomials. The fundamental polynomials Q*_n and P*_n describing the horizontally stratified media are a system of polynomials orthogonal on the unit circle, of first and second order, respectively. The interpretation of geophysical problems corresponds to the identification of the parameters of a system of fundamental orthogonal polynomials. The theorems of orthogonal polynomials are applied to the solution of identification problems. A formula for calculating theoretical curves and direct resistivity interpretation is proposed for the case of arbitrary resistivity of the substratum. The basic equation for synthetic seismograms is reformulated in appendix A. In appendix B a method is indicated for the conversion of the seismic transfer function from arbitrary to perfectly reflective substratum.     

The reservoir-associated earthquakes of April 1983 in western Thailand: Source modeling and implications for induced seismicity

On 22 April 1983, a very large area of Thailand and part of Burma were strongly shaken by a rare earthquake (m _b=5.8,M _s=5.9). The epicenter was located at the Srinagarind reservoir about 190 km northwest of Bangkok, a relatively stable continental region that experienced little previous seismicity. The main shock was preceded by some foreshocks and followed by numerous aftershocks. The largest foreshock ofm _b=5.2 occurred 1 week before the main shock, and the largest aftershock ofm _b=5.3 took place about 3 hours after the main shock. Focal mechanisms of the three largest events in this earthquake sequence have been studied by other seismologists using firts-motion data. However, the solutions for the main shock and the largest aftershock showed significant inconsistency with known surface geology and regional tectonics. We reexamined the mechanisms of these three events by using teleseismicP-andS-waveforms and through careful readings ofP-wave first motions. The directions of theP axes in our study range from NNW-SSE to NNE-SSW, and nodal planes strike in the NW-SE to about E-W in agreement with regional tectonics and surface geology. The main shock mechanism strikes 255°, dips 48°, and slips 63.5°. The fault motions during the main shock and the foreshock are mainly thrust, while the largest aftershock has a large strike-slip component. The seismic moment and the stress drop of the mainshock are determined to be 3.86×10²⁴ dyne-cm and 180 bars, respectively. The occurrence of these thrust events appears to correlate with the unloading of the Srinagarind reservoir. The focal depths of the largest foreshock, the main shock, and the largest aftershock are determined to be 5.4 km, 8 km, and 22.7 km, respectively, from waveform modeling and relative location showing a downward migration of hypocenters of the three largest events during the earthquake sequence. Other characteristics of this reservoir-induced earthquake sequence are also discussed.