Integrating the LPO Constraint into 3D Subduction Dynamics Simulations

俯冲是一个构造板块下沉到另一个构造板块之下进入地幔的过程。俯冲带是地球上构造和地震最活跃的地区。它们有巨大的破坏性地震和丰富的火山活动。这种地质活动是由下沉的构造板块驱动的，它导致周围地幔的变形（应力积累）和熔融。这些地下过程是利用地球物理和地球化学方法间接观测到的。在这些观测中，只有地震各向异性可以直接与俯冲板块下沉引起的变形的空间模式联系起来。为了将地震各向异性的观测结果与地幔变形联系起来，科学家们需要计算变形岩石中形成的一种结构，这种结构被称为晶格择优取向（LPO）。拟议的研究将模拟俯冲，以预测LPO结构。这是首次将LPO计算完全纳入大规模三维模拟。通过分析这些结果，研究人员将从浅地幔三维变形的角度对地震各向异性观测提供更好的解释。更好地了解地幔变形（应力）对于解决俯冲带地震和火山活动造成的自然灾害非常重要。为该项目开发的所有计算模块将公开提供，以推动该项目以外的类似研究。该项目将培养一名博士生，学习计算建模、复杂系统分析以及如何将科学成果传达给广大受众。该项目还将为本科生提供研究机会，培养他们在批判性思维，项目开发，领导能力和科学方法方面的技能。拟议的研究将使用地幔模拟代码ASPECT创建，运行和分析约16个，3-D，随时间变化的俯冲模拟，使用示踪粒子和运动学（D-Rex）LPO计算，对整个模型域中LPO结构的演变进行了完全集成的计算。将对这些模型进行分析，以确定LPO结构中与定义俯冲系统的不同参数（物理性质、几何形状）相关的特征。他们也将进行分析，以确定模式的LPO结构有关的俯冲阶段和状态的变形板。在地面观测到的地震各向异性的预测将由一组设计成与卡斯卡迪亚、S。美洲和中美洲俯冲带。然后，研究人员将比较预测和观测到的地震各向异性，以确定观测结果与地幔流动和板块变形的关系。研究生研究员将主要负责创建，运行和处理模拟，以及俯冲系统参数和俯冲阶段/状态的结果分析。PI将主要负责编写后处理代码，从ASPECT获取LPO结构，计算弹性张量和由此产生的地震各向异性，并将这些预测与三个地点的观测结果进行比较。他们预计至少出版六份载有该项目主要调查结果的手稿。

Subduction is the process by which one tectonic plate sinks beneath another tectonic plate into the Earth's mantle. Subduction zones are some of the most tectonically and seismically active areas on the planet. They have large destructive earthquakes and abundant volcanic activity. This geologic activity is driven by the sinking tectonic plate, which causes deformation (build up of stresses) and melting in the surrounding mantle. These subsurface processes are observed indirectly using geophysical and geochemical methods. Of these observations, only seismic anisotropy can be directly linked to the spatial pattern of deformation caused by sinking of the subducted plate. To make the link from observations of seismic anisotropy to deformation in the mantle scientists need to calculate a fabric that develops in the deformed rocks called lattice-preferred orientation (LPO). The proposed research will simulate subduction to predict the LPO structure. This is the first time LPO calculations will be fully incorporated into large-scale 3-D simulations. By analyzing these results the investigators will provide better interpretations of seismic anisotropy observation in terms of the 3-D deformation in the shallow mantle. A better understanding of the mantle deformation (stresses) is important for addressing natural hazards due to earthquakes and volcanic activity in subduction zones. All of the computational modules developed for the project will be made publicly available to advance similar studies beyond this project. The project will train a PhD student in computational modeling, analysis of complex systems and how to communicate scientific results to a broad audience. The project will also provide research opportunities for undergraduates to develop their skills in critical-thinking, project development, leadership, and the scientific method.The proposed research will create, run and analyze approximately 16, 3-D, time-dependent simulations of subduction using the mantle simulation code ASPECT, with fully-integrated calculations of the evolution of LPO structure throughout the model domain using tracer particles and a kinematic (D-Rex) LPO calculation. The models will be analyzed to identify features in the LPO structure that are related to different parameters (physical properties, geometry) defining the subduction system. They will also be analyzed to identify patterns in the LPO structure related to the phase of subduction and state of deformation in the slab. Predictions of seismic anisotropy observed at the surface will be made for a subset of models designed to match the general characteristics of the Cascadia, S. America and Middle America subduction zones. The investigators will then compare the predicted and observed seismic anisotropy to determine how the observations are related to flow in the mantle and deformation of the slab. The graduate student researcher, will be primarily responsible for creating, running, and processing simulations, as well as the analysis of the results in terms of subduction system parameters and phase/state of subduction. The PI will be primarily responsible for writing the post-processing codes that takes the LPO structure from ASPECT, calculates the elastic tensor and resulting seismic anisotropy, and for the comparison of these predictions to observations from three locations. They anticipate publication of at least six manuscripts with the major findings of the project.