- earthquake and tectonic studies
- High-pressure geophysics and materials
- Seismology and Earthquake Studies
- Geological and Geochemical Analysis
- Earthquake Detection and Analysis
- Seismic Imaging and Inversion Techniques
- Seismic Waves and Analysis
- Geotechnical and Geomechanical Engineering
- Landslides and related hazards
- Dental Implant Techniques and Outcomes
- Hydraulic Fracturing and Reservoir Analysis
- Geological Modeling and Analysis
- Geophysics and Sensor Technology
- Geophysics and Gravity Measurements
- Geological Studies and Exploration
- Rock Mechanics and Modeling
- Methane Hydrates and Related Phenomena
- Geology and Paleoclimatology Research
- Geochemistry and Geologic Mapping
- Seismic Performance and Analysis
- Geomagnetism and Paleomagnetism Studies
- Oral and Maxillofacial Pathology
- Geotechnical Engineering and Underground Structures
- Shoulder Injury and Treatment
- Statistical and numerical algorithms
The University of Tokyo
2001-2024
Tokyo University of Science
2000-2024
Kurume University
2020-2024
Miyazaki Agricultural Research Institute
2024
The Institute of Statistical Mathematics
2013-2022
Research Organization of Information and Systems
2021
Showa University Hospital
2019
Showa University Dental Hospital
2019
Showa University
1994-2012
Showa University Karasuyama Hospital
2012
We developed a new inversion method to reconstruct static images of seismic sources from geodetic data, using Akaike's Bayesian Information Criterion (ABIC). Coseismic surface displacements are generally related with slip distribution on fault by linear integral equations. Parametric expansion the finite number known basis functions yields set observation equations expressed in simple vector form. Incorporating prior constraints smoothness equations, we construct model unknown...
Iterated linear inversion theory can often solve nonlinear inverse problems. Also, provides convenient error estimates and other interpretive measures. But are these measures valid for problems? We address this question in terms of the joint probability density function (pdf) estimated parameters. Briefly, will be if observations functions parameters within a reasonable (say, 95%) confidence region about optimal estimate, estimate is unique. use Bayes' rule to show how prior information...
[1] We can regard the occurrence of earthquakes as partial release tectonic stress by sudden brittle rupture. In framework linear elasticity, any indigenous source including earthquake rupture is represented a moment tensor. The tensor mathematically equivalent to volume integral over whole elastic region surrounding source, and so we quantitatively relate centroid (CMT) seismic events with an unknown field. On basis such idea Bayesian statistical inference theory, developed inversion method...
A kinematic model for the earthquake cycle at convergent plate boundaries has been constructed on basis of dislocation theory. We lithosphere-asthenosphere system by a stratified semi-infinite medium under gravity, consisting an elastic surface layer, intervening layer with Maxwell viscoelasticity, and substratum. The steady motion convergence is naturally represented uniform slip constant rate upper boundary descending oceanic plate. occurrence interplate earthquakes regarded as...
Geodetic observations at many active plate margins reveal relatively steady aseismic motion during the time between major earthquakes. We model by relative several blocks that suffer frictional resistance in upper crust. Frictional drag is represented uniform dislocation on a fault surface. assume friction occurs only “brittle zone” of each block boundary and frictionless sliding below this zone. Using Bayesian inversion procedure, we invert geodetic data to determine parameters model....
We have obtained general expressions for quasi-static internal deformation fields due to a dislocation source in multilayered elastic/viscoelastic half-space under gravity by applying the correspondence principle of linear viscoelasticity associated elastic solution (Fukahata & Matsu'ura 2005). The use upgoing propagator matrix region below and downgoing above derivation mathematical guarantees numerical stability viscoelastic over whole region. tend certain steady state with progress stress...
An inverse method to find optimum fault parameters from geodetic data with random errors is extended so as be applicable a case of the including systematic error caused by movements reference points in triangulation. Application new static displacement associated Kanto earthquake 1923 yields dislocation source model which adequately explains both seismological and data.From data, it found that motion reverse, right-lateral slip 4.8m slip-angle 140° on plane dips 25° towards N24°E, where...
The Mw 9.0 off northeast Japan (Tohoku-oki) earthquake occurred on 2011 March 11 at the interface between North American (NA) and Pacific (PA) plates. occurrence of interplate earthquakes can be regarded as sudden release tectonic stress accumulated by interseismic gradual increase slip deficit in source regions. inversion analysis GPS velocity data for a seismically calm period (1996-2000) before Tohoku-oki has revealed existence five remarkable slip-deficit zones distributed NA-PA plate...
We developed a new method of geodetic data inversion to estimate slip history at plate interface by using Akaike's Bayesian Information Criterion (ABIC). In this we considered the effects viscoelastic stress relaxation in asthenosphere, which cannot be neglected during one earthquake cycle. also introduced proper formulation incorporate two sorts partially dependent prior information into observed Bayes' rule. By applying levelling for 1893–1983 Shikoku, southwestern Japan, reconstructed...
The process of earthquake generation is governed by a coupled non-linear system consisting the equation motion in elastodynamics and fault constitutive relation. On basis results stick-slip experiments we constructed theoretical source model with slip-dependent law. Using model, simulated transition numerically from quasi-static nucleation to high-speed rupture propagation succeeded quantitatively explaining three phases observed experiments, that very slow (1 cm s−1) preceding onset dynamic...
The spatial distribution of the strength interplate coupling between subducting Philippine Sea and overlying continental plates in Tokai district, central Japan, was investigated detail through inversion analysis geodetic data using Akaike's Bayesian Information Criterion (ABIC). used for are annual rates level changes (1972–1984) horizontal length (1977–1988), which presumably represent average crustal movements during interseismic period. result shows existence a strongly coupled region...
We developed a model which can naturally explain the entire process of nucleation, dynamic propagation, and stop earthquake rupture. In this we represent frictional interaction between fault surfaces by slip‐dependent constitutive relation, use it as fundamental law governing considered broad weak zone with locally strong part (asperity) on plane examined rupture proceeding increase external shear stress through numerical simulations. The nucleation first proceeds quasi‐statically at weakest...
We have obtained general expressions for internal displacement and stress fields due to a point dislocation source in multilayered elastic half-space under gravity. Most previous the deformation were by applying one of two different types Thomson–Haskell propagator matrix, namely up-going matrix proposed Singh (1970) down-going Sato (1971). The solution derived with is stable below source, but becomes unstable above source. In contrast, succeeded unifying matrices into generalized applied it...
This paper presents a method for the computation of quasi-static surface displacements caused by sudden appearance dislocation in stratified elastic half-space with an intervenient, Maxwellian viscoelastic layer. Integral representations are derived from those associated case applying correspondence principle linear viscoelasticity. Evaluation integral is carried out using approximating part integrand analytical function. The motion consists instantaneous rebound at time faulting and...
A quasi-linear inverse problem with a priori information about model parameters is formulated in stochastic framework by using singular value decomposition technique for arbitrary rectangular matrices.In many geophysical problems, we have from which most plausible solution and the statistics of its probable error can be guessed. Starting solution, optimization made successive iteration solving set standardized linear equations corrections at each step. In under-determined cases, depends...
Surface deformations due to the fault spreading over several layers in a multi-layered medium are investigated. Machine program is developed by which contour map of vertical components and vector horizontal around focal area can be drawn, designating para-meters as inputs.Some numerical examples surface presented for purpose comparison with those obtained homogeneous semi-infinite model.For dip-slip fault, example, it known that displacement at upheaval on one side strike subsidence other...
Analytical expressions for the strain and tilt fields due to an inclined rectangular fault in a semi-infinite medium with arbitrary Poisson's ratio are derived. Since only simple algebraic functions involved, numerical computation can be quickly carried out.