月球表面撞击坑亮温模拟.pdfVIP

2．方法与模型 2．1计算方法 图1大起伏月表亮温计算流程图(a)与微波辐射传输示意图(b) 本文提出的大起伏月表亮温计算方法【8】，如图1(a)所示。首先，基于DEM数字高程 的分辨率，将撞击坑划分为相应分辨率的网格单元。其次，利用改进的有效太阳辐照度模型 (包括遮蔽效应及倾斜平面的影响，如图1(b)所示)，计算得到每个网格单元的有效太阳 辐照度。再次，通过求解热传导方程得到每一个网格单元的温度剖面分布。接着，基于HUST 多层平面亮温模型【9】(包括倾斜平面的非零度观测角对亮温的影响，如图l(b)所示)，计 算得到每个网格单元的亮温。最后，在一个微波分辨单元内，结合嫦娥天线方向图【lo】对所 有的网格单元亮温进行积分，进而得到辐射计分辨单元内天线接收到的亮温。 2．2物理参数模型 鉴于篇幅，关于月表的物理模型包括物理温度和亮温模型可参见文献嘲，这里主要介绍 改进的太阳辐照度模型和介电常数。 2．2．1太阳辐照度模型 本文考虑起伏月表对当地太阳入射角的影响，采用余弦校正修正李雄耀等人【1l】提出的 月表有效太阳辐照度实时模型，即： ESI=kcos(／) (1) =E,m(cos口s证9+siIl口Cos臼cos国) =k(coso：sinO+sin￡zcosOcos(缈一y)) 式中，k为李等人计算的有效太阳辐照度，口为坡度角，9为太阳入射角，缈为太阳方 位角9与坡向角y的差。 2．2．2介电常数模型 认为月表以下6cm处为月尘和月壤的分界线，即热导率与密度的分界线，这种分层结构和 A15、A17的实测热流数据吻合的很好。本章中月球表层4cm为月尘层。 由于月壤的介电常数对频率非常敏感0s]，Montopoli用非线性重新拟合Apollo的样本并 提出了一种和频率有关的月壤介电常数模型【16，lrl： {sEI=1901pp(眇z)拂№Is“=g’．109uM2卜蜘 鳓。’ I,／=0．027，v=3．058和 g杪)=ql·f+q2 (3) 其中，q．=0．0272，f(GHz)是频率，q，=0．2967。 2．3辐射传输模型 本文中的辐射传输模型如图l(b)所示，在计算嫦娥天线接收的亮温时，既要考虑地 形的遮蔽和倾斜效应。地形的遮蔽主要是对月表物理温度分布产生很大的影响，而倾斜效应 既影响物理温度也影响月表的反射率。更多的细节，请参见文献is]。 3．结果与讨论 为了确定方法的适用性，分别对不同的月表地形进行仿真，包括，起伏较大的月球坑(年 轻撞击坑、古老撞击坑)和Apollo地区(A15、A17)。鉴于嫦娥二号的分辨率更高，因此， 本节用来验证的实测微波数据为嫦娥二号的2C级微波数据【lo】。仿真中，采用的DEM数据 Reconnaissance 为美国探测者轨道器Lunar Orbiter(LRO)上的LOLA激光高度计获取的 DEM数据【·8】。其中，模拟亮温的月球地方时采用实测微波数据的实际观测时间。仿真中， 反照率为0．1，红外发射率为O．94，月尘厚度为4cm。 3．1仿真结果 3．1．1古老撞击坑 一D虫nleadE∞J．∞n∞ulLJoc∞ ．54 ．53 —52 ．51 —50 latitude(。) ∞c∞d‘ng (e)Helicon坑的一号影像图 (f)Schiller坑的一号影像图 分别为0556和0572。 这里选择四个不同直径的撞击坑作为古老撞击坑的实例。从24kin．180km，分别选取了 和Schiller(51．9。S，