SItomo - 用于东阿尔卑斯山分裂强度断层摄影和应用的工具箱

IF 2.1 3区 地球科学 Q2 GEOCHEMISTRY & GEOPHYSICS
Frederik Link, Maureen D. Long
{"title":"SItomo - 用于东阿尔卑斯山分裂强度断层摄影和应用的工具箱","authors":"Frederik Link,&nbsp;Maureen D. Long","doi":"10.1016/j.jog.2024.102018","DOIUrl":null,"url":null,"abstract":"<div><p>The tomographic inversion of shear wave splitting data for upper mantle anisotropy has been a longstanding challenge. This is due to the ray-based approximation of classical approaches and the near-vertical incidence of the core-mantle converted phases such as SKS that are often used. Recent developments include the calculation of finite-frequency sensitivity kernels for SKS splitting intensity observations, which allows us to accurately take into account the sensitivity to anisotropic structure with depth. A requirement of this tomographic technique is a dense station spacing, which results in overlapping sensitivity kernels at depth and allows for the localization of anisotropic structure. This is satisfied by a growing number of temporary seismic deployments, which motivates the desire to image anisotropic complexities with depth. Here, we introduce and make available a toolbox for the MATLAB environment that facilitates the application of finite-frequency splitting intensity tomography to dense seismic arrays. Our implementation includes several key features, including: 1) A forward calculation of splitting intensities and sensitivity kernels for a complex anisotropic model space. 2) Consideration of the dominant period of the wave, allowing for multiple-frequency analysis, as well as the incoming wave’s non-vertical incidence. 3) The inversion can be based on a classical gradient descent, on a form of the conjugate gradient method known as the BFGS algorithm, or on a gradient-informed stochastic reversible jump algorithm, allowing for a data-driven parametrization of the model space. 4) Importing splitting intensity measurements from waveforms processed in SplitRacer allows for fast pre-processing of large data sets due to its fully automatic design. To illustrate our method, we present both synthetic tests and an application to real data. We apply our inversion procedure to data from the Swath-D network, which densely covers the transition of the Central to the Eastern Alps. Previous studies showed evidence for an abrupt lateral change of layered seismic anisotropy that had been attributed to an opening for channeled asthenospheric flow. Using an SKS splitting intensity tomography approach, we can confirm previous inferences while providing additional constraints on the distribution of anisotropy laterally and with depth.</p></div>","PeriodicalId":54823,"journal":{"name":"Journal of Geodynamics","volume":null,"pages":null},"PeriodicalIF":2.1000,"publicationDate":"2024-01-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"SItomo – A toolbox for splitting intensity tomography and application in the Eastern Alps\",\"authors\":\"Frederik Link,&nbsp;Maureen D. Long\",\"doi\":\"10.1016/j.jog.2024.102018\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"<div><p>The tomographic inversion of shear wave splitting data for upper mantle anisotropy has been a longstanding challenge. This is due to the ray-based approximation of classical approaches and the near-vertical incidence of the core-mantle converted phases such as SKS that are often used. Recent developments include the calculation of finite-frequency sensitivity kernels for SKS splitting intensity observations, which allows us to accurately take into account the sensitivity to anisotropic structure with depth. A requirement of this tomographic technique is a dense station spacing, which results in overlapping sensitivity kernels at depth and allows for the localization of anisotropic structure. This is satisfied by a growing number of temporary seismic deployments, which motivates the desire to image anisotropic complexities with depth. Here, we introduce and make available a toolbox for the MATLAB environment that facilitates the application of finite-frequency splitting intensity tomography to dense seismic arrays. Our implementation includes several key features, including: 1) A forward calculation of splitting intensities and sensitivity kernels for a complex anisotropic model space. 2) Consideration of the dominant period of the wave, allowing for multiple-frequency analysis, as well as the incoming wave’s non-vertical incidence. 3) The inversion can be based on a classical gradient descent, on a form of the conjugate gradient method known as the BFGS algorithm, or on a gradient-informed stochastic reversible jump algorithm, allowing for a data-driven parametrization of the model space. 4) Importing splitting intensity measurements from waveforms processed in SplitRacer allows for fast pre-processing of large data sets due to its fully automatic design. To illustrate our method, we present both synthetic tests and an application to real data. We apply our inversion procedure to data from the Swath-D network, which densely covers the transition of the Central to the Eastern Alps. Previous studies showed evidence for an abrupt lateral change of layered seismic anisotropy that had been attributed to an opening for channeled asthenospheric flow. Using an SKS splitting intensity tomography approach, we can confirm previous inferences while providing additional constraints on the distribution of anisotropy laterally and with depth.</p></div>\",\"PeriodicalId\":54823,\"journal\":{\"name\":\"Journal of Geodynamics\",\"volume\":null,\"pages\":null},\"PeriodicalIF\":2.1000,\"publicationDate\":\"2024-01-19\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"Journal of Geodynamics\",\"FirstCategoryId\":\"89\",\"ListUrlMain\":\"https://www.sciencedirect.com/science/article/pii/S0264370724000012\",\"RegionNum\":3,\"RegionCategory\":\"地球科学\",\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q2\",\"JCRName\":\"GEOCHEMISTRY & GEOPHYSICS\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Geodynamics","FirstCategoryId":"89","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0264370724000012","RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q2","JCRName":"GEOCHEMISTRY & GEOPHYSICS","Score":null,"Total":0}
引用次数: 0

摘要

上地幔各向异性剪切波分裂数据的层析反演是一项长期挑战。这是由于经典方法基于射线的近似,以及经常使用的核幔转换相(如 SKS)的近垂直入射。最近的发展包括计算 SKS 分裂强度观测的有限频率灵敏度核,这使我们能够准确地考虑到各向异性结构随深度变化的灵敏度。这种层析技术的一个要求是密集的站点间距,这将导致深度上的灵敏度核重叠,并允许对各向异性结构进行定位。越来越多的临时地震部署满足了这一要求,从而激发了对各向异性复杂性进行深度成像的愿望。在此,我们介绍并提供了一个 MATLAB 环境工具箱,该工具箱有助于将有限频率分裂强度层析成像技术应用于密集地震阵列。我们的实现包括以下几个主要特点1) 复杂各向异性模型空间的分裂强度和灵敏度核的前向计算。2) 考虑波的主导周期,允许多频率分析,以及入射波的非垂直入射。3) 反演可以基于经典的梯度下降法、共轭梯度法的一种形式(即 BFGS 算法),也可以基于梯度信息随机可逆跃迁算法,从而对模型空间进行数据驱动参数化。4)从 SplitRacer 处理过的波形中导入分裂强度测量值,由于其全自动设计,可对大型数据集进行快速预处理。我们将反演程序应用于 Swath-D 网络的数据,该网络密集覆盖了中阿尔卑斯山向东阿尔卑斯山的过渡区域。之前的研究显示,有证据表明层状地震各向异性发生了突然的横向变化,而这种变化被归因于星体流的通道开口。利用 SKS 分裂强度层析成像方法,我们可以证实之前的推论,同时为各向异性的横向分布和深度分布提供额外的约束条件。
本文章由计算机程序翻译,如有差异,请以英文原文为准。
SItomo – A toolbox for splitting intensity tomography and application in the Eastern Alps

The tomographic inversion of shear wave splitting data for upper mantle anisotropy has been a longstanding challenge. This is due to the ray-based approximation of classical approaches and the near-vertical incidence of the core-mantle converted phases such as SKS that are often used. Recent developments include the calculation of finite-frequency sensitivity kernels for SKS splitting intensity observations, which allows us to accurately take into account the sensitivity to anisotropic structure with depth. A requirement of this tomographic technique is a dense station spacing, which results in overlapping sensitivity kernels at depth and allows for the localization of anisotropic structure. This is satisfied by a growing number of temporary seismic deployments, which motivates the desire to image anisotropic complexities with depth. Here, we introduce and make available a toolbox for the MATLAB environment that facilitates the application of finite-frequency splitting intensity tomography to dense seismic arrays. Our implementation includes several key features, including: 1) A forward calculation of splitting intensities and sensitivity kernels for a complex anisotropic model space. 2) Consideration of the dominant period of the wave, allowing for multiple-frequency analysis, as well as the incoming wave’s non-vertical incidence. 3) The inversion can be based on a classical gradient descent, on a form of the conjugate gradient method known as the BFGS algorithm, or on a gradient-informed stochastic reversible jump algorithm, allowing for a data-driven parametrization of the model space. 4) Importing splitting intensity measurements from waveforms processed in SplitRacer allows for fast pre-processing of large data sets due to its fully automatic design. To illustrate our method, we present both synthetic tests and an application to real data. We apply our inversion procedure to data from the Swath-D network, which densely covers the transition of the Central to the Eastern Alps. Previous studies showed evidence for an abrupt lateral change of layered seismic anisotropy that had been attributed to an opening for channeled asthenospheric flow. Using an SKS splitting intensity tomography approach, we can confirm previous inferences while providing additional constraints on the distribution of anisotropy laterally and with depth.

求助全文
通过发布文献求助,成功后即可免费获取论文全文。 去求助
来源期刊
Journal of Geodynamics
Journal of Geodynamics 地学-地球化学与地球物理
CiteScore
4.60
自引率
0.00%
发文量
21
审稿时长
6-12 weeks
期刊介绍: The Journal of Geodynamics is an international and interdisciplinary forum for the publication of results and discussions of solid earth research in geodetic, geophysical, geological and geochemical geodynamics, with special emphasis on the large scale processes involved.
×
引用
GB/T 7714-2015
复制
MLA
复制
APA
复制
导出至
BibTeX EndNote RefMan NoteFirst NoteExpress
×
提示
您的信息不完整,为了账户安全,请先补充。
现在去补充
×
提示
您因"违规操作"
具体请查看互助需知
我知道了
×
提示
确定
请完成安全验证×
copy
已复制链接
快去分享给好友吧!
我知道了
右上角分享
点击右上角分享
0
联系我们:info@booksci.cn Book学术提供免费学术资源搜索服务,方便国内外学者检索中英文文献。致力于提供最便捷和优质的服务体验。 Copyright © 2023 布克学术 All rights reserved.
京ICP备2023020795号-1
ghs 京公网安备 11010802042870号
Book学术文献互助
Book学术文献互助群
群 号:481959085
Book学术官方微信