Earthquake rupture complexity due to dynamic nucleation and interaction of subsidiary faults

We numerically study the dynamic interaction of propagating cracks. It is assumed that propagating cracks can nucleate and drive subsidiary cracks because of shear strain enhancement near the propagating crack tips. The critical strain fracture criterion is assumed in the analysis. Intense interaction is expected to occur among the cracks. All the cracks are assumed to be parallel and antiplane strain deformation is assumed in the computation.In the interaction of two non-coplanar cracks, a strain shadow is formed in the neighborhood of each crack because of the strain release by the introduction of the crack. The growth of each crack is accelerated when the propagating tips of each crack are outside of the strain shadow of the other crack. In general, the crack tips enter the strain shadow, and the crack tips decelerate. The calculation shows that only one of the two cracks can continue to grow, and the other's growth is decelerated and arrested. If we can assume that the suite of cracks interact in a pairwise manner only, then this may suggest that only a limited number of cracks can continue to grow during the final stage of the rupture process. Hence the crack interaction causes complexity in dynamic earthquake faulting. The concepts of barrier and asperity have been employed by many researchers for the interpretation of complex seismic wave data. However, the physical realities of such concepts are obscure. Our calculations show that dynamic crack interactions can produce barriers and asperities in some cases; the crack tip deceleration or arrest due to the interactions among non-coplanar cracks can be interpreted as being due to a barrier. The dynamic coalescence among the coplanar cracks can be regarded as an asperity.Umeda found a localized area that strongly radiates high-frequency seismic waves in the epicentral areas of some large shallow earthquakes. He defined this as an earthquake bright spot. Our analysis implies that only a limited number of cracks continue to grow when many interactive cracks nucleate, and that all other cracks stop extending soon after nucleation. Hence, if the nucleation and termination of several cracks occur in a localized area, it will be observed seismologically as an earthquake bright spot. This is because it is theoretically known that the sudden termination of crack growth and dynamic crack coalescence efficiently emits high-frequency elastic waves.     

Scattering of elastic waves by a fracture zone containing randomly distributed cracks

We theoretically study the scattering ofP, SV andSH waves by a zonal distribution of cracks, which simulates a fault fracture zone. An investigation is conducted how the geometrical properties of the crack distribution and the frictional characteristics of the crack surface are reflected in the attenuation and dispersion of incident waves, as well as in the amplitudes of the transmitted and reflected waves from the zone. If the crack distribution within the fault zone changes temporally during the preparation process of the expected earthquake, it will be important for earthquake prediction to monitor it, utilizing the scattering-induced wave phenomena.We consider the two-dimensional problem. Aligned cracks with the same length are assumed to be randomly distributed in a zone with a finite width, on which elastic waves are assumed to be incident. The distribution of cracks is assumed to be homogeneous and sparse. The crack surface is assumed to be stress-free, or to undergo viscous friction; the latter case simulates fluid-filled cracks. The opening displacement of the crack is assumed to be negligibly small. The idea of the mean wave formalism is employed in the analysis, and Foldy's approximation is assumed.When the crack surface is stress-free, it is commonly observed for every wave mode (P, SV andSH) that the attenuation coefficientQ ^–1 peaks aroundka1, the phase velocity is almost independent ofk in the rangeka<1 and it increases monotonically withk in the rangeka>1, wherek is the intrinsicS wavenumber anda is the half length of the crack. The effect of the friction is to shift the peak ofQ ^–1 and the corner of the phase velocity curve to the low wavenumber range. The high wavenumber asymptote ofQ ^–1 is proportional tok ^–1 independently of model parameters and the wave modes. If the seismological observation thatQ ^–1 ofS waves has a peak at around 0.5 Hz in the earth's crust is combined with our results, the upper limit of crack size within the crust is estimated about 4 km. The information regarding the transmitted and reflected waves, such as the high wavenumber limit of the amplitude of the transmitted wave etc., allows estimation of the strength of the friction.     

Energy transfer among Pb, Ce and Mn in fluorescent calcite from Kuerle, Xinjiang, China

Natural calcite from Kuerle, Xinjiang, China, shows orange-red fluorescence when exposed to short-wave ultraviolet (UV) light (Hg 253.7 nm). Photoluminescence (PL) emission and excitation spectra of the calcite are observed at room temperature in detail. The PL emission spectrum under 208 nm excitation consists of three bands: two UV bands at 325 and 355 nm and an orange-red band at 620 nm. The three bands are ascribed to Pb²⁺, Ce³⁺ and Mn²⁺, respectively, as activators. The Pb²⁺ excitation band is observed at 243 nm, and the Ce³⁺ excitation band at 295 nm. The Pb²⁺ excitation band is also observed by monitoring the Ce³⁺ fluorescence, and the Pb²⁺ and Ce³⁺ excitation bands, in addition to six Mn²⁺ excitation bands, are also observed by monitoring the Mn²⁺ fluorescence. These indicate that four types of the energy transfer can occur in calcite through the following processes: (1) Pb²⁺ → Ce³⁺, (2) Pb²⁺ → Mn²⁺, (3) Ce³⁺ → Mn²⁺ and (4) Pb²⁺ → Ce³⁺ → Mn²⁺.     

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.     

Dynamic nucleation process of shallow earthquake faulting in a fault zone

Two distinct phases are commonly observed at the initial part of seismograms of large shallow earthquakes: low-frequency and low-amplitude waves following the onset of a P wave ( P ₁) are interrupted by the arrival of the second impulsive phase P₂ enriched with high-frequency components. This observation suggests that a large shallow earthquake involves two qualitatively different stages of rupture at its nucleation.
We propose a theoretical model that can naturally explain the above nucleation behaviour. The model is 2-D and the deformation is assumed to be anti-plane. A key clement in our model is the assumption of a zone in which numbers of pre-existing cracks are densely distributed; this cracked zone is a model for the fault zone. Dynamic crack growth nucleated in such a zone is intensely affected by the crack interactions, which exert two conflicting effects: one tends to accelerate the crack growth, and the other tends to decelerate it. The accelerating and decelerating effects are generally ascribable to coplanar and non-coplanar crack interactions, respectively. We rigorously treat the multiple interactions among the cracks, using the boundary integral equation method (BIEM), and assume the critical stress fracture criterion for the analysis of spontaneous crack propagation.
Our analysis shows that a dynamic rupture nucleated in the cracked zone begins to grow slowly due to the relative predominance of non-coplanar interactions. This process radiates the P₁ phase. If the crack continues to grow, coalescence with adjacent coplanar cracks occurs after a short time. Then, coplanar interactions suddenly begin to prevail and crack growth is accelerated; the P₂ phase is emitted in this process. It is interpreted that the two distinct phases appear in the process of the transition from non-coplanar to coplanar interaction predominance.     

Attenuation of coda waves and Q c value beneath the Chengdu telemetered seismic network

Based on the single scattering model of coda power spectrum analysis, digital waveform data of 50 events recorded by the real-time processing system of the Chengdu telemetry network are analyzed to estimate the Q _c values of earth medium beneath the Chengdu telemetry network for several specified frequencies. It is found that the Q _c shows the frequency dependency in the form of Q _c = Q ₀ f ⁿ in the range of 1.0 to 20.0Hz. Estimated Q ₀ ranges from 60.83 to 178.05, and n is found to be 0.713 to 1.159. The average value of Q ₀ and n are 117 and 0.978 respectively. This result indicates the strong frequency dependency of the attenuation of coda waves beneath the Chengdu telemetry network. Comparing with the results obtained in other regions of the world, it is found that Q ₀ ⁻¹ value and its change with frequency are similar to those in regions with strong tectonic activity. This subject is supported by the Ministry of Personnel, China for partly sponsoring.     

Primary Ore Mineral Assemblage and Fluid Inclusion Study of the Batu Hijau Porphyry Cu-Au Deposit, Sumbawa, Indonesia

Abstract. The Batu Hijau porphyry Cu‐Au deposit, Sumbawa Island, Indonesia, is associated with a tonalitic intrusive complex. The temperature‐pressure condition of mineralization at the Batu Hijau deposit is discussed on the basis of fluid inclusion microthermometry. Then, the initial Cu‐Fe sulfide mineral assemblage is discussed. Bornite and chalcopyrite are major copper ore minerals associated with quartz veinlets. The quartz veinlets have been classified into ‘A’ veinlets associated with bornite, digenite, chalcocite and chalcopyrite, ‘B’ veinlets having chalcopyrite bornite along vuggy center‐line, rare ‘C’ chalcopyrite‐quartz veinlets, and late ‘D’ veinlets consisting of massive pyrite and quartz (Clode et al., 1999). Copper and gold mineralization is associated with abundant ‘A’ quartz veinlets. Abundant fluid inclusions are found in veinlet quartz consisting mainly of gas‐rich inclusions and polyphase inclusions throughout the veinlet types. The hydrothermal activity occurred in temperature‐pressure conditions of aqueous fluid immiscibility into hypersaline brine and dilute vapor. The halite dissolution (Tm[halite]) and liquid‐vapor homogenization (Th) temperatures of the polyphase inclusions in veinlet quartz range from 270 to 472d?C and from 280 to 454d?C, respectively. The estimated salinity ranges from 36 to 47 wt% (NaCl equiv.). The apparent pressures lower than 300 bars are estimated to have been along the liquid‐vapor‐halite curve for the fluid inclusions having the Th lower than the Tm that trapped the brine saturated with halite, or at slightly higher pressure relative to liquid‐vapor‐halite curve for the fluid inclusions having the Th higher than the Tm that trapped the brine unsaturated with halite. The actual temperature and pressure during the hydrothermal activity at the Batu Hijau deposit are estimated to have been around 300d?C and 50 bars. At such temperature‐pressure conditions, the principal and initial Cu‐Fe sulfide mineral assemblages are thought to be chalcopyrite + bornite solid solution (bnss) for the chalcopyrite‐bearing assemblage, and chalcocite‐digenite solid solution and bnss for the chalcopyrite‐free assemblage.     

Deep structure of the northeastern Japan arc and its implications for crustal deformation and shallow seismic activity

Seismic tomography studies in the northeastern Japan arc have revealed the existence of an inclined sheet-like seismic low-velocity and high-attenuation zone in the mantle wedge at depths shallower than about 150 km. This sheet-like low-velocity, high-attenuation zone is oriented sub-parallel to the subducted slab, and is considered to correspond to the upwelling flow portion of the subduction-induced convection. The low-velocity, high-attenuation zone reaches the Moho immediately beneath the volcanic front (or the Ou Backbone Range) running through the middle of the arc nearly parallel to the trench axis, which suggests that the volcanic front is formed by this hot upwelling flow. Aqueous fluids supplied by the subducted slab are probably transported upward through this upwelling flow to reach shallow levels beneath the Backbone Range where they are expelled from solidified magma and migrate further upward. The existence of aqueous fluids may weaken the surrounding crustal rocks, resulting in local contractive deformation and uplift along the Backbone Range under the compressional stress field of the volcanic arc. A strain-rate distribution map generated from GPS data reveals a notable concentration of east–west contraction along the Backbone Range, consistent with this interpretation. Shallow inland earthquakes are also concentrated in the upper crust of this locally large contraction deformation zone. Based on these observations, a simple model is proposed to explain the deformation pattern of the crust and the characteristic shallow seismic activity beneath the northeastern Japan arc.