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.     

胶东海阳所超镁铁岩中强亲铁元素的地球化学特征及其构造环境探讨

强亲铁元素与亲石元素具有不同的地球化学行为，因此能够从不同的角度为造山带中超镁铁岩的成因及演化提供重要信息。位于苏鲁造山带东北端的胶东海阳所超镁铁岩主要由橄榄岩和辉石岩组成，它们常以团块状赋存于花岗质片麻岩中。虽然前人对这些超镁铁岩已经开展大量岩石学研究，但关于其成因及构造属性仍存在较大争议。本文开展了海阳所超镁铁岩的全岩主微量元素、强亲铁元素及Re-Os同位素的分析工作，结果显示蛇纹石化橄榄岩具有较高的MgO和Fe₂O₃^T含量，较低的Al₂O₃、TiO₂和CaO含量，明显富集流体迁移元素（U、Pb），亏损高场强元素（Zr、Hf），强亲铁元素没有发生明显分异，但Ru显示正异常，表明海阳所蛇纹石化橄榄岩是经历了低-中等程度部分熔融及熔/流体交代作用影响的残余地幔橄榄岩。海阳所辉石岩的主量元素表现出明显的结晶分异特征，稀土元素较原始地幔富集，铂族元素（PGEs）含量较低且发生了明显的分异，表明辉石岩的地幔源区经历过高程度的部分熔融和硫化物的分离。海阳所蛇纹石化橄榄岩的Os同位素地球化学特征表现出大洋亲和性，与辉石岩不具有熔体-残留体的关系。由于该地区发育较深层次的韧性剪切带，蛇纹石化橄榄岩中的橄榄石与辉石表现出韧性变形的特征，同时有辉石岩侵入到橄榄岩的现象，表明该地区的蛇纹石化地幔橄榄岩与辉石岩既不同时，也不同源，因此，暗示了该套岩石组合可能形成于大洋核杂岩（OCC）与洋脊型蛇绿岩（MOR）堆晶岩交互发育环境。     

Possible creep-related tilt precursors obtained in the Central Apennines (Italy) and in the Southern Caucasus (Georgia)

Daily averaged tilt component data from two sites of the Central Apennines (Italy) and of the Southern Caucasus (Georgia), respectively, revealed intermediate-term tilts as possible precursors to earthquakes (M=3.0÷4.7) which occurred in the above-mentioned seismic areas within a distance of 50 km from the sites. A good temporal correlation as well as a fair spatial correspondence between these residual tilts (with amplitude and duration of some microradians and months, respectively) and main shocks were pointed out, by removing both secular trends and seasonal thermoelastic effects from the raw tilts. An attempt was made to justify the above-mentioned results, based on the assumption that the observed intermediate-term preseismic tilts are the manifestation of aseismic creep episodes of comparable duration in the fault materials of thrust faults close to the tilt sites. The mechanism refers to a strain field slowly propagating from the preparation focal area to the tilt site, through crustal blocks separated by weak transition zones. This propagation is thought to be the cause of the local aseismic fault slip recorded by the tiltmeters. Previously, both discrete structures and strain propagation effects were revealed in the Central Apennines and are thought also to exist in the Southern Caucasus. As in the past, the rheological properties of fault materials are revealed as viscoelastic ones. In fact, creep equations obtained by applying several viscoelastic models on our data, proved to fit quite well some of the observed tilt precursors, producing viscosity and rigidity values very similar to those reported in literature.Professor Petr Viktorovich Manjgaladze died during the writing of this paper     

Seasonal variation of atmospheric dimethylsulfide at Amsterdam Island in the southern Indian Ocean

Daily measurements of atmospheric concentrations of dimethylsulfide (DMS) were carried out for two years in a marine site at remote area: the Amsterdam Island (37°50S–77°31E) located in the southern Indian Ocean. DMS concentrations were also measured in seawater. A seasonal variation is observed for both DMS in the atmosphere and in the sea-surface. The monthly averages of DMS concentrations in the surface coastal seawater and in the atmosphere ranged, respectively, from 0.3 to 2.0 nmol l^-1 and from 1.4 to 11.3 nmol m^-3 (34 to 274 pptv), with the highest values in summer. The monthly variation of sea-to-air flux of DMS from the southern Indian Ocean ranges from 0.7 to 4.4 mol m^-2 d^-1. A factor of 2.3 is observed between summer and winter with mean DMS fluxes of 3.0 and 1.3 mol m^-2 d^-1, respectively.     

Fault zone geometry of a mature active normal fault: A potential high permeability channel (Pirgaki fault, Corinth rift, Greece)

We present results from petrophysical analysis of a normal fault zone with the aim of defining the flow pathways and their behavior during seismic and interseismic periods. Data are obtained on porosity geometry, strain structure and mineralogy of different domains of a normal fault zone in the Corinth rift. Data point out a close relationship between mineralogy of the clayey minerals, porosity network and strain structures and allow definition of a macroscopic anisotropy of the flow parameters with a strong control by microscopic ultracataclasite structures. The Pirgaki fault zone, developed within pelagic limestone, has a sharp asymmetric porosity profile, with a high porosity volume in the fault core and in the damage zone of the hanging wall. From porosity volumes and threshold measurements, a matrix permeability variation of 6 orders of magnitude could be expected between the protolith and the fault core. Modifications of this pathway during seismic and interseismic phases are depicted. Healing of cracks formed during seismic slip events occurred in the fault core zone and the porous network in the damage zone is sealed in a second step. The lens geometry of the fault core zone is associated with dissolution surfaces and open conduits where dissolved matter could move out of the fault core zone. These elementary processes are developed in particular along Riedel's structures and depend on the orientation of the strain surfaces relative to the local stress and depend also on the roughness of each surface type. P-surfaces are smooth and control shearing process. R-surfaces are rough and present two wavelengths of roughness. The long one controls localization of dissolution surfaces and conduits; the short one is characteristic of dissolution surfaces. The dissolved matter can precipitate in the open structures of the hanging wall damage zone, decreasing the connectivity of the macroscopic conduit developed within this part of the fault zone.     

Fracture analysis in the south-western Corinth rift (Greece) and implications on fault hydraulic behavior

This paper reviews the data concerning the fracture network and the hydraulic characteristics of faults in an active zone of the Gulf of Corinth. Pressure gap measured through fault planes shows that in this area the active normal faults (Aigion, Helike) act, at least temporarily and locally, as transversal seal. The analysis of the carbonate cements in the fractures on both the hangingwall and the footwall of the faults also suggests that they have acted as local seals during the whole fault zone evolution. However, the pressure and the characteristics of the water samples measured in the wells indicate that meteoric water circulates from the highest part of the relief to the coast, which means it goes through the fault zones. Field quantitative analysis and core studies from the AIG-10 well have been performed to define both regional and fault-related fracture networks. Then laboratory thin section observations have been done to recognize the different fault rocks characterizing the fault zone components. These two kinds of approach give information on the permeability characteristics of the fault zone. To synthesize the data, a schematic conceptual 3D fluid flow modeling has been performed taking into account fault zone permeability architecture, sedimentation, fluid flow, fault vertical offset and meteoric water influx, as well as compaction water flow. This modeling allows us to fit all the data with a model where the fault segments act as a seal whereas the relays between these segments allow for the regional flow from the Peloponnese topographic highs to the coast.     

Fault slip analysis of Quaternary faults in southeastern Korea

The Quaternary stress field has been reconstructed for southeast Korea using sets of fault data. The subhorizontal direction of the maximum principal stress (σ₁) trended ENE and the direction of the minimum principal stress (σ₃) was nearly vertical. The stress ratio (Φ = (σ₂ − σ₃) / (σ₁ − σ₃)) value was 0.65. Two possible interpretations for the stress field can be made in the framework of eastern Asian tectonics; (1) The σ_Hmax trajectory for southeast Korea fits well with the fan-shaped radial pattern of maximum principal stress induced by the India–Eurasia collision. Thus, we suggest that the main source for this recent stress field in southeast Korea is related to the remote India–Eurasia continental collision. (2) The stress field in Korea shows a pattern similar to that in southwestern Japan. The origin for the E–W trending σ_Hmax in Japan is known to be related to the mantle upwelling in the East China Sea. Thus, it is possible that Quaternary stress field in Korea has evolved synchronously with that in Japan. We suggest further studies (GPS and in situ stress measurement) to test these hypotheses.     

A laboratory experiment to monitor the contact state of a fault by transmission waves

We performed a series of laboratory experiments in which elastic waves were transmitted across a simulated fault. Two types of experiments were carried out: (1) Normal Stress Holding Test (NSHT): normal stress was kept constant for about 3 h without shear stress and transmission waves were observed. (2) Shear Stress Increasing Test (SSIT): shear stress was gradually increased until a stick-slip event occurred. Transmission waves were continuously observed throughout the process of stress accumulation. We focused on the change in transmission waves during the application of shear stress and especially during precursory slips.It was found in NSHT that the amplitude of transmission waves linearly increased with the logarithm of stationary contact time. The increase amounted to a few percent after about 3 h. Creep at asperity contacts is responsible for this phenomenon. From a theoretical consideration, it was concluded that the real contact area increased with the logarithm of stationary contact time.We observed in SSIT a significant increase in wave amplitude with shear stress application. This phenomenon cannot be attributed to the time effect observed in NSHT. Instead, it can be explained by the mechanism of “junction growth” proposed by Tabor. Junction growth yields an increase in real contact area. It is required for junction growth to occur that the material in contact is already plastic under a purely normal loading condition. A computer simulation confirmed that this requirement was satisfied in our experiments. We also found that the rate at which the amplitude increased was slightly reduced prior to a stick-slip event. The onset time of the reduction well coincides with the onset of precursory slip. The cause of the reduction is attributed to the reset of stationary contact time due to displacement. This interpretation is supported by the result of NSHT. Taking the time of stationary contact in SSIT into account, we may expect the change in wave amplitude to be, at most, only a few percent. The observed slight reduction in increasing rate is, in this sense, reasonable. The static stiffness of the fault also decreases with precursory slip. It was also found that low frequency waves are a better indicator of precursory slip than high frequency waves. This might suggest that low frequency waves with longer wavelength are a better indicator of average behavior of faults. The problem, however, merits a further investigation. The shifts in phase were also found to be a good indicator of the change in contact state of the fault. The changes in both amplitude and phase of transmission waves are unifyingly understood through the theory of transmission coefficient presented by Pyrak-Nolte et al. Rough surfaces have a tendency to give larger stick-slips than smooth surfaces. The amount of precursory slip is larger for rough surfaces than for smooth surfaces. Although it was confirmed by a computer simulation that rough surfaces have larger contact diameters than smooth surfaces, the rigorous relationship between the surface roughness (contact diameter) and the amount of precursory slips was not established.     

Deep structure of the Nojima Fault, southwest Japan, estimated from borehole observations of fault-zone trapped waves

To estimate the deep structure of the southern part of the Nojima Fault, southwest Japan without the influence of near-surface structures, we analyzed the Love-wave-type fault-zone trapped waves (LTWs) recorded by a borehole seismometer at 1800 m depth. We examined the polarization, dispersion, and dominant frequency of the wavetrain following the direct S-wave in each seismogram to identify the LTW. We selected eight candidates for typical LTWs from 462 records. Because the duration of the LTW increases with hypocentral distance, we infer that the low velocity fault-zone of the Nojima Fault continues towards the seismogenic depth. In addition, since the duration of the LTW increases nonlinearly with hypocentral distance, we infer that the S-wave velocity of the fault-zone increases with depth. The location of events showing the LTW indicates that the fault-zone dips to the southeast at 75° and continues to a depth of approximately 10 km. We assumed a uniform low-velocity waveguide to estimate the average structure of the fault-zone. We estimated the average width, S-wave velocity, and Q_s of the fault-zone by comparing an analytical solution of the LTW with measured data. The average width, S-wave velocity, and Q_s of the fault-zone are 150 to 290 m, 2.5 to 3.2 km/s, and 40 to 90, respectively. Hence the fault-zone structure with a larger width and smaller velocity reduction than the fault-zone model estimated by previous surface observation is more suitable to represent the average fault-zone structure of the Nojima fault. The present study also indicated that the shallow layers and/or a shallow fault-zone structure drastically changes the characteristics of the LTW recorded at the surface, and therefore cause a discrepancy in the fault-zone model between the borehole observation and surface observation.