Complex geometry and segmentation of the surface rupture associated with the 14 November 2001 great Kunlun earthquake, northern Tibet, China

The 14 November 2001 Kunlun, China, earthquake with a moment magnitude (M_w) 7.8 occurred along the Kusai Lake–Kunlun Pass fault of the Kunlun fault system. We document the spatial distribution and geometry of surface rupture zone produced by this earthquake, based on high-resolution satellite (Landsat ETM, ASTER, SPOT and IKONOS) images combined with field measurements. Our results show that the surface rupture zone can be divided into five segments according to the geometry of surface rupture, including the Sun Lake, Buka Daban–Hongshui River, Kusai Lake, Hubei Peak and Kunlun Pass segments from west to east. These segments, each 55 to 130 km long, are separated by step-overs. The Sun Lake segment extends about 65 km with a strike of N45° 75°W (between 90°05′E 90°50′E) along the previously unrecognized West Sun Lake fault. A gap of about 30 km long exists between the Sun Lake and Buka Daban Peak where no obvious surface ruptures can be observed either from the satellite images or field observations. The Buka Daban–Hongshui River, Kusai Lake, Hubei Peak and Kunlun Pass segments run about 365 km striking N75° 85°W along the southern slope of the Kunlun Mountains (between 91°07′E 94°58′E). This segmentation of the surface rupture is well correlated with the pattern of slip distribution measured in the field. Detailed mapping suggest that these five first-order segments can be further separated into over 20 second-order segments with a length of 10–30 km, linked by smaller scale step-overs or bends.Our result also shows that the total coseismic surface rupture length produced by the 2001 Kunlun earthquake is about 430 km (excluding the 30-km-long gap), which is the longest coseismic surface rupture for an intracontinental earthquake ever recorded.Finally, we suggest a multiple bilateral rupture propagation model that shows the rupture process of the 2001 M_w 7.8 earthquake is complex. It consists of westward and eastward rupture propagations and interaction of these bilateral rupture processes.     

Fault system of the 2004 Mid Niigata Prefecture Earthquake and its aftershocks

Data recorded by a seismic network deployed the day after the 2004 Mid Niigata Prefecture Earthquake (M6.8) in central Japan are used to determine the major source faults responsible for the mainshock and major aftershocks. Using this high-resolution seismic data, three major source faults are identified: two parallel faults dipping steeply to the west located 5 km apart, and the other dipping eastward and oriented perpendicular to the west-dipping faults. The analysis also reveals that the lateral variation in seismic velocity observed at the surface extends to a depth of 15 km, encompassing the source area of the mainshock. This strong heterogeneity of the crust, related to the complex geological and tectonic evolution of the area, is considered to be responsible for the prominent aftershock activity following the 2004 Niigata event.     

Tidal triggering of earthquakes in the subducting Philippine Sea plate beneath the locked zone of the plate interface in the Tokai region, Japan

We found a characteristic space–time pattern of the tidal triggering effect on earthquake occurrence in the subducting Philippine Sea plate beneath the locked zone of the plate interface in the Tokai region, central Japan, where a large interplate earthquake may be impending. We measured the correlation between the Earth tide and earthquake occurrence using microearthquakes that took place in the Philippine Sea plate for about two decades. For each event, we assigned the tidal phase angle at the origin time by theoretically calculating the tidal shear stress on the fault plane. Based on the distribution of the tidal phase angles, we statistically tested whether they concentrate near some particular angle or not by using Schuster's test. In this test, the result is evaluated by p-value, which represents the significance level to reject the null hypothesis that earthquakes occur randomly irrespective of the tidal phase angle. As a result of analysis, no correlation was found for the data set including all the earthquakes. However, we found a systematic pattern in the temporal variation of the tidal effect; the p-value significantly decreased preceding the occurrence of M ≥ 4.5 earthquakes, and it recovered a high level afterwards. We note that those M ≥ 4.5 earthquakes were considerably larger than the normal background seismicity in the study area. The frequency distribution of tidal phase angles in the pre-event period exhibited a peak at the phase angle where the tidal shear stress is at its maximum to accelerate the fault slip. This indicates that the observed small p-value is a physical consequence of the tidal effect. We also found a distinctive feature in the spatial distribution of p-values. The small p-values appeared just beneath the strongly coupled portion of the plate interface, as inferred from the seismicity rate change in the past few years.     

Fracture Toughness Evaluation Based on Tension-softening Model and its Application to Hydraulic Fracturing

This paper discusses the applicability of the tension-softening model in the determination of the fracture toughness of rocks, where the fracture toughness evaluated based on the tension-softening model is compared with the crack growth resistance deduced from laboratory-scale hydraulic fracturing tests. It is generally accepted that the fracture process is dominated by the growth of a fracture process zone for most types of rocks. In this study, the J-integral based technique is employed to determine the fracture toughness of Iidate granite on the basis of the tension-softening model, where compact tension specimens of different dimensions were tested in order to examine the specimen size effect on the measured fracture toughness. It was shown that the tension-softening relation deduced from the J-integral based technique allowed us to determine the specimen size independent fracture toughness K_c of Iidate granite. Laboratory-scale hydraulic fracturing tests were performed on cubic specimens (up to a 10 m sized specimen), where cyclic pressurization was conducted using a rubber-made straddle packer to observe the extent of the hydraulically induced crack. The experimental results of pressure and crack length were then used to construct the crack growth resistance curve based on the stress intensity factor K. The crack growth resistance obtained from the hydraulic fracturing tests was observed to initially increase and then level off, giving a constant K value for a long crack extension stage. The plateau K value in the crack growth resistance curve was found to be in reasonable agreement with the fracture toughness K_c deduced from the tension-softening relation. It was demonstrated that the tension-softening model provides a useful tool to determine the appropriate fracture toughness of rocks, which may be applicable for the analysis of the process of large-scale crack extension in rock masses.     

Three-dimensional Mapping of Magma Source and Transport Regions from Seismic Data: The Mantle Wedge beneath Northeastern Japan

—We investigate the distribution of partial melt in island arc using the seismic velocity structure of the mantle wedge beneath northeastern Japan. The comparison of the seismic tomography with laboratory velocity data on a partially-molten mantle rock yields estimates of melting zones in three dimensions. We employ experimental data on the degree of partial melt in hydrous peridotite to give constraints on the melt fraction and temperature. Melting and magma-rich zones derived from the velocity structure coincide with observed low Q zones. The results of the three-dimensional mapping indicate that the source of magma in island arc is diapir-like melting patches localized within the low velocity zones of the mantle wedge. Extensive volcanic activity along the volcanic front is due to the presence of vast magma-rich zones just beneath the Moho. Those melting zones in the uppermost mantle may, in turn, cause melting of lower crustal materials and produce felsic magma. Melt appears to stay at and beneath the Moho, where crystallization fractionation may proceed. Melt exists at greater depths in the back-arc region, which may correlate with across-arc variations of chemical compositions of the volcanic rocks observed in northeastern Japan. We suggest that magma migration in the ductile lower crust may cause low-frequency microearthquakes, and magma penetration into the brittle upper crust may produce mid-crustal S-wave reflectors.     

The Holocene development of Kushu Lake on Rebun Island in Hokkaido, Japan

Diatom assemblages and sulfur content in sediments were analyzed to clarify changes in the sedimentary environment of Kushu Lake, a coastal lake on Rebun Island in Hokkaido, Japan. Salinity variations were assessed by means of a diatom-based index of paleosalinity and the sedimentary sulfur content. This paper discusses the Holocene development of the lake, in relation to Holocene relative sea-level change. For paleoenvironmental interpretation of the lake development, the rationale of the threshold method (Anundsen et al., 1994) was applied.At ca. 8000 yr BP, a coastal embayment (paleo-Kushu Bay) resulted from marine ingression. The threshold elevation at the mouth of the paleo-Kushu Bay kept pace with the rising sea-level, resulting in its enclosure at the culmination of Holocene marine transgression (ca. 6500–5000 yr BP). From predicted relative sea-level at ca. 6000 yr BP for Rebun Island (Nakada et al., 1991), the threshold may have been at least above –3 to –5 m altitude. A freshwater lake environment with strongly anoxic bottom conditions may have occurred from ca. 5500 to 5100 yr BP. After an important episode of marine ingression, the lake was isolated completely from the open sea at ca. 4900 yr BP. The diatom record suggests that the maximum lacustrine extent occurred at ca. 4900–3100 yr BP. Thereafter, water depth decreased at the lake margins.In Kushu Lake, the threshold elevation, due to a build-up of a coastal barrier, prevents us from determining the amplitude of sea-level changes, even though the age of isolation contacts corresponds to periods of regression and climatic deterioration. In spite of isostatic subsidence, the effective protection provided by the well-developed barrier did not allow registration of any relative sea-level fluctuations since its isolation.     

Crustal movement along a collision boundary of plates (case of eastern Taiwan)

Recent crustal movements in eastern Taiwan were investigated based on geodetic and seismological data as a case study of the tectonics along a colliding plate boundary. Secular horizontal strain along the Longitudinal Valley, a probable arc-continent collision boundary between the Philippine Sea and Eurasion plates, is compressive with a large rate of 2 microstrains/year, and the maximum compression axis trends in NW-SE direction. This means that nearly two-thirds of the convergence rate between the two plates is consumed in the narrow zone along the valley and confirms the view that the valley is the plate boundary between the two. The amount of uplift in the Central Range, deduced from trigonometric leveling, reaches more than one meter during the past 60 years. Such active orogeny is caused by the high strain rate in this area. Horizontal displacement of the Coastal Range, a leading edge of the Philippine Sea plate, has a left-lateral component, suggesting that the boundary is a region of oblique collision and the range moves toward the north, producing subduction of the Philippine Sea plate off the Coast of northeastern Taiwan.An extensional strain pattern is found northeastern Taiwan near Hualien, in northern Taiwan near Ilan, and in southeastern Taiwan south Taitung. In northeastern Taiwan, the direction of the maximum extension axis is nearly parallel to the isodipth contours of deep-focus earthquakes that occurred under the southern part of the Ryukyu Trench. From a model calculation using a finite-element method, this pattern of strain is interpreted as an edge effect of the collision-subduction junction in northeastern Taiwan. The direction of the maximum extension in northern Taiwan forms high angles with the direction of the tectonic line near Ilan stretching from the Okinawa Trough. This fact is harmonious with the view that the Okinawa Trough is spreading. Extensional strain observed in southeastern Taiwan may be related to the eastward subduction since the area is situated on a transition zone, from collision to subduction.     

Wind fluctuations near a cold vortex-tropopause funnel system observed by the MU radar

Vertical and temporal variations of three-dimensional wind velocity associated with an upper-tropospheric cold vortex-tropopause funnel system were observed by an MST radar in Japan (the MU radar). Marked changes of vertical velocity and horizontal wind direction were found between the inside and outside of the cold vortex. The vertical velocity activity outside the vortex was asymmetric; it was most active in a sector before the vortex. Unsaturated internal gravity waves in their generation stage contribute predominantly to the vertical velocity activity, suggesting that tropospheric occluded cyclones may be a possible source of middle-atmospheric gravity waves through the geostrophic adjustment process.