Subduction zones are places where one hydrated oceanic plate goes underneath another plate, and releases its fluids into the overlying mantle wedge. Slab-derived fluids play a key role in subduction zone processes. They serpentinize the cold forearc mantle at shallow depths; and deeper, they trigger hydrous mantle melting beneath the arc volcanoes and sometimes at backarc basin (BAB) spreading center. Examining the composition of arc and BAB magmas helps understanding genesis of subduction-related magmas, nature and composition of their mantle sources and slab-derived fluids. However, investigating such processes at shallow subduction zones is challenging, because the cold forearc mantle generally does not melt. Here, I investigate an unusual region in the southernmost Mariana convergent margin in the Western Pacific, near the Challenger Deep. The SE Mariana forearc stretched to accommodate opening of the southernmost Mariana Trough ~5Ma ago, opening the SE Mariana forearc rift (SEMFR) and causing seafloor spreading ~2.7-3.7Ma ago. The subducted slab beneath SEMFR deepens from <50km to ~100km, thus studying SEMFR lavas provides a unique opportunity to understand shallow subduction processes. By examining the major and trace element composition, the Pb-Nd-Sr isotopic ratios and the volatile contents (H2O, CO2, Cl, S, F) of SEMFR basalts, associated glassy rinds and olivine-hosted melt inclusions (Ol-MI) collected during three cruises (YK08-08, YK10-12, TN273), I show that: (i) SEMFR lavas were produced by adiabatic decompression melting of depleted asthenospheric BAB-like mantle at ~30±6.6 km depth and 1224±40oC; (ii) Ol-MI represent hydrous melts trapped by forearc mantle olivines. Xenocrysts were entrained with SEMFR basalts during ascent; (iii) SEMFR mantle flowed from the forearc towards the arc volcanoes and was metasomatized by shallow aqueous fluids; (iv) SEMFR shallow fluids are more aqueous than the fluids released beneath the Mariana arc and Mariana BAB; (v) the aqueous slab-derived fluids and the volatile fluxes are greatest at ~50-100km slab depth, suggesting that the minerals from the subducting plate mostly broke down beneath the arc to release their fluids. Such results provide new insights into shallow subduction processes, as previous studies showed that volatile fluxes and aqueous slab-derived fluids should increase toward the trench.

Subduction zone megathrusts produce the majority of the world's largest earthquakes. To understand the processes that control seismicity here, it is important to improve our knowledge on the subduction interface characteristics and its spatial variations. Nicoya Peninsula, Costa Rica, extends the continental landmass ~50 km towards the trench, making it a very suitable place to study interface activity from right on the top of the seismogenic zone of the Middle America Subduction Zone (MASZ). We contribute to and utilize an earthquake catalog of 8765 analyst-picked events to determine the spatial variability in the earthquake frequency-magnitude distribution (FMD) in this region. After initial detection, magnitude determination and location, the events are precisely relocated using a locally derived 3-D seismic compressional and shear wave velocity model (DeShon et al., 2006). After restricting the dataset to events nearest the interface and with low formal error (horizontal location error < 5 km), we retain a subset of 3226 events that best resolves interface activity. Beneath Nicoya, we determine the spatial variability and mean FMD of the interface, and focus on the relative relationship of small-to-large earthquakes, termed b-value. Across the region, the overall b-value (1.18 ± 0.04) is higher than the global average (b~1), and much larger than the global subduction zone average (b~0.6). Significant variation in b-value is observed along the active plate interface. A well resolved zone of lower b is observed at and offshore central Nicoya coast, in a previously determined locked patch using deformation observed from Global Positioning System (GPS). Conversely, high b-values prevail over the subducted portion of the Fisher ridge, which likely ruptured in the 1990 Gulf of Nicoya Mw 7.0 earthquake. Observed regions of low b-value approximately corresponds to more strongly-locked segments of the subduction interface resulting in higher differential stress, which may be released in the next large interface earthquake in this part of the MASZ. Across the region the b-value is found to vary inversely with the degree of interface locking. Thus, it is proposed that if sufficient data exist, spatial b-value mapping can be used as a proxy to determine interface locking. This method is especially useful along the subduction megathrust, which is generally offshore making geodetic measurements difficult.

3D seismic reflection data were acquired by the R/V Langseth and used to extract heat flow information using bottom-simulating reflector (BSR) depths across the southern Costa Rica convergent margin. These data are part of the CRISP Project, which will seismically image the Middle America subduction zone in 3D. The survey was conducted in an area approximately 55x11 km, northwest of the Osa Peninsula, Costa Rica. For the analysis presented here, seismic data were processed using a post-stack time migration. The BSR—a reverse polarity seismic reflection indicating the base of the gas hydrate phase boundary—is imaged clearly within the slope-cover sediments of the margin wedge. If pressure is taken into account, in deep water environments the BSR acts as a temperature gauge revealing subsurface temperatures across the margin. Two heat flow models were used in this analysis. In the Hornbach model BSR depth is predicted using a true 3D diffusive heat flow model combined with Integrated Ocean Drilling Program (IODP) thermal conductivity data and results are compared with actual BSR depth observations to constrain where heat flow anomalies exist. In the second model heat flow values are estimated using the heat flow equation. Uniform heat flow in the region should result in a deeper BSR downslope toward the trench due to higher pressure; however results indicate the BSR is deepest at over 325 meters below the seafloor (mbsf) further landward and shoals near the trench to less than 100 mbsf, suggesting elevated heat flow towards the toe of the accretionary prism. Heat flow values also reflect this relation. In addition to this survey-wide trend, local heat flow anomalies appear in the form of both circular patterns and linear trends extending across the survey, which can be related to mounds, thrust faults, folds, double BSRs, and seafloor erosion imaged in the seismic data. I suggest that these areas of higher local heat flow represent sites where advection of heat from deep, upward-migrating, thermogenically-sourced fluids and/or gases may be taking place. These heat flow trends have implications for not only earthquake nucleation, but also methane hydrate reserve stability. / text

The 2004 Sumatra-Andaman and the 2011 Tohoku-Oki earthquakes demonstrate the importance of understanding subduction zone earthquakes and the faults that produce them. Faults that produce earthquakes and/or tsunamis in these systems include plate boundary megathrusts, splay faults (out of sequence thrusts), and strike-slip faults from strain partitioning. Offshore Japan, IODP Exp. 314 collected logging while drilling (LWD) data across several seismically-imaged fault splays in the Nankai Trough accretionary prism. I combine LWD resistivity data with a model of fluid invasion to compare the permeabilities of sands. My results indicate that sands within faulted zones are 2-3 orders of magnitude more permeable than similar undisturbed sands. Therefore fault zones are likely to be fluid conduits within the accretionary wedge. Fluids can affect the physical and chemical properties of the faulted material, increasing pore pressures and effectively lubricating the faults. Fluids play an important role in fault slip, but hazard analysis also requires an understanding of fault geometry and slip direction. Both Japan and Sumatra exhibit strain partitioning, where oblique convergence between tectonic plates is partitioned between the megathrust and strike-slip faults proximal to the arc. Offshore Sumatra, I combine profiles from a 2D seismic survey (SUMUT) with previous bathymetry and active seismic surveys to characterize the West Andaman Fault adjacent to the Aceh forearc Basin. Along this fault I interpret transpressional flower structures that cut older thrust faults. These flower structures indicate that the modern West Andaman Fault is a right lateral strike-slip fault and thus helps to accommodate the translational component of strain in this highly oblique subduction zone. Offshore the Kii Peninsula, Japan, I analyze a trench-parallel depression that forms a notch in the seafloor just landward of the megasplay fault system, along the seaward edge of the forearc Kumano Basin. Using a 12 km wide, 3D seismic volume, I observe vertical faults and faults which dip toward the central axis of the depression, forming apparent flower structures. The along-strike geometry of the vertical faults makes predominantly normal or thrust motion unlikely. I conclude, therefore, that this linear depression is the bathymetric expression of a transtensional fault system. While the obliquity of convergence in the Nankai Trough is small (~15 degrees), this Kumano Basin Edge Fault Zone could be due to partitioning of the plate convergent strain. The location of the West Andaman Fault and KBEFZ within the forearc may be controlled by the rheology contrast between active accretionary wedges and the more stable crust beneath forearc basins. / text

No description available.

No description available.

This study applies Bayesian inversion to receiver functions (RF) to estimate local shear wave velocity (Vs) structure of the crust and upper mantle beneath two ocean bottom seismometers (OBS) offshore, and two land-based seismometers onshore Vancouver Island, British Columbia, Canada. We use passive seismic data recorded on NC89, a permanent NEPTUNE (North-east Pacific Time-series Undersea Networked Experiments) OBS located on the continental slope, and on a temporary autonomous KECK foundation OBS, KEBB, located at the Endeavour segment of the Juan de Fuca Ridge (JdFR). The two land based seismometers (OZB and PGC) are located on Vancouver Island and are part of the Canadian National Seismograph Network (CNSN). The introduction of NEPTUNE has helped to fill a gap in offshore seismic monitoring, however; due to high noise levels and a relatively short deployment time, few useful events have been recorded (to date) for RF analysis. In this study, we utilize three-component, broadband recordings of large (M6+), distant (30 -100 degrees) earthquakes to compute RFs due to locally generated P (compressional) to S (shear) converted waves. RFs are then inverted using a non-linear Bayesian approach which yields optimal profiles of Vs, Vp (compressional wave velocity), and strike and dip angles, as well as rigorous uncertainty estimates for these parameters. Near the JdFR a thin sediment layer (<1 km) is resolved overlying a 2 km thick oceanic crust. The crust contains a large velocity contrast at the depth of an expected axial magma chamber. The oceanic crust thickens to 10 km at the continental slope where it is overlain by 5 km of sediments. At the coastal station (OZB) a low velocity zone is imaged at 16 km depth dipping approximately 12 degrees NE. Evidence for this low velocity zone is also seen beneath southern Vancouver Island (PGC) at a depth consistent with previous studies. Determining such models at a number of locations (from the spreading ridge to the coast) provides new information regarding local structure and can aid in seismic hazard analysis. / Graduate

Communities on Vancouver Island, British Columbia face significant exposure to damaging earthquakes. This seismic risk arises not only from the Island’s proximity to crustal, sub-crustal and subduction earthquake sources in the Cascadia Subduction Zone and from their associated aftershock sequences, but also from environmental (natural and human-made) and social vulnerabilities in Vancouver Island communities and their current capacities to respond and recover from a large seismic event. Seeking to 1) assist community officials and the general public to better understand the scope of the earthquake risk on Vancouver Island; 2) raise awareness of the gaps in Vancouver Island’s risk assessment; 3) encourage and facilitate comprehensive seismic risk discussions at all levels of governance; and 4) offer quantitative data on which to base sound funding and policy decisions, this dissertation offers three new studies, presented in paper format, toward the comprehensive management of seismic risk on Vancouver Island. The first paper, reviews the components of risk and, building on international risk management standards and best practices, develops a new, comprehensive Disaster Risk Management (DRM) Framework for practitioners. This DRM Framework is then used to review existing knowledge of Vancouver Island’s seismic risk. A number of information gaps are identified, and two in particular, mainshock and aftershock hazard assessment, are targeted for further analysis. / Graduate

京都大学 / 新制・課程博士 / 博士(理学) / 甲第23713号 / 理博第4803号 / 新制||理||1687(附属図書館) / 京都大学大学院理学研究科地球惑星科学専攻 / (主査)准教授 河上 哲生, 教授 田上 高広, 教授 下林 典正 / 学位規則第4条第1項該当 / Doctor of Science / Kyoto University / DFAM

L'existence de glissements lents a été observé pour la première fois en Californie, sur une portion de la faille de San Andreas (Steinbrugge et al., 1960,Tocher 1960). Ils ont ensuite été détéctés dans les années 90s avec l'avènement des GPS. Les interactions entre glissement lents, ou glissements asismiques (qui n'émettent pas d'ondes élastiques) et les glissements sismiques est cependant mal connue.Pourtant, cette question est fondamentale puisque des glissements lents ont été observés avant plusieurs séismes.Par exemple, l’installation d’un glissement lent avant le séisme d’Izmit de 1999 en Turquie a été mise en évidence grâce à l’analyse de signaux répétitifs (Bouchon et al., 2011). De même, en étudiant le comportement des séismes répétitifs avant le méga-séisme de Tohoku de 2011, deux séquences de glissement lent en direction du point d’initiation ont été observées (Kato et al., 2012). D’autres glissements lents précédents des séismes de grandes ampleurs ont été reportés, notamment au Mexique, avant le séisme de Papanoa (Radiguet et al., 2016) de magnitude 7.3 et au Chili, avant le séisme d’Iquique de 2014 de magnitude 8.2 (Ruiz et al., 2014).L'objectif de ce travail de thèse est de mieux caractériser les interactions qui existent entre glissements sismiques et asismiques dans une zone de subduction très largement instrumentée: le Japon.La très grande densité des réseaux Japonais (sismique courte et longue période, GPS) autorise des seuils de détection des glissements sismiques ou lents très bas, et permet d'utiliser des méthodes de réseaux afin de maximiser le rapport signal sur bruit. La recherche d'épisodes de déformation lente est basée sur plusieurs types d'observables et de méthodes, et est complétée par l'analyse des changements de taux de sismicité accompagnant ces épisodes. Cette thèse vise ainsi à mieux comprendre comment du glissement lent peut -ou non- accélérer l'occurrence de grands séismes, et les conditions requises pour se faire. / Transient aseismic slip events (that do not emit elastic waves) were first discovered on the San Andreas fault in central California in 1960 (Steinbrugge et al., 1960, Tocher 1960) and were later confirmed by the development and installation of GPS stations. Aseismic slip can occur on continental faults as well as on subduction zones. However, the interactions between aseismic and seismic slips are not fully captured yet. Understanding the mechanisms at stake on fault planes is fundamental since several large earthquakes were preceded by aseismic slip episodes.For example, the setting of a slow slip event before the 1999 Izmit earthquake in Turkey was evidenced by the observation of repeating signals at one station (Bouchon et al., 2011). Also, by studying repeating earthquakes before the 2011 Tohoku earthquake, two slow slip episodes that migrated towards the rupture intitiation were observed (Kato et al., 2012). Other transient slips preceding large earthquakes occurred, including in Mexico, before the M7.3 Papanoa earthquake (Radiguet et al., 2016) and in Chile, before the 2014 Iquique earthquake (Ruiz et al., 2014).The aim of this thesis is to better characterize the interactions between seismic and aseismic slip that can occur in a subduction zone largely instrumented: Japan.The very high density of the japanese monitoring networks(seismic networks, both at short and long periods, and GPS network), allows the detection ofseismic and aseismic slip events with low intensity / size, and is amenable to the use of arraymethods to improve the signal-to-noise ratio. The search for episodes of transient deformation will be based on several observables andtechniques, and is complemented by the analysis of the changes in earthquake rates concomittent to these episodes. This work thus gives new clues on how aseismic deformation can -or not- accelerate the occurrence of strong earthquakes, and on what are the conditions (tectonic, dynamic) for this to happen.