Study of the Nankai seismogenic fault using dynamic wave propagation modelling of digital rock from the Nobeoka Fault
Chandoeun Eng 1 Tatsunori Ikeda 2 Takeshi Tsuji 1 2 3
+ Author Affiliations
- Author Affiliations
1 Department of Earth Resources Engineering, Kyushu University, 744 Motooka Nishi-ku, Fukuoka 819-0395, Japan.
2 International Institute for Carbon-Neutral Energy Research (I2CNER), Kyushu University, 744 Motooka, Nishi-ku, Fukuoka 819-0395, Japan.
3 Corresponding author. Email: tsuji@mine.kyushu-u.ac.jp
Exploration Geophysics 49(1) 11-20 https://doi.org/10.1071/EG17129
Submitted: 27 September 2017 Accepted: 2 October 2017 Published: 25 October 2017
Abstract
To understand the characteristics of the Nankai seismogenic fault in the plate convergent margin, we calculated the P- and S-wave velocities (VP and VS) of digital rock models constructed from core samples of an ancient plate boundary fault at Nobeoka, Kyushu Island, Japan. We first constructed 3D digital rock models from microcomputed tomography images and identified their heterogeneous textures such as cracks and veins. We replaced the cracks and veins with air, water, quartz, calcite and other materials with different bulk and shear moduli. Using the Rotated Staggered Grid Finite-Difference Method, we performed dynamic wave propagation simulations and quantified the effective VP, VS and the ratio of VP to VS (VP/VS) of the 3D digital rock models with different crack-filling minerals. Our results demonstrate that the water-saturated cracks considerably decreased the seismic velocity and increased VP/VS. The VP/VS of the quartz-filled rock model was lower than that in the water-saturated case and in the calcite-filled rock model. By comparing the elastic properties derived from the digital rock models with the seismic velocities (e.g. VP and VP/VS) around the seismogenic fault estimated from field seismic data, we characterised the evolution process of the deep seismogenic fault. The high VP/VS and low VP observed at the transition from aseismic to coseismic regimes in the Nankai Trough can be explained by open cracks (or fractures), while the low VP/VS and high VP observed at the deeper coseismic fault zone suggests quartz-filled cracks. The quartz-rich fault zone characterised as low VP/VS and high VP in this study could partially relate to the coseismic behaviour as suggested by previous studies, because quartz exhibits slip-weakening behaviour (i.e. unstable coseismic slip).
Key words: earthquake, fractures, modelling, rock physics, wave propagation.
References
Ando, M., 1975, Source mechanisms and tectonic significance of historical earthquake derived from tsunami data: Tectonophysics, 27, 119–140| Source mechanisms and tectonic significance of historical earthquake derived from tsunami data:Crossref | GoogleScholarGoogle Scholar |
Andra, H., Combaret, N., Dvorkin, J., Glatt, E., Han, J., Kabel, M., Keehm, Y., Krzikalla, F., Lee, M., Madonna, C., and Marsh, M., 2013, Digital rock physics benchmarks part II: computing effective properties: Computers & Geosciences, 50, 33–43
| Digital rock physics benchmarks part II: computing effective properties:Crossref | GoogleScholarGoogle Scholar |
Audet, P., and Burgmann, R., 2014, Possible control of subduction zone slow-earthquake periodicity by silica enrichment: Nature, 510, 389–392
| Possible control of subduction zone slow-earthquake periodicity by silica enrichment:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BC2cXhtVSgsL%2FJ&md5=d86005108127bbbc9d367140f2b4ed60CAS |
Cerjan, C., Kosloff, D., Kosloff, R., and Reshef, M., 1985, A nonreflecting boundary for discrete acoustic and elastic wave equation: Geophysics, 50, 705–708
| A nonreflecting boundary for discrete acoustic and elastic wave equation:Crossref | GoogleScholarGoogle Scholar |
Di Toro, G., Goldsby, D. L., and Tullis, T. E., 2004, Friction falls to towards zero in quartz rock as slip velocity approaches: Nature, 427, 436–439
| 1:CAS:528:DC%2BD2cXntFWhsQ%3D%3D&md5=634710659c0f9dd8323fa8fb60cef4bdCAS |
DiTullio, L., Laughland, M., Byrne, T., and Underwood, M., 1990, Deformation paths in the shallow levels of an accretionary prism: the Eocene Shimanto belt of southwest Japan: Geological Society of America Bulletin, 102, 1420–1438
| Deformation paths in the shallow levels of an accretionary prism: the Eocene Shimanto belt of southwest Japan:Crossref | GoogleScholarGoogle Scholar |
Fisher, D., and Byrne, T., 1992, Strain variations in an ancient accretionary complex: implications for forearc evolution: Tectonics, 11, 330–347
| Strain variations in an ancient accretionary complex: implications for forearc evolution:Crossref | GoogleScholarGoogle Scholar |
Fukuchi, R., Fujimoto, K., Kameda, J., Hamahashi, M., Yamaguchi, A., Kimura, G., Hamada, Y., Hashimoto, Y., Kitamura, Y., and Saito, S., 2014, Changes in illite crystallinity within an ancient tectonic boundary thrust caused by thermal, mechanical, and hydrothermal effects: an example from the Nobeoka Thrust, southwest Japan: Earth, Planets and Space, 66, 116
| Changes in illite crystallinity within an ancient tectonic boundary thrust caused by thermal, mechanical, and hydrothermal effects: an example from the Nobeoka Thrust, southwest Japan:Crossref | GoogleScholarGoogle Scholar |
Gardner, G. H. F., Gardner, L. W., and Gregory, A. R., 1974, Formation velocity and density - the diagnostic basics for stratigraphic traps: Geophysics, 39, 770–780
| Formation velocity and density - the diagnostic basics for stratigraphic traps:Crossref | GoogleScholarGoogle Scholar |
Goldsby, D. L., and Tullis, T. E., 2002, Low frictional strength of quartz rocks at subseismic slip rates: Geophysical Research Letters, 29, 25-1–25-4
| Low frictional strength of quartz rocks at subseismic slip rates:Crossref | GoogleScholarGoogle Scholar |
Hamahashi, M., Saito, S., Kimura, G., Yamaguchi, A., Fukuchi, R., Kameda, J., Hamada, Y., Kitamura, Y., Fujimoto, K., Hashimoto, Y., Hina, S., and Eida, M., 2013, Contrasts in physical properties between the hanging wall and footwall of an exhumed seismogenic megasplay fault in a subduction zone-an example from the Nobeoka thrust drilling project: Geochemistry, Geophysics, Geosystems, 14, 5354–5370
| Contrasts in physical properties between the hanging wall and footwall of an exhumed seismogenic megasplay fault in a subduction zone-an example from the Nobeoka thrust drilling project:Crossref | GoogleScholarGoogle Scholar |
Hamahashi, M., Hamada, Y., Yamaguchi, A., Kimura, G., Fukuchi, R., Saito, S., Kameda, J., Kitamura, Y., Fujimoto, K., and Hashimoto, Y., 2015, Multiple damage zone structure of an exhumed seismogenic megasplay fault in a subduction zone – a study from the Nobeoka Thrust Drilling Project: Earth, Planets and Space, 67, 30
| Multiple damage zone structure of an exhumed seismogenic megasplay fault in a subduction zone – a study from the Nobeoka Thrust Drilling Project:Crossref | GoogleScholarGoogle Scholar |
Hyndman, R. D., Wang, K., and Yamano, M., 1995, Thermal constraints on the seismogenic portion of the southwestern Japan subduction thrust: Journal of Geophysical Research, 100, 15373–15392
| Thermal constraints on the seismogenic portion of the southwestern Japan subduction thrust:Crossref | GoogleScholarGoogle Scholar |
Jiang, F., and Tsuji, T., 2015, Impact of interfacial tension on residual CO2 clusters in porous sandstone: Water Resources Research, 51, 1710–1722
| Impact of interfacial tension on residual CO2 clusters in porous sandstone:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BC2MXms1OjsbY%3D&md5=8d6ca4297ef77956a0a43ed7678885e5CAS |
Jiang, F., and Tsuji, T., 2016, Estimation of three-phase relative permeability by simulating fluid dynamics directly on rock-microstructure images: Water Resources Research, 53, 11–32
| Estimation of three-phase relative permeability by simulating fluid dynamics directly on rock-microstructure images:Crossref | GoogleScholarGoogle Scholar |
Kamei, R., Pratt, R. G., and Tsuji, T., 2012, Waveform tomography imaging of a megasplay fault system in the seismogenic Nankai subduction zone: Earth and Planetary Science Letters, 317–318, 343–353
| Waveform tomography imaging of a megasplay fault system in the seismogenic Nankai subduction zone:Crossref | GoogleScholarGoogle Scholar |
Kamei, R., Pratt, R. G., and Tsuji, T., 2013, On acoustic waveform tomography of wide-angle OBS data – strategies for pre-conditioning and inversion: Geophysical Journal International, 194, 1250–1280
| On acoustic waveform tomography of wide-angle OBS data – strategies for pre-conditioning and inversion:Crossref | GoogleScholarGoogle Scholar |
Kamei, R., Pratt, R. G., and Tsuji, T., 2014, Misfit functionals in Laplace-Fourier domain waveform inversion, with application to wide-angle ocean bottom seismograph data: Geophysical Prospecting, 62, 1054–1074
| Misfit functionals in Laplace-Fourier domain waveform inversion, with application to wide-angle ocean bottom seismograph data:Crossref | GoogleScholarGoogle Scholar |
Kimura, G., Hina, S., Hamada, Y., Kameda, J., Tsuji, T., Kinoshita, M., and Yamaguchi, A., 2012, Runaway slip to the trench due to rupture of highly pressurized megathrust beneath the middle trench slope: the tsunamigenesis of the 2011 Tohoku earthquake off the east coast of northern Japan: Earth and Planetary Science Letters, 339–340, 32–45
Kodaira, S., Lidaka, T., Kato, A., Park, J. O., Iwasaki, T., and Kaneda, Y., 2004, High pore fluid pressure may cause silent slip in the Nankai Trough: Science, 304, 1295–1298
| High pore fluid pressure may cause silent slip in the Nankai Trough:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BD2cXkt1Ggurc%3D&md5=ed819b04955f1290935d01e6308d61c6CAS |
Kondo, H., Kimura, G., Masago, H., Ohmori-Ikehara, K., Kitamura, Y., Ikesawa, E., Sakaguchi, A., Yamaguchi, A., and Okamoto, O., 2005, Deformation and fluid flow of a major out-of-sequence thrust located at seismogenic depth in an accretionary complex: Nobeoka Thrust in the Shimanto Belt, Kyushu, Japan: Tectonics, 24, TC6008
| Deformation and fluid flow of a major out-of-sequence thrust located at seismogenic depth in an accretionary complex: Nobeoka Thrust in the Shimanto Belt, Kyushu, Japan:Crossref | GoogleScholarGoogle Scholar |
Kumagai, H., 1996, Time sequence and the recurrence models for large earthquakes along the Nankai trough revisited: Geophysical Research Letters, 23, 1139–1142
| Time sequence and the recurrence models for large earthquakes along the Nankai trough revisited:Crossref | GoogleScholarGoogle Scholar |
Madonna, C., Almqvist, B. S. G., and Saenger, E. H., 2012, Digital rock physics: numerical prediction of pressure-dependent ultrasonic velocities using micro-CT imaging: Geophysical Journal International, 189, 1475–1482
| Digital rock physics: numerical prediction of pressure-dependent ultrasonic velocities using micro-CT imaging:Crossref | GoogleScholarGoogle Scholar |
Martin, K. M., Gulick, S. P. S., Bangs, N. L. B., Moore, G. F., Ashi, J., Park, J. O., Kuramoto, S., and Taira, A., 2010, Possible strain partitioning structure between the Kumano fore-arc basin and the slope of the Nankai Trough accretionary prism: Geochemistry, Geophysics, Geosystems, 11, Q0AD02
| Possible strain partitioning structure between the Kumano fore-arc basin and the slope of the Nankai Trough accretionary prism:Crossref | GoogleScholarGoogle Scholar |
Mavko, G., Mukerji, T., and Dvorkin, J., 1998, The rock physics handbook: tools for seismic analysis in porous media: Cambridge University Press.
Miyazaki, S. I., and Heki, K., 2001, Crustal velocity field of southwest Japan: subduction and arc-arc collision: Journal Geophysics Research, 106, 4305–4326
| Crustal velocity field of southwest Japan: subduction and arc-arc collision:Crossref | GoogleScholarGoogle Scholar |
Moore, J. C., and Allwardt, A., 1980, Progressive deformation of a Tertiary trench slope, Kodiak Islands, Alaska: Journal of Geophysical Research, 85, 4741–4756
| Progressive deformation of a Tertiary trench slope, Kodiak Islands, Alaska:Crossref | GoogleScholarGoogle Scholar |
Moore, J. C., and Saffer, D., 2001, Updip limit of the seismogenic zone beneath the accretionary prism of southwest Japan: an effect of diagenetic to low-grade metamorphic processes and increasing effective stress: Geology, 29, 183–186
| Updip limit of the seismogenic zone beneath the accretionary prism of southwest Japan: an effect of diagenetic to low-grade metamorphic processes and increasing effective stress:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BD3MXhtlynurY%3D&md5=1f338aff5c62186a4428fcf0bd73b0acCAS |
Moore, J. C., and Vrolijk, P., 1992, Fluids in accretionary prisms: Reviews of Geophysics, 30, 113–135
| Fluids in accretionary prisms:Crossref | GoogleScholarGoogle Scholar |
Obana, K., and Kodaira, S., 2009, Low-frequency tremors associated with reverse faults in a shallow accretionary prism: Earth and Planetary Science Letters, 287, 168–174
| Low-frequency tremors associated with reverse faults in a shallow accretionary prism:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BD1MXht1Wqt7vF&md5=40ed0929349c9eb52fea92df9837df85CAS |
Obara, K., and Ito, Y., 2005, Very low frequency earthquakes excited by the 2004 off the Kii peninsula earthquakes: a dynamic deformation process in the large accretionary prism: Earth, Planets and Space, 57, 321–326
| Very low frequency earthquakes excited by the 2004 off the Kii peninsula earthquakes: a dynamic deformation process in the large accretionary prism:Crossref | GoogleScholarGoogle Scholar |
Park, J. O., Tsuru, T., Kodaira, S., Cummins, P. R., and Kaneda, Y., 2002, Splay fault branching along the Nankai subduction zone: Science, 297, 1157–1160
| Splay fault branching along the Nankai subduction zone:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BD38XmsVOhsrk%3D&md5=6971785fe74bcbab20b6dc5839dbe2d1CAS |
Park, J. O., Naruse, H., and Bangs, N. L., 2014, Along-strike variations in the Nankai shallow décollement properties and their implications for tsunami earthquake generation: Geophysical Research Letters, 41, 7057–7064
| Along-strike variations in the Nankai shallow décollement properties and their implications for tsunami earthquake generation:Crossref | GoogleScholarGoogle Scholar |
Research Group for Active Faults of Japan, 1991, Active faults in Japan: sheet maps and inventories: University of Tokyo Press.
Saenger, E. H., 2008, Numerical methods to determine effective elastic properties: International Journal of Engineering Science, 46, 598–605
| Numerical methods to determine effective elastic properties:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BD1cXkvV2ltrk%3D&md5=6220e4f03d0df322dbdc8682bee8a418CAS |
Saenger, E. H., Gold, N., and Shapiro, S. A., 2000, Modeling the propagation of elastic waves using a modified finite-difference grid: Wave Motion, 31, 77–92
| Modeling the propagation of elastic waves using a modified finite-difference grid:Crossref | GoogleScholarGoogle Scholar |
Saenger, E. H., Kruger, O., and Shapiro, S., 2004, Effective elastic properties of randomly fractured soils: 3D numerical experiments: Geophysical Prospecting, 52, 183–195
| Effective elastic properties of randomly fractured soils: 3D numerical experiments:Crossref | GoogleScholarGoogle Scholar |
Saenger, E. H., Enzmann, F., Keehm, Y., and Steeb, H., 2011, Digital rock physics: effect of fluid viscosity on effective elastic properties: Journal of Applied Geophysics, 74, 236–241
| Digital rock physics: effect of fluid viscosity on effective elastic properties:Crossref | GoogleScholarGoogle Scholar |
Saffer, D. M., and Tobin, H. J., 2011, Hydrogeology and mechanics of subduction zone forearcs: fluid flow and pore pressure: Annual Review of Earth and Planetary Sciences, 39, 157–186
| Hydrogeology and mechanics of subduction zone forearcs: fluid flow and pore pressure:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BC3MXnsFOnurw%3D&md5=eec289921d6b60707f32e558307d85b4CAS |
Scholz, C. H., 1998, Earthquakes and friction lows: Nature, 391, 37–42
| Earthquakes and friction lows:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DyaK1cXjvFWnug%3D%3D&md5=87a3c6c7415f426f348e0e4583dfbb33CAS |
Screaton, E., Saffer, D., Henry, P., Hunze, S., and Leg 190 Shipboard Scientific Party 2002, Porosity loss within the underthrust sediments of the Nankai accretionary complex: implications for overpressures: Geology, 30, 19–22
| Porosity loss within the underthrust sediments of the Nankai accretionary complex: implications for overpressures:Crossref | GoogleScholarGoogle Scholar |
Seno, T., Stein, S., and Gripp, A. E., 1993, A model for the motion of the Philippine Sea plate consistent with Nuvel-1 and geological data: Journal of Geophysical Research, 98, 17941–17948
| A model for the motion of the Philippine Sea plate consistent with Nuvel-1 and geological data:Crossref | GoogleScholarGoogle Scholar |
Tobin, H. J., and Saffer, D. M., 2009, Elevated fluid pressure and extreme mechanical weakness of a plate boundary thrust, Nankai Trough subduction zone: Geology, 37, 679–682
| Elevated fluid pressure and extreme mechanical weakness of a plate boundary thrust, Nankai Trough subduction zone:Crossref | GoogleScholarGoogle Scholar |
Tsuji, T., Kimura, G., Okamoto, S., Kono, F., Mochinaga, H., Saeki, T., and Tokuyama, H., 2006, Modern and ancient seismogenic out-of-sequence thrusts in the Nankai accretionary prism: comparison of laboratory-derived physical properties and seismic reflection data: Geophysical Research Letters, 33, L18309
| Modern and ancient seismogenic out-of-sequence thrusts in the Nankai accretionary prism: comparison of laboratory-derived physical properties and seismic reflection data:Crossref | GoogleScholarGoogle Scholar |
Tsuji, T., Tokuyama, H., Costa Pisani, P., and Moore, G., 2008, Effective stress and pore pressure in the Nankai accretionary prism off the Muroto Peninsula, southwestern Japan: Journal of Geophysical Research, 113, B11401
| Effective stress and pore pressure in the Nankai accretionary prism off the Muroto Peninsula, southwestern Japan:Crossref | GoogleScholarGoogle Scholar |
Tsuji, T., Dvorkin, J., Mavko, G., Nakata, N., Matsuoka, T., Nakanishi, A., Kodaira, S., and Nishizawa, O., 2011, V p/V s ratio and shear-wave splitting in the Nankai Trough seismogenic zone: insights into effective stress, pore pressure and sediment consolidation: Geophysics, 76, WA71–WA82
| V p/V s ratio and shear-wave splitting in the Nankai Trough seismogenic zone: insights into effective stress, pore pressure and sediment consolidation:Crossref | GoogleScholarGoogle Scholar |
Tsuji, T., Kamei, R., and Pratt, R. G., 2014a, Pore pressure distribution of a mega-splay fault system in the Nankai Trough subduction zone: insight into up-dip extent of the seismogenic zone: Earth and Planetary Science Letters, 396, 165–178
| Pore pressure distribution of a mega-splay fault system in the Nankai Trough subduction zone: insight into up-dip extent of the seismogenic zone:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BC2cXotFyhtL4%3D&md5=bf1deeeb25f7f6464b067b8baa55d3abCAS |
Tsuji, T., Ashi, J., and Ikeda, Y., 2014b, Strike-slip motion of a mega-splay fault system in the Nankai oblique subduction zone: Earth, Planets and Space, 66, 120
| Strike-slip motion of a mega-splay fault system in the Nankai oblique subduction zone:Crossref | GoogleScholarGoogle Scholar |
Tsuji, T., Jiang, F., and Christensen, K., 2016, Characterization of immiscible fluid displacement processes with various capillary numbers and viscosity ratios in 3D natural sandstone: Advances in Water Resources, 95, 3–15
| Characterization of immiscible fluid displacement processes with various capillary numbers and viscosity ratios in 3D natural sandstone:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BC28XkvVeltrc%3D&md5=26ff21cd0c8090d543f2c988b1a8ed5cCAS |
Wang, K., and Hu, Y., 2006, Accretionary prisms in subduction earthquake cycles: the theory of dynamic Coulomb wedge: Journal of Geophysical Research, 111, B06410
| Accretionary prisms in subduction earthquake cycles: the theory of dynamic Coulomb wedge:Crossref | GoogleScholarGoogle Scholar |
Yamabe, H., Tsuji, T., Liang, Y., and Matsuoka, T., 2016, Influence of fluid displacement patterns on seismic velocity during supercritical CO2 injection: Simulation study for evaluation of the relationship between seismic velocity and CO2 saturation: International Journal of Greenhouse Gas Control, 46, 197–204
| Influence of fluid displacement patterns on seismic velocity during supercritical CO2 injection: Simulation study for evaluation of the relationship between seismic velocity and CO2 saturation:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BC28XitFChtrs%3D&md5=5afd71284f89ea0e78de829c378fbbb1CAS |
Yamaguchi, A., Cox, S. F., Kimura, G., and Okamoto, S., 2011, Dynamic changes in fluid redox state associated with episodic fault rupture along a megasplay fault in a subduction zone: Earth and Planetary Science Letters, 302, 369–377
| Dynamic changes in fluid redox state associated with episodic fault rupture along a megasplay fault in a subduction zone:Crossref | GoogleScholarGoogle Scholar | 1:CAS:528:DC%2BC3MXhtlKms7s%3D&md5=fd830093a5e59ca0b8e03e21458014d2CAS |