• No results found

Synthetic inversions for density in the Earth’s interior from seismic and geodetic data

N/A
N/A
Protected

Academic year: 2022

Share "Synthetic inversions for density in the Earth’s interior from seismic and geodetic data"

Copied!
1
0
0

Bezig met laden.... (Bekijk nu de volledige tekst)

Hele tekst

(1)

Synthetic inversions for density in the Earth’s interior from seismic and geodetic data

Nienke Blom (Universiteit Utrecht - n.a.blom@uu.nl); Christian Böhm, Andreas Fichtner (ETH Zürich)

AGU paper nr.:

S23C-2741

Questions?

find Nienke

Methods - synthetic experiments

We perform synthetic experiments where the target model is known.

P- and S-wave emitting sources (×) lie at 56 km depth, and receivers (o) are located at the top of the domain. L+R boundaries are absorbing.

In this target model density, S-velocity and P-velocity are uncorrelated by design. This is because we want to image density independently without any prior constraints about its geometry and distribution.

We test the following things:

• What is the effect of including prior information about the seismic velocities inside the model?

• To what extent does gravity information help the recovery of density?

Why density?

In plate tectonics and mantle dynamics, density plays a major role in deter- mining the forcings on the systems. Along with knowledge of local seismic velocities, density can help to determine whether tomographic anomalies are of thermal or compositional nature. However, density has thus far only been studied as an independent parameter on the very largest scales us- ing normal modes. Here we try to answer the question: how can we image density?

Figure 1. Tomography traditionally images variations in seismic velocities inside the Earth. However, links between tomography and geodynamics are based on assumptions which do not always hold. In order to understand the dynamics of the Earth’s interior, it is therefore necessary to image density variations directly.

?

tomography geodynamics

-60 -40 -20 0 20 40 60

x [km]

0 1000 2000 3000 4000 5000 6000

z [km]

0

500

1000

1500

2000

2500

S-wave velocity [km/s]

S-wave velocity [km/s]

x [km]

0 1000 2000 3000 4000 5000 6000

z [km]

0

500

1000

1500

2000

2500

-50 -40 -30 -20 -10 0 10 20 30 40 50

density [kg/m3]

density [kg/m3]

-100 -50 0 50 100

x [km]

0 1000 2000 3000 4000 5000 6000

z [km]

0

500

1000

1500

2000

2500

P-wave velocity [km/s]

P-wave velocity [km/s]

target model

Results - how can we best image density?

It is possible to image density using waveform inversion on a global scale, without the use of normal modes. Best re- sults are obtained with an L-BFGS descent method where the lowest frequencies are inverted for first.

-60 -40 -20 0 20 40 60

x [km]

0 1000 2000 3000 4000 5000 6000

z [km]

0

500

1000

1500

2000

2500

S-wave velocity [km/s]

S-wave velocity [km/s]

x [km]

0 1000 2000 3000 4000 5000 6000

z [km]

0

500

1000

1500

2000

2500

-50 -40 -30 -20 -10 0 10 20 30 40 50

density [kg/m3]

density [kg/m3]

-100 -50 0 50 100

x [km]

0 1000 2000 3000 4000 5000 6000

z [km]

0

500

1000

1500

2000

2500

P-wave velocity [km/s]

P-wave velocity [km/s]

starting model

Why does gravity not help?

Gravity measurements are non-unique: infinitely many different structures can yield the same grav- ity response. The sensitivity of gravity falls off with 1/r2. The algorithm will thus put most of the density structure near the top of the domain - which is in- correct here.

Unless prior informa- tion as to the relative distribution of density anomalies is supplied, this cannot be solved.

gravity response

possible subsurface density anomaly

weak

moderate strong

point mass

What if the wrong velocities are fixed?

iteration number

0 20 40 60 80 100 120 140

|current - target model | / | target - background model|

0.7 0.8 0.9 1 1.1 1.2 1.3 1.4

velocities fixed - 50% wrong velocities fixed - 15% wrong velocities fixed - 100% correct

Fixing the seismic velocities to predetermined val- ues improves the density inversion result, but if these velocities are wrong, un-accounted for veloc- ity structure will map into density:

x [km]

0 1000 2000 3000 4000 5000 6000

z [km]

0

500

1000

1500

2000

2500

-50 -40 -30 -20 -10 0 10 20 30 40 50

density [kg/m3]

density [kg/m3]

Figure 7. velocity model 15%

wrong

x [km]

0 1000 2000 3000 4000 5000 6000

z [km]

0

500

1000

1500

2000

2500

-50 -40 -30 -20 -10 0 10 20 30 40 50

density [kg/m3]

density [kg/m3]

Figure 8. velocity model 50%

wrong

Is it beneficial to include gravity into the inver- sion scheme?

No. The use of gravity information only works to de- teriorate the inversion result. This is the case both if gravity is included in the misfit functional (and thus if gravity partial derivatives with respect to density are calculated), and if gravity is used as a so-called “hard constraint”, i.e. if every update is forced to satisfy the gravity data.

x [km]

0 1000 2000 3000 4000 5000 6000

z [km]

0

500

1000

1500

2000

2500

-50 -40 -30 -20 -10 0 10 20 30 40 50

density [kg/m3]

density [kg/m3 ]

Figure 6: density result of an inversion run (121 iterations with frequencies in- creasing every 20 iterations), where no gravity information is used.

data: seis only

prior information:

vs, vpcorrect + fixed

x [km]

0 1000 2000 3000 4000 5000 6000

z [km]

0

500

1000

1500

2000

2500

-50 -40 -30 -20 -10 0 10 20 30 40 50

density [kg/m3]

density [kg/m3 ]

Figure 4. density result of an inversion run (131 iterations at frequencies in- creasing every 20 iterations), where the velocities are correct and fixed.

data: seis + grav prior information:

vs, vpcorrect +fixed

What is the effect of including prior information about the seismic velocities?

Because S- and P-velocities are well known inside the Earth, we can fix these values using a subspace meth- od, and only update density. In this case, the inversion is much more efficient.

iteration number

0 50 100 150

| (model - target model) / target model |1.1 1.2 1.3 1.4 1.5 1.6 1.7

Development of difference from target model

velocities not fixed velocities fixed

Can waveform inversion identify density anom- alies on a global scale?

Yes. The inversion scheme correctly identifies the waveform differences as being caused by density, with only minor contamination in S- and P-velocity struc- ture. This is without the use of normal modes.

x [km]

0 1000 2000 3000 4000 5000 6000

z [km]

0

500

1000

1500

2000

2500

-50 -40 -30 -20 -10 0 10 20 30 40 50

density [kg/m3]

density [kg/m3 ]

Figure 2. density result of an inversion run (144 iterations at frequencies in- creasing every 20 iterations), where the velocities are correct but not fixed.

data: seis + grav prior information:

vs, vpcorrect + free

-60 -40 -20 0 20 40 60

x [km]

0 1000 2000 3000 4000 5000 6000

z [km]

0

500

1000

1500

2000

2500

S-wave velocity [km/s]

S-wave velocity [km/s]

-100 -50 0 50 100

x [km]

0 1000 2000 3000 4000 5000 6000

z [km]

0

500

1000

1500

2000

2500

P-wave velocity [km/s]

P-wave velocity [km/s]

Figure 3. S- and P-velocity results only show minor contaminations from density structure

x [km]

0 1000 2000 3000 4000 5000 6000

z [km]

0

500

1000

1500

2000

2500

typical 1st iteration gravity update

(max: 4.15 kg/m3)

Figure 5. The lower the graph, the closer the result is to the tar- get model. If seismic velocities are fixed, the inversion scheme is much more efficient at recover- ing density.

Referenties

GERELATEERDE DOCUMENTEN

The knowledge which is needed has been described in Chapter 1.3.2 by the following two research questions: “What kind of business process modelling technique is suitable for

Figure 4 (up): a seismogram (z component) for a chosen station for the medium complexity structure experiment (200 km correlation length), black – 3D velocity and 3D

We show that waveform inversion in alternative para- metrisations, using gravity measurements and a fixed v s , v p model, results in significant density recovery in synthetic

A systematic comparison between model results and the seismic strain rates obtained from the National Earthquake Information Center catalog, the geodetic velocity field and

Note the clear quadratic convergence in the iterations – the last column settles to approximately −0.11546 independent of the number of iterations.. This behaviour would not be

In the rendering, the large kernel density field is used for color mapping and the aggregated density fields are used for the illumination.. In the final density map, the color

Next to specifying the energy densities of the non-interacting and physical systems, the energy density of the strong-interaction limit is derived, and implications for local ISI

This is very important because it means that, at least for narrow streams (but see also Sect. 4.1), for the slow observed acceleration ap of DACs (compared to the mean wind