2.5D real waveform and real noise simulation of receiver functions in 3D models
There are several reasons why a real-data receiver function differs from the theoretical receiver function in a 1D model representing the stratification under the seismometer. Main reasons are ambient noise, spectral deficiencies in the impinging P-waveform, and wavefield propagation in laterally va...
Main Authors: | , , |
---|---|
Format: | Conference Object |
Language: | English |
Published: |
2014
|
Subjects: | |
Online Access: | https://pure.au.dk/portal/da/publications/25d-real-waveform-and-real-noise-simulation-of-receiver-functions-in-3d-models(76431297-8e7b-4041-a0c5-b936b2abd52b).html |
id |
ftuniaarhuspubl:oai:pure.atira.dk:publications/76431297-8e7b-4041-a0c5-b936b2abd52b |
---|---|
record_format |
openpolar |
spelling |
ftuniaarhuspubl:oai:pure.atira.dk:publications/76431297-8e7b-4041-a0c5-b936b2abd52b 2023-05-15T16:03:56+02:00 2.5D real waveform and real noise simulation of receiver functions in 3D models Schiffer, Christian Jacobsen, B. H. Balling, N. 2014-04 https://pure.au.dk/portal/da/publications/25d-real-waveform-and-real-noise-simulation-of-receiver-functions-in-3d-models(76431297-8e7b-4041-a0c5-b936b2abd52b).html eng eng info:eu-repo/semantics/restrictedAccess Schiffer , C , Jacobsen , B H & Balling , N 2014 , ' 2.5D real waveform and real noise simulation of receiver functions in 3D models ' , EGU General Assembly 2014 , Vienna , Austria , 27/04/2014 - 02/05/2014 . Seismology Receiver Functions Numerical modelling conferenceObject 2014 ftuniaarhuspubl 2020-07-18T21:33:04Z There are several reasons why a real-data receiver function differs from the theoretical receiver function in a 1D model representing the stratification under the seismometer. Main reasons are ambient noise, spectral deficiencies in the impinging P-waveform, and wavefield propagation in laterally varying velocity variations. We present a rapid "2.5D" modelling approach which takes these aspects into account, so that a given 3D velocity model of the crust and uppermost mantle can be tested more realistically against observed recordings from seismometer arrays. Each recorded event at each seismometer is simulated individually through the following steps: A 2D section is extracted from the 3D model along the direction towards the hypocentre. A properly slanted plane or curved impulsive wavefront is propagated through this 2D section, resulting in noise free and spectrally complete synthetic seismometer data. The real vertical component signal is taken as a proxy of the real impingent wavefield, so by convolution and subsequent addition of real ambient noise recorded just before the P-arrival we get synthetic vertical and horizontal component data which very closely match the spectral signal content and signal to noise ratio of this specific recording. When these realistic synthetic data undergo exactly the same receiver function estimation and subsequent graphical display we get a much more realistic image to compare to the real-data receiver functions. We applied this approach to the Central Fjord area in East Greenland (Schiffer et al., 2013), where a 3D velocity model of crust and uppermost mantle was adjusted to receiver functions from 2 years of seismometer recordings and wide angle crustal profiles (Schlindwein and Jokat, 1999; Voss and Jokat, 2007). Computationally this substitutes tens or hundreds of heavy 3D computations with hundreds or thousands of single-core 2D computations which parallelize very efficiently on common multicore systems. In perspective, the comparison can even be made quantitative, and an iterative inverse 3D model updating would be possible. Furthermore the "2.5D" modelling approach of the in reality 3D problem must be investigated in terms of accuracy of the approximation, in particular with focus on highly 3D structures and multiple phases. Conference Object East Greenland Greenland Aarhus University: Research Greenland |
institution |
Open Polar |
collection |
Aarhus University: Research |
op_collection_id |
ftuniaarhuspubl |
language |
English |
topic |
Seismology Receiver Functions Numerical modelling |
spellingShingle |
Seismology Receiver Functions Numerical modelling Schiffer, Christian Jacobsen, B. H. Balling, N. 2.5D real waveform and real noise simulation of receiver functions in 3D models |
topic_facet |
Seismology Receiver Functions Numerical modelling |
description |
There are several reasons why a real-data receiver function differs from the theoretical receiver function in a 1D model representing the stratification under the seismometer. Main reasons are ambient noise, spectral deficiencies in the impinging P-waveform, and wavefield propagation in laterally varying velocity variations. We present a rapid "2.5D" modelling approach which takes these aspects into account, so that a given 3D velocity model of the crust and uppermost mantle can be tested more realistically against observed recordings from seismometer arrays. Each recorded event at each seismometer is simulated individually through the following steps: A 2D section is extracted from the 3D model along the direction towards the hypocentre. A properly slanted plane or curved impulsive wavefront is propagated through this 2D section, resulting in noise free and spectrally complete synthetic seismometer data. The real vertical component signal is taken as a proxy of the real impingent wavefield, so by convolution and subsequent addition of real ambient noise recorded just before the P-arrival we get synthetic vertical and horizontal component data which very closely match the spectral signal content and signal to noise ratio of this specific recording. When these realistic synthetic data undergo exactly the same receiver function estimation and subsequent graphical display we get a much more realistic image to compare to the real-data receiver functions. We applied this approach to the Central Fjord area in East Greenland (Schiffer et al., 2013), where a 3D velocity model of crust and uppermost mantle was adjusted to receiver functions from 2 years of seismometer recordings and wide angle crustal profiles (Schlindwein and Jokat, 1999; Voss and Jokat, 2007). Computationally this substitutes tens or hundreds of heavy 3D computations with hundreds or thousands of single-core 2D computations which parallelize very efficiently on common multicore systems. In perspective, the comparison can even be made quantitative, and an iterative inverse 3D model updating would be possible. Furthermore the "2.5D" modelling approach of the in reality 3D problem must be investigated in terms of accuracy of the approximation, in particular with focus on highly 3D structures and multiple phases. |
format |
Conference Object |
author |
Schiffer, Christian Jacobsen, B. H. Balling, N. |
author_facet |
Schiffer, Christian Jacobsen, B. H. Balling, N. |
author_sort |
Schiffer, Christian |
title |
2.5D real waveform and real noise simulation of receiver functions in 3D models |
title_short |
2.5D real waveform and real noise simulation of receiver functions in 3D models |
title_full |
2.5D real waveform and real noise simulation of receiver functions in 3D models |
title_fullStr |
2.5D real waveform and real noise simulation of receiver functions in 3D models |
title_full_unstemmed |
2.5D real waveform and real noise simulation of receiver functions in 3D models |
title_sort |
2.5d real waveform and real noise simulation of receiver functions in 3d models |
publishDate |
2014 |
url |
https://pure.au.dk/portal/da/publications/25d-real-waveform-and-real-noise-simulation-of-receiver-functions-in-3d-models(76431297-8e7b-4041-a0c5-b936b2abd52b).html |
geographic |
Greenland |
geographic_facet |
Greenland |
genre |
East Greenland Greenland |
genre_facet |
East Greenland Greenland |
op_source |
Schiffer , C , Jacobsen , B H & Balling , N 2014 , ' 2.5D real waveform and real noise simulation of receiver functions in 3D models ' , EGU General Assembly 2014 , Vienna , Austria , 27/04/2014 - 02/05/2014 . |
op_rights |
info:eu-repo/semantics/restrictedAccess |
_version_ |
1766399634099404800 |