Shin’ichi Sakai

Variations of fluid pressure within the subducting oceanic crust and slow earthquakes

[1] We show fine-scale variations of seismic velocities and converted teleseismic waves that reveal the presence of zones of high-pressure fluids released by progressive metamorphic dehydration reactions in the subducting Philippine Sea plate in Tokai district, Japan. These zones have a strong correlation with the distribution of slow earthquakes, including long-term slow slip (LTSS) and low-frequency earthquakes (LFEs). Overpressured fluids in the LTSS region appear to be trapped within the oceanic crust by an impermeable cap rock in the fore-arc, and impede intraslab earthquakes therein. In contrast, fluid pressures are reduced in the LFE zone, which is deeper than the centroid of the LTSS, because there fluids are able to infiltrate into the narrow corner of the mantle wedge, leading to mantle serpentinization. The combination of fluids released from the subducting oceanic crust with heterogeneous fluid transport properties in the hanging wall generates variations of fluid pressures along the downgoing plate boundary, which in turn control the occurrence of slow earthquakes.

Multi-fault system of the 2004 Mid-Niigata Prefecture Earthquake and its aftershocks

A seismic network was deployed the day after the main shock of the 2004 Mid-Niigata Prefecture Earthquake to determine the major source faults responsible for the main shock and large aftershocks. Using the high-resolution seismic data for five days, three major source faults were identified: two parallel faults dipping steeply to the west located 5 km apart, and another dipping eastward and oriented perpendicular to the west-dipping faults. Strong lateral changes in the velocity of the source area resulted in the locations of the epicenters determined in this study being located approximately 4.3 km west-north-west of those reported by the JMA routine catalogue. The strong heterogeneity of the crust is related to the complex geological and tectonic evolution of the area and therefore the relatively large aftershocks followed around the main shock. This is considered to be responsible for the prominent aftershock activity following the 2004 Niigata event.

Urgent aftershock observation of the 2004 off the Kii Peninsula earthquake using ocean bottom seismometers

The 2004 off the Kii Peninsula earthquake occurred on September 5, 2004. Knowing the precise aftershock distribution is important for understanding the mechanism of this earthquake. However, the hypocenter of the main shock was located more than 100 km offshore from the nearest station of the land observation network. In the three days after the main shock, we started ocean bottom seismometer (OBS) observation in order to determine the precise distribution of the aftershocks. We assumed a seismic velocity structure for the hypocenter calculation, based on the results of previous seismic refraction study. The station corrections were incorporated to locate the hypocenter precisely. The hypocenters located within an area covered by five OBSs show relatively small errors. It is found that the OBS-located hypocenters are located about 5.5 km east-southeast from those by JMA and the depth range of the aftershocks is about 5–25 km just beneath the Nankai trough axis. The aftershock hypocenters can be grouped into two clusters at different depths of about 10 km and about 20 km. It is inferred that the main shock also has a depth of 5–25 km. Since this extent of the main shock was larger than one of the oceanic crust of the Philippine Sea plate, the fault plane of the main shock extended at the upper most mantle of the Philippine Sea plate. Although we cannot assign the actual fault plane of the main shock form our observation results, it is clarified that intra-plate earthquakes occurred near the trench region. Our OBS result supports that the main shock was the earthquake not at the plate boundary but within the bending Philippine Sea plate near the trough axis.

Precise aftershock distribution of the 2007 Chuetsu-oki Earthquake obtained by using an ocean bottom seismometer network

The Chuetsu-Oki Earthquake occurred on July 16, 2007. To understand the mechanism of earthquake generation, it is important to obtain a detailed seismic activity. Since the source region of the 2007 Chuetsu-oki Earthquake lies mainly offshore of Chuetsu region, a central part of Niigata Prefecture, it is difficult to estimate the geometry of faults using only the land seismic network data. A precise aftershock distribution is essential to determine the fault geometry of the mainshock. To obtain the detailed aftershock distribution of the 2007 Chuetsu-oki Earthquake, 32 Ocean Bottom Seismometers (OBSs) were deployed from July 25 to August 28 in and around the source region of the mainshock. In addition, a seismic survey using airguns and OBSs was carried out during the observation to obtain a seismic velocity structure below the observation area for precise hypocenter determination. Seven hundred and four aftershocks were recorded with high spatial resolution during the observation period using OBSs, temporally installed land seismic stations, and telemetered seismic land stations and were located using the double-difference method. Most of the aftershocks occurred in a depth range of 6–15 km, which corresponds to the 6-km/s layer. From the depth distribution of the hypocenters, the aftershocks occurred along a plane dipping to the southeast in the whole aftershock region. The dip angle of this plane is approximately 40°. This single plane with a dip to the southeast is considered to represent the fault plane of the mainshock. The regions where few aftershocks occurred are related to the asperities where large slip is estimated from the data of the mainshock. The OBS observation is indispensable to determine the precise depths of events which occur in offshore regions even close to a coast.

Imaging heterogeneous velocity structures and complex aftershock distributions in the source region of the 2007 Niigataken Chuetsu-oki Earthquake by a dense seismic observation

The velocity structure and accurate aftershock distributions in the source region of the 2007 Niigataken Chuetsu-oki Earthquake (thrust type) are obtained by inverting the arrival times from 848 aftershocks observed by a dense seismic network deployed immediately after the mainshock (8 h later). Both the detailed velocity structure and the accurate aftershock distribution show lateral heterogeneity along the fault strike. In the northeast area, aftershocks are aligned along both the NW- and SE-dipping planes. These planes are conjugate to each other. The mainshock hypocenter is located close to the bottom of an approximately 50° NW-dipping plane, which indicates that the mainshock rupture could have initiated on the NW-dipping plane. The high-Vp body beneath this aftershock alignment shows a convex upward shape. In contrast, from the center to the southwest area, most of the aftershocks are aligned along SE-dipping planes. The high-Vp body beneath this aftershock alignment shows a convex downward shape. Based on these results, we suggest that the crustal structure in the source region is divided into two segments by a boundary zone situated between the northeast and southwest areas. It should be noted that this segment boundary zone is coincident with the complex aftershock zone where numerous conjugate fault planes exist. We propose that the mainshock rupture initiated near the bottom of the NW-dipping fault plane and ran to the southwest, then transferred at the segment boundary zone which has numerous conjugate fault planes to the SE-dipping plane.

Crustal structure in the northern Fossa Magna region, central Japan, modeled from refraction/wide-angle reflection data

The northern Fossa Magna (NFM) is a back-arc rift basin filled with thick Tertiary sediments, which show strong NW-SE shortening deformation. In the NFM, there exist two major active fault systems, the Itoigawa-Shizuoka Tectonic Line active fault system (ISTL) and the Western Nagano Basin active fault system (WNB), both of which have great potentials to cause destructive earthquakes. By reanalyzing five sets of refraction/wide-angle reflection data, we successfully obtained detailed and consistent models of the crustal structure in the NFM region. It was a very effective modeling procedure to incorporate vicinal seismic reflection data and geologic information. The geometries of the active faults in the NFM region were revealed. The ISTL is east dipping, and the WNB is northwest dipping. The Tertiary sedimentary layer (<4.0 km/sec) west and adjacent to the ISTL extends to a depth of 4–5 km. The basement rocks below the Central Uplift Belt (CUB) form a wedge structure, which suggests the westward movement of the CUB basement rocks.

The 1998 Miyako fireball’s trajectory determined from shock wave records of a dense seismic array

A high velocity passage of a meteoroid through the atmosphere generates a shock wave with a conical front. When the shock front arrives at the surface, it causes high frequency ground motions that are registered on the seismograms. We can use seismological data to determine the trajectory of the meteoroid in the atmosphere. A strong shock wave from the 1998 Miyako fireball is recorded by more than 20 stations in a dense array of seismographs installed in the northeastern region of Honshu Island, Japan. We determine the velocity and the trajectory of the fireball in the upper atmosphere using the arrival times of the shock wave at the stations.

Estimation and correction for the effect of sound velocity variation on GPS/Acoustic seafloor positioning : An experiment off Hawaii Island

A GPS/Acoustic experiment on the southeastern slope of Hawaii Island presented precise seafloor positioning in the condition of large water depth (2.5—4.5 km) and large velocity variations. We estimated sound velocity variations from acoustic ranging, and found that temperature variation can well explain the velocity variation. The effect of daily variation in the sound velocity amounted to +/- 0.7 m on acoustic ranging of 4—7 km with a fixed velocity structure. CTD data observed about every 3 hours could decrease the range residuals to +/- 0.4 m. These large residuals were fairly well canceled in the positioning of the array center of three acoustic transponders. The estimated precision of the array center positioning was about 3 cm in latitude and longitude.

Precise hypocenter locations of midcrustal low-frequency earthquakes beneath Mt. Fuji, Japan

Midcrustal low-frequency earthquakes (MLFs) have been observed at seismic stations around Mt. Fuji, Japan. In September–December 2000 and April–May 2001, abnormally high numbers of MLFs occurred. We located hypocenters for the 80 MLFs during 1998–2003 by using the hypoDD earthquake location program (Waldhauser and Ellsworth, 2000). The MLF hypocenters define an ellipsoidal volume some 5 km in diameter ranging from 11 to 16 km in focal depth. This volume is centered 3 km northeast of the summit and its long axis is directed NW-SE. The direction of the axis coincides with the major axis of tectonic compression around Mt. Fuji. The center of the MLF epicenters gradually migrated upward and 2–3 km from southeast to northwest during 1998–2001. We interpret that the hypocentral migration of MLFs reflects magma movement associated with a NW-SE oriented dike beneath Mt. Fuji.

Rapid ground deformation of the Miyakejima volcano on 26-27 June 2000 detected by kinematic GPS analysis

A kinematic GPS analysis of data from the Miyakejima volcano captured a fast developing deformation event on 26–27 June 2000 in unprecedented spatial and temporal detail. Initial ground deformation toward east and upward was observed in the southeastern part of the volcano at 18:00 on 26 June 2000, almost simultaneous with earthquake swarms. Some time after 21:30 on 26 June 2000 the displacements at these sites turned from eastward to westward. Three hours later the displacement rates increased gradually at GPS sites in the western part of Miyakejima as the seismicity migrated and approached the west coast, and reached a climax with submarine eruption at 09:00 on 27 June 2000. A Genetic Algorithm was used to explore the parameter space and to find the best fitting source geometry. This analysis leads to an interpretation that the 18:00 26 June earthquake swarm was caused by a dike intrusion near the Oyama crater. Starting from 21:30 this dike deflated and a new dike intruded near the west coast. Following the propagation of this dike to the offshore, a spherical source began deflating in the southwest of Oyama crater.