准确定位火星地震（火震）能够为理解火星内部地震活动性、地下结构、物质组分以及深部动力学过程提供关键信息。通常，描述火震的位置需要两个基本参数：震中距和反方位角。在美国InSight（洞察号）任务的单台定位研究中，地震波走时分析和波形极化分析分别用于确定震中距和反方位角。然而，相比于震中距，反方位角的有效测定更具挑战性，因为InSight记录的地震波振幅比走时更容易受到噪声干扰。前人的研究已使用多种方法来确定火震的反方位角，但对于每一个火震，通常需要人为选择最优的地震记录时间窗和频率范围来进行分析。由于对这些参数的选择具有较强的主观性，分析结果往往存在较大的不确定度，因此目前仅获得了少量低频火震位置信息，而且不同研究对部分火震的定位结果也存在明显差异，从而制约了对火震活动性的深入认识和对火星内部结构的有效约束。 

　　针对火震反方位角定位难题，中国科学院地质与地球物理研究所王旭博士后、陈凌研究员和王新副研究员经过反复试验测试，提出了一种基于变参数极化分析的单台站事件反方位角有效测定方法，用以确定InSight低频火震的反方位角。该方法采用滑动时间和频率窗来处理波形数据，基于主成分分析分别获取每个参数组合下的火震事件波形极化特征和噪声极化特征，并通过统计分析和与噪声对比确定事件反方位角的最优估计和不确定度。为了验证方法的有效性，研究团队将其应用于两个已知真实位置（由轨道图像确认）的火星陨石撞击事件S1094b和S1000a波形数据。采用新方法约束的位置比已报道的InSight团队MQS（MarsQuake Service）定位结果更接近真实的陨石撞击点，两个事件反方位角测定值与真实值的差距分别缩小了约5度和14度（图1）。这表明新方法能够有效压制参数主观选择所带来的影响，在反方位角测定精度方面具有显著优势。 

图1 两个火星陨石撞击事件的真实位置（黑色五角星）、本研究测定的撞击位置（黄色圆圈）和MQS发布的撞击位置（白色圆圈）的对比。（a）撞击事件在火星表面的位置及其反方位角偏差；（b）S1094b事件反方位角统计测定结果；（c）S1000a事件反方位角统计测定结果

　　研究团队采用新方法进一步深入分析了 44个低频火震，其中39个获得了可靠的反方位角测定结果。结合MQS发布的震中距信息，对这些火震进行了重新定位。新的定位结果显示（图2a），低频火震主要分布在InSight着陆点以东约1500公里处的科柏洛斯槽沟（Cerberus Fossae）区域，同时在北部低地、南部高地以及沿火星二分边界附近也存在火震分布。这一新的研究结果提供了目前为止最为完整的低频火震分布图像，揭示科柏洛斯槽沟区域的火震更靠近大约5万年前仍有火山活动的断裂尾端（图2b），而在其他区域火震的分布也更为广泛，特别是首次在南部赫斯珀利亚高原（Hesperia Planum）发现了4个火震事件，表明火星的构造活动可能比以往所认为的更为活跃。然而，要深入理解不同区域低频火震的成因及其反映的火星构造特征，还需要开展更多的火震定位和震源机制研究。 

图2 重新定位的低频火星震分布

　　研究成果发表于国际学术期刊GRL（王旭，陈凌^*，王新. Renewed epicentral distribution of Low Frequency marsquakes by varying-parameter polarization analysis of InSight data[J]. Geophysical Research Letters, 50: e2023GL103896. DOI: 10.1029/2023GL103896）。研究得到国家重点研发计划（2022YFF0503203）、国家自然科学基金（42004041, 42288201）和中国博士后科学基金（2021T140659）共同资助。