SEARCH
Search Details
YOSHIOKA ShoichiResearch Center for Urban Safety and SecurityProfessor
Research activity information
■ Paper- Elsevier BV, Aug. 2025, Tectonophysics, 911, 230849 - 230849[Refereed]Scientific journal
- Tectonic plates bend and deform when approaching a subduction zone, creating intense faulting and highly variable stress and strain fields across short distances inside the slab. In September 2017, a large Mw8.2 intraplate normal fault occurred in southern Mexico, with an epicentral area located within a seismic gap where no megathrust had struck in more than a century. Despite the relatively young and hot Cocos plate, this seismic episode ruptured almost the entire slab below the brittle–ductile transition zone that normally limits the depth extent of such events. Here, we present a high-resolution thermomechanical model of spontaneous subduction for this area, where bending-induced brittle and ductile deformation and grain plate damage are considered. Modeling results show that the 2017 Mw8.2 Tehuantepec normal fault earthquake occurred due to the reactivation of one of the outer-rise-formed abyssal faults. In addition, the hypocenter was located in a stable, hydrated region of the lithospheric mantle at the transition limit between the elastic and ductile regimes. We found that earthquake rupture orientation is consistent with a region where a clear localized shear band of reduced effective viscosity is predicted. We propose that the rupture of this large intraslab event propagated in the ductile portion of the slab initially by a transformational faulting process, followed by a thermal runaway mechanism at greater depths and higher temperature.Frontiers Media SA, Jul. 2025, Frontiers in Earth Science, 13[Refereed]Scientific journal
- Springer Science and Business Media LLC, Apr. 2025, Scientific Reports, 15(1) (1)[Refereed]Scientific journal
- Abstract The Philippine Sea (PHS) plate embracing the Shikoku Basin subducts with a relatively young age varying between ~ 15 and 26 Ma. It exhibits a diversity of subduction styles along the Nankai Trough, from shallow beneath Shikoku and Chugoku regions to steep subduction beneath the Kii Peninsula. Although the onset of subduction along the entire Nankai Trough is synchronous at 15 Ma, present-day subduction geometries differ greatly between Shikoku and the Kii Peninsula. Seismic tomography beneath the Kii Peninsula reveals several high-velocity zones associated with the presence of a high density block. Some studies suggest that the location of this high density block coincides with the acidic rock group represented by the Kumano pluton, which is probably responsible for the high dip angle subduction of the Philippine Sea (PHS) slab beneath the Kii Peninsula. This study analyzes the effect of a high density and viscosity block placed in the vicinity of the Nankai Trough beneath the Kii Peninsula on the subduction dynamics of the PHS slab. We develop a series of 2D subduction numerical models using a visco-elasto-plastic rheology with spontaneous slab bending where we incorporate a high density and viscosity block in the proximity of the convergent margin. Our numerical results are compared with seismic tomography, present-day slab geometry, and seismicity distribution along an across-arc profile passing through the Kii Peninsula. The results show a good relationship between the slab geometry dynamics and high density and viscosity materials localized near the subduction zone. The high density and also possible high viscosity block identified in seismic tomography may hamper subduction of the PHS plate with a low dip angle beneath the Kii Peninsula, resulting in steep subduction.Springer Science and Business Media LLC, Apr. 2025, Progress in Earth and Planetary Science, 12(1) (1)[Refereed]Scientific journal
- Abstract We analyzed the spatiotemporal changes in the displacement and strain fields in the Niigata–Kobe tectonic zone (NKTZ), central Japan, with temporal resolutions of one and 3 years, respectively, using Global Navigation Satellite System (GNSS) time series data provided by the Geospatial Information Authority of Japan (GSI) over a long period of time from 1997 to 2022. We also attempted to extract steady-state variations without assuming linearity, which have never been obtained. The results show that the contraction rates of the dilatation rate, maximum shear strain rate and principal strain rate when postseismic deformations associated with large earthquakes were removed were consistently greater in the northern part of the NKTZ than in the central and southern parts of the NKTZ, reaching a maximum absolute value of 5.1 × 10–6 of contraction in the accumulated dilatation and an accumulated maximum shear strain of 3.0 × 10–6 for the 26 year period except for the years 2011 and 2012. These results are considered to result from the continuous expression of the elastic deformation of the thick sedimentary layer in the Echigo Plain and the viscous deformation of the lower crust, in addition to the effects of subduction of the oceanic Pacific plate. For the maximum shear strain, the central part of the NKTZ also has relatively large values, which may be caused by lateral deformation due to creep along the Atotsugawa and Ushikubi Fault Zones. Unlike before the Mw 9.0 2011 Tohoku-Oki earthquake, after the event, the dilatation rate was lower in the east‒west direction, and the variation in the dilatation rate appeared to return to the preearthquake dilatation rates both in terms of the spatial patterns and amplitudes over time within several years. The short-wavelength component of the dilatation rate in the east‒west direction in the northern part of the NKTZ shows that the spatial variations in the dilatation rate were similar before and after the 2011 Tohoku-Oki earthquake at 37.0° N, indicating continuous contraction. However, at 37.5° N, the spatial distribution of the dilatation rate changed from before to after the occurrence of the earthquake, and the characteristic large contractions that were identified before the earthquake were not observed after the event. Furthermore, Mj 4.0 or greater earthquakes were relatively uniformly distributed in the southern and central parts of the NKTZ. On the other hand, in the northern part of the NKTZ, there was a seismic gap, especially in areas with high strain concentrations in the Echigo Plain, where only large deformations could be caused by the low elasticity in the thick sedimentary layer. These results might indicate that a large earthquake will occur in the seismic gap in the future. In addition, a relatively large number of past major destructive earthquakes have occurred in the northern part of the NKTZ. In this area, since the Eastern Margin Fault Zone of the Takada Plain has a relatively high probability of an earthquake occurring in the next 30 years, additional detailed investigations of this area are needed.Springer Science and Business Media LLC, Mar. 2025, Progress in Earth and Planetary Science, 12(1) (1)[Refereed]Scientific journal
- Springer Science and Business Media LLC, Mar. 2025, Scientific Reports, 15(1) (1)[Refereed]Scientific journal
- Earth and Planetary Physics, 2025, Earth and Planetary Physics, 9(2) (2), 266 - 278[Refereed]Scientific journal
- Abstract The Philippine Sea (PHS) plate is currently subducting beneath the Kyushu region in southwestern Japan and exhibits a steep dip angle at depths of approximately 70–100 km. It has been suggested that the increase in the dip angle may be attributed to a hot plume, which was purportedly present in the northwestern part of Kyushu approximately 10–6 Ma, flowing southeastward toward southern Kyushu and causing steep shallow bending of the PHS plate. Using numerical simulations, this study quantitatively evaluates the impact of an upwelling hot plume from the northwestern part of Kyushu on the abrupt bending of the PHS plate at shallow depths. We use a time-dependent two-dimensional thermal convection model that incorporates the subduction history of the PHS plate along a vertical plane passing through a hot plume from northwestern Kyushu to the Nankai Trough. Numerical simulations were conducted to investigate the behavior of a slab when colliding with an upwelling hot plume, considering different values of slab viscosity. The results demonstrate that when a hot plume impacts a slab with a viscosity lower than that of a normal strong slab, the bending of the slab is reproduced at a depth of approximately 100 km, which is consistent with the observed PHS slab bending at depths of 70–100 km.Springer Science and Business Media LLC, Dec. 2024, Pure and Applied Geophysics, 182(4) (4), 1861 - 1888[Refereed]Scientific journal
- Elsevier BV, Aug. 2024, Environmental Modelling & Software, 179, 106129 - 106129[Refereed]Scientific journal
- Abstract Using horizontal and vertical GNSS time series data from the GSI in sotheastern Kyushu from January 1, 2017 to June 30, 2022, we detected a recent long-term slow slip event (L-SSE) that occurred in the Hyuga-nada region, southwest Japan, and estimated its spatiotemporal slip distribution. We performed such analysis considering the piecewise evolution of slip over time with time windows of 0.8 years to ensure a good signal-to-noise ratio in the horizontal displacements in each time window. The results showed that a slip of more than 2 cm occurred in the central part of Miyazaki Prefecture from 2018.5 to 2019.3 (the units were 0.1 years = 36.5 days). Then, the amount of slip increased and expanded slightly to the southern part of Miyazaki Prefecture from 2019.3 to 2020.1, and the amount of slip reached a maximum of 3.9 cm in the subsequent period (2020.1–2020.9). A smaller slip occurred at almost the same location in the following 0.8 years. Therefore, the duration of the L-SSE was approximately 3.2 years from approximately 2018.5 to 2021.7. The annual average maximum slip rate was approximately 4.9 cm/yr (3.9 cm/0.8 yr) during the period from 2020.1 to 2020.9. The maximum total slip was estimated to be approximately 12.9 cm, and the equivalent release moment was 4.9 × 1019 Nm, corresponding to Mw7.1. Compared to previous L-SSEs that occurred in the Hyuga-nada region, the annual average maximum slip rate was relatively low, and the main total slip was estimated at almost the same location on the plate interface. On the other hand, the slip duration of this L-SSE was the longest, and the release moment and moment magnitude were the greatest. The total slip area of more than approximately 10 cm in the Hyuga-nada L-SSE estimated in this study almost overlapped with the afterslip area of the December 3, 1996 Hyuga-nada earthquake, with a depth range of approximately 30–40 km. Short-term SSEs also occurred in the same depth range. These coincidental same depth ranges of interplate seismic events stem from large thermal gradients. The large thermal gradients play an important role in narrowing the depth range of frictional parameters (a–b) at the plate interface. In addition, the pore pressure and normal stress are closely related to the critical stiffness. In particular, contrasting values of pore pressure and normal stress at the shallow and deep sides of the mantle wedge corner may also contribute to transient aseismic slip within the same depth range.Springer Science and Business Media LLC, Apr. 2024, Progress in Earth and Planetary Science, 11(1) (1)[Refereed]Scientific journal
- Abstract In southern Chile, the Nazca plate is subducting beneath the South American plate. This region was struck by megathrust earthquakes in 1960 and 2010 and is characterized by the existence of a volcanic chain. In this region, we modeled a three-dimensional thermal structure associated with the subduction of the Nazca plate by using numerical simulations. Based on the obtained temperature distribution, we determined the updip and downdip limit temperatures for the region ruptured by these two megathrust earthquakes. In addition, the distributions of water content and dehydration gradient were calculated by using appropriate phase diagrams and compared with the location of the volcanic chain. As a result, we infer that the coseismic slip of the 2010 Mw8.8 Maule earthquake occurred only at temperatures lower than and around the 350 °C isotherm that resembles the beginning of the brittle‒ductile transition. We also deduce that the rupture of the 1960 Mw9.5 Valdivia earthquake propagated up to the 450 °C isotherm because the magnitude was considerably large and the young hot plate subducted near the Chile Ridge. In addition, the hydrous minerals in the turbidites, MORB and ultramafic rocks released fluids via dehydration reactions, and dehydrated water migrated upward almost vertically, decreasing the melting point of the mantle wedge and contributing to the formation of the volcanic chain.Springer Science and Business Media LLC, Jan. 2024, Geoscience Letters, 11(1) (1)[Refereed]Scientific journal
- Abstract We used Global Navigation Satellite System (GNSS) time series data to estimate the spatiotemporal slip distribution for a long-term slow slip event (L-SSE) that occurred in the Tokai region, central Japan, from 2012 to 2016. Since all the used GNSS data were affected by the postseismic deformation associated with the 2011 Mw9.0 Tohoku-Oki earthquake, we removed such postseismic signal from the time series of three components at each of the stations. The minimal time window for an inversion analysis was set to 0.5 years (6 months), taking into account the signal-to-noise ratio of displacements for each time window. In the horizontal displacement fields, displacements were observed in the south‒southeast and southeast directions on the west and east sides of Lake Hamana, respectively, with temporal changes in their amounts and directions. In the vertical displacement fields, uplift was observed on the east side of Lake Hamana. From these data, we estimated the L-SSE initiated in approximately 2012.5 and ended by 2017.0, indicating the duration time is 4.5 years and the duration was much longer than that obtained in a previous study. Using these data, we performed the inversion analysis, in which three a priori information were assumed, i.e., the spatial distribution of slip is smooth, slip mainly occurs in the direction of plate convergence, and the temporal variation in the slip is smooth, to obtain the spatiotemporal slip distribution on a plate boundary with 3-D geometry. As a result, we identified that the L-SSE consisted of two subevents. The first subevent initiated on the southwest side of Lake Hamana and expanded during the period from 2013.0 to 2014.5. The maximum slip velocity during the period from 2012.5 to 2017.0 was estimated to be approximately 3.5 cm/year there for 2013.5–2014.0. The second subevent took place on the west side of Lake Hamana gradually from 2015.0 to 2015.5, continued, and expanded from 2015.5 to 2016.5. From the cumulative slip distribution, we found that its shape spread in the dip direction and obtained a maximum slip of approximately 10.6 cm, a moment release of 2.7 × 1019 Nm, and an equivalent moment magnitude of 6.9. Comparing our results with the L-SSE that occurred in the Tokai region between 2000 and 2005, we found that the slip initiation location was almost the same, but the subsequent slip location was more southerly for the 2012–2016 Tokai L-SSE. Additionally, the maximum slip velocity and moment magnitude were smaller for the 2012–2016 L-SSE.Springer Science and Business Media LLC, Dec. 2023, Geoscience Letters, 10(1) (1)[Refereed]Scientific journal
- Abstract Compared to normal arc-related volcanic eruptions, the formation of a volcanic caldera is a relatively atypical event. During caldera formation a series of large volumes of magma are erupted, reducing the structural support for the rock above the magma chamber and creating a large depression at the surface called caldera. Los Humeros volcanic field (LHVF) represents one of the largest volcanic calderas in Mexico. It is located some 400 km from the trench at the eastern edge of the Trans Mexican Volcanic Belt where the depth to the Cocos slab is more than 300 km. In this study we employ high-resolution two-dimensional thermomechanical numerical simulations of magma intrusions and a horizontal tectonic strain rate to better understand the influence of crustal deformation for the formation of Los Humeros caldera. A minimum number of three thermal anomaly pulses of hydrated mantle material (with diameter of 15 km or more) and a regional strain rate of 7.927 × 10–16 s−1 are required for magma to reach the surface. Modeling results show that regional extension coupled with deep thermal anomalies (with a temperature excess of ΔT ≥ 100 °C) that come in a specific chain-type sequence produce surface deformation patterns similar to LHVF. We propose an asthenospheric sub-slab deep source (> 300 km depth) for the thermal anomalies where previous studies showed the existence of a gap or tear in the Cocos slab.Springer Science and Business Media LLC, Dec. 2023, Scientific Reports, 13(1) (1)[Refereed]Scientific journal
- Abstract Changing thermal regime is one of the key mechanisms driving seismogenic behaviors at cold megathrusts, but it is difficult to interpret warm subduction zones such as Vanuatu for the temperatures are higher than that accommodates shallow brittle failures. We construct a 3-D thermomechanical model to clarify the thermal structure that controls tectonic seismicity in Vanuatu and predict a warm circumstance associated with abundant seismicity. Results reveal a heterogeneous slab ranging from 300 °C to over 900 °C from the Moho to subvolcanic depth. The subduction seismicity corresponds well to the plate interface where dynamic thermal dehydration is focused. The transformation from hydrated basalts to eclogites along the slab facilitates the occurrence of intense earthquakes and slips. Multistage mineralogical metamorphism affects the dynamic stability of megathrusts and favors the generation of active interplate large events. Therefore, slab thermal dehydration plays a greater role than slab temperature condition in influencing the subduction earthquake distribution in warm subduction systems.Springer Science and Business Media LLC, Nov. 2023, Scientific Reports, 13(1) (1)[Refereed]Scientific journal
- Abstract Large subduction earthquakes induce complex postseismic deformation, primarily driven by afterslip and viscoelastic relaxation, in addition to interplate relocking processes. However, these signals are intricately intertwined, posing challenges in determining the timing and nature of relocking. Here, we use six years of continuous GNSS measurements (2015–2021) to study the spatiotemporal evolution of afterslip, seismicity and locking after the 2015 Illapel earthquake ($$M_w$$ 8.3). Afterslip is inverted from postseismic displacements corrected for nonlinear viscoelastic relaxation modeled using a power-law rheology, and the distribution of locking is obtained from the linear trend of GNSS stations. Our results show that afterslip is mainly concentrated in two zones surrounding the region of largest coseismic slip. The accumulated afterslip (corresponding to $$M_w$$ 7.8) exceeds 1.5 m, with aftershocks mainly occurring at the boundaries of the afterslip patches. Our results reveal that the region experiencing the largest coseismic slip undergoes rapid relocking, exhibiting the behavior of a persistent velocity weakening asperity, with no observed aftershocks or afterslip within this region during the observed period. The rapid relocking of this asperity may explain the almost regular recurrence time of earthquakes in this region, as similar events occurred in 1880 and 1943.Springer Science and Business Media LLC, Nov. 2023, Scientific Reports, 13(1) (1)[Refereed]Scientific journal
- Abstract The subduction of the Philippine Sea (PHS) plate along the Nankai Trough in in southwest Japan is a relatively recent process compared with subduction along the Japan Trench in northeast Japan. However, the tectonic evolution of the PHS plate along the Nankai Trough is still controversial and not fully understood. There are several competing hypotheses based on different estimates for the time variations of convergence rate and plate age. Our study employs numerical modelling of subduction in order to evaluate the slab evolution for the last 15 Myr and aims to evaluate each tectonic scenario against the present-day slab geometry along a profile passing through the Shikoku and Chugoku regions. The modelling strategy involves a parameter study where subduction initiation and various subduction parameters are analyzed in terms of subduction geometry evolution. Two-dimensional visco-elasto-plastic numerical simulations of spontaneous bending subduction predict that convergence rate and plate age variations play an important role in the evolution of subduction geometry. Modeling results after 15 Myr of evolution reveal that the tectonic model based on a high convergence rate between ~ 15 Ma and ~ 3 Ma produces a slab geometry that agrees well with the observed present-day slab shape specific for the Shikoku and Chugoku regions.Springer Science and Business Media LLC, Oct. 2023, Scientific Reports, 13(1) (1)[Refereed]Scientific journal
- Elsevier BV, Oct. 2023, iScience, 26(10) (10), 107936 - 107936[Refereed]Scientific journal
- Abstract The southern Chile subduction zone is a complex tectonic environment, where the Chile Ridge, the Nazca (NZ) and Antarctic (AN) plates subduct underneath the South American (SA) plate. The intersection between the NZ, AN and SA plates is referred to as the Chile Triple Junction (CTJ). In this region, a gap, often referred to as a slab window, has been formed between the NZ and AN slabs due to the divergence in their plate motion velocities, with volcanoes existing mainly above the subducted NZ and AN plates. In this study, we constructed a three-dimensional thermomechanical model associated with simultaneous subduction of the NZ and AN plates near the CTJ. The results show that the current temperature distributions on the upper surface of the slabs are higher closer to the Chile Ridge, and the AN plate has a distribution of elevated temperatures relative to the NZ plate at the same depth due to the northward migration of the CTJ and the slower convergence rate of the AN plate. Moreover, we calculated the water content and dehydration gradient from the temperature distribution near the upper surface of the slab and discussed their relationship to the distribution of volcanoes. In the northern part of the model domain, high dehydration gradients were obtained below the volcanic chain. Therefore, we suggest that the water released from the slab and the mantle wedge decreased the melting point of the mantle wedge just above the slab and produced melts, which may have contributed to form the overlying volcanoes.Springer Science and Business Media LLC, Jul. 2023, Geoscience Letters, 10(1) (1)[Refereed]Scientific journal
- Abstract In this study, we investigated the relationship between the strain accumulation before a slow slip event (SSE) and the strain release during the SSE for three recent SSEs along the Suruga Trough, Sagami Trough, and Nankai Trough, which are subduction zones in central to southwest Japan. The three analysed SSEs were the 2013–2016 Tokai long-term SSE (L-SSE), the 2018 Boso-Oki short-term SSE (S-SSE), and the 2019–2021 Central Shikoku L-SSE. We applied exponential and logarithmic functions to remove the postseismic deformations caused by the Mw 9.0 2011 Tohoku-Oki earthquake. We discovered a strong negative correlation between strain accumulation and strain release in the dilatation of all three SSEs and in the maximum shear strain of the 2018 Boso-Oki S-SSE. A comparison of the amount of strain accumulation with that of strain release revealed that approximately 30% of the strain was released in the 2013–2016 Tokai L-SSE, that 40% of the strain was released in the 2019–2021 Central Shikoku L-SSE, and that approximately 60% of the strain was released in the 2018 Boso-Oki S-SSE. This finding suggests that all of the accumulated strains are not necessarily released by the SSEs.Springer Science and Business Media LLC, Jun. 2023, Progress in Earth and Planetary Science, 10(1) (1)[Refereed]Scientific journal
- Abstract The geometry of the Philippine Sea slab (PHS) subducting beneath the Japanese islands has been imaged to 400 km depth beneath the Kyushu and Chugoku regions, whereas the PHS slab geometry beneath the Hokuriku region has only been determined to ~ 140 km depth, thereby indicating a large east–west asymmetry in the slab subduction. However, geologic evidence suggests that there was symmetrical east–west seafloor spreading along the axis of the Kinan seamount chain when the Shikoku basin was an active spreading center in the PHS plate. This inconsistency suggests that the PHS slab should be present beneath the Hokuriku region. Here we perform P-wave travel-time tomography across central Japan and conduct a two-dimensional plate subduction numerical simulation that reproduces the dual subduction of the PHS and Pacific (PAC) plates to elucidate the PHS slab geometry beneath central Japan. The tomography results reveal a high-velocity anomaly at ~ 150–250 km depth that extends from Wakasa bay to Noto peninsula and a slab window beneath the Hokuriku region. The numerical simulation results suggest that the PHS slab may have torn when it collided with the PAC slab, with the once leading edge of the PHS slab now present along the upper surface of the PAC slab beneath Noto peninsula. These results indicate that the PHS slab exists at ~ 250 km depth beneath the Hokuriku region, although it has been torn owing to its collision with the PAC slab, with this tear propagating westward to form a triangular slab window beneath the Hokuriku region. Graphical AbstractSpringer Science and Business Media LLC, May 2023, Earth, Planets and Space, 75(1) (1)[Refereed]Scientific journal
- Jan. 2023, scientific reports, 13, 1379, English[Refereed]
- Jan. 2023, Physics of the Earth and Planetary Interiors, 334, 106971, English[Refereed]
- Nov. 2022, Tectonophysics, 843, 229604, English[Refereed]
- Abstract We consider a Bayesian multi‐model fault slip estimation (BMMFSE), which incorporates many candidates of the underground structure (Earth structure and plate boundary geometry) model characterized by a prior probability density function (PDF). The technique is used to study long‐term slow slip events (L‐SSEs) that occurred beneath the Bungo Channel, southwest Japan, in around 2010 and 2018. We here focus on the two advantages of BMMFSE: First, it allows for estimating slip distribution without introducing relatively strong prior information such as smoothing constraints, by combining a fully Bayesian inference and better consideration of model uncertainty to avoid overfitting. Second, the posterior PDF for the underground structure is also obtained during the fault slip estimation, which can be used as priors for the estimation of slip distribution for recurring events. The estimated slip distribution obtained using BMMFSE agreed better with the distribution of deep tectonic tremors at the down‐dip side of the main rupture area than those based on stronger prior constraints when the corresponding Coulomb failure stress changes are compared. This finding suggests a mechanical relationship between the L‐SSE and the synchronized tremors. The use of the posterior PDF of the underground structure estimated for the 2010 L‐SSE as prior PDF for the 2018 event resulted in more consistent estimation with the data, indicated by a smaller value of an information criterion.American Geophysical Union (AGU), Aug. 2022, Journal of Geophysical Research: Solid Earth, 127(8) (8), e2021JB023712, English[Refereed]Scientific journal
- Jul. 2022, Tectonophysics, 838, 229479, English[Refereed]
- Jul. 2022, Pure and Applied Geophysics, 179, 4245 - 4265, English[Refereed]
- Abstract The Alaska subduction zone is characterized by a subducting oceanic plateau, which is referred to as the Yakutat terrane. Tectonic tremors occur in this zone, and there are few volcanoes above the subducted Yakutat terrane. In this study, we performed a 3-D numerical simulation of a thermal structure associated with the simultaneous subduction of the Yakutat terrane and Pacific plate to elucidate the mechanism of tectonic tremors, which typically involve the presence of water. We calculated the water content distribution near the slab surface by using the thermal structure obtained from our simulation and phase diagrams of the hydrous minerals included in the slab. As a result, dehydration from the marine sedimentary layer and oceanic crust was observed near the area where tectonic tremors occurred. Tectonic tremors occur only in the Yakutat terrane because the marine sedimentary layer and oceanic crust are thicker there, and the total amount of water content in these layers is higher; therefore, the amount of dehydration is also higher there than in the Pacific plate. Additionally, there are few volcanoes above the subducted Yakutat terrane because little water remains within the slab beneath the volcanic chain where magma is produced.Springer Science and Business Media LLC, Apr. 2022, Scientific Reports, 12(1) (1), 6234[Refereed]Scientific journal
- Jan. 2022, Chinese Journal of Geophysics (in Chinese), 65(1) (1), 1 - 11[Refereed]Scientific journal
- Abstract Long-term slow slip events (L-SSEs) have repeatedly occurred beneath the Bungo Channel in southwestern Japan with durations of several months to a couple of years, with a recurrence interval of approximately 6 years. We estimated the spatiotemporal slip distributions of the 2018–2019 Bungo Channel L-SSE by inverting processed GNSS time series data. This event was divided into two subevents, with the first on the southwest side of the Bungo Channel from 2018.3 to 2018.7 and the second beneath the Bungo Channel from 2018.8 to 2019.4. Tectonic tremors became active on the downdip side of the L-SSE occurrence region when large slow slips took place beneath the Bungo Channel. Compared with the previous Bungo Channel L-SSEs, this spatiotemporal slip pattern and amount were similar to those of the 2002–2004 L-SSE. However, the slip expanded in the northeast and southwest directions in the latter half of the second subevent. The maximum amount of slip, the maximum slip velocity, the total released seismic moment, and the moment magnitude of the 2018–2019 L-SSE were estimated to be 28 cm, 54 cm/year, $$4.4 \times 10^{19}$$ Nm, and 7.0, respectively, all of which were the largest among the 1996–1998, 2002–2004, 2009–2011, and 2018–2019 L-SSEs.Springer Science and Business Media LLC, Jan. 2022, Scientific Reports, 12(1) (1), 343[Refereed]Scientific journal
- Wiley, Nov. 2021, Terra Nova, 34(2) (2), 103 - 112[Refereed]Scientific journal
- Water circulation, along with plate subduction, is considered based on the stabilities of hydrous phases and pressure–temperature profiles of the sinking oceanic plate. Water in a rather hot slab like the present one may be largely liberated at shallow depths (< 150 km) and return to the ocean via. arc magmatism. On the other hand, stabilization of dense hydrous minerals under cooler conditions, which current subduction zones will soon experience, causes the transportation or reflux of seawater to the deep mantle, which reduces the total mass of surface seawater. Simple calculations accepting water contents in the subducting slab suggested by a recent seismic velocity structure model indicate that the Earth's oceans are likely to disappear ∼80 million years hence. Significant changes may happen such as the end of plate tectonics and the onset of snowball Earth, with associated catastrophes affecting life. The only way to confirm this picture of the future of the ocean planet Earth is to examine deep hydration taking place along the outer rise through direct analyses of the upper mantle across the Moho.Tokyo Geographical Society, Aug. 2021, Journal of Geography (Chigaku Zasshi), 130(4) (4), 585 - 597, Japanese[Refereed]Scientific journal
- Several interplate seismic events, such as short-term slow slip events (S-SSEs) and low-frequency earthquakes (LFEs), have been identified in the Ryukyu Trench, southwestern Japan. As one of the specific characteristics of this seismicity, the depths at which S-SSEs occur at the plate interface beneath Okinawa Island are approximately 5-10 km shallower than those beneath the Yaeyama Islands. To elucidate the cause of this difference in depth, we constructed a three-dimensional, Cartesian thermomechanical subduction model and applied the subduction history of the Philippine Sea (PHS) plate in the model region. As a result, the interplate temperatures at which S-SSEs take place were estimated to range from 350 to 450 °C beneath Okinawa Island and from 500 to 600 °C beneath the Yaeyama Islands. The former temperature range is consistent with previous thermal modelling studies for the occurrence of slow earthquakes, but the latter temperature range is by approximately 150 °C higher than the former. Therefore, explaining how the depth difference in S-SSEs could be caused from the aspect of only the thermal regime is difficult. Using phase diagrams for hydrous minerals in the oceanic crust and mantle wedge, we also estimated the water content distribution on and above the plate interface of the PHS plate. Near the S-SSE fault planes, almost the same amount of dehydration associated with phase transformations of hydrous minerals from blueschist to amphibolite and from amphibolite to amphibole eclogite within the oceanic crust were inferred along Okinawa Island and the Yaeyama Islands, respectively. On the other hand, the phase transformations within the mantle wedge were inferred only beneath the Yaeyama Islands, whereas no specific phase transformation was inferred beneath Okinawa Island around the S-SSE occurrence region. Therefore, we conclude that dehydrated fluid derived from the oceanic crust at the plate interface would play a key role in the occurrence of S-SSEs.May 2021, Scientific reports, 11(1) (1), 11251 - 11251, English, International magazine[Refereed]Scientific journal
- American Geophysical Union (AGU), Feb. 2021, Geophysical Research Letters, 48(3) (3)[Refereed]Scientific journal
- Springer Science and Business Media LLC, Dec. 2020, Geoscience Letters, 8(1) (1)
Abstract In repeating seismic event sequences within a specialized horizontal area, the moment magnitude is usually scaled with the recurrence interval. In addition to two horizontal dimensions, the vertical dimension plays a certain role in affecting the scaling law. However, whether and how the changing depth influences the scaling law remain enigmatic. Based on a large number of earthquake records with high-resolution epicenter locations in recent decades in Japan, we focus on a comparison between the 3-D seismic moment and seismic interval, which recognize the vertical dimension as the same dimension as the horizontal distances. The results show that (1) the seismic moment scaling law is applicable in the multiparameter 3-D models by visiting the 1.8 million events collected during a period of 15 years; (2) the vertical dimension of depth plays an important role in the Mo–SI relationship as well as in the variance in the 3-D seismic moment–interval magnitudes; and (3) the seismic moment rate, attributable to the plate convergence rate, varies with area and depth in influencing the regional earthquake recurrence frequency.[Refereed]Scientific journal - Springer Science and Business Media LLC, Dec. 2020, Scientific Reports, 10(1) (1), 15005 - 15005, English, International magazine
Abstract The SW Japan arc built by subduction of the Philippine Sea (PHS) plate exhibits uneven distribution of volcanoes: thirteen Quaternary composite volcanoes form in the western half of this arc, Kyushu Island, while only two in the eastern half, Chugoku district. Reconstruction of the PHS plate back to 14 Ma, together with examinations based on thermal structure models constrained by high-density heat flow data and a petrological model for dehydration reactions suggest that fluids are discharged actively at depths of 90–100 km in the hydrous layer at the top of the old (> 50 Ma), hence, cold lithosphere sinking beneath Kyushu Island. In contrast, the young (15–25 Ma) oceanic crust downgoing beneath Chugoku district releases fluids largely at shallower depths, i.e. beneath the non-volcanic forearc, to cause characteristic tectonic tremors and low-frequency earthquakes (LFEs) and be the source of specific brine springs. Much larger amounts of fluids supplied to the magma source region in the western SW Japan arc could build more densely-distributed volcanoes.[Refereed]Scientific journal - Sep. 2019, Journal of Geodynamics, 129, 299 - 312, EnglishSubduction thermal structure, metamorphism and seismicity beneath northcentral Chile[Refereed]Scientific journal
- Aug. 2019, Scientific Reports, 9, 11987, English[Refereed]Scientific journal
- 2019, Journal of Geophysical Research: Solid Earth, 124(7) (7), 6848 - 6865, English[Refereed]Scientific journal
- Mar. 2018, Journal of Geophysical Research: Solid Earth, 123(4) (4), 3080 - 3097, EnglishSubduction thermal regime, slab dehydration, and seismicity distribution beneath Hikurangi based on 3-D simulations[Refereed]Scientific journal
- Jan. 2018, Tectonophysics, 723, 288 - 296, English[Refereed]Scientific journal
- 2018, Journal of Geodynamics, EnglishSubduction thermal structure, metamorphism and seismicity beneath northcentral Chile[Refereed]Scientific journal
- 2018, Journal of Physics: Conference Series, 1036, 012018, English[Refereed]Scientific journal
- Dec. 2017, SCIENTIFIC REPORTS, 7, 16864, English[Refereed]Scientific journal
- Oct. 2017, TECTONICS, 36(10) (10), 1934 - 1946, English[Refereed]Scientific journal
- Sep. 2017, GEOPHYSICAL RESEARCH LETTERS, 44(18) (18), 9242 - 9252, English[Refereed]Scientific journal
- Jul. 2017, SCIENTIFIC REPORTS, 7, 6129, English[Refereed]Scientific journal
- Mar. 2017, GEOPHYSICAL RESEARCH LETTERS, 44(6) (6), 2679 - 2686, English[Refereed]Scientific journal
- Jan. 2017, JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH, 122(1) (1), 332 - 353, English[Refereed]Scientific journal
- Nov. 2016, PHYSICS OF THE EARTH AND PLANETARY INTERIORS, 260, 44 - 52, English[Refereed]Scientific journal
- Jun. 2016, JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH, 121(6) (6), 4458 - 4482, English[Refereed]Scientific journal
- May 2016, GEOPHYSICAL RESEARCH LETTERS, 43(10) (10), 4905 - 4912, English[Refereed]Scientific journal
- Jun. 2015, GEOPHYSICAL JOURNAL INTERNATIONAL, 201(3) (3), 1437 - 1455, English[Refereed]Scientific journal
- May 2015, GEOPHYSICAL JOURNAL INTERNATIONAL, 201(2) (2), 878 - 890, English[Refereed]Scientific journal
- Mar. 2015, PHYSICS OF THE EARTH AND PLANETARY INTERIORS, 240, 70 - 81, English[Refereed]Scientific journal
- Jan. 2015, GEOSCIENCE FRONTIERS, 6(1) (1), 61 - 78, English[Refereed]Scientific journal
- Springer New York LLC, Oct. 2014, Academic Psychiatry, 66(1) (1), English[Refereed]Scientific journal
- Oct. 2014, EARTH PLANETS AND SPACE, 66, English[Refereed]Scientific journal
- Oct. 2014, EARTH PLANETS AND SPACE, 66, English[Refereed]Scientific journal
- 2014, Earth, Planets and Space, 66(1) (1), English[Refereed]Scientific journal
- Nov. 2013, TECTONOPHYSICS, 608, 1094 - 1108, English[Refereed]Scientific journal
- Nov. 2013, JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH, 118(11) (11), 5838 - 5855, English[Refereed]Scientific journal
- Jul. 2013, TECTONOPHYSICS, 600, 165 - 174, English[Refereed]Scientific journal
- Jan. 2013, BULLETIN OF VOLCANOLOGY, 75(1) (1), 1 - 14, English[Refereed]Scientific journal
- Apr. 2012, PHYSICS OF THE EARTH AND PLANETARY INTERIORS, 196, 23 - 31, English[Refereed]Scientific journal
- The Volcanological Society of Japan, 2012, PROGRAMME AND ABSTRACTS THE VOLCANOLOGICAL SOCIETY OF JAPAN, 2012, 71 - 71, Japanese
- Feb. 2011, EARTH AND PLANETARY SCIENCE LETTERS, 303(1-2) (1-2), 1 - 10, English[Refereed]Scientific journal
- Nov. 2010, PHYSICS OF THE EARTH AND PLANETARY INTERIORS, 183(1-2) (1-2), 321 - 329, English[Refereed]Scientific journal
- Sep. 2009, EARTH AND PLANETARY SCIENCE LETTERS, 287(1-2) (1-2), 12 - 23, English[Refereed]Scientific journal
- Seismological Society of Japan, 2009, Zisin (Journal of the Seismological Society of Japan. 2nd ser.), 61(Supplement) (Supplement), 265 - 271[Refereed]Scientific journal
- Jul. 2008, EARTH AND PLANETARY SCIENCE LETTERS, 272(1-2) (1-2), 189 - 198, English[Refereed]Scientific journal
- Jun. 2008, TECTONOPHYSICS, 454(1-4) (1-4), 1 - 13, English[Refereed]Scientific journal
- Dec. 2007, TECTONOPHYSICS, 445(3-4) (3-4), 200 - 209, English[Refereed]Scientific journal
- Nov. 2007, EARTH AND PLANETARY SCIENCE LETTERS, 263(3-4) (3-4), 309 - 322, English[Refereed]Scientific journal
- Oct. 2007, GEOPHYSICAL JOURNAL INTERNATIONAL, 171(1) (1), 302 - 315, English[Refereed]Scientific journal
- Aug. 2025, 月刊技術士, 704, 4 - 7大地震と地震災害の発生メカニズム
- Mar. 2022, 東濃地震科学研究所報告, 47, 85 - 1022018〜2019年豊後水道長期的SSEのすべりの時空間分布の推定
- Mar. 2021, 東濃地震科学研究所報告, 46, 75 - 91津波波形とGNSSデータを用いた2015年イヤペル地震のすべり分布のインヴァージョン
- Mar. 2019, 都市安全研究センター 研究報告, 23, 1 - 17, Japanese想定海溝型巨大地震に伴う津波の数値シミュレーションReport research institution
- Mar. 2019, 東濃地震科学研究所報告, Seq. no.42, 37 - 53, Japanese関東下における海洋プレートの温度・脱水分布と微小地震の発生の関連性についてReport research institution
- Mar. 2019, 都市安全研究センター 研究報告, 23, 160 - 177, EnglishSurface Deformation Prior to the Imminent 2011 mw 9.0 Tohoku-Oki Earthquake and its Geodynamic ImplicationsReport research institution
- Mar. 2019, 都市安全研究センター 研究報告, 23, 149 - 159, EnglishRevisiting Inter-seismic Coupling in Southern Japan using Spatially Variable Slip Smoothing ConstraintsReport research institution
- Mar. 2018, 都市安全研究センター 研究報告, 22, 1 - 19, Japanese津波地震の津波数値シミュレーションReport research institution
- 地震予知総合研究振興会東濃地震科学研究所, Mar. 2017, 東濃地震科学研究所報告, (38) (38), 107 - 125, Japanese西南日本におけるフィリピン海プレートの沈み込みに伴う温度場・流れ場の3次元数値モデリング
- 神戸大学都市安全研究センター, Mar. 2017, 都市安全研究センター 研究報告, (21) (21), 1 - 16, JapaneseEstimation of the dynamic rupture parameters for the 2016 Tottoriken-chubu earthquakeReport research institution
- 神戸大学都市安全研究センター, Mar. 2017, 都市安全研究センター 研究報告, (21) (21), 184 - 194, EnglishSubduction of oceanic plate irregularities in south-central Mexico and the influence on subduction seismicityReport research institution
- 地震予知総合研究振興会東濃地震科学研究所, Mar. 2015, 東濃地震科学研究所報告, (34) (34), 61 - 77, Japanese太平洋プレートとフィリピン海プレートの沈み込みによる関東地方下の温度分布の数値シミュレーション
- 地震予知総合研究振興会東濃地震科学研究所, Mar. 2012, 東濃地震科学研究所報告, (29) (29), 53 - 68, Japaneseフィリピン海プレートの沈み込みに伴う温度分布の数値シミュレーション(続報)
- 地震予知総合研究振興会東濃地震科学研究所, Mar. 2011, 東濃地震科学研究所報告, (27) (27), 127 - 145, Japaneseフィリピン海プレートの沈み込みに伴う温度分布の数値シミュレーション
- 地震予知総合研究振興会東濃地震科学研究所, Mar. 2010, 東濃地震科学研究所報告, (25) (25), 11 - 20, Japanese九州東部下における深部低周波地震
- Jan. 2009, EARTH AND PLANETARY SCIENCE LETTERS, 277(1-2) (1-2), 292 - 292, EnglishOthers
- Mar. 2005, EARTH AND PLANETARY SCIENCE LETTERS, 231(3-4) (3-4), 239 - 247, English
- 2005, 地震, 58,121-141
- 2005, NONLINEAR PROCESSES IN GEOPHYSICS, 12(1) (1), 117 - 128, EnglishNonlinear multidimensional scaling and visualization of earthquake clusters over space, time and feature space
- Sep. 2004, PHYSICS OF THE EARTH AND PLANETARY INTERIORS, 146(3-4) (3-4), 513 - 530, English
- 2004, U.S. Geological Survey, Open File Report, 2004-1010, 15ppEffect of structural heterogeneity and slip distribution on coseismic vertical displacement from rupture on the Seattle fault
- Dec. 2002, JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH, 107(B12) (B12), English
- Dec. 2002, JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH, 107(B12) (B12), English
- 22 Nov. 2002, Tectonophysics, 359(1-2) (1-2), 171 - 187, English
- Nov. 2002, TECTONOPHYSICS, 359(1-2) (1-2), 171 - 187, English
- Oct. 2002, GEOPHYSICAL JOURNAL INTERNATIONAL, 151(1) (1), 69 - 82, English
- Oct. 2001, BULLETIN OF THE SEISMOLOGICAL SOCIETY OF AMERICA, 91(5) (5), 1182 - 1189, English
- Oct. 2001, Bulletin of the Seismological Society of America, 91(5) (5), 1182 - 1189, English
- Oct. 2001, PHYSICS OF THE EARTH AND PLANETARY INTERIORS, 125(1-4) (1-4), 1 - 17, English
- 2001, EARTH PLANETS AND SPACE, 53(4) (4), 261 - 273, English
- 28 Feb. 2000, Earth and Planetary Science Letters, 176(1) (1), 117 - 130, EnglishBook review
- Feb. 2000, EARTH AND PLANETARY SCIENCE LETTERS, 176(1) (1), 117 - 130, English
- 日本地震学会, 2000, Zisin : Journal of the Seismological Society of Japan, 52(4) (4), 469 - 485, Japanese
- 2000, 地震, 53, 25 - 35
- Sep. 1999, JOURNAL OF GEODYNAMICS, 28(2-3) (2-3), 175 - 192, English
- Aug. 1999, Physics of the Earth and Planetary Interiors, 115(1) (1), 17 - 34, English
- Aug. 1999, Physics of the Earth and Planetary Interiors, 115(1) (1), 17 - 34, English
- Apr. 1999, PHYSICS OF THE EARTH AND PLANETARY INTERIORS, 112(3-4) (3-4), 137 - 157, English
- Apr. 1999, PURE AND APPLIED GEOPHYSICS, 154(2) (2), 307 - 328, English
- Apr. 1999, PURE AND APPLIED GEOPHYSICS, 154(2) (2), 307 - 328, English
- Apr. 1999, PHYSICS OF THE EARTH AND PLANETARY INTERIORS, 112(3-4) (3-4), 137 - 157, English
- 海洋出版, 1999, 月刊地球, 21(4) (4), 197 - 204, JapaneseGPSデータを用いた西南日本におけるプレート間カップリングの推定
- 1999, 月刊地球, 25, 158 - 165GPSデータのインヴァージョン解析から推定した東北日本におけるプレート間カップリングの推定
- 09 Oct. 1998, 日本地震学会講演予稿集 = Programme and abstracts, the Seismological Society of Japan, 1998(2) (2), JapaneseNumerical simulation of temperature distribution and velocity anomalies in the horizontally lying slabs just above the lower mantle
- Oct. 1998, PURE AND APPLIED GEOPHYSICS, 152(3) (3), 443 - 464, English
- Oct. 1998, PURE AND APPLIED GEOPHYSICS, 152(3) (3), 443 - 464, English
- Dec. 1997, PHYSICS OF THE EARTH AND PLANETARY INTERIORS, 104(4) (4), 345 - 361, English
- Dec. 1997, PHYSICS OF THE EARTH AND PLANETARY INTERIORS, 104(4) (4), 331 - 344, English
- Dec. 1997, PHYSICS OF THE EARTH AND PLANETARY INTERIORS, 104(4) (4), 345 - 361, English
- Dec. 1997, PHYSICS OF THE EARTH AND PLANETARY INTERIORS, 104(4) (4), 331 - 344, English
- 15 Sep. 1997, 日本地震学会講演予稿集 = Programme and abstracts, the Seismological Society of Japan, 1997(2) (2), JapaneseVelocity anomalies and spatial distribution of physical properties in a horizontally lying slab in the Northwestern Pacific region
- 日本地震学会, 1997, 地震, 50(3) (3), 277 - 289, Japanese
- 1997, 月刊地球, 19, 699 - 703横たわるスラブの温度分布と速度異常
- Ehime University, 1997, Memoirs of the Faculty of Science, Ehime University, 3, 1-20, 1 - 20, EnglishStress analysis in the vicinity of a crack using finite element method
- 1997, 3, 1 - 20Stress analysis in the vicinity of a crack using finite element method
- Mar. 1996, NATURAL HAZARDS, 13(2) (2), 151 - 177, English
- Mar. 1996, NATURAL HAZARDS, 13(2) (2), 151 - 177, English
- 海洋出版, 1996, 月刊地球, 18(12) (12), 841 - 849, Japanese矩形断層面上のディスロケーションによる表面変位に及ぼす3次元不均質粘弾性構造の影響
- 1995, Island Arc, 4(2) (2), 89 - 103, English
- 1995, Journal of Geophysical Research, 100(B10) (B10), 20223 - 20244
- Jan. 1994, TECTONOPHYSICS, 229(3-4) (3-4), 181 - 200, EnglishINTERPLATE COUPLING IN THE KANTO DISTRICT, CENTRAL JAPAN, DEDUCED FROM GEODETIC DATA INVERSION AND ITS TECTONIC IMPLICATIONS
- Jun. 1993, GEOPHYSICAL JOURNAL INTERNATIONAL, 113(3) (3), 607 - 621, English
- Jul. 1992, PHYSICS OF THE EARTH AND PLANETARY INTERIORS, 72(1-2) (1-2), 83 - 98, English
- 1992, Geophysical Research Letters, 19, 313 - 316
- CRC Press, 2001Dictionary of Geophysics, Astrophysics, and Astronomy
- CRC Press, 2001Dictionary of Geophysics, Astrophysics, and Astronomy
- International Joint Workshop on Slow-to-Fast Earthquakes 2025, Sep. 20253D thermal structural and dehydration modeling in the southern Chile subduction zone and its relationship to interplate earthquakes and the volcanic chainOral presentation
- IAGA /IASPEI Joint Scientific Meeting 2025,, Sep. 20253-D thermal structure and dehydration near the Chile Triple Junction and its relation to slab window, tectonic tremors, and volcanoesOral presentation
- Japan Geoscience Union Meeting 2025, May 2025Three-dimensional thermomechanical modelling beneath Guerrero, Mexico, and its relation to the occurrence of interplate seismic phenomenaPoster presentation
- Japan Geoscience Union Meeting 2025, May 20253-D thermal structure and dehydration near the Chile Triple Junction and its relation to slab window, tectonic tremors, and volcanoesOral presentation
- SZNet Ocean Floor Observational Technologies Workshop, Jan. 2025Bayesian Selection of APG Sensor Locations off Taltal (∼25°S), Chile: Enhancing Imaging of Shallow Megathrust Slip BehaviorPoster presentation
- Slow-to-Fast Earthquake Workshop in Chile, Jan. 2025A MUltiscale Sparse Estimation (MUSE) approach for quasi-static slip inversionOral presentation
- Slow-to-Fast Earthquake Workshop in Chile, Jan. 20253D thermal structural and dehydration modeling in the southern Chile subduction zone and its relationship to interplate earthquakes and the volcanic chainOral presentation
- AGU fall meating 2024, Dec. 2024A MUltiscale Sparse Estimation (MUSE) Approach for Quasi-Static Slip InversionPoster presentation
- 日本地震学会2024年度秋季大会, Oct. 20242D collision models with visco-elasto-plastic rheology applied to WNW-ESE shortening along the Niigata-Kobe Tectonic zoneOral presentation
- 日本地震学会2024年度秋季大会, Oct. 2024最近日向灘で発生した大規模長期的スロースリップイベントの検出とその時空間滑り分布の推定Oral presentation
- International Joint Workshop on Slow-to-Fast Earthquakes 2024, Sep. 2024Detection of a recent large Hyuga-nada long-term slow slip event and estimation of its spatiotemporal slip distributionsPoster presentation
- 日本地球惑星科学連合2024年大会, May 2024WNW-ESE shortening in the Niigata-Kobe Tectonic zone explained by 2D visco-elasto-plastic collision numerical modelsPoster presentation
- 日本地球惑星科学連合2024年大会, May 2024Subduction thermal state, slab metamorphism and seismicity in Makran subduction zoneOral presentation
- International Joint Workshop on Slow-to-Fast Earthquakes 2024, Mexico, Feb. 2024Numerical modeling of subduction of the Philippine Sea plate tectonic history along the Nankai TroughPoster presentation
- 世界の沈み込み帯におけるSlowとFastの破壊現象の実像に関する研究集会:富山2024, Feb. 20242D visco-elasto-plastic subduction numerical models applied for the Kii Peninsula and Shikoku-Chugoku regionsOral presentation
- 日本地震学会2023年度秋季大会, Nov. 2023, 横浜, Domestic conferenceレシーバー関数を用いたオリビン準安定領域のイメージングPoster presentation
- 日本地震学会2023年度秋季大会, Oct. 2023, 横浜, Domestic conference九州下におけるホットプルームがフィリピン海スラブの屈曲に及ぼす影響 (3)Poster presentation
- 日本地震学会2023年度秋季大会, Oct. 2023, 横浜, Domestic conference新潟-神戸歪集中帯における歪場の時空間変化 -GNSSを用いた長期的解析-Poster presentation
- 日本地震学会2023年度秋季大会, Oct. 2023, 横浜, Domestic conference2D subduction models in the presence of a high-density rigid continental block: a case study for the slab geometry beneath the Kii PeninsulaOral presentation
- SF地震学B03班定例ミーティング, Aug. 2023, オンライン, Domestic conferenceアラスカ沈み込み帯におけるプレート境界地震イベントと3次元温度構造・脱水分布との関連性Oral presentation
- IUGG2023, Jul. 2023, Berlin, Domestic conference3-D thermal structure and dehydration in the southern Chile subduction zone and relation to interplate earthquakes and the volcanic chainOral presentation
- 日本地球惑星科学連合2023年大会, May 2023, 千葉Slab diversity and segmentation along the Middle America Trench: a framework for understanding geodynamic processes of subductionPoster presentation
- 日本地球惑星科学連合2023年大会, May 2023, 千葉, Domestic conferenceThermal regime and slab dehydration beneath the Izu-Bonin arc: Implications for fast and slow subduction earthquakeOral presentation
- 日本地球惑星科学連合2023年大会, May 2023, 千葉, Domestic conferenceMosaicking Crustal Surface Morphostructure and Seismic Cycle Deformation Patterns in Chile and Japan using GNSS Velocities and Machine LearningOral presentation
- 日本地球惑星科学連合2023年大会, May 2023, 千葉, Domestic conferenceNumerical modeling of subduction and evaluation of the Philippine Sea plate tectonic history along the Nankai TroughOral presentation
- 日本地球惑星科学連合2023年大会, May 2023, 千葉, Domestic conferenceRheological segmentation of the Cocos slab and its relation with the 2017 M8.2 Tehuantepec earthquakePoster presentation
- SF地震学A03班グループミーティング, Feb. 2023, English, Beppu, Domestic conferenceSubduction zone diversity along Middle America TrenchOral presentation
- SF地震学A03班グループミーティング, Feb. 2023, English, Beppu, Domestic conferenceNumerical modeling of subduction and evaluation of the Philippine Sea plate tectonic history along the Nankai TroughOral presentation
- 日本地震学会2022年度秋季大会, Oct. 2022, Japanese, 札幌, Domestic conference九州下におけるホットプルームがフィリピン海スラブの屈曲に及ぼす影響(2)Oral presentation
- 日本地震学会2022年度秋季大会, Oct. 2022, Japanese, 札幌, Domestic conferenceチリ沈み込み帯南部における3次元温度構造・脱水分布と海溝型巨大地震・通常地震・火山列との関連性Poster presentation
- 日本地球惑星科学連合大会2022, Jun. 2022, Japanese, オンライン・千葉, Domestic conference九州下におけるホットプルームがフィリピン海スラブの屈曲に及ぼす影響Poster presentation
- 日本地球惑星科学連合大会2022, May 2022, English, オンライン・千葉, Domestic conferenceIlluminating slow earthquakes: an approach from thermal modeling[Invited]Oral presentation
- 日本地球惑星科学連合大会2022, May 2022, English, オンライン・千葉, Domestic conferenceStrain accumulation and release processes associated with theoccurrence of SSEs in the Japanese IslandsOral presentation
- 日本地球惑星科学連合大会2022, May 2022, Japanese, オンライン・千葉, Domestic conference中部日本下深部のフィリピン海プレート形状とスラブウィンドウの形成Oral presentation
- 日本地球惑星科学連合大会2022, May 2022, Japanese, オンライン・千葉, Domestic conference地下構造の不確かさを考慮したベイズマルチモデル断層すべり推定について[Invited]Oral presentation
- 地殻変動研究のこれまでとこれから, Mar. 2022豊後水道下で発生する長期的スロースリップイベントについてOral presentation
- 都市地震リスク軽減に向けたデータ駆動型研究の最前線, Mar. 2022地下構造の不確かさを考慮したベイズマルチモデル断層すべり推定Oral presentation
- JST-JICA-UNAM Final Workshop, Mar. 20226-year achievement of the Japanese B1 groupOral presentation
- SF地震学B03班キックオフミーティング, Dec. 2021本プロジェクトにおける研究計画Oral presentation
- 学術変革領域研究A03班キックオフミーティング, Nov. 2021, オンライン今年度及び5カ年の活動計画Oral presentation
- 本地震学会2021年度秋季大会, Oct. 2021, オンラインアラスカ沈み込み帯における深部低周波微動と3次元温度構造・脱水分布との関連性Poster presentation
- 本地震学会2021年度秋季大会, Oct. 2021, オンライン北陸地方下のフィリピン海プレート形状を探るOral presentation
- 本地震学会2021年度秋季大会, Oct. 2021, オンライン豊後水道長期的スロースリップにおけるすべり分布のベイズマルチモデル推定:先験的拘束条件の推定結果に対する影響Oral presentation
- International Joint Workshop on Slow Earthquakes 2021, Sep. 2021, オンラインA Bayesian multi-model inference for fault slip distribution: the effect of prior constraints in the estimation for long-term slow slip events beneath the Bungo ChannelPoster presentation
- Japan Geoscience Union Meeting 2021, Jun. 2021, オンラインSpatiotemporal slip distributions of the 2018-2019 Bungo Channel long-term slow slip eventOral presentation
- 日本地球惑星科学連合2021年大会, Jun. 2021, オンラインアラスカ沈み込み帯における3次元温度構造モデリングOral presentation
- 日本地球惑星科学連合2021年大会, Jun. 2021, オンライン沈み込むスラブの年代の違いによる西南日本の火山分布の特徴~熱構造の観点から~[Invited]Invited oral presentation
- 日本地球惑星科学連合2021年大会, Jun. 2021, オンライン中部日本下に沈み込んだフィリピン海プレートの行方Oral presentation
- Online Workshop in February 2021 for “The Project for Hazard Assessment of Large Earthquakes and Tsunamis in the Mexican Pacific Coast for Disaster Mitigation”, Feb. 2021Three-dimensional thermal modeling associated with subduction of the Cocos plate beneath southwestern Mexico, and its relation to the occurrence of interplate seismic eventsOral presentation
- JpGU-AGU Joint Meeting 20203-D thermal modeling of generation mechanisms for short-term slowslip events and low-frequency earthquakes along the Ryukyu trenchPoster presentation
- JpGU-AGU Joint Meeting 2020, Jul. 20202011年東北地方太平洋沖地震前15年間のプレート間カップリングの時空間分布の推定Oral presentation
- AGU fall meating 2019, Dec. 2019, English, International conferenceFault stress inversion reveals seismogenic asperity of the 2011 Mw 9.0 Tohoku-Oki earthquakePoster presentation
- International Joint Workshop on Slow Earthquakes 2019, Sep. 2019, 仙台, International conferenceSpatial distribution of long-term slow slip event from 2018 to 2019 beneath the Bungo Channel under sparsity constraintsPoster presentation
- 日本地震学会2019年度秋季大会, Sep. 2019, 京都, Domestic conference琉球海溝におけるフィリピン海プレートの沈み込みに伴う3次元温度構造モデリングOral presentation
- 日本地震学会2019年度秋季大会, Sep. 2019, 京都, Domestic conference理論津波波形を用いた断層すべり分布のインヴァージョンの検討Poster presentation
- 日本地震学会2019年度秋季大会, Sep. 2019, 京都, Domestic conferenceGNSSデータを用いた2000年三宅島・神津島周辺の地震・火山活動に伴う地殻変動のモデル化Oral presentation
- 27th IUGG General Assembly, Jul. 2019, English, Montreal, International conferenceSubduction thermal structure, metamorphism and seismicity beneath northcentral ChileOral presentation
- Japan Geoscience Union Meeting 2019, May 2019, English, Makuhari Messe, Chiba, Domestic conferenceThree-dimensional thermal modeling associated with subduction of the Philippine Sea plate at the Ryukyu Trench, and its relation to the occurrence of interplate seismic eventsPoster presentation
- Japan Geoscience Union Meeting 2019, May 2019, English, Makuhari Messe, Chiba, Domestic conferenceSpatial distribution of slow slip events off the Boso peninsula from 1996 to 2018 under sparsity constraintsPoster presentation
- AGU Fall Meeting 2018, Dec. 2018, English, Washington, D.C., International conferenceSpatiotemporal locking state on the plate interface prior to the 2011 Tohoku-Oki earthquake inverted from GNSS dataPoster presentation
- AGU Fall Meeting 2018, Dec. 2018, English, Washington, D.C., International conferenceSeismogenesis of Dual Subduction Beneath Kanto, Central Japan Controlled by Fluid ReleaseOral presentation
- JST-JICA-UNAM Workshop for the Mid-term Evaluation, Nov. 2018, English, メキシコシティ, International conferenceThree-dimensional thermal modeling associated with subduction of the Cocos plate beneath southern Mexico, and its relation to the occurrence of interplate seismic eventsOral presentation
- 日本地震学会2018年度秋季大会, Oct. 2018, Japanese, 郡山, Domestic conference東北地方太平洋沖地震前 15 年間の地殻変動から推定したプレート間の固着・非地震性すべりの時空間分布Poster presentation
- 日本地震学会2018年度秋季大会, Oct. 2018, Japanese, 郡山, Domestic conference東北地方太平洋域地震の隣接地域における想定海溝型巨大地震に伴う津波の数値シミュレーションPoster presentation
- 10th ACES International Workshop, Sep. 2018, English, 南あわじ, International conferenceSpatiotemporal Distributions of Interplate Coupling and Aseismic Slips Prior to the 2011 Tohoku-Oki Earthquake Inferred From GNSS DataOral presentation
- 10th ACES International Workshop, Sep. 2018, English, 南あわじ, International conferenceSeismogenesis of Double Subduction Beneath Kanto, Central Japan Controlled by Fluid ReleaseOral presentation
- International Joint Workshop on Slow Earthquakes 2018, Sep. 2018, English, 福岡, International conferenceDiscussion based on spatial distribution of long-term slow slip events beneath the Bungo ChannelPoster presentation
- the AOGS 15th Annual Meeting, Jun. 2018, English, Honolulu, International conferenceGeodetic Data Inversion for Spatial Distribution of Long-term Slow Slip Events Beneath the Bungo Channel, Southwest Japan, Using Sparse Modelling,Poster presentation
- 日本地球惑星科学連合2018年大会, May 2018, English, 千葉, Domestic conferenceSpatial distribution of long-term slow slip events beneath the Bungo Channel under sparsity constraints (II)[Invited]Invited oral presentation
- 日本地球惑星科学連合2018年大会, May 2018, English, 千葉, Domestic conferencePostseismic displacement field associated with the 2011 Tohoku-oki earthquake by temperature dependent viscosity modelPoster presentation
- 日本地球惑星科学連合2018年大会, May 2018, Japanese, 千葉, Domestic conferenceGNSSデータを用いた日本列島の歪速度場の時空間変化Poster presentation
- 日本地球惑星科学連合2018年大会, May 2018, English, 千葉, Domestic conference3D Rheology effects on postseismic viscoelastic surface displacement fields in subduction zonesPoster presentation
- 研究集会「夢のある話を」, Mar. 2018, Japanese, 京都, Domestic conference環太平洋沈み込み帯における温度・脱水分布とプレート間地震の解明に向けてOral presentation
- 東京大学大気海洋研究所共同利用研究集会「地球環境と固体地球のかかわり」, Nov. 2017, Japanese, 柏, Domestic conferenceテクトニックな変動から推定した2011年東北地方太平洋沖地震前のプレートの固着状態と非地震性すべりの時空間分布Oral presentation
- スロー地震合同研究集会2017, Sep. 2017, English, 松山, Domestic conferenceSpatiotemporal distribution of interplate locking and aseismic slip prior to the 2011 Tohoku-Oki earthquakePoster presentation
- スロー地震合同研究集会2017, Sep. 2017, English, 松山, Domestic conferenceOn spatial distribution of dynamic rupture parameters for the 2016 Tottoriken-chubu earthquakePoster presentation
- International Meeting on “High-Dimensional Data-Driven Science” (HD3-2017), Sep. 2017, English, Kyoto, International conferenceGeodetic invention for spatial distridution of slow earthquakes under sparsity constraintsOral presentation
- IAG-IASPEI 2017, Aug. 2017, English, 神戸, International conferenceTwo-dimentional thermal modeling associated with subduction of the Philippine Sea Plate beneath southern Kyushu, JapanOral presentation
- IAG-IASPEI 2017, Aug. 2017, English, Kobe, International conferenceSpatiotemporal distribution of locking and aseismic slips prior to the 2011 Tohoku-oki earthquakePoster presentation
- IAG-IASPEI 2017, Aug. 2017, English, Kobe, International conferenceInterplate thermal regime and slab dehydration at the source region of episodic tremor and slow slip events in the Cascadia subduction zoneOral presentation
- IAG-IASPEI 2017, Aug. 2017, English, Kobe, International conferenceEstimation of the dynamic repture parameters for the 2016 Tottoriken-chubu earthquakePoster presentation
- JpGU-AGU Joint Meeting 2017, May 2017, Japanese, 千葉, International conferenceSpatial distribution of long-term slow slip events beneath the Bungo Channel under sparsity constraintsPoster presentation
- JpGU-AGU Joint Meeting 2017, May 2017, Japanese, 千葉, Domestic conferenceThermal regime and slab dehydration in the subducted Juan de Fuca plate beneath the Cascadia subduction zone based on 3D numerical simulationPoster presentation
- JpGU-AGU Joint Meeting 2017, May 2017, Japanese, 千葉, Domestic conference三次元数値シミュレーションによるカスかディア沈み込み帯における?ファンデフカプレートの温度構造とスラブ脱水についてPoster presentation
- JpGU-AGU Joint Meeting 2017, May 2017, English, 千葉, International conference九州南部の非火山地域における2次元熱対流モデリングPoster presentation
- JpGU-AGU Joint Meeting 2017, May 2017, Japanese, 千葉, Domestic conferenceテクトニックな変動から推定した東北地方太平洋沖地震前のプレート間の固着と非地震性すべりの時空間分布Poster presentation
- JpGU-AGU Joint Meeting 2017,, May 2017, English, 千葉, Domestic conferenceThermal regime and slab dehydration in the subducted Juan de Fuca plate beneath the Cascadia subduction zone based on 3D numerical simulation,Poster presentation
- JpGU-AGU Joint Meeting 2017, May 2017, English, 千葉, Domestic conferenceThermal regime and slab dehydration in the subducted Juan de Fuca plate beneath the Cascadia subduction zone based on 3D numerical simulation,Poster presentation
- JpGU-AGU Joint Meeting 2017, May 2017, English, 幕張メッセ, International conferenceSpatial distribution of long-term slow slip events beneath the Bungo Channel under sparsity constraintsPoster presentation
- JpGU-AGU Joint Meeting 2017, May 2017, English, 千葉, Domestic conference2-D thermal modeling along a non-volcanic region in southern Kyushu, Japan,Poster presentation
- JpGU-AGU Joint Meeting 2017, May 2017, English, 千葉, International conference2-D thermal modeling along a non-volcanic region in southern Kyushu, JapanPoster presentation
- JpGU-AGU Joint Meeting 2017, May 2017, Japanese, 千葉, Domestic conference2016年鳥取県中部地震の動的断層破壊シミュレーションPoster presentation
- AGU FALL MEETING 2016, Dec. 2016, English, San Francisco, International conferenceThermal regime, slab dehydration and seismicity distribution beneath the Hikurangi subduction zone based on 3D simulation,Poster presentation
- AGU FALL MEETING 2016, Dec. 2016, English, San Francisco, International conferenceNumerical simulations of temperature, dehydration, and flow fields associated with subduction of the Cocos plate, and its relation to the occurrence of interplate seismic events in southern Mexico,Poster presentation
- Annual Meeting 2016, Nov. 2016, English, Puerto Vallarta, Mexico, International conferenceNumerical simulations of temperature, dehydration, and flow fields associated with subduction of the cocos plate, and its relation to the occurrence of interplate seismic events in southern mexico,Oral presentation
- 日本地震学会秋季大会, Oct. 2016, English, 愛知, Domestic conferenceSlab dehydration, thermal regime, the distribution of tectonic tremors and seismicity beneath HikurangiOral presentation
- 35rd General Assembly of the European Seismological Commission, Sep. 2016, English, Trieste, Italy,, International conferenceThree-dimensional numerical modeling of temperature and dehydration fields associated with subduction of the Philippine Sea plate, southwest Japan,Oral presentation
- The 1st Asia-Pacific Workshop on Lithosphere and Mantle Dynamics, Sep. 2016, English, Taipei, International conferenceThree-dimensional numerical modeling of temperature and dehydration associated with subduction of the Philippine Sea plate, southwest Japan,[Invited]Invited oral presentation
- Joint Workshop on Slow Earthquakes 2016, Sep. 2016, English, 東京, Domestic conference3D thermal convection modeling associated with subduction of the Cocos plate in southern Mexico - Effects of trench retreat on temperature and flow fields -,Poster presentation
- 23rd Electromagnetic Induction Workshop, Aug. 2016, English, ChiangMai, Thailand,, International conferenceGeo-fluids distribution in mantle inferred from the electrical conductivity and simulated thermalOral presentation
- Crustal Dynamics 2016, Jul. 2016, English, 岐阜, Domestic conferenceSlab Dehydration, Interplate Hydrousity and Seismic Distribution in Modeling deformation, faulting, and mountain building in the island-arc crust of northeastern and central Japan considering heterogeneous thermal structurePoster presentation
- Crustal Dynamics 2016, Jul. 2016, English, 岐阜, Domestic conferenceMapping Crustal Structure Change in Tohoku With Time-Varying P-to-S Receiver FunctionsPoster presentation
- 日本地球惑星科学連合2016年大会, May 2016, Japanese, 千葉, Domestic conference九州地方の電気伝導度構造と温度構造から求められたマントルの流体分布Oral presentation
- 日本地球惑星科学連合2016年大会, May 2016, English, 千葉, Domestic conferenceThree-dimensional numerical modeling for subduction thermal regime, slab dehydration, and mantle flow beneath Kanto to Tohoku, Japan,Poster presentation
- 日本地球惑星科学連合2016年大会, May 2016, English, 千葉, Domestic conferenceSubduction of oceanic plate irregularities in Mexico and Japan and the influence on large megathrust earthquakes,Poster presentation
- AGU Chapman Conference on the Slow Slip Phenomena, Feb. 2016, English, Ixtapa, Mexico, International conferenceThree-dimensional thermal modeling associated with subduction of the Philippine Sea plate, southwest JapanPoster presentation
- AGU Chapman Conference on the Slow Slip Phenomena, Feb. 2016, English, Ixtapa, Mexico, International conferenceThree-dimensional thermal modeling associated with subduction of the Cocos plate beneath southern Mexico, and its relationship to the occurrence of interplate seismic eventsPoster presentation
- AGU Chapman Conference on the Slow Slip Phenomena, Feb. 2016, English, Ixtapa, Mexico, International conferenceSeismic distribution and dehydration of MORB associated with subduction of the Pacific and the Philippine Sea plates beneath the Tohoku and Kanto districts, JapanPoster presentation
- AGU Chapman Conference on the Slow Slip Phenomena, Feb. 2016, English, Ixtapa, Mexico, International conferenceA Study of Low-Frequency Earthquake Magnitudes in Northern Vancouver IslandPoster presentation
- AGU Fall Meeting 2015, Dec. 2015, English, San Francisco, International conferenceCrustal structure in the Tohoku region before and after the Tohoku-Oki earthquake with P to S receiver functionsPoster presentation
- AGU Fall Meeting 2015, Dec. 2015, English, San Francisco, International conferenceA Study of Low-Frequency Earthquake Magnitudes in Northern Vancouver IslandPoster presentation
- GeoBerlin 2015, Oct. 2015, English, Berlin, International conferenceImpact of olivine-spinel phase change kinetics on the deformation of the Mariana SlabOral presentation
- 地殻ダイナミクス -東北沖地震後の内陸変動の統一的理解-2015年研究集会, Sep. 2015, Japanese, 宮城, Domestic conference東北地方太平洋沖地震の地震前・地震時・地震後の固着状態とすべり分布の推定Oral presentation
- IUGG XXVI 2015,IUGG-1469, Jun. 2015, English, Prague, International conferenceTemperatures and fault slips on the upper surface of the subducting Philippine Sea plate beneath the Kanto district, central JapanOral presentation
- Central Asian tectonics and Western Pacific Geodynamics international workshop, Jun. 2015, English, Wuhan, International conference3-D thermal modeling associated with subduction of the Philippine Sea plate in southwest Japan.Oral presentation
- 東京大学地震研究所共同利用研究集会, May 2015, Japanese, 東京, Domestic conference2D and 3D thermal modeling associated with subduction of the Philippine Sea plate in central and southwest JapanOral presentation
- Tectonic Tremor and Silent Seismicity - 2015 International Workshop, Feb. 2015, English, Mexico City, International conferenceSpatiotemporal slip distributions of four long-term slow slip events beneath the Bungo Channel, southwest Japan, inferred from inversion analyses of GPS dataOral presentation
- AGU fall meeting, Dec. 2014, English, San Francisco, International conferenceImpact of Phase Change Kinetics on the Mariana Slab Within the Framework of 2-D Mantle ConvectionPoster presentation
- Joint Workshop on Slow earthquakes: "The prospects for studies of slow earthquakes toward Nankai Megaquake predict ions and disaster preventions", Sep. 2014, English, 京都, International conferenceRelations among temperature, dehydration of the PHS plate, and the three seismic events in the Tokai districtOral presentation
- AOGS 2014, Aug. 2014, English, 札幌, International conferenceRelationships among temperature, dehydration, and the occurrence of a megathrust earthquake, LFEs, and a SSE in the Tokai region, central JapanOral presentation
- IASPEI-LACSC, Jul. 2014, English, Bogota, Colombia, International conference3D thermal modeling associated with subduction of the Philippine Sea plate in southwest JapanOral presentation
- 日本地球惑星科学連合2014年大会, May 2014, Japanese, 横浜, Domestic conference東海地域におけるフィリピン海プレートの沈み込みに伴う温度場・脱水と想定東海地震・SSE・LFEとの関連性Poster presentation
- Geoscience Union Meeting 2014, Apr. 2014, English, 横浜, International conferenceThermal modeling associated with subduction of the Philippine Sea plate in southwest JapanOral presentation
- Geoscience Union Meeting 2014, Apr. 2014, English, 横浜, International conference3D numerical modeling of thermal regime and mantle flow associated with subduction of the two oceanic plates, JapanOral presentation
- 平成25年度研究集会, Feb. 2014, Japanese, 東京, Domestic conference豊後水道下で発生する長期的スロースリップイベントについてOral presentation
- 平成25年度研究集会, Feb. 2014, Japanese, 東京, Domestic conference東海地域におけるフィリピン海プレートの沈み込みに伴う温度場・脱水と想定東海地震・SSE・LFEとの関連性Poster presentation
- AGU fall meeting, Dec. 2013, English, San Francisco, International conferenceSlip distributions of three slow slip events beneath the Bungo Channel,southwest Japan,inferred from new inversion analyses of GPS dataPoster presentation
- 日本地震学会2013年度秋季大会, Oct. 2013, Japanese, 横浜, Domestic conference豊後水道下で発生した長期的スロースリップイベントについてOral presentation
- 日本地震学会2013年度秋季大会, Oct. 2013, Japanese, 横浜, Domestic conference東北地方の三次元P波減衰構造と島弧マグマ活動Oral presentation
- 日本地震学会2013年度秋季大会, Oct. 2013, Japanese, 横浜, Domestic conference東海地域におけるフィリピン海プレートの沈み込みに伴う温度場の数値シュミレーションPoster presentation
- 日本地震学会2013年度秋季大会, Oct. 2013, Japanese, 横浜, Domestic conference西日本におけるプレート境界面の摩擦係数と放射性元素の空間分布が地殻熱流量に与える影響Poster presentation
- 日本地震学会2013年度秋季大会, Oct. 2013, Japanese, 横浜, Domestic conference差分法を用いたマルチスケール動的断層破壊シュミレーションの試みPoster presentation
- IAHS-IAPSO-IASPEI, Jul. 2013, English, Gottenberg, International conferenceSlip distributions of three slow slip events beneath the Bungo Channel,southwest Japan,inferred from inversion analyses of GPS dataOral presentation
- AOGS Annual Meeting, Jun. 2013, English, Brisbane, International conferenceThree-dimensional numerical simulations of temperature, fluid flow and heat flow associated with subduction of the Philippine Sea plate beneath southwest Japan,Oral presentation
- 日本地球惑星科学連合2013年大会, May 2013, Japanese, 千葉, Domestic conference断層のステップ部における動的破壊の相互作用の数値シミュレーションPoster presentation
- 日本地球惑星科学連合2013年大会, May 2013, Japanese, 千葉, Domestic conference西南日本におけるフィリピン海プレートの沈み込みに伴う温度・流れ・地殻熱流量の3次元数値モデリングPoster presentation
- 日本地球惑星科学連合2013年大会, May 2013, Japanese, 千葉, Domestic conference曲面スラブの沈み込みに伴う温度場・流れ場・地殻熱流量の3次元数値モデリングPoster presentation
- 日本地球惑星科学連合2013年大会, May 2013, Japanese, 千葉, Domestic conference2012年12月7日三陸沖地震に伴う津波の数値シミュレーションPoster presentation
- AGU 20121 fall meeting, Dec. 2012, English, Americal Geophysical Union, San Francisco, USA, International conferenceStress perturbation given on the Mount Fuji Valcano Magma System caused by the Tohoku Megathrust Earthquake, JapanOral presentation
- AGU fall meeting, Dec. 2012, English, AGU, San Francisco, International conferenceAlong strike variation of tremor activities and thermal structures in various subduction zonesOral presentation
- 日本火山学会2012年度秋季大会, Oct. 2012, Japanese, 日本火山学会, 御代田, Domestic conference東北地方太平洋沖地震による富士山影響評価Oral presentation
- 日本地震学会2012年度秋季大会, Oct. 2012, Japanese, 日本地震学会, 函館, Domestic conference沈み込み帯における深部微動と温度の関係Oral presentation
- 日本地震学会2012年度秋季大会, Oct. 2012, Japanese, 日本地震学会, 函館, Domestic conference曲面スラブの沈み込に伴う温度場・流れ場・地殻熱流量の3次元数値シミュレーションPoster presentation
- 日本地震学会2012年度秋季大会, Oct. 2012, Japanese, 日本地震学会, 函館, Domestic conference関東地方下と東北地方下の温度構造から推定した太平洋プレート上面の温度プロファイルと含水量Poster presentation
- 日本地震学会2012年度秋季大会, Oct. 2012, Japanese, 日本地震学会, 函館, Domestic conference環太平洋地域の沈み込み帯におけるプレートの斜め沈み込みとその地域で発生する逆断層型地震のすべり方向のずれPoster presentation
- ECGS Workshop 2012, Oct. 2012, English, ECGS, Luxembourg, International conferenceTremor activities and thermal structures in various subduction zonesPoster presentation
- 日本地震学会2012年度秋季大会, Oct. 2012, Japanese, 日本地震学会, 函館, Domestic conference3次元スタガード格子差分法による自発的断層破壊と地震動のシミュレーションOral presentation
- The 9th Asian Seismological Commission General Assembly 2012, Sep. 2012, English, Asian Seismological Commission, Ulaanbaatar, International conferenceNumerical simulations of temperature fields associated with subduction of the oceanic plates beneath Tohoku and Kanto districtsOral presentation
- 日本地球惑星科学連合2012年大会, May 2012, Japanese, 日本地球惑星科学連合, 千葉, Domestic conference東北地方と関東地方におけるプレートの沈み込みに伴う温度分布の数値シミュレーションPoster presentation
- 日本地球惑星科学連合2012年大会, May 2012, Japanese, 日本地球惑星科学連合, 千葉, Domestic conference地球物理学と物質科学に基づくポスト・スピネル転移に対する制約:トンガスラブの例Oral presentation
- 日本地球惑星科学連合2012年大会, May 2012, Japanese, 日本地球惑星科学連合, 千葉, Domestic conferenceニュージーランド北島におけるテクトニックな微動活動と温度構造Poster presentation
- 日本地球惑星科学連合2012年大会, May 2012, Japanese, 日本地球惑星科学連合, 千葉, Domestic conferenceココスプレートの沈み込み方向と中央アメリカ海溝で発生する逆断層型地震のすべり角の方向のずれPoster presentation
- International Conference on a New Perspective of Great Earthquakes along Subduction Zones, Feb. 2012, English, 文部科学省科研費新学術領域研究会KANAME, 高知, International conferenceTowards comprehensive modeling of earthquake process in Nankai subduction zoneOral presentation
- International Conference on a New Perspective of Great Earthquakes along Subduction Zones, Feb. 2012, English, 文部科学省科研費新学術領域研究会KANAME, 高知, International conferenceNumerical simulations on temperature fields associated with subduction of the Philippine Sea plate in Southwest JapanOral presentation
- AGU 2011 fall meeting, Dec. 2011, English, Americal Geophysical Union, San Francisco, USA, International conference2-D numerical simulations on temperature fields associated with subduction of the Philippine Sea plate in Southwest JapanOral presentation
- 日本地震学会2011年秋季大会, Oct. 2011, Japanese, 日本地震学会, 静岡, Domestic conference豊後水道で発生する長期的スロースリップイベントについてOral presentation
- 日本地震学会2011年秋季大会, Oct. 2011, Japanese, 日本地震学会, 静岡, Domestic conference過去6000年間の三陸沖超巨大古津波履歴を示す二つの露頭Oral presentation
- 日本地震学会2011年秋季大会, Oct. 2011, Japanese, 日本地震学会, 静岡, Domestic conference伊勢湾付近の深部低周波地震の空白域についてOral presentation
- 日本地震学会2011年秋季大会, Oct. 2011, Japanese, 日本地震学会, 静岡, Domestic conferenceニュージーランド南島で発生した2010年Darfield地震と2011年クライストチャーチ地震による被害Oral presentation
- Asia Oceania Geosciences Society 8th Annual Meeting 2011, Aug. 2011, English, Asia Oceania Geosciences Society, Taipei, Taiwan, International conference2-D Numerical Simulations on Temperature Fields Associated with Subduction of the Philippine Sea Plate in Southwest JapanOral presentation
- 日本地球惑星科学連合2011年大会, May 2011, Japanese, 日本地球惑星科学連合, 千葉, Domestic conference津波はリアス式海岸で増幅したのか? -現地調査による湾内津波の挙動-Oral presentation
- 日本地球惑星科学連合2011年大会, May 2011, Japanese, 日本地球惑星科学連合, 千葉, Domestic conference太平洋プレートとフィリピン海プレートの沈み込みによる関東下の温度分布の数値シミュレーションOral presentation
- 日本地球惑星科学連合2011年大会, May 2011, Japanese, 日本地球惑星科学連合, 千葉, Domestic conference下部マントルへのスタグナントスラブ崩落の数値シミュレーションOral presentation
- 日本地球惑星科学連合2011年大会, May 2011, Japanese, 日本地球惑星科学連合, 千葉, Domestic conferenceいかにして巨大津波から避難したか? ~陸前高田市におけるインタビュー調査からの考察(序報)~Oral presentation
- 日本地球惑星科学連合2011年大会, May 2011, Japanese, 日本地球惑星科学連合, 千葉, Domestic conference2011年3月15日静岡県東部で発生した地震(Mjma6.4)による富士山マグマ溜まりへの影響評価Oral presentation
- EGU2011, Apr. 2011, English, European Geosciences Union, Vienna, Austria, International conferenceNumerical simulations on fall of stagnant slabs into the lower mantleOral presentation
- 第890回地震研究所談話会, Jan. 2011, Japanese, 東京大学地震研究所, 東京, Domestic conferenceDirect imaging of the metastable olivine wedge for a deep dry cold slab beneath southwest JapanOral presentation
- AGU 2010 fall meeting, Dec. 2010, English, Americal Geophysical Union, San Francisco, USA, International conferenceEffects of trench migration on fall of stagnant slabs into the lower mantleOral presentation
- The 8th ASC General Assembly 2010, Nov. 2010, English, Asian Seismological Commission, Hanoi, Vietnam, International conference2-D thermal modeling of the Philippine Sea plate in southwest Japan -effects of hot plume, yield stress, and surface erosion-Oral presentation
- The 8th ASC General Assembly 2010, Nov. 2010, English, Asian Seismological Commission, Hanoi, Vietnam, International conference2D temperature model in deep slabs incorporating kinetics and latent heat release associated with the phase transformationsOral presentation
- 日本地震学会2010年秋季大会, Oct. 2010, Japanese, 日本地震学会, 広島, Domestic conference西南日本におけるフィリピン海プレートの沈み込みの2次元温度モデリング~ホットプリューム、降伏応力、削剥の効果~Oral presentation
- 日本地震学会2010年秋季大会, Oct. 2010, Japanese, 日本地震学会, 広島, Domestic conferenceGPS データを用いた西南日本のすべり欠損分布の推定~プレート形状,観測点配置の影響~Oral presentation
- The Meeting of the Americas, Aug. 2010, English, Americal Geophysical Union, Foz do Iguassu, Brasil, International conference2D temperature model in deep slabs incorporating kinetics of the 410-km and 660-km phase transformationsOral presentation
- 日本地球惑星科学連合2010年大会, May 2010, Japanese, 日本地球惑星科学連合, 千葉, Domestic conference西南日本におけるフィリピン海プレートの沈み込みに伴う温度分布の数値シミュレーションOral presentation
- EGU2010, May 2010, English, European Geosciences Union, Vienna, Austria, International conferenceNumerical simulations of temperature distributions associated with subduction of the Philippine Sea plate, southwest Japan,Oral presentation
■ Works
- 1998 - 2003プレート境界・内陸活断層のモデル化及び3次元粘弾性媒質中 での応力伝播の並列シミュレーション・モデル化に関する研究
- 1996 - 2001南海トラフにおける海溝型巨大地震災害軽減のための地震発生機構のモデル化・観測シス テム高度化に関する総合研究
- Japan Society for the Promotion of Science, Grants-in-Aid for Scientific Research, Grant-in-Aid for Transformative Research Areas (A), Kyoto University, 10 Sep. 2021 - 31 Mar. 2026Slow-to-Fast earthquakes through comparison across global subduction zones
- Japan Society for the Promotion of Science, Grants-in-Aid for Scientific Research, Grant-in-Aid for Scientific Research (B), Kobe University, 01 Apr. 2021 - 31 Mar. 2025Why do intra-slab earthquakes occur ? Relationship with the three-dimensional distribution of olivine metastable phase, temperature, and dehydration研究の目的と研究実施計画は以下の通りである。 ・研究の目的:深さ200~500 kmで発生するスラブ内地震に関して、以下の5項目を本研究の目的とし、吉岡と中島で分担・協力して、研究を遂行する。 ・研究実施計画:① 地震発生メカニズムの包括的理解に不可欠な同一のスラブにおける地震の震源パラメターの系統的な深さ変化の解明(中島) ② 太平洋スラブ内のオリビン準安定相領域の3次元空間分布の特定(中島・吉岡) ③ オリビン準安定相の存在を考慮したスラブ深部の3次元温度分布を計算するコードの開発(吉岡) ④ 太平洋スラブ深部の3次元温度分布・含水量分布の推定(吉岡) ⑤脱水脆性化説、相転移断層説、熱的不安定説の特定(吉岡・中島) 上記①、②に関しては、昨年度までにほぼ実施が完了した。 吉岡は、上記③に関して、2次元鉛直断面内で海洋プレートの沈み込みに伴うオリビンからウォズレアイトの相転移のカイネティクスを取り入れたプログラムのチューニングを行い、中部地方下の太平洋スラブに適用し、数値モデリングを実施する準備を整えた。 中島は、遠地地震を用いたレシーバー関数解析プログラムコードを新規に作成し、日本列島下の太平洋スラブの深部構造およびマントル遷移層不連続面の空間変化を調査した。その結果、太平洋スラブの上部境界面付近に正の振幅をもつ地震波不連続面が連続的に存在すること、その不連続面は410km不連続面に繋がるようにみえること、スラブ下(裏側)の410km、660km不連続面はほぼフラットであることが明らかになった。これらの結果は、含水相の相転移の深さ、マントル対流の温度異常などの情報を含んでいると考えられる。次年度では解析のイメージング精度の向上を目指しつつ結果の解釈を行い、スラブの相転移と深発地震の関係、およびスラブ下に存在すると考えられているマントル上昇流の原因を考察する。
- Japan Society for the Promotion of Science, Grants-in-Aid for Scientific Research, Grant-in-Aid for Scientific Research (B), Kobe University, 01 Apr. 2021 - 31 Mar. 2025Why do intra-slab earthquakes occur ? Relationship with the three-dimensional distribution of olivine metastable phase, temperature, and dehydration
- Japan Society for the Promotion of Science, Grants-in-Aid for Scientific Research, Grant-in-Aid for Scientific Research on Innovative Areas (Research in a proposed research area), The University of Tokyo, 30 Jun. 2016 - 31 Mar. 2021Study on Geoscientific Modeling of Earthquake Phenomena from Low-speed Deformation to High-speed SlipThrough geoscientific modeling, we conducted research to elucidate the principles of slow earthquakes. Specifically, the research was divided into three parts: (A) clarification of the temporal and spatial relationships among slow earthquake phenomena, (B) modeling in realistic plate tectonic systems, and (C) investigation of the predictability of plate motion including giant earthquakes. In (A), we made a major contribution to understanding the relationship between apparently different phenomena, focusing on the clarification of the broadband picture of slow earthquakes. In (B), we developed various models incorporating realistic factors such as heat, water, and multiple faults. In (C), experimental studies were particularly advanced, and a major achievement was the understanding of the two-dimensional growth of the rupture nucleus. We also obtained clues for new research directions using machine learning.
- 科学研究費補助金/基盤研究(B), Apr. 2016 - Mar. 2020, Principal investigatorCompetitive research funding
- 科学研究費補助金/新学術領域研究, Apr. 2015 - Mar. 2017, Principal investigatorCompetitive research funding
- 国立研究開発法人科学技術振興機構, 国際科学技術共同研究推進事業 地球規模課題対応国際科学技術協力プログラム(SATREPS), 2017, Principal investigator【JST-SATREPS】メキシコ沿岸部の巨大地震・津波災害の軽減に向けた総合的研究Competitive research funding
- 国際科学技術共同研究推進事業 地球規模課題対応国際科学技術協力プログラム(SATREPS), 2016, Principal investigator【JST-SATREPS】メキシコ沿岸部の巨大地震・津波災害の軽減に向けた総合的研究Competitive research funding
- 科学研究費補助金/基盤研究(B), Apr. 2010 - Mar. 2015, Principal investigatorCompetitive research funding
- Japan Society for the Promotion of Science, Grants-in-Aid for Scientific Research, Grant-in-Aid for Scientific Research on Innovative Areas (Research in a proposed research area), The University of Tokyo, 23 Jul. 2009 - 31 Mar. 2014Comprehensive modeling of pre- and co-seismic processesWe studied inter-, pre-, and co-seismic process of mega-quake in the Nankai subduction zone. The main topics are (A) slow deformation on the shallow plate interface around accretionary prism, (B) dynamic rupture in geometrically complex fault system, and (C) earthquake cycle including preparation stage of a mega-quake. We successfully constructed a new comprehensive model for physical processes of mega-earthquakes in subduction zone, in a wide frequency range, through the modeling of slow earthquakes and the 2011 Tohoku-Oki earthquake.
- Japan Society for the Promotion of Science, Grants-in-Aid for Scientific Research, Grant-in-Aid for Scientific Research on Priority Areas, Japan Agency for Marine-Earth Science and Technology, 2004 - 2009Stagnant Slab : A Key Word for the Mantle Dynamics本年度は、以下の研究を実施した。 1. 平成21年7月16日に通算第9回目の総括班会議を開き、これまでの各計画研究の成果の状況を再確認し、領域全体の研究のとりまとめの方向および役割分担を議論した。 2. 欧文国際誌Physics of Earth Planetary Interiorに特集号の投稿受け付けを、前年度に開始して平成21年8月31日に締め切った。本研究領域からだけではなく、海外からも含め、37篇の論文が投稿された。平成22年度中に刊行される予定で査読が進められている。特集号「Deep Slab and Mantle Dynamics」のゲスト編集者は、末次大輔(海洋研究開発機構)、井上徹(愛媛大学地球深部ダイナミクス研究センター)、Craig Bina(米国ノースウェスタン大学)、Douglas Wiens(米国ワシントン大学)の4名であり、末次、井上が特定領域研究に参加した研究者である。論文の分野内訳は、マントル対流モデリングの論文が10篇、高温高圧鉱物学の論文が13篇、地震学・地球電磁気学など地球物理学の論文が14篇と、特定領域研究がカバーする広い学問分野にわたっている。学問分野をまたいだ論文も数多く投稿されており、特定領域研究の特徴がよく表れている。現在、15篇の採否が決定されており、残りも含めて平成22年度中には刊行される予定で編集作業が進められている。 3. 平成22年1月24日(日)に、地震研究所を会場にしてサイエンスカフェ「プレートのゆくえ」を開催した。深尾領域代表者による講演に引き続き、40数名の参加者との活発な質疑応答が行なわれ、盛況であった。 4. 平成21年3月に領域の成果をとりまとめた報告書(冊子体)を印刷した。また、領域の成果をわかりやすく書き下ろした一般向け印刷物「プレートのゆくえ」を刊行した。記事の執筆には、早稲田大学大学院政治学研究科の吉戸智明氏および大学院生の協力を得た。
- 2003 - 2009スタグナントスラブの形成及び崩落過程の解明Competitive research funding
- 2003 - 2009Numerical modeling of formation and collapse processes of stagnat slabsCompetitive research funding
- Japan Society for the Promotion of Science, Grants-in-Aid for Scientific Research, Grant-in-Aid for Scientific Research on Priority Areas, Kyushu University, 2004 - 2008Modeling of stagnation/falling processes by downward flow本研究では、主として、3次元箱型、2次元箱型の熱対流を用いたスラブ沈み込みモデルを構築した。3次元箱型モデルでは、高速な数値解法アルゴリズムの開発を進め、海洋プレートの沈み込みを実現し、410km、660kmでの相転移、660kmでの粘性ジャンプ、粘性率の温度・圧力依存性、海溝の後退などをモデルに取り入れた。数値シミュレーションの結果、滞留スラブの形成には海溝の後退が重要であることを示した。2次元モデルでは、プレートに強制的な沈み込み速度を与えるモデルと、自発的にプレートが沈み込むモデルを構築した。前者では、海溝の後退がある場合、660kmでの相転移によってスラブを浮かせるモードと下部マントルの高粘性によってスラブを支えるモードがあることを示した。また、海溝後退によって生成された滞留スラブが、下部マントルへと崩落するメカニズムとして、海溝後退の停止や前進が重要であることを示した。後者のモデルでは、マントル遷移層でのスラブの振る舞いと地表のテクトニクスの相互作用を扱うことを可能にし、様々な構造を持つ滞留スラブを再現した。また、滞留スラブの形成には、海溝の後退が非常に重要であり、滞留・崩落の過程には、2つの様式があることを明らかにした。
- Japan Society for the Promotion of Science, Grants-in-Aid for Scientific Research, Grant-in-Aid for Scientific Research (B), Hiroshima University, 2004 - 2006SEISMOLOGICAL ANALYSIS AND NUMERICAL SIMULATION FOR THE MANTLE STRUCTURE WITH SUBDUCTED SLABBroadband seismic waveform analyses were carried out to study the structure associated with stagnant slab (SS) in the northwestern Pacific. The region of model M3.11 with high velocity anomaly (HVA) and the 660 km discontinuity depth depression is bounded by zones that are delineated by model M2.0 with HVA but without depression of the discontinuity. The structural change is associated with the geochemical properties of SS where the contrast of a hydrous garnet rich layer versus a bulk of peridotite involves. Anomalously broadened P waveforms of deep focus events were observed at some regional stations after sampling the region with SS. The rays of the anomalous waves propagated in the vicinity of the regions which are modeled with M3.11 or M2.0. The broadened P waveforms may be SV converted waves within an anomalous narrow zone near the "660 km" discontinuity. The creation of such heterogeneities may involve in the complex structure of phase transformation with water and cold temperature anomaly with SS. Numerical simulations were also performed to test various physical conditions associated with SS. To explain formation process of SS, effects of adiabatic compression, viscous dissipation, and trench retreat were incorporated into a 2-D thermal-fluid flow numerical model to clarify the roles of various factors on formation process of SS. We also developed a 2-D numerical model incorporating spatio-temporal variations of physical properties such as density and viscosity, and a method of calculation for temperature distribution of subducting slab with arbitrary shape. Using the model, we tested conditions to cause depression of the 660 km discontinuity, by calculating temperature distribution along a profile across the Japanese islands where seismic velocity structure has been well determined using receiver functions and seismic tomography.
- 日本学術振興会, 科学研究費助成事業, 基盤研究(C), 九州大学, 2002 - 2004フィリピン海スラブの沈み込みに伴う温度・速度構造とスラブ内地震・地表変形の地域性1.熱と流れの3次元プログラムの開発・モデル化 熱と流れの2次元差分プログラムを基に、3次元箱型モデルの差分プログラムを完成させた。パラメターの無次元化、プログラムの並列化・高速化、計算結果の描画プログラムの開発などを行った。このプログラムを用いて、西南日本において、フィリピン海プレートの過去の沈み込み様式、斜め沈み込み、紀南海山列からの距離の違いによる異なる温度境界条件などを考慮した熱と流れの3次元沈み込みモデルの構築を行った。観測された地殻熱流量と計算値の比較を行い、プレート上面の温度分布を新たに提案し、東南海・南海地震の地震発生域を推定することに成功した。 2.四国沖〜九州東部沖のプレート間カップリングの推定と地表変形の地域性の検討 九州において、GPSデータを用いたインヴァージョン解析により、日向灘におけるプレート間カップリングと九州地域の地殻変動の関連性を明らかにした。その結果、九州の地殻変動は、プレート間カップリングに加えて、別府-島原地溝帯の深部での右横ずれと、沖縄トラフの拡大に伴うプリュームによる南部での南東方向へのマントルドラッグにより説明できることがわかった。さらに、四国のGPSデータのインヴァージョンにより、四国沖でのプレート間カップリングを推定し、1968年日向灘地震(M7.5)や1946年南海地震(M8.1)の地震時のすべり領域と今回得られたプレート間カップリングの空間分布との関連性やこれらの地震の繰り返し間隔を推定した。
- 日本学術振興会, 科学研究費助成事業, 奨励研究(A), 九州大学, 2000 - 2001背弧下のプリュームの構造・ダイナミクスと背弧拡大のメカニズムの関連性背弧下の上部マントルを数百mごとの格子点に分割し、各格子点での深さからその点での圧力と温度分布を計算した。得られた背弧下での温度分布と、パイロライト組成に対する温度-圧力ダイヤグラムを用いて、すべての格子点での温度・圧力に対する密度、ヤング率、剛性率などの各種物性パラメターの値を計算し、背弧下の上部マントル内での物性パラメターの空間分布を求めた。その結果、背弧下の上昇するプリューム内の高温領域では、これらの物性値の値が小さくなることが見出された。また、相転移境界がプリュームの流れに及ぼす影響を定量的に評価し、相転移境界がある場合には、ない場合に比べて、660km直下でプリュームが妨げられ、横長の塊状になる形状を示すこと、プリューム域を除いては温度分布のパターンは両者で似ているが、相転移がある場合は、ない場合に比べてその時間的進行が遅れること、などが明らかになった。さらに、パイロライト組成を仮定し、背弧下の上部マントル内でP波・S波速度分布、及び地震波速度異常の空間分布を求め、P波、S波とも地震波速度分布と温度分布との対応がよく、低速度域は高温域に、高速度域は低温域に対応していること、P波よりもS波速度異常の方が値が大きいこと、などがわかった。また、背弧拡大モデルに対する流線から背弧下の上部マントル内での流れの速度分布の計算を行い、浅部ではプリューム域から離れるに従って速度ベクトルの向きが鉛直方向から水平方向に向いてくること、相転移境界がない場合でも、660km付近で反流が生じる可能性があることを示した。また、背弧におけるリソスフェア下部に働く剪断応力を計算し、現在の九州の南東向きのGPSデータを説明するためには、10〜20MPa程度の断応力が必要である可能性を示した。
- 日本学術振興会, 科学研究費助成事業, 奨励研究(A), 九州大学, 1997 - 1998スラブ内部の相転移に伴う応力場と深発地震の発生の場との関連性本研究では、物性値の空間分布を考慮して、スラブ内の応力場の数値シミュレーションを行った。スラブ構成物質をオリビン(Mg0.89Fe0.11)2SiO4と仮定し、差分(FD)法を用いてスラブ内の温度分布を計算し、各タイムステップごとに物性値の空間分布を反映させた。スラブ内部の応力場に影響を与える要素として、スラブと周囲のマントルの密度差により生じる正・負の浮力、沈み込むスラブの上面・下面に働く剪断力、スラブ先端がマントルから受ける抵抗力、熱応力及び不均質体積変化により生じる応力を考慮し、有限要素法(FEM)を用いて弾性応力解析を行った。正・負の浮力に起因する応力は、α相内のdown.dip tension、γ相内に卓越するdown-dip compression、γ→Pv+Mw相転移境界付近のdown-dip tensionで特徴づけられ、最大剪断応力はオリビン準安定相境界及びγ相内深部に集中し30MPa程度を示した。剪断力に起因する応力場は、α相内では圧縮場が、γ相内ではdown-dip tensionが現れ、最大剪断応力の大きさは、数MPaと非常に小さかった。スラブ先端が受ける抵抗力に起因する応力場では、down-dip compressionが卓越し、最大剪断応力はスラブ先端部で最大値を示した。熱応力及び不均質体積変化により生じる応力は、各相転移境界で応力集中が見られた。深発地震発生数と解析結果のdown-dip方向の応力を比較した結果、正・負の浮力による応力が深発地震のピークと一致し、700km以深では最大剪断応力がゼロに近づき、サイスミシティが見られないことがうまく説明されることにより、考慮した5つの要素の中で最も良い一致を示すことがわかった。また、正・負の応力に起因する応力場とハーバード大学のグループによるCMT解との比較を試み、これについてもスラブと周囲のマントルの密度差による正・負の浮力が地震学的観測データを最もうまく説明する支配的な要素であると結論づけた。
- Japan Society for the Promotion of Science, Grants-in-Aid for Scientific Research, Grant-in-Aid for Scientific Research (B), Ehime University, 1995 - 1997High Pressure Phase Transformations Down to the Lower Mantle Conditions and Implications for the Nature of the Seismic Discontinuities.High pressure and high temperature experiments have been made on mantle materials using both quench method and in situ X-ray diffraction measurement. The phase transformations in a pyrolite composition were studied by the former method at pressures to 28 GPa and temperatures along a representative mantle geotherm.It was demonstrated that iron partitioning of the coexisting phases in this composition affects significantly the width of the pressure interval between the modified spinel and spinel phases of olivine, which is consistent with the nature of the 410km discontinuity revealed by the recent seismological studies. On the other hand, we presented some data on the partitioning of iron between perovskite and magnesiowustite under the lower mantle conditions, and showed that the iron partition coefficient approaches unity with increasing pressure as a result of the increase of aluminum in the perovskite phase. The sharpness of the 660km seismic is discussed on the basis of these experimental data. The phase boundary between spinel and the postspinel phase has been determined by in situ X-ray measurements using a combination of the synchrotron source (SPring-8) and a multianvil device. The boundary was determined at temperatures up to about 2000゚C,demonstrating that the postspinel phase boundary has a negative Clapeyron slope as estimated by quench experiments and thermodynamic analyzes. The boundary was located at about 21 GPa, at 1600゚C,which is about 2GPa lower than the earlier estimations based on other high pressure studies.
- 日本学術振興会, 科学研究費助成事業, 奨励研究(特別研究員), 京都大学, 1990 - 1990北西太平洋地域の地殼応力場の時間的変化と大地震発生の関連性
- 巨大地震に伴う粘弾性応答Competitive research funding
- Viscoelastic response associated with great earthquakesCompetitive research funding
