冷 佩，廖前瑜，※，任 超，李召良

（1. 中國(guó)農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所，北京100081；2. 桂林理工大學(xué)測(cè)繪地理信息學(xué)院，廣西桂林541004）

0 引言

近地表氣溫是指在近地表（一般指距地面2 m 處高度）觀測(cè)的大氣溫度。近地表氣溫作為描述地表大氣環(huán)境的重要指標(biāo)，控制著自然系統(tǒng)中大多數(shù)的生物和物理過程［1-2］，是各種地表過程模型，如地表蒸散發(fā)模型［3-4］、水文模型［5］、土壤—植被—水分系統(tǒng)動(dòng)力模型［6］等的重要驅(qū)動(dòng)參數(shù)。在全球氣溫變暖的大背景下，區(qū)域氣溫的持續(xù)攀升導(dǎo)致了生態(tài)環(huán)境的改變和氣象災(zāi)害頻發(fā)，尤其是對(duì)農(nóng)業(yè)生產(chǎn)的可持續(xù)發(fā)展產(chǎn)生了重大影響［7-10］。因此，準(zhǔn)確估計(jì)近地表氣溫在地球表面空間上的連續(xù)分布對(duì)全球農(nóng)業(yè)發(fā)展具有重要意義。

氣溫是表征空氣冷熱的物理量。近地表氣溫的變化主要受到輻射、湍流和平流等多種因素影響［11］。目前，分布在世界各地的氣象站能夠提供較為精確而連續(xù)時(shí)間序列的近地表氣溫?cái)?shù)據(jù)，然而在稀疏空間分布的站點(diǎn)氣溫信息無法描述近地表氣溫在連續(xù)空間上的異質(zhì)性。對(duì)于沒有氣象站點(diǎn)的地區(qū)，傳統(tǒng)方法大多通過距離最近的有限的氣象站氣溫觀測(cè)數(shù)據(jù)插值得到［12-18］。相比傳統(tǒng)觀測(cè)，遙感方法可以獲取區(qū)域尺度的近地表氣溫影響因子，具有實(shí)時(shí)性、區(qū)域性、經(jīng)濟(jì)適用性等無可替代的優(yōu)勢(shì)。遙感技術(shù)的廣泛應(yīng)用為獲取區(qū)域尺度近地表氣溫提供了強(qiáng)大的數(shù)據(jù)支撐。因此，眾多學(xué)者圍繞近地表氣溫遙感反演進(jìn)行了大量研究，并已形成了眾多基于遙感數(shù)據(jù)或與遙感數(shù)據(jù)結(jié)合的區(qū)域氣溫反演方法，這些主要可以歸納為4 類：大氣廓線外推方法、地表溫度—植被指數(shù)方法（Temperature Vegetation Index，TVX）、基于數(shù)據(jù)統(tǒng)計(jì)的方法以及基于地表能量平衡的方法。該文系統(tǒng)回顧近地表氣溫遙感反演算法模型及最新研究進(jìn)展，為相關(guān)領(lǐng)域獲取區(qū)域近地表氣溫提供借鑒，對(duì)近地表氣溫遙感反演模型未來發(fā)展方向進(jìn)行展望。

1 近地表氣溫遙感反演方法

1.1 大氣廓線外推方法

大氣廓線外推方法建立在同一垂直面空氣溫度是連續(xù)變化的假設(shè)之上。在此假設(shè)基礎(chǔ)之上，通過絕熱遞減率（Adiabatic Lapse Rate，ALR）插值得到近地表氣溫。近地表氣溫與ALR之間的關(guān)系可以表示為：

式（1）中Ta為近地表氣溫，H為地面海拔高度，H0為大氣廓線離地面最近的氣壓層高度，T0則為H0點(diǎn)對(duì)應(yīng)的氣溫。

目前，MOD07_L2 是較為常用的基于大氣廓線外推方法反演近地表氣溫的數(shù)據(jù)源。MOD07_L2 可提供全球每日分布在20 個(gè)垂直大氣壓水平上（1 000、950、920、850、780、700、620、500、400、300、250、200、150、100、70、50、30、20、10、5 hPa）的大氣溫度剖面數(shù)據(jù)，具有較高水平空間分辨率（5 km），但是這些數(shù)據(jù)僅限于晴空條件下。因此，在使用MOD07_L2 大氣廓線參數(shù)時(shí)，需要對(duì)云干擾像元進(jìn)行掩膜處理和質(zhì)量控制。Bisht 等［19］利用垂直氣壓水平1 000 hPa 下的空氣溫度作為近地表氣溫。然而Tang等［20］認(rèn)為由于海拔高度的變化，利用氣壓水平1 000 hPa 下的空氣溫度替代近地表氣溫顯然是不合適的，并提出了通過大氣流體靜力學(xué)的方法獲取近地表氣溫。Cosgrove 等［21］基于MOD07_L2 產(chǎn)品開發(fā)了一個(gè)簡(jiǎn)單的大氣廓線外推參數(shù)化方法，將ALR 設(shè)定為每公里上升6.5℃或者下降6.5℃。因此，根據(jù)ALR 的假設(shè)似乎能夠獲取晴空下垂直面上任意高度的空氣溫度。然而固定的ALR 假設(shè)并不能準(zhǔn)確地估計(jì)區(qū)域尺度的近地表氣溫，因?yàn)锳LR 在空間分布上是變化的；除此之外，ALR 還存在日變化、顯著的季節(jié)變化和年變化［22］。因此，基于ALR 是常數(shù)的假設(shè)僅僅適用于均一下墊面和特定的氣壓水平區(qū)域［23］。Zhu 等［24］基于MOD07_L2 產(chǎn)品大氣廓線利用最接近地表的兩個(gè)氣壓層和所對(duì)應(yīng)的空氣溫度逐像元計(jì)算ALR，然后通過地表氣壓插值得到一個(gè)氣溫，并利用估算得到的氣溫與地表溫度（Land Surface Temperature，LST）取平均得到晴空下地面2 m 的氣溫，最后將晴空下的氣溫和MOD06_L2 產(chǎn)品的LST 建立線性回歸關(guān)系，有云情況下根據(jù)回歸系數(shù)得到有云條件下地面2 m 的氣溫，全天候氣溫反演結(jié)果均方根誤差范圍為2.4～3.3℃。Famiglietti 等［25］基于MOD07_L2 大氣廓線利用地表大氣壓非線性內(nèi)插得到近地表氣溫，并將結(jié)果與全球范圍內(nèi)109 個(gè)地面氣象站觀測(cè)數(shù)據(jù)進(jìn)行了對(duì)比驗(yàn)證，平均均方根誤差為3.47℃。大氣廓線外推的方法原理簡(jiǎn)單，不需要地面輔助數(shù)據(jù)就能得到精度較高的晴空下近地表氣溫。有云條件下則需要其他的遙感數(shù)據(jù)建模估算。

1.2 地表溫度—植被指數(shù)方法

1989 年Nemani 等［26］人首次提出地表溫度-植被指數(shù)特征空間模型。圖1 為TVX 示意圖。由于濃密植被和周圍空氣的熱容接近，TVX 方法估算近地表氣溫是基于遙感圖像植被全覆蓋時(shí)的地表溫度近似等于該地點(diǎn)空氣溫度的假設(shè)。因此，LST隨著植被覆蓋度的增加趨近于該點(diǎn)的空氣溫度，即使不發(fā)生蒸散發(fā)，全植被覆蓋冠層的輻射溫度也與冠層內(nèi)空氣溫度保持平衡［1，27］?；诖耍琓VX 方法能夠成功地應(yīng)用于空氣溫度的估算有兩個(gè)限制條件［28-31］：（1）在晴空條件下具有較為穩(wěn)定的大氣條件；（2）地表在一定范圍內(nèi)滿足土壤水分均一的條件。TVX 模型由于其物理原理明確，模型簡(jiǎn)單，不需要輔助數(shù)據(jù)的特點(diǎn)，被廣泛應(yīng)用于區(qū)域近地表氣溫的估算。Zhu 等［29］將植被指數(shù)與地表溫度的負(fù)相關(guān)系數(shù)閾值降至0.1后，成功地利用TVX 方法從MODIS 陸地產(chǎn)品中提取出了日最高氣溫和日最低氣溫，均方根誤差為3.79℃。Renaud 等［32］基于TVX 模型估算了城市植被氣溫，均方根誤差范圍在3.34～6.23℃。

圖1 TVX 示意圖Fig.1 TVX configuration

TVX 方法估算近地表氣溫一般是通過建立像元陣列的NDVI 和LST 回歸關(guān)系，然后外推求得理論植被全覆蓋情況下的LST，最后求得各個(gè)像元的近地表氣溫。關(guān)系式如下：

式（2）、（3）中LST為地表溫度，NDVI為歸一化植被指數(shù)，ai和bi分別為每個(gè)移動(dòng)窗口內(nèi)的LST和NDVI線性回歸系數(shù)的截距和斜率，Ta為近地表氣溫，NDVImax為理論植被全覆蓋植被指數(shù)。

1.3 基于數(shù)據(jù)統(tǒng)計(jì)的方法

1.3.1 單因子統(tǒng)計(jì)模型

單因子統(tǒng)計(jì)模型是基于空氣溫度和地表溫度的高度相關(guān)性，通過建立空氣溫度和遙感數(shù)據(jù)地表溫度的線性回歸關(guān)系求得氣溫。Chen 等［33］建立了基于GOES（Geostationary Operational Environmental Satellite）靜止氣象衛(wèi)星反演的地表溫度和地面站點(diǎn)實(shí)測(cè)空氣溫度的線性回歸關(guān)系，統(tǒng)計(jì)結(jié)果顯示地表溫度和空氣溫度之間平均相關(guān)系數(shù)為0.87，估算的空氣溫度標(biāo)準(zhǔn)方差為1.57℃。侯英雨等［34］基于AVHRR（Advanced Very High Resolution Radiometer）衛(wèi)星的地表溫度與觀測(cè)空氣溫度的相關(guān)關(guān)系，建立了稀疏植被下墊面不同高程范圍的空氣溫度單因子統(tǒng)計(jì)模型，經(jīng)檢驗(yàn)反演結(jié)果絕對(duì)誤差在1.5～1.8℃之間。單因子統(tǒng)計(jì)模型可以表述為：

該方法的優(yōu)點(diǎn)在于模型簡(jiǎn)單，輸入?yún)?shù)少。然而其回歸系數(shù)依賴于建模數(shù)據(jù)的獲取時(shí)間和地點(diǎn)，無法推廣到區(qū)域尺度。Lin 等［35］的研究表明在不同地形和不同時(shí)刻條件下，空氣溫度和地表溫度的關(guān)系是變化的，在解釋復(fù)雜地形表面能量通量的時(shí)空變化時(shí)必須要考慮相關(guān)的環(huán)境因子。因此，僅考慮地表溫度對(duì)空氣溫度的影響的統(tǒng)計(jì)模型并不適用于所有時(shí)刻和區(qū)域的空氣溫度估算。

1.3.2 多因子統(tǒng)計(jì)模型

多因子統(tǒng)計(jì)模型考慮了多個(gè)空氣溫度的影響因子建立線性或非線性模型求解空氣溫度。影響空氣溫度的因素主要可歸類為兩類：（1）地理變量包括海拔、緯度和離海岸線距離等；（2）物理變量（遙感變量）包括地表溫度、植被指數(shù)和地表反照率等。多因子統(tǒng)計(jì)模型可以表述為：

式（5）中Albedo為地表反照率，h為高程。

Kawashima 等［36］在單因子統(tǒng)計(jì)模型基礎(chǔ)之上，引入了NDVI，通過平均地表溫度和歸一化植被指數(shù)的二元線性回歸分析得到了氣溫，反演標(biāo)準(zhǔn)誤差在0.47～1.05℃。Cristóbal等［37］發(fā)現(xiàn)結(jié)合遙感變量（地表溫度、地表反照率和植被指數(shù)）和地理變量（海拔、緯度和太陽輻射）的空氣溫度模型要比單純地理模型和單純遙感模型精度要高，混合模型日尺度氣溫反演結(jié)果的平均均方根誤差為1.75℃，月尺度和年尺度的氣溫反演結(jié)果平均均方根誤差為1℃。徐偉燕等［38］考慮了地表溫度、植被指數(shù)、水汽壓、地表反照率及高程5 個(gè)影響因子建立多元線性回歸模型估算近地表氣溫，反演結(jié)果均方根誤差為2.31℃，同時(shí)對(duì)其分白天和夜晚進(jìn)行了檢驗(yàn)，發(fā)現(xiàn)夜間精度比白天高。白琳等［39］利用地表溫度、植被指數(shù)、差異水體指數(shù)、地表反照率、不透水面蓋度和高程6 個(gè)影響因子，將隨機(jī)森林方法應(yīng)用于近地表氣溫的反演，反演結(jié)果的均方根誤差為1.06℃，發(fā)現(xiàn)多因子的隨機(jī)森林方法比單純的多元線性回歸模型精度高。Jang 等［40］將通過AVHRR 的5 個(gè)波段反演得到的反射率和亮溫?cái)?shù)據(jù)，結(jié)合地理變量作為神經(jīng)網(wǎng)絡(luò)的輸入層來反演空氣溫度，結(jié)果顯示隱層節(jié)點(diǎn)數(shù)為22 時(shí)精度最高，均方根誤差為1.79℃，有95%的樣本誤差在3℃以內(nèi)。Chenour 等［41］采用地理加權(quán)回歸方法比較了單變量和多變量模型的精度，結(jié)果顯示采用地表溫度、植被指數(shù)、高程、坡度及反照率的地理加權(quán)回歸方法進(jìn)行大氣估算精度最高，平均均方根誤差為0.62℃，平均相關(guān)系數(shù)為0.99。大量的研究表明多因子的統(tǒng)計(jì)模型精度要比單因子的統(tǒng)計(jì)模型高。但是可移植性仍然較差，而且并不能反映湍流和平流對(duì)空氣溫度的影響。

1.4 基于地表能量平衡方法

忽略平流的影響，根據(jù)能量平衡方程可以建立近地表氣溫的表達(dá)式：

式（6）中，Ta為近地表氣溫，T0為空氣動(dòng)力學(xué)溫度，H為顯熱通量，ρ為空氣密度，Cp為空氣定壓比熱，ra為空氣動(dòng)力學(xué)阻抗。由于空氣動(dòng)力學(xué)溫度在實(shí)際中難以獲取，因此通常假設(shè)地表溫度等于空氣動(dòng)力學(xué)溫度?？諝鈩?dòng)力學(xué)溫度與地表溫度之間的差異會(huì)隨著植被覆蓋度降低而變大，這一假設(shè)應(yīng)用在稀疏植被下墊面和裸土下墊面時(shí)會(huì)帶來較大的誤差［42-43］。

Pape 等［44］基于能量平衡方程，并結(jié)合地面實(shí)測(cè)的氣象數(shù)據(jù)和遙感數(shù)據(jù)求解得到近地表氣溫，均方根誤差在0.37～1.02℃之間，但是其僅僅考慮了能量閉合系統(tǒng)內(nèi)的湍流影響和輻射影響。然而，在實(shí)際情況中地表能量平衡方程通常是不閉合的，平流作用也會(huì)對(duì)氣溫產(chǎn)生影響。Zhang 等［45］提出了一個(gè)基于能量平衡思想估算區(qū)域氣溫的算法（Advection-Energy Balance for Surface Air Temperature Algorithm，ADEBAT），同時(shí)考慮了局部驅(qū)動(dòng)力和平流的影響，對(duì)于區(qū)域內(nèi)每一個(gè)待反演像元，尋找距離其最近的滿足有溫差和風(fēng)速風(fēng)向相近條件的兩個(gè)氣象站點(diǎn)估算平流因子，結(jié)果均方根誤差在0.2～0.42 K。Liu 等［46］通過反距離加權(quán)插值改進(jìn)了ADEBAT 算法的平流因子。Zhu 等［47］提出了一種基于地表能量平衡算法（Surface Energy Balance Algorithm for Land，SEBAL）的晴空氣溫反演方法，通過對(duì)極端像元（水面像元和干燥裸土像元）的迭代計(jì)算，建立地表溫度和氣溫的線性關(guān)系，利用MODIS 數(shù)據(jù)反演得到的結(jié)果均方根誤差為1.7℃，利用Landsat8數(shù)據(jù)反演得到的結(jié)果均方根誤差為2.6℃。

2 近地表氣溫遙感反演方法優(yōu)缺點(diǎn)比較

大氣廓線外推方法簡(jiǎn)單易行，可操作性強(qiáng)，不需要地面輔助數(shù)據(jù)就可以獲取晴天條件下的區(qū)域氣溫，但是其結(jié)果精度相對(duì)于其他3 種方法仍然有待提高。有云條件下，可以通過已獲取的晴空下空氣溫度和地表溫度建立回歸關(guān)系求得云干擾下的空氣溫度［24］。與此同時(shí)，面臨的另外一個(gè)問題是通過熱紅外并不能獲取云干擾下的地表溫度數(shù)據(jù)，實(shí)際上，與空氣溫度建立回歸關(guān)系的地表溫度是使用統(tǒng)計(jì)或時(shí)空插值方法填充獲取的，僅僅是理論無云條件下的地表溫度近似值，而不是實(shí)際觀測(cè)的云下地表溫度［48］。大氣廓線外推方法由于其特點(diǎn)，得到了廣泛應(yīng)用，特別是在復(fù)雜陸地過程模型中，常常將大氣廓線外推方法獲取的空氣溫度作為輸入?yún)?shù)之一［20，23］。

TVX 方法的成功應(yīng)用建立在兩個(gè)假設(shè)條件下：（1）地表溫度和植被指數(shù)之間存在很強(qiáng)的負(fù)相關(guān)關(guān)系；（2）全植被冠層的溫度與冠層內(nèi)的空氣溫度處于平衡狀態(tài)。然而，研究表明TVX 方法雖然依賴于地表溫度和植被指數(shù)之間的負(fù)相關(guān)關(guān)系，但是這種依賴關(guān)系并不十分敏感，相關(guān)系數(shù)閾值可低至0.1［28］。TVX 方法的一個(gè)關(guān)鍵是求得NDVImax，研究表明NDVImax隨著研究區(qū)和衛(wèi)星圖像的不同而變化［1］，一般通過地面輔助數(shù)據(jù)和遙感數(shù)據(jù)之間的經(jīng)驗(yàn)關(guān)系來判定NDVImax。而且當(dāng)該方法應(yīng)用于低植被覆蓋度區(qū)域和裸土區(qū)域時(shí)，由于移動(dòng)窗口內(nèi)的地表溫度和植被指數(shù)變化范圍較小而導(dǎo)致構(gòu)成的線性關(guān)系變得隨機(jī)。因此，TVX 方法不適用于低植被覆蓋度或非植被覆蓋下墊面。盡管TVX 方法得益于簡(jiǎn)單的模型和輸入，但NDVImax的經(jīng)驗(yàn)性大大限制了TVX 模型的估算精度和實(shí)際應(yīng)用。

數(shù)據(jù)統(tǒng)計(jì)類的方法是基于地面觀測(cè)數(shù)據(jù)和遙感數(shù)據(jù)建立空氣溫度與其影響因子的一種線性或非線性的經(jīng)驗(yàn)關(guān)系。該方法原理簡(jiǎn)單，需要大量的輸入數(shù)據(jù)，結(jié)果精度受到輸入數(shù)據(jù)的影響。由于地面氣象站點(diǎn)的疏密程度影響，可移植性差是這類模型的最大特點(diǎn)。隨著計(jì)算機(jī)技術(shù)的發(fā)展和大數(shù)據(jù)時(shí)代的到來，數(shù)據(jù)統(tǒng)計(jì)模型有了較大的發(fā)展，逐漸向更加復(fù)雜的非線性算法發(fā)展［39］。基于數(shù)據(jù)統(tǒng)計(jì)的方法由于其簡(jiǎn)單的原理，常常被引入到其他3 種氣溫反演模型中用來提高氣溫的反演精度。

基于地表能量平衡的方法估算氣溫則具有良好的普適性和物理基礎(chǔ)，并且具有較高的反演精度。但是模型復(fù)雜，參數(shù)過多導(dǎo)致而容易造成誤差傳遞和累積，且具有遙感難以獲取的地面參數(shù)，如風(fēng)速、空氣動(dòng)力學(xué)阻抗和空氣動(dòng)力學(xué)溫度等。在實(shí)際情況中，地表能量平衡通常都是不閉合的，需要考慮湍流和平流的影響。因此，地表能量平衡模型通常需要地面輔助數(shù)據(jù)計(jì)算湍流因子和平流因子對(duì)其結(jié)果進(jìn)行校正，限制了其在大區(qū)域范圍上的應(yīng)用。

表1 給出了上述4 類近地表氣溫遙感方法的對(duì)比?？偟膩碚f，每種方法都有各自的優(yōu)缺點(diǎn)。然而，目前還沒有一種方法既能夠獨(dú)立于（或較少依賴于）地面輔助數(shù)據(jù)而又具有普遍適用性。

表1 不同近地表氣溫遙感反演方法對(duì)比Table 1 Comparison of methods for estimating near-surface air temperature

3 結(jié)論與展望

該文系統(tǒng)回顧和總結(jié)了近地表氣溫遙感反演的研究進(jìn)展，并將這些方法歸納為4 類：大氣廓線外推方法、地表溫度—植被指數(shù)方法（TVX）、基于數(shù)據(jù)統(tǒng)計(jì)的方法以及基于地表能量平衡的方法。前3 種方法是基于空氣溫度和其他變量之間的經(jīng)驗(yàn)關(guān)系，易于理解和使用；而地表能量平衡模型則依賴于空氣溫度的物理形成機(jī)制，較為復(fù)雜，需要大量輸入數(shù)據(jù)，具有良好的可移植性和通用性。綜合目前國(guó)內(nèi)外近地表氣溫遙感反演的發(fā)展現(xiàn)狀，對(duì)近地表氣溫遙感反演算法未來的發(fā)展方向做出展望。

（1）在晴空條件下，上述4 類方法都有一定適用性。然而現(xiàn)有的氣溫遙感反演算法除了大氣廓線外推方法之外，其余3 類方法一般都依賴地面觀測(cè)數(shù)據(jù)與遙感數(shù)據(jù)共同求解。距離通過純遙感數(shù)據(jù)反演得到高精度的全球氣溫還有一定的距離。同時(shí)TVX 方法和基于地表能量平衡的方法不適用于低植被覆蓋區(qū)域和裸土區(qū)域。未來應(yīng)該加強(qiáng)不同地形條件下的大氣輻射傳輸過程研究。此外，晴空條件下，地表溫度的遙感反演技術(shù)已較為成熟，而地表溫度與近地表氣溫密切相關(guān)。考慮到靜止氣象衛(wèi)星具有較高的時(shí)間分辨率，而現(xiàn)有氣溫反演算法較少考慮地表溫度日變化信息。因此，將地表溫度和其他地表參數(shù)的日變化信息引入氣溫反演中，有可能是提高晴空條件下氣溫反演精度的有效途徑。

（2）在云覆蓋情況下，目前僅有數(shù)據(jù)統(tǒng)計(jì)方法適用。面臨的主要困難是在有云條件下，無法獲取地表溫度、地表反照率以及植被指數(shù)等諸多氣溫反演的重要量。這也是定量遙感領(lǐng)域面臨的最大挑戰(zhàn)之一。盡管目前已經(jīng)發(fā)展了大量估算云下地表參數(shù)的算法，但大多數(shù)是基于時(shí)間和空間的插值填充獲取的，并不能反映實(shí)際云覆蓋條件下的地面真值。微波由于能穿透云層的特點(diǎn)，被廣泛應(yīng)用于獲取云下地表參數(shù)。因此，作者認(rèn)為未來應(yīng)該加強(qiáng)有云條件下的輻射傳輸過程研究，發(fā)展利用微波遙感數(shù)據(jù)獲取區(qū)域尺度近地表氣溫的方法。

（3）無論是在晴空還有云覆蓋條件下，現(xiàn)有的近地表氣溫遙感反演算法基本上忽略了平流對(duì)空氣溫度的影響，然而在實(shí)際中平流對(duì)空氣溫度的影響是顯著的。平流因子可以描述為關(guān)于風(fēng)速的函數(shù)，但是目前還沒有可行的風(fēng)速反演方法。現(xiàn)有的算法大多通過地面觀測(cè)氣溫與方程求解氣溫之間的差分估算平流因子。因此大氣流動(dòng)對(duì)氣溫影響的量化建模是今后近地表氣溫遙感反演能否擺脫地面輔助數(shù)據(jù)，建立可靠的普適氣溫反演模型的關(guān)鍵。