地表温度反演的不可能性.docxVIP

地表温度反演的不可能性 0 地表温度反演和地表热辐射传输方程 地表地表温度（lst）是表面能量平衡的重要参数。在表面温度和大气之间的相互作用中起着重要作用。它尤其适用于许多领域，如气、地质、水文、生态等。土壤水分状况,森林火灾的检测,军事伪装的应用,石油和铀矿的寻找等都离不开地表温度。通过反演可以得到地表温度。地表温度反演以地表热辐射传输方程为基础,即通过建立能量平衡方程来反演地表温度。辐射传输方程描述了卫星的微波辐射计所观测到的辐射总强度,不仅有来自地表的辐射 ,而且还有来自大气的向上和向下的路径辐射。这些辐射成分在穿过大气层到达遥感器的过程中 ,还受到大气层吸收作用的影响而削减。同时 ,地表和大气的辐射也在这一过程中产生不可忽略的影响。因此,地表温度的演算实际上是一个复杂的求解问题。 1 太阳光束自适应识别反射模型 辐射传输方程是描述辐射传播通过介质时与介质发生相互作用(吸收、散射、发射等)而使辐射能按一定规律传输的方程。在地球大气条件下,求解非常复杂,只能在一些假定下求得解析解,因此辐射传输方程的求解,一直是大气辐射学研究的重要内容。 辐射传输方程的基本形式为: R(vj?μ)=B(vj?ts)ε(vj?μ)+Ra(vj?μ)+∫2π0∫10μ′fr(μ?-μ′?φ′)Rd(vj?μ?-μ′?φ′)dμ′dφ′+τ*(vj?μ?-μ0?0)(-μ0)E0(vj)fr(μ?-μ0?φ0)。(1) R(vj,μ)是一个波段的平均光谱辐射亮度,vj是波数,B(vj,ts)是地面温度ts的普朗克方程,ε(vj,μ)是地面发射率,τ(vj,μ,Ps)是沿着θ角从地面到大气层顶部的透射率,Ra(vj,μ)表示大气上行辐射,Rd(vj,μ,-μ′,φ′)表示大气下行辐射,它的入射方向由-μ′和φ′表示(负号表示方向是向下的),μ′=cosθ′,φ′表示方位角,最后一项分式表示向下的太阳光束被地面反射,μ0是太阳天顶角的余弦值,φ0是观测方向与太阳光束之间的夹角。fr(μ,-μ′,φ′)是双向反射分布函数,E0(vj)表示从大气顶部入射的太阳光谱辐照度。 直接求解辐射传输方程是不可能的,也就是说它是不适定的。在通常情况下,不适定问题的解是不存在的。当解存在时,解并不稳定,即解不具有对初-边值、时-空几何以及模拟过程中容易产生误差的各种数据的连续性依赖。在研究不适定问题时,通常避开解的不适定性,而假设某些受到适当约束性限制的函数作为方程的解,求出解对前面提到的数据的连续性依赖。因此,对不适定性问题的研究,常常需要对解加上某些约束性限制。 为了反演物理参数的实用性,Xia.L.Ma等人做以下实用性的假设: (1)水平均一的大气层分布状态,大气热辐射不依赖于方位角。 (2)采用镜面反射模型,可假设μ=μ′,因此Rd(vj,μ,-μ′,φ′)也被简化。 (3)不考虑不同层次不同密度的单一材料或表面的双向漫发射因子,用于短波区域的太阳光束的反射。漫反射因子定义为a=πfr(μ,-μ′?φ′)r?r是假定的朗伯体的反射率。根据Kirchhoff定律,r=1-ε,因此方程(1)可以简化为: R(vj?μ)=B(vj?ts)ε(vj?μ)τ(vj?μ?ps)-∫Ρs0{B[vj?t(p)]?τ(vj?μ?p)?pdp+[1-ε(vj?μ)]∫Ρs0{B[vj?t(p)]?τ*(vj?μ?-μ?p)?p}dp+[1-ε(vj?μ)]aτ*(vj?μ?-μ0?0)(-μ0)E0(vj)π。(2) τ*(vj,μ,-μ,p)=τ(vj,μ,ps)τ(vj,-μ,ps)/τ(vj,μ,p),是反射光谱从地表到压力面p的透射率。假设地表是朗伯体,也就是说地表发射率不依赖于天顶角,不考虑太阳光束传播过程中的辐射亮度,那么辐射传输公式可以被近似成以下的扰动形式: δτB(j)=δtsΚts(j)+δε(j)Κε(j)+ls∑i=1δτ(i)Κt(i?j)+ls∑i=1{δln[q(i)]Κq(i?j)。(3) 扰动量δ是关于先验估计或平均值的扰动,τB是单个波段的亮温向量,Kts、Kε、Kt和Kq是地表温度 (ts),地表发射率(ε),大气温度(t)和水汽含量(q)的权重函数。i表示大气积分压力面(i=1,…,ls,从大气顶层到地表),j表示波带数,ls表示地表压力积分面,方程(3)可以表示为矩阵形式: δy=Κδx。(4)Κ=(ΚtΚqΚtsΚε)。(5)δx=(δtδlnqδtsδε)。(6) 其中δy=δtB,即为亮温扰动量,ts是地表温度。 2 辐射方程的数值解 2.1 地表温度反演 研究反演算法的基础是辐射传输方程,由于其解的不适定性,在分裂窗算法中,假定: (1)晴空大气条件; (2)地球表面为朗伯体; (3)忽略大气分子和气溶胶的散射。 则辐射强度表