Mathematical modelling of tsunami wave generation

Size: px
Start display at page:

Download "Mathematical modelling of tsunami wave generation"

Transcription

1 Mathematical modelling of tsunami wave generation DENYS DUTYKH 1 Chargé de Recherche CNRS 1 CNRS-LAMA, Université de Savoie Campus Scientifique Le Bourget-du-Lac, France Institut Jean le Rond d Alembert DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 1 / 53

2 Acknowledgements Collaborators: Frédéric Dias: Professor, University College of Dublin Youen Kervella: PhD student at IFREMER, Brest Special thanks to: Costas Synolakis: Professor, Univeristy of Southern California Emile Okal: Professor, Northwestern University DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 2 / 53

3 Outline 1 Introduction 2 Tsunami generation 3 Dislocation dynamics 4 Sediments effect DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 3 / 53

4 Outline 1 Introduction 2 Tsunami generation 3 Dislocation dynamics 4 Sediments effect DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 4 / 53

5 Why study tsunami generation? Main motivation for this research Objectives: 1 Understand seismology/ hydrodynamics coupling 2 What are main physical effects to take into account? 3 Provide an initial condition for various tsunami propagation codes Figure: Global seismic wave propagation model Difficulties: Almost no information about sea-bed deformations DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 5 / 53

6 Okada s solution Seismic deformation computation Physical assumptions: Isotropic homogeneous half-space Linear elasticity Static solution Fault is a rectangular dislocation Figure: U 1 strike-slip, U 2 dip-slip, U 3 tensile fault Reference Y. Okada. Surface deformation due to shear and tensile faults in a half-space. Bull. Seism. Soc. Am., 1985, 75, DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 6 / 53

7 Okada s solution Three elementary solutions parameter value δ 13 d 40 km L 150 km W 80 km ν 0.27 E Pa U 15 m Figure: Dip-slip deformation Reference: F. Dias, D. Dutykh. Dynamics of tsunami waves. In book: Extreme Man-Made and Natural Hazards in Dynamics of Structures, 2006 DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 7 / 53

8 Okada s solution Three elementary solutions parameter value δ 13 d 40 km L 150 km W 80 km ν 0.27 E Pa U 15 m Figure: Strike-slip deformation Reference: F. Dias, D. Dutykh. Dynamics of tsunami waves. In book: Extreme Man-Made and Natural Hazards in Dynamics of Structures, 2006 DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 7 / 53

9 Okada s solution Three elementary solutions parameter value δ 13 d 40 km L 150 km W 80 km ν 0.27 E Pa U 15 m Figure: Tensile deformation Reference: F. Dias, D. Dutykh. Dynamics of tsunami waves. In book: Extreme Man-Made and Natural Hazards in Dynamics of Structures, 2006 DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 7 / 53

10 Main ingredients of Okada s solution - I Volterra s theory of dislocations Annales scientifiques de l E.N.S., 24, 1907 Displacement field in an elastic medium: u k = 1 b i τ k 8πµ ijν j dσ Σ DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 8 / 53

11 Main ingredients of Okada s solution - II Fundamental solution for an isotropic homogeneous half-space Entire space: Green-Somigliana s tensor Half-space: Mindlin s solution Combination of nuclei of strain in entire space Physics, 7, , 1936 DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 9 / 53

12 Finite fault inversion algorithm Data for December 2004 event: courtesy of Chen Ji, Caltech Figure: Surface projection of the slip distribution Figure: Comparison of data (black) and synthetic seismograms (red) DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 10 / 53

13 Application to the Indian Ocean Tsunami 2004 Initial condition construction (Ref: Grilli et al. (2007)) Figure: Dislocation fault segments Figure: Initial surface elevation DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 11 / 53

14 Indian Ocean Tsunami 2004 Global simulation by V. Titov (NOAA) DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 12 / 53

15 Outline 1 Introduction 2 Tsunami generation 3 Dislocation dynamics 4 Sediments effect DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 13 / 53

16 Traditional approach to tsunami generation - I Passive generation Put coseismic displacements directly on the free surface and let it propagate: DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 14 / 53

17 Traditional approach to tsunami generation - I Passive generation Put coseismic displacements directly on the free surface and let it propagate: DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 14 / 53

18 Tsunami generation process Video illustration by courtesy of J. Sharpe DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 15 / 53

19 Traditional approach to tsunami generation - II A few drawbacks Some obvious shortcomings: Interaction between moving bottom and ocean is neglected Time characteristics of source are not taken into account Initial velocity field is neglected u t=0 = 0 DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 16 / 53

20 Comparison between linear and nonlinear models First minutes of tsunami propagation Main objective: To check the importance of nonlinear effects Frequency dispersion 1 Complete water wave problem 2 Cauchy-Poisson problem 3 Nonlinear Shallow-Water Equations (NSWE) Y. Kervella, D. Dutykh, F. Dias. Comparison between three-dimensional linear and nonlinear tsunami generation models, Theoretical and Computational Fluid Dynamics, 2007 DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 17 / 53

21 Complete water wave problem Free-surface potential flow Governing equations: φ = 0, η t + η φ = φ z, z = η, φ t φ 2 + gη = 0, z = η, φ z + h t + φ h = 0, z = h. Solved by BIEM (PhD thesis of P. Guyenne) Accelerated by Fast Multipole Algorithm Reference: C. Fochesato, F. Dias. (2006), A fast method for nonlinear three-dimensional free-surface waves, Proceedings of the Royal Society of London A 462, DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 18 / 53

22 Cauchy-Poisson problem Linearized water-wave problem Governing equations: φ = 0, η t = φ z, z = η, φ t + gη = 0, z = η, φ z + h t = 0, z = h. Fourier-Laplace transform An analytical solution for separable geometries Reference: D. Dutykh, F. Dias. Water waves generated by a moving bottom. In book Tsunami and nonlinear waves, Kundu Anjan (Editor), 2007 DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 19 / 53

23 Nonlinear Shallow Water Equations (NSWE) Saint-Venant equations Govering equations: η t + ((h+η) u ) = h t u t u 2 + gη = 0 Finite volume scheme FVCF Reference: J.-M. Ghidaglia, A. Kumbaro, G. Le Coq. On the numerical solution to two fluid models via a cell centered finite volume method. Eur. J. Mech. B-Fluids, 20, 2001 D. Dutykh, F. Dias. Tsunami generation by dynamic displacement of sea bed due to dip-slip faulting. Accepted to Mathematics and computers in simulation, 2008 DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 20 / 53

24 Comparison results - I Weakly Dispersive and weakly nonlinear waves ε := a 0 h 0 = , µ 2 := ( h 0 l ) 2 = 10 4, S := ε µ 2 = 5. Figure: Linearized solution, Fully nonlinear, NSWE DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 21 / 53

25 Comparison results - II Dispersive and weakly nonlinear waves ε := a 0 h 0 = , µ 2 := ( h 0 l ) 2 = 0.02, S := ε µ 2 = Tide gauge 1 Tide gauge z, m 0 z, m t, s t, s 0.04 Tide gauge 3 Tide gauge z, m 0 z, m t, s t, s Tide gauge 6 Tide gauge z, m 0 z, m t, s t, s Figure: Linearized solution, Fully nonlinear, NSWE DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 22 / 53

26 Cauchy-Poisson problem Analytical solution for flat bottom General solution by Fourier-Laplace transform: η( x, t) = 1 (2π) 2 R 2 1 cosh( k h) 2πi e i k x µ+i µ i s 2 ζ( k, s) s 2 +ω 2 est ds d k ω 2 = g k tanh( k h) ζ( x, t): is the unknown dynamic sea-bed displacement There is no available analytical solution for fault dynamics Issue: Use static solution and make assumptions about the dynamics Reference: D. Dutykh, F. Dias. Water waves generated by a moving bottom. In book Tsunami and nonlinear waves, Kundu Anjan (Editor), 2007 DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 23 / 53

27 Dynamic sea-bed deformation Main ingredients Variables separation: ζ( x, t) = D( x)t(t) Celebrated Okada s solution provides D( x) Assumptions about time evolution T(t) (dynamic scenarios): Instantaneous Exponential Trigonometric Linear T i (t) = H(t) T e (t) = (1 e αt )H(t) T c (t) = H(t t 0 )+ 1 2 [1 cos(πt t 0 )]H(t 0 t) T l (t) = ( H(t t 0 )+ t t 0 H(t 0 t) ) H(t) Application to linear waves: η( x, t) = 1 (2π) 2 R 2 D( k)e i k x 1 cosh( k h) 2πi µ+i µ i s 2ˆT(s) s 2 +ω 2 est ds d k DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 24 / 53

28 Application to tsunami generation problems - I How large is error in translating sea-bed deformation onto free surface? Passive: deformation translated on free surface η( x, t) = 1 D( (2π) i 2 k)e k x cosωt d k Active: instantaneous scenario η i ( x, t) = 1 D( k)e i k x (2π) 2 cosh( k h) cosωt d k Reference: D. Dutykh, F. Dias, Y. Kervella. Linear theory of wave generation by a moving bottom, CRAS, 343, 2006 R 2 R 2 DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 25 / 53

29 Application to tsunami generation problems - II Numerical results and conclusions z, m Active generation Passive generation time (s) z, m time (s) η active η passive / η active time (s) η active η passive / η active time (s) z, m time (s) z, m time (s) η active η passive / η active time (s) η active η passive / η active time (s) Figure: Tide gauge records Figure: Relative difference between two solutions New drawbacks: wave amplitude is always slightly exceeded water has effect of low-pass filter DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 26 / 53

30 Outline 1 Introduction 2 Tsunami generation 3 Dislocation dynamics 4 Sediments effect DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 27 / 53

31 Towards more realistic dynamic source model Fault dynamics modelling Earth crust is a linear viscoelastic material (Kelvin-Voigt model) Isotropic homogeneous or heterogeneous medium Fault modeled as a Volterra dislocation Displacement field is increased by the amount of the Burgers vector along any loop enclosing the dislocation d u = b C Simplified situation with respect to fracture mechanics: location and displacement jump are known Reference: D. Dutykh, F. Dias. Tsunami generation by dynamic displacement of sea bed due to dip-slip faulting. Accepted to Mathematics and computers in simulation, 2008 DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 28 / 53

32 Haskell s model (1969) Rupture propagation and rise time b( x, t) = 0 t ζ/v < 0 ( b 0 /T)(t ζ/v) 0 t ζ/v T b0 t ζ/v > T z x T is the rise time, V the rupture velocity O y ζ is a coordinate along the fault Front propagates unilaterally along the y axis W L 2 d L 2 δ DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 29 / 53

33 Coupled computation Seismology/hydrodynamics coupling Viscoelastodynamics Space derivatives are discretized by P2 FEM Second-order implicit time stepping Hydrodynamics Governing equations are NSWE: Solved by FVCF scheme η t + ((h+η) v ) = t h, v t v 2 + g η = 0. Coupling with FEM computation is done through the bathymetry h = h( x, t) DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 30 / 53

34 Comparison of two approaches Remarks about simulation Active generation: We simulate only first 10 s of the earthquake and couple it with hydrodynamic solver For t > 10 s we assume that bottom remains in its latest configuration Passive generation: Translate static dislocation solution onto free surface as initial condition Multiscale nature: Two different scales: elastic waves and water gravity waves DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 31 / 53

35 Results of numerical computation Comparison between passive and active generation DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 32 / 53

36 Discrepancy with tide gauges records Chilean 1960 event Reference: J.C. Borrero, B. Uslu, V. Titov, C.E. Synolakis (2006). Modeling tsunamis for California ports and harbors. Proceedings of the thirtieth International Conference on Coastal Engineering, ASCE DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 33 / 53

37 Outline 1 Introduction 2 Tsunami generation 3 Dislocation dynamics 4 Sediments effect DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 34 / 53

38 Sediments Definition and formation Definition: Sediment is any particulate matter that can be transported by fluid flow and which eventually is deposited as a layer of solid particles on the bottom of water. Formation mechanisms: Material washed into the sea from the land Remains of tiny forms of sea life Layers of sand settle on the ocean floor DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 35 / 53

39 DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 36 / 53

40 Global Seismic Hazard Map Figure: The map shows the peak ground acceleration (PGA) that a site can expect during the next 50 years with 10 percent probability DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 37 / 53

41 Previous studies: plate movement modelling Ref: C.W. Fuller, S.D. Willett and M.T. Brandon. Geology (2006) It appears that the most severe subduction zone earthquakes occur in areas where sediment-filled basins are found Method: Thermo-mechanical modelling of plate movement Frictional plastic rheology and thermally-activated viscous rheology Numerics: 2D FEM DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 38 / 53

42 Previous studies: simulation results Ref: Fuller et al. Geology (2006) Figure: No sedimentation Figure: Sedimentation DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 39 / 53

43 Previous studies: simulation results Ref: Fuller et al. Geology (2006) Figure: No sedimentation Figure: Sedimentation DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 39 / 53

44 Previous studies: simulation results Ref: Fuller et al. Geology (2006) Figure: No sedimentation Figure: Sedimentation DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 39 / 53

45 Previous studies: conclusions Reference: C.W. Fuller, S.D. Willett and M.T. Brandon. Formation of forearc basins and their influence on subduction zone earthquakes. Geology, 34, no. 2, pp (2006) Conclusion: Our numerical simulations demonstrate that sedimentation stabilizes the underlying wedge, preventing internal deformation beneath the basin [... ] The lack of deformation in these stable regions increases the likelihood of thermal pressurization of the subduction thrust, allows the fault to load faster, and allows greater healing of the fault between rupture events. DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 40 / 53

46 Numerical scheme Solving elastodynamic governing equations Variational form of governing equations: Ω ρ 2 u 2 t v dω+ Ω σ( u) : v dω = Ω ρ g v dω+ f v ds Γ N Numerics: Unstructured triangular mesh (can be adapted) Discretization in space by a P2 FEM Fully implicit discretization scheme in time Programming: FreeFem++ : program takes 200 lines with comments! DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 41 / 53

47 Two test-cases description Parameters used in computations and geometry of the domain. parameter d δ L b G g G s value 4000 m m 10 m Pa Pa ν g 0.27 ν s 0.3 (c s ) g (c s ) s 3600 m s 330 m s Table: From Mei and Foda (1981) Figure: Test-case with homogeneous medium y O d L x DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 42 / 53

48 Two test-cases description Parameters used in computations and geometry of the domain. parameter d δ L b G g G s value 4000 m m 10 m Pa Pa ν g 0.27 ν s 0.3 (c s ) g (c s ) s 3600 m s 330 m s Table: From Mei and Foda (1981) h s y O d Figure: Test-case with sediment layer L x DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 42 / 53

49 Static solution Comparison with/without sediment layer computations Static sea bed deformation Homogeneous half space With sand layer y, m x, km Figure: Volterra dislocation source: h s = 600 m DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 43 / 53

50 Dynamic computation - I Sand layer thickness h s = 600 m DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 44 / 53

51 Dynamic computation - II Sand layer thickness h s = 100 m DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 45 / 53

52 Sediment amplification factor - I How to quantify the influence of the sediment layer on maximal sea-bed amplitude? Definition: S a = max v s(x, t) (x,t) max v 0(x, t) 1 (x,t) v 0 (x, t): Vertical displacement at the sea-bed in homogeneous case v s (x, t): The same quantity in presence of the sediment layer Question: How does it depend on sediment layer relative thickness h s /d? DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 46 / 53

53 Sediment amplification factor - II Dependence on sediment layer relative thickness 0.8 Sediment Impact factor Author: Denys Dutykh, CMLA, ENS de Cachan S impact h /d s DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 47 / 53

54 Application to the Tsunami Boxing Day 2004 Seismic parameters and sediments thickness estimation Figure: Sediment thickness map DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 48 / 53

55 Application to the Tsunami Boxing Day 2004 Seismic parameters and sediments thickness estimation Figure: Seismic parameters DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 48 / 53

56 Application to the Tsunami Boxing Day 2004 Seismic parameters and sediments thickness estimation 0.8 Sediment Impact factor Author: Denys Dutykh, CMLA, ENS de Cachan S impact h s /d Figure: Sediment impact factor DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 48 / 53

57 New Scientist Phone conversation with Stephen Battersby... D. Dutykh, F. Dias. Influence of sedimentary layering on tsunami generation. Computer Methods in Applied Mechanics and Engineering. In Press DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 49 / 53

58 Further study of sediments effects In collaboration with E. Okal (Northwestern University) Normal modes framework : Numerous cases were considered with real-world parameters Sediments effect is computed in the near field and in the far field Some examples are shown below Conclusion: Results sometimes show amplification, but not in the most conditions! DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 50 / 53

59 Conclusions Conclusions: Sediment layer has no effect on static deformation, but... It makes big difference on dynamics: Slows down elastic and Rayleigh wave propagation Has strong effect on the amplitude: amplification factor up to 70% We should (maybe) revise initial condition for tsunami wave modelling in some subduction zones according to Sediment Thickness Map No information about sediments in Mediterranean region Perspectives: Perform these computations in 3D Take into account porosity in the mud layer DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 51 / 53

60 Other research directions Tsunami wave modelling (with F. Dias and R. Poncet) Tsunami energy estimation Visco-potential flows (with F. Dias) Spectral methods for water wave problem (with D. Clamond, LJAD) Two-phase flows (with LAMA, CMLA, UCD) Powder snow avalanches (in collaboration with D. Bresch & C. Acary-Robert) Modelling (Kazhikov-Smagulov, Brenner-NS,... ) Numerical simulations Protecting structures efficiency estimation DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 52 / 53

61 Thank you for your attention! dutykh DENYS DUTYKH (CNRS LAMA) Tsunami wave generation Institut d Alembert 53 / 53

Numerical simulation of dispersive waves

Numerical simulation of dispersive waves Numerical simulation of dispersive waves DENYS DUTYKH 1 Chargé de Recherche CNRS 1 LAMA, Université de Savoie 73376 Le Bourget-du-Lac, France Colloque EDP-Normandie DENYS DUTYKH (CNRS LAMA) Dispersive

More information

arxiv: v5 [physics.class-ph] 12 Jul 2012

arxiv: v5 [physics.class-ph] 12 Jul 2012 ON THE CONTRIBUTION OF THE HORIZONTAL SEA-BED DISPLACEMENTS INTO THE TSUNAMI GENERATION PROCESS DENYS DUTYKH, DIMITRIOS MITSOTAKIS, LEONID B. CHUBAROV, AND YURIY I. SHOKIN arxiv:111.1741v5 [physics.class-ph]

More information

Resonant wave run-up on sloping beaches and vertical walls

Resonant wave run-up on sloping beaches and vertical walls Resonant wave run-up on sloping beaches and vertical walls DENYS DUTYKH 1 Chargé de Recherche CNRS 1 Université de Savoie Mont Blanc Laboratoire de Mathématiques (LAMA) 73376 Le Bourget-du-Lac France Seminar:

More information

Three Dimensional Simulations of Tsunami Generation and Propagation

Three Dimensional Simulations of Tsunami Generation and Propagation Chapter 1 Earth Science Three Dimensional Simulations of Tsunami Generation and Propagation Project Representative Takashi Furumura Authors Tatsuhiko Saito Takashi Furumura Earthquake Research Institute,

More information

Energy of tsunami waves generated by bottom motion arxiv: v2 [physics.ao-ph] 23 Oct 2008

Energy of tsunami waves generated by bottom motion arxiv: v2 [physics.ao-ph] 23 Oct 2008 Energy of tsunami waves generated by bottom motion arxiv:0808.2861v2 [physics.ao-ph] 23 Oct 2008 By Denys Dutykh, Frédéric Dias CMLA, ENS Cachan, CNRS, PRES UniverSud, 61, avenue du President Wilson, 94230

More information

RELATION BETWEEN RAYLEIGH WAVES AND UPLIFT OF THE SEABED DUE TO SEISMIC FAULTING

RELATION 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 information

Earthquake and Volcano Deformation

Earthquake and Volcano Deformation Earthquake and Volcano Deformation Paul Segall Stanford University Draft Copy September, 2005 Last Updated Sept, 2008 COPYRIGHT NOTICE: To be published by Princeton University Press and copyrighted, c

More information

Algebraic geometry for shallow capillary-gravity waves

Algebraic geometry for shallow capillary-gravity waves Algebraic geometry for shallow capillary-gravity waves DENYS DUTYKH 1 Chargé de Recherche CNRS 1 Université de Savoie Laboratoire de Mathématiques (LAMA) 73376 Le Bourget-du-Lac France Seminar of Computational

More information

Kinematics of the Southern California Fault System Constrained by GPS Measurements

Kinematics of the Southern California Fault System Constrained by GPS Measurements Title Page Kinematics of the Southern California Fault System Constrained by GPS Measurements Brendan Meade and Bradford Hager Three basic questions Large historical earthquakes One basic question How

More information

Seth Stein and Emile Okal, Department of Geological Sciences, Northwestern University, Evanston IL USA. Revised 2/5/05

Seth Stein and Emile Okal, Department of Geological Sciences, Northwestern University, Evanston IL USA. Revised 2/5/05 Sumatra earthquake moment from normal modes 2/6/05 1 Ultra-long period seismic moment of the great December 26, 2004 Sumatra earthquake and implications for the slip process Seth Stein and Emile Okal,

More information

Earthquakes and Earthquake Hazards Earth - Chapter 11 Stan Hatfield Southwestern Illinois College

Earthquakes and Earthquake Hazards Earth - Chapter 11 Stan Hatfield Southwestern Illinois College Earthquakes and Earthquake Hazards Earth - Chapter 11 Stan Hatfield Southwestern Illinois College What Is an Earthquake? An earthquake is the vibration of Earth, produced by the rapid release of energy.

More information

Influence of a sedimentary basin infilling description on the 2D P-SV wave propagation using linear and nonlinear constitutive models

Influence of a sedimentary basin infilling description on the 2D P-SV wave propagation using linear and nonlinear constitutive models Influence of a sedimentary basin infilling description on the 2D P-SV wave propagation using linear and nonlinear constitutive models C. Gélis IRSN, France L.F. Bonilla Université Paris Est - IFSTTAR,

More information

LECTURE 6 EARTHQUAKES AS TSUNAMI SOURCES

LECTURE 6 EARTHQUAKES AS TSUNAMI SOURCES LECTURE 6 EARTHQUAKES AS TSUNAMI SOURCES Northwestern University, 2007 TSUNAMI GENERATION by EARTHQUAKE SOURCES CLASSICAL APPROA CH 1. Consider model of EarthquakeRupture 2. Compute Static Deformation

More information

Differentiating earthquake tsunamis from other sources; how do we tell the difference?

Differentiating earthquake tsunamis from other sources; how do we tell the difference? Differentiating earthquake tsunamis from other sources; how do we tell the difference? David Tappin (1), Stephan Grilli (2), Jeffrey Harris (2), Timothy Masterlark (3), James Kirby (4), Fengyan Shi Shi

More information

Surface changes caused by erosion and sedimentation were treated by solving: (2)

Surface changes caused by erosion and sedimentation were treated by solving: (2) GSA DATA REPOSITORY 214279 GUY SIMPSON Model with dynamic faulting and surface processes The model used for the simulations reported in Figures 1-3 of the main text is based on two dimensional (plane strain)

More information

Inversion of tsunami data. A. Sladen CNRS, Géoazur 1/35

Inversion of tsunami data. A. Sladen CNRS, Géoazur 1/35 Inversion of tsunami data A. Sladen CNRS, Géoazur 1/35 DEFINITION Tsunami waves are gravity wave with a long period need a BIG source! 2/35 DEFINITION Krakatoa, 1883 Summer 2015, E.T. pers. comm. Lituya

More information

Earthquakes. Earthquake Magnitudes 10/1/2013. Environmental Geology Chapter 8 Earthquakes and Related Phenomena

Earthquakes. Earthquake Magnitudes 10/1/2013. Environmental Geology Chapter 8 Earthquakes and Related Phenomena Environmental Geology Chapter 8 Earthquakes and Related Phenomena Fall 2013 Northridge 1994 Kobe 1995 Mexico City 1985 China 2008 Earthquakes Earthquake Magnitudes Earthquake Magnitudes Richter Magnitude

More information

Earthquakes. Building Earth s Surface, Part 2. Science 330 Summer What is an earthquake?

Earthquakes. Building Earth s Surface, Part 2. Science 330 Summer What is an earthquake? Earthquakes Building Earth s Surface, Part 2 Science 330 Summer 2005 What is an earthquake? An earthquake is the vibration of Earth produced by the rapid release of energy Energy released radiates in all

More information

Research Article Linearized Shallow-water Wave Theory of Tsunami Generation and Propagation by Three-dimensional Stochastic Seismic Bottom Topography

Research Article Linearized Shallow-water Wave Theory of Tsunami Generation and Propagation by Three-dimensional Stochastic Seismic Bottom Topography Research Journal of Applied Sciences, Engineering and Technology 7(19): 4035-4055, 2014 DOI:10.19026/rjaset.7.765 ISSN: 2040-7459; e-issn: 2040-7467 2014 Maxwell Scientific Publication Corp. Submitted:

More information

Modified Serre Green Naghdi equations with improved or without dispersion

Modified Serre Green Naghdi equations with improved or without dispersion Modified Serre Green Naghdi equations with improved or without dispersion DIDIER CLAMOND Université Côte d Azur Laboratoire J. A. Dieudonné Parc Valrose, 06108 Nice cedex 2, France didier.clamond@gmail.com

More information

Materials and Methods The deformation within the process zone of a propagating fault can be modeled using an elastic approximation.

Materials and Methods The deformation within the process zone of a propagating fault can be modeled using an elastic approximation. Materials and Methods The deformation within the process zone of a propagating fault can be modeled using an elastic approximation. In the process zone, stress amplitudes are poorly determined and much

More information

UGRC 144 Science and Technology in Our Lives/Geohazards

UGRC 144 Science and Technology in Our Lives/Geohazards UGRC 144 Science and Technology in Our Lives/Geohazards Session 3 Understanding Earthquakes and Earthquake Hazards Lecturer: Dr. Patrick Asamoah Sakyi Department of Earth Science, UG Contact Information:

More information

I. Locations of Earthquakes. Announcements. Earthquakes Ch. 5. video Northridge, California earthquake, lecture on Chapter 5 Earthquakes!

I. Locations of Earthquakes. Announcements. Earthquakes Ch. 5. video Northridge, California earthquake, lecture on Chapter 5 Earthquakes! 51-100-21 Environmental Geology Summer 2006 Tuesday & Thursday 6-9:20 p.m. Dr. Beyer Earthquakes Ch. 5 I. Locations of Earthquakes II. Earthquake Processes III. Effects of Earthquakes IV. Earthquake Risk

More information

Co-seismic Gravity Changes Computed for a Spherical Earth Model Applicable to GRACE Data

Co-seismic Gravity Changes Computed for a Spherical Earth Model Applicable to GRACE Data Chapter 2 Co-seismic Gravity Changes Computed for a Spherical Earth Model Applicable to GRACE Data W.Sun,G.Fu,andSh.Okubo Abstract Dislocation theories were developed conventionally for a deformed earth

More information

Earthquakes Chapter 19

Earthquakes Chapter 19 Earthquakes Chapter 19 Does not contain complete lecture notes. What is an earthquake An earthquake is the vibration of Earth produced by the rapid release of energy Energy released radiates in all directions

More information

2C09 Design for seismic and climate changes

2C09 Design for seismic and climate changes 2C09 Design for seismic and climate changes Lecture 10: Characterisation of seismic motion Aurel Stratan, Politehnica University of Timisoara 07/04/2017 European Erasmus Mundus Master Course Sustainable

More information

Hitoshi Hirose (1), and Kazuro Hirahara (2) Abstract. Introduction

Hitoshi Hirose (1), and Kazuro Hirahara (2) Abstract. Introduction Three dimensional simulation for the earthquake cycle at a subduction zone based on a rate- and state-dependent friction law: Insight into a finiteness and a variety of dip-slip earthquakes Hitoshi Hirose

More information

Contribution of HPC to the mitigation of natural risks. B. Feignier. CEA-DAM Ile de France Département Analyse, Surveillance, Environnement

Contribution of HPC to the mitigation of natural risks. B. Feignier. CEA-DAM Ile de France Département Analyse, Surveillance, Environnement Contribution of HPC to the mitigation of natural risks B. Feignier CEA-DAM Ile de France Département Analyse, Surveillance, Environnement Introduction Over the last 40 years, the increase in computational

More information

Supporting Information for Inferring field-scale properties of a fractured aquifer from ground surface deformation during a well test

Supporting Information for Inferring field-scale properties of a fractured aquifer from ground surface deformation during a well test GEOPHYSICAL RESEARCH LETTERS Supporting Information for Inferring field-scale properties of a fractured aquifer from ground surface deformation during a well test Jonathan Schuite 1, Laurent Longuevergne

More information

Preliminary Study of Possible Tsunami Hazards in Taiwan Region

Preliminary Study of Possible Tsunami Hazards in Taiwan Region Preliminary Study of Possible Tsunami Hazards in Taiwan Region Xiaoming Wang and Philip L.-F. Liu Cornell University (First Draft on May 25 2006) (Second Draft on June 1 2006) (Final Update on June 8 2006)

More information

Elastoplastic Deformation in a Wedge-Shaped Plate Caused By a Subducting Seamount

Elastoplastic Deformation in a Wedge-Shaped Plate Caused By a Subducting Seamount Elastoplastic Deformation in a Wedge-Shaped Plate Caused By a Subducting Seamount Min Ding* 1, and Jian Lin 2 1 MIT/WHOI Joint Program, 2 Woods Hole Oceanographic Institution *Woods Hole Oceanographic

More information

Name Date Class. radiate in all directions, carrying some of the. of plate boundaries have different usual patterns of.

Name Date Class. radiate in all directions, carrying some of the. of plate boundaries have different usual patterns of. Chapter Outline Earthquakes CHAPTER 6 Lesson 1: Earthquakes and Plate Boundaries A. What is an earthquake? 1. A(n) is the rupture and sudden movement of rocks along a fault. A fault is a fracture surface

More information

ON NEAR-FIELD GROUND MOTIONS OF NORMAL AND REVERSE FAULTS FROM VIEWPOINT OF DYNAMIC RUPTURE MODEL

ON 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 information

Very basic tsunami physics...

Very basic tsunami physics... Very basic tsunami physics... Energy Bottom uplift & Waterberg formation E R 4.8 + 1.5M E T = 1 2 ρglλ(δh)2 L ~ 1 6 m; λ ~ 1 4 m; δh ~ 5m E R 1 18 J 1 2 E T Center of mass falls... Wavelength λ H ~ 4;

More information

Ground displacement in a fault zone in the presence of asperities

Ground displacement in a fault zone in the presence of asperities BOLLETTINO DI GEOFISICA TEORICA ED APPLICATA VOL. 40, N. 2, pp. 95-110; JUNE 2000 Ground displacement in a fault zone in the presence of asperities S. SANTINI (1),A.PIOMBO (2) and M. DRAGONI (2) (1) Istituto

More information

TSUNAMI HAZARD ASSESSMENT IN NORTHERN EGYPT USING NUMERICAL SIMULATION

TSUNAMI HAZARD ASSESSMENT IN NORTHERN EGYPT USING NUMERICAL SIMULATION TSUNAMI HAZARD ASSESSMENT IN NORTHERN EGYPT USING NUMERICAL SIMULATION Abutaleb Ali Supervisor: Bunichiro SHIBAZAKI MEE16717 Yushiro FUJII ABSTRACT To investigate the tsunami hazard along the northern

More information

7 SEISMIC LOADS. 7.1 Estimation of Seismic Loads. 7.2 Calculation of Seismic Loads

7 SEISMIC LOADS. 7.1 Estimation of Seismic Loads. 7.2 Calculation of Seismic Loads 1 7 SEISMIC LOADS 7.1 Estimation of Seismic Loads 7.1.1 Seismic load and design earthquake motion (1) For ordinary buildings, seismic load is evaluated using the acceleration response spectrum (see Sec.7.2)

More information

Synthetic Seismicity Models of Multiple Interacting Faults

Synthetic Seismicity Models of Multiple Interacting Faults Synthetic Seismicity Models of Multiple Interacting Faults Russell Robinson and Rafael Benites Institute of Geological & Nuclear Sciences, Box 30368, Lower Hutt, New Zealand (email: r.robinson@gns.cri.nz).

More information

Earthquakes Earth, 9th edition, Chapter 11 Key Concepts What is an earthquake? Earthquake focus and epicenter What is an earthquake?

Earthquakes Earth, 9th edition, Chapter 11 Key Concepts What is an earthquake? Earthquake focus and epicenter What is an earthquake? 1 2 3 4 5 6 7 8 9 10 Earthquakes Earth, 9 th edition, Chapter 11 Key Concepts Earthquake basics. "" and locating earthquakes.. Destruction resulting from earthquakes. Predicting earthquakes. Earthquakes

More information

Slip 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 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 information

Earthquake and Volcano Clustering at Mono Basin (California)

Earthquake and Volcano Clustering at Mono Basin (California) Excerpt from the Proceedings of the COMSOL Conference 2010 Paris Earthquake and Volcano Clustering at Mono Basin (California) D. La Marra *,1, A. Manconi 2,3 and M. Battaglia 1 1 Dept of Earth Sciences,

More information

Simulation of earthquake rupture process and strong ground motion

Simulation 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 information

Gravity Tectonics Volcanism Atmosphere Water Winds Chemistry. Planetary Surfaces

Gravity Tectonics Volcanism Atmosphere Water Winds Chemistry. Planetary Surfaces Gravity Tectonics Volcanism Atmosphere Water Winds Chemistry Planetary Surfaces Gravity & Rotation Polar flattening caused by rotation is the largest deviation from a sphere for a planet sized object (as

More information

Elastic Rebound Theory

Elastic Rebound Theory Earthquakes Elastic Rebound Theory Earthquakes occur when strain exceeds the strength of the rock and the rock fractures. The arrival of earthquakes waves is recorded by a seismograph. The amplitude of

More information

Derivation of Table 2 in Okada (1992)

Derivation of Table 2 in Okada (1992) Derivation of Table 2 in Okada (1992) [ I ] Derivation of Eqs.(4) through (6) Eq.(1) of Okada (1992) can be rewritten, where is the displacement at due to a j-th direction single force of unit magnitude

More information

The Mechanics of Earthquakes and Faulting

The Mechanics of Earthquakes and Faulting The Mechanics of Earthquakes and Faulting Christopher H. Scholz Lamont-Doherty Geological Observatory and Department of Earth and Environmental Sciences, Columbia University 2nd edition CAMBRIDGE UNIVERSITY

More information

1. classic definition = study of deformed rocks in the upper crust

1. classic definition = study of deformed rocks in the upper crust Structural Geology I. Introduction 1. classic definition = study of deformed rocks in the upper crust deformed includes translation, rotation, and strain (change of shape) All rocks are deformed in some

More information

Tsunami waveform analyses of the 2006 underthrust and 2007 outer-rise Kurile earthquakes

Tsunami waveform analyses of the 2006 underthrust and 2007 outer-rise Kurile earthquakes Author(s) 2008. This work is licensed under a Creative Commons License. Advances in Geosciences Tsunami waveform analyses of the 2006 underthrust and 2007 outer-rise Kurile earthquakes Y. Tanioka 1, Y.

More information

Scientific Research on the Cascadia Subduction Zone that Will Help Improve Seismic Hazard Maps, Building Codes, and Other Risk-Mitigation Measures

Scientific Research on the Cascadia Subduction Zone that Will Help Improve Seismic Hazard Maps, Building Codes, and Other Risk-Mitigation Measures Scientific Research on the Cascadia Subduction Zone that Will Help Improve Seismic Hazard Maps, Building Codes, and Other Risk-Mitigation Measures Art Frankel U.S. Geological Survey Seattle, WA GeoPrisms-Earthscope

More information

Defmod, an earthquake simulator that adaptively switches between quasi-static and dynamic states

Defmod, an earthquake simulator that adaptively switches between quasi-static and dynamic states Defmod, an earthquake simulator that adaptively switches between quasi-static and dynamic states C. Meng (cmeng@mit.edu) B. Hager ERL, MIT May 16, 2016 Reservoir production/injection induced seismicity

More information

Rotation of the Principal Stress Directions Due to Earthquake Faulting and Its Seismological Implications

Rotation of the Principal Stress Directions Due to Earthquake Faulting and Its Seismological Implications Bulletin of the Seismological Society of America, Vol. 85, No. 5, pp. 1513-1517, October 1995 Rotation of the Principal Stress Directions Due to Earthquake Faulting and Its Seismological Implications by

More information

Important Concepts. Earthquake hazards can be categorized as:

Important Concepts. Earthquake hazards can be categorized as: Lecture 1 Page 1 Important Concepts Monday, August 17, 2009 1:05 PM Earthquake Engineering is a branch of Civil Engineering that requires expertise in geology, seismology, civil engineering and risk assessment.

More information

log 4 0.7m log m Seismic Analysis of Structures by TK Dutta, Civil Department, IIT Delhi, New Delhi. Module 1 Seismology Exercise Problems :

log 4 0.7m log m Seismic Analysis of Structures by TK Dutta, Civil Department, IIT Delhi, New Delhi. Module 1 Seismology Exercise Problems : Seismic Analysis of Structures by TK Dutta, Civil Department, IIT Delhi, New Delhi. Module Seismology Exercise Problems :.4. Estimate the probabilities of surface rupture length, rupture area and maximum

More information

High-Resolution Finite Volume Methods and Adaptive Mesh Refinement

High-Resolution Finite Volume Methods and Adaptive Mesh Refinement High-Resolution Finite Volume Methods and Adaptive Mesh Refinement Randall J. LeVeque Department of Applied Mathematics University of Washington CLAWPACK and TsunamiClaw Software http://www.amath.washington.edu/~claw

More information

TSUNAMI PROPAGATION AND INUNDATION MODELINGS ALONG SOUTH-EAST COAST OF PAPUA NEW GUINEA

TSUNAMI PROPAGATION AND INUNDATION MODELINGS ALONG SOUTH-EAST COAST OF PAPUA NEW GUINEA TSUNAMI PROPAGATION AND INUNDATION MODELINGS ALONG SOUTH-EAST COAST OF PAPUA NEW GUINEA Martin WAREK Supervisor: Yushiro FUJII MEE12620 Bunichiro SHIBAZAKI ABSTRACT This study covers tsunami generation,

More information

Dynamic Earthquake Triggering Due to Stress from Surface Wave Particle Displacement

Dynamic Earthquake Triggering Due to Stress from Surface Wave Particle Displacement Emily Morton Geop 523: Theoretical Seismology 29 April 2011 Introduction Dynamic Earthquake Triggering Due to Stress from Surface Wave Particle Displacement Earthquakes can be triggered by other earthquakes.

More information

Research Article Theoretical Study on the Flow Generated by the Strike-Slip Faulting

Research Article Theoretical Study on the Flow Generated by the Strike-Slip Faulting ISRN Geophysics Volume 23, Article ID 947672, 7 pages http://dx.doi.org/.55/23/947672 Research Article heoretical Study on the Flow Generated by the Strike-Slip Faulting Chi-Min Liu,,2 Ray-eng ang, 2,3

More information

Salmon: Introduction to ocean waves

Salmon: Introduction to ocean waves 9 The shallow-water equations. Tsunamis. Our study of waves approaching the beach had stopped at the point of wave breaking. At the point of wave breaking, the linear theory underlying Propositions #1

More information

RELAXED VARIATIONAL PRINCIPLE FOR WATER WAVE MODELING

RELAXED VARIATIONAL PRINCIPLE FOR WATER WAVE MODELING RELAXED VARIATIONAL PRINCIPLE FOR WATER WAVE MODELING DENYS DUTYKH 1 Senior Research Fellow UCD & Chargé de Recherche CNRS 1 University College Dublin School of Mathematical Sciences Workshop on Ocean

More information

News Release December 30, 2004 The Science behind the Aceh Earthquake

News Release December 30, 2004 The Science behind the Aceh Earthquake News Release December 30, 2004 The Science behind the Aceh Earthquake PASADENA, Calif. - Kerry Sieh, the Robert P. Sharp Professor of Geology at the California Institute of Technology and a member of Caltech's

More information

STUDY ON APPROPRIATE MODELING OF TSUNAMIS IN MALAYSIA FOR RISK EVALUATION

STUDY ON APPROPRIATE MODELING OF TSUNAMIS IN MALAYSIA FOR RISK EVALUATION STUDY ON APPROPRIATE MODELING OF TSUNAMIS IN MALAYSIA FOR RISK EVALUATION Zaty Aktar binti Mokhtar* Supervisor: Fumihiko Imamura** MEE06025 Shunichi Koshimura** ABSTRACT In order to design a tsunami warning

More information

Effect of the Emperor seamounts on trans-oceanic propagation of the 2006 Kuril Island earthquake tsunami

Effect 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 information

Spectral Element simulation of rupture dynamics

Spectral 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 information

Widespread 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 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 information

Rheology III. Ideal materials Laboratory tests Power-law creep The strength of the lithosphere The role of micromechanical defects in power-law creep

Rheology III. Ideal materials Laboratory tests Power-law creep The strength of the lithosphere The role of micromechanical defects in power-law creep Rheology III Ideal materials Laboratory tests Power-law creep The strength of the lithosphere The role of micromechanical defects in power-law creep Ideal materials fall into one of the following categories:

More information

Deformation cycles of great subduction earthquakes in a viscoelastic Earth

Deformation cycles of great subduction earthquakes in a viscoelastic Earth Deformation cycles of great subduction earthquakes in a viscoelastic Earth Kelin Wang Pacific Geoscience Centre, Geological Survey of Canada School of Earth and Ocean Science, University of Victoria????

More information

Instabilities and Dynamic Rupture in a Frictional Interface

Instabilities and Dynamic Rupture in a Frictional Interface Instabilities and Dynamic Rupture in a Frictional Interface Laurent BAILLET LGIT (Laboratoire de Géophysique Interne et Tectonophysique) Grenoble France laurent.baillet@ujf-grenoble.fr http://www-lgit.obs.ujf-grenoble.fr/users/lbaillet/

More information

GEOL 10: Environmental Geology Mid Term II Study Guide

GEOL 10: Environmental Geology Mid Term II Study Guide (1) 10 pts. Use the seismographs below, along with the travel time plot for California, and the location map for seismographs in this region to do the following: a. Determine the origin time for the earthquake

More information

2. Tsunami Source Details

2. Tsunami Source Details 2. Tsunami Source Details The Northland area faces a range of potential tsunamigenic sources that include several local and distant fault systems and underwater landslides. A NIWA study (Goff et al. 2006)

More information

Elastic rebound theory

Elastic rebound theory Elastic rebound theory Focus epicenter - wave propagation Dip-Slip Fault - Normal Normal Fault vertical motion due to tensional stress Hanging wall moves down, relative to the footwall Opal Mountain, Mojave

More information

NUMERICAL SIMULATION AS GUIDANCE IN MAKING TSUNAMI HAZARD MAP FOR LABUAN ISLAND

NUMERICAL SIMULATION AS GUIDANCE IN MAKING TSUNAMI HAZARD MAP FOR LABUAN ISLAND NUMERICAL SIMULATION AS GUIDANCE IN MAKING TSUNAMI HAZARD MAP FOR LABUAN ISLAND MOHD RIDZUAN bin Adam Supervisor: Fumihiko IMAMURA MEE09199 ABSTRACT At the northeast end of the South China Sea, tsunamis

More information

Seismic hazard analysis of Tianjin area based on strong ground motion prediction

Seismic hazard analysis of Tianjin area based on strong ground motion prediction Earthq Sci (2010)23: 369 375 369 Doi: 10.1007/s11589-010-0734-6 Seismic hazard analysis of Tianjin area based on strong ground motion prediction Zhao Boming School of Civil Engineering, Beijing Jiaotong

More information

Tsunami potential and modeling

Tsunami potential and modeling Tsunami potential and modeling GEORGE PRIEST OREGON DEPT. OF GEOLOGY AND MINERAL INDUSTRIES NEWPORT COASTAL FIELD OFFICE April 7, 2012 GeoPRISMS Cascadia Workshop, Portland, Oregon What creates most uncertainty

More information

Comparison between three-dimensional linear and nonlinear tsunami generation models

Comparison between three-dimensional linear and nonlinear tsunami generation models Please note that this is an author-produced PDF of an article accepted for publication following peer review. The definitive publisher-authenticated version is available on the publisher Web site Theoretical

More information

3D Finite Element Modeling of fault-slip triggering caused by porepressure

3D Finite Element Modeling of fault-slip triggering caused by porepressure 3D Finite Element Modeling of fault-slip triggering caused by porepressure changes Arsalan Sattari and David W. Eaton Department of Geoscience, University of Calgary Suary We present a 3D model using a

More information

Fault-Zone Properties arid Earthquake Rupture Dynamics

Fault-Zone Properties arid Earthquake Rupture Dynamics Fault-Zone Properties arid Earthquake Rupture Dynamics Eiichi Fukuyama National Research Institute for Earth Science,and Disaster Prevention Tsukuba, Japan i ':-> i' ' -':'.." \,' " '' L VH ELSEVIER AMSTERDAM

More information

Effects of Surface Geology on Seismic Motion

Effects 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 VELOCITY STRUCTURE INVERSIONS FROM HORIZONTAL TO VERTICAL

More information

Section Forces Within Earth. 8 th Grade Earth & Space Science - Class Notes

Section Forces Within Earth. 8 th Grade Earth & Space Science - Class Notes Section 19.1 - Forces Within Earth 8 th Grade Earth & Space Science - Class Notes Stress and Strain Stress - is the total force acting on crustal rocks per unit of area (cause) Strain deformation of materials

More information

TSUNAMI HAZARD ASSESSMENT FOR THE CENTRAL COAST OF PERU USING NUMERICAL SIMULATIONS FOR THE 1974, 1966 AND 1746 EARTHQUAKES

TSUNAMI HAZARD ASSESSMENT FOR THE CENTRAL COAST OF PERU USING NUMERICAL SIMULATIONS FOR THE 1974, 1966 AND 1746 EARTHQUAKES TSUNAMI HAZARD ASSESSMENT FOR THE CENTRAL COAST OF PERU USING NUMERICAL SIMULATIONS FOR THE 1974, 1966 AND 1746 EARTHQUAKES Sheila Yauri Supervisor: Yushiro FUJII MEE10521 Bunichiro SHIBAZAKI ABSTRACT

More information

Tsunamis and ocean waves

Tsunamis and ocean waves Department of Mathematics & Statistics AAAS Annual Meeting St. Louis Missouri February 19, 2006 Introduction Tsunami waves are generated relatively often, from various sources Serious tsunamis (serious

More information

Source 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 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 information

Earthquakes. Pt Reyes Station 1906

Earthquakes. Pt Reyes Station 1906 Earthquakes Pt Reyes Station 1906 Earthquakes Ground shaking caused by the sudden release of accumulated strain by an abrupt shift of rock along a fracture in the earth. You Live in Earthquake Country

More information

Linear potential theory for tsunami generation and propagation

Linear potential theory for tsunami generation and propagation Linear potential theory for tsunami generation and propagation 2 3 Tatsuhiko Saito 4 5 6 7 8 9 0 2 3 4 5 National Research Institute for Earth Science and Disaster Prevention, Tsukuba, Ibaraki, Japan Running

More information

UNIT - 7 EARTHQUAKES

UNIT - 7 EARTHQUAKES UNIT - 7 EARTHQUAKES WHAT IS AN EARTHQUAKE An earthquake is a sudden motion or trembling of the Earth caused by the abrupt release of energy that is stored in rocks. Modern geologists know that most earthquakes

More information

Effect of an outer-rise earthquake on seismic cycle of large interplate earthquakes estimated from an instability model based on friction mechanics

Effect of an outer-rise earthquake on seismic cycle of large interplate earthquakes estimated from an instability model based on friction mechanics Effect of an outer-rise earthquake on seismic cycle of large interplate earthquakes estimated from an instability model based on friction mechanics Naoyuki Kato (1) and Tomowo Hirasawa (2) (1) Geological

More information

Stress equilibrium in southern California from Maxwell stress function models fit to both earthquake data and a quasi-static dynamic simulation

Stress equilibrium in southern California from Maxwell stress function models fit to both earthquake data and a quasi-static dynamic simulation Stress equilibrium in southern California from Maxwell stress function models fit to both earthquake data and a quasi-static dynamic simulation Peter Bird Dept. of Earth, Planetary, and Space Sciences

More information

What Is an Earthquake? What Is an Earthquake? Earthquake

What Is an Earthquake? What Is an Earthquake? Earthquake Earthquakes Earth,, Chapter 11 Built Earthquakes: on shaky summary ground, in haiku form Memphis, south of New Madrid - whole lotta shakin'... What Is an Earthquake? An earthquake is ground shaking caused

More information

FRICTIONAL HEATING DURING AN EARTHQUAKE. Kyle Withers Qian Yao

FRICTIONAL HEATING DURING AN EARTHQUAKE. Kyle Withers Qian Yao FRICTIONAL HEATING DURING AN EARTHQUAKE Kyle Withers Qian Yao Temperature Change Along Fault Mode II (plain strain) crack rupturing bilaterally at a constant speed v r Idealize earthquake ruptures as shear

More information

Detailed analysis of tsunami waveforms generated by the 1946 Aleutian tsunami earthquake

Detailed analysis of tsunami waveforms generated by the 1946 Aleutian tsunami earthquake Detailed analysis of tsunami waveforms generated by the 1946 Aleutian tsunami earthquake Y. Tanioka, T. Seno To cite this version: Y. Tanioka, T. Seno. Detailed analysis of tsunami waveforms generated

More information

Mechanics of Earthquakes and Faulting

Mechanics of Earthquakes and Faulting Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Overview Milestones in continuum mechanics Concepts of modulus and stiffness. Stress-strain relations Elasticity Surface and body

More information

Data Repository Hampel et al., page 1/5

Data Repository Hampel et al., page 1/5 GSA DATA REPOSITORY 2138 Data Repositor Hampel et al., page 1/5 SETUP OF THE FINITE-ELEMENT MODEL The finite-element models were created with the software ABAQUS and consist of a 1-km-thick lithosphere,

More information

Geophysics Departmental Exam: 2004 Part 1

Geophysics Departmental Exam: 2004 Part 1 2004 Geophysics Departmental Exam: 2004 Part 1 This section is 90 minutes, closed book, and consists of questions designed to test your knowledge of facts and figures in the geosciences. The focus will

More information

Earthquakes. Forces Within Eartth. Faults form when the forces acting on rock exceed the rock s strength.

Earthquakes. Forces Within Eartth. Faults form when the forces acting on rock exceed the rock s strength. Earthquakes Vocabulary: Stress Strain Elastic Deformation Plastic Deformation Fault Seismic Wave Primary Wave Secondary Wave Focus Epicenter Define stress and strain as they apply to rocks. Distinguish

More information

Originally published as:

Originally published as: Originally published as: Lorenzo Martín, F., Wang, R., Roth, F. (2002): The effect of input parameters on visco-elastic models of crustal deformation. - Física de la Tierra, 14, 33-54 The effect of input

More information

TSUNAMI CHARACTERISTICS OF OUTER-RISE EARTHQUAKES ALONG THE PACIFIC COAST OF NICARAGUA - A CASE STUDY FOR THE 2016 NICARAGUA EVENT-

TSUNAMI CHARACTERISTICS OF OUTER-RISE EARTHQUAKES ALONG THE PACIFIC COAST OF NICARAGUA - A CASE STUDY FOR THE 2016 NICARAGUA EVENT- TSUNAMI CHARACTERISTICS OF OUTER-RISE EARTHQUAKES ALONG THE PACIFIC COAST OF NICARAGUA - A CASE STUDY FOR THE 2016 NICARAGUA EVENT- Amilcar Cabrera Supervisor: Yuichiro TANIOKA MEE16718 ABSTRACT Nicaragua

More information

Numerical Modeling for the Propagation of Tsunami Wave and Corresponding Inundation

Numerical Modeling for the Propagation of Tsunami Wave and Corresponding Inundation IOSR Journal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 78-1684,p-ISSN: 30-334X, Volume 1, Issue Ver. IV (Mar - Apr. 015), PP 55-6 www.iosrjournals.org Numerical Modeling for the Propagation

More information

A METHOD FOR DETERMINING ASPERITY PARAMETERS PRODUCING SPECIFIC MAXIMUM GROUND MOTION

A METHOD FOR DETERMINING ASPERITY PARAMETERS PRODUCING SPECIFIC MAXIMUM GROUND MOTION 13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 24 Paper No. 395 A METHOD FOR DETERMINING ASPERITY PARAMETERS PRODUCING SPECIFIC MAXIMUM GROUND MOTION Masayuki YOSHIMI

More information

An Earthquake is a rapid vibration or shaking of the Earth s crust created by a release in energy from sudden movement of a part of a plate along a

An Earthquake is a rapid vibration or shaking of the Earth s crust created by a release in energy from sudden movement of a part of a plate along a An Earthquake is a rapid vibration or shaking of the Earth s crust created by a release in energy from sudden movement of a part of a plate along a fault. Energy released radiates in all directions from

More information

Geo736: Seismicity along mid-ocean ridges

Geo736: Seismicity along mid-ocean ridges Geo736: Seismicity along mid-ocean ridges Course Notes: S. G. Wesnousky Spring 2018 Bathymetric maps show the ocean basins of the world are characteristically divided by a bathymetric ridge. The bathymetric

More information

Three Fs of earthquakes: forces, faults, and friction. Slow accumulation and rapid release of elastic energy.

Three Fs of earthquakes: forces, faults, and friction. Slow accumulation and rapid release of elastic energy. Earthquake Machine Stick-slip: Elastic Rebound Theory Jerky motions on faults produce EQs Three Fs of earthquakes: forces, faults, and friction. Slow accumulation and rapid release of elastic energy. Three

More information