Collaborative Research: Incorporating SPECFEM3D numerical seismograms in the Global CMT Project

合作研究：将 SPECFEM3D 数值地震图纳入全球 CMT 项目

• 批准号: 2218859
• 负责人: Jeroen Tromp
• 金额: $ 22.66万
• 依托单位: Princeton University
• 依托单位国家: 美国
• 项目类别: Standard Grant
• 财政年份: 2022
• 资助国家: 美国
• 起止时间: 2022-09-01 至 2024-08-31
• 项目状态: 已结题
• 来源: https://www.nsf.gov/awardsearch/showAward?AWD_ID=2218859&HistoricalAwards=false
• 关键词: Collaborative; Research; Incorporating; SPECFEM3D; numerical

This research will lay the groundwork for improvements in the scientific description of earthquakes, benefitting researchers in earthquake statistics, seismic tomography, seismic hazard, and seismic discrimination, as well as tectonic interpretation and analysis. This is the first step towards making the results of complex and expensive numerical calculations available for multiple applications and a wide range of users. This work focuses on demonstrating the feasibility of the methodology, but the database and method will be valuable for any research requiring accurate predictions of the global seismic wavefield from an arbitrary earthquake at any location without expensive computation. Two graduate students will be trained in research at the nexus of earthquake science, theoretical seismology, and computational Earth science.This two-year, focused effort will (1) develop a new method for calculating, storing, and accessing high-fidelity long-period synthetic seismograms for state-of-the-art 3D tomographic models of the Earth, and (2) incorporate these seismograms in the earthquake analysis of the Global Centroid-Moment-Tensor (CMT) Project. Currently, the CMT synthetic seismograms are calculated using modern 3D Earth models, but accuracy is limited by the validity of the path-average approximation for mode summation and surface-wave ray theory, inexact predictions of the amplitude and polarization of ground motion, and other unmodeled effects, bias retrieved earthquake parameters. The incorporation of higher-fidelity synthetic seismograms in the CMT analysis will improve the characterization of seismic sources and remove concerns about a key source of uncertainty and bias. The team will adapt the spectral-element wave-equation solver SPECFEM3D_GLOBE to generate a database of kernel seismograms on a global grid of hypocenters, for a large set of station locations, using source-receiver reciprocity to speed up the calculation. Kernel seismograms on the grid will be organized and stored in a format that facilitates rapid access to a particular source region and the stations of the Global Seismographic Network. Kernel seismograms for an arbitrary centroid location will be efficiently calculated by spatial interpolation, in a manner that matches the accuracy of the full forward calculation. The CMT code will be modified to ingest the interpolated SPECFEM3D_GLOBE seismograms and testing will allow the assessment of success of the approach and method. The Princeton numerical seismology group and Lamont earthquake-analysis group will jointly evaluate the approach and fidelity of the waveform interpolation, develop practical formats for accessing the (massive) database of global waveforms, and assess the success of these developments.This project is jointly supported by the Geophysics and Instrumentation and Facilities programs in the Division of Earth Sciences. It is also co-funded by a collaboration between the Directorate for Geosciences and Office of Advanced Cyberinfrastructure to support AI/ML and open science activities in the geosciences.This award reflects NSF's statutory mission and has been deemed worthy of support through evaluation using the Foundation's intellectual merit and broader impacts review criteria.

这项研究将为改进地震的科学描述奠定基础，使地震统计、地震层析成像、地震灾害、地震判别以及构造解释和分析方面的研究人员受益。这是使复杂且昂贵的数值计算结果可用于多种应用程序和广泛用户的第一步。这项工作的重点是证明该方法的可行性，但该数据库和方法对于任何需要准确预测任何位置的任意地震的全球地震波场而无需昂贵的计算的研究都是有价值的。两名研究生将接受地震科学、理论地震学和计算地球科学相关研究的培训。这项为期两年的重点工作将（1）开发一种用于计算、存储和访问高保真长期地震数据的新方法。用于最先进的地球 3D 断层扫描模型的周期合成地震图，以及 (2) 将这些地震图纳入全球质心矩张量 (CMT) 项目的地震分析中。目前，CMT 合成地震图是使用现代 3D 地球模型计算的，但精度受到模态求和和表面波射线理论的路径平均近似的有效性、地震动幅度和极化的不精确预测等因素的限制。未建模的影响，偏差检索地震参数。在 CMT 分析中纳入更高保真度的合成地震图将改善震源的表征，并消除对不确定性和偏差的关键来源的担忧。该团队将采用谱元波动方程求解器 SPECFEM3D_GLOBE，在全球震源网格上生成大量台站位置的核心地震图数据库，并使用源接收器互易性来加速计算。网格上的核心地震图将以一种便于快速访问特定震源区域和全球地震网络台站的格式进行组织和存储。任意质心位置的核地震图将通过空间插值有效地计算，其方式与完全正演计算的精度相匹配。 CMT 代码将被修改以获取插值的 SPECFEM3D_GLOBE 地震图，并且测试将允许评估该方法和方法的成功。普林斯顿数值地震学小组和拉蒙特地震分析小组将共同评估波形插值的方法和保真度，开发用于访问（海量）全球波形数据库的实用格式，并评估这些开发的成功。该项目得到共同支持由地球科学部的地球物理学和仪器仪表和设施项目负责。它还由地球科学理事会和先进网络基础设施办公室之间的合作共同资助，以支持地球科学领域的人工智能/机器学习和开放科学活动。该奖项反映了 NSF 的法定使命，并通过使用基金会的评估进行评估，被认为值得支持。智力价值和更广泛的影响审查标准。