On-chip earthquake simulation model using potentials

A two-dimensional (2-D) Cellular Automata (CA) dynamic system constituted of cells-charges has been proposed for the simulation of the earthquake process. The CA model has been calibrated with the use of real data. The calibration incorporates major seismic characteristics of the area under test. The simulation results are found in good quantitative and qualitative agreement with the recorded Gutenberg–Richter (GR) scaling relations. The model is enriched with a powerful multi-parameter interface that enables the user to load real data from different regions. This paper examines the on-chip realisation of the model and its instrumentation. The CA model hardware implementation is based on Field Programmable Gate Array (FPGA) logic. It utilises an array of 32 × 32 cells. Parameters that construct the local CA rule constitute the input data. The initial seed, which to some extent corresponds to the seismic features of the area under test, is loaded in a semi-parallel way and the process is completed in a certain number of time steps. The automatic response of the processor provides the corresponding GR scaling law of the area under study. The hardware implementation of the CA-based earthquake simulation model is advantageous in terms of low-cost, high-speed, compactness and portability features. It can operate as a preliminary data-acquisition filter that accelerates the evaluation of recorded data as far as its origin time, spatial and magnitude completeness and quality are concerned. Software that performs reliable automatic phase picking, as well as data elaboration, can be assembled next to the earthquake recording instruments (the whole network) output to assure a quick and reliable iteration of the focal parameters of a recorded earthquake (epicentre coordinates, focal depth and magnitude). The dedicated processor can be accommodated right after this stage (before any manual elaboration) focusing on the near real-time development of a reliable qualitative dynamical seismic record and a mapping of the seismic characteristics of the area.     

Vp/Vs Anisotropy and Implications for Crustal Composition Identification and Earthquake Prediction

Abstract: The ratio of P- to S-wave velocities (Vp/Vs) is regarded as one of the most diagnostic properties of natural rocks. It has been used as a discriminant of composition for the continental crust and provides valuable constraints on its formation and evolution processes. Furthermore, the spatial and temporal changes in Vp/Vs before and after earthquakes are probably the most promising avenue to understanding the source mechanics and possibly predicting earthquakes. Here we calibrate the variations in Vp/Vs in dry, anisotropic crustal rocks and provide a set of basic information for the interpretation of future seismic data from the Wenchuan earthquake Fault zone Scientific Drilling (WFSD) project and other surveys. Vp/Vs is a constant (Ф0) for an isotropic rock. However, most of crustal rocks are anisotropic due to lattice-preferred orientations of anisotropic minerals (e.g., mica, amphibole, plagioclase and pyroxene) and cracks as well as thin compositional layering. The Vp/Vs ratio of an anisotropic rock measured along a selected pair of propagation-vibration directions is an apparent value (Фij) that is significantly different from the value for its isotropic counterpart (Ф0). The usefulness of apparent Vp/Vs ratios as a diagnostic of crustal composition depends largely on rock seismic anisotropy. A 5% of P- and S-wave velocity anisotropy is sufficient to make it impossible to determine the crustal composition using the conventional criteria (Vp/Vs≤1.756 for felsic rocks, 1.7561.944 fluid-filled porous/fractured or partially molten rocks) if the information about the wave propagation-polarization directions with respect to the tectonic framework is unknown. However, the variations in Vp/Vs measured from borehole seismic experiments can be readily interpreted according to the orientations of the ray path and the polarization of the shear waves with respect to the present-day principal stress directions (i.e., the orientation of cracks) and the frozen fabric (i.e., foliation and lineation).     

龙池乡高位滑坡抛射机理分析

5.12汶川大地震龙门山地区发育地震高位滑坡,通常发生在10~11度的极震区,这种灾害难预测,发生速度快,容易造成惨重的生命和财产损失。本文通过在龙池乡川主坪的实地调查,对典型地震抛射滑坡形成的抛射堆积体进行剖析,探究地震加速度受地形放大效应的影响效应及崩塌被抛射的运动程式。     

The spatial and temporal changes inV p /V s before strong earthquakes in Kamchatka

A method was developed to obtain from a signal station the spatial and temporal distribution ofV _p /V _s ratios before earthquakes of magnitude>6. It was shown thatV _p /V _s values strongly depend upon the relative positions of the stations, the future large earthquake and the foci of the smaller earthquakes used forV _p /V _s determination. The appearance of a zone of anomalousV _p /V _s values with linear dimensions of the order of 100 km was noted at least 4 years before a deep earthquake of magnitude 7. Similar size anomalous zones were detected one year before some magnitude 6 earthquakes. V _p /V _s values decreased by a small but distinct amount during this time. Additionally, local short term minima inV _p /V _s ratios were observed some months before the major event. The epicenters of the large earthquakes were located within the 100 km size zone where the gradients of theV _p /V _s field were largest.     

Influence of earthquakes on the stability of slopes

Earthquakes are a major trigger for instability of natural and man-made slopes. Often the instability of slopes due to an earthquake causes more destruction and kills more people than the actual earthquake itself. A comparison is made between different methodologies to analyze the potential stability of slopes during earthquakes. Theoretically, it seems simple to calculate the stability of a slope during an earthquake. In reality, however, the stability is influenced by so many parameters that are either not known or which influence is so poorly known that a decent estimation of stability cannot be made. Offshore the situation is worse because proper data required for stability calculations are even less available than onshore. On- and offshore, erosion and weathering create continuously slopes that may become unstable during a future earthquake, offshore also sedimentation creates continuously new slopes. Another fundamental problem in stability analysis is the complicated and largely unknown behavior of seismic waves in three-dimensions in natural materials. The lack of accurate data and the unknown behavior of seismic waves in three-dimensions make estimations of slope stability during an earthquake unreliable.     

6600 years of earthquake record in the Julian Alps (Lake Bohinj,Slovenia)

Sequences of lake sediments often form long and continuous records that may be sensitive recorders of seismic shaking. A multi‐proxy analysis of Lake Bohinj sediments associated with a well‐constrained chronology was conducted to reconstruct Holocene seismic activity in the Julian Alps (Slovenia). A seismic reflection survey and sedimentological analyses identified 29 homogenite‐type deposits related to mass‐wasting deposits. The most recent homogenites can be linked to historical regional earthquakes (i.e. 1348 ad , 1511 ad and 1690 ad ) with strong epicentral intensity [greater than ‘damaging’ (VIII ) on the Medvedev–Sponheuer–Karnik scale]. The correlation between the historical earthquake data set and the homogenites identified in a core isolated from local stream inputs, allows interpretation of all similar deposits as earthquake related. This work extends the earthquake chronicle of the last 6600 years in this area with a total of 29 events recorded. The early Holocene sedimentary record is disturbed by a seismic event (6617 ± 94 cal yr bp ) that reworked previously deposited sediment and led to a thick sediment deposit identified in the seismic survey. The period between 3500 cal yr bp and 2000 cal yr bp is characterized by a major destabilization in the watershed by human activities that led to increases in erosion and sedimentation rates. This change increased the lake's sensitivity to recording an earthquake (earthquake‐sensitivity threshold index) with the occurrence of 72 turbidite‐type deposits over this period. The high turbidite frequency identified could be the consequence of this change in lake earthquake sensitivity and thus these turbidites could be triggered by earthquake shaking, as other origins are discarded. This study illustrates why it is not acceptable to propose a return period for seismic activity recorded in lake sediment if the sedimentation rate varies significantly.     

Process of Housing Damage Assessment: The 1995 Hanshin-Awaji Earthquake Disaster Case

Rebuilding of victims' livelihoods was a crucial issue in the restoration process in the1995 Hanshin-Awaji Earthquake Disaster. Housing damage assessment influencedmost of the rebuilding of the livelihood in the long term, because the Victim Certificatesissued by the local governments based on the results of the Housing damage assessmentwas required to receive most of the individual assistance measures. In the process ofHousing damage assessment, many complex problems arose, leading to extensivework on the part of the disaster responders. Consequently, a considerable number ofvictims were dissatisfied with the assessment and applied for a resurvey. Due to a floodof requests for resurvey, disaster responders had to work on damage assessment, leavingrelief activities aside.In order to facilitate Housing damage assessment, this paper discusses thefollowing five points: (1) the processes and the problems of assessments performedin the Hanshin-Awaji Earthquake Disaster, (2) the changes in the nature of informationneeded by the victims, (3) the improvements over the present damage assessment, (4)the housing situation in Japan, and (5) the international situation on damage assessment.It is obvious from the results that a poor damage assessment system and the size ofthe disaster produced a very large work load. Differences in appreciation among theinvestigators also contributed to unfair assessments and led to the victims beingincreasingly dissatisfied by the survey results. Finally, a design concept for acomprehensive damage assessment system, which has been derived from theabove five points, is proposed for post-disaster management.     

Holocene beach deposits for assessing coastal uplift of the northeastern Boso Peninsula, Pacific coast of Japan

This paper presents a case study that assessed spatial variations in the tectonic uplift rates of beach deposits in the relict Kujukuri strand plain, situated on the northeastern coast of the Boso Peninsula, eastern Japan. The southern Boso Peninsula is tilted downward to the northeast due to plate subduction along the Sagami Trough. However, the cause of the northeastern coast uplift creating the relict strand plain is unclear, due to the absence of a Holocene raised marine terrace sequence. Elevations and ages of beach deposits were collected from drilled cores and ground-penetrating radar profiles along three shore-normal sections in the southern Kujukuri strand plain. From this, alongshore variations in the relative sea level since the mid-Holocene could be seen. These corresponded to north-to-northeast downward tilting at a rate of 0.4 m/ka for an interval 10 km and are concordant with the longer term tilting of the last interglacial marine terrace surrounding the plain. Although it is difficult to assess shore-normal variations of uplift based on the present dataset, the recognized tilting apparently continues to the tilting of the southern Boso Peninsula, implying the Sagami Trough probably affects the uplift of the Kujukuri coast.     

SOURCE RADIATION AND RESPONSES OF WAVE PROPAGATION

Recordings of seismic waves propagating from earthquake source to a station at the earth's surface are a system response function.The convolution operator in time domain can be simplified as a multiplication operator in frequency domain.We discuss in frequency domain the separation of source,path and site effects for global scaling of earthquake source radiation.Also discussed are source scaling model,faulting mechanism,and the H/V inversion problems with crustal and near surface structures.Gross features of apparent source spectra appear to be not much region-dependent although there may be difference between tectonic styles within a region of tectonic mixture for which we need further study as data accumulate.Vertical spectra may be a better approach to approximate source radiation,as it has less crustal amplification effects than horizontal spectra.The H/V ratio is evidently a comprehensive indicator of amplification effects from near surface to deep structure.This gives it potential as an inversion tool to deduce site crustal structure.