Seismic structure and seismicity of the cooling lava lake of Kilauea Iki, Hawaii

The use of multiple methods is indispensable for the determination of the seismic properties of a complex body such as a partially solidified lava lake. The combined results from various active and passive experiments conducted on Kilauea Iki in March 1976 suggest a rather thin (less than 10 m) zone of residual melt, the lateral extent of which is delineated by a sharp drop in the activity of seismic events occurring within the crust of the lake. The heavily fractured crust is characterized by P velocities ranging from 0.8 to 2.2 km/s and the P velocity could be even lower in the melt possibly because of the presence of gas bubbles. The present seismicity of the lake is compatible with the formation of tensile cracks caused by contraction due to cooling.     

Investigation on the Load/Unload Response Ratio With Coda Q-1 as Response

The load/unload response ratio YQ with the geophysical parameter coda Q-1 of the crust as response is denned in this study.The variation in YQ-1 before and after the Northridge earthquake of January 17,1994(California)has been investigated by using the data of coda Q-1 with frequencies of 1.5,3.0,6.0,12.0,and 24.0 Hz in the Southern California from 1987 to 1994.It can be found that YQ-1 for coda waves with all frequencies,the frequency of 12.0 Hz excluded,ascended to a certain extent prior to the occurrence of the rnainshock and returned to normality after the main shock.     

Seismicity simulation with a rate- and state-dependent friction law

The dynamic motions and stabilities of a single-degree-of-freedom elastic system controlled by different friction laws are compared. The system consists of a sliding block connected to an elastic spring, driven at a constant velocity. The friction laws are a laboratory-inferred friction law called the rate-and-state-dependent friction law, proposed by Dieterich and Ruina, and a simple friction law described by dynamic and static frictions. We further extend the solution to a one-dimensional mass-spring model which is an analog of a fault controlled by the rate-and-state-dependent friction law. This model predicts non uniform slip and stress drop along the rupture length of a heterogeneous fault. This result is very different from some earlier modelings based on the simple friction law and a slip weakening friction law. In those earlier modelings the stress and slip functions become smoother with time along the length of the fault rupture, owing to the interactions between fault segments during slip. Because of this smoothing process the number of small events will decrease with time, and the universilly observed stationary magnitude-frequency relation cannot be explained. The interaction between a fault segment and its neighboring segments can be measured when the post-slip stress on this segment is compared with the stress on an identical segment (represented by a block in this modeling) without neighboring segments. If the post-slip stress of the former is much higher than that of the latter, strong interaction exists; if the two are close, only weak interaction exists. The interaction is determined by the relative motion between fault segments and the time duration of interaction. Our new modeling with the rate-and-state-dependent friction law appears to show no such smoothing effect and provides a physical mechanism for the roughening process in the difference between the fault strength and stress that is necessary to explain the observed stationary magnitude-frequency relation. The noninstantaneous healing predicted by the rate-and-state-dependent friction law may be repsonsible for the recurring nonuniform slip and stress drop, and may be explained by the reduction of interaction among fault segments due to the low frictional strength during the fault stopping. The very low friction during slip stopping allows much longer times than does the higher friction due to instantaneous healing for the fault segments to adjust their motions from an upper-limit slip velocity to almost rest. According to newton's second law, a process with fixed masses and constant velocity changes involves low forces and weak interactions if it is accomplished in a long time period, and vice versa. Our modeling also indicates that the existence of strong patches with higher effective stress on a fault is needed for the occurrence of major events. The creeping section of a fault, such as the one along the San Andreas fault in central California, on the other hand, can be simulated with the rate-and-state-dependent friction law by certain model parameters, which, however, must not include strong patches. In this case small earthquakes and aseismic creep relieve the accumulating strain without any large events.     

The fractal nature of the inhomogeneities in the lithosphere evidenced from seismic wave scattering

In this paper we show evidences of the fractal nature of the 3-D inhomogeneities in the lithosphere from the study of seismic wave scattering and discuss the relation between the fractal dimension of the 3-D inhomogeneities and that of the fault surfaces. Two methods are introduced to measure the inhomogeneity spectrum of a random medium: 1. the coda excitation spectrum method, and 2. the method of measuring the frequency dependence of scattering attenuation. The fractal dimension can be obtained from the inhomogeneity spectrum of the medium. The coda excitation method is applied to the Hindu-Kush data. Based on the observed coda excitation spectra (for frequencies 1–25 Hz) and the past observations on the frequency dependence of scattering attenuation, we infer that the lithospheric inhomogeneities are multiple scaled and can be modeled as a bandlimited fractal random medium (BLFRM) with an outer scale of about 1 km. The fractal dimension of the 3-D inhomogeneities isD ₃=31/2–32/3, which corresponds to a scaling exponent (Hurst number)H=1/2–1/3. The corresponding 1-D inhomogeneity spectra obey the power law with a powerp=2H+1=2–5/3. The intersection between the earth surface and the isostrength surface of the 3-D inhomogeneities will have fractal dimensionD ₁=1.5–1.67. If we consider the earthquake fault surface as developed from the isosurface of the 3-D inhomogeneities and smoothed by the rupture dynamics, the fractal dimension of the fault trace on the surface must be smaller thanD ₁, in agreement with recent measurements of fractal dimension along the San Andreas fault.     

Bias in the estimate of seismic moment tensor by the linear inversion method

Summary. We investigate the effects of various sources of error on the estimation of the seismic moment tensor using a linear least squares inversion on surface wave complex spectra. A series of numerical experiments involving synthetic data subjected to controlled error contamination are used to demonstrate the effects. Random errors are seen to enter additively or multiplicitively into the complex spectra. We show that random additive errors due to background recording noise do not pose difficulties for recovering reliable estimates of the moment tensor. On the other hand, multiplicative errors from a variety of sources, such as focusing, multipathing, or epicentre mislocation, may lead to significant overestimation or underestimation of the tensor elements and in general cause the estimates to be less reliable.     

Observations experimentales de la constance de la pesanteur au Bureau International des Poids et Mesures

Sans résumé Présenté à la Commission Gravimétrique Internationale, 7–11 Septembre 1970, Paris.