Researcher Database

Researcher Profile and Settings

Master

Affiliation (Master)

  • Faculty of Science Institute of Seismology and Volcanology

Affiliation (Master)

  • Faculty of Science Institute of Seismology and Volcanology

researchmap

Profile and Settings

Degree

  • Doctor of Science(The University of Tokyo)
  • Master of Science(Kyoto University)

Profile and Settings

  • Name (Japanese)

    Murai
  • Name (Kana)

    Yoshio
  • Name

    200901038400320476

Alternate Names

Achievement

Research Interests

  • 地球惑星内部物理学   テクトニクス   地震学   Physics of the Earth and Planetary Interiors   Tectonics   Seismology   

Research Areas

  • Natural sciences / Solid earth science

Research Experience

  • 2007 - 2009 Hokkaido University Institute of Seismology and Volcanology Faculty of Science
  • 2009 - 北海道大学大学院理学研究院附属地震火山研究観測センター 准教授
  • 2006 - 2007 Hokkaido University Institute of Seismology and Volcanology Faculty of Science
  • 2006 - 2007 Research Associate,Institute of Seismology and Volcanology, Faculty of Science, Hokkaido University
  • 1998 - 2006 北海道大学大学院理学研究科附属地震火山研究観測センター 助手
  • 1998 - 2006 Research Associate,Institute of Seismology and Volcanology, Graduate School of Science, Hokkaido University
  • 1997 - 1998 Hokkaido University School of Science
  • 1997 - 1998 Research Associate,Laboratory for Ocean Bottom Seismology, Faculty of Science, Hokkaido University
  • 1995 - 1997 Japan Meteorological Agency
  • 1995 - 1997 Technical officer,Japan Meteorological Agency

Committee Memberships

  • 2002 -2008   The Seismological Society of Japan   delegate   The Seismological Society of Japan

Published Papers

  • Evgeny A. Podolskiy, Yoshio Murai, Naoya Kanna, Shin Sugiyama
    The Journal of the Acoustical Society of America 151 (1) 6 - 16 0001-4966 2022/01/03 [Refereed]
  • Evgeny A. Podolskiy, Yoshio Murai, Naoya Kanna, Shin Sugiyama
    SEISMOLOGICAL RESEARCH LETTERS 92 (5) 2850 - 2865 0895-0695 2021/09 [Refereed]
     
    About 70% of Earth's surface is covered by ocean, for which seismic observations are challenging. Seafloor seismology overcame this fundamental difficulty and radically transformed the earth sciences, as it expanded the coverage of seismic networks and revealed otherwise inaccessible features. At the same time, there has been a recent increase in the number of studies on cryoseismology. These have yielded multiple discoveries but are limited primarily to land and ice-surface receivers. Near ice calving fronts, such surface stations are noisy, primarily due to crevassing and wind, are hazardous to maintain, and can be lost due to iceberg calving. To circumvent these issues, we have applied ocean-bottom seismology to the calving front of a tidewater glacier in northwest Greenland. We present details of this experiment, and describe the technical challenges, noise analysis, and examples of recorded data. This includes tide-modulated seismicity with thousands of icequakes per day and the first near-source (similar to 200-640 m) underwater record of a major kilometer-scale calving event in Greenland, which generated a glacial earthquake that was detectable similar to 420 km away. We also identified a decrease in bottom-water temperature, presumably due to modified water stratification driven by extreme Greenland glacial melting, at the end of July 2019. Importantly, we identify glacial sediments as the key reason for the anomalously long (similar to 9.7 hr) delay in the sensor release from the fjord seafloor. Our study demonstrates a methodology to undertake innovative, interdisciplinary, near-source studies on glacier basal sliding, calving, and marine-mammal vocalizations.
  • Evgeny A Podolskiy, Yoshio Murai, Naoya Kanna, Shin Sugiyama
    Nature communications 12 (1) 3929 - 3929 2021/06/24 [Refereed]
     
    Shearing along subduction zones, laboratory experiments on analogue faults, and sliding along glacier beds are all associated with aseismic and co-seismic slip. In this study, an ocean-bottom seismometer is deployed near the terminus of a Greenlandic tidewater glacier, effectively insulating the signal from the extremely noisy surface seismic wavefield. Continuous, tide-modulated tremor related to ice speed is recorded at the bed of the glacier. When noise interference (for example, due to strong winds) is low, the tremor is also confirmed via analysis of seismic waveforms from surface stations. The signal resembles the tectonic tremor commonly observed during slow-earthquake events in subduction zones. We propose that the glacier sliding velocity can be retrieved from the observed seismic noise. Our approach may open new opportunities for monitoring calving-front processes in one of the most difficult-to-access cryospheric environments.
  • Asbjørn Johan Breivik, Jan Inge Faleide, Rolf Mjelde, Ernst R. Flueh, Yoshio Murai
    Tectonophysics 776 0040-1951 2020/02/05 [Refereed]
     
    The Norwegian continental shelf has been through several rift phases since the Caledonian orogeny. Early Cretaceous rifting created the largest sedimentary basins, and Early Cenozoic continental breakup between East Greenland and Europe affected the continental shelf to various degrees. The Lofoten/Vesterålen shelf is located off Northern Norway, bordering the epicontinental Barents Sea to the northeast, and the deep-water Lofoten Basin to the west. An ocean bottom seismometer/hydrophone (OBS) survey was conducted over the shelf and margin areas in 2003 to constrain crustal structure and margin development. This study presents Profile 8-03, located between the islands of Lofoten/Vesterålen and the shelf edge. The wide-angle seismic data were modeled using forward/inverse raytracing to build a crustal velocity-depth transect. Gravity modeling was used to resolve an ambiguity in seismic Moho identification in the southwestern part. Results show a crustal thickness of ~31 km, significantly thicker than what a vintage land station based study suggested. Profile 8-03 and other OBS profiles to the southwest show high sedimentary velocities at or near the seafloor, increasing rapidly with depth. Sedimentary velocities were compared to the velocity-depth function derived from an OBS profile at the Barents Sea margin, tied to a coincident well log, where there is little erosion. Results from this profile and the crossing Profile 6-03 (Breivik et al. 2017) indicate three major erosion episodes; Late Triassic-Early Jurassic, tentatively mid-Cretaceous, Late Cretaceous–early Cenozoic, and a minor late glacial erosion episode off Vesterålen.
  • Miyamachi H., Uchida K., Miyamachi R., Isoda K., Triahadini A, Teguri Y., Kamiya Y., Shimizu H., Katao H., Shibutani T, Miura T., Kobayashi R., Nakagawa J., Yoneda I., Kato S., Tameguri T., Takishita K., Nakai K., Yamashita Y., Maeda Y., Watanabe T., Horikawa S., Yakiwara H., Matsuhiro K., Okuda T., Tsuji S., Sogawa N., Hasegawa D., Nakahigashi K., Kurashimo E., Iwasaki T., Shinohara M., Yamada T., Hirano S., Abe H., Ando M., Tanaka S., Ikezawa S., Sato T., Yamamoto M., Hirahara S., Nakayama T., Azuma R., Hino R., Kubo T., Ohtomo S., Suzuki S., Tsutsui T., Inoue Y., Takei R., Tada Y., Takahashi H., Aoyama H., Ohzono M., Shiina T., Soda M., Takada M., Ichiyanagi M., Yamaguchi T., Ono N., Saito K., Ito C., Susukida Y., Murai Y., Nakagaki T., Tanaka Y., Sakao K., Unno N., Matsushima T.
    PROGRAMME AND ABSTRACTS THE VOLCANOLOGICAL SOCIETY OF JAPAN 特定非営利活動法人 日本火山学会 2019 (0) 44 - 44 2433-5320 2019 [Refereed][Not invited]
  • R. Azuma, R. Hino, Y. Ohta, Y. Ito, K. Mochizuki, K. Uehira, Y. Murai, T. Sato, T. Takanami, M. Shinohara, T. Kanazawa
    Journal of Geophysical Research: Solid Earth 123 (6) 5249 - 5264 2169-9313 2018/06/01 [Refereed][Not invited]
     
    ©2018. The Authors. The 2011 Tohoku-oki earthquake occurred with a rupture length of 500 km along the Japan Trench, causing a large slip (>30 m) at the shallowest portion of the plate boundary fault south of 39°N off Miyagi and a smaller slip (~10 m) north of 39°N off Sanriku, the northern part of the source area. We estimated the P wave velocity (Vp) structure around the shallowest portion of the plate boundary along the trench to investigate the spatial correlation between the structural variation and coseismic shallow slip distribution of the 2011 earthquake. For this purpose, we analyzed data from an air gun-ocean bottom seismometer survey on an along-arc line on the lower part of the landward slope of the Japan Trench. We detected a high Vp zone in the hanging wall of the plate interface off Miyagi, which corresponds to the Cretaceous backstop. The Vp model also showed the absence of the high Vp zone off Sanriku. This suggests that the low-Vp deformed prism is wider off Sanriku than off Miyagi. Therefore, the backstop is close to the trench particularly off Miyagi. We suggest a correlation between this along-arc variation of the hanging wall side structure and the extent of the coseismic shallow slip in the trench axial region; the large slip reached the trench off Miyagi but not off Sanriku.
  • Asbjorn Johan Breivik, Jan Inge Faleide, Rolf Mjelde, Ernst R. Flueh, Yoshio Murai
    TECTONOPHYSICS 718 25 - 44 0040-1951 2017/10 [Refereed][Not invited]
     
    The Early Eocene continental breakup was magma-rich and formed part of the North Atlantic Igneous Province. Extrusive and intrusive magmatism was abundant on the continental side, and a thick oceanic crust was produced up to a few m.y. after breakup. However, the extensive magmatism at the Voring Plateau off mid-Norway died down rapidly northeastwards towards the Lofoten/Vesteralen Margin. In 2003 an Ocean Bottom Seismometer profile was collected from mainland Norway, across Lofoten, and into the deep ocean. Forward/inverse velocity modeling by raytracing reveals a continental margin transitional between magma-rich and magma-poor rifting. For the first time a distinct lower-crustal body typical for volcanic margins has been identified at this outer margin segment, up to 3.5 km thick and similar to 50 km wide. On the other hand, expected extrusive magmatism could not be clearly identified here. Strong reflections earlier interpreted as the top of extensive lavas may at least partly represent high-velocity sediments derived from the shelf, and/or fault surfaces. Early post-breakup oceanic crust is moderately thickened (similar to 8 km), but is reduced to 6 km after 1 m.y. The adjacent continental crystalline crust is extended down to a minimum of 4.5 km thickness. Early plate spreading rates derived from the Norway Basin and the northern Voring Plateau were used to calculate synthetic magnetic seafloor anomalies, and compared to our ship magnetic profile. It appears that continental breakup took place at similar to 53.1 Ma, similar to 1 m.y. later than on the Voring Plateau, consistent with late strong crustal extension. The low interaction between extension and magmatism indicates that mantle plume material was not present at the Lofoten Margin during initial rifting, and that the observed excess magmatism was created by late lateral transport from a nearby pool of plume material into the lithospheric rift zone at breakup time. (C) 2017 Elsevier B.V. All rights reserved.
  • Trond Kvarven, Rolf Mjelde, Berit Oline Hjelstuen, Jan Inge Faleide, Hans Thybo, Ernst R. Flueh, Yoshio Murai
    TECTONOPHYSICS 666 144 - 157 0040-1951 2016/01 [Refereed][Not invited]
     
    The inner part of the volcanic, passive More Margin, mid-Norway, expresses an unusual abrupt thinning from high onshore topography with a thick crust to an offshore basin with thin crystalline crust. Previous P-wave modeling of wide-angle seismic data revealed the presence of a high-velocity (7.7-8.0 km/s) body in the lower crust in this transitional region. These velocities are too high to be readily interpreted as Early Cenozoic intrusions, a model often invoked to explain lower crustal high-velocity bodies in the region. We present a Vp/Vs model, derived from the modeling of wide-angle seismic data, acquired by use of Ocean Bottom Seismograph horizontal components. The modeling suggests dominantly felsic composition of the crust. An average Vp/Vs value for the lower crustal body is modeled at 1.77, which is compatible with a mixture of continental blocks and Caledonian eclogites. The results are compiled with earlier results into a transect extending from onshore Norway to onshore Greenland. Back-stripping of the transect to Early Cenozoic indicates asymmetric conjugate magmatism related to the continental break-up. Further back-stripping to the time when most of the Caledonian mountain range had collapsed indicates that the thinning during the first phase of extension was about 25% higher for proto Norway than proto Greenland. (C) 2015 Elsevier B.V. All rights reserved.
  • Bryndis Brandsdottir, Emilie E. E. Hooft, Rolf Mjelde, Yoshio Murai
    GEOCHEMISTRY GEOPHYSICS GEOSYSTEMS 16 (3) 612 - 634 1525-2027 2015/03 [Refereed][Not invited]
     
    Variations in crustal structure along the 700 km long KRISE7 refraction/reflection and gravity profile, straddling 66.5 degrees N across the Iceland Shelf, Iceland Plateau and western Norway Basin confirm that extinct spreading centers coexisted with the now extinct aegir Ridge prior to the initiation of the Kolbeinsey Ridge at 26 Ma. The western 300 km of the profile, across the Iceland shelf, formed by rifting at the Kolbeinsey Ridge, whereas the eastern 400 km, across the Iceland Plateau and the western Norway Basin, formed by earlier rifting, possibly containing slivers of older oceanic or continental crust rifted off the central E-Greenland margin along with the Jan Mayen Ridge. Crustal thickness increases gradually across the Iceland shelf, from 12 to 13 km near the Kolbeinsey Ridge to 24-28 km near the eastern shelf edge, decreasing abruptly across the shelf edge, to 12-13 km. The Iceland Plateau has crustal thickness ranging from 12 to 15 km decreasing to 5-8 km across the western Norway Basin and 4-5 km at the aegir Ridge. We suggest that high-velocity lower crustal domes and corresponding gravity highs across the Iceland plateau mark the location of extinct rift axes that coexisted with the aegir Ridge. Similar lower crustal domes are associated with the currently active rift segments within Iceland and the Kolbeinsey Ridge.
  • Yukihiro Nakatani, Kimihiro Mochizuki, Masanao Shinohara, Tomoaki Yamada, Ryota Hino, Yoshihiro Ito, Yoshio Murai, Toshinori Sato
    GEOPHYSICAL RESEARCH LETTERS 42 (5) 1384 - 1389 0094-8276 2015/03 [Refereed][Not invited]
     
    The southern limit of the 2011 Tohoku earthquake is considered to be located around off Ibaraki. However, it is not well constrained how far south the large slip extended. To give better constraints, we investigated seismicity including small earthquakes before and after the Tohoku earthquake around off Ibaraki using dense ocean bottom seismic array data. We automatically identified epicenters by backprojecting semblance values resulting in a considerable increase in the number of detected events compared with those listed in the catalog based on onshore observation. The results revealed a couple of seismicity-activated region. The largest aftershock also occurred similar to 30min after the main shock. Our detailed results suggest that this highly activated seismicity was initiated by the largest aftershock instead of the main shock. It, then, suggests that the large coseismic slip zone of the Tohoku earthquake may not have extended off Ibaraki.
  • Asbjorn Breivik, Jan Inge Faleide, Rolf Mjelde, Ernst Flueh, Yoshio Murai
    JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH 119 (9) 6733 - 6755 2169-9313 2014/09 [Refereed][Not invited]
     
    The Voring Plateau off mid-Norway is a volcanic passive margin, located north of the East Jan Mayen Fracture Zone (EJMFZ). Large volumes of magmatic rocks were emplaced during Early Eocene margin formation. In 2003, an ocean bottom seismometer survey was acquired over the margin. One profile crosses from the VOring Plateau to the VOring Spur, a bathymetric high north of the EJMFZ. The P wave data were ray traced into a 2-D crustal velocity model. The velocity structure of the VOring Spur indicates up to 15 km igneous crustal thickness. Magmatic processes can be estimated by comparing seismic velocity (V-P) with igneous thickness (H). This and two other profiles show a positive H-V-P correlation at the VOring Plateau, consistent with elevated mantle temperature at breakup. However, during the first 2 Ma magma production was augmented by a secondary process, possibly small-scale convection. From approximate to 51.5 Ma excess melting may be caused by elevated mantle temperature alone. Seismic stratigraphy around the VOring Spur shows that it was created by at least two uplift events, with the main episode close to the Miocene/Pliocene boundary. Low H-V-P correlation of the spur is consistent with renewed igneous growth by constant, moderate-degree mantle melting, not related to the breakup magmatism. The admittance function between bathymetry and free-air gravity shows that the high is near local isostatic equilibrium, precluding that compressional flexure at the EJMFZ uplifted the high. We find a proposed Eocene triple junction model for the margin to be inconsistent with observations.
  • Trond Kvarven, Jorg Ebbing, Rolf Mjelde, Jan Inge Faleide, Audun Libak, Hans Thybo, Ernst R. Flueh, Yoshio Murai
    TECTONOPHYSICS 626 21 - 40 0040-1951 2014/06 [Refereed][Not invited]
     
    The More Margin in the NE Atlantic represents a dominantly passive margin with an unusual abrupt transition from alpine morphology onshore to a deep sedimentary basin offshore. In order to study this transition in detail, three ocean bottom seismometer profiles with deep seismic reflection and refraction data were acquired in 2009; two dip-profiles which were extended by land stations, and one tie-profile parallel to the strike of the More-Trondelag Fault Complex. The modeling of the wide-angle seismic data was performed with a combined inversion and forward modeling approach and validated with a 3D-density model. Modeling of the geophysical data indicates the presence of a 12-15 km thick accumulation of sedimentary rocks in the More Basin. The modeling of the strike profile located closer to land shows a decrease in crustal velocity from north to south. Near the coast we observe an intra-crustal reflector under the Trondelag Platform, but not under the Slorebotn Sub-basin. Furthermore, two lower crustal high-velocity bodies are modeled, one located near the More Marginal High and one beneath the Slorebotn Sub-basin. While the outer lower crustal body is modeled with a density allowing an interpretation as magmatic underplating, the inner body has a density close to mantle density which might suggest an origin as an eclogized body, formed by metamorphosis of lower crustal gabbro during the Caledonian orogeny. The difference in velocity and extent of the lower crustal bodies seems to be controlled by the Jan Mayen Lineament, suggesting that the lineament represents a pre-Caledonian structural feature in the basement. (C) 2014 Elsevier B.V. All rights reserved.
  • 村井芳夫, 東龍介, 篠原雅尚, 町田祐弥, 山田知朗, 中東和夫, 真保敬, 望月公廣, 日野亮太, 伊藤喜宏, 佐通利典, 塩原肇, 植平賢司, 八木原寛, 尾鼻浩一郎, 高橋成実, 小平秀一, 平田賢治, 対馬弘晃, 岩崎貴哉
    北海道大学地球物理学研究報告 北海道大学大学院理学研究院 (76) 147 - 158 0439-3503 2013/03/19 [Not refereed][Not invited]
     
    2011年3月11日に,太平洋プレートと日本列島を乗せた陸側のプレートとの境界で2011年東北地方太平洋沖地震が発生した.この地震は,日本周辺では観測史上最大のマグニチュード9という巨大地震だった.本震発生後には多数の余震が発生するが,大地震発生のメカニズムを解明するためには,正確な余震分布を調べることが重要である.全国の6つの大学と海洋研究開発機構,気象庁気象研究所は,本震発生直後から共同で100台以上の海底地震計を用いて余震観測を行った.2011年6月中旬までのデータから,震源域全体で約3か月間の精度の良い震源分布が得られた.余震の震源の深さは,全体的に陸に近づくにつれて深くなっていた.震源分布からは,本震時に大きくすべったプレート境界では余震活動が低いことがわかった.上盤の陸側プレート内では余震活動が活発で,正断層型と横ずれ型が卓越していた.太平洋プレート内の余震も多くが正断層型か横ずれ型だった.このことから,日本海溝付近の太平洋プレート内の深部と上盤の陸側プレート内では,本震の発生によって応力場が圧縮場から伸張場に変化したことが示唆される.
  • Gou Fujie, Seiichi Miura, Shuichi Kodaira, Yoshiyuki Kaneda, Masanao Shinohara, Kimihiro Mochizuki, Toshihiko Kanazawa, Yoshio Murai, Ryota Hino, Toshinori Sato, Kenji Uehira
    Earth, Planets and Space 65 (2) 75 - 83 1880-5981 2013 [Refereed][Not invited]
     
    Large destructive interplate earthquakes, such as the 2011 Mw 9.0 Tohoku-oki earthquake, have occurred repeatedly in the northern Japan subduction zone. The spatial distribution of large interplate earthquakes shows distinct along-trench variations, implying regional variations in interplate coupling. We conducted an extensive wide-angle seismic survey to elucidate the along-trench variation in the seismic structure of the forearc and to examine structural factors affecting the interplate coupling beneath the forearc mantle wedge. Seismic structure models derived from wide-angle traveltimes showed significant along-trench variation within the overlying plate. In a weakly coupled segment, (i) the sediment layer was thick and flat, (ii) the forearc upper crust was extremely thin, (iii) the forearc Moho was remarkably shallow (about 5 km), and (iv) the P-wave velocity within the forearc mantle wedge was low, whereas in the strongly coupled segments, opposite conditions were found. The good correlation between the seismic structure and the segmentation of the interplate coupling implies that variations in the forearc structure are closely related to those in the interplate coupling. © The Society of Geomagnetism and Earth, Planetary and Space Sciences (SGEPSS) The Seismological Society of Japan The Volcanological Society of Japan The Geodetic Society of Japan The Japanese Society for Planetary Sciences TERRAPUB.
  • Azuma, Ryosuke, Murai, Yoshio, Katsumata, Kei, Nishimura, Yuichi, Yamada, Takuji, Mochizuki, Kimihiro, Shinohara, Masanao
    Geochemistry Geophysics Geosystems 13 (8) Q08015 - 10.1029/2012GC004135 1525-2027 2012/08 [Refereed][Not invited]
     
    The Kuril Trench subduction zone is one of the most seismogenic regions, where underthrust earthquakes with M > 8 recur along the trench. The seismic gap between the source areas of the 1973 Nemuro-oki and 2003 Tokachi-oki earthquakes, which are typical underthrust earthquakes faulting with rupture velocities of similar to 3 km/s, has been ruptured by the 1952 Tokachi-oki earthquake. The seismic gap has also slipped incidental to neighboring asperities. The difference in slip pattern on the plate interface generally appears as a spatial difference in seismic structure on the plate interface, such as a reflectivity of the plate interface. We estimated the crustal velocity structure and analyzed the reflectivity of the plate interface to investigate the physical properties of the plate interface by performing an air gun-ocean bottom seismometer experiment on the along-trench profile across the seismic gap. Strong reflections from the plate interface were observed in the 1952 Tokachi-oki source area including the seismic gap, rather than in the 1973 Nemuro-oki source area. The strong reflectivity of the plate interface in such the seismic gap with an incidental slip suggests that a slip pattern in the corresponding seismic gap would be conditionally stable. The coupling condition in the source areas of the eastern part of the source area of the 1952 earthquake is different from that in source areas of typical underthrust earthquakes, such as the 2003 Tokachi-oki and 1973 Nemuro-oki earthquakes. Our results suggest that the 1952 Tokachi-oki earthquake was a complex earthquake with the characteristic of a tsunami earthquake.
  • Libak, Audun, Mjelde, Rolf, Keers, Henk, Faleide, Jan Inge, Murai, Yoshio
    Marine Geophysical Research 33 (2) 185 - 207 0025-3235 2012/06 [Refereed][Not invited]
     
    This paper describes results from a geophysical study in the Vestbakken Volcanic Province, located on the central parts of the western Barents Sea continental margin, and adjacent oceanic crust in the Norwegian-Greenland Sea. The results are derived mainly from interpretation and modeling of multichannel seismic, ocean bottom seismometer and land station data along a regional seismic profile. The resulting model shows oceanic crust in the western parts of the profile. This crust is buried by a thick Cenozoic sedimentary package. Low velocities in the bottom of this package indicate overpressure. The igneous oceanic crust shows an average thickness of 7.2 km with the thinnest crust (5-6 km) in the southwest and the thickest crust (8-9 km) close to the continent-ocean boundary (COB). The thick oceanic crust is probably related to high mantle temperatures formed by brittle weakening and shear heating along a shear system prior to continental breakup. The COB is interpreted in the central parts of the profile where the velocity structure and Bouguer anomalies change significantly. East of the COB Moho depths increase while the vertical velocity gradient decreases. Below the assumed center for Early Eocene volcanic activity the model shows increased velocities in the crust. These increased crustal velocities are interpreted to represent Early Eocene mafic feeder dykes. East of the zone of volcanoes velocities in the crust decrease and sedimentary velocities are observed at depths of more than 10 km. The amount of crustal intrusions is much lower in this area than farther west. East of the Knolegga Fault crystalline basement velocities are brought close to the seabed. This fault marks the eastern limit of thick Cenozoic and Mesozoic packages on central parts of the western Barents Sea continental margin.
  • Breivik, Asbjorn Johan, Mjelde, Rolf, Faleide, Jan Inge, Murai, Yoshio
    Geophysical Journal International 188 (3) 798 - 818 0956-540X 2012/03 [Refereed][Not invited]
     
    The Jan Mayen microcontinent (JMMC) in the NE Atlantic was created through two Cenozoic rift episodes. Originally part of East Greenland, the JMMC rifted from NW Europe during the Early Eocene under extensive magmatism. The eastern margin is conjugate to the MoreFaeroes volcanic margin. The western JMMC margin underwent prolonged extension before it finally separated from East Greenland during the Late Oligocene. Here we present the modelling by forward/inverse ray tracing of two wide-angle seismic profiles acquired using Ocean Bottom Seismometers, across the northern and the southern JMMC. Early Eocene breakup magmatism at the eastern JMMC produced an igneous thickness of 79 km in the north, and 1214 km in the south. While the continent is clear in the north, the southern JMMC appears to be affected by later Icelandic magmatism. Reduced seismic velocity and increased crustal thickness are compatible with continental crust adjacent to the volcanic margin in the south, but the continental presence towards the Iceland shelf is less clear. Our magnetic track off the southern JMMC gives seafloor spreading rates comparable to that of the conjugate More Margin. Transition to ultraslow seafloor spreading occurs at similar to 43 Ma, indicating onset of major deformation of the JMMC. Calculating the igneous thicknessmean V P relationship at the eastern volcanic margin gives the typical positive correlation seen elsewhere on the NE Atlantic margins. The results indicate temperature driven breakup magmatism under passive mantle upwelling, with a maximum mantle temperature anomaly of similar to 50 degrees C in the north and 90150 degrees C in the south.
  • Kazuo Nakahigashi, Masanao Shinohara, Kimihiro Mochizuki, Tomoaki Yamada, Ryota Hino, Toshinori Sato, Kenji Uehira, Yoshihiro Ito, Yoshio Murai, Toshihiko Kanazawa
    EARTH PLANETS AND SPACE 64 (12) 1149 - 1156 1343-8832 2012 [Refereed][Not invited]
     
    We present the result of a seismic experiment conducted using ocean bottom seismometers and controlled sources in the region off Ibaraki and the Boso Peninsula. This region is the southern edge of the rupture zone of the 2011 off the Pacific coast of Tohoku Earthquake. We estimated the P-wave seismic velocity structure beneath the profile using a 2-D ray-tracing method. The crustal structure in the southern area is more heterogeneous than that of the northern area. This heterogeneity is thought to be related with subducting the Philippine Sea plate (PHS). The plate boundary between the landward plate and the Pacific plate (PAC) is positioned at depths of 20 km at a distance of 170 km from the southern end of the profile. The subducting PHS is imaged on the southern part of the profile. However, we could not obtain a distinct image of the contact zone of PHS and PAC. The contact zone of PHS and PAC is estimated to have a large heterogeneity resulting from strong deformation due to the collision of the two plates. We infer that the termination of the rupture, and the large afterslip in the collision region, are caused by this strong heterogeneity.
  • Masanao Shinohara, Yuya Machida, Tomoaki Yamada, Kazuo Nakahigashi, Takashi Shinbo, Kimihiro Mochizuki, Yoshio Murai, Ryota Hino, Yoshihiro Ito, Toshinori Sato, Hajime Shiobara, Kenji Uehira, Hiroshi Yakiwara, Koichiro Obana, Narumi Takahashi, Shuichi Kodaira, Kenji Hirata, Hiroaki Tsushima, Takaya Iwasaki
    EARTH PLANETS AND SPACE 64 (12) 1137 - 1148 1343-8832 2012 [Refereed][Not invited]
     
    The 2011 off the Pacific coast of Tohoku Earthquake occurred at the plate boundary between the Pacific plate and the landward plate on March 11, 2011, and had a magnitude of 9. Many aftershocks occurred following the mainshock. Obtaining a precise aftershock distribution is important for understanding the mechanism of earthquake generation. In order to study the aftershock activity of this event, we carried out extensive sea-floor aftershock observations using more than 100 ocean-bottom seismometers just after the mainshock. A precise aftershock distribution for approximately three months over the whole source area was obtained from the observations. The aftershocks form a plane dipping landward over the whole area, nevertheless the epicenter distribution is not uniform. Comparing seismic velocity structures, there is no aftershock along the plate boundary where a large slip during the mainshock is estimated. Activity of aftershocks in the landward plate in the source region was high and normal fault-type, and strike-slip-type, mechanisms are dominant. Within the subducting oceanic plate, most earthquakes have also a normal fault-type, or strike-slip-type, mechanism. The stress fields in and around the source region change as a result of the mainshock.
  • Czuba, Wojciech, Grad, Marek, Mjelde, Rolf, Guterch, Aleksander, Libak, Audun, Krueger, Frank, Murai, Yoshio, Schweitzer, Johannes, IPY Project Grp
    Geophysical Journal International 184 (2) 541 - 554 0956-540X 2011/02 [Refereed][Not invited]
     
    P>A 410 km long Ocean Bottom Seismometer profile spanning from the Bear Island, Barents Sea to oceanic crust formed along the Mohns Ridge has been modelled by use of ray-tracing with regard to observed P-waves. The northeastern part of the model represents typical continental crust, thinned from ca. 30 km thickness beneath the Bear Island to ca. 13 km within the Continent-Ocean-Transition. Between the Hornsund FZ and the Kn circle divide legga Fault, a 3-4 km thick sedimentary basin, dominantly of Permian/Carboniferous age, is modelled beneath the ca. 1.5 km thick layer of volcanics (Vestbakken Volcanic Province). The P-wave velocity in the 3-4 km thick lowermost continental crust is significantly higher than normal (ca. 7.5 km s-1). We interpret this layer as a mixture of mafic intrusions and continental crystalline blocks, dominantly related to the Paleocene-Early Eocene rifting event. The crystalline portion of the crust within the south-western part of the COT consists of a ca. 30 km wide and ca. 6 km thick high-velocity (7.3 km s-1) body. We interpret the body as a ridge of serpentinized peridotites. The magmatic portion of the ocean crust accreted along the Knipovich Ridge from continental break-up at ca. 35 Ma until ca. 20 Ma is 3-5 km thicker than normal. We interpret the increased magmatism as a passive response to the bending of this southernmost part of the Knipovich Ridge. The thickness of the magmatic portion of the crust formed along the Mohns Ridge at ca. 20 Ma decreases to ca. 3 km, which is normal for ultra slow spreading ridges.
  • Yamada, Tomoaki, Nakahigashi, Kazuo, Kuwano, Asako, Mochizuki, Kimihiro, Sakai, Shin'ichi, Shinohara, Masanao, Hino, Ryota, Murai, Yoshio, Takanami, Tetsuo, Kanazawa, Toshihiko
    Earth Planets and Space 63 (7) 841 - 845 1343-8832 2011 [Refereed][Not invited]
     
    The 2011 Tohoku earthquake hit the forearc region of the Japan Trench on March 11, 2011. The rupture zone seemed to reach off the coast of the Kanto region. We had conducted ocean bottom seismographic observations off the coast of the Kanto in 2008 and estimated 851 hypocenter locations around the south part of the 2011 off the Pacific coast of Tohoku Earthquake by using over 50 ocean bottom seismometers (OBSs) and routine data jointly. The hypocenters distributed some clusters, and we found a few seismic gaps at the boundary zones of the clusters. The most remarkable seismic gap was positioned at the edge of intervened Philippine Sea plate (PHS) between the North American plate (NA) and Pacific plate (PAC). We compare the epicenter distributions with the 2011 aftershock distribution determined by routine data. The aftershocks are segmented spatially and there are some seismic gaps among the segments. The remarkable low places are consistent with the boundary zones of each cluster we estimated from the 2008 data. We infer that those regions have strong heterogeneity resulting from strong deformations caused by various subduction processes, such as intervening PHS between NA and PAC, seamount chains and changes in physical properties.
  • Shinohara, Masanao, Yamada, Tomoaki, Nakahigashi, Kazuo, Sakai, Shin'ichi, Mochizuki, Kimihiro, Uehira, Kenji, Ito, Yoshihiro, Azuma, Ryusuke, Kaiho, Yuka, No, Tetsuo, Shiobara, Hajime, Hino, Ryota, Murai, Yoshio, Yakiwara, Hiroshi, Sato, Toshinori, Machida, Yuya, Shinbo, Takashi, Isse, Takehi, Miyamachi, Hiroki, Obana, Koichiro, Takahashi, Narumi, Kodaira, Shuichi, Kaneda, Yoshiyuki, Hirata, Kenji, Yoshikawa, Sumio, Obara, Kazushige, Iwasaki, Takaya, Hirata, Naoshi
    Earth Planets and Space 63 (7) 835 - 840 1343-8832 2011 [Refereed][Not invited]
     
    The 2011 off the Pacific coast of Tohoku Earthquake occurred offshore of northeast Japan region on March 11th, 2011. In order to study the aftershock activity of this event, we started deployment of seventy-two ocean bottom seismometers (OBSs) four days after the mainshock. In the south of the source region, thirty-four long-term OBSs (LT-OBSs) had been deployed before the occurrence of the mainshock, and we recovered three LT-OBSs to clarify the depth distribution of aftershocks. Using the data of OBSs, ninety-nine aftershocks were located. Most of the aftershocks were located in a depth range of 5-30 km and concentrate in the plate boundary region. In addition, aftershocks occurred within the subducting oceanic crust and the 6.2-km/s layer of the landward plate. No aftershocks were found in the mantle of the subducting plate. From the results of a previous seismic survey using OBSs and controlled sources, the subducting Philippine Sea plate is estimated to be in contact with the subducting Pacific plate. The southern end of the seismic activity region of the aftershocks corresponds to the contact region of two subducting plates. We infer that the rupture of the mainshock sequence was terminated at the oceanic plate contact region.
  • Kandilarov, Aleksandre, Landa, Hildegunn, Mjelde, Rolf, Pedersen, Rolf B., Okino, Kyoko, Murai, Yoshio
    Marine Geophysical Research 31 (3) 173 - 195 0025-3235 2010/09 [Refereed][Not invited]
     
    A combined ocean bottom seismometer, multichannel seismic reflection and gravity study has been carried out along the spreading direction of the Knipovich Ridge over a topographic high that defines a segment center. The youngest parts of the crust in the immediate vicinity of the ridge reveal fractured Oceanic Layer 2 and thermally expanded and possibly serpentinized Oceanic Layer 3. The mature part of the crust has normal thickness and seismic velocities with no significant crustal thickness and seismic velocity variations. Mature Oceanic Layer 2 is in addition broken into several rotated fault blocks. Comparison with a profile acquired similar to 40 km north of the segment center reveals significant differences. Along this profile, reported earlier, periods of slower spreading led to generation of thin crust with a high P-wave velocity (Vp), composed of a mixture of gabbro and serpentinized mantle, while periods of faster spreading led to generation of more normal gabbroic crust. For the profile across the segment center no clear relation exists between spreading rate and crustal thickness and seismic velocity. In this study we have found that higher magmatism may lead to generation of oceanic crust with normal thickness even at ultra-slow spreading rates.
  • Breivik, Asbjorn Johan, Mjelde, Rolf, Raum, Thomas, Faleide, Jan Inge, Murai, Yoshio, Flueh, Ernst R.
    Norwegian Journal of Geology 90 (4) 141 - 161 0029-196X 2010 [Refereed][Not invited]
     
    The outer mid-Norwegian margin is characterized by strong breakup magmatism and has been extensively surveyed. The crustal structure of the inner continental shelf, however, is less studied, and its relation to the onshore geology, Caledonian structuring, and breakup magmatism remains unclear. Two Ocean Bottom Seismometer profiles were acquired across the Trondelag Platform in 2003, as part of the Euromargins program. Additional- land stations recorded the marine shots. The P-wave data were modeled by ray-tracing, supported by gravity modeling. Older multi-channel seismic data allowed for interpretation of stratigraphy down to the top of the Triassic. Crystalline basement velocity is similar to 6 km s(-1) onshore. Top basement is difficult to identify offshore, as velocities (5.3-5.7 km s(-1)) intermediate between typical crystalline crust and Mesozoic sedimentary strata appear 50-80 km from the coast. This layer thickens towards the Klakk-Ytreholmen Fault Complex and predates Permian and later structur-ing. The velocities indicate sedimentary rocks, most likely Devonian. Onshore late- to post-Caledonian detachments have been proposed to extend offshore, based on the magnetic anomaly pattern. We do not find the expected correlation between upper basement velocity structure and detachments. However, there is a distinct, dome-shaped lower-crustal body with a velocity of 6.6-7.0 km s(-1). This is thickest under the Froan Basin, and the broad magnetic anomaly used to delineate the detachments correlates with this. The proposed offshore continuation of the detachments thus appears- unreliable. While we find indications of high density and velocity (similar to 7.2 km s(-1)) lower crust under the Ras Basin, similar to the proposed igneous underplating of the outer margin, this is poorly constrained near the end of our profiles. The gravity field indicates that this body may be continuous from the pre-breakup basement structures of the Utgard High to the Froya High, suggesting that it could be an island arc or oceanic terrane- accreted during the Caledonian orogeny. Thus, we find no clear evidence of early Cenozoic igneous underplating of the inner part of the shelf.
  • Mjelde, R., Raum, T., Kandilarov, A., Murai, Y., Takanami, T.
    Tectonophysics 468 (1-4) 224 - 243 0040-1951 2009/04 [Refereed][Not invited]
     
    Five regional Ocean Bottom Seismic profiles from the northwestern part of the More Margin, NE Atlantic, have been modeled with regard to crustal structure, and the interpretations have been supported by the use of multi-channel seismic and potential field data. The chi-squared values (uncertainty) for the studied profiles are generally in the range 1-3, and the uncertainty in the Moho depth is estimated to 1 km. The Continent-Ocean-Transition is found to be relatively sharp (10-40 km wide) expressing a change in physical properties corresponding to felsic continental crust along the landward part of the marginal high (V(p)=6.26.5 km/s), and thick, mafic oceanic crust further northwestwards. The Early Eocene breakup was associated with excessive magmatism, and the models document direct continuity between this magmatism and a lower crustal high-velocity layer in the More Basin (V(p)=7.1-7.3 km/s), suggesting a genetic link. The Jan Mayen Lineament is observed across the entire study area as northeastward termination of the lower crustal body, as well as structuring in the top of the crystalline basement and the Moho. We present a speculative model involving lower crustal bodies of eclogites (V(p)>8 km/s) inferred in the easternmost part of the More Basin linking the Jan Mayen Lineament to the post-Caledonian Nordfjord-Sogn Detachment. (C) 2008 Elsevier B.V. All rights reserved.
  • Becel, Anne, Laigle, Mireille, de Voogd, Beatrice, Hirn, Alfred, Taymaz, Tuncay, Galve, Audrey, Shimamura, Hideki, Murai, Yoshio, Lepine, Jean-Claude, Sapin, Martine, Ozalaybey, Serdar
    Tectonophysics 467 (1-4) 1 - 21 0040-1951 2009/03 [Refereed][Not invited]
     
    Understanding further the nature and evolution at lithospheric scale of the Sea of Marmara on the North Anatolian Fault needs constraints on the deep crustal and Mono spatial variation. This has been probed here with offshore-onshore and OBS, Ocean Bottom Seismometer refraction seismics, in addition to coincident MCS, marine multichannel reflection seismic profiles over the whole North Marmara Trough. The diverse strikes of MCS profiles in a dense grid allow to avoid misinterpretation of late echoes in the deep basin as Moho reflections and attribute them to sidesweeps. Moho is instead positively identified from reversed observations of first-arrival head and refracted waves at the top of the mantle obtained at large offset by land stations. A significant and sharp reduction in its depth, on the order of 5 km occurs beneath both the eastern and western rims of the North Marmara Trough, with a more progressive crustal thinning from the south. The wide-angle reflections on OBS and land stations document in addition to Moho the top of a lower crustal reflective layer, which is also sampled by MCS profiles, and appears to follow Moho topography. The dense grid of MCS profiles along the southwestern margin of the North Marmara Trough reveals a dipping reflector through the upper crust with tilted basement blocks on top. This low-angle fault is suggested as a normal sense detachment extending in depth towards the reflective lower crust. The upwarp of the Moho and lower crustal layer towards the North Marmara Trough suggests that crustal thinning occurs mostly in the upper crustal part, with lateral transport of the material towards WSW in the footwall of the detachment, and possibly other features to the south, in the motion of Anatolia with respect to stable Eurasia oblique to the North Marmara Trough. Thinning can be accommodated in an asymmetric partitioning of the displacement on several branching faults at lithospheric scale. (C) 2008 Elsevier B.V. All rights reserved.
  • Machida, Yuya, Shinohara, Masanao, Takanami, Tetsuo, Murai, Yoshio, Yamada, Tomoaki, Hirata, Naoshi, Suyehiro, Kiyoshi, Kanazawa, Toshihiko, Kaneda, Yoshiyuki, Mikada, Hitoshi, Sakai, Shin'ichi, Watanabe, Tomoki, Uehira, Kenji, Takahashi, Narumi, Nishino, Minoru, Mochizuki, Kimihiro, Sato, Takeshi, Araki, Ei'ichiro, Hino, Ryota, Uhira, Kouichi, Shiobara, Hajime, Shimizu, Hiroshi
    Tectonophysics 465 (1-4) 164 - 176 0040-1951 2009/02 [Refereed][Not invited]
     
    Large earthquakes have repeatedly occurred in the area off southeastern Hokkaido Island, Japan, as the Pacific Plate subducts beneath the island, which is on the North American Plate. The most recent large earthquake in this area, the 2003 Tokachi-oki earthquake (Mw=8.0), occurred on September 26, 2003. In order to investigate aftershock activity in the rupture area, 47 Ocean Bottom Seismometers (OBSs) were quickly deployed after the main shock. In the present study, we simultaneously estimate the hypocenters and 3-D seismic velocity models from the P- and S-wave arrivals of the aftershocks recorded by OBSs. The subducting plate is clearly imaged as a northwest dipping zone in which Vp is greater than 7 km/s, and the relocated hypocenters also show the subducting Pacific Plate. The aftershock distribution reveals that the dip angle of the plate boundary increases abruptly around 90 km from the Kuril Trench. The bending of the subducting plate corresponds to the southeastern edge of the rupture area. The island arc crust on the overriding plate has P-wave velocities of 6-7 km/s and a Vp/Vs of 1.73. A region of Vp/Vs greater than 1.88 was found north of the epicenter of the main shock. The depth of the high Vp/Vs region extends about 10 km upward from the plate interface. The plate boundary just below the high Vp/Vs region has the largest slip at the main rupture. A high Vp anomaly (similar to 7.5 km/s) is found in the island arc crust in northeast part of the study area, which we interpret as a structural boundary related to the arc-arc collisional tectonics of the Hokkaido region, as the rupture of the main shock terminated at this high Vp region. We suggest that the plate interface geometry and the trench-parallel velocity heterogeneity in the landward plate are principal factors in controlling the rupture area of the main shock. (C) 2008 Elsevier B.V. All rights reserved.
  • 北海道大学地球物理学研究報告 (72) 23 - 35 2009 [Not refereed][Not invited]
  • Kandilarov, Aleksandre, Mjelde, Rolf, Okino, Kyoko, Murai, Yoshio
    Marine Geophysical Research 29 (2) 109 - 134 0025-3235 2008/06 [Refereed][Not invited]
     
    The ultra-slow, asymmetrically-spreading Knipovich Ridge is the northernmost part of the Mid Atlantic ridge system. In the autumn of 2002 a combined ocean-bottom seismometer multichannel seismic (OBS/MCS) and gravity survey along the spreading direction of the Knipovich Ridge was carried out. The main objective of the study was to gain an insight into the crustal structure and composition of what is assumed to be an amagmatic segment of oceanic crust. P-wave velocity and Vp/Vs models were built and complemented by a gravity model. The 190 km long transect reveals a much more complex crustal structure than anticipated. The magmatic crust is thinner than the global average of 7.1 +/- 1.0 km. The young fractured portion of Oceanic Layer 2 has low seismic velocities while the older part has normal seismic velocities and is broken into several rotated fault blocks seen as thickness variations of Layer 2. The youngest part of Oceanic Layer 3 is also dominated by low velocities, indicative of fracturing, seawater circulation and thermal expansion. The remaining portion of Layer 3 exhibits inverse variations in thickness and seismic velocity. This is explained by a sequence of periods of faster spreading (estimated to be up to 8 mm/year from interpretation of magnetic anomalies) when more normal gabbroic crust was being generated and periods of slower spreading (5.5 mm/year) when amagmatic stretching and serpentinization of the upper mantle occurred, and crust composed of mixed gabbro and serpentinized mantle was generated. The volumetric changes and upward fluid migration, associated with the process of serpentinization in this part of the crust, caused disruption to the overlying sedimentary layers.
  • Nakamura, Masaki, Yoshida, Yasuhiro, Zhao, Dapeng, Takayama, Hiroyuki, Obana, Koichiro, Katao, Hiroshi, Kasahara, Junzo, Kanazawa, Toshihiko, Kodaira, Shuichi, Sato, Toshinori, Shiobara, Hajime, Shinohara, Masanao, Shimamura, Hideki, Takahashi, Narumi, Nakanishi, Ayako, Hino, Ryota, Murai, Yoshio, Mochizuki, Kimihiro
    Physics of the Earth and Planetary Interiors 168 (1-2) 49 - 70 0031-9201 2008/05 [Refereed][Not invited]
     
    We determined the three-dimensional Vp and Vs structures beneath Japan by applying seismic tomography to a large number of arrival times recorded at temporary stations in the Japan Sea and the Pacific Ocean, as well as those at permanent stations on the Japan Islands. As a result, we obtained more precise seismic images than previous studies. In the crust and the uppermost mantle, southwestern Honshu exhibited weaker heterogeneity than the other areas in Japan, corresponding to the distribution of active volcanoes. Stripe-like heterogeneities exist in the subducting Pacific slab. Relatively low-velocity zones correspond to low-seismicity areas in the Pacific slab, suggesting that the slab is possibly torn or thin around the areas. The fact that nonvolcanic deep tremors associated with the subducting Philippine Sea slab beneath Shikoku, Kii, and Tokai do not occur in zones of high Vp, high Vs, and low Vp/Vs ratio may reflect the existence of fluids generated by the dehydration processes of the slab. Prominent and wide low Vp and Vs zones exist beneath central Honshu at the depth range of 30-60 km, where the volcanic front related to the subducting Pacific plate is located and seismicity around the Philippine Sea plate is very low. This condition may exist because magma genesis processes related to the subducting Pacific plate activate the same processes around the Philippine Sea plate. 2008 Elsevier B.V. All rights reserved.
  • Shinohara, Masanao, Kanazawa, Toshihiko, Yamada, Tomoaki, Nakahigashi, Kazuo, Sakai, Shin'ichi, Hino, Ryota, Murai, Yoshio, Yamazaki, Akira, Obana, Koichiro, Ito, Yoshihiro, Iwakiri, Kazuhiro, Miura, Ryo, Machida, Yuya, Mochizuki, Kimihiro, Uehira, Kenji, Tahara, Michitaka, Kuwano, Asako, Amamiya, Shin'ichiro, Kodaira, Shuichi, Takanami, Tetsuo, Kaneda, Yoshiyuki, Iwasaki, Takaya
    Earth Planets and Space 60 (11) 1121 - 1126 1880-5981 2008 [Refereed][Not invited]
     
    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 degrees. 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.
  • Mjelde, R., Raum, T., Murai, Y., Takanami, T.
    Journal of Geodynamics 43 (3) 374 - 392 0264-3707 2007/04 [Refereed][Not invited]
     
    A tectono-magmatic model across the Voring Plateau, NE Atlantic, has been derived from a densely sampled ocean bottom seismic profile, constrained by multichannel seismic and gravity data. The last phase of rifting, breakup at ca. 55 Ma and first phase of drifting was associated with abundant magmatism, probably related to the arrival of the Icelandic hot-spot. The continent-ocean-transition is modelled over a ca. 25 km wide zone, from an area with clearly defined stretched continental crust, to an area with oceanic crustal velocities and densities. The last phase of rifting is interpreted to be closely related to the development of crustal-scale detachment faults. The initiation of the episode may be referred to as passive rifting, whereas the breakup was actively related to the crustal detachment fault that was closest to the local injection centre con the Voring Plateau. This detachment fault is defined as the continent-ocean-boundary (COB), since it delimits intruded continental crust to the southeast from anomalous oceanic crust to the northwest. The modelled P-wave velocities and crustal thickness indicates active upwelling of mantle and elevated mantle potential temperatures from breakup to ca. 5 Ma spreading. (C) 2006 Elsevier Ltd. All rights reserved.
  • Murai, Yoshio
    Geophysical Journal International 168 (1) 211 - 223 0956-540X 2007/01 [Refereed][Not invited]
     
    We compute the synthetic seismograms of multiply scattered SH waves in 2-D elastic media with densely distributed parallel cracks. We assume randomly distributed cracks in a rectangular-bounded region, which simulate a cracked zone. The crack surfaces are assumed to be stress-free. When the incident wavelength is longer than the crack size, the delay in the arrival of the primary wave is observed at stations beyond the cracked zone and the amplitude of the primary wave is amplified in the cracked zone in the synthetic seismograms. This is because the cracked zone behaves as a low velocity and soft material to the incident long-wavelength wave due to the crack distribution. When the half-wavelength of the incident wave is shorter than the crack length, the scattered waves are clearly observed in the synthetic seismograms and the amplitude of the primary wave is largely attenuated beyond the cracked zone. The calculated attenuation coefficient Q(-1) of the primary wave is directly proportional to the crack density in the range of nu a(2) <= 0.1, where nu and alpha are the number density and half the length of cracks, respectively. This is consistent with that obtained by a stochastic analysis based on Foldy's approximation. A periodic distribution of cracks in a zone is considered as an utterly different model in order to investigate the effect of spatial distributions on the attenuation and dispersion of seismic waves. When cracks are distributed densely, the values of Q(-1) for the periodic crack distribution appear to differ from those for the random distribution of cracks in the low wavenumber range. This suggests that the effect of multiple interactions among densely distributed cracks depends on not only the density but also the spatial distribution of cracks at low wavenumbers. The calculated phase velocity of the primary wave is consistent with that from the stochastic analysis in the range of nu alpha(2) = 0.1 and does not depend on the spatial distribution of cracks and does not depend on the spatial distribution of cracks. This suggests that the multiple crack interactions have a smaller effect to the phase velocity. Therefore, the crack density can be estimated from the values of the phase velocity for the cases of densely distributed cracks even if the effect of the multiple crack interactions is not considered. We can clearly observe the reflected waves in the synthetic seismograms. The elastic constant of a single anisotropic layer equivalent to the cracked zone is derived from the crack density at the long-wavelength limit. The reflection coefficients calculated from the synthetic seismograms are consistent with those of the anisotropic layer calculated from its elastic constants and thickness in low wavenumber range. This means that a fracture zone distributed parallel cracks is considered as an anisotropic layer for long incident wavelengths. Therefore, the elastic constants, crack density and the thickness of the fracture zone can be estimated from the frequency dependence of the reflection coefficients for long incident wavelengths. On the contrary the wavenumber dependence of the reflection coefficients cannot be explained theoretically in the high wavenumber range.
  • Breivik, Asbjorn Johan, Mjelde, Rolf, Faleide, Jan Inge, Murai, Yoshio
    Journal of Geophysical Research-Solid Earth 111 (B7) B07102 - B07102 0148-0227 2006/07 [Refereed][Not invited]
     
    [ 1] In year 2000, an ocean bottom seismometer (OBS) profile was acquired across the More margin to the Aegir Ridge, an extinct seafloor spreading axis. The margin is an early Eocene volcanic passive margin, located between the Faeroe-Iceland Ridge ( FIR) and the East Jan Mayen Fracture Zone (EJMFZ). The P wave data were modeled by ray tracing to give a crustal transect showing a 10 - 11 km thick igneous crust created by breakup magmatism, tapering off to magma-starved seafloor spreading by C23 time (51.4 Ma). The location of the EJMFZ was reinterpreted from a satellite derived gravity map, and spreading direction in the Norway Basin reevaluated. No other fracture zones were confirmed, and both thin oceanic crust ( 4 - 5 km) and lack of fracture zones resemble ultraslow spreading on the Arctic Gakkel Ridge. Magnetic seafloor spreading anomalies were identified from the magnetic track recorded with the OBS profile, and half spreading rates were derived. Early seafloor spreading was slow ( 15 - 32 mm yr(-1)), approaching ultraslow ( 6 - 8 mm yr(-1)) by C20 time (42.7 Ma). AV- shaped pattern seen in the gravity field located only around the northern part of the Aegir Ridge corresponds to increased crustal thickness in the seismic model, recording northeast transport ( 3 - 6 mm yr(-1)) of more melt-fertile asthenosphere zones. The magma-starved character of the Norwegian Basin seen also during slow seafloor spreading may be the result of depletion of the asthenosphere when the Iceland plume constructed the FIR to the south, as the asthenosphere is subsequently transported into the Norway Basin.
  • Hooft, EEE, Brandsdottir, B, Mjelde, R, Shimamura, H, Murai, Y
    Geochemistry Geophysics Geosystems 7 10.1029/2005GC001123  1525-2027 2006/05 [Refereed][Not invited]
     
    [1] We present the results of a seismic refraction experiment that constrains crustal structure and thickness along 225 km of the Kolbeinsey Ridge and Tjornes Fracture Zone and thus quantifies the influence of the Iceland hot spot on melt flux at the spreading center north of Iceland. North of the Iceland shelf, crustal thickness is relatively constant over 75 km, 9.4 +/- 0.2 km. Along the southern portion of the Kolbeinsey Ridge, on the Iceland shelf, crustal thickness increases from 9.5 +/- 0.1 km to 12.1 +/- 0.4 km over 90 km. Gravity inversion indicates a residual crustal gravity anomaly that decreases by about 30 - 40 mGal toward Iceland. We infer that the variations in crustal thickness and gravity are accompanied by mantle temperature changes of 40 degrees to 50 degrees C. At similar distances from the Iceland hot spot, crustal thickness along the Kolbeinsey Ridge is 2 - 2.5 km less than at the Reykjanes Ridge, consistent with the asymmetry in plume-ridge interaction that has been inferred from the axial depth and geochemistry of these ridges. Average lower crustal velocities are also higher along the Kolbeinsey Ridge, consistent with a lower degree of active upwelling than along the Reykjanes Ridge. Topography and crustal thickness patterns at the spreading centers around Iceland are consistent with isostatic support for normal crustal and mantle densities. However, we infer that the lower crust beneath central Iceland is considerably denser than that beneath the adjacent ridges. Crustal thickness and geochemical patterns suggest that deep melting is spatially limited and asymmetric about Iceland while shallow melting is enhanced over a broad region. This asymmetry may be due to a mantle plume that is tilted from south to north in the upper mantle and preferentially melts deeper enriched material beneath the Reykjanes Ridge.
  • Masanao Shinohara, Ei'ichiro Araki, Toshihiko Kanazawa, Kiyoshi Suyehiro, Masashi Mochizuki, Tomoaki Yamada, Kazuo Nakahigashi, Yuka Kaiho, Yoshio Fukao
    ANNALS OF GEOPHYSICS 49 (2-3) 625 - 641 1593-5213 2006/04 [Refereed][Not invited]
     
    Seismological networks provide critical data for better Understanding the dynamics of the Earth; however, a great limitation on existing networks is the uneven distribution of stations. In order to achieve a more uniform distribution of seismic stations, observatories must be Constructed in marine areas. The best configuration for oceanic seismic observatories is thought to be placement of seismometers in deep boreholes. Two deep-sea borehole seismological observatories (WP-1 and WP-2) were Constructed in the Western Pacific and form the initial installations of a 1000 kin span network. At present, seismic records of more than 400 total days were retrieved from both the WP-1 and WP-2. Long-term variations in broadband seismic noise spectra (3 mHz-10 Hz) in the Western Pacific were revealed from these records, and the data showed that ambient seismic noise levels in borehole observatories are comparable to those of the quietest land seismic stations. In addition, there is little temporal variation of noise levels in periods greater than 10 s. Due to this low seismic noise environment, many teleseismic events with magnitudes greater than 5 were recorded. It is confirmed that seismic observation in deepsea borehole gives the best environment for earthquake observation in marine areas.
  • Raum, T, Mjelde, R, Shimamura, H, Mural, Y, Brastein, E, Karpuz, RM, Kravlk, K, Kolsto, HJ
    Tectonophysics 415 (1-4) 167 - 202 0040-1951 2006/03 [Refereed][Not invited]
     
    The dominantly passive volcanic Voring and More Margins, NE Atlantic, are separated by the 200km long Voring Transform Margin (VTM). The southern Voring Basin and the VTM have been studied by use of four regional Ocean Bottom Seismograph (OBS) profiles, combined by gravity modelling. The models demonstrate a complex pattern of magmatism along the transform margin. The distribution of magmatism seems to be related to the existence and trend of a lower crustal 8+ km/s body, interpreted as eclogitized rocks, present in the southern Voring Basin. Early Tertiary breakup related magmatic 'leakage' across the Continent-Ocean-Transition (COT) appears to be facilitated where this layer is absent. These results Support earlier workers who have concluded that the Jan Mayen Fracture Zone originated from a Caledonian zone of weakness. We propose that partly eclogitized rocks were uplifted into the lower crust close to this zone during the Caledonian orogeny and that this body acted as a barrier to magma emplacement during the Late Cretaceous-Early Eocene phase of rifting/breakup. The eclogitized terrain also appears to have caused northeastward channeling of the Late Cretaceous-Early Tertiary intrusions within the Voring Basin. An up to 10km thick pre-Cretaceous sedimentary basin in the southern Voring Basin may be genetically related to the NS-trending Late Paleozoic and Mesozoic rift basins in North-East Greenland. (c) 2006 Elsevier B.V. All rights reserved.
  • Rouzo, Stephane, Klingelhoefer, Frauke, Jonquet-Kolsto, Helene, Karpuz, Ridvan, Kravik, Karl, Mjelde, Rolf, Murai, Yoshio, Raum, Thomas, Shimamura, Hideki, Williamson, Paul, Geli, Louis
    Marine Geophysical Researches 27 (3) 181 - 199 0025-3235 2006 [Refereed][Not invited]
     
    This study presents the modelling of 2-D and 3-D wide-angle seismic data acquired on the complex, volcanic passive margin of the Voring Plateau, off Norway. Three wide-angle seismic profiles were shot and recorded simultaneously by 21 Ocean Bottom Seismometers, yielding a comprehensive 3-D data set, in addition to the three in-line profiles. Coincident multi-channel seismic profiles are used to better constrain the modelling, but the Mesozoic and deeper structures are poorly imaged due to the presence of flood basalts and sills. Velocity modelling reveals an unexpectedly large 30 km basement high hidden below the flood basalt. When interpreted as a 2-D structure, this basement high produces a modelled gravity anomaly in disagreement with the observed gravity. However, both the gravity and the seismic data suggest that the structure varies in all three directions. The modelling of the entire 3-D set of travel times leads to a coherent velocity structure that confirms the basement high; it also shows that the abrupt transition to the slower Cretaceous basin coincides in position and orientation with the fault system forming the Ran Ridge. The positive gravity anomaly over the Ran Ridge originates from the focussed and coincident elevation of the high velocity lower crust and pre-Cretaceous basement. Although the Moho is not constrained by the seismic data, the gravity modelled from the 3-D velocity model shows a better fit along the profiles. This study illustrates the interest of a 3-D acquisition of wide-angle seismic over complex structures and the benefit of the subsequent integrated interpretation of the seismic and gravity data.
  • AJ Breivik, R Mjelde, P Grogan, H Shimamura, Y Murai, Y Nishimura
    TECTONOPHYSICS 401 (1-2) 79 - 117 0040-1951 2005/05 [Refereed][Not invited]
     
    The western Barents Sea and the Svalbard archipelago share a common history of Caledonian basement formation and subsequent sedimentary deposition. Rock formations from the period are accessible to field study on Svalbard, but studies of the near offshore areas rely on seismic data and shallowdrilling. Offshore mapping is reliable down to the Permian sequence, but multichannel reflection seismic data do not give a coherent picture of older stratigraphy. A survey of 10 Ocean Bottom Seismometer profiles was collected around Svalbard in 1998. Results show a highly variable thickness of pre-Permian sedimentary strata, and a heterogeneous crystalline crust tied to candidates for continental sutures or major thrust zones. The data shown in this paper establish that the observed gravity in some parts of the platform can be directly related to velocity variations in the crystalline crust, but not necessarily to basement or Moho depth. The results from three new models are incorporated with a previously published profile, to produce depth-to-basement and -Moho maps south of Svalbard. There is a 14 km deep basement located approximately below the gently structured Upper Paleozoic Sorkapp Basin, bordered by a 7 km deep basement high to the west, and 7-9 km depths to the north. Continental Moho-depth range from 28 to 35 km, the thickest crust is found near the island of Hopen, and in a NNW trending narrow crustal root located between similar to 19 degrees E and 20 degrees E, the latter is interpreted as a relic of westward dipping Caledonian continental collision or major thrusting. There is also a basement high on this trend. Across this zone, there is an eastward increase in the V(P),V(P)/V(S) ratio, and density, indicating a change towards a more mafic average crustal composition. The northward basement/Moho trend projects onto the Billefjorden Fault Zone (BFZ) on Spitsbergen. The eastern side of the BFZ correlates closely with coincident linear positive gravity and magnetic anomalies on western Ny Friesland, apparently originating from an antiform with high-grade metamorphic Caledonian terrane. A double linear magnetic anomaly appears on the BFZ trend south of Spitsbergen, sub-parallel to and located 10-50 km west of the crustal root. Based on this correlation, it is proposed that the suture or major thrust zone seen south of Svalbard correlates to the BFZ. The preservation of the relationship between the crustal suture, the crustal root, and upper mantle reflectivity, challenges the large-offset, post-collision sinistral transcurrent movement on the BFZ and other trends proposed in the literature. In particular, neither the wide-angle seismic data, nor conventional deep seismic reflection data south of Svalbard show clear signs of major lateral offsets, as seen in similar data around the British Isles. (c) 2005 Elsevier B.V. All rights reserved.
  • MJELDE Rolf, RAUM Thomas, MYHREN Bjarte, MYHREN Bjarte, SHIMAMURA Hideki, MURAI Yoshio, TAKANAMI Tetsuo, KARPUZ Ridvan, NAESS Unni, NAESS Unni
    J Geophys Res 110 (B5) B05101 - 10.1029/2004JB003026 0148-0227 2005/05 [Refereed][Not invited]
     
    [1] A 180 km long seismic wide-angle profile was acquired across the Voring Plateau, NE Atlantic, using a tuned air gun array and three-component ocean bottom seismometers deployed with similar to 5 km spacing. The seismic P wave data have been modeled by ray tracing/inversion, and the model has been constrained by S wave and gravity modeling. The data and modeling have allowed us to depict the crustal structure and nature of the continent-ocean transition ( COT) in a classical volcanic margin case. The P wave velocity near the top of the main crustal layer is estimated to similar to 6.0 km/s landward of the similar to 25 km wide COT. The seaward increase to similar to 6.5 km/s in the COT is conformable with heavily intruded continental crust within this zone. Farther seaward, the velocity increase to similar to 6.9 km/s in the same layer suggests the presence of oceanic crust. The abundant magmatism landward of, and within, the COT is primarily observed as extrusives forming wedges of seaward dipping reflectors and mafic lower crustal intrusions/underplating. The maximum thickness of the oceanic crust is measured to similar to 23.5 km at the seaward termination of the COT.
  • Murai, Y, Akiyama, S, Yamashina, T, Kuwano, A
    北海道大学地球物理学研究報告 北海道大学大学院理学研究科地球惑星科学専攻(地球物理学) (68) 219 - 231 0439-3503 2005/03 [Not refereed][Not invited]
  • YAMADA TOMOAKI, KANEDA YOSHIYUKI, TAKANAMI TETSUO, MIKADA HITOSHI, SAKAI SHIN'ICHI, UEHIRA KENJI, MURAI YOSHIO, NISHINO MINORU, UHIRA KOICHI
    地震 57 (3) 281 - 290 0037-1114 2005/02 [Refereed][Not invited]
  • TAKANAMI TETSUO, MURAI YOSHIO, MACHIDA YUYA, SAITO ICHISUKE, MAKINO YUMI, KATSUMATA KEI, YAMAGUCHI TERUHIRO, NISHINO MINORU
    地震 日本地震学会 57 (3) 291 - 303 0037-1114 2005/02 [Refereed][Not invited]
  • 高波鉄夫, 村井芳夫, 町田祐弥, 斉藤市輔, 牧野由美, 勝俣啓, 山口照寛, 西野実
    月刊地球 (号外49) 80 - 91 0387-3498 2005/01 [Not refereed][Not invited]
  • R. Mjelde, T. Raum, A. Breivik, H. Shimamura, Y. Murai, T. Takanami, J. I. Faleide
    Petroleum Geology Conference Proceedings 6 803 - 813 2047-9921 2005 [Refereed][Not invited]
     
    Modelling of extensive seismic datasets recorded on Ocean Bottom Seismographs (OBS) on the outer Vøring Margin, NE Atlantic, has provided significant new insights into deeper sedimentary structures, distribution of sill-intrusions in the sedimentary section, top of the crystalline crust, the lower crust and Moho. Primarily based on the modelling of S-waves, it is concluded that the high-velocity lower crust most likely consists of a mixture of plume-related Late Cretaceous/Early Tertiary mafic intrusions mixed with older continental blocks. Northeastwards in the Vøring Basin, the landward limit of the lower crustal high-velocity layer steps gradually seawards, closely related to five crustal scale lineaments. Evidence for an interplay between active and passive rifting components is found on regional and local scales on the margin. The active component is evident through the decrease in magmatism with increased distance from the Iceland plume, and the passive component is illustrated by the fact that all resolved crustal lineaments to a certain degree acted as barriers to magma emplacement. A lithospheric delamination model is invoked to explain the observed variations in crustal velocities and thickness. The location of six Tertiary domal structures in the Vøring Basin is between, or in the vicinity of, pre-breakup high-velocity structures, which may act as rigid blocks during compression. It is proposed that the existence and trend of these high-velocity structures, subject to mild NW-SE compression, is the most important factor controlling the formation, spatial distribution and trend of the domes. Structures in the high-velocity lower crust may be the single most important element in controlling the formation of the domes all modelled highs in the lower crustal Early Tertiary intrusive layer seem to be related to the formation of domes. © 2005 Petroleum Geology Conferences Ltd.
  • F Ljones, A Kuwano, R Mjelde, A Breivik, H Shimamura, Y Murai, Y Nishimura
    TECTONOPHYSICS 378 (1-2) 17 - 41 0040-1951 2004/01 [Refereed][Not invited]
     
    The crustal structure along a 312 km transect, stretching from the axial mountains of the North Atlantic Knipovich Ridge to the continental shelf of Svalbard, has been obtained using seismic reflection data and wide angle OBS data. The resulting seismic V-p and V-s models are further constrained by a 2-D-gravity model. The principal objective of this study is to describe and resolve the physical and compositional properties of the crust in order to understand the processes and creation of oceanic crust in this extremely slow-spreading counterpart of the North Atlantic Ridge Systems. V, is estimated to be 3.50-6.05 km/s for the upper oceanic crust (oceanic layer 2), with a marked increase away from the ridge. The measured Vp of 6.55-6.95 km/s for oceanic layer 3A and 7.10-7.25 km/s for layer 3B, both with a V-p/V-s ratio of 1.81, except for slightly higher values at the ridge axis, does not allow a clear distinction between gabbro and mantle-derived peridotite (10-40% serpentized). The thickness of the oceanic crust varies a lot along the transect from the minimum of 5.6 km to a maximum of 8.1 km. The mean thickness of 6.7 km for the oceanic crust is well above the average thickness for slow-spreading ridges (< 10 mm/year half-spreading rate). The areas of increased thickness could be explained by large magma production-rates found in the zones of axial highs at the ridge axis, which also have generated the off-axial highs adjacent the ridge. We suggest that these axial and off-axial highs along the ridge control the lithological composition of the oceanic crust. This approach suggests normal gabbroic oceanic crust to be found in the areas bound by the active magma segments (the axial and off-axial highs) and mantle-derived peridotite outside these zone. (C) 2003 Elsevier B.V. All rights reserved.
  • M Shinohara, T Yamada, T Kanazawa, N Hirata, Y Kaneda, T Takanami, H Mikada, K Suyehiro, S Sakai, T Watanabe, K Uehira, Y Murai, N Takahashi, M Nishino, K Mochizuki, T Sato, E Araki, R Hino, K Uhira, H Shiobara, H Shimizu
    EARTH PLANETS AND SPACE 56 (3) 295 - 300 1343-8832 2004 [Refereed][Not invited]
     
    The Tokachi-Oki earthquake occurred on September 26, 2003. Precise aftershock distribution is important to understand the mechanism of this earthquake generation. To study the aftershock activity, we deployed forty-seven ocean bottom seismometers (OBSs) and two ocean bottom pressure meters (OBPs) at thirty-eight sites in the source region. We started the OBS observation four days after the mainshock for an observation period of approximately two months. In the middle of the observation period, nine OBSs near the epicenter of the mainshock were recovered to clarify the depth distribution of aftershocks near the mainshock. From the data overall OBS, seventy-four aftershocks were located with high spatial resolution. Most of the aftershocks were located in a depth range of 15-20 km and occurred within the subducting oceanic crust, the 5.5-km/s layer of the landward plate and the plate boundary. No aftershocks were found in the mantle of the subducting plate. The low seismic activity beneath the trench area where the water depth is greater than about 2000 m suggests a weak coupling between the two plates. The depth of the mainshock is inferred to be 15-20 km from the aftershock distribution.
  • AJ Breivik, R Mjelde, P Grogan, H Shimamura, Y Murai, Y Nishimura
    TECTONOPHYSICS 369 (1-2) 37 - 70 0040-1951 2003/07 [Refereed][Not invited]
     
    The Barents Sea is located in the northwestern corner of the Eurasian continent, where the crustal terrain was assembled in the Caledonian orogeny during Late Ordovician and Silurian times. The western Barents Sea margin developed primarily as a transform margin during the early Tertiary. In the northwestern part south of Svalbard, multichannel reflection seismic lines have poor resolution below the Permian sequence, and the early post-orogenic development is not well known here. In 1998, an ocean bottom seismometer (OBS) survey was collected southwest to southeast of the Svalbard archipelago. One profile was shot across the continental transform margin south of Svalbard, which is presented here. P-wave modeling of the OBS profile indicates a Caledonian suture in the continental basement south of Svalbard, also proposed previously based on a deep seismic reflection line coincident with the OBS profile. The suture zone is associated with a small crustal root and westward dipping mantle reflectivity., and it marks a boundary between two different crystalline basement terrains. The western terrain has low (6.2-6.45 km s(-1)) P-wave velocities, while the eastern has higher (6.3-6.9 kin s(-1)) velocities. Gravity modeling agrees with this, as an increased density is needed in the eastern block. The S-wave data predict a quartz-rich lithology compatible with felsic gneiss to granite within and west of the suture zone, and an intermediate lithological composition to the east. A geological model assuming westward dipping Caledonian subduction and collision can explain the missing lower crust in the western block by subduction erosion of the lower crust, as well as the observed structuring. Due to the transform margin setting, the tectonic thinning of the continental block during opening of the Norwegian-Greenland Sea is restricted to the outer 35 km of the continental block, and the continent-ocean boundary (COB) can be located to within 5 kin in our data. Distinct from the outer high commonly observed on transform margins, the upper part of the continental crust at the margin is dominated by two large, rotated down-faulted blocks with throws of 2-3 kin on each fault, apparently formed during the transform margin development. Analysis of the gravity field shows that these faults probably merge to one single fault to the south of our profile, and that the downfaulting dominates the whole margin segment from Spitsbergen to Bjornoya. South of Bjornoya, the faulting leaves the continental margin to terminate as a graben 75 km south of the island. Adjacent to the continental margin, there is no clear oceanic layer 2 seismic signature. However, the top basement velocity of 6.55 kin s(-1) is significantly lower than the high (7 km s(-1)) velocity reported earlier from expanding spread profiles (ESPs), and we interpret the velocity structure of the oceanic crust to be a result of a development induced by the 7-8-km-thick sedimentary overburden. (C) 2003 Elsevier Science B.V. All rights reserved.
  • Y Murai, S Akiyama, K Katsumata, T Takanami, T Yamashina, T Watanabe, Cho, I, M Tanaka, A Kuwano, N Wada, H Shimamura, Furuya, I, DP Zhao, R Sanda
    GEOPHYSICAL RESEARCH LETTERS 30 (9) 10.1029/2002GL016459  0094-8276 2003/05 [Refereed][Not invited]
     
    [1] The Kuril arc collides with the northeast Japan arc in the southern part of Hokkaido, Japan. 3-D tomographic inversion of data from a dense network of sensitive ocean-bottom seismographs and land stations has allowed imaging of previously unseen details of the arc-arc collision structure. A low velocity body dips gently southwestward, at depths of 35 to 45 km, from east of the Hidaka Mountains to the source area of the 1982 Urakawa-oki destructive earthquake (Ms 6.8). The low velocity body is the lower half of the lower crust of the Kuril arc, which must have been delaminated by the collision. We believe that the continuing collision of the delaminated lower crust with the northeast Japan arc resulted as an episode of aseismic slow slip prior to the 1982 Urakawa-oki earthquake as well being the reason for the high seismic activity in this region.
  • TAKANAMI TETSUO, MURAI YOSHIO, HONDA RYO, NISHIMURA YUICHI, KATSUMATA KEI, SHIMAMURA HIDEKI, HASEGAWA SEIZO, UKI NAGAHISA
    北海道大学地球物理学研究報告 北海道大学大学院理学研究科地球惑星科学専攻(地球物理学) (66) 63 - 75 0439-3503 2003/03/20 [Not refereed][Not invited]
  • AJ Breivik, R Mjelde, P Grogan, H Shimamura, Y Murai, Y Nishimura, A Kuwano
    TECTONOPHYSICS 355 (1-4) 67 - 97 0040-1951 2002/09 [Refereed][Not invited]
     
    The assembly of the crystalline basement of the western Barents Sea is related to the Caledonian orogeny during the Silurian. However, the development southeast of Svalbard is not well understood, as conventional seismic reflection data does not provide reliable mapping below the Permian sequence. A wide-angle seismic survey from 1998, conducted with ocean bottom seismometers in the northwestern Barents Sea, provides data that enables the identification and mapping of the depths to crystalline basement and Moho by ray tracing and inversion. The four profiles modeled show pre-Permian basins and highs with a configuration distinct from later Mesozoic structural elements, Several strong reflections from within the crystalline crust indicate an inhomogeneous basement terrain, Refractions from the top of the basement together with reflections from the Moho constrain the basement velocity to increase from 6.3 km s(-1) at the top to 6.6 km s(-1) at the base of the crust. On two profiles, the Moho deepens locally into root structures, which are associated with high top mantle velocities of 8.5 km s(-1). Combined P- and S-wave data indicate a mixed sand/clay/carbonate lithology for the sedimentary section, and a predominantly felsic to intermediate crystalline crust. In general, the top basement and Moho surfaces exhibit poor correlation with the observed gravity field, and the gravity models required high-density bodies in the basement and upper mantle to account for the positive gravity anomalies in the area. Comparisons with the Ural suture zone suggest that the Barents Sea data may be interpreted in terms of a proto-Caledonian subduction zone dipping to the southeast, with a crustal root representing remnant of the continental collision, and high mantle velocities and densities representing eclogitized oceanic crust. High-density bodies within the crystalline crust may be accreted island arc or oceanic terrain. The mapped trend of the suture resembles a previously published model of the Caledonian orogeny. This model postulates a separate branch extending into central parts of the Barents Sea coupled with the northerly trending Svalbard Caledonides, and a microcontinent consisting of Svalbard and northern parts of the Barents Sea independent of Laurentia and Baltica at the time. Later, compressional faulting within the suture zone apparently formed the Sentralbanken High. (C) 2002 Elsevier Science B.V. All rights reserved.
  • 村井芳夫, 山品匡史, 長郁夫, 田中昌之, 桑野亜佐子, 和田直人, 島村英紀, 古屋逸夫, 三田亮平
    月刊地球 海洋出版 24 (7) 495 - 498 0387-3498 2002/07 [Not refereed][Not invited]
  • Genshiro Kitagawa, Tetsuo Takanami, Asako Kuwano, Yoshio Murai, Hideki Shimamura
    Progress in Discovery Science 449 - 458 2002 [Not refereed]
  • R Mjelde, AJ Breivik, H Elstad, AE Ryseth, Skilbrei, JR, JG Opsal, H Shimamura, Y Murai, Y Nishimura
    NORWEGIAN JOURNAL OF GEOLOGY 82 (3) 183 - 202 0029-196X 2002 [Refereed][Not invited]
     
    The study integrates seismic wide-angle data from Ocean Bottom Seismographs (OBS), multichannel reflection seismics (MCS), and gravity- and aeromagnetic data, acquired along a 195 km profile from the Lofoten Basin, across the Sorvestsnaget Basin and the Veslemoy High. The study also includes OBS and MCS data from a 77 km long crossing profile, and available data from boreholes have been utilized. The MCS data generally provide the best imaging of the shallow to intermediate deep sedimentary levels (down to ca. 5 km), whereas the OBS- and gravity data are mainly used to constrain the geometries and velocity/density of the deeper sedimentary section, the crystalline crust and upper mantle layers. In the Lofoten Basin, an approximately 6 km wedge of Cenozoic sediments overlies ca. 6 km thick 'normal' oceanic crust. The continent-ocean-transition crossed by the profile occurs over a ca. 20 km wide zone, where a pronounced landward thickening of the crust is associated with the presence of igneous intrusions. Interpretational models with large amounts of salt in the Sorvestsnaget Basin have been tested, but none of the datasets used supported the presence of these. The modelling suggests that an interface at 7-8 km depth, which originally was interpreted as the base Cretaceous level, rather represents mid-Cretaceous levels in the Sorvestsnaget Basin, and that the base Cretaceous section is located as deep as 9-10 km here. A deeper interface (11-12 km) is interpreted as the mid-Jurassic level, based on similarities with observations in the Tromso Basin. The depth to the crystalline basement is estimated to be approximately 17 km which implies the presence of a ca. 5 km thick section of late Paleozoic to early Mesozoic sedimentary strata in the basin. Beneath the westernmost part of the Veslemoy High igneous intrusions are thought to be present at both upper and lower crystalline crustal levels, i.e. at about 8 and 15 km depth, respectively. Igneous rocks are also inferred to be present locally at ca. 6 km depth, beneath the base Cretaceous interface further east on the high. The depth to the Moho is estimated to be ca. 28 km beneath the Veslemoy High.
  • KATSUMATA KEI, WADA NAOTO, MURAI YOSHIO, TAKANAMI TETSUO, KASAHARA MINORU, SHIMAMURA HIDEKI
    月刊地球 (号外27号) 200 - 203 0387-3498 1999/12 [Not refereed][Not invited]
  • Y Murai, T Yamashita
    GEOPHYSICAL JOURNAL INTERNATIONAL 134 (3) 677 - 688 0956-540X 1998/09 [Refereed][Not invited]
     
    We have developed a new, simple method to calculate elastic waves scattered by zonally distributed cracks, taking account of multiple crack interactions. In this paper, our method is applied to SH-wave propagation in 2-D elastic media. An array of cracks is represented by an imperfectly bonded interface with a spatially inhomogeneous strength. A section of the interface with vanishing strength represents a stress-free crack, while a section with a significantly high elastic stiffness approximates an intact section. A cracked zone is therefore represented by a multilayered zone, partitioned by many such interfaces. In this formation, the interactions among all cracks can be considered implicitly and a mixed boundary-value problem can be circumvented. In our method, the elastic waves scattered by an interface are expanded into discrete reflected and transmitted plane-wave components. The wavefield in a cracked zone can be easily calculated by the wave propagator method or by the reflection and transmission operator method, on the basis of this expansion. The accuracy of our method is shown to be quite satisfactory for some examples. The zeroth-order reflection coefficient, which is defined as the coefficient for a reflected-wave component that has a reflection angle equal to the incidence angle, is specifically studied for some crack arrays, as an example of the application of our method of analysis.
  • Y MURAI, J KAWAHARA, T YAMASHITA
    GEOPHYSICAL JOURNAL INTERNATIONAL 122 (3) 925 - 937 0956-540X 1995/09 [Refereed][Not invited]
     
    We compute synthetic seismograms of SH waves that are multiply scattered by randomly distributed cracks. All the cracks are assumed to have the same length and strike direction; the crack surfaces are assumed to be stress-free, or to undergo viscous friction. We analyse the deterministic wave equation, and rigorously treat multiple crack interactions. We first calculate the wavefield in the wavenumber domain, and then we obtain the time-domain solution by its Fourier transform. A plane wave whose time dependence is described by the Ricker wavelet is assumed to be incident upon the region of crack distribution. The scattered waves are efficiently excited when the half-wavelength of the incident wave is close to or shorter than the crack length. High-wavenumber components are shown to be more abundant in the scattered waves when the crack distribution is denser. The time delay of the arrival of the primary wave, due to crack scattering, is shown to be prominent when the wavelength of the incident wave is much longer than the crack length. When the crack surfaces are subject to viscous friction, both the amplitudes of the scattered waves and the time delay of the primary-wave arrivals are smaller than those for the case of stress-free crack surfaces. When the crack distribution is statistically homogeneous, the calculated attenuation coefficient Q(-1) and phase velocity upsilon of the primary wave are generally consistent with those obtained by a stochastic analysis based on Foldy's approximation. A short analysis on the effect of inhomogeneous crack distribution shows that the wavenumber at which Q(-1) is at its peak value is smaller than that expected from the stochastic analysis for homogeneous crack distribution.
  • MURAI YOSHIO, KOBAYASHI YOSHIMASA
    地震 45 (2) 129 - 143 0037-1114 1992/09 [Refereed][Not invited]
  • Murai, Y, Matsunami, K, Kobayashi, Y
    京都大学防災研究所年報 京都大学防災研究所 (33B-1) 103 - 111 0386-412X 1990/04 [Not refereed][Not invited]
  • TSUTSUI Tomoki, KOBAYASHI Yoshimasa, SHIBA Yoshiaki, SUDA Yoshiyuki, MURAI Yoshio, IWATA Tomotaka, FUJIWARA Hiroyuki, MATSUI Ichiro
    Zisin (Journal of the Seismological Society of Japan. 2nd ser.) 公益社団法人 日本地震学会 42 (4) 405 - 418 0037-1114 1989 [Refereed][Not invited]
     
    Shallow subsurface structure around Hino River in the south-east shore of Lake Biwa was investigated by means of seismic reflection survey. The target area on the western margin of a seismic quiescent zone between the lake and the Suzuka range includes a belt of steep gradient of Bouguer anomaly which inclines west.
    Survey lines were designed, one on the east and another on the west of Hino River, in order to determine the subsurface structure beneath the steep gradient belt of gravity. Time sections were constructed by the Common Depth Point (CDP) method and followed by the frequency-wavenumber (F-K) migration to discuss the subsurface structure in detail.
    In the time sections of two survey lines, a horizon “A”, which is smoothly continuous, and another horizon “B”, which is inclined to south-west are recognized. The “A” horizon is considered to be a reflective plane in the Biwako-group, Quaternary deposits, by referring to the deep drilling log available in the neighborhood of the survey line. The “B” horizon is inferred to be an uncomformity boundary between the Biwako-group and the basement rocks because no continuous reflector can be found thereunder. Since the “B” horizon exibits a buried valley in the center of the western survey line and an south-westward inclination through the eastern survey line, this structure may be a remnant channel which used to run SE to NW during the erosive period of this area.
    It is noticeable that the “A” horizon shows onlap boundary relation against “B” horizon in the center of the eastern survey line. This suggests that this area turned into subsidence after the erosive period.
    No active faults were found above the “B” horizon through this survey. Hence low seismicity around the target area is likely inherent in character.

MISC

Books etc

  • Seismic waves, research and analysis
    MURAI Yoshio (Joint workWave propagation from a line source embedded in a fault zone)
    InTech 2012/01
  • Progress in discovery science
    Kitagawa, G, Takanami, T, Kuwano, A, Murai,Y, Shimamura, H (Joint workExtraction of signal from high dimensional time series: Analysis of ocean bottom seismograph data)
    Springer 2002
  • 震度を知る-基礎知識とその活用-
    気象庁監修 (Joint work現地調査等による震度の推定)
    株式会社ぎょうせい 1996/09

Presentations

  • Wave propagation from a line source embedded in a fault zone containing densely distributed parallel cracks  [Not invited]
    MURAI Yoshio
    Workshop on Recent developments in seismic wave scattering and heterogeneities in the Earth  2011/10
  • Estimate of the length and density of cracks distributed in a fault zone based on the results by DPEM  [Not invited]
    Murai Y, Yamashita T
    3rd SEGJ (The Society of Exploration Geophysicists of Japan)/SEG (The Society of Exploration Geophysicists) International Symposium on Geotomography -Fracture Imaging-  1995/11
  • Multiple scattering of SH waves in 2-D media with cracks  [Not invited]
    Murai Y, Kawahara J, Yamashita T
    2nd workshop on characterization of subsurface cracks  1992/07

Association Memberships

  • 日本地球惑星科学連合   米国地球物理学連合(American Geophysical Union)   日本地震学会   Japan Geoscience Union   American Geophysical Union   The Seismological Society of Japan   

Works

  • ノルウェー・トロムソ沖における地下構造探査
    村井 芳夫 2013/09 -2013/09
  • 村井 芳夫 2013/06 -2013/09
  • 北海道・十勝沖から根室沖における地下構造探査
    2010
  • 伊豆大島における海底地下構造探査
    2009
  • ノルウェー・モール海盆における地下構造探査
    2009
  • 茨城沖・房総沖における地下構造探査のための海底地震計組立
    2008
  • ノルウェー・ベア島沖における地下構造探査
    2008
  • 北海道積丹半島沖における自然地震観測
    2008
  • 2007年新潟県中越沖地震の余震観測
    2007
  • 北海道・えりも沖における地下構造探査
    2007
  • 茨城沖における地下構造探査
    2006
  • 十勝沖から根室沖における地下構造探査
    2006

Research Projects

  • 不均質媒質における波動伝播の物理学の深化-理論・実験・観測・シミュレーションからの複合的アプローチ-
    The Earthquake Research Institute, the University of Tokyo:cooperative research programs
    Date (from‐to) : 2012/04 -2013/03 
    Author : NAKAHARA Hisashi
  • 切迫度の高い震源域の先行調査観測
    文部科学省:海洋研究開発機構からの再委託事業
    Date (from‐to) : 2008/04 -2013/03
  • Promotion of studies on imaging and monitoring of small-scale heterogeneities in the lithosphere - Toward understanding the spatiotemporal variation of the structure of active regions in the crust
    The Earthquake Research Institute, the University of Tokyo:cooperative research programs
    Date (from‐to) : 2011/04 -2012/03 
    Author : MURAI Yoshio
  • 2011年東北地方太平洋沖地震に関する総合調査
    科研費補助金:特別研究促進費
    Date (from‐to) : 2011/04 -2012/03
  • 2011年東北地方太平洋沖地震に関する総合調査
    科研費補助金:特別研究促進費
    Date (from‐to) : 2011/03 -2011/03
  • Promotion of studies on imaging and monitoring of small-scale heterogeneities in the lithosphere -Toward understanding the spatiotemporal variation of the structure of seismogenic zones
    The Earthquake Research Institute, the University of Tokyo:cooperative research programs
    Date (from‐to) : 2010/04 -2011/03 
    Author : MURAI Yoshio
  • ノルウェー・モール海盆における地下構造探査
    ベルゲン大学(ノルウェー):奨学寄附金
    Date (from‐to) : 2009/04 -2010/03 
    Author : 村井 芳夫
  • 地震波による断層面の不均質性の解明に向けた理論的研究
    科研費補助金:基盤研究(C)
    Date (from‐to) : 2007/04 -2010/03 
    Author : 村井 芳夫
  • Japan Society for the Promotion of Science:Grants-in-Aid for Scientific Research
    Date (from‐to) : 2010 -2010 
    Author : SHINOHARA Masanao, MURAI Yoshio, FUJIMOTO Hiromi, HINO Ryota, SATO Toshinori, HIRATA Naoshi, OBARA Kazunari, SHIOBARA Hajime, IIO Yoshihisa, UEHIRA Kenji, MIYAMACHI Hiroki, KANEDA Yoshiyuki, KODAIRA Syuichi, MATSUZAWA Toru, OKADA Tomomi, YAGI Yuji, KOKETSU Kazuki, YAMANAKA Yoshiko, HIRAHARA Kazuro, TANIOKA Yuichiro, IMAMURA Fumihiko, SATAKE Kenji, TANAKA Jun, TAKAHASHI Tomoyuki, OKAMURA Makoto, YASUDA Susumu, KABEYASAWA Toshimi, HORI Muneo, HIRATA Kenji, TSUJI Yoshinobu, TAKAHASHI Yoshikazu, GOTO Hiroyuki, MORIKAWA Hitoshi
     
    The 2011 off the Pacific coast of Tohoku earthquake occurred offshore of northeast Japan region on March 11, 2011, and large ground motion with long duration and huge tsunami damaged human society. To clarify the generation mechanism of this large earthquake, we carried out various observations. Seismic observations in marine and land areas revealed variation of seismic activities in time and space. Geometry of the source was estimated by seismic surveys and ocean bottom crustal movement measurements. In addition, the slip distribution during the mainshock was obtained using various geophysical datasets. Field surveys revealed tsunamis and the generation process of the tsunamis was estimated. Damage of constructions was also clarified and useful information for disaster mitigation in future was obtained.
  • ノルウェー・ベア島沖における地下構造探査
    ベルゲン大学(ノルウェー):奨学寄附金
    Date (from‐to) : 2008/04 -2009/03 
    Author : 村井 芳夫
  • アスペリティ周辺の地震活動特性の研究
    文部科学省:東京大学地震研究所からの再委託事業
    Date (from‐to) : 2004/04 -2009/03
  • リソスフェアの短波長不均質構造の物理的解釈
    The Earthquake Research Institute, the University of Tokyo:cooperative research programs (B)
    Date (from‐to) : 2003/04 -2006/03
  • 日本学術振興会:科学研究費助成事業
    Date (from‐to) : 2003 -2004 
    Author : 西村 裕一, 島村 英紀, 高波 鐵夫, 村井 芳夫, 知北 和久
     
    前年度から今年度にかけて,水中という過酷な観測条件でも十分に使用でき,かつ測定機器や得られたデータの取り扱いが容易な湖底地震観測システムの構築を目指した.富士山の山頂での観測に使用された白山工業株式会社製のデータロガーMS7000を,センサーである地震計とともにガラス球に入れて密封する方式が最もよさそうであるとの結論に達した.最大の課題であった,ガラス球を介しての外部との通信や時計の校正については,いくつかの関係業者とともに検討し,ガラス球の内と外をモデムおよびイーサケーブルで繋ぐ方針でシステムを構築した.また,内蔵する電池もこのケーブルを介して充電する方式も可能とした.この設計図を元に,6ピンのケーブル4本がガラス球から外に出て来て,それぞれがモールドされている形で,ガラス球を加工して制作した.しかしながら,今年度は最終年度であるため機器開発を終え観測実験を開始する予定であったものの,紀伊半島沖地震津波(2004年9月5日)やスマトラ島沖地震津波(2004年12月26日)が発生し,研究代表者は緊急に結成された津波の調査チームに加わり,現地調査やデータ解析に従事することになった.緊急を要するこれらの不測の事態のため,機器の製作や観測実験の計画は大きく遅れてしまい,年度内完了が困難となり,結局は研究廃止申請を提出して終了することとした.開発がほぼ終わっている機器であるので,次年度以降,何らかの形で実験観測を行い,さらには実際の観測にも役立てるようにしたいと考えている.
  • 世界各地の特異な海溝-海嶺の深部地下構造と、微小地震から見た「活動のいま」の解明
    科研費補助金:基盤研究(A)
    Date (from‐to) : 2001/04 -2003/03
  • 日本学術振興会:科学研究費助成事業
    Date (from‐to) : 1999 -2000 
    Author : 村井 芳夫
     
    本研究では、断層破砕帯を平行な亀裂が密分布する領域としてモデル化し、そこでの波動伝播の理論計算を行い、断層破砕帯の速度構造、卓越する亀裂長と分布密度を推定する方法の確立を目的としている。今年度は、前年度までに開発された数値計算プログラムに、亀裂長を明確に定義できるよう修正を加え、亀裂長が散乱波の発生にどのように寄与しているか再検討した。その際、亀裂分布が断層トラップ波にどのように影響するか調べるために、震源と観測点は断層破砕帯の中心に置き、震源から射出される波は等方的であるとした。まず、断層破砕帯を亀裂が密に分布するだけの領域と仮定して計算したところ、震源から離れた観測点では、断層にトラップされた大振幅の長周期後続波をシミュレートすることはできなかった。そこで、亀裂の存在する領域が、周囲の媒質より低速度の帯状の領域として計算を行ったところ、直達S波の直後に大振幅の長周期後続波が現れた。したがって、実際の断層破砕帯は、断層面に平行な亀裂が非常に密に分布している領域というだけではなく、低速度の物質になっていると考えられる。次に、直達波の部分のスペクトルをとると、断層トラップ波に対応した長周期のピークが現れた。この長周期のピークは断層破砕帯の低速度層の速度と幅に対応している。一方、大振幅の長周期後続波より後のコーダ波の部分のスペクトルをとると、亀裂からの散乱波に対応した短周期のピークが現れた。この短周期のピークは、波数をk、亀裂長の半分をaとしたとき、亀裂間隔を変化させても、ka=1〜1.5で現れることから、亀裂長と入射波の半波長がほぼ等しい時に現れることがわかり、また、亀裂分布密度が高いほど、ピーク振幅が大きくなることがわかった。以上より、断層トラップ波を含む直達波の部分から断層破砕帯の速度構造が、コーダ波の部分から卓越する亀裂長と分布密度を推定できる可能性が得られた。
  • Japan Society for the Promotion of Science:Grants-in-Aid for Scientific Research
    Date (from‐to) : 1997 -1999 
    Author : SHIMAMURA Hideki, TAKANAMI Tetsuo, NISHIMURA Yuichi, TYO Ikuo, WATANABE Tomoki
     
    By the use of Ocean Bottom Seismographs (OBS), which was developed by the head investigator and his colleagues in Japan, we investigated the subbottom structure of the region from Mid Atlantic Ridge through Norwegian Atlantic coast with detail by means of controlled source seismology. Also, microearthquake survey by the use of OBS has revealed the detailed seismicity which is a direct results from the activity and movement of oceanic lithosphere which was born at the Mid Atlantic Ridge.
  • Study on structure of oceanic lithosphere
    Date (from‐to) : 1997
  • Study on tectonics of Northern Mid Atlantic Ridge
    International Joint Research Projects
    Date (from‐to) : 1997
  • Study on subsurface structure and seismicity at subduction zones
    Cooperative Research
    Date (from‐to) : 1997
  • Theoretical study on elastic wave propagation in a fault zone
    Date (from‐to) : 1995


Copyright © MEDIA FUSION Co.,Ltd. All rights reserved.