我们通常将距离地表90km以上的高空大气称为高层大气,在这个高度以上,大气密度随着高度的增加急剧下降,到200km高度时,大气密度已经约为地面的五亿分之一!
我们通常将距离地表90km以上的高空大气称为高层大气,在这个高度以上,大气密度随着高度的增加急剧下降,到200km高度时,大气密度已经约为地面的五亿分之一!尽管那里的大气密度已经如此稀薄,但卫星高速地绕地球运行(约7.9km/s),所受到的大气阻力仍然不可小视。就像当我们在公路上骑行,经常因来自空气的阻力感到吃力一样,卫星也会因为大气阻力的影响而慢慢感到“劳累”从而发生轨道衰减。现有的科学技术已经可以非常准确地模拟出除大气阻力之外卫星所受的各种扰动力,包括地球非球形扰动、日月扰动和辐射光压扰动等,因此大气阻力产生的扰动已经成为卫星定轨和预报的最大误差来源。为了计算出卫星受到的大气阻力,需要尽可能准确地知道那里的大气密度到底是多少。那么如何测量出那里的大气密度呢?下面本文列举出目前几种常见的测量方法。
星载加速仪数据具有精度高、分辨率高的特点。通过表1可以看出CHAMP加速仪精度达10-9m/s2,GRACE卫星精度更是达到10-10m/s2,换算成大气密度,精度可达10-15~10-14kg/m3。而在CHAMP和GRACE的高度,大气密度量级约为10-13~10-11kg/m3。因此通过加速仪数据可以较为准确地测量出卫星轨道处的大气密度。
表1 加速仪基本信息
从表1中我们能够发现目前的高精度加速仪大部分都是法国生产的,遗憾的是,我国目前尚不具备生产如此高精度加速仪的能力。我国航天工程中的加速度计多用于载体的微重力测量系统和高精度惯导系统中,也可用于高精度的静态角度测量系统中,将加速度用于大气反演方面尚无工程先例。如图5是国内工业级的用于惯导方面的加速度计,其分辨率量级约为10-6m/s2。
与加速仪原理类似的一种仪器是“阻力平衡仪”,它由弹簧系统取代静电系统,通过弹簧位移得到阻力,进而得到大气密度。该方法最初由意大利罗马大学率先提出,主要搭载于微小卫星平台。该方法曾为MSIS90模式的建立提供了阻力平衡器数据。
轨道数据反演
利用卫星轨道数据反演大气密度也是测量大气密度的一项重要来源,许多已知的半经验大气模型都是以这种方法计算得到的大气密度作为数据基础而建立的,如Jacchia和CIRA系列模型。其基本原理(图7)是:卫星在大气阻力的作用下会偏离预定轨道,利用轨道参数的变化与大气密度的关系即可反演出大气密度。
近十几年出现了利用卫星两行轨道根数(TLE)反演大气密度的方法。顾名思义,TLE数据是由两行轨道根数组成的,包含了卫星的编号、时间和六个轨道参数(轨道倾角、偏心率、每日绕行圈数、升交点赤经、近地点角距、平近点角)等信息。TLE数据从上世纪50年代起就开始记录,截止2017年6月底,在编目标已达42000多个,共累计数据1亿多组,因此TLE数据具有良好的空间覆盖率和时空连续性,这为大气密度的反演和研究大气密度的长期变化趋势提供了宝贵的数据来源。
TLE中的每日绕行圈数可以提供轨道衰减信息。卫星在大气阻力的作用下轨道逐渐衰减,导致每日绕行圈数逐渐增加,利用这一参数,并结合其他各项轨道参数与大气密度的关系即可计算出卫星轨道处的大气密度。
由于每颗卫星每天只记录有限的几组TLE数据,因此TLE的时间分辨率并不理想,针对这一问题,近几年又出现了利用卫星高精度GPS数据反演大气密度的方法,该方法原理与TLE类似,但由于GPS数据记录频率高,因此具有非常高的时间分辨率。该方法的缺点是对GPS数据的质量要求非常高,需要GPS精度达到厘米量级,而目前达到这一要求的卫星数据非常少,因此所能反演的大气密度十分有限。
以上是目前测量大气密度的几种常见方法,所得大气密度即可用来建立新的密度模型,又可对已有的模型进行修正。在实际应用过程中,卫星定轨和预报大部分是用大气模型来模拟大气密度的影响,而现有的模型普遍存在15%以上的误差,因此如何测量出更加准确的大气密度仍然任重而道远。