日本南海海槽发生罕遇地震情况下我国华东沿海的海啸危险性研究

基于日本南海海槽地震活动性和历史海啸事件记载的分析,本文对日本南海海槽发生M_W9.1罕遇地震情况下的海啸进行了数值模拟研究.结果表明:该地震可引发初始波幅约10 m的海啸,6个小时后传至浙江沿海,近岸各处波幅为1—2 m;8个小时后靠近上海海岸线,最大波幅约2 m,受地形影响局地爬高至近3 m;11个小时后抵达苏北黄海沿岸,预计波幅普遍在1 m左右.海啸的上岸高度与海岸附近的海深和海岸线的形态密切相关.我国近岸海域地形变化复杂,海湾众多,对海啸波有放大作用,该模拟结果可能比实际传播到近岸时偏小,因此综合评估日本海啸影响我国华东地区的规模m可达1—2级左右.一旦日本南海发生罕遇地震对我国的影响不容忽视,尤其遇上风暴潮与天文大潮叠加,则可能会造成一定程度的海啸灾害.     

Structure and evolution of the backstop in the eastern Nankai Trough area (Japan): Implications for the soon-to-come Tokai earthquake

Abstract We present data showing that the intra-oceanic shortening now occurring south of the eastern Nankai Trough that has produced the Zenisu Ridge has also been responsible for the formation of a previous ridge now buried below the continental margin. This ridge, that we refer to as Paleo-Zenisu, is presently adjacent to the backstop and its location coincides with the outer limit of the seismogenic decollement. The subduction of the paleo-Zenisu ridge below the wedge has led to its complete reorganization and has given its identity to the Great Tokai earthquake rupture zone. The formation of paleo-Zenisu and its consequent subduction has induced the tilting of the backstop toward the northwest since ca 2 Ma. This model suggests that the backstop and possibly the wedge are dextrally sheared because they are extruded southwestward in relation to the collision of the Izu-Bonin Ridge with Japan. We use the finite motion from Zenisu to paleo-Zenisu to derive both the subduction vectors along the Nankai Trough and the shortening vectors within Zenisu-Izu. The amount of shortening absorbed within Zenisu-Izu increases toward the northeast. The corresponding subduction vectors of the Zenisu platelet below the wedge decrease accordingly to the northeast from 50 to less than 20 mm/year and the Zenisu body rotates clockwise with a pole near 36° North, 139° East. This might explain the apparent longer repetition time of great earthquakes in the Tokai area. On the other hand, the 25-35 mm/year obtained for the rate of shortening along the Zenisu thrust indicates a high seismic potential there.     

Variations in sediment thickness and type along the northern Philippine Sea Plate at the Nankai Trough

Abstract   We documented regional and local variations in basement relief, sediment thickness, and sediment type in the Shikoku Basin, northern Philippine Sea Plate, which is subducting at the Nankai Trough. Seismic reflection data, tied with ocean drilling program drill cores, reveal that variations in the incoming sediment sequences are correlated with basement topography. We mapped the three-dimensional seismic facies distribution and measured representative seismic sequences and units. Trench-parallel seismic profiles show three regional provinces in the Shikoku Basin that are distinguished by the magnitude of basement relief and sediment thickness: Western (<200–400 m basement relief, >600 m sediment thickness), Central (>1500 m relief, ∼2000 m sediments), and Eastern (<600 m relief, ∼1200 m sediments) provinces. The total thickness of sediment in basement lows is as much as six times greater than that over basement highs. Turbidite sedimentation in the Shikoku Basin reflects basement control on deposition, leading to the local presence or absence of turbidite units deposited during the middle Oligocene to the middle Miocene. During the first phase of sedimentation, most basement lows were filled with turbidites, resulting in smooth seafloor morphology that does not reflect basement relief. A second phase of turbidite deposition in the Eastern Province was accompanied by significant amounts of hemipelagic sediments interbedded with turbidite layers compared to the other provinces because of its close proximity to the Izu–Bonin Island Arc. Both regional and local variations in basement topography and sediment thickness/type have caused lateral heterogeneities on the underthrusting plate that will, in turn, influence lateral fluid flow along the Nankai accretionary prism.     

Tsunami deposits from the 1993 Southwest Hokkaido earthquake and the 1640 Hokkaido Komagatake eruption,northern Japan

The southwest Hokkaido tsunami of July 12th, 1993, left continuous onshore sand deposits along the west coast of Oshima Peninsuka, Hokkaido, northern Japan. We investigated spatial distribution and lithofacies of the new tsunami deposits for its identification of ancient tsunami deposits. An eyewitness acount and bent plants helped our interpretation of the onshore tsunami behavior. We regard the following properties as typical of the coastal tsunami sand deposits: (1) The deposits cover the surface almost continuously on gentle topography. (2) Deposit thicknesses and mean grain sizes descrease with distance from the sea. (3) Deposit thicknesses and lithofacies vary greatly across local surface undulation. (4) Graded bedding reflecting tsunami runup and backwash is present in thick deposits. (5) The deposits are widely distributed along the coast and extend inland several tens of meters to 100 m. We examined a candidate for the paleo-tsunami deposits associated with the 1640 Komagatake eruption, and confirmed that the similar patterns are typical of ancient tsunami deposits.     

Millennium‐scale recurrent uplift inferred from beach deposits bordering the eastern Nankai Trough,Omaezaki area,central Japan

A flight of Holocene marine terraces on the southwestern coast of Cape Omaezaki of central Japan provides evidence of recurrent millennium‐scale uplift events. We reconstructed the uplift history of these terraces by using facies analysis of drill core and geoslicer samples, environmental analysis of trace fossils, and ¹⁴C age determinations. Coastal uplift can be identified by the displacement of beach deposits such as foreshore deposits, which represent the intertidal swash zone of a wave‐dominated sandy coast. Three levels of former beach deposits facing the Nankai Trough were identified near the coast in the Omaezaki area. The highest of these, dated at about 3020–2880 BC, records a maximum of 2.2–2.7 m of emergence. The middle beach surface, of minimum age 370–190 BC, shows 1.6–2.8 m of emergence. The lowest beach surface, which is older than 1300–1370 AD, records 0.4–1.6 m of emergence. Our analysis of vertical crustal deformation data during the Holocene in this region suggests that rapid and strong uplift was restricted to the southwestern coast of the Omaezaki area and was probably caused by high‐angle thrusting on subsidiary faults branching from the underlying plate boundary megathrust.     

Tectonics and sedimentation around Kashinosaki Knoll: A subducting basement high in the eastern Nankai Trough

Abstract   When seamounts and other topographic highs on an oceanic plate are subducted, they cause significant deformation of the overriding plate and may act as asperities deeper in the seismogenic zone. Kashinosaki Knoll (KK) is an isolated basement high of volcanic origin on the subducting Philippine Sea Plate that will soon be subducted at the eastern Nankai Trough. Seismic reflection imaging reveals a thick accumulation of sediments (∼1200 m) over and around the knoll. The lower portion of the sedimentary section has a package of high-amplitude, continuous reflections, interpreted as turbidites, that lap onto steep basement slopes but are parallel to the gentler basement slopes. Total sediment thickness on the western and northern slopes is approximately 40–50% more than on the summit and southeastern slopes of KK. These characteristics imply that the basal sedimentary section northwest of KK was deposited by infrequent high-energy turbidity currents, whereas the area southeast of KK was dominated by hemipelagic sedimentation over asymmetric basement relief. From the sediment structure and magnetic anomalies, we estimate that the knoll likely formed near the spreading center of the Shikoku Basin in the early Miocene. Its origin differs from that of nearby Zenisu Ridge, which is a piece of the Shikoku Basin crust uplifted along a thrust fault related to the collision of the Izu–Bonin arc and Honshu. KK has been carried into the margin of the Nankai Trough, and its high topography is deflecting Quaternary trench turbidites to the south. When KK collides with the accretionary prism in about 1 My, the associated variations in sediment type and thickness around the knoll will likely result in complex local variations in prism deformation.     

Hydrothermal sediments associated with a relict back-arc spreading center in the Shikoku Basin, recovered from the Nankai accretionary prism, Japan

The hemipelagic mudrocks of the Nankai accretionary prism, Japan, contain hydrothermal deposits associated with a relict spreading center in the Shikoku Basin. Initial work on core samples from Ocean Drilling Program site 808 found several samples with elevated concentrations of calcium, magnesium, iron and manganese, at depths of between 1060 and 1111 m below sea floor. However, the origin of these sediments was uncertain, due to a lack of data. There was no recorded evidence of whether these elevated concentrations were present throughout this interval of core, or if they were present as discrete layers with the background hemipelagic mudrocks in between. In the present study the core was resampled, and the sediments with anomalous chemical compositions were found to be present in discrete layers. This fact, along with a detailed interpretation of their geochemistry, has allowed them to be identified as hydrothermal sediments, associated with the relict spreading center in the Shikoku Basin. The lower (older) two layers display a chemical composition typical of umbers, while the upper (younger) two layers are metalliferous mudrocks typical of deposits found further from the spreading center.     

Stratigraphic records of tsunamis along the Japan Sea,southwest Hokkaido,northern Japan

The stratigraphy of tsunami deposits along the Japan Sea, southwest Hokkaido, northern Japan, reveals tsunami recurrences in this particular area. Sandy tsunami deposits are preserved in small valley plains, whereas gravelly deposits of possible tsunami origin are identified in surficial soils covering a Holocene marine terrace and a slope talus. At least five horizons of tsunami events can be defined in the Okushiri Island, the youngest of which immediately overlies the Ko‐d tephra layer (1640 AD) and was likely formed by the historical Oshima‐Ohshima tsunami in 1741 AD. The four older tsunami deposits, dated using accelerator mass spectrometry ¹⁴C, were formed at around the 12th century, 1.5–1.6, 2.4–2.6, and 2.8–3.1 ka, respectively. Tsunami sand beds of the 1741 AD and circa 12th century events are recognized in the Hiyama District of Hokkaido Island, but the older tsunami deposits are missing. The deposits of these two tsunamis are found together at the same sites and distributed in regions where wave heights of the 1993 tsunami (Hokkaido Nansei‐oki earthquake, Mw = 7.7) were less than 3 m. Thus, the 12th century tsunami waves were possibly generated near the south of Okushiri Island, whereas the 1993 tsunami was generated towards the north of the island. The estimated recurrence intervals of paleotsunamis, 200–1100 years with an average of 500 years, likely represents the recurrence interval of large earthquakes which would have occurred along several active faults offshore of southwest Hokkaido.     

Structure of the Nankai accretionary prism as revealed from IZANAGI sidescan imagery and multichannel seismic reflection profiling

Abstract The Nankai accretionary prism, off southwest Japan represents one of the best developed clastic prisms in the world. A combination of swath mapping including Sea Beam and 'IZANAGI' sidescan sonar and closely spaced seismic reflection data was used to investigate the relationship between the progressive landward change in surface morphology and the internal structural evolution of the prism. The prism surface is divided into three zones sub-parallel to the trough axis on the basis of the IZANAGI backscattering image. The frontal part of the prism is characterized by several continuous lineaments that are approximately perpendicular to the plate convergence direction. These lineaments correspond to anticlinal ridges caused by active imbricate thrusting. Landward, these anticlinal ridges become progressively masked by fine-grained hemipelagic slope sediments that are constantly supplied to the entire prism slope. However, these overlying sediments show little deformation. This implies a change in deformation style from frontal thrusting with fault-bend folds to internal refolding of thrust sheets. In the middle to upper prism slope, the IZANAGI image shows numerous landslide features and large fault scarps, suggesting that exposed sediments are lithified enough to fail in brittle mode compared with the wet sediment deformation at the prism toe. Prism evolution is strongly affected by the decollement depth which may be indirectly controlled by oceanic basement relief; a topographic embayment coincides with a regional minimum of sediment offscraping where a basement high has been subducted. The small tapered prism observed in the embayment may be due to lateral supply of overpressured pore fluids from the adjacent prism. Strain caused by the differential rate of prism growth across the basement relief forms faults trending at high angles to the trough axis.     

Simple Spring-mass Model Simulation of Earthquake Cycle along the Nankai Trough in Southwest Japan

— Along the Nankai trough in southwest Japan, due to the subduction of the Philippine Sea plate, great earthquakes have occurred repeatedly. The rupture zone is divided into five segments. Historical documents show the characteristic features of the past several earthquake cycles, such as almost simultaneous rupture occurrence in segments in spite of different convergence rates, the recurrence time of 90 to 150 years, the existence of segment pairing in earthquake ruptures and the different coseismic slip behaviors in respective segments. Based on the rate and state friction law, we simulate these features with a simple block-spring model to investigate the physical mechanisms of earthquake cycle. Considering the actual fault parameters related to the geometry and the kinematics of the convergent plate in five segments, we calculate the model parameters for the corresponding blocks in the simulation. The features of the observed earthquake cycle are successfully reproduced by assigning the other following model parameters; (1) the stick-slip periods are the same (more than 150 years) for non-interacting blocks, (2) the different pairs of frictional parameters a-b and D _c are assigned in each segment, (3) the interactions between segments are large, (4) the convergence rate in the eastern Tokai segment is about half of those in the other segments.     

Source Fault Model of the 1771 Yaeyama Tsunami, Southern Ryukyu Islands, Japan, Inferred from Numerical Simulation

The 1771 Yaeyama tsunami is successfully reproduced using a simple faulting model without submarine landslide. The Yaeyama tsunami (M 7.4), which struck the southern Ryukyu Islands of Japan, produced unusually high tsunami amplitudes on the southeastern coast of Ishigaki Island and caused significant damage, including 12,000 casualties. Previous tsunami source models for this event have included both seismological faults and submarine landslides. However, no evidence of landslides in the source has been obtained, despite marine surveying of the area. The seismological fault model proposed in this study, describing a fault to the east of Ishigaki Island, successfully reproduces the distribution of tsunami runup on the southern coast of the Ryukyu Islands. The unusual runup heights are found through the numerical simulation attributable to a concentration of tsunami energy toward the southeastern coast of Ishigaki Island by the effect of the shelf to the east. Thus, the unusual runup heights observed on the southeastern coast of Ishigaki Island can be adequately explained by a seismological fault model with wave-ray bending on the adjacent shelf.     

Oblique and near collision subduction, Sagami and Suruga Troughs —preliminary results of the French-Japanese 1984 Kaiko cruise, Leg 2

Leg 2 of the French-Japanese 1984 Kaiko cruise has surveyed the Suruga and the Sagami Troughs, which lie on both sides of the northwestward moving and colliding Izu-Bonin Ridge, the northernmost part of the Philippine Sea plate. The transition from the Nankai Trough to the Suruga Trough is characterized by northward decrease in width of the accretionary prism, in good agreement with the increasing obliquity between the through axis and the direction of the convergence, as the strike of the convergent boundary changes from ENE-NNE to south-north. South of the area, the southern margin of the Zenisu Ridge shows contractional deformations. This supports the interpretation made by the team of Leg 1 who studied the western extension of the area we studied, that it is an intra-oceanic thrusting of the ridge over the Shikoku Basin. In the Sagami Trough, where the relative motion is highly oblique to the plate boundary, active subduction is mostly confined in the east-west trending portions of the trough located south of the Boso Peninsula and along the lower Boso Canyon, near the TTT triple junction. In between, the present motion is mainly right-lateral along the northwest trending Boso escarpment. However, an inactive but recent (Pliocene to lower Pleistocene) accretionary prism exists south of the Boso escarpment, which suggests that the relative motion was more northerly than at present before about 1 Ma ago.     

Determination and application of path duration of seismic ground motions based on the K-NET data in Sagami Bay,Japan

Duration models are one of the important parameters in ground-motion simulations. This model varies in different study areas, and plays a critical role in nonlinear structural response analysis. Currently, available empirical models are being globally used in ground-motion simulations, with limited research focusing on path duration in specific regions. In this study, we collected 6,486 sets of three-component strong-motion records from 29 K-NET stations in the Sagami Bay, Japan, and its surrounding areas between January 2000 to October 2018. We extracted the effective duration of 386 pieces of ground-motion records by manually picking up the S-wave arrival time and calculating the significant duration. We then obtained the path duration model of the study area based on the empirical equation of dynamic corner frequency and source duration of [7]. Compared with the results of the available empirical models, the Fourier spectrum of the simulated ground motion from our effective duration model showed higher accuracy in the long-term range, with less fitting residuals. This path duration model was then applied to simulate two earthquakes of M_W5.4 and M_W6.2, respectively, in the region using the stochastic finite-fault method with a set of reliable source, path, and site parameters determined for the study area. The simulation results of most stations fit well with observation records in the 0–30 Hz frequency band. For the M_W5.4 earthquake, the simulated ground motions at KNG005/KNG010/SZO008 stations were relatively weak in the mid to high frequency band (1–30 Hz) because the quality factor and geometric diffusion model used in the simulation were the averages of the entire Sagami Bay region, causing a bias in the results of a few stations owing to local crustal velocity anomalies and topographic effects. For the M_W6.2 earthquake, the simulated ground motions were relatively weak at all SZO and TKY stations, mainly because of the close distance from these stations to the epicenter and the complex seismic-wave propagation paths. The analysis suggests that the differences between the simulation results of the two earthquakes were mainly related to complex geological conditions and seismic-wave propagation paths.     

Great earthquake at 7.3 ka inferred from tsunami deposits in the Sukumo Bay area,Southwestern Japan

Tsunami deposits in Kyushu Island, Southwestern Japan, have been attributed to the 7.3 ka Kikai caldera eruption, but their origin has not been confirmed. We analyzed an 83-cm-thick Holocene event deposit in the SKM core, obtained from incised valley fill in the coastal lowlands near Sukumo Bay, Southwestern Shikoku Island. We confirmed that the event deposit contains K-Ah volcanic ash from the 7.3 ka eruption. The base of the event deposit erodes the underlying inner-bay mud, and the deposit contains material from outside the local terrestrial and marine environment, including angular quartz porphyry from a small inland exposure, oyster shell debris, and a coral fragment. Benthic foraminifers and ostracods in the deposit indicate various habitats, some of which are outside Sukumo Bay. The sand matrix contains low-silica volcanic glass from the late stage of the Kikai caldera eruption. We also documented the same glass in an event deposit in the MIK1 core, from the incised Oyodo River valley in the Miyazaki Plain on Southeastern Kyushu. These two 7.3 ka tsunami deposits join other documented examples that are widely distributed in Southwestern Japan including the Bungo Channel and Beppu Bay in Eastern Kyushu, Tachibana Bay in Western Kyushu, and Zasa Pond on the Kii Peninsula as well as around the caldera itself. The tsunami deposits near the caldera have been divided into older and younger 7.3 ka tsunami deposits, the younger ones matching the set of widespread deposits. We attribute the younger 7.3 ka tsunami deposits to a large tsunami generated by a great interplate earthquake in the Northern part of the Ryukyu Trench and (or) the Western Nankai Trough just after the late stage of the Kikai caldera eruption and the older 7.3 ka tsunami deposits to a small tsunami generated by an interplate earthquake or Kikai caldera eruption.     

Dewatering and extensional deformation of the Shikoku Basin hemipelagic sediments in the Nankai Trough

The Shikoku Basin hemipelagic sequence, which underlies the Nankai Trough wedge, S.W. Japan, thins by 50% between the outer edge of the trench wedge and DSDP Site 582, 14 km arcward. A sedimentation model, which incorporates changes in sedimentation rates with time and with distance from the trench wedge toe, indicates that 57% of the total thinning occurs as a result of temporally varying sedimentation rates and a time transgressive facies boundary between the trench wedge turbidites and the underlying hemipelagites. Burial-induced consolidation beneath the wedgeshaped turbiditic overburden, accounts for the remaining 43% of arcward thinning within the hemipelagic unit. Rapid dewatering, modeled as one-dimensional consolidation suggests that the excess pore water pressures are quite low during this progressive dewatering. Thus, high pore water pressures should not be assumed to occur universally in convergent margin settings. Normal faults and vein structures in the hemipelagites suggest near-horizontal extension in addition to vertical consolidation. The estimates of excess pore water pressures together with fault geometries and horizontal extensional strains, determined from the geometry of the subducting oceanic plate, can be used to constrain the stress conditions at failure. The expulsion of hot water from the rapidly dewatering sediments in the Nankai Trough may help to explain anomalously high heat flow in the central part of the trough.     

The origin of strata within the inner accretionary prism of Nankai Trough: Evidence from clay mineral assemblages along the NanTroSEIZE transect

One of the more prominent architectural elements of the Nankai subduction margin, offshore southwest Japan, is an out‐of‐sequence thrust fault (megasplay) that separates the inner accretionary prism from the outer prism. The inner prism (hanging wall of the megasplay) is dominated by mudstone, which is enigmatic when the sedimentary facies is compared to coeval deposits in the Shikoku Basin (i.e. inputs from the subducting Philippine Sea plate) and to coarser‐grained turbidite sequences from the Quaternary trench wedge. Clay mineral assemblages amplify the mismatches of sedimentary facies. Mudstones from the inner prism are uniformly depleted in smectite, with average bulk values of 23–24 wt%, whereas the Shikoku Basin deposits show progressive decreases in proportions of smectite over time, from averages of 46–48 wt% at 10 Ma to 17–21 wt% at 1 Ma. Plate‐boundary reconstructions for the Philippine Sea region provide one solution to the conundrum. Between 15 Ma and 10 Ma, the Pacific plate subducted near the NanTroSEIZE transect, and a trench‐trench‐trench triple junction migrated to the northeast. Accretion during that period involved sediments that had been deposited on the Pacific plate. Motion of the Philippine Sea plate changed from 10 Ma to 6 Ma, resulting in sinistral slip along the proto‐Nankai Trough. Sediments accreted during that period probably had been deposited near the triple junction, with a hybrid detrital provenance. Renewed subduction of the Philippine Sea plate at 6 Ma led to reorganization of watersheds near the Izu–Honshu collision zone and gradual incision of large submarine canyons on both sides of the colliding Izu arc. Accreted Pliocene mudstones share more of an affinity to the triple junction paleoenvironment than they do to Shikoku Basin. These differences between subducting Shikoku Basin strata and accreted Pacific plate sediments have important implications for interpretations of frictional properties, structural architecture, and diagenetic fluid production.     

An integrated geodynamic model of the Nankai subduction zone and neighboring regions from geophysical inversion and modeling

I investigate large-scale deep crustal structures of the Nankai subduction zone and neighboring region using regional magnetic and gravity anomalies, heat flow measurements, and earthquake hypocenters. It is found that ages, dip angles, and geothermal states of the subducting slab have direct influences on mantle wedge serpentinization. The weakest serpentinization observed in the Nankai forearc region is associated with the youngest downgoing plate of the Shikoku Basin. Conspicuous gravity anomalies identified in the forearc region are coincidental spatially with magnetic anomalies after the reduction to the pole, a mathematical procedure that helps relocate magnetic sources and boundaries, and allows us to more easily interpret magnetic data. It is argued that these patches of magnetic and gravity anomalies are caused by the same sources of anomalous density and magnetization, and are linked directly to preexisting structures such as magnetic anomalies and their boundaries in the subducting oceanic crust. Since the gravity and magnetic anomaly patches are found to be closely related to interplate seismogenic behaviors in the Nankai subduction zone, I suggest that major magnetic boundaries in the Shikoku Basin are likely weak places for slab tears that trigger seismic segmentations along the subduction zone.     

Deep-sea submersible survey in the Suruga, Sagami and Japan Trenches: preliminary results of the 1985 Kaiko cruise, Leg 2

Nine submersible dives were made in three trenches off central Japan, between 2990 and 5900 m of water depth. Our observations confirm the interpretation that Daiichi-Kashima Seamount is a Cretaceous guyot formed on the Pacific plate that has traveled into the Japan Trench. We also confirmed the previous interpretation of a large normal fault that splits the seamount in two halves, the lower one being now subducting beneath the Japan margin. Compressional deformation was identified within the lower part of the inner slope in front of the seamount. The pattern of deformation that affects Quaternary sediments is in agreement with the present kinematics of the convergence between the Pacific plate and Japan. Deep-water (5700 m) clam colonies are associated with advection of fluids, driven by the subduction-related overpressures. In the northern slope of the Boso Canyon, along the Sagami Trough system (Philippine Sea plate-Japan boundary), the deformation affecting a thick upper Miocene to lower Pliocene sequence indicates two directions of shortening: a N175°E direction which is consistent with the present relative motion along the Sagami Trough (N285–N300°E) and a N30°E direction which could be related to a more northerly direction of convergence that occured during the early Quaternary and earlier.     

Rejuvenated extension of the Philippine Sea plate and its effect on subduction dynamics in the Nankai Trough

The Nankai Trough, Japan, is a subduction zone characterized by the recurrence of disastrous earthquakes and tsunamis. Slow earthquakes and associated tremor also occur intermittently and locally in the Nankai Trough and the causal relationship between slow earthquakes and large earthquakes is important to understanding subduction zone dynamics. The Nankai Trough off Muroto, Shikoku Island, near the southeast margin of the rupture segment of the 1946 Nankai earthquake, is one of three regions where slow earthquakes and tremor cluster in the Nankai Trough. On the Philippine Sea plate, the rifting of the central domain of the Shikoku Basin was aborted at ~15 Ma and underthrust the Nankai forearc off Muroto. Here, the Tosa-Bae seamount and other high-relief features, which are northern extension of the Kinan Seamount chain, have collided with and indented the forearc wedge. In this study, we analyzed seismic reflection profiles around the deformation front of accretionary wedge and stratigraphically correlated them to drilling sites off Muroto. Our results show that the previously aborted horst-and-graben structures, which were formed around the spreading center of the Shikoku Basin at ~15 Ma, were rejuvenated locally at ~6 Ma and more regionally at ~3.3 Ma and have remained active since. The reactivated normal faulting has enhanced seafloor roughness and appears to affect the locations of slow earthquakes and tremors. Rejuvenated normal faulting is not limited to areas near the Nankai Trough, and extends more than 200 km into the Shikoku Basin to the south. This extension might be due to extensional forces applied to the Philippine Sea plate, which appear to be driven by slab-pull in the Ryukyu and Philippine trenches along the western margin of the Philippine Sea plate.     

The Tsunami of 26 December, 2004: Numerical Modeling and Energy Considerations

A numerical model for the global tsunamis computation constructed by Kowalik et al. (2005), is applied to the tsunami of 26 December, 2004 in the World Ocean from 80°S to 69°N with spatial resolution of one minute. Because the computational domain includes close to 200 million grid points, a parallel version of the code was developed and run on a Cray X1 supercomputer. An energy flux function is used to investigate energy transfer from the tsunami source to the Atlantic and Pacific Oceans. Although the first energy input into the Pacific Ocean was the primary (direct) wave, reflections from the Sri Lankan and eastern shores of Maldives were a larger source. The tsunami traveled from Indonesia, around New Zealand, and into the Pacific Ocean by various routes. The direct path through the deep ocean to North America carried miniscule energy, while the stronger signal traveled a considerably longer distance via South Pacific ridges as these bathymetric features amplified the energy flux vectors. Travel times for these amplified energy fluxes are much longer than the arrival of the first wave. These large fluxes are organized in the wave-like form when propagating between Australia and Antarctica. The sources for the larger fluxes are multiple reflections from the Seychelles, Maldives and a slower direct signal from the Bay of Bengal. The energy flux into the Atlantic Ocean shows a different pattern since the energy is pumped into this domain through the directional properties of the source function. The energy flow into the Pacific Ocean is approximately 75% of the total flow to the Atlantic Ocean. In many locations along the Pacific and Atlantic coasts, the first arriving signal, or forerunner, has lower amplitude than the main signal which often is much delayed. Understanding this temporal distribution is important for an application to tsunami warning and prediction.