Collaborative Research: Wave Equation Tomography And Data Assimilation: A New Approach To Estimating P And S Speed Variations In Earth's Lower Mantle

主要研究人员寻求合作研究的资金，旨在通过明确整合地震数据分析和解释的观测，理论和计算方面来开发一种新的层析成像模型。这一努力扩展了他们以前在层析成像方面的工作，涉及新数据集的测量以及波传播理论和多重网格技术的发展，这种技术允许联合解释对地球结构具有不同敏感性的数据。在我们看来，对这些组分进行定量积分，是对弹性性质、温度和组成的空间变化进行精确映射的可行且必要的一步。在这里，PI的将集中在VP和VS的相对变化，在底部~1000公里的地球下地幔，其中，他们认为，包含关键线索tunderstanding地幔对流和地球的热化学演化的地质time.Intellectual优点：在过去的几十年中，全球层析成像产生了壮观的图像，例如，地幔流动轨迹和结构复杂的地幔底部附近。然而，不均匀的数据覆盖面和异构的数据质量呈现非唯一的，模糊的图像，在可靠性方面有很大的空间差异。正则化和使用不准确的波传播理论可能会产生不正确的估计弹性参数，即使在采样似乎足够，和波速变化的幅度通常约束不良。此外，基于不同数据集的结果通常在重要方面不一致，并且对对地球结构具有不同敏感性的数据的正确联合解释（例如，体波和表面波，在不同频率下测量的P波或S波）仍然是一个重大挑战。"红蓝"图像的近似性质阻碍了定量解释和与其他地球物理约束的整合，并使层析成像无法充分发挥其作为地球深部定量探测的潜力。这是PI想要改变的。为了更好地估计参数，他们需要利用丰富的宽带波形，他们需要更强大的理论框架，整合和联合解释不同的数据集。我们对全球层析成像的多分辨率数据融合方法的最终目标是通过改进，在一系列长度尺度上，从各种地震学数据中产生更好的地球深部内部三维模型。（并明确联系）成像的三个基本方面：数据质量和覆盖范围：使用自动化程序和多分辨率概念（如时间频率小波）它们将通过从大量可用的波形中提取相速度和到达时间信息来增强空间和频谱数据的覆盖范围通过国际数据中心。波传播理论：认识到需要考虑所考虑的数据的不同采样特性（并从中受益），并受到最近在理解有限频率效应方面的进展的启发，他们将计算新测量数据的反投影的精确灵敏度核。参数化和正则化：为了保持和利用三维灵敏度核的局部化特性，我们将采用自适应多重网格参数化和正则化技术进行联合反演.本文的研究重点是：（i）测量三维P型和S型体波走时，（ii）反演. lnVS/. lnVP的三维变化（或相关参数）在地球的地幔，和（iii）精炼-或驳斥-现有的观点组成的不均匀性在最低的地幔。他们可以在观测地震学和层析成像（货车der Hessel，麻省理工学院）以及波传播和反演理论（De Hoop，CSM）方面积累经验，在自动化数据处理方面，他们将与Ritsema（IPGP，法国）合作，并涉及一名博士后助理（为此寻求一些资金匹配）。更广泛的影响：沿着矿物物理学数据，需要精确估计弹性参数，以限制成分和温度的空间变化，从而限制地幔动力学和矿物学模型。此外，这里开发的数据融合的概念和工具准备的USARAY数据的大数据集的处理和解释。拟议的工作构成了麻省理工学院博士项目的第一部分，但麻省理工学院和CSM的学生将参与研究的各个方面，无论是作为本科研究机会（UROP）还是满足一般考试要求。

The Principal Investigators seek funding for collaborative research that aims to develop a new class of tomographic model by explicit integration of observational, theoretical, and computational aspects of seismic data analysis and interpretation. This effort extends their previous work in tomography and involves the measurement of new data sets and the development of wave propagation theory and multi-grid technology that allowsjoint interpretation of data with different sensitivities to Earth's structure. Quantitative integration of thesecomponents is, in our view, a viable and essential step toward the accurate mapping of spatial variations in elastic properties, temperature, and composition. Here the PI's will focus on relative variations in VP and VS in the bottom ~1000 km of Earth's lower mantle, which, they believe, contains critical clues tounderstanding mantle convection and Earth's thermo-chemical evolution over geological time.Intellectual Merit: Over the past decades global tomography has produced spectacular images of, for instance, mantle flow trajectories and structural complexity near the base of the mantle. However, uneven data coverage and heterogeneous data quality render non-unique, fuzzy images, with substantial spatial variations in reliability. Regularization and the use of inaccurate wave propagation theory probably produce incorrect estimates of elastic parameters even where sampling seems adequate, and the magnitude of wavespeed variations is usually poorly constrained. Moreover, results based on different data sets often disagree in important aspects, and correct joint interpretation of data with different sensitivities to Earth's structure (e.g., body- and surface waves, P or S waves measured at different frequencies) remains a major challenge. The approximate nature of the "red and blue" images impedes quantitative interpretation and integration with other geophysical constraints and keeps tomography from reaching its full potential as a quantitative probe of Earth's deep interior. This the PI's seek to change. For better parameter estimation they need to exploit the richness of broad-band waveforms and they need more powerful theoretical frameworks for integration and joint interpretation of diverse data sets. The ultimate objective of our approach toward multi-resolution data fusion for global tomography is to produce better 3-D models of Earth's deep interior - on a range of length scales and from a variety of seismological data - by improving (and explicitly linking) three essential aspects of imaging: Data quality and coverage: using automated procedures and multi-resolution concepts (such as time frequency wavelets) they will enhance spatial and spectral data coverage by extracting phase velocity and arrival time information from the vast number of waveforms available through international data centers. Wave propagation theory: recognizing the need to account for (and benefit from) the different sampling properties of the data considered, and inspired by recent advances in understanding finite frequency effects, they will compute accurate sensitivity kernels for the back-projection of the newly measured data. Parameterization and regularization: to preserve and exploit the localization properties of 3-D sensitivity kernels we will use adaptive multi-grid parameterization and regularization techniques for joint inversion.The research proposed here focuses on (i) measuring teleseismic P and S type body-wave travel times, (ii) inversion for 3-D variations in .lnVS/.lnVP (or related parameters) in Earth's mantle, and (iii) refining - or refuting - existing views on compositional heterogeneity in the lowermost mantle. They can build on experience in observational seismology and tomography (Van der Hilst, MIT) and wave propagation and inversion theory (De Hoop, CSM), and for the automated data processing they will collaborate with Ritsema (IPGP, France) and involve a postdoctoral associate (for which some fund matching is sought). Broader Impact: Along with mineral physics data, accurate estimates of elastic parameters are needed to constrain spatial variations in compositon and temperature and, thus, models of mantle dynamics and mineralogy. Furthermore, the concept of and tools for data fusion developed here prepare for the handling and interpretation of large data sets of USARRAY data. The proposed work constitutes the first part of a PhD project at MIT, but students at MIT and CSM will be involved in aspects of the research,either as a Undergraduate Research OPportunity (UROP) or in fulfillment of General Exam requirements.