Joint inversion of InSAR and broadband teleseismic waveform data with ABIC: application to the 1997 Manyi, Tibet earthquake
|
|
- Dora Little
- 5 years ago
- Views:
Transcription
1 Joint inversion of InSAR and broadband teleseismic waveform data with ABIC: application to the 1997 Manyi, Tibet earthquake Gareth Funning 1, Yukitoshi Fukahata 2, Yuji Yagi 3 & Barry Parsons 4 1 University of California, Berkeley, USA 2 University of Tokyo, Japan 3 University of Tsukuba, Japan 4 COMET, University of Oxford, UK
2 Overview Introduction Rationale for joint inversions The 1997 Manyi, Tibet earthquake Properties of the earthquake source model Inverting for source parameters with ABIC Brief description of the method InSAR-only, seismology-only and joint inversion results
3 Why attempt a joint inversion? InSAR has excellent spatial resolution (< 100 m) can map and model deformation in fine detail
4 Example: the Bam earthquake Deformation pattern consistent with slip on two parallel faults. Bam, Iran: track 385 (ascending), frame 0575 (detail) Funning et al., 2005, JGR
5 Why attempt a joint inversion? InSAR has excellent spatial resolution (< 100 m) can map deformation in fine detail but temporal resolution is poor ( 35 days) no information on timing.
6 Why attempt a joint inversion? InSAR has excellent spatial resolution (< 100 m) can map deformation in fine detail but temporal resolution is poor ( 35 days) no information on timing. Body wave seismology has good temporal resolution, and is also sensitive to location.
7 Deformation pattern consistent with slip on two parallel faults. Jackson et al., 2005, submitted to GJI
8 Deformation pattern consistent with slip on two parallel faults. Jackson et al., 2005, submitted to GJI
9 Why attempt a joint inversion? InSAR has excellent spatial resolution (< 100 m) can map deformation in fine detail but temporal resolution is poor ( 35 days) no information on timing or sequence of events. Body wave seismology has good temporal resolution, and is also sensitive to location. Aim is to generate a solution that is well constrained in space and time.
10 The 1997 Manyi, Tibet earthquake 11th November 1997, M w 7.5 largest event of that year Occurred in a remote area of N Tibet; no casualties reported, no field analysis. Ideal remote earthquake to study in this way both InSAR and seismology require no field visits.
11 The 1997 Manyi, Tibet earthquake 36 8th November 1997, M w 7.5 Aftershock pattern spans 3 satellite tracks 35 SRTM topography in 34 shaded relief
12 The 1997 Manyi, Tibet earthquake Peak LOS offset is 2.4 m ( 7 m of left-lateral offset)
13 Earthquake source model properties Define a fault geometry using the InSAR data.
14 The fault geometry Assumed fault geometry: Strike 258 Dip 90 Rake -5 Length 180 km Width 18 km 180 km 18 km (schematic)
15 The fault model Assumed fault geometry: Strike 258 Dip 90 Rake -5 Length 180 km Width 18 km 180 km Fault is divided into 6 6 km patches. 18 km (schematic)
16 Earthquake source model properties Define a fault geometry using the InSAR data. Define a source model with multiple time steps.
17 The source model The amount of slip is solved for at each of five time steps. t 1 t 2 Commencement of slip on each patch occurs after a rupture front has reached the patch. t 3 t 4 t 5 Each time step has a duration of 2 seconds. increasing time
18 Earthquake source model properties Define a fault geometry using the InSAR data. Define a source model with multiple time steps. Model must be smooth in space oscillations imply unphysical strains on the fault.
19 Spatial smoothing The total slip on a given patch a 1 a 4 a 7 a M a 2 a 5 a 8 a 3 a 6 a 9 (schematic)
20 Spatial smoothing The total slip on a given patch is related to that of patches adjacent to it, by a finitedifference Laplacian approximation. (a 2 a 5 ) (a 5 a 8 )+(a 4 a 5 ) (a 5 a 6 ) = 0 a 2 + a 4 4a 5 + a 6 + a 8 = 0 a 1 a 4 a 7 a 2 a 5 a 8 a M a 3 a 6 a 9 (schematic)
21 Earthquake source model properties Define a fault geometry using the InSAR data. Define a source model with multiple time steps. Model must be smooth in space oscillations imply unphysical strains on the fault Model must be smooth in time oscillations imply changes in acceleration.
22 Temporal smoothing The slip on a given patch, at a given time step t 2 a 95 (schematic)
23 Temporal smoothing The slip on a given patch, at a given time step t 1 t 2 a 5 is related to that for the same patch at the previous and subsequent time steps. a 95 a 185 t 3 increasing time (schematic)
24 Temporal smoothing The slip on a given patch, at a given time step is related to that for the same patch at the previous and subsequent time steps. a 5 a 95 a 185 t 1 t 2 t 3 Laplacian smoothing in 1D is used here (a 5 a 95 ) (a 95 a 185 ) = 0 a 5 2a 95 + a 185 = 0
25 Earthquake source model properties Define a fault geometry using the InSAR data. Define a source model with multiple time steps. Model must be smooth in space oscillations imply unphysical strains on the fault Model must be smooth in time oscillations imply changes in acceleration Model must also be consistent with the available data, of course
26 Inverting data using ABIC Akaike s Bayesian Information Criterion (Akaike, 1980) What level of resolution is appropriate, given the input data and fault geometry? Calculation of ABIC allows the optimal estimation of the relative weighting of different constraints on the inversion.
27 Inversions using ABIC Model, a, has two types of constraints: Observations: d = Ha + e o d = observed data H = data kernels a = fault slip e o = errors in observations Smoothing: Sa + e s = 0 Ta + e t = 0 S = 2D Laplacian operator e s = errors in smoothing T = 1D Laplacian operator e t = errors in smoothing
28 Inversions using ABIC Observation equation: Smoothing equations: d = Ha + e o Sa + e s = 0 Ta + e t = 0
29 Inversions using ABIC Data distribution: Prior distribution: p(d a;σ 2,γ 2 ) = (2πσ 2 ) -N/2 E -1/2 exp[-(1/2σ 2 )(d-ha) T E -1 (d-ha)] p(a;ρ s2, ρ t2 ) = (2π) -M/2 (1/ρ s2 )G s + (1/ρ t2 )G t 1/2 exp{-a T [(1/2ρ s2 )G s + (1/2ρ t2 )G t ]a}
30 Inversions using ABIC Data distribution: Prior distribution: p(d a;σ 2,γ 2 ) = (2πσ 2 ) -N/2 E -1/2 exp[-(1/2σ 2 )(d-ha) T E -1 (d-ha)] p(a;ρ s2,ρ t2 ) = (2π) -M/2 (1/ρ s2 )G s + (1/ρ t2 )G t 1/2 exp{-a T [(1/2ρ s2 )G s + (1/2ρ t2 )G t ]a} Combined (posterior) distribution: p(a;σ 2,α 2,β 2,γ 2 d) = cp(d a;σ 2,γ 2 ) p(a;ρ s2,ρ t2 )
31 Inversions using ABIC Data distribution: Prior distribution: p(d a;σ 2,γ 2 ) = (2πσ 2 ) -N/2 E -1/2 exp[-(1/2σ 2 )(d-ha) T E -1 (d-ha)] p(a;ρ s2,ρ t2 ) = (2π) -M/2 (1/ρ s2 )G s + (1/ρ t2 )G t 1/2 exp{-a T [(1/2ρ s2 )G s + (1/2ρ t2 )G t ]a} Combined (posterior) distribution: p(a;σ 2,α 2,β 2,γ 2 d) = cp(d a;σ 2,γ 2 ) p(a;ρ s2,ρ t2 ) ABIC: ABIC(α 2,β 2,γ 2 )= -2 log p(a;σ 2,α 2,β 2,γ 2 d) da
32 Inverting data using ABIC ABIC is a function of (in this case) up to three hyperparameters, α 2, β 2 and γ 2. Evaluate ABIC numerically for a range of values. The hyperparameter values that give a minimum are the optimum values. Recover best-fitting set of model parameters, a*, for each model, as part of this process.
33 Modelling the Manyi earthquake I calculate three models for the earthquake, using the ABIC technique: 1. InSAR data only 2. Seismic data only 3. Both datasets jointly In all cases I use the same simplified fault geometry
34 InSAR data inversion Simplest case 1 time step. No temporal smoothing. Only 1 hyperparameter α 2, the relative weighting of data and smoothing.
35 ABIC for the InSAR model Evaluate ABIC for a range of values of α 2 until a minimum is obtained
36 InSAR model results Majority of significant slip occurs in the upper 12 km of the fault
37 InSAR model fit to data Fit to data is reasonable for this very crude geometry.
38 Seismic data inversion 5 time steps, temporal smoothing is required. 2 hyperparameters α 2 and β 2, controlling spatial and temporal smoothing, respectively. Seismic kernels calculated by Yuji Yagi using the method of Kikuchi and Kanamori (1991).
39 Broadband teleseismic data Broadband vertical component data used from 10 GSN stations
40 Broadband teleseismic data 80 s time window, Band-pass filtered between 0.01 and 0.8 Hz Sample interval 0.25 s.
41 ABIC for the seismic model Minimum at α 2 = 0.008, β 2 =0.8, again obtained numerically
42 Seismic model results Greater slip in the lower 6 km of the fault, and less at the upper 6 km than for the InSAR model.
43 Fit to data of the seismic model The fit to the waveforms is generally good.
44 Fit to data of the seismic model Total slip gives a poor fit to the InSAR data
45 Joint inversion 5 time steps; temporal smoothing required. 3 hyperparameters α 2, β 2 and γ 2 controlling spatial smoothing, temporal smoothing and relative weighting of the datasets.
46 ABIC for the joint model Minimum at α 2 = , β 2 =0.0006, γ 2 =117
47 Joint model results Less slip in the lower 6 km than for the seismic-only inversion it is relocated to the surface.
48 Fit to data of the joint model Total slip gives a fit comparable to that for the InSAR-only inversion
49 Fit to data of the joint model The fit to the waveforms is very similar to that of the seismic-only inversion. (degradation of 2% in misfit)
50 Slip history of the Manyi earthquake
51 Comparison of model results InSAR Seis. Joint
52 Conclusions It is possible to find a solution with ABIC which satisfies both InSAR and seismic datasets. InSAR is the main control on the spatial location of slip, which occurred mainly at shallow depths. This is a viable technique to apply to other large continental earthquakes where near-source information is not available.
53
Buried Strike Slip Faults: The 1994 and 2004 Al Hoceima, Morocco Earthquakes.
Buried Strike Slip Faults: The 1994 and 2004 Al Hoceima, Morocco Earthquakes. Juliet Biggs 1, Eric Bergman 2, Brian Emmerson 3, Gareth Funning 4, James Jackson 3, Barry Parsons 1,Tim Wright 1. 1 COMET,
More informationHybrid Back-Projection: Theory and Future Direction of Seismic Array Analysis
Hybrid Back-Proection: Theory and Future Direction of Seismic Array Analysis YAGI Yui (Univ. of Tsukuba) Collaborators: FUKAHATA Yukitoshi (DPRI, Kyoto University) KASAHARA Amato, NAKAO Atsushi, EZAKI
More informationSOURCE PROCESS OF THE 2003 PUERTO PLATA EARTHQUAKE USING TELESEISMIC DATA AND STRONG GROUND MOTION SIMULATION
Synopses of Master Papers Bulletin of IISEE, 47, 19-24, 2013 SOURCE PROCESS OF THE 2003 PUERTO PLATA EARTHQUAKE USING TELESEISMIC DATA AND STRONG GROUND MOTION SIMULATION Fabricio Moquete Everth* Supervisor:
More informationTwo Contrasting InSAR Studies of Recent Earthquakes in Tibet
Two Contrasting InSAR Studies of Recent Earthquakes in Tibet Barry Parsons Department of Earth Sciences University of Oxford John Elliott, Wanpeng Feng,, James Jackson, Zhenhong Li, Xinjian Shan, Alastair
More informationJournal of Geophysical Research Letters Supporting Information for
Journal of Geophysical Research Letters Supporting Information for InSAR observations of strain accumulation and fault creep along the Chaman Fault system, Pakistan and Afghanistan H. Fattahi 1, F. Amelung
More information27th Seismic Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
GROUND TRUTH LOCATIONS USING SYNERGY BETWEEN REMOTE SENSING AND SEISMIC METHODS-APPLICATION TO CHINESE AND NORTH AFRICAN EARTHQUAKES C. K. Saikia 1, H. K. Thio 2, D. V. Helmberger 2, G. Ichinose 1, and
More informationInversion of Earthquake Rupture Process:Theory and Applications
Inversion of Earthquake Rupture Process:Theory and Applications Yun-tai CHEN 12 * Yong ZHANG 12 Li-sheng XU 2 1School of the Earth and Space Sciences, Peking University, Beijing 100871 2Institute of Geophysics,
More informationPost-seismic motion following the 1997 Manyi (Tibet) earthquake: InSAR observations and modelling
Geophys. J. Int. (7) 69, 9 7 doi:./j.365-46x.6.33.x Post-seismic motion following the 997 Manyi (Tibet) earthquake: InSAR observations and modelling Isabelle Ryder, Barry Parsons, Tim J. Wright and Gareth
More informationRELOCATION OF THE MACHAZE AND LACERDA EARTHQUAKES IN MOZAMBIQUE AND THE RUPTURE PROCESS OF THE 2006 Mw7.0 MACHAZE EARTHQUAKE
RELOCATION OF THE MACHAZE AND LACERDA EARTHQUAKES IN MOZAMBIQUE AND THE RUPTURE PROCESS OF THE 2006 Mw7.0 MACHAZE EARTHQUAKE Paulino C. FEITIO* Supervisors: Nobuo HURUKAWA** MEE07165 Toshiaki YOKOI** ABSTRACT
More informationMoment tensor inversion of near source seismograms
Moment tensor inversion of near source seismograms Yuji Yagi and Naoki Nishimura ABSTRACT We construct a program set for estimating moment tensor solution using near source seismograms. We take the effect
More informationOn the Limitation of Receiver Functions Method: Beyond Conventional Assumptions & Advanced Inversion Techniques
On the Limitation of Receiver Functions Method: Beyond Conventional Assumptions & Advanced Inversion Techniques Hrvoje Tkalčić RSES, ANU Acknowledgment: RSES Seismology & Mathematical Geophysics Group
More information27th Seismic Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
GROUND TRUTH OF AFRICAN AND EASTERN MEDITERRANEAN SHALLOW SEISMICITY USING SAR INTERFEROMETRY AND GIBBS SAMPLING INVERSION Benjamin A. Brooks 1, Francisco Gomez 2, Eric A. Sandvol 2, and Neil Frazer 1
More informationWidespread Ground Motion Distribution Caused by Rupture Directivity during the 2015 Gorkha, Nepal Earthquake
Widespread Ground Motion Distribution Caused by Rupture Directivity during the 2015 Gorkha, Nepal Earthquake Kazuki Koketsu 1, Hiroe Miyake 2, Srinagesh Davuluri 3 and Soma Nath Sapkota 4 1. Corresponding
More informationEmpirical Green s Function Analysis of the Wells, Nevada, Earthquake Source
Nevada Bureau of Mines and Geology Special Publication 36 Empirical Green s Function Analysis of the Wells, Nevada, Earthquake Source by Mendoza, C. 1 and Hartzell S. 2 1 Centro de Geociencias, Universidad
More informationSource rupture process of the 2003 Tokachi-oki earthquake determined by joint inversion of teleseismic body wave and strong ground motion data
LETTER Earth Planets Space, 56, 311 316, 2004 Source rupture process of the 2003 Tokachi-oki earthquake determined by joint inversion of teleseismic body wave and strong ground motion data Yuji Yagi International
More informationSUPPLEMENTARY INFORMATION
doi:10.1038/nature11492 Figure S1 Short-period Seismic Energy Release Pattern Imaged by F-net. (a) Locations of broadband seismograph stations in Japanese F-net used for the 0.5-2.0 Hz P wave back-projection
More informationOpen Access FULL PAPER. Yukitoshi Fukahata * and Manabu Hashimoto
DOI 1.116/s63--5- FULL PAPER Open Access Simultaneous estimation of the dip angles and slip distribution on the faults of the Kumamoto earthquake through a weak nonlinear inversion of InSAR data Yukitoshi
More informationDisplacement field and slip distribution of the 2005 Kashmir earthquake from SAR imagery
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 33,, doi:10.1029/2006gl027193, 2006 Displacement field and slip distribution of the 2005 Kashmir earthquake from SAR imagery E. Pathier, 1
More informationSeismic Source Mechanism
Seismic Source Mechanism Yuji Yagi (University of Tsukuba) Earthquake Earthquake is a term used to describe both failure process along a fault zone, and the resulting ground shaking and radiated seismic
More informationComplicated repeating earthquakes on the convergent plate boundary: Rupture processes of the 1978 and 2005 Miyagi-ken Oki earthquakes
Complicated repeating earthquakes on the convergent plate boundary: Rupture processes of the 1978 and 2005 Miyagi-ken Oki earthquakes Changjiang Wu 1 and Kazuki Koketsu Earthquake Research Institute, University
More informationDear editors and reviewer(s), thank for your comments and suggestions. Replies as follows:
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 Dear editors and reviewer(s), thank for your comments and suggestions. Replies as follows:
More informationGeodesy (InSAR, GPS, Gravity) and Big Earthquakes
Geodesy (InSAR, GPS, Gravity) and Big Earthquakes Mathew Pritchard Teh-Ru A. Song Yuri Fialko Luis Rivera Mark Simons UJNR Earthquake Research Panel, Morioka, Japan - Nov 6, 2002 Goals Accurate and high
More informationReceived 17 June 2011; revised 20 March 2013; accepted 7 June 2013; published 27 August 2013.
JOURNAL OF GEOPHYSICAL RESEARCH: SOLID EARTH, VOL. 118, 4599 4627, doi:1.12/jgrb.5244, 213 A three-step maximum a posteriori probability method for InSAR data inversion of coseismic rupture with application
More informationData Repository: Seismic and Geodetic Evidence For Extensive, Long-Lived Fault Damage Zones
DR2009082 Data Repository: Seismic and Geodetic Evidence For Extensive, Long-Lived Fault Damage Zones Fault Zone Trapped Wave Data and Methods Fault zone trapped waves observed for 2 shots and 5 local
More informationStrong-Motion and Teleseismic Waveform Inversions for the Source Process of the 2003 Bam, Iran, Earthquake
Bulletin of the Seismological Society of America, Vol. 12, No. 4, pp. 1477 1496, August 212, doi: 1.178/1211198 Strong-Motion and Teleseismic Waveform Inversions for the Source Process of the 23 Bam, Iran,
More informationSOURCE MODELING OF RECENT LARGE INLAND CRUSTAL EARTHQUAKES IN JAPAN AND SOURCE CHARACTERIZATION FOR STRONG MOTION PREDICTION
SOURCE MODELING OF RECENT LARGE INLAND CRUSTAL EARTHQUAKES IN JAPAN AND SOURCE CHARACTERIZATION FOR STRONG MOTION PREDICTION Kimiyuki Asano 1 and Tomotaka Iwata 2 1 Assistant Professor, Disaster Prevention
More informationBROADBAND STRONG MOTION SIMULATION OF THE 2004 NIIGATA- KEN CHUETSU EARTHQUAKE: SOURCE AND SITE EFFECTS
Third International Symposium on the Effects of Surface Geology on Seismic Motion Grenoble, France, 30 August - 1 September 2006 Paper Number: 105 BROADBAND STRONG MOTION SIMULATION OF THE 2004 NIIGATA-
More information14 S. 11/12/96 Mw S. 6/23/01 Mw S 20 S 22 S. Peru. 7/30/95 Mw S. Chile. Argentina. 26 S 10 cm 76 W 74 W 72 W 70 W 68 W
175 Chapter 5 Comparision of co-seismic and post-seismic slip from the November 12, 1996, M w 7.7 and the June 23, 2001, M w 8.4 southern Peru subduction zone earthquakes 176 Abstract We use InSAR and
More informationGeodetic data inversion using ABIC to estimate slip history during one earthquake cycle with viscoelastic slip-response functions
Geophys. J. Int. () 15, 1 153 doi: 1.1111/j.135-X..1.x Geodetic data inversion using ABIC to estimate slip history during one earthquake cycle with viscoelastic slip-response functions Yukitoshi Fukahata,
More informationJoint inversion of strong motion, teleseismic, geodetic, and tsunami datasets for the rupture process of the 2011 Tohoku earthquake
GEOPHYSICAL RESEARCH LETTERS, VOL. 38,, doi:10.1029/2011gl050098, 2011 Joint inversion of strong motion, teleseismic, geodetic, and tsunami datasets for the rupture process of the 2011 Tohoku earthquake
More informationDynamic source inversion for physical parameters controlling the 2017 Lesvos earthquake
Dynamic source inversion for physical parameters controlling the 2017 Lesvos earthquake F. Kostka 1, E. Sokos 2, J. Zahradník 1, F. Gallovič 1 1 Charles University in Prague, Faculty of Mathematics and
More informationJournal of Geophysical Research - Solid Earth
Journal of Geophysical Research - Solid Earth Supporting Information for Transpressional Rupture Cascade of the 2016 M w 7.8 Kaikoura Earthquake, New Zealand Wenbin Xu 1*, Guangcai Feng 2*, Lingsen Meng
More informationSlip distributions of the 1944 Tonankai and 1946 Nankai earthquakes including the horizontal movement effect on tsunami generation
Slip distributions of the 1944 Tonankai and 1946 Nankai earthquakes including the horizontal movement effect on tsunami generation Toshitaka Baba Research Program for Plate Dynamics, Institute for Frontier
More informationBasics of the modelling of the ground deformations produced by an earthquake. EO Summer School 2014 Frascati August 13 Pierre Briole
Basics of the modelling of the ground deformations produced by an earthquake EO Summer School 2014 Frascati August 13 Pierre Briole Content Earthquakes and faults Examples of SAR interferograms of earthquakes
More informationCoseismic slip distribution of the 2005 off Miyagi earthquake (M7.2) estimated by inversion of teleseismic and regional seismograms
Coseismic slip distribution of the 2005 off Miyagi earthquake (M7.2) estimated by inversion of teleseismic and regional seismograms Tadashi Yaginuma 1, Tomomi Okada 1, Yuji Yagi 2, Toru Matsuzawa 1, Norihito
More informationJCR (2 ), JGR- (1 ) (4 ) 11, EPSL GRL BSSA
Dun Wang ( ) In collaboration with: Hitoshi Kawakatsu, Jim Mori, Kazuki Koketsu, Takuto Maeda, Hiroshi Tsuroka, Jiancang Zhunag, Lihua Fang, and Qiang Yao School of Geosciences, China University of Geosciences
More informationEARTHQUAKE SOURCE PARAMETERS OF MODERATELY EARTHQUAKE IN THE SOUTH EASTERN IRAN BASED ON TELESEISMIC AND REGIONAL DISTANCES
EARTHQUAKE SOURCE PARAMETERS OF MODERATELY EARTHQUAKE IN THE SOUTH EASTERN IRAN BASED ON TELESEISMIC AND REGIONAL DISTANCES M. Mostafazadeh 1 and S.Ashkpour Motlagh 2 1 Assistant Professor,Seismology Research
More informationDETERMINATION OF SLIP DISTRIBUTION OF THE 28 MARCH 2005 NIAS EARTHQUAKE USING JOINT INVERSION OF TSUNAMI WAVEFORM AND GPS DATA
Synopses of Master Papers Bulletin of IISEE, 47, 115-10, 013 DETERMINATION OF SLIP DISTRIBUTION OF THE 8 MARCH 005 NIAS EARTHQUAKE USING JOINT INVERSION OF TSUNAMI WAVEFORM AND GPS DATA Tatok Yatimantoro
More informationFault identification for buried strike-slip earthquakes using InSAR: The 1994 and 2004 Al Hoceima, Morocco earthquakes
Geophys. J. Int. (2006) 166, 1347 1362 doi: 10.1111/j.1365-246X.2006.03071.x Fault identification for buried strike-slip earthquakes using InSAR: The 1994 and 2004 Al Hoceima, Morocco earthquakes Juliet
More informationSource process of the 2011 off the Pacific coast of Tohoku Earthquake with the combination of teleseismic and strong motion data
LETTER Earth Planets Space, 63, 565 569, 2011 Source process of the 2011 off the Pacific coast of Tohoku Earthquake with the combination of teleseismic and strong motion data Yasuhiro Yoshida 1, Hiroshi
More informationThe April 6 th 2009, L Aquila (Italy) earthquake: DInSAR analysis and seismic source model inversion
ESA ESRIN 30th November - 4th December 2009 Frascati, Italy The April 6 th 2009, L Aquila (Italy) earthquake: DInSAR analysis and seismic source model inversion Simone Atzori, Christian Bignami, Marco
More informationMeasurement of differential rupture durations as constraints on the source finiteness of deep-focus earthquakes
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 111,, doi:10.1029/2005jb004001, 2006 Measurement of differential rupture durations as constraints on the source finiteness of deep-focus earthquakes Linda M. Warren
More informationSurface displacements and source parameters of the 2003 Bam (Iran) earthquake from Envisat advanced synthetic aperture radar imagery
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 110,, doi:10.1029/2004jb003338, 2005 Surface displacements and source parameters of the 2003 Bam (Iran) earthquake from Envisat advanced synthetic aperture radar imagery
More information7 Ground Motion Models
7 Ground Motion Models 7.1 Introduction Ground motion equations are often called attenution relations but they describe much more than just the attenutation of the ground motion; they describe the probability
More informationMethod to Determine Appropriate Source Models of Large Earthquakes Including Tsunami Earthquakes for Tsunami Early Warning in Central America
Pure Appl. Geophys. 174 (2017), 3237 3248 Ó 2017 The Author(s) This article is published with open access at Springerlink.com DOI 10.1007/s00024-017-1630-y Pure and Applied Geophysics Method to Determine
More information3.3. Waveform Cross-Correlation, Earthquake Locations and HYPODD
3.3. Waveform Cross-Correlation, Earthquake Locations and HYPODD 3.3.1 Method More accurate relative earthquake locations depend on more precise relative phase arrival observations so I exploit the similarity
More informationLocation and mechanism of the Little Skull Mountain earthquake as constrained by satellite radar interferometry and seismic waveform modeling
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 107, NO. B6, 2118, 10.1029/2001JB000627, 2002 Location and mechanism of the Little Skull Mountain earthquake as constrained by satellite radar interferometry and seismic
More informationSource rupture process inversion of the 2013 Lushan earthquake, China
Geodesy and Geodynamics 2013,4(2) :16-21 http://www. jgg09. com Doi:10.3724/SP.J.1246.2013.02016 Source rupture process inversion of the 2013 Lushan earthquake, China Zhang Lifen 1 ' 2, lman Fatchurochman
More informationDepth (Km) + u ( ξ,t) u = v pl. η= Pa s. Distance from Nankai Trough (Km) u(ξ,τ) dξdτ. w(x,t) = G L (x,t τ;ξ,0) t + u(ξ,t) u(ξ,t) = v pl
Slip history during one earthquake cycle at the Nankai subduction zone, inferred from the inversion analysis of levelling data with a viscoelastic slip response function Mitsuhiro Matsu'ura, Akira Nishitani
More informationJOURNAL OF GEOPHYSICAL RESEARCH, VOL. 116, B01406, doi: /2010jb007849, 2011
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 116,, doi:10.1029/2010jb007849, 2011 Mechanical constraints on inversion of coseismic geodetic data for fault slip and geometry: Example from InSAR observation of
More informationDistribution of slip from 11 M w > 6 earthquakes in the northern Chile subduction zone
Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 111,, doi:10.1029/2005jb004013, 2006 Distribution of slip from 11 M w > 6 earthquakes in the northern Chile subduction zone M. E. Pritchard,
More informationIntroduction to Displacement Modeling
Introduction to Displacement Modeling Introduction Deformation on the Earth surface informs us about processes and material properties below surface Observation tools: GPS (static-dynamic) InSAR (static)
More informationFULL MOMENT TENSOR ANALYSIS USING FIRST MOTION DATA AT THE GEYSERS GEOTHERMAL FIELD
PROCEEDINGS, Thirty-Eighth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 11-13, 2013 SGP-TR-198 FULL MOMENT TENSOR ANALYSIS USING FIRST MOTION DATA AT
More information28th Seismic Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
GROUND TRUTH LOCATIONS USING SYNERGY BETWEEN REMOTE SENSING AND SEISMIC METHODS: SSSC AT IMS STATIONS FOR TIBETAN PLATEAU EARTHQUAKES Gene A. Ichinose 1, Chandan K. Saikia 2*, Donald V. Helmberger 3, and
More informationThe IISEE earthquake catalog, Catalog of Damaging Earthquakes in the World, IISEE-NET,, and BRI strong motion observation
The IISEE earthquake catalog, Catalog of Damaging Earthquakes in the World, IISEE-NET,, and BRI strong motion observation Tatsuhiko Hara International Institute of Seismology and Earthquake Engineering,
More informationVertical to Horizontal (V/H) Ratios for Large Megathrust Subduction Zone Earthquakes
Vertical to Horizontal (V/H) Ratios for Large Megathrust Subduction Zone Earthquakes N.J. Gregor Consultant, Oakland, California, USA N.A. Abrahamson University of California, Berkeley, USA K.O. Addo BC
More informationCHARACTERIZING EARTHQUAKE SLIP MODELS FOR THE PREDICTION OF STRONG GROUND MOTION
CHARACTERIZING EARTHQUAKE SLIP MODELS FOR THE PREDICTION OF STRONG GROUND MOTION P G SOMERVILLE 1, K IRIKURA, N ABRAHAMSON 3, S SAWADA 4, T KAGAWA 5 And Y TATSUMI 6 SUMMARY Over the past fifteen years,
More informationSpatio-temporal variation in slip rate on the plate boundary off Sanriku, northeastern Japan, estimated from small repeating earthquakes
Spatio-temporal variation in slip rate on the plate boundary off Sanriku, northeastern Japan, estimated from small repeating earthquakes T. Matsuzawa, N. Uchida, T. Igarashi *, N. Umino, and A. Hasegawa
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2012) doi: 10.1111/j.1365-246X.2012.05602.x Low-frequency and broad-band source models for the 2009 L Aquila, Italy, earthquake Natalia Poiata, 1 Kazuki
More informationJournal of Geophysical Research (Solid Earth) Supporting Information for
Journal of Geophysical Research (Solid Earth) Supporting Information for Postseismic Relocking of the Subduction Megathrust Following the 2007 Pisco, Peru earthquake D.Remy (a), H.Perfettini (b), N.Cotte
More informationShort Note Source Mechanism and Rupture Directivity of the 18 May 2009 M W 4.6 Inglewood, California, Earthquake
Bulletin of the Seismological Society of America, Vol. 100, No. 6, pp. 3269 3277, December 2010, doi: 10.1785/0120100087 Short Note Source Mechanism and Rupture Directivity of the 18 May 2009 M W 4.6 Inglewood,
More informationCase Study 1: 2014 Chiang Rai Sequence
Case Study 1: 2014 Chiang Rai Sequence Overview Mw 6.1 earthquake on 5 May 2014 at 11:08:43 UTC Largest recorded earthquake in Thailand Fault Orientation How does the orientation of the fault affect the
More informationSource of the July 2006 West Java tsunami estimated from tide gauge records
GEOPHYSICAL RESEARCH LETTERS, VOL. 33, L24317, doi:10.1029/2006gl028049, 2006 Source of the July 2006 West Java tsunami estimated from tide gauge records Yushiro Fujii 1 and Kenji Satake 2 Received 13
More informationThe combined inversion of seismic and geodetic data for the source process of the 16 October,
1 The combined inversion of seismic and geodetic data for the source process of the October, 1999, Mw7.1 Hector Mine, California, earthquake Asya Kaverina, Douglas Dreger and Evelyn Price Berkeley Seismological
More informationConstraining earthquake source inversions with GPS data: 2. A two-step approach to combine seismic and geodetic data sets
Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 114,, doi:10.1029/2008jb005746, 2009 Constraining earthquake source inversions with GPS data: 2. A two-step approach to combine seismic
More informationAutomatic Moment Tensor Analyses, In-Situ Stress Estimation and Temporal Stress Changes at The Geysers EGS Demonstration Project
PROCEEDINGS, 42nd Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 13-15, 2017 SGP-TR-212 Automatic Moment Tensor Analyses, In-Situ Stress Estimation and
More informationSonglin Li 1, Xiaoling Lai 1 Zhixiang Yao 2 and Qing Yang 1. 1 Introduction
Earthq Sci (2009)22: 417 424 417 Doi: 10.1007/s11589-009-0417-3 Fault zone structures of northern and southern portions of the main central fault generated by the 2008 Wenchuan earthquake using fault zone
More informationTsunami Simulation of 2009 Dusky Sound Earthquake in New Zealand
Tsunami Simulation of 2009 Dusky Sound Earthquake in New Zealand Polina Berezina 1 Institute of Geology, Taras Shevchenko National University of Kyiv, Kyiv, Ukraine Supervisor: Prof. Kenji Satake Earthquake
More informationSpatial and Temporal Distribution of Slip for the 1999 Chi-Chi, Taiwan, Earthquake
Bulletin of the Seismological Society of America, 91, 5, pp. 1069 1087, October 2001 Spatial and Temporal Distribution of Slip for the 1999 Chi-Chi, Taiwan, Earthquake by Kuo-Fong Ma, Jim Mori, Shiann-Jong
More informationAn intermediate deep earthquake rupturing on a dip-bending fault: Waveform analysis of the 2003 Miyagi-ken Oki earthquake
GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L24619, doi:10.1029/2004gl021228, 2004 An intermediate deep earthquake rupturing on a dip-bending fault: Waveform analysis of the 2003 Miyagi-ken Oki earthquake Changjiang
More informationThe March 11, 2011, Tohoku-oki earthquake (Japan): surface displacement and source modelling
The March 11, 2011, Tohoku-oki earthquake (Japan): surface displacement and source modelling Salvatore Stramondo Bignami C., Borgstrom S., Chini M., Guglielmino F., Melini D., Puglisi G., Siniscalchi V.,
More informationRupture Process of the Great 2004 Sumatra-Andaman Earthquake
Rupture Process of the Great 2004 Sumatra-Andaman Earthquake Supporting Online Materials Submitted to Science, March 12, 2005 Charles J. Ammon 1, Ji Chen 2, Hong-Kie Thio 3, David Robinson 5, Sidao Ni
More informationEffect of the Emperor seamounts on trans-oceanic propagation of the 2006 Kuril Island earthquake tsunami
GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L02611, doi:10.1029/2007gl032129, 2008 Effect of the Emperor seamounts on trans-oceanic propagation of the 2006 Kuril Island earthquake tsunami S. Koshimura, 1 Y.
More informationDynamic Triggering Semi-Volcanic Tremor in Japanese Volcanic Region by The 2016 Mw 7.0 Kumamoto Earthquake
Dynamic Triggering Semi-Volcanic Tremor in Japanese Volcanic Region by The 016 Mw 7.0 Kumamoto Earthquake Heng-Yi Su 1 *, Aitaro Kato 1 Department of Earth Sciences, National Central University, Taoyuan
More informationEffects of Surface Geology on Seismic Motion
4 th IASPEI / IAEE International Symposium: Effects of Surface Geology on Seismic Motion August 23 26, 2011 University of California Santa Barbara LONG-PERIOD (3 TO 10 S) GROUND MOTIONS IN AND AROUND THE
More informationBuilding up Seismsic Models for Ground Motion Prediction of Taiwan: Problems and Challenges
Building up Seismsic Models for Ground Motion Prediction of Taiwan: Problems and Challenges Kuo-Fong Ma 馬國鳳 Institute of Geophysics National Central University What approaches we can make toward a reliable
More informationRupture Characteristics of Major and Great (M w 7.0) Megathrust Earthquakes from : 1. Source Parameter Scaling Relationships
Journal of Geophysical Research Solid Earth Supporting Information for Rupture Characteristics of Major and Great (M w 7.0) Megathrust Earthquakes from 1990-2015: 1. Source Parameter Scaling Relationships
More informationSource studies of the ongoing ( ) sequence of recent large earthquakes in Canterbury
Source studies of the ongoing (00-0) sequence of recent large earthquakes in Canterbury C. Holden & J. Beavan GNS Science New Zealand) SUMMARY: On September, 00, a surface rupturing crustal earthquake
More informationSimulation of Strong Ground Motions for a Shallow Crustal Earthquake in Japan Based on the Pseudo Point-Source Model
6 th International Conference on Earthquake Geotechnical Engineering -4 November 25 Christchurch, New Zealand Simulation of Strong Ground Motions for a Shallow Crustal Earthquake in Japan Based on the
More informationEXAMINATION ON CONSECUTIVE RUPTURING OF TWO CLOSE FAULTS BY DYNAMIC SIMULATION
EXAMINATION ON CONSECUTIVE RUPTURING OF TWO CLOSE FAULTS BY DYNAMIC SIMULATION M. Muto 1, K. Dan 1, H. Torita 1, Y. Ohashi 1, and Y. Kase 2 1 Ohsaki Research Institute, Inc., Tokyo, Japan 2 National Institute
More informationto: Interseismic strain accumulation and the earthquake potential on the southern San
Supplementary material to: Interseismic strain accumulation and the earthquake potential on the southern San Andreas fault system by Yuri Fialko Methods The San Bernardino-Coachella Valley segment of the
More informationMeasuring seismicity in the Groningen Field. Bernard Dost, Elmer Ruigrok, Jesper Spetzler, Gert-Jan van den Hazel, Jordi Domingo
Measuring seismicity in the Groningen Field Bernard Dost, Elmer Ruigrok, Jesper Spetzler, Gert-Jan van den Hazel, Jordi Domingo Monitoring induced seismicity in the Netherlands. Instrumentation and network
More informationSimulation of earthquake rupture process and strong ground motion
Simulation of earthquake rupture process and strong ground motion Takashi Miyatake (1) and Tomohiro Inoue (2) (1) Earthquake Research Institute, University of Tokyo, Yayoi, Bunkyo, Tokyo, 113-0032, Japan
More informationby A. Ozgun Konca, Sebastien Leprince, Jean-Philippe Avouac, and Don V. Helmberger
Bulletin of the Seismological Society of America, Vol. 100, No. 1, pp. 267 288, February 2010, doi: 10.1785/0120090072 E Rupture Process of the 1999 M w 7.1 Duzce Earthquake from Joint Analysis of SPOT,
More informationEarthquake patterns in the Flinders Ranges - Temporary network , preliminary results
Earthquake patterns in the Flinders Ranges - Temporary network 2003-2006, preliminary results Objectives David Love 1, Phil Cummins 2, Natalie Balfour 3 1 Primary Industries and Resources South Australia
More informationThe source process of the 2001 July 26 Skyros Island (Greece) earthquake
Geophys. J. Int. (2004) 156, 541 548 doi: 10.1111/j.1365-246X.2004.02124.x The source process of the 2001 July 26 Skyros Island (Greece) earthquake Zafeiria Roumelioti, 1 Anastasia Kiratzi 1 and Douglas
More informationRAPID SOURCE PARAMETER DETERMINATION AND EARTHQUAKE SOURCE PROCESS IN INDONESIA REGION
RAPI SOURCE PARAETER ETERINATION AN EARTHQUAKE SOURCE PROCESS IN INONESIA REGION Iman Suardi Supervisor: Yui YAGI EE06012 Akio KATSUATA Osamu KAIGAICHI ABSTRACT We proposed a magnitude determination method
More informationExtending the magnitude range of seismic reservoir monitoring by Utilizing Hybrid Surface Downhole Seismic Networks
Extending the magnitude range of seismic reservoir monitoring by Utilizing Hybrid Surface Downhole Seismic Networks Gisela Viegas*, ESG, Kingston, Ontario, Canada Gisela.Fernandes@esgsolutions.com and
More informationThe Combined Inversion of Seismic and Geodetic Data for the Source Process of the 16 October 1999 M w 7.1 Hector Mine, California, Earthquake
Bulletin of the Seismological Society of America, Vol. 92, No. 4, pp. 12 128, May 22 The Combined Inversion of Seismic and Geodetic Data for the Source Process of the October 1999 M w 7.1 Hector Mine,
More informationNumerical Modeling for Earthquake Source Imaging: Implications for Array Design in Determining the Rupture Process
TAO, Vol. 15, No. 2, 133-150, June 2004 Numerical Modeling for Earthquake Source Imaging: Implications for Array Design in Determining the Rupture Process Yi-Ling Huang 1, *, Bor-Shouh
More informationCo-seismic slip from the July 30, 1995, M w 8.1 Antofagasta, Chile, earthquake as constrained by InSAR and GPS observations
98 Chapter 3 Co-seismic slip from the July 30, 1995, M w 8.1 Antofagasta, Chile, earthquake as constrained by InSAR and GPS observations Published by Blackwell Publishing Ltd. in Geophysical Journal International
More informationWas the February 2008 Bukavu seismic sequence associated with magma intrusion?
Was the February 2008 Bukavu seismic sequence associated with magma intrusion? N. d Oreye, C. Lucas National Museum of Natural History, Luxembourg nicolas.doreye@ecgs.lu P. Gonzalez, J. Fernández Univ.
More informationA non-linear geodetic data inversion using ABIC for slip distribution on a fault with an unknown dip angle
Geophys. J. Int. () 173, 353 36 doi: 1.1111/j.1365-6X.7.3713.x A non-linear geodetic data inversion using ABIC for slip distribution on a fault with an unknown dip angle Yukitoshi Fukahata 1 and Tim J.
More informationInterseismic strain accumulation across the Manyi fault (Tibet) prior to the 1997 M w 7.6 earthquake
GEOPHYSICAL RESEARCH LETTERS, VOL. 38,, doi:10.1029/2011gl049762, 2011 Interseismic strain accumulation across the Manyi fault (Tibet) prior to the 1997 M w 7.6 earthquake M. A. Bell, 1 J. R. Elliott,
More informationSource model of the 2005 Miyagi-Oki, Japan, earthquake estimated from broadband strong motions
Earth Planets Space, 59, 1155 1171, 2007 Source model of the 2005 Miyagi-Oki, Japan, earthquake estimated from broadband strong motions Wataru Suzuki and Tomotaka Iwata Disaster Prevention Research Institute,
More informationSendai Earthquake NE Japan March 11, Some explanatory slides Bob Stern, Dave Scholl, others updated March
Sendai Earthquake NE Japan March 11, 2011 Some explanatory slides Bob Stern, Dave Scholl, others updated March 14 2011 Earth has 11 large plates and many more smaller ones. Plates are 100-200 km thick
More informationA systematic measurement of shear wave anisotropy in southern California Report for SCEC Award #15081 Submitted March 15, 2016
A systematic measurement of shear wave anisotropy in southern California Report for SCEC Award #15081 Submitted March 15, 2016 Investigators: Zhigang Peng (Georgia Tech) I. Project Overview... i A. Abstract...
More informationSource Characteristics of Large Outer Rise Earthquakes in the Pacific Plate
Source Characteristics of Large Outer Rise Earthquakes in the Pacific Plate T. Sasatani, N. Takai, M. Shigefuji, and Y. Miyahara Hokkaido University, Sapporo, Japan W. Kawabata Electric Power Development
More informationRapid magnitude determination from peak amplitudes at local stations
Earth Planets Space, 65, 843 853, 2013 Rapid magnitude determination from peak amplitudes at local stations Akio Katsumata 1, Hiroshi Ueno 1, Shigeki Aoki 1, Yasuhiro Yoshida 2, and Sergio Barrientos 3
More information2008 Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
STRUCTURE OF THE KOREAN PENINSULA FROM WAVEFORM TRAVEL-TIME ANALYSIS Roland Gritto 1, Jacob E. Siegel 1, and Winston W. Chan 2 Array Information Technology 1 and Harris Corporation 2 Sponsored by Air Force
More information