[論文レビュー] 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は、主小惑星帯の小惑星動力学をシミュレートし、カタリナ・スカイ・サーベイ（CSS）観測結果と照合することで、近地小惑星（NEO）の新しい軌道分布モデルを提示する。このモデルは、主にν₆および3:1共振に起因する小惑星源からのサイズ依存的サンプリングを明らかにし、小規模なNEO（H ≈ 25）はν₆に支配され、大規模なNEO（H ≈ 15）は3:1共振に由来する。これは、最大のヤルコフスキー効果による移動と、共振近傍での特定の小惑星の傾き（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.

研究の動機と目的

• 近地小惑星（NEO）の軌道および絶対視等級分布の物理的に整合性があり、観測的に補正されたモデルを構築し、衝突リスクの評価およびサーベイ計画の改善を図ること。
• サイズ依存的動力学、近日点付近での破壊、ヤルコフスキー効果を組み込むことで、NEOのフラックス推定値の不一致を解消すること。
• 理論的小惑星源モデル（特に3:1およびν₆共振）と観測されたNEOフラックスの整合性を検証すること。
• 将来のサーベイ（例：ベルラ・C・ルビン天文台、NEOサーベイター）のデータ統合を想定した柔軟かつアップデート可能なフレームワークを提供すること。

提案手法

• 3:1やν₆などの源共振からの主小惑星帯の小惑星10⁵個をN体数値積分によりシミュレートし、内太陽系へのNEO注入を再現する。
• マルチネストアルゴリズムを用いたベイズ推論により、シミュレートされたNEO検出結果を実際のカタリナ・スカイ・サーベイ（CSS）観測結果と照合する。
• 絶対視等級H = 15からH = 28の範囲で柔軟かつ正確に表現可能な、立方スプラインを用いてNEOの絶対視等級分布をモデル化する。
• Granvikら（2016）に基づき、近日点距離が小さい場所でのNEOのサイズ依存的破壊を組み込み、観測の不完全性を補正する。
• マルチネストによる最適化により、観測データから源貢献、ヤルコフスキー効果の速度、および傾き制約を推定する。
• 補正済みモデルからユーザー定義のNEOサンプルを生成可能な公開済みのNEOMODシミュレータの開発。

実験結果

リサーチクエスチョン

• RQ1 異なる主小惑星帯共振（例：3:1、ν₆）が、さまざまな絶対視等級においてNEO集団に及ぼす相対的貢献は何か？
• RQ2 3:1共振へのH < 18のNEOのフラックスはヤルコフスキー効果によって持続可能か？これは小惑星の傾きにどのような意味を持つのか？
• RQ3 近日点付近でのサイズ依存的破壊効果は、小規模と大規模NEOの軌道分布にどのように影響を与えるか？
• RQ4 CSS検出統計に基づくと、活動を示さないコメット（JFC）は現在のNEO集団にどの程度貢献しているか？
• RQ5 モデル予測の地球、火星、月への衝突フラックスは、以前の推定値およびクレーター年代測定データと比較してどうなるか？

主な発見

• ν₆および3:1共振は、それぞれ絶対視等級H = 15およびH = 25のNEOの約30%および80%を生じさせ、サイズ依存的源貢献が顕著であることが示された。
• H < 18のNEOが3:1共振に流入するためには、主小惑星帯の小惑星が最大のヤルコフスキー効果速度（D = 1 km、a = 2.5 auの場合で約2×10⁻⁴ au Myr⁻¹）で移動する必要があり、a < 2.5 auでは約0°、a > 2.5 auでは約180°の傾きを示すことが示唆された。
• 10 mボリデの地上衝突の平均間隔は約30年と推定され、以前の推定値の3倍以上に上昇し、H ≈ 25–28の範囲でCSSが不足している可能性を示唆している。
• ν₆および3:1源からのメテオライト落下のPM超過は、それぞれ0.47±0.02および0.50±0.05であり、モルビデリとグレッドマン（1998）の予測と整合性が高く、ウィズダム（2017, 2020）の予測よりは低いが、観測値の0.63±0.02よりは低い。
• 活動を示さないコメット（JFC）のNEO集団への貢献はα_JFC < 0.017（68.3%信用区間）に制限され、ベイズ因子はコメット寄与なしのモデルをわずかに支持した。
• 火星/月衝突フラックス比（Ma/Mo）は、H = 15では7.1、H = 25では4.7であり、それぞれR_b = 2.0および1.2を示した。これは、クレーター年代測定研究で従来用いられてきたR_b ≈ 2.6よりも顕著に低い。