On the Limitation of Receiver Functions Method: Beyond Conventional Assumptions & Advanced Inversion Techniques
|
|
- Annabel Warren
- 6 years ago
- Views:
Transcription
1 On the Limitation of Receiver Functions Method: Beyond Conventional Assumptions & Advanced Inversion Techniques Hrvoje Tkalčić RSES, ANU Acknowledgment: RSES Seismology & Mathematical Geophysics Group members
2 Some passive-source receiverbased methods Teleseismic travel times Teleseismic receiver functions (P and S) Body wave earthquake interferometry Ambient noise dispersion Interstation method surface wave dispersion..
3 Basin structure and geometry from nuclear blasts waveforms and teleseismic travel times Rodgers et al. PAGEOPH 26; Tkalčić et al., BSSA 28
4 The benefits of RFs + SWD Earth Vs structure can be inverted using: 1. Receiver functions (RF) 2. Surface wave dispersion curves 3. RF + dispersion curves (jointly) or other datasets Different approaches to modeling Forward modeling Linearized inversion Grid-search Non-linear inversion with optimization Multi-step approach
5 IRFFM (Interactive RF Forward Modeling) Multi-step approach Linearized Inversion Tkalčić et al., JGR 26 Chen et al., JGR 21 Stipčević et al., GJI 211 Tkalčić et al., GJI 211; GJI 212
6 Lithospheric structure of Saudi Arabia, China, Australia & Croatia from multi-step modeling of RFs and SWs SE China (RFs combined with tomography) Croatia and Adriatic Sea SE Australia (RFs combined with ambient noise)
7 Advantages and limitations of Advantages RFs A way to invert for Vs structure under a single station Sensitive to gradients (discontinuities) in Vs velocity A needed complement to crustal tomography RF + SW dispersion curves (jointly) or other datasets Limitations of conventional methods Information limited to a volume beneath a single station Insensitive to absolute velocity unless SW are added Simplifications/assumptions often cannot explain real Earth (1. lack of data, 2. anisotropy, 3. dipping layers, 4. non-uniqueness and noise in the data)
8 1. Exploiting seismic signal and noise in an aseismic environment to constrain crustal structure c. CAPRA station CP8 Initial IRFFM2 model RF only inversion Best models from joint inversion Preferred model from joint inversion Observed RF Time (s) Depth (km) Velocity (km/s) Group velocity data from noise Phase velocity data from interstation method Group velocity dispersion curve from RF modelling Phase velocity dispersion curve from RF modelling Group velocity dispersion curve from joint modelling Phase velocity dispersion curve from joint modelling Velocity (km/s) d. CAPRA station CP Period (s) Initial IRFFM model average model: RF analysis Best models: joint inversion Final average model: joint inversion time Young et al., GJI 212 depth velocity velocity Group velocity data from noise Phase velocity data from interstation method Group velocity dispersion curve from RF modelling Phase velocity dispersion curve from RF modelling Group velocity dispersion curve from joint modelling Phase velocity dispersion curve from joint modelling period
9 2. Dipping Moho From Stipčević, PhD Thesis, paper in preparation Moho depth determined using H κ (above) and NA (right) method
10 From Stipčević, PhD Thesis, paper in preparation
11 Stations with similar results obtained using H-κ and NA methods Stipčević et al., in preparation
12 Stations for which there is a large difference between the H-κ i NA results Stipčević et al., in preparation
13 Moho dip determined using NA algorithm Stipčević et al., in preparation
14 Dipping Moho synthetic experiment Starting model : Two 2km thick layers in the crust Moho at 4 km with 2 dip & 27 strike Synthetics are calculated using Fredrickson and Bostock method assuming an isotropic medium, and synthetic RFs are determined by deconvolution. These synthetic RF data are then linearly stacked and inverted for Earth structure using NA method introduced in Exercise 6.
15 Dipping Moho synthetic experiment Starting model: Two 2km thick layers in the crust Moho at 4 km with 2 dip & 27 strike Synthetics are calculated using the Fredrickson and Bostock method assuming an isotropic medium and synthetic RFs are determined by deconvolution.
16 Dipping Moho Synthetic Experiment 2 km 36 km Earth model This is a result of the NA inversion when laterally homogeneous horizontal layers are assumed. It is also assumed that the Moho is horizontal (not dipping). (Exercise 1 in synopsis)
17 Dipping Moho Synthetic Experiment 2 km 4km Earth model Moho: D~17 ; S~262 Now inverting for the Moho dip and orientation
18 3. A multi-step approach including polarization anisotropy Tkalčić et al., JGR 26
19 4. Non-uniqueness etc. Different approaches to inverse problems Forward modeling Linearized inversion Grid-search Non-linear inversion with optimization Multi-step approach Non-linear inversion with the Bayesian framework Transdimensional Bayes framework hierarchical Bayes theorem: p(m d obs ) p(d obs m) p(m)
20 The importance of knowing the data noise in trans-dimensional formulation! est = x ! est = 3 Number of samples Number of cells x ! est = 8 Number of Samples Number of cells 1 x Number of samples Number of cells
21 Hierarchical Models Relationship between data noise and model complexity Treating data noise σ as an unknown in the problem Data noise is uncorrelated Likelihood function unknown p(d m) 1 ( 2πσ 2 ) & exp ( d g(m) N '( 2σ 2 2 ) + * + unknown unknown Data noise = measurement uncertainty + modeling uncertainty
22 Covariance Matrix of Noise in Data Data noise is correlated Misfit Φ(m) = [d g(m)] T C D 1 [d g(m)] Likelihood p(m d) = 1 (2π) N C D % exp Φ(m) ( & ' 2 ) *
23 Noise Parameterization How do we parameterize C D? $ 1 r r 2... r N 1 ' & ) & r 1 r r N 2 ) C D = σ 2 & r 2 r 1 r N 3 ) & ) & ) %& r N 1 r N 2 r N () 2 parameters : Magnitude of noise σ Correlation of noise r
24 Synthetic experiment Receiver Function !.1 Noisy Receiver Func8on!.2! Time (s) Magnitude and correlation of noise are treated as unknowns Depth (km) Depth(km) True model Solution is a large ensemble of models distributed according to the target distribution Vs (km/s) Vs(km/s) Bodin et al., JGR 212
25 Synthetic experiment Receiver Function !.1!.2! True model Mean Noisy Receiver Func8on Time (s) Number of models Magnitude and correlation of noise are unknown Depth(km) Different ways to look at the solution Vs(km/s) Vs(km/s) 6 p(transition) Bodin et al., JGR 212
26 Synthetic experiment Receiver Function !.1!.2! True model Mean Noisy Receiver Func8on Time (s) Solu8on for maximum number of models Magnitude and correlation of noise are unknown 1 1 Probability distribution at 31 km Depth(km) Depth(km) Ensemble Solution True Model Average Solution Maximun Solution Vs(km/s) Vs(km/s) Vs (Km/s)
27 Synthetic experiment Depth(km) Receiver Function !.1!.2! True model 1 Noisy Receiver Func8on Time (s) Number of layers Number of models Prior Posterior True Model Algorithm is able to recover the complexity of the model and the level of data noise Vs(km/s).1.2.3! r
28 Synthetic experiment Receiver Function !.1!.2! True model Noisy Receiver Func8on Time (s) Magnitude and correlation of noise are unknown Trade-offs between parameters Depth(km) Moho Depth Moho depth (km) Vs(km/s) Vs(km/s) Vs in the last layer of the crust
29 Application to field data Bodin et al., JGR 212
30 Application to Joint Inversion Dispersion curve for Rayleigh waves Group Velocity Receiver Function Period (s) Receiver Function!.2! Time(s) d = d RF + d dis " C = C RF % D D $ ' C dis # D & Algorithm naturally weights the information brought by each data type.
31 Application to Joint Inversion (synthetic test) Receiver Function Dispersion Joint Inversion Depth(km) 3 Depth(km) 3 4 Depth(km) Vs(km/s) p(discontinuity) Vs(km/s) p(discontinuity) Vs(km/s) Bodin et al., JGR 212
32 Probability Application to Joint Inversion WOMBAT data from SE Australia: RFs + ambient noise dispersion 1 Only RF Only SWD Joint (RF+SWD) Average models Depth(km) 15 Depth(km) 15 Depth(km) 15 Depth(km) Vs(km/s) Vs(km/s) Vs(km/s) RF SWD Joint Vs(km/s) Bodin et al., JGR 212
33 Future research Modelling multiple geophysical datasets will be approached through the transdimensional hieararchical Bayesian framework, where the number of free parameters and the data noise will be treated as unknowns in the inversion. Various simplifications that hinder the progress in crustal and lithospheric imaging using passive-source data and permanent/temporary seismic receivers will be gradually incorporated in the Bayesian inversion strategy this includes, but is not limited to: anisotropy, dipping layers and 3D structure, noisy data, etc. This is a general strategy that can be applied to other types of inverse problems in Earth Science and for imaging various parts of the Earth s crust.
Transdimensional Markov Chain Monte Carlo Methods. Jesse Kolb, Vedran Lekić (Univ. of MD) Supervisor: Kris Innanen
Transdimensional Markov Chain Monte Carlo Methods Jesse Kolb, Vedran Lekić (Univ. of MD) Supervisor: Kris Innanen Motivation for Different Inversion Technique Inversion techniques typically provide a single
More informationBayesian Monte Carlo Inversion of Multiple Data Sets to Image the Crust and Uppermost Mantle Beneath the Continental United States
University of Colorado, Boulder CU Scholar Physics Graduate Theses & Dissertations Physics Spring 1-1-2014 Bayesian Monte Carlo Inversion of Multiple Data Sets to Image the Crust and Uppermost Mantle Beneath
More informationGeophysical Journal International. Joint inversion of surface wave dispersion and receiver functions: A Bayesian Monte-Carlo approach
Joint inversion of surface wave dispersion and receiver functions: A Bayesian Monte-Carlo approach Journal: Geophysical Journal International Manuscript ID: GJI-S--0.R Manuscript Type: Research Paper Date
More informationJoint inversion of InSAR and broadband teleseismic waveform data with ABIC: application to the 1997 Manyi, Tibet earthquake
Joint inversion of InSAR and broadband teleseismic waveform data with ABIC: application to the 1997 Manyi, Tibet earthquake Gareth Funning 1, Yukitoshi Fukahata 2, Yuji Yagi 3 & Barry Parsons 4 1 University
More information2008 Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
STRUCTURE OF THE KOREAN PENINSULA FROM WAVEFORM TRAVEL-TIME ANALYSIS Roland Gritto 1, Jacob E. Siegel 1, and Winston W. Chan 2 Array Information Technology 1 and Harris Corporation 2 Sponsored by Air Force
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2012) 189, 1536 1556 doi: 10.1111/j.1365-246X.2012.05414.x Transdimensional tomography with unknown data noise Thomas Bodin, 1, Malcolm Sambridge, 1
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2012) 189, 1681 1700 doi: 10.1111/j.1365-246X.2012.05442.x Multistep modelling of receiver-based seismic and ambient noise data from WOMBAT array: crustal
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2012) 189, 1681 1700 doi: 10.1111/j.1365-246X.2012.05442.x Multistep modelling of receiver-based seismic and ambient noise data from WOMBAT array: crustal
More informationOn the Nature of the P-Wave Velocity Gradient in the Inner Core beneath Central America
Journal of Earth Science, Vol. 24, No. 5, p. 699 705, October 2013 ISSN 1674-487X Printed in China DOI: 10.1007/s12583-013-0365-7 On the Nature of the P-Wave Velocity Gradient in the Inner Core beneath
More informationContents of this file
Geophysical Research Letters Supporting Information for Intraplate volcanism controlled by back-arc and continental structures in NE Asia inferred from trans-dimensional ambient noise tomography Seongryong
More information27th Seismic Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
SEISMIC SOURCE AND PATH CALIBRATION IN THE KOREAN PENINSULA, YELLOW SEA, AND NORTHEAST CHINA Robert B. Herrmann 1, Young-Soo Jeon 1, William R. Walter 2, and Michael E. Pasyanos 2 Saint Louis University
More informationData Repository Item
Data Repository Item 2009003 An abrupt transition from magma-starved to magma-rich rifting in the eastern Black Sea Donna J. Shillington, Caroline L. Scott, Timothy A. Minshull, Rosemary A. Edwards, Peter
More informationRegional 3D velocity structure
Seismic imaging: Regional 3D velocity structure 3 Seismic anisotropy courtesy of Ed Garnero This is birefringence: tells us about the preferential orientation of minerals 1 Anisotropy beneath the East
More informationTaiwan Moho Discontinuity Reference Model
Taiwan Moho Discontinuity Reference Model Receiver Function Studies How-Wei Chen ( 陳浩維 ) Colleagues: H. L. Wang, L. Zhu, John Y.H. Duan and K. J. Shyu RF Based TMDRM Studies Related WPGM Presentations
More informationInferring the Orientation of Anisotropy from Receiver Functions
GA DATA RPOITORY irth et al. Methodology: Receiver Functions To interrogate the anisotropic character of the MLD beneath the central U.., we computed receiver functions at four long running seismic stations
More informationImaging Moho topography beneath the Alps by multdisciplinary seismic tomography
Imaging Moho topography beneath the Alps by multdisciplinary seismic tomography Edi Kissling ETH Zürich SPP short course February 1+2, 218, Berlin, Germany Alpine Moho map from CSS Moho uncertainty derived
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2012) 188, 144 164 doi: 10.1111/j.1365-246X.2011.05249.x Estimating crustal seismic anisotropy with a joint analysis of radial and transverse receiver
More informationAmbient Noise Tomography in the Western US using Data from the EarthScope/USArray Transportable Array
Ambient Noise Tomography in the Western US using Data from the EarthScope/USArray Transportable Array Michael H. Ritzwoller Center for Imaging the Earth s Interior Department of Physics University of Colorado
More information2008 Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
STRUCTURE OF THE KOREAN PENINSULA FROM WAVEFORM TRAVEL-TIME ANALYSIS Roland Gritto 1, Jacob E. Siegel 1, and Winston W. Chan 2 Array Information Technology 1 and Harris Corporation 2 Sponsored by Air Force
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2017) 208, 10 23 Advance Access publication 2016 October 8 GJI Seismology doi: 10.1093/gji/ggw383 A Markov chain Monte Carlo with Gibbs sampling approach
More informationGlobal 1-D Earth models
Topic Global 1-D Earth models Compiled Peter Bormann (formerly GeoForschungsZentrum Potsdam, Telegrafenberg, by D-14473 Potsdam, Germany); E-mail: pb65@gmx.net Version March 2002 Below, data and background
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (014) 198, 1431 1437 GJI Seismology doi: 10.1093/gji/ggu13 Estimation of shear velocity contrast from transmitted P s amplitude variation with ray-parameter
More informationUncertainty analysis for the integration of seismic and CSEM data Myoung Jae Kwon & Roel Snieder, Center for Wave Phenomena, Colorado School of Mines
Myoung Jae Kwon & Roel Snieder, Center for Wave Phenomena, Colorado School of Mines Summary Geophysical inverse problems consist of three stages: the forward problem, optimization, and appraisal. We study
More informationEvidence of an axial magma chamber beneath the ultraslow spreading Southwest Indian Ridge
GSA Data Repository 176 1 5 6 7 9 1 11 1 SUPPLEMENTARY MATERIAL FOR: Evidence of an axial magma chamber beneath the ultraslow spreading Southwest Indian Ridge Hanchao Jian 1,, Satish C. Singh *, Yongshun
More information29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
MODELING TRAVEL-TIME CORRELATIONS BASED ON SENSITIVITY KERNELS AND CORRELATED VELOCITY ANOMALIES William L. Rodi 1 and Stephen C. Myers 2 Massachusetts Institute of Technology 1 and Lawrence Livermore
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (215) 23, 1832 1846 GJI Seismology doi: 1.193/gji/ggv41 A Bayesian approach to assess the importance of crustal corrections in global anisotropic surface
More information29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies MODELING P WAVE MULTIPATHING IN SOUTHEAST ASIA
MODELING P WAVE MULTIPATHING IN SOUTHEAST ASIA Ali Fatehi and Keith D. Koper Saint Louis University Sponsored by the Air Force Research Laboratory ABSTRACT Contract No. FA8718-06-C-0003 We have used data
More informationSUPPLEMENTARY INFORMATION
Inability of additional parameters to resolve the Rayleigh-Love discrepancy Radial anisotropy is introduced to resolve the Rayleigh-Love misfit discrepancy that exists across large regions of the western
More informationInverting Ground Motion from a Seismometer Array to Obtain the Vertical Component of Rotation: A Test Using Data from Explosions
Inverting Ground Motion from a Seismometer Array to Obtain the Vertical Component of Rotation: A Test Using Data from Explosions Wu-Cheng Chi 1, W. H. K. Lee 2, C. J. Lin 1, John Aston 3, and C. C. Liu
More informationThe Earth s Structure from Travel Times
from Travel Times Spherically symmetric structure: PREM - Crustal Structure - Upper Mantle structure Phase transitions Anisotropy - Lower Mantle Structure D D - Structure of of the Outer and Inner Core
More informationReceiver Function Inversion
Receiver Function Inversion Advanced Studies Institute on Seismological Research Kuwait City, Kuwait - January 19-22, 2013 Jordi Julià Universidade Federal do Rio Grande do Norte, Brasil Outline Introduction
More informationDR
DR2003071 0 0 270 0 30 0 90 0 60 0 120 0 150 0 90 0 180 0 180 0 A) RadialReceiverFunctions B ackazimuth (in degrees relative to north) -135-90 -45 0 45 90 135 180-5.0-2.5 Tangential R eceiver Functions
More informationSurface Wave Tomography
Surface Wave Tomography Alexandra Mauerberger gassner@gfz-potsdam.de GEO-DEEP9300 course at CEED, University of Oslo Nov. 6 - Nov. 10 2017 Outline Recap General Methods Techniques to Derive Tomography
More information29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
SASIA3D: A 3-D CRUST AND UPPER-MANTLE VELOCITY MODEL OF SOUTH ASIA DERIVED FROM JOINT INVERSION OF P-WAVE TRAVEL TIMES AND SURFACE-WAVE DISPERSION DATA William L. Rodi 1 and Delaine T. Reiter 2 Massachusetts
More informationTomographic imaging of P wave velocity structure beneath the region around Beijing
403 Doi: 10.1007/s11589-009-0403-9 Tomographic imaging of P wave velocity structure beneath the region around Beijing Zhifeng Ding Xiaofeng Zhou Yan Wu Guiyin Li and Hong Zhang Institute of Geophysics,
More informationCURRICULUM VITAE WEISEN SHEN EDUCATION
CURRICULUM VITAE WEISEN SHEN Center for Imaging the Earth s Interior, Department of Physics, CU Boulder Boulder, CO, 80309 Research ID: J-3969-2012 http://www.researcherid.com/rid/j-3969-2012 Citizenship:
More informationConstraints on Shallow Low-Viscosity Earth Layers from Future GOCE Data
Constraints on Shallow Low-Viscosity Earth Layers from Future GOCE Data Hugo Schotman 1,2, Bert Vermeersen 2, Pieter Visser 2 1 2 3 rd International GOCE User Workshop, ESA Esrin, 7 November 2006 glacial-isostatic
More informationSeismogram Interpretation. Seismogram Interpretation
Travel times in the Earth Ray paths, phases and their name Wavefields in the Earth: SH waves, P-SV waves Seismic Tomography Receiver Functions Seismogram Example Long-period transverse displacement for
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2014) 197, 443 457 Advance Access publication 2014 January 24 doi: 10.1093/gji/ggt515 Estimating sedimentary and crustal structure using wavefield continuation:
More informationWhen one of things that you don t know is the number of things that you don t know
When one of things that you don t know is the number of things that you don t know Malcolm Sambridge Research School of Earth Sciences, Australian National University Canberra, Australia. Sambridge. M.,
More informationThe crustal structure beneath SE Romania from teleseismic receiver functions
Geophys. J. Int. (25) 163, 238 251 doi: 1.1111/j.1365-2X.25.2715.x The crustal structure beneath SE Romania from teleseismic receiver functions T. Diehl, J. R. R. Ritter and the CALIXTO Group Geophysical
More informationSeismic tomography: Art or science?
Seismic tomography: Art or science? Frederik J Simons Princeton University What s inside the Earth? www.usgs.gov www.glencoe.com Dalton, Nature 2003 Only seismic waves have actually been there, done that
More informationChangbaishan volcanism in northeast China linked to subduction-induced mantle upwelling
SUPPLEMENTARY INFORMATION DOI: 10.1038/NGEO2166 Changbaishan volcanism in northeast China linked to subduction-induced mantle upwelling Youcai Tang, Masayuki Obayashi, Fenglin Niu, Stephen P. Grand, Yongshun
More informationAPPLICATION OF RECEIVER FUNCTION TECHNIQUE TO WESTERN TURKEY
APPLICATION OF RECEIVER FUNCTION TECHNIQUE TO WESTERN TURKEY Timur TEZEL Supervisor: Takuo SHIBUTANI MEE07169 ABSTRACT In this study I tried to determine the shear wave velocity structure in the crust
More informationReconstruction of Subsurface Structures from Ambient Seismic Noise in Jakarta Area
Reconstruction of Subsurface Structures from Ambient Seismic Noise in Jakarta Area Sri Widiyantoro 1 Andri Dian Nugraha 1 Zulfakriza 2 Bayu Pranata 3 1 Global Geophysics Research Group, Faculty of Mining
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2015) 203, 1302 1318 GJI Seismology doi: 10.1093/gji/ggv368 Ultra-low velocity zones beneath the Philippine and Tasman Seas revealed by a trans-dimensional
More informationWhat does Seismic Anisotropy tell us about the Lithosphere-Asthenosphere Boundary?
What does Seismic Anisotropy tell us about the Lithosphere-Asthenosphere Boundary? Jean-Paul Montagner (1), Gael Burgos (1), Eric Beucler (2), Antoine Mocquet (2) and Yann Capdeville (2), Mathias Obrebski
More information29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies ADVANCED WAVEFORM SIMULATION FOR SEISMIC MONITORING EVENTS
ABSTRACT ADVANCED WAVEFORM SIMULATION FOR SEISMIC MONITORING EVENTS Don V. Helmberger 1, Jeroen Tromp 1, and Arthur J. Rodgers 2 California Institute of Technology 1 and Lawrence Livermore National Laboratory
More informationPure and Applied Geophysics Calculation of station-representative isotropic receiver functions
Pure and Applied Geophysics Calculation of station-representative isotropic receiver functions --Manuscript Draft-- Manuscript Number: Full Title: Article Type: Keywords: Corresponding Author: PAAG-D--00R
More informationHigh Resolution Imaging of Fault Zone Properties
Annual Report on 1998-99 Studies, Southern California Earthquake Center High Resolution Imaging of Fault Zone Properties Yehuda Ben-Zion Department of Earth Sciences, University of Southern California
More information28th Seismic Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
CRUSTAL AND UPPER-MANTLE P- AND S-VELOCITY STRUCTURE IN CENTRAL AND SOUTHERN ASIA FROM JOINT BODY- AND SURFACE-WAVE INVERSION Delaine Reiter 1 and William Rodi 2 Weston Geophysical Corporation 1 and Earth
More informationP191 Bayesian Linearized AVAZ Inversion in HTI Fractured Media
P9 Bayesian Linearized AAZ Inversion in HI Fractured Media L. Zhao* (University of Houston), J. Geng (ongji University), D. Han (University of Houston) & M. Nasser (Maersk Oil Houston Inc.) SUMMARY A new
More informationSurface wave array tomography in SE Tibet from ambient seismic noise and two-station analysis II. Crustal and upper-mantle structure
Geophys. J. Int. (2008) doi: 10.1111/j.1365-246X.2007.03696.x Surface wave array tomography in SE Tibet from ambient seismic noise and two-station analysis II. Crustal and upper-mantle structure Huajian
More informationBayesian ISOLA: new tool for automated centroid moment tensor inversion
Bayesian ISOLA: new tool for automated centroid moment tensor inversion A CMT catalog automatically generated for Swiss data Jiří Vackář 1,2 Jan Burjánek 2 František Gallovič 1 Jiří Zahradník 1 John Clinton
More informationThe stratification of seismic azimuthal anisotropy in the western US
The stratification of seismic azimuthal anisotropy in the western US Fan-Chi Lin 1, Michael H. Ritzwoller 1, Yingjie Yang 1, Morgan P. Moschetti 1, and Matthew J. Fouch 2 1 - Center for Imaging the Earth's
More informationGlobal geophysics and wave propagation
Global geophysics and wave propagation Reading: Fowler p76 83 Remote sensing Geophysical methods Seismology Gravity and bathymetry Magnetics Heat flow Seismology: Directly samples the physical properties
More informationCatherine A. Rychert University of Bristol
Catherine A. Rychert University of Bristol Thanks to: James O. S. Hammond, J. Michael Kendall, Nicholas Harmon, Derek Keir, Cindy Ebinger, Atalay Ayele, Ian Bastow, Graham Stuart, Manahloh Belachew What
More informationImproving Regional Seismic Event Location Through Calibration of the International Monitoring System
UCRL-JC-133019 PREPRINT Improving Regional Seismic Event Location Through Calibration of the International Monitoring System C.A. Schultz S.C. Myers M. Flanagan J. Swenson W. Hanley M. Pasyanos This paper
More informationMultiple Scenario Inversion of Reflection Seismic Prestack Data
Downloaded from orbit.dtu.dk on: Nov 28, 2018 Multiple Scenario Inversion of Reflection Seismic Prestack Data Hansen, Thomas Mejer; Cordua, Knud Skou; Mosegaard, Klaus Publication date: 2013 Document Version
More information27th Seismic Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
P AND S WAVE VELOCITY STRUCTURE OF THE CRUST AND UPPER MANTLE UNDER CHINA AND SURROUNDING AREAS FROM BODY AND SURFACE WAVE TOMOGRAPHY M. Nafi Toksöz, Robert D. Van der Hilst, Youshun Sun, and Chang Li
More information29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies ENHANCEMENTS OF GEOPHYSICAL MODELS FOR MONITORING
ENHANCEMENTS OF GEOPHYSICAL MODELS FOR MONITORING Michael E. Pasyanos, William R. Walter, Megan P. Flanagan Rengin Gok, Stephen C. Myers, Kathleen M. Dyer Lawrence Livermore National Laboratory Sponsored
More informationMeasuring seismicity in the Groningen Field. Bernard Dost, Elmer Ruigrok, Jesper Spetzler, Gert-Jan van den Hazel, Jordi Domingo
Measuring seismicity in the Groningen Field Bernard Dost, Elmer Ruigrok, Jesper Spetzler, Gert-Jan van den Hazel, Jordi Domingo Monitoring induced seismicity in the Netherlands. Instrumentation and network
More informationEdinburgh Research Explorer
Edinburgh Research Explorer Global crustal thickness from neural network inversion of surface wave data Citation for published version: Meier, U, Curtis, A & Trampert, J 27, 'Global crustal thickness from
More informationAuxiliary Material. Subduction of oceanic asthenosphere: evidence from sub-slab seismic anisotropy. Teh-Ru Alex Song, Hitoshi Kawakatsu
Auxiliary Material Subduction of oceanic asthenosphere: evidence from sub-slab seismic anisotropy Teh-Ru Alex Song, Hitoshi Kawakatsu correspondence to: tehrusong@gmail.com 1 1. The effect of anisotropy
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2011) 185, 1103 1119 doi: 10.1111/j.1365-246X.2011.05004.x Crustal and uppermost mantle structure beneath the External Dinarides, Croatia, determined
More informationCURRICULUM VITAE Weisen Shen 沈伟森
CURRICULUM VITAE Weisen Shen 沈伟森 BORN: August 8, 1986, Meishan, China; CITIZENSHIP: China ADDRESS: Department of Earth and Planetary Sciences Washington University in St Louis One Brookings Drive St Louis,
More informationGEOPHYSICAL RESEARCH LETTERS, VOL. 37, L02304, doi: /2009gl041835, 2010
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 37,, doi:10.1029/2009gl041835, 2010 Seismic structure of the Longmen Shan region from S wave tomography and its relationship with the Wenchuan
More informationCrustal shear velocity structure of the western US inferred from ambient seismic noise and earthquake data
JOURNAL OF GEOPHYSICAL RESEARCH, VOL.???, XXXX, DOI:10.1029/, 1 2 3 Crustal shear velocity structure of the western US inferred from ambient seismic noise and earthquake data M. P. Moschetti, 1 M. H. Ritzwoller,
More informationSeismic tomography: Art or science? Frederik J Simons Princeton University
Seismic tomography: Art or science? Frederik J Simons Princeton University What s inside the Earth? www.usgs.gov www.glencoe.com Dalton, Nature 2003 Only seismic waves have actually been there, done that
More informationSlabs, plumes and their interaction: new insights from global anisotropy tomography
Slabs, plumes and their interaction: new insights from global anisotropy tomography Ana M G Ferreira Seismological Laboratory, Department of Earth Sciences University College London, UK Sung-Joon Chang,
More informationBasic Ray Tracing. Rick Aster and Sue Bilek. October 3, 2003
Basic Ray Tracing Rick Aster and Sue Bilek October 3, 3 A key observation that we can make about a seismic signal is its arrival time. From systematic observations of arrival times, we can deduce useful
More informationCrustal Anisotropy Across Eastern Tibet and Surroundings Modeled as a Depth-Dependent Tilted Hexagonally Symmetric Medium
Crustal Anisotropy Across Eastern Tibet and Surroundings Modeled as a Depth-Dependent Tilted Hexagonally Symmetric Medium Journal: Manuscript ID Draft Manuscript Type: Research Paper Date Submitted by
More informationInterval anisotropic parameters estimation in a least squares sense Case histories from West Africa
P-263 Summary Interval anisotropic parameters estimation in a least squares sense Patrizia Cibin*, Maurizio Ferla Eni E&P Division (Milano, Italy), Emmanuel Spadavecchia - Politecnico di Milano (Milano,
More information2008 Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies ADVANCED WAVEFORM SIMULATION FOR SEISMIC MONITORING
ADVANCED WAVEFORM SIMULATION FOR SEISMIC MONITORING Donald V. Helmberger 1, Arthur J. Rodgers 2, Sidao Ni 1,3, Shengji Wei 1, and Jeroen Tromp 1 California Institute of Technology 1, Lawrence Livermore
More information28th Seismic Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
GROUND TRUTH LOCATIONS USING SYNERGY BETWEEN REMOTE SENSING AND SEISMIC METHODS: SSSC AT IMS STATIONS FOR TIBETAN PLATEAU EARTHQUAKES Gene A. Ichinose 1, Chandan K. Saikia 2*, Donald V. Helmberger 3, and
More informationTravel time tomography of the uppermost mantle beneath Europe
Chapter 3 Travel time tomography of the uppermost mantle beneath Europe We have obtained a detailed P and S model of the uppermost mantle beneath Europe using regional travel time data based on the ISC
More informationImagerie de la Terre profonde avec le bruit sismique. Michel Campillo (ISTERRE, Grenoble)
Imagerie de la Terre profonde avec le bruit sismique Michel Campillo (ISTERRE, Grenoble) Body waves in the ambient noise: microseisms (Gutenberg, Vinnik..) The origin of the noise in the period band 5-10s:
More informationCentral Coast Seismicity Locations. Jeanne Hardebeck US Geological Survey Menlo Park, CA
Central Coast Seismicity Locations Jeanne Hardebeck US Geological Survey Menlo Park, CA 11/29/2011 Seismic Network Issues for Earthquake Location: - Many stations are single-component, difficult to identify
More informationState-of-the-art in passive seismic array methodology, And useful application for mineral exploration. Nick Smith PassiveX Pty. Ltd.
State-of-the-art in passive seismic array methodology, And useful application for mineral exploration Nick Smith PassiveX Pty. Ltd. l Microtremor wavefield l Why use arrays? l Array techniques l State
More informationData Repository Item For: Kinematics and geometry of active detachment faulting beneath the TAG hydrothermal field on the Mid-Atlantic Ridge
GSA Data Repository Item: 2007183 Data Repository Item For: Kinematics and geometry of active detachment faulting beneath the TAG hydrothermal field on the Mid-Atlantic Ridge Brian J. demartin 1*, Robert
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION Seismic evidence for a global low velocity layer within the Earth s upper mantle SUPPLEMENTARY MATERIAL Benoît Tauzin 1, Eric Debayle 2 & Gérard Wittlinger 3 1 Department of Earth
More information28th Seismic Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
DETERMINATION OF CRUST AND UPPER MANTLE STRUCTURE BENEATH AFRICA USING A GLOBAL OPTIMIZATION BASED WAVEFORM MODELING TECHNIQUE Jay Pulliam, Mrinal K. Sen, and Abhijit Gangopadhyay University of Texas at
More informationStructural Cause of Missed Eruption in the Lunayyir Basaltic
GSA DATA REPOSITORY 2015140 Supplementary information for the paper Structural Cause of Missed Eruption in the Lunayyir Basaltic Field (Saudi Arabia) in 2009 Koulakov, I., El Khrepy, S., Al-Arifi, N.,
More informationStrong, Wen (Shih Chung Wen, 溫士忠 ) TEL: ext FAX:
Strong, Wen (Shih Chung Wen, 溫士忠 ) TEL: +886-5-2720411 ext. 61212 FAX: +886-6-2720807 E-mail: strong@eq.ccu.edu.tw strong6212@gmail.com [Education] Ph.D., Institute of Seismology, National Chung Cheng
More informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (2016) 204, 215 235 GJI Seismology doi: 10.1093/gji/ggv441 Seismic evidence for lithospheric modification associated with intracontinental volcanism in
More informationDRAGON ADVANCED TRAINING COURSE IN ATMOSPHERE REMOTE SENSING. Inversion basics. Erkki Kyrölä Finnish Meteorological Institute
Inversion basics y = Kx + ε x ˆ = (K T K) 1 K T y Erkki Kyrölä Finnish Meteorological Institute Day 3 Lecture 1 Retrieval techniques - Erkki Kyrölä 1 Contents 1. Introduction: Measurements, models, inversion
More informationSCEC Community Fault and Velocity Models and how they might contribute to the DCPP seismic hazard assessment. Andreas Plesch Harvard University
SCEC Community Fault and Velocity Models and how they might contribute to the DCPP seismic hazard assessment Andreas Plesch Harvard University SCEC Unified Structural Representation (USR) The SCEC USR
More informationWidespread 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 informationoverlie the seismogenic zone offshore Costa Rica, making the margin particularly well suited for combined land and ocean geophysical studies (Figure
Chapter 1 Introduction Historically, highly destructive large magnitude (M w >7.0) underthrusting earthquakes nucleate along the shallow segment of subduction zone megathrust fault, and this region of
More information29th Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
TRANSITION ZONE WAVE PROPAGATION: CHARACTERIZING TRAVEL-TIME AND AMPLITUDE INFORMATION Peter M. Shearer and Jesse F. Lawrence University of California San Diego, Institute of Geophysics and Planetary Physics
More informationThe Effect of Sedimentary Basins on Surface Waves That Pass Through Them
The Effect of Sedimentary Basins on Surface Waves That Pass Through Them Journal: Manuscript ID GJI-S--00 Manuscript Type: Research Paper Date Submitted by the Author: 0-Mar-0 Complete List of Authors:
More informationStress 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 informationGeophysical Journal International
Geophysical Journal International Geophys. J. Int. (016) 07, 101 109 Advance Access publication 016 August 5 GJI Seismology doi: 10.1093/gji/ggw30 Fast and cheap approximation of Green function uncertainty
More informationCharacterization of Induced Seismicity in a Petroleum Reservoir: A Case Study
Characterization of Induced Seismicity in a Petroleum Reservoir: A Case Study Edmond Sze, M. Nafi Toksöz, and Daniel R. Burns Earth Resources Laboratory Dept. of Earth, Atmospheric and Planetary Sciences
More informationVariation in the crustal structure across central Iceland
Geophys. J. Int. (1) 145, 246 264 Variation in the crustal structure across central Iceland Zhijun Du 1 and G. R. Foulger 2 1 Institute of Theoretical Geophysics, University of Cambridge, Downing Street,
More informationThe Math, Science and Computation of Hydraulic Fracturing
The Math, Science and Computation of Hydraulic Fracturing 03/21/2013 LLNL-PRES-XXXXXX This work was performed under the auspices of the U.S. Department of Energy by under contract DE-AC52-07NA27344. Lawrence
More informationTomography of the 2011 Iwaki earthquake (M 7.0) and Fukushima
1 2 3 Auxiliary materials for Tomography of the 2011 Iwaki earthquake (M 7.0) and Fukushima nuclear power plant area 4 5 6 7 8 9 Ping Tong 1,2, Dapeng Zhao 1 and Dinghui Yang 2 [1] {Department of Geophysics,
More information28th Seismic Research Review: Ground-Based Nuclear Explosion Monitoring Technologies
GROUND TRUTH HYPOCENTERS AND 3D CRUSTAL VELOCITY STRUCTURE IN CENTRAL ASIA FROM IN-COUNTRY NETWORKS Thomas de la Torre 1, Gaspar Monsalve-Mejia 1, Anne F. Sheehan 1, Charlotte Rowe 2, and Michael Begnaud
More informationPorosity Bayesian inference from multiple well-log data
Porosity inference from multiple well log data Porosity Bayesian inference from multiple well-log data Luiz Lucchesi Loures ABSTRACT This paper reports on an inversion procedure for porosity estimation
More informationRelevance Vector Machines
LUT February 21, 2011 Support Vector Machines Model / Regression Marginal Likelihood Regression Relevance vector machines Exercise Support Vector Machines The relevance vector machine (RVM) is a bayesian
More information