密度、强度、反照率等物理特性参数是评估近地小行星撞击风险与危害、制定小行星防御处置策略的基本依据。然而，小行星密度、强度等特性参数一般仅能通过飞行器抵近探测获取，迄今为止人类仅掌握极少数实施过飞行探测任务的近地小行星的这一类物理特性参数。

　　近日，中国科学院国家空间科学中心复杂航天系统电子信息技术重点实验室小行星防御课题组博士研究生耿淑娟与研究员周炳红、李明涛，提出了利用小行星进入大气层的火流星观测数据约束其物理性质和轨道的方法，同时，利用该方法对2022EB5近地小行星的物理性质进行约束，推断其或是一颗低密度、低强度、低反照率的碳质小行星。该研究为未来通过火流星观测数据获取小行星密度、强度、反照率等物理特性的分布特征提供了新途径。

　　小行星2022EB5是人类第5颗提前预警的撞击地球的近地小行星，于格林尼治时间2022年3月11日由匈牙利天文学家Krisztián Sárneczky首次发现，并在约两小时后进入地球大气层，于挪威海上空发生解体。NASA近地小天体研究中心（CNEOS）利用美国卫星传感器网络（USGS）记录了该小行星在地球大气内飞行过程的光变曲线以及峰值亮度时刻的速度、位置，并根据光能推断该小行星的撞击能量（约4千吨TNT）。美国宇航局和欧空局等通过进入大气层前的光学观测资料确定了该小行星的轨道，利用小行星的平均反照率估计了小行星的尺寸，而该小行星其他物理性质尚不清楚。

　　对于已知峰值亮度时刻状态矢量的小行星，传统的推断物理性质的方法是将其峰值亮度时刻的速度视为进入大气时刻的速度，将峰值亮度时刻的动压视为其解体强度的最大值，从而根据估算的撞击能量进一步计算小行星的质量和轨道，推测小行星的材质。对于有光学观测信息的小行星，通常假设平均反照率，根据绝对星等、直径、反照率间的经验关系，对小行星的尺寸进行计算。然而，这种传统的计算方式相当于跳过小行星进入大气过程的计算，忽略了大气层对小行星轨迹的影响。同时，平均反照率估算的小行星尺寸往往与实际情况存在较大偏差。

　　然而，对小行星的大气进入过程建模，需要输入小行星密度、结构强度、烧蚀热等物理特性参数。这些物理特性参数恰是需要估计的未知量。如何在未知物理特性参数的条件下对小行星的大气进入过程进行建模并同步估计未知物理特性参数是难题之一。

　　该研究将计算分为四组，分别假设小行星为C类、S类、M类及彗星，并根据经验参数在每组中假设小行星的密度、烧蚀热为对应材质类型的经典值。研究将其他未知的物理量视为待优化变量，包括小行星峰值亮度时刻的半径、半径变化率、质量以及小行星的解体强度。以煎饼模型的解体条件作为优化的约束条件，研究以此计算小行星的尺寸、质量、解体强度。结合小行星各组是否有解以及解体强度的范围，研究对小行星的材质进行初步判断。例如，若仅有S类的组别有解，则该小行星为S类小行星；若两组以上有解，则根据强度对小行星类型进行判断。通常而言，彗星的强度<0.1MPa，M类小行星强度>10MPa，但C类与S类小行星的强度往往不能严格区分，需要其他辅助信息进一步约束。

　　该课题组利用Chelyabinsk流星体、Kosice流星体对此方法进行验证。结果显示，使用该方法计算得到的小行星物理性质以及轨道根数，均较传统方法计算得到的的数值更符合观测结果。

　　研究显示，对于2022EB5小行星，使用传统方法计算，其反照率将被假设为0.147，得到的小行星直径约1.6米；其解体强度约3.5MPa，若仅根据强度推断，2022 EB5可能被确定为S类小行星。利用这一新方法则能够对2022EB5的物理性质进行进一步约束，如表1所示。其中，Ryugu/Bennu组是日本和美国小行星取样返回任务的目标小行星，密度约为1200kg/m³，用于模拟低密度C类小行星的撞击场景。　　　　　　

2022 EB5能量沉积曲线