RECEIVER FUNCTIONS analysis for Vs profiles and Moho

  • Slides: 17
Download presentation
RECEIVER FUNCTIONS analysis for Vs profiles and Moho depth estimations Andrea Licciardi (1, 2)

RECEIVER FUNCTIONS analysis for Vs profiles and Moho depth estimations Andrea Licciardi (1, 2) Nicola Piana Agostinetti (1) Geophysics Section, Dublin Institute for Advanced Studies (2) School of Geological Sciences, University College Dublin TABOO Meeting, 23 rd – 24 th April 2013, Ancona

RFs data-set Andrea Licciardi Dublin Institute for Advanced Studies Ancona 23 rd April 2013

RFs data-set Andrea Licciardi Dublin Institute for Advanced Studies Ancona 23 rd April 2013

RFs data-set Num. of recorded earthquakes Num. of selected earthquakes Num. of computed RFs

RFs data-set Num. of recorded earthquakes Num. of selected earthquakes Num. of computed RFs Andrea Licciardi Dublin Institute for Advanced Studies Ancona 23 rd April 2013

1) RFs: Moho depth. The Zhu and Kanamori (2000) method • Based on a

1) RFs: Moho depth. The Zhu and Kanamori (2000) method • Based on a weighted stack of RRFs at the predicted arrival times of Ps, Pp. Ss+Ps. Ps phases • Grid search in the H-Vp/Vs domain • For the correct values of H and Vp/Vs the stacking function is max • RFs are computed with a frequency domain deconvolution (Di Bona, 1998) • Fq ≈ 1 Hz. • Ver. Res ≈ 1 km Andrea Licciardi Dublin Institute for Advanced Studies Ancona 23 rd April 2013

1) RFs: Moho depth Andrea Licciardi Dublin Institute for Advanced Studies Ancona 23 rd

1) RFs: Moho depth Andrea Licciardi Dublin Institute for Advanced Studies Ancona 23 rd April 2013

1) RFs: Moho depth Piana Agostinetti et al. 2011 Pietralunga Umbertide Promano Andrea Licciardi

1) RFs: Moho depth Piana Agostinetti et al. 2011 Pietralunga Umbertide Promano Andrea Licciardi Trasimeno Lake Dublin Institute for Advanced Studies Piccione Sigillo Ancona 23 rd April 2013

2) RFs: mapping the ATF Andrea Licciardi Dublin Institute for Advanced Studies Ancona 23

2) RFs: mapping the ATF Andrea Licciardi Dublin Institute for Advanced Studies Ancona 23 rd April 2013

2) RFs: mapping the ATF Direct P polarity reversal in the very first second

2) RFs: mapping the ATF Direct P polarity reversal in the very first second of the Radial and Transverse Receiver Functions • DIPPING INTERFACE • HIGH VELOCITY CONTRAST Andrea Licciardi Dublin Institute for Advanced Studies Ancona 23 rd April 2013

Evidence of ATF from RFs BAZ Dist. The BAZ of the min amplitude on

Evidence of ATF from RFs BAZ Dist. The BAZ of the min amplitude on the RADIAL and of the zero energy on the TRANSVERSE gives the DIP-DIRECTION of the interface Andrea Licciardi Dublin Institute for Advanced Studies Ancona 23 rd April 2013

2) RFs: mapping the ATF ATMC ATMI DIP-DIRECTION STRIKE N 002 BAZ Dist. Andrea

2) RFs: mapping the ATF ATMC ATMI DIP-DIRECTION STRIKE N 002 BAZ Dist. Andrea Licciardi Dublin Institute for Advanced Studies AT 01 BAZ Dist. Ancona 23 rd April 2013

2) RFs: mapping the ATF AT 01 N 002 Andrea Licciardi Dublin Institute for

2) RFs: mapping the ATF AT 01 N 002 Andrea Licciardi Dublin Institute for Advanced Studies Ancona 23 rd April 2013

3) High Frequency RF: comparison with borehole data Freq ~ 1 Hz = V-Res:

3) High Frequency RF: comparison with borehole data Freq ~ 1 Hz = V-Res: 1 km AT 05: same location as Pozzo San Donato ATVO: Monte Valentino, ~1 -2 km from Monte Civitello Focus on Radial RF for isotropic S-vel profile Andrea Licciardi Dublin Institute for Advanced Studies Ancona 23 rd April 2013

3) High Frequency RF: comparison with borehole data Harmonic decomposition of the RF data-set

3) High Frequency RF: comparison with borehole data Harmonic decomposition of the RF data-set (decoupling of isotropic/anisotropic components) K=0 analysis (i. e. isotropic) Max freq. ~ 4 Hz , V-res ~ 0. 25 km (? ) Sharp negative arrival = velocity inversion at depth (within the first 5 km) Andrea Licciardi Dublin Institute for Advanced Studies Ancona 23 rd April 2013

3) High Frequency RF: comparison with borehole data Posterior mean Vp and Vs Maximum

3) High Frequency RF: comparison with borehole data Posterior mean Vp and Vs Maximum probability for an interface at depth 1 D Posterior distribution of seismic interfaces at depth Andrea Licciardi Fit between mean synthetic and observed RF 1 D Posterior distribution for S-velocity Dublin Institute for Advanced Studies Ancona 23 rd April 2013

3) High Frequency RF: comparison with borehole data Litho-discontinuities from Borehole data Interpretation (e.

3) High Frequency RF: comparison with borehole data Litho-discontinuities from Borehole data Interpretation (e. g. Mirabella et al 2011) Vs Discontinuities from RF data (1) (2) (3) (5) (4) (6) (1) ~200 m between top Evaporites (2) Intra-Evaporites discontinuities retrieved (3) Top of the Acoustic Basement missed by ~200 m, but presence of fine -layering (i. e. Ligurids) on top. Basement is not low Vs as expected (4) Intra-basement discontinuity retrieved (5) No discontinuity between Basement and Evaporites (i. e. Evaporites have low-Vs) (6) Very low-Vs top of the Acoustic Basement Andrea Licciardi Dublin Institute for Advanced Studies Ancona 23 rd April 2013

3) High Frequency RF: comparison with borehole data - AT 05 0 m Careful

3) High Frequency RF: comparison with borehole data - AT 05 0 m Careful inspection of borehole data should be used to interpret the seismic discontinuities retrieved by RF analysis For example, for AT 05, the Vs discontinuity at about 3800 m might be associated to a sharp resistivity variation “Argilloscisti” sequence ~3030 m (4) ~4483 m

CONCLUSIONS 1. Results from RFs analysis show a complex Moho topography in between the

CONCLUSIONS 1. Results from RFs analysis show a complex Moho topography in between the Tyrrhenian and Adriatic domains; 2. RFs can be used to map faults with high velocity contrast in the shallow crust; 3. Velocitity profiles from high frequency RFs can be directly compared with borehole data providing additional information for geological interpretations; OUTLOOK: 1. Results from station ATVO will be compared with the available data from Monte Civitello borehole; 2. 1 D Vs velocity profiles will be computed at each station. Andrea Licciardi Dublin Institute for Advanced Studies Ancona 23 rd April 2013

MODEL OF HUMAN OCCUPATION MOHO THE BEGINNING MOHOMODEL OF HUMAN OCCUPATION MOHO THE BEGINNING MOHO
What a Receiver Really Sees Receiver 1 ReceiverWhat a Receiver Really Sees Receiver 1 Receiver
Profiles and Assessment Profiles Profiles New Link MenuProfiles and Assessment Profiles Profiles New Link Menu
Exercise 04 Receiver functions will be determined ReceiverExercise 04 Receiver functions will be determined Receiver
PODREDUMBRE O MOHO GRIS Botrytis cinerea FITOPATOLOGIA AGRICOLAPODREDUMBRE O MOHO GRIS Botrytis cinerea FITOPATOLOGIA AGRICOLA
Moho The plate boundaries Plate Boundary Types DivergentMoho The plate boundaries Plate Boundary Types Divergent
The GALILEO GPS Software Receiver Company Software ReceiverThe GALILEO GPS Software Receiver Company Software Receiver
Review of Important AM Receiver Concepts 1 ReceiverReview of Important AM Receiver Concepts 1 Receiver
Cosc 54730 Broadcast Receiver Broadcast receiver A broadcastCosc 54730 Broadcast Receiver Broadcast receiver A broadcast
Intent Broadcast Receiver l Intent l Broadcast ReceiverIntent Broadcast Receiver l Intent l Broadcast Receiver
Tecom IP Receiver Overview The Tecom IP ReceiverTecom IP Receiver Overview The Tecom IP Receiver
Conventional Eddycurrent NDE Exciter Receiver Electronics Exciter ReceiverConventional Eddycurrent NDE Exciter Receiver Electronics Exciter Receiver
Using the Receiver Workbench Concept Using the ReceiverUsing the Receiver Workbench Concept Using the Receiver
IEEE 1394 Internet encoder transmitter receiver NAIST receiverIEEE 1394 Internet encoder transmitter receiver NAIST receiver
UniDirectionalLink UDL Receiver Feeder BiDirectionalLink BDL UDL ReceiverUniDirectionalLink UDL Receiver Feeder BiDirectionalLink BDL UDL Receiver
Fiber Optic Receiver A fiber optic receiver isFiber Optic Receiver A fiber optic receiver is