Search

Dense strong motion observation networks provided us with valuable data for studying strong motion generation from large earthquakes. From kinematic waveform inversion of seismic data, the slip distribution on the fault surface of large earthquakes is known to be spatially heterogeneous. Because heterogeneities in the slip and stress drop distributions control the generation of near-source ground motion, it is important to characterize these heterogeneities for past earthquakes in constructing a source model for reliable prediction of strong ground motion. The stress changes during large earthquakes on the faults recently occurring in Japan are estimated from the detailed slip models obtained by the kinematic waveform inversion. The stress drops on and off asperities are summarized on the basis of the stress change distributions obtained here. In this paper, we define the asperity to be a rectangular area whose slip is 1.5 or more times larger than the average slip over the fault according to the previous study for inland crustal earthquakes. The average static stress drops on the asperities of the earthquakes studied here are in the range 6-23 MPa, whereas those off the asperities are below 3 MPa. We compiled the stress drop on the asperities together with a data set from previous studies of other inland earthquakes in Japan and elsewhere. The static stress drop on the asperity depends on its depth, and we obtained an empirical relationship between the static stress drop and the asperity's depth. Moreover, surface-breaking asperities seemed to have smaller stress drops than buried asperities. Simple ground motion simulations using the characterized asperity source models reveal that deep asperities generate larger ground motion than shallow asperities, because of the different stress drops of the asperities. These characteristics can be used for advanced source modeling in strong ground motion prediction for inland crustal earthquakes. [ABSTRACT FROM AUTHOR]

Izmir, the third largest city and one of the major economic centers in Turkey, has more than three million residents and one-half million buildings. The city, located in a seismically active region in western Anatolia, was a subject of the 1997 RADIUS (Risk Assessment Tools for Diagnosis of Urban Areas against Seismic Disaster) project. In this paper, the seismic hazard of Izmir is investigated through probabilistic seismic hazard assessment. First, the seismic setting of Izmir is presented. Considering the statistics of earthquakes that took place in the region during the period 1900-2005, a simple seismic hazard model is used to facilitate the assessment. To account for modeling uncertainties associated with the values of seismicity parameters, a logic tree procedure is employed in carrying out the seismic hazard computations. The resulting weighted average seismic hazard, presented in terms of peak ground acceleration and associated probability of exceedence, could be considered the 'best estimate' of seismic hazard for Izmir. Accordingly, for a return period of 475 years, for rock sites, a PGA value of 0.34 g is calculated. This PGA hazard estimate is close to the current code-recommended design acceleration level for Izmir. [ABSTRACT FROM AUTHOR]

Seismic noise is a fundamental part of seismic data which cannot be avoided when conducting any seismic survey. It consists of coherent and random noise. Noise removal or filtering is one of the major concerns in the field of seismic processing. In this paper, we introduce an image filtering technique based on a detection-estimation algorithm for Gaussian and random noise removal in seismic data, namely the trilateral filter, based on a statistic called rank-ordered absolute differences. The non-linear and adaptive behaviour of this filter makes it very robust in the presence of random and coherent noise, in addition to its computational simplicity and its ability to automatically identify noise in data. We have modified the strategy of trilateral filtering by adapting the rank-ordered absolute differences formula in order to extract the signal component. We have successfully used this filter for the removal of surface waves and random spiky noise from synthetic and field data. Results are very encouraging and show the superiority of this filter compared with other filters, particularly when used recursively. [ABSTRACT FROM AUTHOR]

In this paper, we present a model for studying aftershock sequences that integrates Coulomb static stress change analysis, seismicity equations based on rate-state friction nucleation of earthquakes, slip of geometrically complex faults, and fractal-like, spatially heterogeneous models of crustal stress. In addition to modeling instantaneous aftershock seismicity rate patterns with initial clustering on the Coulomb stress increase areas and an approximately 1/t diffusion back to the pre-mainshock background seismicity, the simulations capture previously unmodeled effects. These include production of a significant number of aftershocks in the traditional Coulomb stress shadow zones and temporal changes in aftershock focal mechanism statistics. The occurrence of aftershock stress shadow zones arises from two sources. The first source is spatially heterogeneous initial crustal stress, and the second is slip on geometrically rough faults, which produces localized positive Coulomb stress changes within the traditional stress shadow zones. Temporal changes in simulated aftershock focal mechanisms result in inferred stress rotations that greatly exceed the true stress rotations due to the main shock, even for a moderately strong crust (mean stress 50 MPa) when stress is spatially heterogeneous. This arises from biased sampling of the crustal stress by the synthetic aftershocks due to the non-linear dependence of seismicity rates on stress changes. The model indicates that one cannot use focal mechanism inversion rotations to conclusively demonstrate low crustal strength (≤10 MPa); therefore, studies of crustal strength following a stress perturbation may significantly underestimate the mean crustal stress state for regions with spatially heterogeneous stress. [ABSTRACT FROM AUTHOR]

The ability to successfully predict the future behavior of a system is a strong indication that the system is well understood. Certainly many details of the earthquake system remain obscure, but several hypotheses related to earthquake occurrence and seismic hazard have been proffered, and predicting earthquake behavior is a worthy goal and demanded by society. Along these lines, one of the primary objectives of the Regional Earthquake Likelihood Models (RELM) working group was to formalize earthquake occurrence hypotheses in the form of prospective earthquake rate forecasts in California. RELM members, working in small research groups, developed more than a dozen 5-year forecasts; they also outlined a performance evaluation method and provided a conceptual description of a Testing Center in which to perform predictability experiments. Subsequently, researchers working within the Collaboratory for the Study of Earthquake Predictability (CSEP) have begun implementing Testing Centers in different locations worldwide, and the RELM predictability experiment—a truly prospective earthquake prediction effort—is underway within the U.S. branch of CSEP. The experiment, designed to compare time-invariant 5-year earthquake rate forecasts, is now approximately halfway to its completion. In this paper, we describe the models under evaluation and present, for the first time, preliminary results of this unique experiment. While these results are preliminary—the forecasts were meant for an application of 5 years—we find interesting results: most of the models are consistent with the observation and one model forecasts the distribution of earthquakes best. We discuss the observed sample of target earthquakes in the context of historical seismicity within the testing region, highlight potential pitfalls of the current tests, and suggest plans for future revisions to experiments such as this one. [ABSTRACT FROM AUTHOR]

The objective of this paper is to quantify the use of past seismicity to forecast the locations of future large earthquakes and introduce optimization methods for the model parameters. To achieve this the binary forecast approach is used where the surface of the Earth is divided into l° × l° cells. The cumulative Benioff strain of m ≥ m c earthquakes that occurred during the training period, Δ Ttr, is used to retrospectively forecast the locations of large target earthquakes with magnitudes ≥ m T during the forecast period, Δ Tfor. The success of a forecast is measured in terms of hit rates (fraction of earthquakes forecast) and false alarm rates (fraction of alarms that do not forecast earthquakes). This binary forecast approach is quantified using a receiver operating characteristic diagram and an error diagram. An optimal forecast can be obtained by taking the maximum value of Pierce’s skill score. [ABSTRACT FROM AUTHOR]

A brief account is given of the principles of stochastic modelling in seismology, with special regard to the role and development of stochastic models for seismicity. Stochastic models are seen as arising in a hierarchy of roles in seismology, as in other scientific disciplines. At their simplest, they provide a convenient descriptive tool for summarizing data patterns; in engineering and other applications, they provide a practical way of bridging the gap between the detailed modelling of a complex system, and the need to fit models to limited data; at the most fundamental level they arise as a basic component in the modelling of earthquake phenomena, analogous to that of stochastic models in statistical mechanics or turbulence theory. As an emerging subdiscipline, statistical seismology includes elements of all of these. The scope for the development of stochastic models depends crucially on the quantity and quality of the available data. The availability of extensive, high-quality catalogues and other relevant data lies behind the recent explosion of interest in statistical seismology. At just such a stage, it seems important to review the underlying principles on which statistical modelling is based, and that is the main purpose of the present paper. [ABSTRACT FROM AUTHOR]

We develop a constitutive model for rocks that are constituted from brittle particles, based on the theory of breakage mechanics. The model connects between the energetics and the micromechanics that drive the process of confined comminution. Given this ability, our model not only describes the entire stress-strain response of the material, but also connects this response to predicting the evolution of the grain size distribution. The latter fact enables us to quantify how the permeability reduces within cataclasite zones, in relation to aspects of grain crushing. Finally, our paper focuses on setting a framework for quantifying how the energy budget of earthquakes is expensed in relation to dissipation events in cataclasis. We specifically distinguish between the dissipation directly from the creation of new surface area, which causes further breakage dissipation from the redistribution of locked-in stored energy from surrounding particles, dissipations from friction and from the configurational reorganisation of particles. [ABSTRACT FROM AUTHOR]

In this paper, ground motion during the Independence Day earthquake of August 15, 1950 (Mw 8.6, B en-M enahem et al., 1974) in the northeastern part of India is estimated by seismological approaches. A hybrid simulation technique which combines the low frequency ground motion simulated from an analytical source mechanism model with the stochastically simulated high-frequency components is used for obtaining the acceleration time histories. A series of ground motion simulations are carried out to estimate the peak ground acceleration (PGA) and spectral accelerations at important cities and towns in the epicentral region. One sample PGA distribution in the epicentral region encompassing the epicenter is also obtained. It is found that PGA in the epicentral region has exceeded 1 g during this earthquake. The estimated PGA’s are validated to the extent possible using the MMI values. The simulated acceleration time histories can be used for the assessment of important engineering structures in northeastern India. [ABSTRACT FROM AUTHOR]

There are two fundamentally different approaches to assessing the probabilistic risk of earthquake occurrence. The first is fault based. The statistical occurrence of earthquakes is determined for mapped faults. The applicable models are renewal models in that a tectonic loading of faults is included. The second approach is seismicity based. The risk of future earthquakes is based on the past seismicity in the region. These are also known as cluster models. An example of a cluster model is the epidemic type aftershock sequence (ETAS) model. In this paper we discuss an alternative branching aftershock sequence (BASS) model. In the BASS model an initial, or seed, earthquake is specified. The subsequent earthquakes are obtained from statistical distributions of magnitude, time, and location. The magnitude scaling is based on a combination of the Gutenberg-Richter scaling relation and the modified Båth’s law for the scaling relation of aftershock magnitudes relative to the magnitude of the main earthquake. Omori’s law specifies the distribution of earthquake times, and a modified form of Omori’s law specifies the distribution of earthquake locations. Unlike the ETAS model, the BASS model is fully self-similar, and is not sensitive to the low magnitude cutoff. [ABSTRACT FROM AUTHOR]

The purpose of this paper is to discuss the statistical distributions of recurrence times of earthquakes. Recurrence times are the time intervals between successive earthquakes at a specified location on a specified fault. Although a number of statistical distributions have been proposed for recurrence times, we argue in favor of the Weibull distribution. The Weibull distribution is the only distribution that has a scale-invariant hazard function. We consider three sets of characteristic earthquakes on the San Andreas fault: (1) The Parkfield earthquakes, (2) the sequence of earthquakes identified by paleoseismic studies at the Wrightwood site, and (3) an example of a sequence of micro-repeating earthquakes at a site near San Juan Bautista. In each case we make a comparison with the applicable Weibull distribution. The number of earthquakes in each of these sequences is too small to make definitive conclusions. To overcome this difficulty we consider a sequence of earthquakes obtained from a one million year “Virtual California” simulation of San Andreas earthquakes. Very good agreement with a Weibull distribution is found. We also obtain recurrence statistics for two other model studies. The first is a modified forest-fire model and the second is a slider-block model. In both cases good agreements with Weibull distributions are obtained. Our conclusion is that the Weibull distribution is the preferred distribution for estimating the risk of future earthquakes on the San Andreas fault and elsewhere. [ABSTRACT FROM AUTHOR]

The LURR theory is a new approach for earthquake prediction, which achieves good results in earthquake prediction within the China mainland and regions in America, Japan and Australia. However, the expansion of the prediction region leads to the refinement of its longitude and latitude, and the increase of the time period. This requires increasingly more computations, and the volume of data reaches the order of GB, which will be very difficult for a single CPU. In this paper, a new method was introduced to solve this problem. Adopting the technology of domain decomposition and parallelizing using MPI, we developed a new parallel tempo-spatial scanning program. [ABSTRACT FROM AUTHOR]

In this paper the statistical properties of solar flares and earthquakes are compared by analyzing the energy distributions, the time series of energies and interevent times, and, above all, the distributions of interevent times per se. It is shown that the two phenomena have different statistics of scaling, and even the same phenomenon, when observed in different periods or at different locations, is characterized by different statistics that cannot be uniformly rescaled onto a single, universal curve. The results indicate apparent complexity of impulsive energy release processes, which neither follow a common behaviour nor could be attributed to a universal physical mechanism. [ABSTRACT FROM AUTHOR]

The December 26, 2004 Sumatra-Andaman earthquake that registered a moment magnitude ( M w ) of 9.1 was one of the largest earthquakes in the world since 1900. The devastating tsunami that resulted from this earthquake caused more casualties than any previously reported tsunami. The number of fatalities and missing persons in the most seriously affected countries were Indonesia - 167,736, Sri Lanka - 35,322, India - 18,045 and Thailand - 8,212. This paper describes two field visits to assess tsunami effects in Sri Lanka by a combined team of Japanese and Sri Lankan researchers. The first field visit from December 30, 2004 – January 04, 2005 covered the western and southern coasts of Sri Lanka including the cities of Moratuwa, Beruwala, Bentota, Pereliya, Hikkaduwa, Galle, Talpe, Matara, Tangalla and Hambantota. The objectives of the first field visit were to investigate the damage caused by the tsunami and to obtain eyewitness information about wave arrival times. The second field visit from March 10–18, 2005 covered the eastern and southern coasts of Sri Lanka and included Trincomalee, Batticaloa, Arugam Bay, Yala National Park and Kirinda. The objectives of the second visit were mainly to obtain eyewitness information about wave arrival times and inundation data, and to take relevant measurements using GPS instruments. [ABSTRACT FROM AUTHOR]

We develop stochastic approaches to determine the potential for tsunami generation from earthquakes by combining two interrelated time series, one for the earthquake events, and another for the tsunami events. Conditional probabilities for the occurrence of tsunamis as a function of time are calculated by assuming that the inter-arrival times of the past events are lognormally distributed and by taking into account the time of occurrence of the last event in the time series. An alternative approach is based on the total probabilitiy theorem. Then, the probability for the tsunami occurrence equals the product of the ratio, r (= tsunami generating earthquakes/total number of earthquakes) by the conditional probability for the occurrence of the next earthquake in the zone. The probabilities obtained by the total probability theorem are bounded upwards by the ratio r and, therefore, they are not comparable with the conditional probabilities. The two methods were successfully tested in three characteristic seismic zones of the Pacific Ocean: South America, Kuril-Kamchatka and Japan. For time intervals of about 20 years and over the probabilities exceed 0.50 in the three zones. It has been found that the results depend on the approach applied. In fact, the conditional probabilities of tsunami occurrence in Japan are slightly higher than in the South America region and in Kuril-Kamchatka they are clearly lower than in South America. Probabilities calculated by the total probability theorem are systematically higher in South America than in Japan while in Kuril-Kamchatka they are significantly lower than in Japan. The stochastic techniques tested in this paper are promising for the tsunami potential assessment in other tsunamigenic regions of the world. [ABSTRACT FROM AUTHOR]

The concept of state vector stems from statistical physics, where it is usually used to describe activity patterns of a physical field in its manner of coarsegrain. In this paper, we propose an approach by which the state vector was applied to describe quantitatively the damage evolution of the brittle heterogeneous systems, and some interesting results are presented, i.e., prior to the macro-fracture of rock specimens and occurrence of a strong earthquake, evolutions of the four relevant scalars time series derived from the state vectors changed anomalously. As retrospective studies, some prominent large earthquakes occurred in the Chinese Mainland (e.g., the M 7.4 Haicheng earthquake on February 4, 1975, and the M 7.8 Tangshan earthquake on July 28, 1976, etc) were investigated. Results show considerable promise that the time-dependent state vectors could serve as a kind of precursor to predict earthquakes. [ABSTRACT FROM AUTHOR]

21 earthquakes recorded by a temporary seismic network in the Changbaishan Tianchi volcanic area in Northeast China operated during the summer of 2002 and 2003 were analyzed to estimate the S coda attenuation. The attenuation quality factor Q c was estimated using the single scattering attenuation model of Sato (1977) in the frequency band from 4 to 24 Hz. All the events studied in this paper occurred at depths from 2 to 6 km with ML of 1.4–2.8. The epicentral distances are less than 25 km. For all events which occurred near the Tianchi Lake (caldera), the Q c patterns obtained at the stations near the lake are similar, and the Q c values are relatively small. At the stations located about 15 km east of the Tianchi Lake, however, the average Q c is significantly higher. For an event which occurred 25km from the lake to the west, Q c patterns derived at the stations near the lake are quite similar to the above mentioned Q c for stations located in the east. Further study shows that Q c value in the north and central areas of the volcano is relatively lower than that in the surrounding area. Compared to other volcanic areas in the world, the average Q c of the Changbaishan Tianchi volcanic area is obviously lower. The deep seismic sounding and teleseismic receiver function studies indicated more than one lower velocity layer in the crust. The MT studies suggested the presence of high conductive bodies beneath the area. We interpret the strong attenuation of coda waves near the Changbaishan Tianchi volcano as being possibly related to high temperature medium caused by shallow magma chambers. [ABSTRACT FROM AUTHOR]

Time delays associated with processes leading to a failure or stress relaxation in materials and earthquakes are studied in terms of continuum damage mechanics. Damage mechanics is a quasi-empirical approach that describes inelastic irreversible phenomena in the deformation of solids. When a rock sample is loaded, there is generally a time delay before the rock fails. This period is characterized by the occurrence and coalescence of microcracks which radiate acoustic signals of broad amplitudes. These acoustic emission events have been shown to exhibit power-law scaling as they increase in intensity prior to a rupture. In case of seismogenic processes in the Earth's brittle crust, all earthquakes are followed by an aftershock sequence. A universal feature of aftershocks is that their rate decays in time according to the modified Omori's law, a power-law decay. In this paper a model of continuum damage mechanics in which damage (microcracking) starts to develop when the applied stress exceeds a prescribed yield stress (a material parameter) is introduced to explain both laboratory experiments and systematic temporal variations in seismicity. [ABSTRACT FROM AUTHOR]

Aki (1969) first modeled coda waves of a local earthquake as a superposition of scattered surface waves. This paper attempts to clarify the constituents of surface-wave coda at long periods at very long lapse times. For a large earthquake of magnitude 7 or larger, vertical component oscillation in periods from 90 to 180 s persists for more than 20 hours from the earthquake origin time. Although the early portion of the coda envelope is successfully modeled by assuming incoherent scattered Rayleigh waves by heterogeneities distributed all over the Earth, the later potion of the observed coda envelope (roughly later than 35,000 s) has systematically larger amplitude than theoretical prediction. To clarify the cause of this discrepancy, we studied the constituents of vertical-component seismograms of three large earthquakes recorded by the F-net in Japan using the f- k power spectral analysis. We found that the direct and scattered fundamental-mode Rayleigh waves of velocity about 3.7 km/s are dominant in the earlier part of each envelope. It justifies the use of a scattering model of the fundamental Rayleigh waves for synthesizing the envelope. At lapse times later than 20,000 s–35,000 s, higher modes with phase velocities around 20 km/s become dominant. The transition time to the dominance of higher modes is found to become earlier for a deeper focus earthquake. The small coda attenuation factor from (1.90±0.23) × 10−3 to (2.38±0.32) × 10−3 estimated from later coda envelopes recorded at IRIS stations distributed worldwide also agrees with the attenuation factor of spheroidal modes according to PREM. We may interpret that higher mode waves are uniformly distributed at large lapse time due to large velocity dispersion and/or scattering and they dominate over the fundamental mode waves because of smaller attenuation in the lower mantle. The coda attenuation measurement proposed by Aki is found to be useful even for long periods and at very large lapse times. [ABSTRACT FROM AUTHOR]

The article presents a list of published papers authored or co-authored by geophysicist Keiiti Aki. Periodical articles the list include "Quantitative prediction of earthquake occurrence as stochastic phenomena," "Thermal processes near the earth's crust," "Relay computer for correlogram analysis," and "Correlogram analysis of seismograms."

Seismic hazard in terms of spectral acceleration ( SA) has been estimated for the first time in northern Algeria. For this purpose, we have used the spatially-smoothed seismicity approach. The present paper is intended to be a continuation of previous work in which we have evaluated the seismic hazard in terms of peak ground acceleration ( PGA) using the same methodology. To perform these evaluations, four complete and Poissonian seismic models have been used. One of them considers earthquakes with magnitudes above M S 6.5 in the last 300 years, that is, the most energetic seismicity in the region. Firstly, seismic hazard maps in terms of SA, at periods of 0.1, 0.2, 0.3, 0.4, 0.5, 1.0, 1.5 and 2.0 sec, with 39.3% and 10% probability of exceedance in 50 years, have been obtained. Therefore, uniform hazard spectra ( UHS) are computed and examined in detail for twelve of the most industrial and populated cities in northern Algeria. All the reported results in this study are for rock soil and 5% of damping. It is noteworthy that, in the seismic hazard maps as well as in the UHS plots, we observe maximum SA values in the central area of the Tell. The higher values are reached in the Chleff region (previously El Asnam), specifically around the location of the destructive earthquakes of September 9, 1954 ( M S 6.8), and October 10, 1980 ( M S 7.3). These maximum values, 0.4 g and 1.0 g, are associated with periods of about 0.2 and 0.3 sec for return periods of 100 and 475 years, respectively. [ABSTRACT FROM AUTHOR]

The occurrence of the Algiers earthquake ( M 6.8) of May 21, 2003, has motivated the necessity to reassess the probabilistic seismic hazard of northern Algeria. The fact that this destructive earthquake took place in an area where there was no evidence of previous significant earthquakes, neither instrumental nor historical, strongly encourages us to review the seismic hazard map of this region. Recently, the probabilistic seismic hazard of northern Algeria was computed using the spatially smoothed seismicity methodology. The catalog used in the previous computation was updated for this review, and not only includes information until June 2003, but also considers a recent re-evaluation of several historical earthquakes. In this paper, the same methodology and seismicity models are utilized in an effort to compare this methodology against an improved and updated seismic catalog. The largest mean peak ground acceleration (PGA) values are obtained in northernmost Algeria, specifically in the central area of the Tell Atlas. These values are of the order of 0.48 g for a return period of 475 years. In the City of Algiers, the capital of Algeria, and approximately 50 km from the reported epicenter of this latest destructive earthquake, a new mean PGA value of 0.23 g is obtained for the same return period. This value is 0.07 g greater than that obtained in the previous computation. In general, we receive greater seismic hazard results in the surrounding area of Algiers, especially to the southwest. The main reason is not this recent earthquake by itself, but the significant increase in the mmax magnitude in the seismic source where the city and the epicenter are included. [ABSTRACT FROM AUTHOR]

The statistical properties of aftershock sequences are associated with three empirical scaling relations: (1) Gutenberg-Richter frequency-magnitude scaling, (2) Båth’s law for the magnitude of the largest aftershock, and (3) the modified Omori’s law for the temporal decay of aftershocks. In this paper these three laws are combined to give a relation for the aftershock decay rate that depends on only a few parameters. This result is used to study the temporal properties of aftershock sequences of several large California earthquakes. A review of different mechanisms and models of aftershocks are also given. The scale invariance of the process of stress transfer caused by a main shock and the heterogeneous medium in which aftershocks occur are responsible for the occurrence of scaling laws. We suggest that the observed partitioning of energy could play a crucial role in explaining the physical origin of Båth’s law. We also study the stress relaxation process in a simple model of damage mechanics and find that the rate of energy release in this model is identical to the rate of aftershock occurrence described by the modified Omori’s law. [ABSTRACT FROM AUTHOR]

The development of fault interaction models has triggered the need for an accurate estimation of seismicity rate changes following the occurrence of an earthquake. Several statistical methods have been developed in the past to serve this purpose, each relying on different assumptions (e.g., stationarity, gaussianity) pertaining to the seismicity process.In this paper we review these various approaches, discuss their limitations, and propose further improvements. The feasibility of mapping robust seismicity rate changes, and more particularly rate decreases (i.e., seismicity shadows), in the first few days of an aftershock sequence, is examined. To this aim, the hypothesis of large numbers of earthquakes, hence the use of Gaussian statistics, as is usually assumed, must be dropped.Finally, we analyse the modulation in seismicity rates following the 1992, June 28M w 7.3 Landers earthquake in the region of the 1992, April 22M w 6.1 Joshua Tree earthquake. Clear instances of early triggering (i.e., in the first few days) followed by a seismicity quiescence, are observed. This could indicate the existence of two distinct interaction regimes, a first one caused by the destabilisation of active faults by the travelling seismic waves, and a second one due to the remaining static stress perturbation. [ABSTRACT FROM AUTHOR]

This paper presents an analysis of the distribution of the time t between two consecutive events in a stationary point process. The study is motivated by the discovery of unified scaling laws for t for the case of seismic events. We demonstrate that these laws cannot exist simultaneously in a seismogenic area. Under very natural assumptions we show that if, after rescaling to ensureEt =1, the interevent time has a universal distributionF, thenFmust be exponential. In other words, Corral’s unified scaling law cannot exist in the whole range of time. In the framework of a general cluster model we discuss the parameterization of an empirical unified law and the physical meaning of the parameters involved. [ABSTRACT FROM AUTHOR]

This paper introduces a statistical technique, based on the recently developed Multiscale Trend Analysis (MTA), for quantifying correlations between non-stationary processes observed at irregular non-coincident time grids. We apply this technique to studying the temporal correlation between the dynamics of the ductile and brittle layers in the lithosphere. Our results confirm the previously reported strong positive correlation between the codaQ-1 and seismicity and its drop before major earthquakes observed in California. The proposed technique has significant advantages over the conventional correlation analysis: (1) MTA allows one to work directly with non-coincident time series without preliminary resampling the data; (2) the correlation is defined via the stable objects-trends-rather than noisy individual observations, hence it is highly robust; (3) the correlations are quantified at different time scales. The suggested technique seems promising for the wide range of applied problems dealing with coupled time series. [ABSTRACT FROM AUTHOR]

The recent earthquakes in California and Japan have shown the fundamental role that the road infrastructure plays in emergencies. In fact, only the maintenance of a sufficient level of efficiency can help to quickly reach the affected areas and thus avoid further serious consequences. The necessity of guaranteeing the functionality of the transport network during seismic events therefore requires seismic risk planning extended also to the road infrastructures in order to support the management of post-earthquake emergencies. Analogously it is fundamentally important to have analysis instruments of the road system able to preventatively evaluate the effects of earthquakes in order to identify possible emergencies, therefore preparing a program of intervention to reduce seismic risk on road networks. This paper proposes a methodology for the evaluation of seismic risk of road infrastructures according to the following points:•Study of seismic hazard of the site for the definition of a seismic scenario using attenuation models in relation to historical seismology and the geological and tectonic characteristics of the territory;•Analysis of the direct exposure connected to the probability of the presence of road users on the different parts of the network directly exposed to the seismic event;•Analysis of the indirect exposure relative to the distribution of the population and the infrastructures for which post-earthquake accessibility must be guaranteed;•Evaluation of the functional vulnerability in relation to the potential replaceability of damaged stretches considering network configuration and geometrical characteristics;•Evaluation of structural vulnerability of the stretch correlated to the characteristics (structural, mechanical, technological, etc.) of the different components (bridges, embankments, trenches, tunnels) that make up the stretches obtained by the use of correctly elaborated tables for each component.The determination of global risk indexes of the single stretches and of the network, evaluated by means of a relationship between the ascertained parameters derived from the investigation of the previous points, provides the necessary information for the definition of mitigation measures to reduce the risk and for management planning before and after disaster. The proposed methodology, which has already been applied to a restricted area, is currently being applied to the province of Catania (Sicily, Italy), which is one of the geographical regions of highest seismic risk in Europe, and its future extension to all of eastern Sicily is foreseeable. [ABSTRACT FROM AUTHOR]

- In this paper, the site characteristics of the Dahan downhole array are studied by analyzing the September 21, 1999M7.3 Chi-Chi earthquake sequence including the main shock and some aftershocks. The four-level array (0 m, 50 m, 100 m and 200 m) is located to the north of Hualien City in eastern Taiwan. Polarization analysis is used to check the orientation errors of the seismometers at different levels of depth. If the surface instrument is chosen as reference, the angle between the major polarization axes of the surface and any downhole records is the orientation error that must be corrected for the downhole accelerographs. The orientation errors at depths of 50 m, 100 m and 200 m are 32°, 120° and -84°. After the corrections, the coherency between the surface and downhole records is substantially improved. Spectral ratio analysis shows that the predominant frequency of the Chi-Chi main shock shifts to a lower frequency. We also simulate ground motions at different depths by using the Haskell method with a linear velocity structure model. The record at surface is chosen as the input motion. Compared with the observed data, ground acceleration can be well reproduced for the aftershocks (weak-motion events) of the September 21, 1999M7.3 Chi-Chi earthquake. However, for the Chi-Chi main shock, the synthetic waveform cannot match well with the observation neither in amplitude nor in phase. This indicates that large ground shaking probably induced the nonlinear site effect at that time, and the model used cannot support it. [ABSTRACT FROM AUTHOR]

- An algorithm recently developed by RUNDLEet al. (2002) to find regions of anomalous seismic activity associated with large earthquakes identified the location of anM w = 5.6 earthquake near Calexico, Mexico. In this paper we analyze the regional seismicity before this event, and a nearbyM w = 5.7 event, using time-to-failure algorithms developed by BOWMANet al. (1998) and BOWMAN and KING (2001a,b). The former finds the radius of a circular region surrounding the epicenter that optimizes the time-to-failure acceleration of seismic release. The latter optimizes acceleration based on the expected stress accumulation pattern for a dislocation source. Both methods found a period of accelerating seismicity in an optimal region, the size of which agrees with previously proposed scaling relations. This positive result suggests that the Rundle algorithm may provide a useful technique to identify regions of accelerating seismicity, which can then be analyzed using signal optimization time-to-failure techniques. [ABSTRACT FROM AUTHOR]

This paper provides a generic equation for the evaluation of the maximum earthquake magnitude mmax for a given seismogenic zone or entire region. The equation is capable of generating solutions in different forms, depending on the assumptions of the statistical distribution model and/or the available information regarding past seismicity. It includes the cases (i) when earthquake magnitudes are distributed according to the doubly-truncated Gutenberg-Richter relation, (ii) when the empirical magnitude distribution deviates moderately from the Gutenberg-Richter relation, and (iii) when no specific type of magnitude distribution is assumed. Both synthetic, Monte-Carlo simulated seismic event catalogues, and actual data from Southern California, are used to demonstrate the procedures given for the evaluation of mmax. The three estimates of mmax for Southern California, obtained by the three procedures mentioned above, are respectively: 8.32 ± 0.43, 8.31 ± 0.42 and 8.34 ± 0.45. All three estimates are nearly identical, although higher than the value 7.99 obtained by Field et al. (1999). In general, since the third procedure is non-parametric and does not require specification of the functional form of the magnitude distribution, its estimate of the maximum earthquake magnitude mmax is considered more reliable than the other two which are based on the Gutenberg-Richter relation. [ABSTRACT FROM AUTHOR]

A simplified technique for simulation of wide-band strong motion based on simple regression relations and Empirical Green’s Function (EGF) technique by Irikura (1986) is presented in this paper. The method uses the acceleration envelope as a shaping window for a filtered white Gaussian noise, to get the synthetic accelerogram from each subfault. Correction factors for slip of large and small events and transmission factors at each boundary of different layers are included in this synthetic accelerogram. The synthetic accelerogram obtained from each subfault is used as the Green’s function to get resultant records. Simulations are made for the confirmed models of the Uttarkashi and the Chamoli earthquakes at a number of stations to get wide-band strong ground motion. The comparison of synthetic with the observed records over a wide range of frequencies for two different Himalayan earthquakes establishes the efficacy of the present technique. [ABSTRACT FROM AUTHOR]

This paper provides a generic equation for the evaluation of the maximum earthquake magnitude mmax for a given seismogenic zone or entire region. The equation is capable of generating solutions in different forms, depending on the assumptions of the statistical distribution model and/or the available information regarding past seismicity. It includes the cases (i) when earthquake magnitudes are distributed according to the doubly-truncated Gutenberg-Richter relation, (ii) when the empirical magnitude distribution deviates moderately from the Gutenberg-Richter relation, and (iii) when no specific type of magnitude distribution is assumed. Both synthetic, Monte-Carlo simulated seismic event catalogues, and actual data from Southern California, are used to demonstrate the procedures given for the evaluation of mmax. The three estimates of mmax for Southern California, obtained by the three procedures mentioned above, are respectively: 8.32 ± 0.43, 8.31 ± 0.42 and 8.34 ± 0.45. All three estimates are nearly identical, although higher than the value 7.99 obtained by Field et al. (1999). In general, since the third procedure is non-parametric and does not require specification of the functional form of the magnitude distribution, its estimate of the maximum earthquake magnitude mmax is considered more reliable than the other two which are based on the Gutenberg-Richter relation. [ABSTRACT FROM AUTHOR]

— The paper discusses the very recent seismotectonic and geological investigations for the Sofia valley. It also reviews available historical and recent seismological data. One possible variant of seismic zoning (interpretation of available geological and seismological data) is used to evaluate uniform probability Fourier spectra of strong ground motion at Sofia City. The distribution of the maximum expected displacements, using 1-D models and SHAKE’91 techniques for 50 years of exposure time is mapped. The maps of eigensolution results for 30 m, 50 m depths to the rock are illustrated. These issues are important for definition of design parameters for evaluation of seismic safety of underground infrastructure, as well as for tall structures. [ABSTRACT FROM AUTHOR]

— The 1940 (Mw=7.7) and 1977 (Mw=7.4) Vrancea earthquakes (Romania) inflicted heavy damage and casualties in Bucharest and the statistics indicate a recurrence interval of 25 years for Mw ≥ 7.0 events. Under these circumstances, the seismic microzonation represents important information for detailed urban planning that establishes an appropriate level of preparedness to the earthquake threat. This paper reviews the main studies concerning the seismicity of the Vrancea region, the site conditions of the city, the characterization of the building stock, and the codes of practice that regulate the antiseismic design. The first-order microzonation of Bucharest was performed starting from the existing database of structural and geotechnical parameters. New insights originating from direct instrumental observation and interpretation of the local effects as well as realistic numerical modeling that update and improve the input data necessary for a detailed microzoning map of the city are also discussed. [ABSTRACT FROM AUTHOR]

— The study of surface deformation due to seismic activity is often made using dislocations with uniform slip and simple geometries. A better modeling of coseismic and postseismic surface displacements can be obtained by using dislocations with variable slip and nonregular shapes. This is consistent with the asperity model of fault surfaces, assuming a friction distribution on faults made of locked zones with much higher friction than surrounding zones. In this paper we consider the 1997–1998 Colfiorito seismic sequence. The coseismic surface displacements in the Colfiorito zone are used in order to infer the slip distribution on the fault surface at different stages of the sequence. The displacement field has been modeled varying the slip distribution on the fault, and comparing the deformation observed by SAR and GPS techniques with model results. The slip distribution is calculated by Monte Carlo simulations on a normal fault with the dip angle equal to 40°. A good approximation is obtained by using square asperity units of 1.5×1.5 km2. In the first stage, we employed a simplified model with uniform slip, in which each asperity unit is allowed to slip a constant amount or not to slip at all, and in the second stage, we evaluate the slip distribution in the dislocation area determined by the Monte Carlo inversion: in this case we allow unit cells to undergo different values of slip in order to refine the initial dislocation model. The results show that the 1997 seismic events of the sequence can be modeled by irregular dislocations, obtaining a good fit to the DInSAR and GPS observations. The model also confirms the results of previous studies by a different methodology, defining the distribution of asperities on the fault plane using the fault geometry, the geodetic data and the seismic moment of the 1997–1998 Colfiorito seismic sequence. Furthermore, the analysis of 1997 aftershocks in the seismogenic region shows a strong correlation between most events and the asperity distribution, which can be considered as an independent test of the validity of the model. [ABSTRACT FROM AUTHOR]

Southwest Portugal, the Gulf of Cadiz and Morocco are under the potential threat of natural hazards linked to seismicity and tsunami generation. We report the results of two multi-channel seismic (MCS) surveys carried out in 1992 and 1998 along the continental margin and oceanic crust of SW Iberia. This MCS data set shows the evidence of the compressional deformation which involves both the continental and the oceanic crust of the study area. The area of deformation extends from the southern border of the Tagus Abyssal Plain to the Seine Abyssal Plain, encompassing the continental margin of SW Portugal. Most of the structures observed are probably related to a Mid-Miocene phase of Africa-Europe plate convergence. In this paper we discuss the recent advances on the identification of the tectonic structures that are still active and that may generate great earthquakes and tsunamis. The tectonic structures identified are located respectively at the Guadalquivir Bank, along the eastern border of the Horseshoe Abyssal Plain and along the southern continental slope of SW Portugal. [ABSTRACT FROM AUTHOR]

We analyse the seismic catalogue of the local earthquakes which occurred at Somma-Vesuvius volcano in the past three decades (1972–2000). The seismicity in this period can be described as composed of a background level, characterised by a low and rather uniform rate of energy release and by sporadic periods of increased seismic activity. Such relatively intense seismicity periods are characterised by energy rates and magnitudes progressively increasing in the critical periods. The analysis of the b value in the whole period evidences a well-defined pattern, with values of b progressively decreasing, from about 1.8 at the beginning of the considered period, to about 1.0 at present. This steady variation indicates an increasing dynamics in the volcanic system. Within this general trend it is possible to identify a substructure in the time sequence of the seismic events, formed by the alternating episodes of quiescence and activity. The analysis of the source moment tensor of the largest earthquakes shows that the processes at the seismic source are generally not consistent with simple double-couples, but that they are compatible with isotropic components, mostly indicating volumetric expansion. These components are shown to be statistically significant for most of the analysed events. Such focal mechanisms can be interpreted as the effect of explosion phenomena, possibly related to volatile exsolution from the crystallising magma. The availability of a reduced amount of high quality data necessary for the inversion of the source moment tensor, the still limited period of systematic observation of Vesuvius micro-earthquakes and, above all, the absence of eruptive events during such interval of time, cannot obviously permit the outlining of any formal premonitory signal. Nevertheless, the analysis reported in this paper indicates a progressively evolving dynamics, characterised by a generally increasing trend in the seismic activity in the volcanic system and by a significant volumetric component of recent major events, thus posing serious concern for a future evolution towards eruptive activity. [ABSTRACT FROM AUTHOR]

— Submarine and shoreline slope failures that accompany large earthquakes and large tsunamis are triggered by several mechanisms. Triggering mechanisms range from direct effects, such as inertial forces from earthquake shaking, to indirect effects, such as rapid drawdown that occurs when an earthquake-generated tsunami first approaches a shoreline. Soil shear strength also plays an important role in earthquake-related slope failures. Earthquakes change the shear strength of the soil by inducing excess pore water pressures. These excess pore water pressures change with time after the earthquake, resulting in changes in shear strength and slope stability with time. This paper reviews earthquake-related triggering mechanisms for submarine and shoreline slope failures. The variation in shear strength with time following an earthquake is examined and it is shown that delayed slope failures after an earthquake can occur as a result of changes in earthquake-induced excess pore water pressures and shear strength with time. [ABSTRACT FROM AUTHOR]

— On the evening of 17 July 1998, on th e Aitape Coast of Papua New Guinea, a strongly felt earthquake was followed some 10–25 minutes later by a destructive tsunami. The tsunami comprised three waves, each estimated to be about 4 m high. The second of the three waves rose to a height of 10–15 m above sea level after it had crossed the shoreline and caused most damage. Maximum wave heights and the greatest damage were recorded along a 14-km sector of coast centered on Sissano Lagoon. In this sector the wavefronts moved from east to west along the coast; all structures were destroyed, and in the two main villages 20–40 percent of the population was killed. Partial destruction extended 23 km to the southeast and 8 km to the northwest, and effects of the tsunami were felt as far as 250 km to the west–northwest, beyond the international border. More than 1600 people are known to have died, with some estimates as high as 2200; 1000 were seriously injured, and 10,000 survivors were displaced. This paper presents information from interviews with eye-witnesses and from mapping of damage and inundation, and includes new information on the height, shape and timing of the waves; on the possible escape of petroleum and other gases from beneath the seafloor before and during the tsunami; on unusual sound effects that preceded the waves, and lighting effects that followed; on possible deep circulation (to 250 m) of sea water in the waves; on subsidence of the order of 50–70 cm at the coastal sand barrier; and on the resilience and potential protective capacity of certain species of trees. Eye-witness accounts indicate that the tsunami reached the shore at between 09:00 and 09:08 UT, which is earlier than is proposed in published models of the timing and location of the source of the tsunami. [ABSTRACT FROM AUTHOR]

— Some underwater landslides are triggered by strong ground motions caused by earthquakes. This paper reviews current concepts and trends in the characterization of strong ground motion. Improved empirical ground motion models have been derived from a strong motion data set that has grown markedly over the past decade. However, these empirical models have a large degree of uncertainty because the magnitude-distance-soil category parameterization of these models often oversimplifies reality. This reflects the fact that other conditions that are known to have an important influence on strong ground motions, such as near-fault rupture directivity effects, crustal waveguide effects, and basin response effects, are not treated as parameters of these simple models. Numerical ground motion models based on seismological theory that include these additional effects have been developed and extensively validated against recorded ground motions, and used to estimate the ground motions of past earthquakes and predict the ground motions of future scenario earthquakes. [ABSTRACT FROM AUTHOR]

— The historical development of spontaneous rupture propagation, starting from the landmark paper of Griffith in 1920, through to the late 1980s is traced, with particular emphasis on the work carried out at MIT in the 1970s by K. Aki and his co-workers. Numerical applications of Kostrov's method for planar shear cracks were developed by Hamano, Das and Aki. Simultaneously at MIT, Madariaga considered the radiated field of a dynamic shear crack. The further development of these ideas, for example, three-dimensional spontaneous planar faulting models, continued through the 1980s. Major insight into the maximum possible rupture speeds for earthquakes developed, with the acceptance of the theoretical possibility of supersonic rupture speeds for faults with cohesion and friction, the theoretical developments spurring the search for such observations for earthquake ruptures. Possible mechanisms by which faults stop were elucidated. It was shown that a propagating rupture can jump over barriers for cracks with a cohesive zone at its tip. Complex faulting models, namely the barrier and asperity models, and their associated radiated field developed. In the late 1980s, it was shown that “dynamic” or transient asperities can develop during the complex rupturing process. Even seemingly relatively simple physical situations, can lead to such complex rupturing processes that the usual idea of “rupture velocity” needs to be abandoned in those cases. Some of the work initiated by Aki and his co-workers, such as the details of the transition from sub-Rayleigh to super-shear speeds in inplane shear mode, and the behavior of the cohesive zone size as the crack extends, still remains the subject of research today. [ABSTRACT FROM AUTHOR]

— An overview of the most important developments in Engineering (or Strong Motion) Seismology is presented alongside Professor Keiiti Aki's contributions, who is one of the founders of this field. The mechanics of earthquake rupture are discussed with due emphasis on the various physical phenomena. The presentation is made in a tutorial manner, borrowing freely from Keiiti Aki's papers, and endeavoring to emulate his unique style of clarity, simplicity and synthetic ability. [ABSTRACT FROM AUTHOR]

— A simple and powerful method for simulating ground motions is to combine parametric or functional descriptions of the ground motion's amplitude spectrum with a random phase spectrum modified such that the motion is distributed over a duration related to the earthquake magnitude and to the distance from the source. This method of simulating ground motions often goes by the name “the stochastic method.” It is particularly useful for simulating the higher-frequency ground motions of most interest to engineers (generally, f>0.1 Hz), and it is widely used to predict ground motions for regions of the world in which recordings of motion from potentially damaging earthquakes are not available. This simple method has been successful in matching a variety of ground-motion measures for earthquakes with seismic moments spanning more than 12 orders of magnitude and in diverse tectonic environments. One of the essential characteristics of the method is that it distills what is known about the various factors affecting ground motions (source, path, and site) into simple functional forms. This provides a means by which the results of the rigorous studies reported in other papers in this volume can be incorporated into practical predictions of ground motion. [ABSTRACT FROM AUTHOR]

— The phenomenon of “ Lg blockage,” where Lg is strongly attenuated by crustal heterogeneities, poses a serious problem to CTBT monitoring because Lg is an important seismic phase for discrimination. This paper examines blockage in three continental regions where the Lg blockages may be caused by large, enclosed sedimentary basins along the propagation path. The Barents Sea Basin blocks Lg propagation across the Barents Sea from the Russian nuclear test sites at Novaya Zemlya to Scandinavian stations. Also, “early Lg” waves are observed in Sn codas on NORSAR, NORESS, and ARCESS recordings of Novaya Zemlya explosions where direct Lg is blocked. Early Lg waves may have resulted from Sn-to- Lg mode conversion at the contact between the Barents Basin and the Kola Peninsula. The Northern and Southern Caspian Sea Basins also block Lg waves from PNEs and earthquakes, perhaps due to thick, low-velocity, low- Q sediments replacing the granitic layer rocks in the crust. Lg blockage has also been observed in the Western Mediterranean/Levantine Basin due to low- Q sediments and crustal thinning. A “basin capture” model is proposed to explain Lg blockage in sedimentary basins. In this model, shear waves that reverberate in the crust and constitute the Lg wave train are captured, delayed, and attenuated by thick, low-velocity sediments that replace the “granitic” layer rocks of the upper crust along part of the propagation path. Sn waves, which propagate below the basin, would not be blocked and in fact, the blocked Lg waves may be diverted downward into Sn waves by the low velocity sediments in the basin. [ABSTRACT FROM AUTHOR]

A correction to this paper has been published: https://doi.org/10.1007/s00024-021-02692-4 [ABSTRACT FROM AUTHOR]

This paper selectively reviews physical models of earthquake instability. In these models, instability arises as a result of interaction of a fault constitutive relation with deformation of the surrounding material that occurs in response to remote tectonic loading. In contrast to kinematic models in which the fault slip is imposed, it is calculated in physical models and, consequently, these models are essential for understanding precursory processes. Some kind of 'weakening' behavior for the fault constitutive relation is required to produce an instability analogous to an earthquake. Two commonly employed idealizations discussed here are rate-independent slip weakening and rate/state-dependent friction. When these constitutive models are employed on surfaces embedded in elastic half-spaces or layers, possibly coupled to a viscoelastic substrate, the results are capable of simulating realistically some aspects of earthquake occurrence. Common to all models is the prediction that earthquake instability is preceded by precursory slip which produces a departure of surface strain-rate from the background level. Near the epicenter of a moderate to large earthquake, the magnitude of this departure appears to be well within the range of current geodetic measurement accuracy, and its duration is of the order of months to years. However, details depend on a variety of factors, including the modelling of the constitutive relation near peak stress, coupling of elastic crust to the asthenosphere, and coupling of deformation with pore fluid diffusion. [ABSTRACT FROM AUTHOR]

This paper reports the present state of seismological research at Mt. Etna. A schematic classification of the earthquakes that occur on the volcano is proposed, based on both seismogram and spectrum features. We have made both focal solutions and estimates of earthquake source parameters (stress drop values between 2 and 20 bars and small source dimensions). The crust of Etna thus appears as an extremely heterogeneous medium that does not permit great stress accumulation. The coexistence of an extensional regime with an older and deeper compressive one seems confirmed at depths greater than about 7 km. Eruptive and seismic phenomena occur mainly along the principal structural trends of the volcano, but often the directions of the eruptive fractures and the earthquake concentration during the same eruption do not coincide. Tectonics seem to play an important role in controlling seismo-volcanic behaviour. [ABSTRACT FROM AUTHOR]

The physical basis for delayed aftereffects of large crustal earthquakes is analyzed. These effects apparently range in scale from typical aftershocks near the source region of a prior earthquake to the triggering of groundwater and deformation related anomalies, possibly including other large earthquakes, at locations as much as a few hundred kilometers from the source of a great earthquake. Observational reports of such aftereffects are reviewed and possible mechanisms for their occurrence are discussed. It is noted that distant coseismic stress changes associated with representative M=7 to M=8 earthquakes seem far too small to induce aftereffects of the type reported. It is shown, however, that a substantial increase of distant stresses in the crustal seismogenic layer can occur as relaxation processes take place, and approximate stress estimates are presented. These relaxation processes include, on a presumably short time scale possibly on the order of a year or less, deep aseismic viscous slip or concentrated shear flow on downward lithospheric extensions of the rupture zone and, on a longer time scale of possibly tens of years, viscous flow in the asthenosphere. The higher stresses occur in the relaxed condition first because lithospheric stress changes are then carried predominantly by the seismogenic layer, and second because the alleviation of base shear stress causes stress alterations to be channeled in the lithospheric plate, consequently causing less rapid stress attenuation with distance. The conclusion is that a fault segment in the vicinity of a large earthquake experiences a sudden coseismic stress change and, as seems more important at greater distances, an alteration in local tectonic loading rate due to the relaxation processes mentioned. In this sense a large earthquake can trigger delayed seismic or aseismic phenomena elsewhere. Using a specific experimentally motivated rate and evolving state constitutive description for frictional slip, the paper closes with an elementary analysis of the transient slippage response of a single-degree-of-freedom fault model to stress perturbations of this type. Depending on material properties, especially on whether the fault surface exhibits long-term velocity strengthening or weakening, on a normal stress-dependent measure of fault stiffness in interaction with its surrounding, and on the size of the stress perturbations and preperturbation state of the fault surface, the response may involve slip motions that accelerate to a delayed seismic instability or that transiently accelerate and then decelerate in rate, accomplishing stress relaxation in an aseismic manner. The concept of a stability domain is introduced in a phase plane whose points locate all possible conditions of the surface, and the dependence of the boundary of this domain on stiffness and constitutive properties is reported. [ABSTRACT FROM AUTHOR]

The phenomenological theory of premonitory sequence prior to an earthquake is considered. Formation of stresses related to cracks and the law governing fracture propagation lead to the evolution equation of internal stresses. A consistent description involves redistribution of the internal stresses with concentration effect before an earthquake and diffusion-like outflow after it. Formation of cracks and penetration of pore fluid cause changes in the mechanical, electrical, and other properties of rocks. Variations in the electrical resistivity of rocks are related to the fracturing processes with interaction of an interstitial fluid. Changes in rock resistivity are described here by a new parameter called a connection factor; it depends on crack density, and especially on its threshold value related to the percolation process. The influence of percolation on rock resistivity is here related to a formal analytical theory given by Kirkpatrick (1973). The author also refers to results of the time-continuous rock resistivity measurements carried out in mines in Poland. Some examples indicate a role of rock resistivity observations in situ as a precursory parameter related to rockbursts. The paper also deals with some secondary effects, like those related to a frequence response and to influence of electrokinetic potentials on an effective rock resistivity. [ABSTRACT FROM AUTHOR]

In this paper we propose a two-stage model of rock fracture. In the first stage, cracks or local regions of failure are uncorrelated and occur randomly throughout the rock in response to loading of pre-existing flaws. As damage accumulates in the rock, there is a gradual increase in the probability that large clusters of closely spaced cracks or local failure sites will develop. Based on statistical arguments, a critical density of damage will occur where clusters of flaws become large enough to lead to larger-scale failure of the rock (stage two). While crack interaction and cooperative failure is expected to occur within clusters of closely spaced cracks, the initial development of clusters is predicted based on the random variation in pre-existing flaw populations. Thus the onset of the unstable second stage in the model can be computed from the generation of random, uncorrelated damage. The proposed model incorporates notions of the kinetic (and therefore time-dependent) nature of the strength of solids as well as the discrete hierarchic structure of rocks and the flaw populations that lead to damage accumulation. The advantage offered by this model is that its salient features are valid for fracture processes occurring over a wide range of scales including earthquake processes. A notion of the rank of fracture (fracture size) is introduced, and criteria are presented for both fracture nucleation and the transition of the failure process from one scale to another. [ABSTRACT FROM AUTHOR]