Search

The spatial and temporal variation of LURR (Load/Unload Respond Ratio) in California during April 2002 to June 2004 was studied in this paper. The result shows that before the San Simeon earthquake (35.7 N, 121.1 W) on Dec. 22, 2003, Y/ Y c anomalous region occurred successively near the epicenter from April 2002 to June 2002, and the maximum anomaly of Y/ Y c occurred in May, 2002. The published research work pointed out that the Y/ Y c anomaly near the San Simeon earthquake appeared from March, 2002. Compared with the five earthquake cases out of the six with M ≥ 6.5 in California during the period from 1980 to 2001, the maximum Y/ Y c and duration of Y/ Y c anomaly before this earthquake are among the normal ranges, but the time delay from the maximum anomaly time to the occurrence time of this earthquake is the longest one. The result also shows that two areas with Y/ Y c anomalies occurred from Oct. 2002 and Dec. 2002, respectively. According to statistical characteristics of the relationship between Y/ Y c anomalies and the coming earthquakes, the seismic tendency in California was discussed in this paper. [ABSTRACT FROM AUTHOR]

Two large earthquakes occurred in the western part of China in 2008, one of them being the Yutian (35.6°N, 81.6°E) M7.3 earthquake that occurred on March 21 (BJT) and the other the Wenchuan (31.0°N, 103.4°E) M8.0 earthquake that occurred on May 12 (BJT). In this paper, the West Continental China (included in 20.0°-50.0°N, 70.0°-110.0°E region) was the study region for verifyong the predictability of the pattern informatics (PI) method using the receiver-operating characteristic curve (ROC) test and R score test. Different forecasting maps with different calculating parameters were obtained. The calculating parameters were the grid size Δ x, base time t, reference interval t to t, change interval t to t, and forecasting interval t to t. In this paper, the base time t fixed to June 1, 1971, the ending forecast time t fixed to June 1, 2008, and the forecasting interval t to t changed from 1 to 10 years, and the grid sizes were chosen as 1° × 1° and 2° × 2°, respectively. The results show that the PI method could forecast the Yutian M7.3 and Wenchuan M8.0 earthquakes only using suitable parameters. Comparing the forecast results of grid sizes 1° × 1° and 2° × 2°, the models with 2° × 2° grids were better. Comparing the forecast results with different forecasting windows from 1 to 10 years, the models with forecasting windows of 4-8 years were better using the ROC test, and the models with forecasting windows of 7-10 years were better using the R score test. The forecast efficiency of the model with a grid size of 2° × 2° and forecast window of 8 years was the best one using either the ROC test or the R score test. [ABSTRACT FROM AUTHOR]

There are many studies on the anomalous variations of the ionospheric TEC prior to large earthquakes. However, whether or not the morphological characteristics of the TEC anomalies in the daytime and at night are different is rarely studied. In the present paper, based on the total electron content (TEC) data from the global ionosphere map (GIM), we carry out a statistical survey on the spatial-temporal distribution of TEC anomalies before 1339 global M6.0+ earthquakes during 2003-2014. After excluding the interference of geomagnetic disturbance, the temporal and spatial distributions of ionospheric TEC anomalies prior to the earthquakes in the daytime and at night are investigated and compared. Except that the nighttime occurrence rates of the pre-earthquake ionospheric anomalies (PEIAs) are higher than those in the daytime, our analysis has not found any statistically significant difference in the spatial-temporal distribution of PEIAs in the daytime and at night. Moreover, the occurrence rates of pre-earthquake ionospheric TEC both positive anomalies and negative anomalies at night tend to increase slightly with the earthquake magnitude. Thus, we suggest that monitoring the ionospheric TEC changes at night might be a clue to reveal the relation between ionospheric disturbances and seismic activities. [ABSTRACT FROM AUTHOR]

Tidal stresses are generated in any three-dimensional body influenced by an external inhomogeneous gravity field of rotating planets or moons. In this paper, as a special case, stresses caused within the solid Earth by the body tides are discussed from viewpoint of their influence on seismic activity. The earthquake triggering effects of the Moon and Sun are usually investigated by statistical comparison of tidal variations and temporal distribution of earthquake activity, or with the use of mathematical or experimental modelling of physical processes in earthquake prone structures. In this study, the magnitude of the lunisolar stress tensor in terms of its components along the latitude of the spherical surface of the Earth as well as inside the Earth (up to the core-mantle boundary) were calculated for the PREM (Dziewonski and Anderson in Phys Earth Planet Inter 25(4):297-356, 1981). Results of calculations prove that stress increases as a function of depth reaching a value around some kPa at the depth of 900-1500 km, well below the zone of deep earthquakes. At the depth of the overwhelming part of seismic energy accumulation (around 50 km) the stresses of lunisolar origin are only (0.0-1.0)·103 Pa. Despite the fact that these values are much smaller than the earthquake stress drops (1-30 MPa) (Kanamori in Annu Rev Earth Planet Sci 22:207-237, 1994) this does not exclude the possibility of an impact of tidal forces on outbreak of seismic events. Since the tidal potential and its derivatives are coordinate dependent and the zonal, tesseral and sectorial tides have different distributions from the surface down to the CMB, the lunisolar stress cannot influence the break-out of every seismological event in the same degree. The influencing lunisolar effect of the solid earth tides on earthquake occurrences is connected first of all with stress components acting parallel to the surface of the Earth. The influence of load tides is limited to the loaded area and its immediate vicinity. [ABSTRACT FROM AUTHOR]

In this paper we consider the statistics of the aftershock sequence of the m = 7.65 20 September 1999 Chi-Chi, Taiwan earthquake. We first consider the frequency-magnitude statistics. We find good agreement with Gutenberg-Richter scaling but find that the aftershock level is anomalously high. This level is quantified using the difference in magnitude between the main shock and the largest inferred aftershock $$ {{\Updelta}}m^{ *}. $$ Typically, $$ {{\Updelta}}m^{ *} $$ is in the range 0.8-1.5, but for the Chi-Chi earthquake the value is $$ {{\Updelta}}m^{ *} $$ = 0.03. We suggest that this may be due to an aseismic slow-earthquake component of rupture. We next consider the decay rate of aftershock activity following the earthquake. The rates are well approximated by the modified Omori's law. We show that the distribution of interoccurrence times between aftershocks follow a nonhomogeneous Poisson process. We introduce the concept of Omori times to study the merging of the aftershock activity with the background seismicity. The Omori time is defined to be the mean interoccurrence time over a fixed number of aftershocks. [ABSTRACT FROM AUTHOR]

There are limitations in conventional earthquake loss procedures if attempts are made to apply these to assess the social and economic impacts of recent disastrous earthquakes. This paper addresses the need to develop an applicable model for estimating the significant increases of earthquake loss in mainland China. The casualties of earthquakes were studied first. The casualties of earthquakes are strongly related to earthquake strength, occurrence time (day or night) and the distribution of population in the affected area. Using data on earthquake casualties in mainland China from 1980 to 2000, we suggest a relationship between average losses of life and the magnitude of earthquakes. Combined with information on population density and earthquake occurrence times, we use these data to give a further relationship between the loss of life and factors like population density, intensity and occurrence time of the earthquake. Earthquakes that occurred from 2001 to 2003 were tested for the given relationships. This paper also explores the possibility of using a macroeconomic indicator, here GDP (Gross Domestic Product), to roughly estimate earthquake exposure in situations where no detailed insurance or similar inventories exist, thus bypassing some problems of the conventional method. [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]

— 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]

Summarizes the main results achieved in the framework of the UNESCO-IUGS-IGCP project 414 "Realistic Modelling of Seismic Input for Megacities and Large Urban Areas," centered at the Abdus Salam International Center for Theoretical Physics. Attempt to address the problems of pre-disaster orientation like hazard prediction, risk assessment, and hazard mapping, in connection with seismic activity. Unequal distribution of damage and casualties from earthquakes; Studies on the estimation of the seismic ground motion before the occurrence of a damaging earthquake.

The paper studies the statistical properties of deep earthquakes around North Island, New Zealand. We first evaluate the catalogue coverage and completeness of deep events according to cusum (cumulative sum) statistics and earlier literature. The epicentral, depth, and magnitude distributions of deep earthquakes are then discussed. It is worth noting that strong grouping effects are observed in the epicentral distribution of these deep earthquakes. Also, although the spatial distribution of deep earthquakes does not change, their occurrence frequencies vary from time to time, active in one period, relatively quiescent in another. The depth distribution of deep earthquakes also hardly changes except for events with focal depth less than 100 km. On the basis of spatial concentration we partition deep earthquakes into several groups-the Taupo-Bay of Plenty group, the Taranaki group, and the Cook Strait group. Second-order moment analysis via the two-point correlation function reveals only very small-scale clustering of deep earthquakes, presumably limited to some hot spots only. We also suggest that some models usually used for shallow earthquakes fit deep earthquakes unsatisfactorily. Instead, we propose a switching Poisson model for the occurrence patterns of deep earthquakes. The goodness-of-fit test suggests that the time-varying activity is well characterized by a switching Poisson model. Furthermore, detailed analysis carried out on each deep group by use of switching Poisson models reveals similar time-varying behavior in occurrence frequencies in each group. [ABSTRACT FROM AUTHOR]

The first part of the paper defines the terms and classifications common in earthquake prediction research and applications. This is followed by short reviews of major earthquake prediction programs initiated since World War II in several countries, for example the former USSR, China, Japan, the United States, and several European countries. It outlines the underlying expectations, concepts, and hypotheses, introduces the technologies and methodologies applied and some of the results obtained, which include both partial successes and failures. Emphasis is laid on discussing the scientific reasons why earthquake prediction research is so difficult and demanding and why the prospects are still so vague, at least as far as short-term and imminent predictions are concerned. However, classical probabilistic seismic hazard assessments, widely applied during the last few decades, have also clearly revealed their limitations. In their simple form, they are time-independent earthquake rupture forecasts based on the assumption of stable long-term recurrence of earthquakes in the seismotectonic areas under consideration. Therefore, during the last decade, earthquake prediction research and pilot applications have focused mainly on the development and rigorous testing of long and medium-term rupture forecast models in which event probabilities are conditioned by the occurrence of previous earthquakes, and on their integration into neo-deterministic approaches for improved time-variable seismic hazard assessment. The latter uses stress-renewal models that are calibrated for variations in the earthquake cycle as assessed on the basis of historical, paleoseismic, and other data, often complemented by multi-scale seismicity models, the use of pattern-recognition algorithms, and site-dependent strong-motion scenario modeling. International partnerships and a global infrastructure for comparative testing have recently been developed, for example the Collaboratory for the Study of Earthquake Predictability (CSEP) with test regions in California, Italy, Japan, New Zealand, and the Western Pacific. Algorithms and data bases are operated in a permanently learning and upgrading mode. Future perspectives and research requirements and the feasibility and possible problems encountered with the implementation of earthquake predictions in practice are briefly summarized. [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]

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 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]

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 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]

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 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]

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 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]

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 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]

— 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]

The Region–Time–Length (RTL) algorithm has been applied to different instrumental catalogues to detect seismic quiescence before medium-to-large earthquakes in Italy in the last two decades. RTL performances are sensitive to the choice of spatial and temporal parameters. The method for automatic parameters selection developed by Chen and Wu has been applied to twelve Italian earthquakes with magnitude greater than 5. The limits of the method in constructing maps of seismic quiescence before the earthquake are demonstrated, and a simple improvement is proposed. Then a new technique, namely RTL, is proposed for routine surveys of the Italian seismicity. RTL has been applied to all the earthquakes with magnitude greater than 4 in the Italian area in the time interval 1994–2004; four different sub-areas have been identified, with different characteristics in the level of recorded seismicity. One subarea—Tyrrhenian Sea—was characterized by a too low level of recorded seismicity for the application of the method. In the other three subareas a seismic quiescence was detected before at least the 66% of the earthquakes with magnitude greater or equal to 4 and all the earthquakes with magnitude greater than 5. [ABSTRACT FROM AUTHOR]

The general characteristics of the recent horizontal crustal movement in the Chinese Mainland and the horizontal crustal movements before the Kunlun Mountain earthquake of M8.1 on Nov. 14, 2001 and two earthquakes of M 8.0 and M 7.9 occurred in Sept. 2003 around the Chinese Mainland are analyzed with GPS data obtained in the Crustal Movement Observation Network of China (CMONOC). Relative vertical displacements observed in the last five years at fiducial stations show that there was a significant correlation between the length of day and the vertical displacements at most stations in western China. [ABSTRACT FROM AUTHOR]

Temporal characteristics of the famous Matsushiro earthquake swarm were investigated quantitatively using point-process analysis. Analysis of the earthquake occurrence rate revealed not only the precise and interesting process of the swarm, but also the relation between pore water pressure and the strength of the epidemic effect, and the modified Omori-type temporal decay of earthquake activity. The occurrence rate function ?(t) for this swarm is represented well aswhere f(t) represents the contribution of the swarm driver, which was the erupting water from the deep in this case, and the second term represents an epidemic effect of the modified Omori type. Based on changes in the form of f(t), this two-year long swarm was divided into six periods and one short transitional epoch. The form of f(t) in each period revealed the detail of the water erupting process. In the final stage, f (t) decayed according to the modified Omori-formula form, while it decayed exponentially in the brief respite of the water eruption in the fourth period. When an exponential decay of swarm activity is observed, we have to be cautious of a sudden restart of the violent activity. The epidemic effect is stronger when the pressure of the pore water is higher. Even when the pressure is not high, thepvalue in the epidemic effect is small, when there is plenty of pore water. However, the epidemic effect produced about a quarter of the earthquakes even though there was not much pore water in the rocks. [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]

Strong mining-induced earthquakes exhibit various aftershock patterns. The aftershock productivity is governed by the geomechanical properties of rock in the seismogenic zone, mining-induced stress and coseismic stress changes related to the main shock's magnitude, source geometry and focal mechanism. In order to assess the quantitative aftershock productivity potential in the mining environment we apply a forecast model based on natural seismicity properties, namely constant tectonic loading and the Gutenberg-Richter frequency-magnitude distribution. Although previous studies proved that mining-induced seismicity does not obey the simple power law, here we apply it as an approximation of seismicity distribution to resolve the number of aftershocks, not considering their magnitudes. The model used forecasts the aftershock productivity based on the background seismicity level estimated from an average seismic moment released per earthquake and static stress changes caused by a main shock. Thus it accounts only for aftershocks directly triggered by coseismic process. In this study we use data from three different mines, Mponeng (South Africa), Rudna and Bobrek (Poland), representing different geology, exploitation methods and aftershock patterns. Each studied case is treated with individual parameterization adjusted to the data specifics. We propose the modification of the original model, i.e. including the non-uniformity of M , resulting from spatial correlation of mining-induced seismicity with exploitation. The results show that, even when simplified seismicity distribution parameters are applied, the modified model predicts the number of aftershocks for each analyzed case well and accounts for variations between these values. Such results are thus another example showing that coseismic processes of mining-induced seismicity reflect features of natural seismicity and that similar models can be applied to study the aftershock rate in both the natural and the mining environment. [ABSTRACT FROM AUTHOR]

In this work, we present a first report on magnetic disturbances over the Brazilian sector that occurred on September 16-17, 2015. These disturbances are associated with the Illapel earthquake. Using seven magnetometers installed across Brazil, and distributed within 4500 km from the epicenter, the presence of magnetic disturbances is noted ~15 min after the time of earthquake. The seismometers registered the arrival of the first seismic waves ~10 min prior to the observed signatures on magnetic disturbances, suggesting that such waves are possibly from seismogenic origin, which involves the lithosphere-atmosphere-ionosphere coupling. [ABSTRACT FROM AUTHOR]

The recent earthquakes caused severe damages on the existing buildings. By this motivation, an important amount of research work has been conducted to determine the seismic risk of seismically active regions. For an accurate seismic risk assessment, processing of ground motions would provide an advantage. Using the current technology, it is not possible to precisely predict the future earthquakes. Therefore, most of the current seismic risk assessment methodologies are based on statistical evaluation by using recurrence and magnitude of the earthquakes hit the specified region. Because of the limited number of records on earthquakes, the quality of definitions is questionable. Fuzzy logic algorithm can be used to improve the quality of the definition. In the present study, ground motion data profile of western Turkey is defined using an intelligent hybrid processing. The approach is given in a practical way for an easier and faster calculation. Earthquake data between 1970 and 1999 from western part of Turkey have been used for training. The results are tested and validated with the earthquake data between 2000 and 2015 of the same region. Enough approximation was validated between calculated values and the earthquake data by using the intelligent hybrid processing. [ABSTRACT FROM AUTHOR]

On 16 September 2015, a great ( M 8.3) interplate thrust earthquake ruptured offshore Illapel, Chile, producing a 4.7-m local tsunami. The last major rupture in the region was a 1943 M 7.9 event. Seismic methods for rapidly characterizing the source process, of value for tsunami warning, were applied. The source moment tensor could be obtained robustly by W-phase inversion both within minutes (Chilean researchers had a good solution using regional data within 5 min) and within an hour using broadband seismic data. Short-period teleseismic P wave back-projections indicate northward rupture expansion from the hypocenter at a modest rupture expansion velocity of 1.5-2.0 km/s. Finite-fault inversions of teleseismic P and SH waves using that range of rupture velocities and a range of dips from 16°, consistent with the local slab geometry and some moment tensor solutions, to 22°, consistent with long-period moment tensor inversions, indicate a 180- to 240-km bilateral along-strike rupture zone with larger slip northwest to north of the epicenter (with peak slip of 7-10 m). Using a shallower fault model dip shifts slip seaward toward the trench, while a steeper dip moves it closer to the coastline. Slip separates into two patches as assumed rupture velocity increases. In all cases, localized ~5 m slip extends down-dip below the coast north of the epicenter. The seismic moment estimates for the range of faulting parameters considered vary from 3.7 × 10 Nm (dip 16°) to 2.7 × 10 Nm (dip 22°), the static stress drop estimates range from 2.6 to 3.5 MPa, and the radiated seismic energy, up to 1 Hz, is about 2.2-3.15 × 10 J. [ABSTRACT FROM AUTHOR]

Seismic coda represents a valuable source of information about the attenuation of the high-frequency waves in the studied region. The quality factor Q derived from coda is an integral parameter of the volume surrounding the hypocenter and seismic station and, according to the applied method, represents the total attenuation or the intrinsic and scattering parts. We analyzed records of 13 selected earthquakes in the magnitude range 1.7-2.9 of the 2011 swarm from West Bohemian/Vogtland area (central Europe), which were recorded at epicentral distances from 7 to 50 km. Two methods were applied: coda method for estimation of the Q and the Multiple Lapse Time Windows Analysis for separation of the scattering and intrinsic loss by estimation of Q and Q. Careful selection of the analyzed events was necessary due to the frequent contamination of coda decays by the running seismic swarm activity. The resulting coda Q is relatively high with respect to the geodynamic activity and varies between 100 and 2500 within the analyzed frequency range of 1-18 Hz. The intrinsic loss dominates over scattering attenuation with Q increasing from 100 and 1850 and Q from 300 to 3400 in the same frequency range, which is consistent with the geodynamic activity of the region. We find that the intrinsic attenuation in West-Bohemia/Vogtland is higher that in neighboring Germany, which could be attributed to the heterogeneity of the crust in central Europe. [ABSTRACT FROM AUTHOR]

At The Geysers geothermal reservoir in northern California, evidence strongly suggests that activities associated with production of electric power cause an increase in the number of small earthquakes. First-degree dynamic moment tensors are used to investigate the relationship between induced earthquakes and injection of water into a well as part of a controlled experiment in the northwest Geysers. The estimation of dynamic moment tensors in the complex shallow crust at The Geysers is challenging, so the method is described in detail with particular attention given to the uncertainty in the results. For seismic events in the moment magnitude range of 0.9-2.8, spectral moduli of dynamic moment tensors are reliably recovered in the frequency range of 1-100 Hz, but uncertainty in the associated spectral phases limits their use to a few simple results. A number of different static moment tensors are investigated, with the preferred one obtained from parameters of a model fitted to the spectral modulus of the dynamic moment tensor. Moment tensors estimated for a group of 20 earthquakes exhibit a range of source mechanisms, with over half having significant isotropic parts of either positive or negative sign. Corner frequencies of the isotropic part of the moment tensor are about 40 % larger than the average of the deviatoric moment tensor. Some spatial patterns are present in source mechanisms, with earthquakes closely related in space tending to have similar mechanisms, but at the same time, some nearby earthquakes have very different mechanisms. Tensional axes of displacement in the source regions are primarily horizontal, while the pressure axes range from near horizontal to vertical. Injection of water into the well in the center of the study area clearly causes an increase in the number of earthquakes per day, but an effect upon source mechanisms is not evident. [ABSTRACT FROM AUTHOR]

In this article, a stochastic finite fault source model is calibrated to estimate ground motion in northeastern India for intermediate depth events originating in the Indo-Burmese tectonic domain. A total of 47 three-component accelerograms from eight events with magnitudes ranging from M 4.8-6.4 are used to estimate the input source and site parameters of the finite fault source model. Key seismic parameters such as stress drop (Δσ) and site amplification function are determined from the recorded strong motion data. The obtained stress drop of the eight recorded events lies in between 105 and 165 bars. [ABSTRACT FROM AUTHOR]

The collision between the Arabian and Eurasian plates in eastern Turkey causes the Anatolian block to move westward. The North Anatolian Fault (NAF) is a major strike-slip fault that forms the northern boundary of the Anatolian block, and the Erzincan Basin is the largest sedimentary basin on the NAF. In the last century, two large earthquakes have ruptured the NAF within the Erzincan Basin and caused major damage ( M = 8.0 in 1939 and M = 6.8 in 1992). The seismic hazard in Erzincan from future earthquakes on the NAF is significant because the unconsolidated sedimentary basin can amplify the ground motion during an earthquake. The amount of amplification depends on the thickness and geometry of the basin. Geophysical constraints can be used to image basin depth and predict the amount of seismic amplification. In this study, the basin geometry and fault zone structure were investigated using broadband magnetotelluric (MT) data collected on two profiles crossing the Erzincan Basin. A total of 24 broadband MT stations were acquired with 1-2 km spacing in 2005. Inversion of the MT data with 1D, 2D and 3D algorithms showed that the maximum thickness of the unconsolidated sediments is ~3 km in the Erzincan Basin. The MT resistivity models show that the northern flanks of the basin have a steeper dip than the southern flanks, and the basin deepens towards the east where it has a depth of 3.5 km. The MT models also show that the structure of the NAF may vary from east to west along the Erzincan Basin. [ABSTRACT FROM AUTHOR]

Finite difference simulations of seismic wave propagation are performed in the Niigata area, Japan, for the 2007 Mw 6.6 Niigata-ken Chuetsu-Oki earthquake at low frequencies. We test three 3D structural models built independently in various studies. First aftershock simulations are carried out. The model based on 3D tomography yields correct body waves in the near field, but later phases are imperfectly reproduced due to the lack of shallow sediment layers; other models based on various 1D/2D profiles and geological interpretation provide good site responses but generate seismic phases that may be shifted from those actually observed. Next, for the mainshock simulations, we adopt two different finite source models that differ in the near-field ground motion, especially above the fault plane (but under the sea) and then along the coastline. Each model is found to be calibrated differently for the given stations. For engineering purposes, the variations observed in simulated ground motion are significant, but for seismological purposes, additional parameter calibrations would be possible for such a complex 3D case. [ABSTRACT FROM AUTHOR]

Although the evidence for complexity is overwhelming, the dynamics of faulting is still poorly understood. Whilst it has long been known that discreteness in numerical earthquake models produces complexity, the mathematical structure and form of this complexity has never been fully established. Using a simple 1D nonlinear fault model we show how complexity can arise in discrete models through the presence of nonlinear, scale-dependent (or mesh-dependent) terms. We show that scale-dependencies may be a significant factor in the generation of slip complexity and pulse-like rupture over multiple earthquake cycles. We demonstrate that the introduction of length scales in discrete earthquake models implies that both strongly weakening friction and scale-dependent processes may be necessary in generating the pulse-like rupture mode and earthquake complexity over multiple earthquake cycles. [ABSTRACT FROM AUTHOR]

We investigated the development of a distinct later phase observed at stations near the Japan Trench associated with shallow, outer-rise earthquakes off the coast of Sanriku, northern Japan based on the analysis of three-component broadband seismograms and FDM simulations of seismic wave propagation using a heterogeneous structural model of the Japan Trench subduction zone. Snapshots of seismic wave propagation obtained through these simulations clearly demonstrate the complicated seismic wavefield constructed by a coupling of the ocean acoustic waves and the Rayleigh waves propagating within seawater and below the sea bottom by multiple reflections associated with shallow subduction zone earthquakes. We demonstrated that the conversion to the Rayleigh wave from the coupled ocean acoustic waves and the Rayleigh wave as they propagate upward along the slope of seafloor near the coast is the primary cause of the arrival of the distinct later phase at the station near the coast. Through a sequence of simulations using different structural models of the Japan Trench subduction zone, we determined that the thick layer of seawater along the trench and the suddenly rising sea bottom onshore of the Japanese island are the major causes of the distinct later phase. The results of the present study indicate that for realistic modeling of seismic wave propagation from the subduction zone earthquakes, a high-resolution bathymetry model is very crucial, although most current simulations do not include a water column in their simulation models. [ABSTRACT FROM AUTHOR]

Prior to an earthquake, natural seismicity is correlated across multiple spatial and temporal scales. Many studies have indicated that an earthquake is hard to accurately predict by a single time-dependent precursory method. In this study, we attempt to combine four earthquake prediction methods, i.e. the Pattern Informatics (PI), Load/Unload Response Ratio (LURR), State Vector (SV), and Accelerating Moment Release (AMR) to estimate future earthquake potential. The PI technique is founded on the premise that the change in the seismicity rate is a proxy for the change in the underlying stress. We first use the PI method to quantify localized changes surrounding the epicenters of large earthquakes to objectively quantify the anomalous areas (hot spots) of the upcoming events. Next, we delineate the seismic hazard regions by integrating with regional active fault zones and small earthquake activities. Then, we further evaluate the earthquake potential in the seismic hazard regions using the LURR, SV and AMR methods. Retrospective tests of this new approach on the large earthquakes ( M > 6.5) which have occurred in western China over the last 3 years show that the LURR and SV time series usually climb to an anomalously high peak months to years prior to occurrence of a large earthquake. And, the asymptote time, t, 'predicted' by the AMR method correspond to the time of the actual events. The results may suggest that the multi-methods combined approach can be a useful tool to provide stronger constraints on forecasts of the time and location of future large events. [ABSTRACT FROM AUTHOR]

We explore fractal properties of two observed seismicity distributions prior to the 2003 M 7.4 Colima, Mexico and 1992 M 7.3 Landers, USA earthquakes, together with several mathematical fractal distributions and two non-fractal ones, in order to estimate minimum reliable sample sizes, determine whether fractality for observed seismicity is essentially different from random uniform distributions, and explore the possibility of extracting premonitory information from fractal characteristics of seismicity before large earthquakes. Sample sizes above 800 events for whole catalogs appear to be sufficient to maintain ordered multifractality and to yield dimension estimates that vary smoothly and reliably. Fractal estimates appear to be best for whole catalogs that include aftershocks. The fractal characteristics of spatial distributions of seismicity are essentially different from those of the uniform random distribution, which is the null hypothesis of a non-fractal distribution with minimum information. The fractal dimensions and afractality measures of seismicity distributions change with time and show distinctive behaviors associated with foreshocks and main events, although these behaviors are different for each example. Results suggest the possibility of a priori identification of foreshocks to large earthquakes. A combination of fractal dimension and afractality measures over time may be helpful in large earthquake premonitory studies. [ABSTRACT FROM AUTHOR]

We characterize the heterogeneous source slip model of intraslab earthquakes to compare source scaling properties with those of inland crustal and subduction-zone plate-boundary earthquakes. We extracted rupture area ( S), total area of asperity ( S), average slip ( D) and average slip on asperity ( D) of eleven intraslab earthquakes following the procedure proposed by S omerville et al. (Seism Res Lett 70:59-80, ) and proposed the empirical scaling relationship formula of S, S, and D for intraslab earthquakes. Under the same seismic moment, an intraslab earthquake has a smaller rupture area and total area of asperity, and smaller average slip than an inland crustal earthquake. The area ratio of asperity area and total rupture area of intraslab earthquakes are similar to those of inland crustal earthquakes. The strong motion generation area (SMGA) scaling of intraslab earthquakes appears self-similar, and those results support the idea the characterized source model of intraslab earthquakes can be modeled in a manner similar to that of inland crustal earthquakes. [ABSTRACT FROM AUTHOR]

wavetrain of high-frequency (HF) P waves from a large earthquake, when recorded at a distant station, looks like a segment of modulated noise, with its duration close to the duration of rupture. These wavetrains, with their bursts and fadings, look much more intermittent than a segment of common stationary random noise. We try to describe quantitatively this bursty behavior. To this end, variogram and spectral analyses are applied to time histories of P-wave envelopes (squared-amplitude or instant-power signals) in six HF bands of 1-Hz width. Nine M = 7.6-9.2 earthquakes were examined, using, in total, 232 records and 992 single-band traces. Variograms of integrated instant power are approximately linear on a log-log scale, indicating that the correlation structure of the instant-power signal is approximately self-similar. Also, estimates of the power spectrum of the instant-power signal look approximately linear on a log-log scale. Log-log slopes of the variograms and spectra deliver estimates of the Hurst exponent H that are mostly in the range 0.6-0.9, markedly above the value H = 0.5 of uncorrelated (white-noise) signals. The preferred estimate over the entire data set is H = 0.83, still, this estimate may include some bias, and must be treated as preliminary. The inter-event scatter of H estimates is about 0.04, reflecting individual event-to-event variations of H. Many of the average log-log spectral plots show slight concavity that perturbs the approximately linear slope; this is a secondary effect that seems to be mostly related to the limited bandwidth of the data. Evidence is given in support of the idea that the observed approximately self-similar correlation structure of the P-wave envelope originates in a similar structure of the body wave instant-power signal radiated by the source, so that the propagation-related distortions can be regarded as limited. The facts presented suggest that the space-time organization of the earthquake rupture process is multiscaled and bears significant fractal features; it deviates from the brittle-crack model with its two well-separated characteristic scales. Phenomenologically, the high-frequency body-wave radiation from an earthquake source can be thought of as a product of stationary noise and the square root of a positive random envelope function with a power-law spectrum. From the viewpoint of applications, the self-similarity of body wave envelopes provides a useful constraint for earthquake source models used to simulate strong ground motions. [ABSTRACT FROM AUTHOR]

In weather forecasting, current and past observational data are routinely assimilated into numerical simulations to produce ensemble forecasts of future events in a process termed “model steering”. Here we describe a similar approach that is motivated by analyses of previous forecasts of the Working Group on California Earthquake Probabilities (WGCEP). Our approach is adapted to the problem of earthquake forecasting using topologically realistic numerical simulations for the strike-slip fault system in California. By systematically comparing simulation data to observed paleoseismic data, a series of spatial probability density functions (PDFs) can be computed that describe the probable locations of future large earthquakes. We develop this approach and show examples of PDFs associated with magnitude M > 6.5 and M > 7.0 earthquakes in California. [ABSTRACT FROM AUTHOR]

We employ a computationally efficient fault system earthquake simulator, RSQSim, to explore effects of earthquake nucleation and fault system geometry on earthquake occurrence. The simulations incorporate rate- and state-dependent friction, high-resolution representations of fault systems, and quasi-dynamic rupture propagation. Faults are represented as continuous planar surfaces, surfaces with a random fractal roughness, and discontinuous fractally segmented faults. Simulated earthquake catalogs have up to 10 earthquakes that span a magnitude range from ∼M4.5 to M8. The seismicity has strong temporal and spatial clustering in the form of foreshocks and aftershocks and occasional large-earthquake pairs. Fault system geometry plays the primary role in establishing the characteristics of stress evolution that control earthquake recurrence statistics. Empirical density distributions of earthquake recurrence times at a specific point on a fault depend strongly on magnitude and take a variety of complex forms that change with position within the fault system. Because fault system geometry is an observable that greatly impacts recurrence statistics, we propose using fault system earthquake simulators to define the empirical probability density distributions for use in regional assessments of earthquake probabilities. [ABSTRACT FROM AUTHOR]