Dynamic source inversion for physical parameters controlling the 2017 Lesvos earthquake
|
|
- Evelyn Snow
- 5 years ago
- Views:
Transcription
1 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 Physics, Dept. of Geophysics, Czech Republic 2 University of Patras, Department of Geology, Laboratory of Seismology, Greece
2 The 12 th June, 2017 Lesvos earthquake Aegean sea, south of Lesvos Normal faulting mechanism, Mw 6.3 Hypocentral depth (GI NOA): (12.0 ± 1.7) km Centroidal depth: 8 km (3D Grid search, ISOLA, regional data) SW-dipping fault (40 ) consisent with hypocenter and centroid Multiple point source inversion dominated by a single subevent. The closest village to the epicenter was Plomari, but the most struck was Vrisa (1 casualty, heavy structural damage). USGS Red diamond: hypocenter Beachball at the centroid Blue stars: major aftershocks (M>3.9)
3 Kinematic inversion of the Lesvos earthquake (Kiratzi, 2018) Broad band and strong motion data, Hz Slip confined in a shallow patch, 20 km 10 km along strike and dip, respectively. Mainly unilateral propagation towards WNW Rupture speed: 3.1 km/s Rupture duration: ~6.5 s. M 0 : Nm (MW 6.35) Stress drop: 3.6 MPa Kiratzi, 2018
4 Method dynamic inversion
5 Dynamic source inversion Slip Slip rate Green s functions Hypocenter Fault plane Distribution of prestress frictional parameters Compare with observed data and generate new models (Monte Carlo) Dynamic rupture simulation Friction law E.g. Peyrat and Olsen (2004), Ruiz and Madariaga (2011), Twardzik et al. (2014) Evolution of slip
6 Linear slip-weakening friction (Ida, 1972) 1. Zero slip as long as traction T is smaller than the strength, T u (rupture criterion) 2. When the traction reaches T u, it starts evolving according to: T D = T u 1 D D c for 0 < D < D C T D = 0 otherwise. where D Slip T- Friction T u - Upper yield stress (strength) D C critical slip distance G C = T 0 D C /2 fracture energy density
7 Linear slip-weakening friction (Ida, 1972) 1. Zero slip as long as traction T is smaller than the strength, T u (rupture criterion) 2. When the traction reaches T u, it starts evolving according to: T D = T u 1 D D c for 0 < D < D C T D = 0 otherwise. where D Slip T- Friction T u - Upper yield stress (strength) D C critical slip distance G C = T 0 D C /2 fracture energy density
8 Linear slip-weakening friction (Ida, 1972) 1. Zero slip as long as traction T is smaller than the (strength), T u (rupture criterion) 2. When the traction reaches T u, it starts evolving according to: T D = T u 1 D D c for 0 < D < D C T D = 0 otherwise. where D Slip T- Friction T u - Upper yield stress (strength) D C critical slip distance G C = T 0 D C /2 fracture energy density
9 Linear slip-weakening friction (Ida, 1972) 1. Zero slip as long as traction T is smaller than the strenght, T u (rupture criterion) 2. When the traction reaches T u, it starts evolving according to: T D = T u 1 D D c for 0 < D < D C T D = 0 otherwise. where D Slip T- Friction T u - Upper yield stress (strength) D C critical slip distance G C = T 0 D C /2 fracture energy density
10 Forward solver (need for speed) Dynamic rupture by FD3D code (Madariaga et al., 1998): Finite differences (4 th order) on a Cartesian box 200m grid Vertical fault reaching a free surface Boundary conditions Initial conditions: displacement and traction zero everywhere, except for the fault (which has nonzero prestress and a nucleation zone with T=T nucl >T u ). Symmetry conditions permit to solve the problem on half of the domain Box covers only the fault only slip rates are saved Wave propagation by Axitra (Bouchon, 1981; Coutant, 1989) 1D layered velocity model Pre-calculated Green s functions on a coarser grid, respecting the true fault geometry Representation theorem is used to obtain station waveforms
11 Parametrization : a single elliptic patch (barrier) 13 inverted parameters: 8 geometric parameters (5 for the ellipse, 3 for the nucleation zone) 3 parameters to determine initial stress: T 0, T avg prestress inside the ellipse (gaussian shape) T nucl initial stress inside the nucl. zone 2 parameters governing friction: T u Upper yield stress constant D C critical slip consant The fault plane
12 Bayesian inversion by Monte Carlo approach Data errors: Gaussian distribution with constant variance Prior information (constraints) on model parameters: Homogeneous prior PDF D c > 0.10m (for numerical reasons) Nucl. zone not fixed Posterior PDF is sampled using Monte Carlo Markov Chain by the Parallel tempering method (Sambridge, 2013)
13 Application to the Lesvos earthquake
14 Data, the model fault plane and the velocity model 21 accelerograms from Greek and Turkish national seismic networks. Distance of stations from the centroid: km. Integrated to displacements, processed with causal Butterworth s filter of fourth-order in the range of Hz. The velocity model 0 5 Plot of Vp, Vs Karagianni 5 th=4thlayer Vp Vs Depth (km) Velocity (km/sec) Karagianni et al., 2002
15 Results: the best found model (out of ~10 6 ) Quantity Value T u 13 Mpa T nucl 17 Mpa R nucl D C Stress drop 1.3 km 0.25 m 6.2 MPa M 0 3 x Nm Duration of rupture 8.2 s Avg. rupture velocity 1.1 km/s Variance Reduction 37.2 %
16 Optimal model: evoluton of slip rate
17 Comparison with kinematic inversion Slip from optimal dynamic model Kinematic model (Kiratzi, 2018)
18 Comparison with kinematic inversion Rupture velocity: 1200 m/s vs 3000 m/s Stress drop: 6 MPa vs 3.6 MPa Slip from optimal dynamic model Kinematic model (Kiratzi, 2018)
19 esults : fits to data Variance reduction: VR 1 i 2 s i o i 2 2 = 37 % o i Frequency range: Hz 2 Observed data Simulated data
20 Uncertainties: the stress and friction parameters (1-D marginals)
21 Correlations between dynamic parameters (2-D marginals)
22 Correlations between dynamic parameters (2-D marginals)
23 Correlations between dynamic parameters (2-D marginals)
24 Uncertainties in geometric parameters Fault plane Blue : contours of slip Yellow-orange: Nucleation zones Green star: centroid Red star: hypocenter
25 Histograms of various rupture properties (1-D marginals)
26 Histograms of various rupture properties (1-D marginals) M 0 = 3.5 ± N.m Kiratzi: N.m Papadimitriou et. al.: N.m
27 Histograms of various rupture properties (1-D marginals) V r = 1.1 ± 0.08 km/s Kiratzi: 3.2 km/s Δσ = 5.8 ± 0.8 MPa Kiratzi: 3. 6 MPa
28 V r = 1.1 ± 0.08 km/s Kiratzi: 3.2 km/s Δσ = 5.8 ± 0.8 MPa Kiratzi: 3. 6 MPa
29 Conclusion We use Parallel tempering to obtain a dynamic model along with estimates of uncertainties and correlations. Our dynamic model fits data with 37 % variance reduction in the frequency range of Hz M 0 consistent (within infered std) with other researchers Slow (1 km/s) rupture propagation to the northwest High (6 MPa) stress drop Centroid location within nucleation zone uncertainty (unlike the hypocenter) Some parameters are correlated (e.g. T 0 and Dc, R and Dc) and can be varied together to obtain the same fit to data.
Lesvos June 12, 2017, Mw 6.3 event, a quick study of the source
Lesvos June 12, 2017, Mw 6.3 event, a quick study of the source E. Sokos 1) and J. Zahradník 2) 1) University of Patras, Greece 2) Charles University, Prague, Czech Republic Report sent to EMSC on 21/06/2017
More informationPreliminary slip model of M9 Tohoku earthquake from strongmotion stations in Japan - an extreme application of ISOLA code.
Preliminary slip model of M9 Tohoku earthquake from strongmotion stations in Japan - an extreme application of ISOLA code. J. Zahradnik 1), F. Gallovic 1), E. Sokos 2) G-A. Tselentis 2) 1) Charles University
More informationThe Mw 6.2 Leonidio, southern Greece earthquake of January 6, 2008: Preliminary identification of the fault plane.
The Mw 6.2 Leonidio, southern Greece earthquake of January 6, 28: Preliminary identification of the fault plane. J. Zahradnik 1, E. Sokos 2, A.Serpetsidaki 2, and G A.Tselentis 2 1 Charles University in
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 informationSCENARIO MODELING OF THE 2014 Mw6.0 SOUTH NAPA, CALIFORNIA, EARTHQUAKE USING AN ADVANCED BROADBAND KINEMATIC SOURCE MODEL
1 Best Practices in Physics-based Fault Rupture Models for SCENARIO MODELING OF THE 2014 Mw6.0 SOUTH NAPA, CALIFORNIA, EARTHQUAKE USING AN ADVANCED BROADBAND KINEMATIC SOURCE MODEL F. GALLOVIČ 1 1 Charles
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 informationInverting full waveforms into 1D seismic velocity model of the upper crust by neighborhood algorithm - Corinth Gulf, Greece
Inverting full waveforms into 1D seismic velocity model of the upper crust by neighborhood algorithm - Corinth Gulf, Greece VLADIMÍR PLICKA AND JIŘÍ ZAHRADNÍK Department of Geophysics, Faculty of Mathematics
More information19 years of Research Infrastructure PSLNET, significant influence to study the source parameters of Greece earthquakes.
19 years of Research Infrastructure PSLNET, significant influence to study the source parameters of Greece earthquakes. V. Plicka Charles university Prague Outline History Projects Publications Software
More informationZakynthos 25/10/2018, Mw 6.8 earthquake: Superposition of strike-slip and thrust?
Zakynthos 25/10/2018, Mw 6.8 earthquake: Superposition of strike-slip and thrust? J. Zahradník 1), E. Sokos 2), V. Plicka 1) 1) Charles University, Prague, Czech Republic 2) University of Patras, Greece
More informationTeleseismic waveform modelling of the 2008 Leonidio event
The 6 January 2008 (Mw6.2) Leonidio (southern Greece) intermediate depth earthquake: teleseismic body wave modelling Anastasia Kiratzi and Christoforos Benetatos Department of Geophysics, Aristotle University
More informationON NEAR-FIELD GROUND MOTIONS OF NORMAL AND REVERSE FAULTS FROM VIEWPOINT OF DYNAMIC RUPTURE MODEL
1 Best Practices in Physics-based Fault Rupture Models for Seismic Hazard Assessment of Nuclear ON NEAR-FIELD GROUND MOTIONS OF NORMAL AND REVERSE FAULTS FROM VIEWPOINT OF DYNAMIC RUPTURE MODEL Hideo AOCHI
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 informationTitle from Slip-Velocity Functions on Ear.
Title Stress-Breakdown Time and Slip-Weak from Slip-Velocity Functions on Ear Author(s) Mikumo, Takeshi; Olsen, Kim Yuji B.; Fuk Citation Bulletin of the Seismological Socie 9(): 6-8 Issue Date - URL http://hdl.handle.net//9
More informationDynamic stress field of a kinematic earthquake source model with k-squared slip distribution
Geophys. J. Int. (2007) doi: 10.1111/j.1365-246X.2007.03548.x Dynamic stress field of a kinematic earthquake source model with k-squared slip distribution Jan Burjánek and Jiří Zahradník Charles University
More informationSTRIKE SLIP SPLAY USING DYNAMIC RUPTURE MODELS
STRIKE SLIP SPLAY USING DYNAMIC RUPTURE MODELS Julian Lozos (PEER, UC Berkeley) SWUS GMC Workshop #2, Berkeley, CA, 24 October 2013 Outline Issues in the dynamics of fault branches. Advantages and disadvantages
More informationChallenges in earthquake physics and source imaging
Challenges in earthquake physics and source imaging Jean-Paul Ampuero and Nadia Lapusta (Caltech Seismolab) Main goals and current issues in earthquake dynamics The source imaging inverse problem Parallels
More informationISOLA code for multiple-point source modeling review
In: Moment Tensor Solutions - A Useful Tool for Seismotectonics (edited by Sebastiano D'Amico), Springer ISOLA code for multiple-point source modeling review Jiří Zahradník 1 and Efthimios Sokos 2 1 Charles
More informationDistinguishing barriers and asperities in near-source ground motion
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 110,, doi:10.1029/2005jb003736, 2005 Distinguishing barriers and asperities in near-source ground motion Morgan T. Page, Eric M. Dunham, 1 and J. M. Carlson Department
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (016) 07, 101 109 Advance Access publication 016 August 5 GJI Seismology doi: 10.1093/gji/ggw30 Fast and cheap approximation of Green function uncertainty
More informationSTUDYING THE IMPORTANT PARAMETERS IN EARTHQUAKE SIMULATION BASED ON STOCHASTIC FINITE FAULT MODELING
STUDYING THE IMPORTANT PARAMETERS IN EARTHQUAKE SIMULATION BASED ON STOCHASTIC FINITE FAULT MODELING H. Moghaddam 1, N. Fanaie 2* and H. Hamzehloo 1 Professor, Dept. of civil Engineering, Sharif University
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 informationBulletin of the Geological Society of Greece
Bulletin of the Geological Society of Greece Vol. 43, 2010 STUDY OF THE 2ND DECEMBER 2002 VARTHOLOMIO EARTHQUAKE (WESTERN PELOPONNESE) M5.5 AFTERSHOCK SEQUENCE Serpetsidaki A. University of Patras, Department
More informationEstimation of Peak Ground Acceleration for Delhi Region using Finsim, a Finite Fault Simulation Technique
215 Estimation of Peak Ground Acceleration for Delhi Region using Finsim, a Finite Fault Simulation Technique NEELIMA SATYAM. D* and K. S. RAO** * Earthquake Engineering Research Centre, International
More informationOutstanding Problems. APOSTOLOS S. PAPAGEORGIOU University of Patras
NEAR-FAULT GROUND MOTIONS: Outstanding Problems APOSTOLOS S. PAPAGEORGIOU University of Patras Outline Characteristics of near-fault ground motions Near-fault strong ground motion database A mathematical
More informationBayesian ISOLA: new tool for automated centroid moment tensor inversion
Bayesian ISOLA: new tool for automated centroid moment tensor inversion A CMT catalog automatically generated for Swiss data Jiří Vackář 1,2 Jan Burjánek 2 František Gallovič 1 Jiří Zahradník 1 John Clinton
More informationDynamic stress field of a kinematic earthquake source model
Geophys. J. Int. (2000) 142, 000 000 Dynamic stress field of a kinematic earthquake source model with k-squared slip distribution Jan Burjánek 1 and Jiří Zahradník 1 1 Charles University in Prague, Faculty
More informationHeterogeneous Coulomb stress perturbation during earthquake cycles in a 3D rate-and-state fault model
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L21306, doi:10.1029/2008gl035614, 2008 Heterogeneous Coulomb stress perturbation during earthquake cycles in a 3D rate-and-state fault
More informationBroadband Strong Motion Simulation Based on a Dynamic Fault Rupture Model: Application to the 2000 Tottori, Japan earthquake
Workshop on Numerical odeling of Earthquake Source Dynamics, Smolenice Castle, Slovak Republic Broadband Strong otion Simulation Based on a Dynamic Fault Rupture odel: Application to the 2000 Tottori,
More informationSource Process and Constitutive Relations of the 2011 Tohoku Earthquake Inferred from Near-Field Strong-Motion Data
Source Process and Constitutive Relations of the 2011 Tohoku Earthquake Inferred from Near-Field Strong-Motion Data Kunikazu Yoshida, Anatoly Petukhin & Ken Miyakoshi Geo-Research Institute, Japan Koji
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting Lecture 20, 30 Nov. 2017 www.geosc.psu.edu/courses/geosc508 Seismic Spectra & Earthquake Scaling laws. Seismic Spectra & Earthquake Scaling laws. Aki, Scaling law
More informationGeometry of co-seismic surface ruptures and tectonic meaning of the 23 October 2011 Mw 7.1 Van earthquake (East Anatolian Region, Turkey)
Geometry of co-seismic surface ruptures and tectonic meaning of the 23 October 2011 Mw 7.1 Van earthquake (East Anatolian Region, Turkey) Bülent Dogan & Ahmet Karakaș. Journal of Structural Geology, 46(2013),
More informationResolvability of Isotropic Component in Regional Seismic Moment Tensor Inversion
Bulletin of the Seismological Society of America, Vol. 103, No. 4, pp. 2460 2473, August 2013, doi: 10.1785/0120120097 Resolvability of Isotropic Component in Regional Seismic Moment Tensor Inversion by
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 informationSurface Rupture in Kinematic Ruptures Models for Hayward Fault Scenario Earthquakes
Surface Rupture in Kinematic Ruptures Models for Hayward Fault Scenario Earthquakes Brad Aagaard May 21, 2009 Surface Rupture from Kinematic Rupture Models Objective Holistic approach to calculating surface
More informationTOWARD A JOINT CATALOGUE OF RECENT SEISMICITY IN WESTERN GREECE: PRELIMINARY RESULTS
Δελτίο της Ελληνικής Γεωλογικής Εταιρίας τομ. XXXVII, 2007 Πρακτικά 11 ου Διεθνούς Συνεδρίου, Αθήνα, Μάιος 2007 Bulletin of the Geological Society of Greece vol. XXXVII, 2007 Proceedings of the 11 th International
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 informationFault Representation Methods for Spontaneous Dynamic Rupture Simulation. Luis A. Dalguer
Fault Representation Methods for Spontaneous Dynamic Rupture Simulation Luis A. Dalguer Computational Seismology Group Swiss Seismological Service (SED) ETH-Zurich July 12-18, 2011 2 st QUEST Workshop,
More informationModeling the M L 4.7 mainshock of the February-July 2001 earthquake sequence in Aegion, Greece
Journal of Seismology 8: 247 257, 2004. 2004 Kluwer Academic Publishers. Printed in the Netherlands. 247 Modeling the M L 4.7 mainshock of the February-July 2001 earthquake sequence in Aegion, Greece J.
More informationJoint inversion of InSAR and broadband teleseismic waveform data with ABIC: application to the 1997 Manyi, Tibet earthquake
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
More informationArthur Frankel, William Stephenson, David Carver, Jack Odum, Robert Williams, and Susan Rhea U.S. Geological Survey
Probabilistic Seismic Hazard Maps for Seattle: 3D Sedimentary Basin Effects, Nonlinear Site Response, and Uncertainties from Random Velocity Variations Arthur Frankel, William Stephenson, David Carver,
More informationChapter 2. Earthquake and Damage
EDM Report on the Chi-Chi, Taiwan Earthquake of September 21, 1999 2.1 Earthquake Fault 2.1.1 Tectonic Background The island of Taiwan is located in the complex junction where the Eurasian and Philippine
More informationVariability of Near-Field Ground Motion from Dynamic Earthquake Rupture Simulations
Variability of Near-Field Ground Motion from Dynamic Earthquake Rupture Simulations J. Ripperger, P.M. Mai, J.-P. Ampuero Institute of Geophysics, ETH Zurich, 893 Zurich, Switzerland Abstract. This study
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 informationTwo ways to think about the dynamics of earthquake ruptures
Two ways to think about the dynamics of earthquake ruptures (1) In terms of friction (2) In terms of fracture mechanics Scholz describes conditions for rupture propagation (i.e. instability) via energy
More informationTeleseismic Body Waves from Dynamically Rupturing Shallow Thrust Faults: Are They Opaque for Surface-Reflected Phases?
Bulletin of the Seismological Society of America, Vol. 95, No. 3, pp. 8 87, June 25, doi:.785/237 Teleseismic Body Waves from Dynamically Rupturing Shallow Thrust Faults: Are They Opaque for Surface-Reflected
More informationHybrid k-squared source model for strong ground motion simulations: Introduction
Physics of the Earth and Planetary Interiors 160 (2007) 34 50 Hybrid k-squared source model for strong ground motion simulations: Introduction František Gallovič, Johana Brokešová Department of Geophysics,
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 informationModelling Strong Ground Motions for Subduction Events in the Wellington Region, New Zealand
Proceedings of the Ninth Pacific Conference on Earthquake Engineering Building an Earthquake-Resilient Society 14-16 April, 2011, Auckland, New Zealand Modelling Strong Ground Motions for Subduction Events
More informationSeismic slip on an upper-plate normal fault during a large subduction megathrust rupture
SUPPLEMENTARY INFORMATION DOI: 10.1038/NGEO2585 Seismic slip on an upper-plate normal fault during a large subduction megathrust rupture Stephen P. Hicks* and Andreas Rietbrock NATURE GEOSCIENCE www.nature.com/naturegeoscience
More informationWhich Dynamic Rupture Parameters Can Be Estimated from Strong Ground Motion and Geodetic Data?
Pure appl. geophys. 161 (24) 2155 2169 33 4553/4/122155 15 DOI 1.17/s24-4-2555-9 Ó Birkhäuser Verlag, Basel, 24 Pure and Applied Geophysics Which Dynamic Rupture Parameters Can Be Estimated from Strong
More informationAutomated procedure for point and kinematic source inversion at regional distances
Click Here for Full Article JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.1029/2009jb006450, 2010 Automated procedure for point and kinematic source inversion at regional distances Simone Cesca, 1,2
More informationBayesian ISOLA: new tool for automated centroid moment tensor inversion
Bayesian ISOLA: new tool for automated centroid moment tensor inversion Jiří Vackář 1,2,3 Jan Burjánek 1,3 František Gallovič 2 Jiří Zahradník 2 1 Swiss Seismological Service, ETH Zurich 2 Charles University
More informationJournalofGeophysicalResearch: SolidEarth
JournalofGeophysicalResearch: SolidEarth RESEARCH ARTICLE Key Points: Relations between source dimensions and corner frequencies are established Stress drop estimates are affected by differences in source
More informationMagnitude 8.2 NORTHWEST OF IQUIQUE, CHILE
An 8.2-magnitude earthquake struck off the coast of northern Chile, generating a local tsunami. The USGS reported the earthquake was centered 95 km (59 miles) northwest of Iquique at a depth of 20.1km
More informationSynthetic Near-Field Rock Motions in the New Madrid Seismic Zone
Synthetic Near-Field Rock Motions in the New Madrid Seismic Zone Genda Chen*, Ph.D., P.E., and Mostafa El-Engebawy Engebawy,, Ph.D. *Associate Professor of Civil Engineering Department of Civil, Architecture
More informationof other regional earthquakes (e.g. Zoback and Zoback, 1980). I also want to find out
4. Focal Mechanism Solutions A way to investigate source properties of the 2001 sequence is to attempt finding well-constrained focal mechanism solutions to determine if they are consistent with those
More informationCOULOMB STRESS CHANGES DUE TO RECENT ACEH EARTHQUAKES
COULOMB STRESS CHANGES DUE TO RECENT ACEH EARTHQUAKES Madlazim Physics Department, Faculty Mathematics and Sciences of Surabaya State University (UNESA) Jl. Ketintang, Surabaya 60231, Indonesia. e-mail:
More informationCase Study 2: 2014 Iquique Sequence
Case Study 2: 2014 Iquique Sequence Overview Mw 8.2 earthquake on 1 April 2014 at 11:08:43 UTC Pictures of damage Seismicity maps Foreshock seismicity -> Main shock -> Aftershock 1 April 2014 11:08:43
More informationSlip-weakening behavior during the propagation of dynamic ruptures obeying rate- and state-dependent friction laws
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 108, NO. B8, 2373, doi:10.1029/2002jb002198, 2003 Slip-weakening behavior during the propagation of dynamic ruptures obeying rate- and state-dependent friction laws
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 informationReceived 15 January 2010, Revised 23 August 2010, Accepted 17 September 2010
Indonesian Journal of Physics Vol 21 No. 4, October 2010 Seismic Moment Tensors of Padang 20090930 and Jambi 20091001 Events in Jambi using Three Component Local Waveforms: Identification of the Active
More informationThe Mw November 2015 Lefkada (Greece) Earthquake: Structural Interpretation by Means of the Aftershock Analysis
Pure Appl. Geophys. Ó 2017 Springer International Publishing AG DOI 10.1007/s00024-017-1601-3 Pure and Applied Geophysics The Mw6.5 17 November 2015 Lefkada (Greece) Earthquake: Structural Interpretation
More informationWhat happened before the last five strong earthquakes in Greece: Facts and open questions
86 Proc. Jpn. Acad., Ser. B 82 (2006) [Vol. 82, What happened before the last five strong earthquakes in Greece: Facts and open questions By Panayiotis A. VAROTSOS ) Solid State Section and Solid Earth
More informationHigh-Frequency Ground Motion Simulation Using a Source- and Site-Specific Empirical Green s Function Approach
High-Frequency Ground Motion Simulation Using a Source- and Site-Specific Empirical Green s Function Approach R. Mourhatch & S. Krishnan California Institute of Technology, Pasadena, CA, USA SUMMARY: A
More informationFault Specific, Dynamic Rupture Scenarios for Strong Ground Motion Prediction
Fault Specific, Dynamic Rupture Scenarios for Strong Ground Motion Prediction H. Sekiguchi Disaster Prevention Research Institute, Kyoto University, Japan Blank Line 9 pt Y. Kase Active Fault and Earthquake
More informationTopography on Earthquake Motions in Sedimentary Basins
TRANSPORTATION RESEARCH RECORD 1411 Effects of Three-Dimensional Bedrock Topography on Earthquake Motions in Sedimentary Basins ARTHUR FRANKEL Work being done at the U.S. Geological Survey on 3-D simulations
More informationCalculation of Focal mechanism for Composite Microseismic Events
Calculation of Focal mechanism for Composite Microseismic Events Hongliang Zhang, David W. Eaton Department of Geoscience, University of Calgary Summary It is often difficult to obtain a reliable single-event
More informationIntegration of Seismic and Seismological Data Interpretation for Subsurface Structure Identification
Integration of Seismic and Seismological Data Interpretation for Subsurface Structure Identification Iftikhar Ahmed Satti, and Wan Ismail Wan Yusoff Abstract The structural interpretation of a part of
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 informationSpectral Element simulation of rupture dynamics
Spectral Element simulation of rupture dynamics J.-P. Vilotte & G. Festa Department of Seismology, Institut de Physique du Globe de Paris, 75252 France ABSTRACT Numerical modeling is an important tool,
More informationRUPTURE MODELS AND GROUND MOTION FOR SHAKEOUT AND OTHER SOUTHERN SAN ANDREAS FAULT SCENARIOS
RUPTURE MODELS AND GROUND MOTION FOR SHAKEOUT AND OTHER SOUTHERN SAN ANDREAS FAULT SCENARIOS ABSTRACT : L.A. Dalguer 1, S.M. Day 2, K. Olsen 3, and V.M. Cruz-Atienza 4 1 Institute of Geophysics, ETH-Zurich,
More informationThe 2016, October 26, Central Italy Earthquake Origin Time 17:10:36 UTC, M L(ISNet) =5.3; M W(ISNet) =5.6
The 2016, October 26, Central Italy Earthquake Origin Time 17:10:36 UTC, M L(ISNet) =5.3; M W(ISNet) =5.6 RISSC-Lab: Laboratorio di RIcerca in Sismologia Sperimentale e Computazionale The event as seen
More informationToward understanding slip inversion uncertainty and artifacts
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 115,, doi:10.1029/2010jb007414, 2010 Toward understanding slip inversion uncertainty and artifacts Jiří Zahradník 1 and František Gallovič 1,2 Received 26 January
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 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 informationTriggering of earthquakes during the 2000 Papua New Guinea earthquake sequence
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 112,, doi:10.1029/2006jb004480, 2007 Triggering of earthquakes during the 2000 Papua New Guinea earthquake sequence Sun-Cheon Park 1 and Jim Mori 1 Received 3 May
More informationEarthquake stress drop estimates: What are they telling us?
Earthquake stress drop estimates: What are they telling us? Peter Shearer IGPP/SIO/U.C. San Diego October 27, 2014 SCEC Community Stress Model Workshop Lots of data for big earthquakes (rupture dimensions,
More informationEarthquake Stress Drops in Southern California
Earthquake Stress Drops in Southern California Peter Shearer IGPP/SIO/U.C. San Diego September 11, 2009 Earthquake Research Institute Lots of data for big earthquakes (rupture dimensions, slip history,
More informationDi#erences in Earthquake Source and Ground Motion Characteristics between Surface and Buried Crustal Earthquakes
Bull. Earthq. Res. Inst. Univ. Tokyo Vol. 2+,**0 pp.,/3,00 Di#erences in Earthquake Source and Ground Motion Characteristics between Surface and Buried Crustal Earthquakes Paul Somerville* and Arben Pitarka
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 informationSeismic Efficiency, Overshoot and Enhanced Dynamic Weaking of Fractures Associated with Stimulation in Heavy Oil Reservoirs
Seismic Efficiency, Overshoot and Enhanced Dynamic Weaking of Fractures Associated with Stimulation in Heavy Oil Reservoirs Lindsey N. Meighan, Ted Urbancic and Adam Baig ESG, Kingston, Ontario Summary
More informationModel Inversion for Induced Seismicity
Model Inversion for Induced Seismicity David Castiñeira Research Associate Department of Civil and Environmental Engineering In collaboration with Ruben Juanes (MIT) and Birendra Jha (USC) May 30th, 2017
More informationSource parameters II. Stress drop determination Energy balance Seismic energy and seismic efficiency The heat flow paradox Apparent stress drop
Source parameters II Stress drop determination Energy balance Seismic energy and seismic efficiency The heat flow paradox Apparent stress drop Source parameters II: use of empirical Green function for
More informationIterative Deconvolution of Regional Waveforms and a Double-Event Interpretation of the 2003 Lefkada Earthquake, Greece
Bulletin of the Seismological Society of America, Vol. 95, No. 1, pp. 159 172, February 2005, doi: 10.1785/0120040035 Iterative Deconvolution of Regional Waveforms and a Double-Event Interpretation of
More informationSupplementary Materials for
advances.sciencemag.org/cgi/content/full/4/3/eaao4915/dc1 Supplementary Materials for Global variations of large megathrust earthquake rupture characteristics This PDF file includes: Lingling Ye, Hiroo
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 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 informationInterpreting lab creep experiments: from single contacts... to interseismic fault models
Interpreting lab creep experiments: from single contacts...... to interseismic fault models Delphine Fitzenz, EOST Strasbourg Steve Hickman, U.S.G.S. Menlo Park Yves Bernabé, EOST, now at MIT Motivation
More informationSupporting information for the main manuscript: Strain budget of the Ecuador-Colombia subduction zone: a stochastic view
Supporting information for the main manuscript: Strain budget of the Ecuador-Colombia subduction zone: a stochastic view B. Gombert 1, Z. Duputel 1, R. Jolivet 2, M. Simons 3, J. Jiang 4, C. Liang 5, E.
More informationModeling velocity recordings of the Mw6.0 South Napa, California, earthquake: unilateral event with weak high-frequency directivity
Modeling velocity recordings of the Mw6.0 South Napa, California, earthquake: unilateral event with weak high-frequency directivity František Gallovič 1 1 Charles University, Faculty of Mathematics and
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 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 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 informationSPATIAL DISTRIBUTION OF STRONG GROUND MOTION CONSIDERING ASPERITY AND DIRECTIVITY OF FAULT
SPATIAL DISTRIBUTION OF STRONG GROUND MOTION CONSIDERING ASPERITY AND DIRECTIVITY OF FAULT Shunroku YAMAMOTO SUMMARY Waveform simulations of the 995 Hyogo-ken Nanbu earthquake were carried out to study
More informationAn analysis of the Kefalonia seismic sequence of Jan Feb. 3, 2014
An analysis of the Kefalonia seismic sequence of Jan. 26 - Feb. 3, 2014 Alessandro Caporali 1),Carine Bruyninx 2),Rui Fernandes 3),Athanassios Ganas 4),Ambrus Kenyeres 5),Martin Lidberg 6), Guenter Stangl
More informationRELATION BETWEEN RAYLEIGH WAVES AND UPLIFT OF THE SEABED DUE TO SEISMIC FAULTING
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 24 Paper No. 1359 RELATION BETWEEN RAYLEIGH WAVES AND UPLIFT OF THE SEABED DUE TO SEISMIC FAULTING Shusaku INOUE 1,
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 informationNumerical study on multi-scaling earthquake rupture
GEOPHYSICAL RESEARCH LETTERS, VOL. 31, L266, doi:1.129/23gl1878, 24 Numerical study on multi-scaling earthquake rupture Hideo Aochi Institut de Radioprotection et de Sûreté Nucléaire, France Satoshi Ide
More informationEffects of Fault Dip and Slip Rake Angles on Near-Source Ground Motions: Why Rupture Directivity Was Minimal in the 1999 Chi-Chi, Taiwan, Earthquake
Bulletin of the Seismological Society of America, Vol. 94, No. 1, pp. 155 170, February 2004 Effects of Fault Dip and Slip Rake Angles on Near-Source Ground Motions: Why Rupture Directivity Was Minimal
More informationSeismic and aseismic processes in elastodynamic simulations of spontaneous fault slip
Seismic and aseismic processes in elastodynamic simulations of spontaneous fault slip Most earthquake simulations study either one large seismic event with full inertial effects or long-term slip history
More information