Toward the next generation of earthquake source models by accounting for model prediction error

Download Presentation

Toward the next generation of earthquake source models by accounting for model prediction error

Loading in 2 Seconds...

- 62 Views
- Uploaded on
- Presentation posted in: General

Toward the next generation of earthquake source models by accounting for model prediction error

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Toward the next generation of earthquake source models by accounting for model prediction error

- Zacharie Duputel
- Seismo Lab, GPS division, Caltech

Acknowledgements: Piyush Agram, Mark Simons, Sarah Minson, James Beck,

Pablo Ampuero, Romain Jolivet, Bryan Riel, Michael Aivasis, Hailiang Zhang.

Project : Toward the next generation of source models including realistic statistics of uncertainties

SIV initiative

- Modeling ingredients
- Data:
- Field observations
- Seismology
- Geodesy
- ...

- Theory:
- Source geometry
- Earth model
- ...

- Data:
- Sources of uncertainty
- Observational uncertainty:
- Instrumental noise
- Ambient seismic noise

- Prediction uncertainty:
- Fault geometry
- Earth model

- Observational uncertainty:
- A posteriori distribution

Izmit earthquake (1999)

Slip, m

Depth, km

Slip, m

Depth, km

Single model

Slip, m

Depth, km

Ensemble of

models

2

A reliable stochastic model for the prediction uncertainty

The forward problem

- posterior distribution:

Exact theory

Stochastic (non-deterministic) theory

p(d|m) = δ(d - g( ,m))

p(d|m) = N(d | g( ,m), Cp)

Calculation of Cp based on the physics of the problem: A perturbation approach

Covariance matrix describing uncertainty

in the Earth model parameters

Partial derivatives w.r.t. the elastic parameters (sensitivity kernel)

3

Prediction uncertainty due to the earth model

1000 stochastic realizations

Covariance

Cμ

Cp

Toy model 1: Infinite strike-slip fault

μ1

- Data generated for a layered half-space (dobs)

- 5mm uncorrelated observational noise (→Cd)

- GFs for an homogeneous half-space (→Cp)

- CATMIP bayesian sampler (Minson et al., GJI 2013):

μ2

Synthetic Data + Noise

shallow fault + Layered half-space

Inversion:

Homogeneous half-space

Slip, m

Slip, m

?

μ1

0.9H

0.9H

H

H

μ2

Depth / H

Depth / H

μ2

μ2/μ1 =1.4

2H

2H

Toy model 1: Infinite strike-slip fault

Posterior Mean Model

Input (target) model

Why a smaller misfit does not necessarily indicate a better solution

No Cp (overfitting)

Cp Included (larger residuals)

Depth / H

Depth / H

Slip, m

Slip, m

Displacement, m

Displacement, m

Distance from fault / H

Distance from fault / H

Toy Model 2: Static Finite-fault modeling

Input (target) model

- Finite strike-slip fault
- Top of the fault at 0 km
- South-dipping = 80°
- Data for a layered half-space

Slip, m

Dist. along Dip, km

Dist. along Strike, km

Earth model

Data

Horizontal Disp., m

Vertical Disp., m

North, km

Depth, km

Shear modulus, GPa

East, km

8

Toy Model 2: Static Finite-fault modeling

Input (target) model

- Finite strike-slip fault
- 65 patches, 2 slip components
- 5mm uncorrelated noise(→Cd)
- GFs for an homogeneous half- space (→Cp)

Slip, m

Dist. along Dip, km

Dist. along Strike, km

Earth model

Data

Horizontal Disp., m

Vertical Disp., m

North, km

Depth, km

Model for

Data

Model for

GFs

Shear modulus, GPa

East, km

9

Toy Model 2: Static Finite-fault modeling

Input (target) model - 65 patches average

- Finite strike-slip fault
- 65 patches, 2 slip components
- 5mm uncorrelated noise(→Cd)
- GFs for an homogeneous half- space (→Cp)

Slip, m

Dist. along Dip, km

Dist. along Strike, km

Earth model

Posterior mean model, No Cp

Slip, m

Dist. along Dip, km

Dist. along Strike, km

Depth, km

Posterior mean model, including Cp

Uncertainty on the

shear modulus

Slip, m

Dist. along Dip, km

Dist. along Strike, km

Shear modulus, GPa

10

Conclusion and Perspectives

- Improving source modeling by accounting for realistic uncertainties
- 2 sources of uncertainty
- Observational error
- Modeling uncertainty

- Importance of incorporating realistic covariance components
- More realistic uncertainty estimations
- Improvement of the solution itself

- Accounting for lateral variations
- Improving kinematic source models

- 2 sources of uncertainty

Application to actual data: Mw 7.7 Balochistan earthquake

Jolivet et al., submitted to BSSA

AGU Late breaking session on Tuesday

Toy Model 2: Static Finite-fault modeling

Posterior mean model, including Cp

- Finite strike-slip fault
- 65 patches, 2 slip components
- 5mm uncorrelated noise(→Cd)
- GFs for an homogeneous half- space (→Cp)

Slip, m

Dist. along Dip, km

Dist. along Strike, km

Earth model

Covariance with respect to xr

CpEast(xr), m2

x 104

Depth, km

North, km

Uncertainty on the

shear modulus

xr

Shear modulus, GPa

East, km

13

Toy Model 2: Static Finite-fault modeling

Posterior mean model, including Cp

- Finite strike-slip fault
- 65 patches, 2 slip components
- 5mm uncorrelated noise(→Cd)
- GFs for an homogeneous half- space (→Cp)

Slip, m

Dist. along Dip, km

Dist. along Strike, km

Earth model

Covariance with respect to xr

CpEast(xr), m2

x 104

Depth, km

North, km

xr

Log(μi / μi+1)

East, km

14

Toy model 1: prior: U(-0.5,20)

Posterior Mean Model

Input (target) model

Toy model 1: prior: U(0,20)

Posterior Mean Model

Input (target) model

Toy model including a slip step

Toy model including a slip step

Evolution of m at each beta step

Evolution of Cp at each beta step

Covariance Cμ

1000 realizations

Covariance Cp

1000 realizations

On the importance of Prediction uncertainty

- Observational error:
- Measurements dobs : single realization of a stochastic variable d* which can be described by a probability density p(d*|d) = N(d*|d, Cd)

- Prediction uncertainty: whereΩ = [ μT , φT ]T
- Ωtrue is not known and we work with an approximation
- The prediction uncertainty:
- scales with the with the magnitude of m
- can be described by p(d|m) = N(d | g( ,m), Cp)

- In the Gaussian case, the solution of the problem is given by:

Measurements

Displacement field

Earth

model

Source

geometry

Prior information

Prediction

errors

Measurement

errors

D: Prediction space