Application of a modified pattern informatics method to forecasting the locations of future large earthquakes in the central Japan

We propose a modification of the Pattern Informatics (PI) method that has been developed for forecasting the locations of future large earthquakes. This forecast is based on analyzing the space–time patterns of past earthquakes to find possible locations where future large earthquakes are expected to occur. A characteristic of our modification is that the effect of errors in the locations of past earthquakes on the output forecast is reduced. We apply the modified and original methods to seismicity in the central part of Japan and compared the forecast performances. We also invoke the Relative Intensity (RI) of seismic activity and randomized catalogs to constitute null hypotheses. We do statistical tests using the Molchan and Relative Operating Characteristic (ROC) diagrams and the log-likelihoods and show that the forecast for using the modified PI method is generally better than the competing original-PI forecast and the forecasts from the null hypotheses. Using the bootstrap technique with Monte-Carlo simulations, we further confirm that earthquake sequences simulated based on the modified-PI forecast can be statistically the same as the real earthquake sequence so that the forecast is acceptable. The main and innovative science in this paper is the modification of the PI method and the demonstration of its applicability, showing a considerable promise as an intermediate-term earthquake forecasting tool.     

Background seismicity in the Central Apennines of Italy: The Abruzzo region case study

We investigate background seismic activity of the Abruzzo region, a 5000 km² area located within the Central Apennines of Italy, where in the past 600 years at least 5 large earthquakes (I = XI–X) have occurred.Between April 2003 and September 2004, a dense temporary seismic network composed of 30 digital three-component seismic stations recorded 850 earthquakes with 0.9 < M_L < 3.7. We present earthquake locations and focal mechanisms obtained by standard procedures and an optimized velocity model computed with a search technique based on genetic algorithms.The seismicity occurs at a low and constant rate of  2.6 e− 04 events/daykm² and is sparsely distributed within the first 15 km of the crust. Minor increases in the seismicity rate are related to the occurrence of small and localised seismic sequences that occur at the tip of major active normal faults along secondary structures.We observe that during the 16 months of study period, the Fucino fault system responsible for the 1915 Fucino earthquake (M_S = 7.0), and the major normal faults of the area, did not produce significant seismic activity.Fault plane solutions evaluated using P-wave polarity data show the predominance of normal faulting mechanisms ( 55%) with NE-trending direction of extension coherent with the regional stress field active in this sector of the Apennines. Around 27% of the focal solutions have pure strike–slip mechanisms and the rest shows transtensional faulting mechanisms that mainly characterise the kinematics of the secondary structures activated by the small sequences.We hypothesize that the largest known NW-trending normal faults are presently locked and we propose that in the case of activation, the secondary structures located at their tips may act as transfer faults accommodating a minor part of the extensional deformation with strike–slip motion.     

Variations of Seismicity in the Avachinsky Gulf (Kamchatka, Russia)

Variations of seismic mode in the region of the Avachinsky Gulf (Kamchatka, Russia) are considered. Observed anomalies (seismic quiescence, the ring seismicity, reduction of the slope of the earthquake recurrence diagram) provide a basis to consider this region as a place of strong earthquake preparation. The Kamchatka regional catalogues of earthquakes between 1962–1995 were used in the analysis. A reduced seismicity rate is observed during 10 years in an area of 150 km × 60 km in size. During the last five years, in the vicinity of the area considered, earthquakes with M > 5 occurred three times more often than the average over thirty years. It is interpreted as ring seismicity. The block of 220 km × 220~km in size, including the quiescence zone, is characterized by a continuous decrease of the recurrence diagram slope, which has reached a minimum value for the last 33 years in this region.     

Asperity distribution of the 1964 Great Alaska earthquake and its relation to subsequent seismicity in the region

The 1964 Alaska earthquake was the second largest seismic events in the 20th century. The aim of this work is the use of surface deformation data to determine asperity and slip distributions on the fault plane of the Alaska earthquake: these distributions are calculated by a Monte Carlo method. To this aim, we decompose the fault plane in a large number of small square asperity units with a side of 25 km; this allows us to obtain plane surfaces with an irregular shape. In the first stage, each asperity unit is allowed to slip a constant amount or not to slip at all, providing the geometry of the dislocation surface that best reproduces the observed displacements. To this purpose, a large number of slip distributions have been tried by the use of the Monte Carlo method. The slip amplitude is the same for all the asperities and is equal to the average fault slip inferred from the seismic moment. In the second stage, we evaluate the slip distribution in the dislocation area determined by the Monte Carlo inversion: in this case, we allow unit cells to undergo different values of slip in order to refine the initial dislocation model. The results confirm the previous finding that the slip distribution of the great Alaska earthquake was essentially made of two dislocation areas with a higher slip, the Prince William Sound and the Kodiak asperities. Analysis of the post-1964 seismicity in the rupture region shows a strong correlation between the larger earthquakes (M_w≥6) and the distribution of locked asperities following the 1964 event, which can be considered as an independent test of the validity of the model. We do not find slip values higher than 25 m for any of the patches, and we determine two separate high-slip zones: one correspondent to the Prince William Sound asperity, and one (18 m slip) to the Kodiak asperity. The slip distribution connected with the 1964 shock appears to be consistent with the following seismicity in the region.     

Paleoseismological evidence for non-characteristic behavior of surface rupture associated with the 2004 Mid-Niigata Prefecture earthquake, central Japan

The 2004 Mid-Niigata Prefecture earthquake sequence (mainshock magnitude, M_JMA 6.8), which occurred in an active fold-and-thrust belt in northern central Japan, generated a small thrust surface rupture (< 20 cm of vertical displacement) along a previously unmapped northern extension of the active Muikamachi–Bonchi–Seien fault zone, on the eastern margin of the epicentral region. To better understand past seismic behavior of the rupture, we conducted a paleoseismic trenching study across the 10-cm-high west-side-up surface rupture at the foot of a pre-existing 1.8-m-high east-facing scarp, which probably resulted from past earthquake(s). A well-defined west-dipping thrust fault zone accompanied by drag folding and displacing the upper Pliocene to lower Pleistocene strata and the unconformably overlying upper Pleistocene (?) to Holocene strata was exposed. The principal fault zone is connected directly to the 2004 surface rupture. From the deformational characteristics of the strata and radiocarbon dating, we inferred that two large paleoseismic events occurred during the past 9000 years prior to the 2004 event. These two pre-2004 events have a nearly identical fault slip (at minimum, 1.5 m), which is ≥ 15 times that of the 2004 event (∼ 10 cm). These paleoseismic data, coupled with the geological and geomorphological features, suggest that the 2004 event represented non-characteristic behavior of the fault, which can potentially generate a more destructive earthquake accompanied by meter-scale surface displacement. This study provides insight into the interpretation of past faulting events and increases our understanding of rupture behavior.     

The role of Pan-African structures in intraplate seismicity near the termination of the Romanche fracture zone, West Africa

The distribution of epicenters of both historic earthquakes and recent seismic events in southeastern Ghana, compiled from local and teleseismic networks, show strong correlation with the Pan-African structures onshore and indicate an alignment with disruptions on seismic sections offshore. The seismic reflection sections reveal basement structures of the external zone of the Pan-African Dahomeyide orogen and these structures can be traced to offsets of shelf strata and seabottom reflectors, providing direct evidence, for the first time, for neotectonic activity that may be responsible for seismicity in the area. The deep structure of the external zone consists of moderately-dipping reflectors inferred to represent high-strain zones in the variably deformed margin of the West African craton. Taken together, the available data suggest that active tectonics in this intraplate environment may involve inversion of the Pan-African thrust structures but that this activity is apparently not related to reactivation of the nearby Romanche Fracture Zone.     

Seismic velocity discontinuities in the crust and uppermost mantle beneath the Kanto district, central Japan, identified from receiver function imaging and repeating earthquake activity

We constructed vertical cross-sections of depth-converted receiver function images to estimate the seismic velocity structure of the crust and uppermost mantle beneath the Kanto district, central Japan. Repeating earthquake data for the plate boundary were also used to estimate geometries of the subducting Philippine Sea plate and the subducting Pacific plate. As a result, we present images of some major seismic discontinuities. The upper boundary of the Pacific plate dips to the northwest in northern Kanto and to the west–southwest in southern Kanto with some undulations. On the other hand, the upper boundary of the Philippine Sea plate as a whole dips to the northwest. However, it is concave to the northeast in the southern Boso peninsula. We suggest that the low-velocity mantle wedge may be indicated on the top of both subducting plates. Plate thickness gradually decreases to the northeast. The northeastern end of the Philippine Sea plate is interpreted to be at depths of 45–90 km. The Moho discontinuity in the overriding plate is deeper than 25 km in the northern Kanto. It contacts the subducting Philippine Sea plate in the southwestern part near 35.8°N.     

Regional variation in exhumation and strain rate of the high-pressure Sambagawa metamorphic rocks in central Shikoku, south-west Japan

Regional variation in the P–T path of the Sambagawa metamorphic rocks, central Shikoku, Japan has been inferred from compositional zoning of metamorphic amphibole. Rocks constituting the northern part (Saruta River area) exhibit a hairpin type P–T path, where winchite/actinolite grew at the prograde stage, the peak metamorphism was recorded by the growth of barroisite to hornblende and sodic amphibole to winchite/actinolite grew at the retrograde stage. In the southern part (Asemi River area), rocks exhibit a clockwise type P–T path, where barroisite to hornblende core is rimmed by winchite to actinolite. The difference in P–T path could suggest a faster exhumation rate (i.e. more rapid decompression) in the southern than in the northern part. On the other hand, physical conditions of deformation during the exhumation stage have been independently inferred from microstructures in deformed quartz. Recrystallized quartz grains in rocks from the low‐grade (chlorite and garnet) zones are much more stretched in the southern part (aspect ratio ≥ 4.0) than in the northern part (aspect ratio< 4.0), indicating a higher strain rate in the former than in the latter. These facts may indicate that the exhumation and strain rates are correlated (i.e. the exhumation rate increases with increasing the strain rate). The difference in the exhumation rate inferred from amphibole zoning between the northern and southern parts could be explained by an extensional model involving normal faulting, where the lower plate can be exhumed faster than the upper plate due to the displacement along the fault. Furthermore, the model may explain the positive correlation between the exhumation and strain rates, because the lower plate tended to support more stress than the upper plate.     

Robustness of the holistic seismic risk evaluation in urban centers using the USRi

The Urban Seismic Risk index (USRi) published in a previous article (Carreño et al., Nat Hazards 40:137–172, 2007) is a composite indicator that measures risk from an integrated perspective and guides decision-making for identifying the main interdisciplinary factors of vulnerability to be reduced or intervened. The first step of the method is the evaluation of the potential physical damage (hard approach) as a result of the convolution of the seismic hazard with the physical vulnerability of buildings and infrastructure. Subsequently, a set of social context conditions that aggravate the physical effects is also considered (soft approach). According to this procedure, the physical risk index is evaluated for each unit of analysis from existing loss scenarios, whereas the total risk index is obtained by multiplying the former index by an impact factor using an aggravating coefficient, based on variables associated with the socio-economic conditions of each unit of analysis. The USRi has been developed using the underlying holistic and multi-hazard approach of the Urban Risk Index framework proposed for the evaluation of disaster risk in different megacities worldwide. This article presents the sensitivity analysis of the index to different parameters such as input data, weights and transformation functions used for the scaling or normalization of variables. This analysis has been performed using the Monte Carlo simulation to validate the robustness of this composite indicator, understanding as robustness how the cities maintain the ranking as well as predefined risk level ranges, when compared with the deterministic results of risk. Results are shown for different cities of the world.     

Finite element modelling of elastic intraplate stresses due to heterogeneities in crustal density and mechanical properties for the Jabalpur earthquake region,central India

Deep lower crustal intraplate earthquakes are infrequent and the mechanism of their occurrence is not well understood. The Narmada-Son-lineament region in central India has experienced two such events, the 1938 Satpura earthquake and the 1997 Jabalpur earthquake, having a focal depth of more than 35 km. We have estimated elastic stresses due to the crustal density and mechanical properties heterogeneities along the Hirapur-Mandla profile passing through the Jabalpur earthquake region to analyse conditions suitable for the concentration of shear stresses in the hypocentral region of this earthquake. Elastic stresses have been computed by a finite element method for a range of material parameters. The results indicate that the shear stresses generated by the density heterogeneities alone are not able to locally enhance the stress concentration in the hypocentral region. The role of mechanical properties of various crustal layers is important in achieving this localization of stresses. Among a range of material parameters analysed, the model with a mechanically strong lower crust overlying a relatively weak sub-Moho layer is able to enhance the stress concentration in the hypocentral region, implying a weaker mantle in comparison to the lower crust for this region of central India.     

Seismic investigation of near-surface geological structure in the Paducah, Kentucky, area: application to earthquake hazard evaluation

Seismics method were used to evaluate shallow geological conditions at 33 sites in the vicinity of Paducah, Kentucky. A combined set of P- and S-wave seismic refraction and reflection soundings were used, in addition to local borehole information, to produce structure maps of (1) a shallow (< 30 m deep) horizon believed to represent an unconformity surface at the top of the Eocene, and (2) the Paleozoic bedrock surface (< 85 to > 160 m deep). Shear-wave velocity contrasts across the shallow unconformity were generally 2-to-1 while the contrast at the top of the Paleozoic bedrock exceeds 5-to-1. These seismic boundaries have been determined to be very important in modelling and interpreting earthquake ground motion amplification in the Paducah area. The quality and accuracy of the data, and the cost effective nature of the methods, suggest that other communities in areas at risk to damage from seismic activity, with foundation conditions comparable to Paducah, might benefit from similar characterization in order (1) to identify seismically hazardous, near-surface, geological conditions, and (2) to develop geological models that could be used in computer simulations of site response.     

Pull-down basin in the central part of Japan due to subduction-induced mantle flow

The driving force for the basin subsiding against isostatic balance in and around Lake Biwa in the Kinki district, Japan is discussed. The lake region is characterized by strong negative Bouguer anomalies, especially by a steep horizontal gradient zone of gravity anomaly running along the western margin of the lake. The large negative anomaly (>50 mgal) cannot be explained by low-density sediments beneath it. A down-warping structure extending to the Moho depth should be taken into account. This conjecture has been strongly supported by a short-period receiver function imaging, which shows a clear offset of about 8 km for the Moho discontinuity under the steep gravity gradient zone.A question arises as to what is the driving force to create such a large down-warping structure. We consider that the subduction of the shallow-dipping slab under the region (Philippine Sea Slab) may cause crustal deformation by dragging the viscous mantle downward. In order to verify this model, we simulated the induced mantle flow due to the subduction of the Philippine Sea Slab and the pressure distribution on the crust–mantle boundary. This numerical experiment showed that the induced flow makes a strong negative pressure zone under the lake region if the slab has a vertical offset along the direction of subduction. This offset of the slab is consistent with plate models deduced from hypocentral distributions and Sp phases of the deep-focus earthquakes.     

Propagation of seismic slip from brittle to ductile crust: Evidence from pseudotachylyte of the Woodroffe thrust, central Australia

The Woodroffe thrust, central Australia, is a > 1.5-km-wide mylonitized shear zone marked by large volumes of mm- to cm-scale pseudotachylyte veins. The pseudotachylytes display typical melt-origin features, including rounded and embayed clasts, spherulitic and dentritic microlites, and flow structures within a fine-grained matrix. Three types of pseudotachylyte are identified on the basis of deformation texture, vein morphology, and host-rock lithology: cataclasite-related (C-Pt), mylonite-related (M-Pt), and ultramylonite-related (Um-Pt). The M-Pt and Um-Pt veins intrude into mylonite and ultramylonite and are themselves overprinted by subsequent mylonitization. These pseudotachylytes contain an internal foliation defined by flattened porphyroclasts and layering of the fine-grained vein matrix, and the foliation is generally oriented parallel to foliation in the surrounding mylonite and ultramylonite. These observations constrain the timing and environment of M-Pt and Um-Pt pseudotachylyte formation to a protracted period of deformation and mylonitization within the ductile regime of the crust. The M-Pt and Um-Pt veins, as well as the host mylonite, are overprinted by cataclasis and multiple generations of late-stage C-Pt veins that were generated in the brittle-dominated regime of the upper crust during uplift and exhumation of the shear zone.The coexistence of multiple generations of voluminous C-Pt, M-Pt, and Um-Pt veins indicates that the pseudotachylyte veins represent a large number of large earthquakes and accompanying seismic slip over an extended period of seismicity on the Woodroffe thrust. The timing and distribution of pseudotachylyte indicate that the earthquakes nucleated at the base of the brittle-dominated seismogenic zone and propagated down through the brittle–ductile transition into the ductile-dominated regime of the crust.     

Shallow landslides in deeply weathered slates of the Sierra de Collcerola (Catalonian Coastal Range, Spain)

Stability conditions in an area located NW of Barcelona (Spain) are discussed. Here, several mass movements were observed, mainly affecting weathered Paleozoic slates. Many of these failures involved slopes cut along recent infrastructures: debris flows, wedge and plane failures, generally surficial, occurred more frequently. After a detailed geological and geomorphologic survey, geomechanic characterization was carried out, according to RMR and SMR classifications. This rating gave a prediction of slope behaviour, in fairly good agreement with the real observed one.

Stability numerical analysis was carried out for the main cut slopes, based upon the Limit Equilibrium Method. First of all, the deterministic factor of safety was computed using the mean values of parameters. After that, a simulation technique based upon the Monte Carlo Method was applied in order to obtain factor of safety distributions. The probability of failure was estimated as P(F<1).

Finally, results from deterministic and probabilistic approaches were compared. The effectiveness of different possible remedial measures was highlighted by means of a sensitivity analysis, which showed that the more important parameters in the study area are the geometrical ones (height, slope and failure plane angles). The final technical solutions adopted are briefly outlined.     

Rheological implications of the strong seismic reflector in the Kakkonda geothermal field, Japan

We re-processed the seismic reflection survey data of the Kakkonda geothermal field. The pre-stack migration delineates a strong and continuous reflector between 1800- and 2800-m depth, below which formations are not reflective. Earthquake data exhibit seismicity in the upper crust. The lower boundary of seismogenic layer is interpreted as the brittle–ductile transition. The thermal structure is thought to be the major factor controlling its depth. We compared the strong reflector with the thermal and rheological structure from drillholes. The depth of the reflector corresponds to the top of the highly–very highly fractured zone observed from formation microscanner imagery (FMI) logging in the Miocene formations. The density of fracture in the Kakkonda granite is very low, suggesting that granite corresponds to the nonreflective zone. The temperature–depth profile of well WD-1a shows that the temperature at the highly–very highly fractured zone is about 350 °C. This corresponds to a hydrothermal convection zone filled with two-phase geothermal fluid. The cut-off depth of seismicity that indicates the brittle–ductile transition lies at the isotherm of 300–350 °C near the reflector. We conclude that the strong seismic reflector is a strong contrast in acoustic impedance at the top of the fractured layer. The fractured layer could be a decoupling plane caused by different tectonic behaviors between the upper brittle and the lower ductile layers or a dehydration front by thermal diffusion. The similarity between the strong reflector and K-horizon, the strong reflector, found in southern Tuscany, Italy suggests that the P-wave reflector at the top of highly fractured zone at the brittle–ductile transition be common in areas with magmatic activity.     

Tectonic implications of damage zone-related fault-fracture networks revealed in drill core through the Nojima fault, Japan

1800 m of drill core through the Nojima fault zone, Japan, reveals subsidiary fault and fracture networks that developed in the fault zone that triggered the 1995 Ms 7.2 Kobe earthquake. The subsidiary fault zones contain a fault gouge of < 1 cm bounded by thin zones of foliated cataclasite or breccia. Fractures are filled with calcite veins, calcite-cemented breccias, clay, and iron-oxide and carbonate alternation of the granitic host rock. These features are typical of extensional fractures that form the conduit network for fluid flux close to a major fault zone. The zone of distributed deformation surrounding the main fault is 50 m in width, and the dip of the Nojima fault at > 1 km depth is 75°. The fault-fracture networks associated with the Nojima fault zone are coseismic and were filled with carbonate and fine-grained material during repeated seismic-related infiltration of the fault zone by carbonate-bearing subsurface water. This study shows that fault-related fracture networks plays an important role as fluid flow conduits within seismically active faults, and can change in character from zones of high permeability to low permeability due to cementation and/or pore collapse.     

Occurrence of lawsonite in pelitic schists from the Sanbagawa metamorphic belt, central Shikoku, Japan

The occurrence of lawsonite is described from pelitic schists of the lower-grade part of the pumpellyite-bearing subzone of the chlorite zone in the Asemi River area of central Shikoku. The lawsonite-bearing parageneses are consistent with the generally accepted view that the Sanbagawa facies series represents higher pressures than the lawsonite-bearing facies series in New Zealand.     

Seismic-reflection profiles of the central part of the Clarendon–Linden fault system of western New York in relation to regional seismicity

Geological and geophysical research in upstate New York, with few exceptions, has not definitively associated seismicity with specific Proterozoic basement or Paleozoic bedrock structures. The central part of the Clarendon–Linden fault system (CLFS) between Batavia and Dale, NY is one of those exceptions where seismicity has been studied and has been spatially associated with structure. The CLFS is either a complex system of long faults with associated shorter branches and parallel segments, or a region of many short faults aligned north–south from the Lake Ontario shore southward to Allegany County, NY. Interpretation of 38 km of Vibroseis and approximately 56 km of conventional seismic-reflection data along 13 lines suggests that the CLFS is a broad zone of small faults with small displacements in the lower Paleozoic bedrock section that is at least 77 km long and 7–17 km wide and spatially coincident with a north-trending geophysical (combined aeromagnetic and gravity) lineament within the basement. The relative offset across the faults of the system is more than 91 m near Attica, NY. The CLFS is the expression of tectonic crustal adjustments within the Paleozoic rock above the boundary of two basement megablocks of differing petrologic provinces and differing earthquake characteristics that forms the eastern side of the Elzevir–Frontenac boundary zone. Deep seismic-reflection profiles display concave-eastward listric faults that probably merge at depth near the mid-crustal boundary layer. An interpretive vertical section provides the setting for refined definitions of the CLFS, its extensions at depth and its relation to seismicity. Most modern seismicity in western New York and the Niagara Peninsula of Ontario occurs in apparent patterns of randomly dispersed activity. The sole exception is a line of seven epicenters of small earthquakes that trend east from Attica, NY into the Rochester basement megablock. Earthquakes may be triggered at the intersections of north- and east-trending brittle faults within the Niagara basement megablock. Current interpretations of the mechanisms for earthquake generation in western New York and the Niagara Peninsula of Ontario require conservative estimates of seismic hazards that assume that an earthquake the size of the 1929 Attica, NY, event (M_b=5.2) or larger could occur anywhere in the Eastern Great Lakes Basin (EGLB). The broad zone of small-displacement faults that marks the CLFS in the lower Paleozoic sedimentary section and the uppermost basement may not provide the structural environment for generation of earthquakes in western New York. If this interpretation is correct, most seismicity is generated within the Niagara basement megablock beneath or west of the CLFS. Consequently, we may have to look to the deeper tectonic regime of basement megablocks to understand the distribution of modern seismicity in the EGLB.     

The higher grade metamorphic zonation of the Sambagawa metamorphic belt in central Shikoku, Japan

The higher grade metamorphic zonation of the Sambagawa (= Sanbagawa) belt is established for the first time for the whole area of central Shikoku. As discontinuous reactions to define the isograd are absent, the metamorphic grade is primarily determined by the Mg-Fe partitioning between garnet and chlorite along representative traverses. However, for regional mapping, mineralogical features of the pelitic schists, such as using mineral assemblages of more than divariant equilibrium, the modal garnet to chlorite ratio, and the optical properties of chlorite, are employed as auxiliary criteria.
The presence of the highest grade mineral zone in the middle of the structural level is confirmed, but its spatial distribution is far more complex than hitherto accepted. Thermal axes are now confirmed at three different structural levels. A model is presented in which the stacking of thrust sheets of different grade took place while metamorphic reactions were in progress. Thermal readjustment brought a continuous metamorphic temperature gradient across and within the thrust sheets. Tectonic blocks of metagabbro and ultramafic rock were emplaced synchronously with thinning and subsequently also re-equilibrated. Local anomalies of metamorphic grade, represented by mixing of schists of different metamorphic grade, exist, but they are due to a later stage event.     

Southward extrusion of eclogite-bearing mafic–ultramafic bodies in the Sanbagawa belt, central Shikoku, Japan

Several mafic rock masses, which have experienced eclogite facies metamorphism, are distributed in flat-lying non-eclogitic schists in an intermediate structural level (thermal core) of the Sanbagawa belt. The largest, Iratsu mass, and an associated peridotite, the Higashi-Akaishi mass, extend E–W for about 8 km, and N–S for about 3 km, and are surrounded by pelitic, basic and quartz schists. The Iratsu mass consists of metabasites of gabbroic and basaltic origin, with intercalations of ultramafic rocks, felsic gneiss, quartz schist and metacarbonate. The Iratsu mass can be divided into two layers along a WNW-trending metacarbonate layer. The Higashi-Akaishi mass consists of peridotite with intercalations of garnet clinopyroxenite. It is situated beneath the western half of the Iratsu mass, and their mutual boundary dips gently or steeply to the N or NE. These masses underwent eclogite, and subsequent epidote-amphibolite facies metamorphism as has been reported elsewhere. The Iratsu–Higashi-Akaishi masses and the surrounding rocks underwent ductile deformation under epidote-amphibolite facies (or lower P–T) metamorphic conditions. Their foliation generally trends WNW and dips moderately to the NNE, and the mineral lineation mostly plunges to the N and NE. In non-eclogitic schists surrounding the Iratsu–Higashi-Akaishi masses, the foliation generally trends WNW and dips gently or steeply to the N or S and the mineral lineation mostly plunges to the NW, N and NE. Kinematic analysis of deformation structures in outcrops and oriented samples has been performed to determine shear senses. Consistent top-to-the-north, normal fault displacements are observed in peridotite layers of the Higashi-Akaishi mass and eclogite-bearing epidote amphibolite layers of the Iratsu mass. Top-to-the-northeast or top-to-the-northwest displacements also occur in non-eclogitic pelitic–quartz schists on the northern side of the Iratsu mass. In the structural bottom of the Iratsu–Higashi-Akaishi masses and to the south, reverse fault (top-to-the-south) movements are recognized in serpentinized peridotite and non-eclogitic schists. These observations provide the following constraints on the kinematics of the rock masses: (1) northward normal displacement of Iratsu relative to Higashi-Akaishi, (2) northward normal displacement of non-eclogitic schists on the north of the Iratsu mass and (3) southward thrusting of the Iratsu–Higashi-Akaishi masses upon non-eclogitic schists in the south. The exhumation process of the Iratsu–Higashi-Akaishi masses can be explained by their southward extrusion.