Source characteristics of the 1992 Nicaragua tsunami earthquake inferred from teleseismic body waves

We analyzed the broadband body waves of the 1992 Nicaragua earthquake to determine the nature of rupture. The rupture propagation was represented by the distribution of point sources with moment-rate functions at 9 grid points with uniform spacing of 20 km along the fault strike. The moment-rate functions were then parameterized, and the parameters were determined with the least squares method with some constraints. The centroid times of the individual moment-rate functions indicate slow and smooth rupture propagation at a velocity of 1.5 km/s toward NW and 1.0 km/s toward SE. Including a small initial break which precedes the main rupture by about 10 s, we obtained a total source duration of 110 s. The total seismic moment isM _o =3.4×10²⁰ Nm, which is consistent with the value determined from long-period surface waves,M _o =3.7×10²⁰ Nm. The average rise time of dislocation is determined to be 10 s. The major moment release occurred along a fault length of 160 km. With the assumption of a fault widthW=50 km, we obtained the dislocationD=1.3 m. From andD the dislocation velocity isD=D/0.1 m/s, significantly smaller than the typical value for ordinary earthquakes. The stress drop =1.1 MPa is also less than the typical value for subduction zone earthquakes by a factor of 2–3. On the other hand, the apparent stress defined by 2E _s /M _o , where andE _s are respectively the rigidity and the seismic wave energy, is 0.037 MPa, more than an order of magnitude smaller than . The Nicaragua tsunami earthquake is characterized by the following three properties: 1) slow rupture propagation; 2) smooth rupture; 3) slow dislocation motion.     

Determination of the geopotential scale factor from satellite altimetry

Summary The geopotential scale factor R ₀ =GM/W ₀ has been determined on the basis of satellite altimetry as R _{0=(6 363 672·5±0·3) m} and/or the geopotential value on the geoid W ₀ =(62 636 256·5±3) m ² s ^–2 . It has been stated that R ₀ and/or W ₀ is independent of the tidal distortion of surface W=W ₀ due to the zero frequency tide.

---

¶rt;a nmu amumuu u ama amnmuaa R ₀ =GM/W ₀ =(6 363 672,5±0,3) m u/uu aunmuaa a nmuu¶rt;a W ₀ =(62 636 256,5±3) m² s^–2. m, m R ₀ u/uu W ₀ auum m nm amu a a nuu ¶rt;au nmu W=W ₀ .

     

Source parameters of some large earthquakes in the eastern Mediterranean region based on an itertive maximum entropy technique

TeleseismicP-waves of some large earthquakes that occurred in the eastern Mediterranean region have been analysed by using an iterative maximum entropy technique in order to obtain the independent spectral parameters, the long-period spectral level ₀ and the corner frequencyf ₀ of the far-field displacement spectra.Based on these parameters, the seismic source parameters seismic momentM ₀, source dimensionr, fault lengthl, average displacement u, shear stress drop , radiated energyE _s and apparent stressn are calculated for the considered earthquakes by using Brune's and Madariaga's models.The striking feature of the source parameters obtained in this study is the low stress drop value which varies between 5 and 15 bars. If Madariaga's model had been used, higher stress drop values would have been obtained.The low stress drop earthquakes in the eastern Mediterranean region might be interpreted either by the possible presence of low strength material near the source or by the partial stress drop model.     

(ϖμ/ϖT) P of the lower mantle

We estimate (/T) _P of the lower mantle at seismic frequencies using two distinct approaches by combining ambient laboratory measurements on lower mantle minerals with seismic data. In the first approach, an upper bound is estimated for |(/T) _P | by comparing the shear modulus () profile of PREM with laboratory room-temperature data of extrapolated to high pressures. The second approach employs a seismic tomography constraint ( lnV _S / lnV _P ) _P =1.8–2, which directly relates (/T) _P with (K _S /T) _P . An average (K _S /T) _P can be obtained by comparing the well-established room-temperature compression data for lower mantle minerals with theK _S profile of PREM along several possible adiabats. Both (K _S /T) and (/T) depend on silicon content [or (Mg+Fe)/Sil of the model. For various compositions, the two approaches predict rather distinct (/T) _P vs. (K _S /T) _P curves, which intersect at a composition similar to pyrolite with (/T) _P =–0.02 to –0.035 and (K _S /T) _P =–0.015 to –0.020 GPa/K. The pure perovskite model, on the other hand, yields grossly inconsistent results using the two approaches. We conclude that both vertical and lateral variations in seismic velocities are consistent with variation due to pressure, temperature, and phase transformations of a uniform composition. Additional physical properties of a pyrolite lower mantle are further predicted. Lateral temperature variations are predicted to be about 100–250 K, and the ratio of ( lnp/ lnV _S ) _P around 0.13 and 0.26. All of these parameters increase slightly with depth if the ratio of ( lnV _S / lnV _P ) _P remains constant throughout the lower mantle. These predicted values are in excellent agreement with geodynamic analyses, in which the ratios ( ln / lnV _S ) _P and ( / lnV _S ) _P are free parameters arbitrarily adjusted to fit the tomography and geoid data.     

Radioactivity and geochemistry of granulites and some related rocks from the saxonian granulite complex

Summary The distribution of radioactive(Th, U, K), major and selected trace(Rb, Sr, Ba, Y, Zr, V, Cr, Ni) elements of granulites from the Saxonian Granulite Complex was studied. Similarly to the South Bohemian granulites, the Saxonian granulites can be divided according to the contents of their major and trace elements into two main groups, groupA containing mostly acid and subacid granulites (K ₂ O>2.5%, SiO ₂ >68%), and groupB containing mostly intermediate and basic granulites (K ₂ O<2.5%, SiO ₂ <68%). Statistically significant differences between groupsA andB were found for all major oxides and several trace elements(Rb, V, Cr, Ni). The Saxonian granulites follow the same calc-alkaline trend as the South Bohemian, granulitesA being placed mostly in the rhyolite field and granulitesB mostly in the dacite, andesite and basalt fields of this trend. The investigated granulites are characterized by a considerable scatter ofTh andU contents accompanied by very variableTh/U ratios; theTh andU concentrations of granulitesA are substantially lower than is usual for rocks of corresponding acidity.

---

¶rt;a an¶rt;u a¶rt;uamu(Th, U, K) u ua ¶rt;u(Rb, Sr, Ba, Y, Zr, V, Cr, Ni) m aum n¶rt;a aaum na. naa, m u¶rt;aum n uu aam n aaum u ¶rt;u am aua, u u uu. aum n u uu ma a¶rt;um ¶rt; ¶rt;nn; nnA nua¶rt;ama a au¶rt; u au¶rt;aum (K ₂ O>2,5%, Si O ₂ >68%), nnB ¶rt;u u aum (K ₂ O<2,5%, SiO ₂ <68%). ¶rt; muunnau mm mamumuu m au ¶rt; a u u ¶rt; m ¶rt;u m(Rb, V, Cr, Ni). auaum n¶rt;¶rt;m um- m¶rt; a u -uaum;aumA a¶rt;ma a uum n, uaumB a a ¶rt;aum, a¶rt;um u aam n m m¶rt;a. ¶rt;aum — u unnA — aamum uu ¶rt;au da¶rt;uamu mTh uU.

     

Source function of stress waves of a spherical explosive source

Summary The source functions of the stress wave patterns at an elastic source of these waves are analysed. The comparison of the properties of the functions with the stress wave records obtained earlier showed that their parameters do not satisfy, to a greater or lesser extent, the stress wave patterns in the neighbourhood of explosive sources. For this reason a new source function (1) was defined, which fully approximates the observed stress wave patterns in gravel sandy soil. The coefficientsP ₀, , , and were experimentally determined as functions of the distance from the source, its size and the radius of the elastic source in the medium considered. The properties of source function (1) are demonstrated on an example.Paper presented at the XIIIth General Assembly of the European Seismological Commission, Braov (Romania), 28 August to 6 September, 1972.     

Some comparisons between mining-induced and laboratory earthquakes

Although laboratory stick-slip friction experiments have long been regarded as analogs to natural crustal earthquakes, the potential use of laboratory results for understanding the earthquake source mechanism has not been fully exploited because of essential difficulties in relating seismographic data to measurements made in the controlled laboratory environment. Mining-induced earthquakes, however, provide a means of calibrating the seismic data in terms of laboratory results because, in contrast to natural earthquakes, the causative forces as well as the hypocentral conditions are known. A comparison of stick-slip friction events in a large granite sample with mining-induced earthquakes in South Africa and Canada indicates both similarities and differences between the two phenomena. The physics of unstable fault slip appears to be largely the same for both types of events. For example, both laboratory and mining-induced earthquakes have very low seismic efficiencies where _a is the apparent stress and is the average stress acting on the fault plane to cause slip; nearly all of the energy released by faulting is consumed in overcoming friction. In more detail, the mining-induced earthquakes differ from the laboratory events in the behavior of as a function of seismic momentM ₀. Whereas for the laboratory events 0.06 independent ofM ₀, depends quite strongly onM ₀ for each set of induced earthquakes, with 0.06 serving, apparently, as an upper bound. It seems most likely that this observed scaling difference is due to variations in slip distribution over the fault plane. In the laboratory, a stick-slip event entails homogeneous slip over a fault of fixed area. For each set of induced earthquakes, the fault area appears to be approximately fixed but the slip is inhomogeneous due presumably to barriers (zones of no slip) distributed over the fault plane; at constant , larger events correspond to larger _a as a consequence of fewer barriers to slip. If the inequality _a / 0.06 has general validity, then measurements of _a =µE _a /M ₀, where is the modulus of rigidity andE _a is the seismically-radiated energy, can be used to infer the absolute level of deviatoric stress at the hypocenter.     

Comparison of short- and long-periodP-wave travel times and arrival angles

Summary Using about 100 earthquakes recorded on both the short-period and long-period matched three-component seismograph systems at Umeå (UME) in northern Sweden, relativeP-wave travel times (LP-SP) are measured and arrival angles (incidence and azimuth) are determined from relative amplitudes ofP. Long-period onsets appear to register consistently later than short-period ones by a couple of tenths of a second. This figure is of the same order as the uncertainty in T _LS (difference between LP and SP onset times), but still is shown to be of some statistical significance. The effect is believed due, at least in part, to the relative inability of the long-period seismometer to visibly register an initial shorter-period lower-amplitude focal phase from some earthquakes.The dip on crustal interfaces implied by the short-period azimuth and incidence-angle anomalies is 7 to 15 times greater than that implied by the long-period anomalies. This effect is believed due to undulation or irregularity of, foremost, the Moho, as well as other crustal boundaries. The wavelength of, say, Moho undulations would seem to be limited above and below by considerations of the dominant long-period wavelengths (100 km) and of the diameter of Moho sampled by the short-period waves (50 km). The possible undesirable effect of the crustal layering on the components of surface motion, and thus on the arrival angles, is discussed and cautiously dismissed, due to the absence of sedimentary layering, in favour of the above interpretation.     

Triggered earthquakes as stress gauge: Implication for the uppercrust behavior in the Grenoble area,France

A sequence of moderate shallow earthquakes (3.5M _L5.3) was located within the Vercors massif (France) in the period 1961–1984. This subalpine massif has been a low seismic area for at least 5 centuries. During the period 1962–1963, 12 shallow earthquakes occurred in the neighborhood (10 km) of the Monteynard reservoir, 30 km south of the city of Grenoble. The latest fourM _L4.0 earthquakes occurred in 1979–1984 either at larger distance (35 km) or greater depth (10 km) from the reservoir. Two triggering mechanisms are suggested for this sequence: (i) the direct effect of elastic loading through either increased shear stress or strength reducing by increased pore pressure at depth; (ii) the pore pressure diffusion induced by poroelastic stress change due to the reservoir filling.The weekly water levels, local balanced geological cross sections, and focal mechanisms argue for two types of mechanical connection between the earthquake sequence and the filling cycles of the Monteynard reservoir. The seismic sequence started with the 1962–1963 shallow earthquakes that occurred during the first filling of the reservoir and are typical of the direct effect of elastic loading. The 1979 deeper earthquake is located at a 10 km depth below the reservoir. This event occurred 16 years after the initial reservoir impoundment, but one month after the previous 1963 maximum water level was exceeded. Moreover the yearly reservoir level increased gradually in the period 1962–1979 and has decreased since 1980. Accordingly we suggest that the gradual diffusion of water from reservoir to hypocentral depths decreases the strength of the rock matrices through increased pore pressure. The transition between the two types of seismic response is supported by the analysis ofM _L3.5 earthquakes which all occurred in the period 1964–1971, ranging between 10 and 30 km distance from the reservoir. The three other delayed earthquakes of the 1961–1984 seismic sequence (M _L4 during the 1979–1984 period) are all located 35 km away from the reservoir. Based on the seismic activity, the estimates for the hydraulic diffusivities range between 0.2–10 m²/s, except for the first event that occurred 30 km north of the reservoir, the filling just started. The lack ofin situ measurements of crustal hydrological properties in the area, shared by most of the Reservoir-Induced-Seismicity cases, prevents us from obtaining absolute evidence for the triggering processes. These observations and conceptual models attest that previous recurrence times for moderate natural shocks (4.5M _L5.5) estimated within this area using historical data, could be modified by 0.1–1 MPa stress changes. These small changes in deviatoric stress suggest that the upper crust is in this area nearly everywhere at a state of stress near failure. Although the paucity of both number and size of earthquakes in the French subalpine massif shows that aseismic displacements prevail, our study demonstrates that triggered earthquakes are important tools for assessing local seismic risk through mapping fault zones and identifying their possible seismic behavior.     

Rayleigh waves in an inhomogeneous medium

Summary Dispersion in Rayleigh waves is discussed for semi-infinite media with = ₁(1 ± cos s z) and = ₁(1 ± cosh s z), being the rigidity of the medium. A few workers tried with the above Fourier type of model but failed to find the dispersive nature. Because they neglected s due to the complexity of the calculation they arrived at a non dispersive frequency equation. This difficulty is removed in this paper and a dispersive frequency equation is obtained which shows both direct and inverse dispersion. The second model leads to non-convergent solution forz but shows many interesting results which are also discussed.     

The gross features of the lithosphere-asthenosphere system in Europe from seismic surface waves and body waves

Long-period recordings of dispersive Rayleigh waves along numerous station lines, or profiles, in Europe have for the first time permitted a uniform inversion of these observations based on a new method of phase velocity regionalization.Regional dispersion relations obtained by this method have then been subjected to a complete inversion procedure commonly known as the hedgehog method. The results are presented in a map outlining the thickness of the lower lithosphere (lid) and the shear (S) velocities in both the lid and the asthenosphere channel.A comparison of these results with the minimum compressional (P) wave velocities in the asthenosphere and their corresponding depths provides an estimate of theV _p /V _s ratio for the asthenosphere in the European area.Contribution No. 314, Institute of Geophysics, ETH-Zürich, Switzerland.     

Generalized thermoelastic longitudinal waves in unbounded medium

Summary In this paper the generalized thermoelastic longitudinal waves and the temperature field set up due to coupling of the displacement and the temperature fields, with heat wave travelling with certain finite velocity, in an unbounded medium are studied. The thermoelastic displacement potential and the temperature field at any point are obtained in terms of the surface integrals involving the potential, the temperature and their normal derivatives.Notation x _i the cartesian coordinate system,i=1,2,3 - n (u _i) the displacement vector - (/x _i) the del operator - _t /t the derivative with respect to time - T _o the temperature corresponding to the natural stat of zero stress and strain - T Absolute temperature - c _e the specific heat - , Lamé's constants - ₀ the density - coefficient of linear thermal expansion - K thermal conductivity coefficient - _kk u - ₀ the relaxation time     

Propagation of surface waves on a nonhomogenous spherically aeolotropic shell

Summary This paper studies the propagation of Surface Waves on a spherically aeolotropic shell surrounded by vacuum. The elastic constantsc _ij and density of the material of the shell are assumed to be of the form _ij r ^l and _o r ^m respectively, where _{ij o} are constants andl, m are any integers.     

On the possibilities of horizontal propagation ofPc1 pulsations in the ionosphere

Summary The vertical distribution of the contribution of the energy flux density due to the Alfvén(ordinary) wave, guided by the geomagnetic field(and propagating through the ionosphere to the Earth's surface) in the horizontal direction is demonstrated in the mechanism of the horizontal propagation of the Pc1 signal. The distribution with height is shown of the variations of the polarization characteristics of the propagating wave(e.g. the rotation of the polarization plane, changes in ellipticity, attenuation, etc.), which are the result of coupling in the denser layers of the low ionosphere in which also suitable isotropic(extraordinary) modes are generated. The results obtained using the method described in[4, 13] are demonstrated on a model of the daytime ionosphere under incidence of ordinaryL-modes, frequency f=0.3 Hz, and various meridional angles at the ionosphere.

---

auauma anmau uaa Pc1 naa m an¶rt;u ¶rt;u nmmu ma uu uma anauu maum n n¶rt; , anma u nmu. naa m an¶rt;u uu aamumu nuauu anma (nauau nmu nuauu, uu unmumu, amau u m.¶rt;.), m m ¶rt;mu au¶rt;mu na uu u . ¶rt; mum n¶rt;¶rt;u umn() ¶rt;. mam num m¶rt; [4, 13] ¶rt;mua ¶rt;u ¶rt; u nu na¶rt;uu a u L-¶rt; amm f=0,3 n¶rt; au u¶rt;uau au.

     

Seismic signals in geothermal areas of active volcanism: a case study from ‘La Fossa’, Vulcano (Italy)

Since the last eruption (1888–1890) volcanism at Vulcano, Aeolian Archipelago, southern Tyrrhenian Sea, has taken the form of persistent fumarolic activity. The gas-vapour phases of the geothermal systems are mainly discharged within two restricted areas about 1 km apart from each other, in the northern part of the island. These areas are La Fossa crater, and the beach fumaroles of the Baia di Levante. Fluids released at the two main fumarolic fields display quite different chemical and temperature characteristics, implying different origins. The local seismicity essentially takes the form of discrete shocks of shallow origin (depth1 km) at La Fossa, usually with energy < 10¹³ ergs. They are thought to be related to the uprise of pressurized hot gases and vapours discharged at the crater fumaroles. The present investigation points to the existence of two principal categories of seismic events (called M-shocks and N-shocks). These are short events (normally < 10 s). M-type shocks are thought to be due to resonance vibrations within the interior of the volcano, probably driven by the excitation of shock-waves within cavities deeply affected by deposition and alteration of self-sealant hydrothermal minerals. N-type events display features that resemble those of volcano-tectonic earthquakes, but have no recognizable S-phases. Here they are tentatively attributed to microfracturing of rocks which have been extensively hydrothermally altered. Results of the present study permit a preliminary conceptual model of the local shallow seismic processes in the framework of geochemical modelling of fumarolic activity and geological inferences from geothermal drilling.     

Intensity of vlf and elf phenomena in the outer ionosphere

Summary Spatial and temporal variability of natural ELF-VLF phenomena field intensity has been studied using data from Interkosmos satellites with apogees below 2000 km. The results show a distinct dependence of the diurnal variation and latitudinal intensity distribution on the magnetic activity. The L-value at which the maximum of plasmaspheric noise occurs, e.g., correlates better with the Dst-index than with other indicies.

---

mama u a auau anmu n mm - u a ua a uu ¶rt;a uu m an u 2000 . mam naam m auum m ¶rt;a u um an¶rt;u uu m aum amumu. uua L-naama, a m nm au naa uu, um Dst-u¶rt; ¶rt;uu u¶rt;auaum amumu.

     

Accuracy of integral migration transformation in dependence on space-time sampling

Summary The accuracy of wave field extrapolation is studied with respect to the discretization of field data and integral extrapolator. Assuming a far-field approximation of the Rayleigh-Sommerfeld solution for a two-dimensional scalar wave equation, the minimum and the maximum transmitted frequency are expressed as functions of the sampling intervals t, x, and the half-width x₀ and angle _a of the migration aperture. The theoretical limitation of the transmitted frequency band is tested on numerical examples.

---

aamuam mm manuu auumu m ¶rt;umuauu n u uma manu nama. ¶rt; u uma u -¶rt;a ¶rt; ¶rt; a au, ¶rt; nuuuu ¶rt;a n mu ¶rt; uua u aua n¶rt;aa amm a uu m a -nmam ¶rt;umuauu t u , nuu ₀ u a _a uau anm. mu n¶rt;u amm ¶rt;uanaa mmua a u nua.

     

Is there any relation between the sun's motion and global seismic activity?

Summary The seismic energy released by global earthquake activity with time was correlated with F=||, where || is the absolute value of the change of the Sun's acceleration with time. For deep earthquakes and probability P=0.95, the coefficient of correlation was found to be around 0.4. For shallow earthquakes, the dependence was not proved.

---

¶rt;a u uu (1900–1982.) ua ¶rt; F=||,¶rt; || am au uu u a. mu (h>60 ) u ¶rt; mmu =0.95 uum uu ¶rt;muam0.45. nm mu auum a n¶rt;m¶rt;a.

     

Palaeointensity of the geomagnetic field during upper cainozoic derived from palaeo-slags and porcellanites in north Bohemia

Summary Porcellanites and palaeo-slags from North Bohemia are natural materials which can be used to derive the palaeomagnetic directions and palaeointensity of the geomagnetic field active at the time of caustic alteration. The origin of these rocks, called erdbrands, was due to the caustic alteration of predominantly pelitic sediments as a result of underground fires conditioned by spontaneous ignition of coal seams. The caustic alteration occurred during the Upper Pliocene to the Quaternary. Three procedure based on the methods by Thellier and Nagata are presented in the paper. The newly developed apparatus MAVACS (Magnetic Vacuum Control System) was used for the thermal demagnetization of samples. A procedure based on multi-component analysis was also proposed and tested. Besides some methodic results, it was found that the geomagnetic field intensity varied during the respective period within the limits of 48%±4% to 154%±32% of the present geomagnetic field intensity.

---

aum u na au a mumuu uu n¶rt;mam nu¶rt; amua, m n¶rt;¶rt;um ¶rt; ¶rt;u naaum anau u naumumuaum n, ¶rt;m amu uu. mu n¶rt;, aa ¶rt;a¶rt;, uu n¶rt; uu amu uu num num a¶rt; n¶rt; ¶rt;mu aau . amu uu u m nu¶rt; m nua ¶rt; mmu nu¶rt;a. am n¶rt; mu umnmau nua, n¶rt;¶rt;u ¶rt; ¶rt;u naumumu, nuau a m¶rt; u aama. a aamaa annaama MAVACS (Magnetic Vacuum Control System) a unaa ¶rt; mu aauuau ¶rt; amu aa. n¶rt; u n n¶rt;¶rt;, a a munm aau amuauu. nu m¶rt;uu au, ma ma, m umumaum n u¶rt; nu¶rt; a n¶rt;a 48%±4%-154%±32% au umumu aum n.

     

On the possibility of long-range forecasts of geomagnetic activity

Summary The relation between geomagnetic activity and solar activity is studied in long intervals and in individual 11-year cycles, characterized by certain basic parameters, in connection with investigating the prognostic significance of long-term fluctuations of geomagnetic activity [1, 2]. A number of properties as well as sufficiently close relations with a prognostic significance have been found, such as the secular variation of geomagnetic activity, the 22-year cycle of change in the form of the behaviour of geomagnetic activity in 11-year cycles, long-term concurrence of solar and geomagnetic activity and its changes within the 11-year cycle, relations between solar and geomagnetic parameters characterizing the 11-year cycle.

---

a auu uu u ¶rt; uaum amum ¶rt;u nma u u m¶rt; 11-mu ua am mu ¶rt; ¶rt;m nuauaum amumu.