[论文解读] NEOMOD: A New Orbital Distribution Model for Near Earth Objects

David Nesvorný, Rogerio Deienno+11 more|arXiv (Cornell University)|Jun 15, 2023

Astro and Planetary Science被引用 4

一句话总结

NEOMOD 提出了一种新的近地小行星（NEOs）轨道分布模型，通过模拟主带小行星动力学并校准卡特琳娜天空调查（CSS）观测结果实现。该模型揭示了主带源区（特别是 ν₆ 和 3:1 共振）的尺寸依赖性采样：小尺寸近地小行星（H ≈ 25）主要源自 ν₆ 共振，而大尺寸近地小行星（H ≈ 15）则来自 3:1 共振，其持续通量需要最大 Yarkovsky 驱动及在共振附近特定的自转轴倾角（0° 或 180°）。

ABSTRACT

Near Earth Objects (NEOs) are a transient population of small bodies with orbits near or in the terrestrial planet region. They represent a mid-stage in the dynamical cycle of asteroids and comets, which starts with their removal from the respective source regions -- the main belt and trans-Neptunian scattered disk -- and ends as bodies impact planets, disintegrate near the Sun, or are ejected from the Solar System. Here we develop a new orbital model of NEOs by numerically integrating asteroid orbits from main belt sources and calibrating the results on observations of the Catalina Sky Survey. The results imply a size-dependent sampling of the main belt with the $ν_6$ and 3:1 resonances producing $\simeq 30$\% of NEOs with absolute magnitudes $H = 15$ and $\simeq 80$\% of NEOs with $H = 25$. Hence, the large and small NEOs have different orbital distributions. The inferred flux of $H<18$ bodies into the 3:1 resonance can be sustained only if the main-belt asteroids near the resonance drift toward the resonance at the maximal Yarkovsky rate ($\simeq 2 imes 10^{-4}$ au Myr$^{-1}$ for diameter $D=1$ km and semimajor axis $a=2.5$~au). This implies obliquities $θ\simeq 0^\circ$ for $a<2.5$~au and $θ\simeq 180^\circ$ for $a>2.5$~au, both in the immediate neighborhood of the resonance (the same applies to other resonances as well). We confirm the size-dependent disruption of asteroids near the Sun found in previous studies. An interested researcher can use the publicly available NEOMOD Simulator to generate user-defined samples of NEOs from our model.

研究动机与目标

• 开发一个物理上一致且经过观测校准的近地小行星（NEOs）轨道与绝对星等分布模型，以改进撞击风险评估与巡天规划。
• 通过整合尺寸依赖的动力学、近日点附近解体效应及 Yarkovsky 驱动影响，解决近地小行星通量估计中的差异。
• 检验观测到的近地小行星通量与理论小行星源模型（特别是 3:1 和 ν₆ 共振）的一致性。
• 为未来与即将开展的巡天任务（如薇拉·C·鲁宾天文台、近地天体侦察者）的数据集成提供一个灵活且可更新的框架。

提出的方法

• 对来自源共振区（如 3:1、ν₆）的 10⁵ 颗主带小行星进行数值 N 体积分，模拟其向内太阳系的近地小行星注入过程。
• 利用 MultiNest 算法进行贝叶斯推断，将模拟的近地小行星探测结果与实际卡特琳娜天空调查（CSS）观测进行校准。
• 采用三次样条插值方法建模近地小行星的绝对星等分布，实现在 H = 15 至 H = 28 范围内的灵活且精确表示。
• 基于 Granvik 等人（2016）的研究，引入在近日点附近的小行星尺寸依赖性解体效应，以反映观测的不完整性。
• 通过 MultiNest 优化模型，从观测数据中推断源区贡献、Yarkovsky 驱动速率及倾角约束。
• 开发公开可用的 NEOMOD 模拟器，可基于校准模型生成用户自定义的近地小行星样本。

实验结果

研究问题

• RQ1 不同主带共振（如 3:1、ν₆）对不同绝对星等下近地小行星种群的相对贡献如何？
• RQ2 H < 18 的近地小行星通量能否由 Yarkovsky 驱动维持进入 3:1 共振？这对小行星倾角意味着什么？
• RQ3 近日点附近尺寸依赖的解体效应如何影响小尺寸与大尺寸近地小行星的轨道分布？
• RQ4 基于 CSS 检测统计数据，惰性彗星（JFCs）对当前近地小行星种群的贡献是多少？
• RQ5 模型预测的地球、火星和月球上的撞击通量与以往估计值及陨石坑年代学数据相比如何？

主要发现

• ν₆ 和 3:1 共振分别产生约 30% 和 80% 的绝对星等 H = 15 的近地小行星，以及 H = 25 的近地小行星，表明存在显著的尺寸依赖性源区贡献。
• H < 18 的近地小行星通量进入 3:1 共振要求主带小行星在共振附近以最大 Yarkovsky 速率漂移（D = 1 km、a = 2.5 au 时约为 ~2×10⁻⁴ au Myr⁻¹），这意味着在 a < 2.5 au 区域倾角约为 ~0°，在 a > 2.5 au 区域约为 ~180°。
• 估算得出，10 米火流星撞击地球的平均时间间隔约为 30 年，比以往估计值长三倍，提示 CSS 对 H ≈ 25–28 的天体可能存在漏报。
• 来自 ν₆ 和 3:1 源区的陨石雨 PM 超额值分别为 0.47±0.02 和 0.50±0.05，与 Morbidelli 与 Gladman（1998）的预测更一致，但低于观测值 0.63±0.02。
• 惰性彗星（JFCs）对近地小行星种群的贡献被限制在 α_JFC < 0.017（68.3% 可信区间），贝叶斯因子略倾向于无彗星贡献的模型。
• 对于 H = 15 和 H = 25，火星与月球的撞击通量比（Ma/Mo）分别为 7.1 和 4.7，对应 R_b = 2.0 和 1.2，显著低于陨石坑年代学研究中曾使用的 R_b ≈ 2.6。