張翠峰，萬(wàn)計(jì)，羅含之，楊碧蕊

(1.長(zhǎng)沙市規(guī)劃勘測(cè)設(shè)計(jì)研究院，湖南 長(zhǎng)沙 417000； 2.湖南省測(cè)繪產(chǎn)品質(zhì)量監(jiān)督檢驗(yàn)授權(quán)站，湖南 長(zhǎng)沙 417000；3.加州大學(xué)圣巴巴拉分校，加利福尼亞州 圣巴巴拉市 93106)

1 引 言

熱島效應(yīng)是指由于人口向城市聚集、城市建設(shè)規(guī)模擴(kuò)大造成熱量在局部空間范圍內(nèi)集中，并由此引起的一種地區(qū)性氣候現(xiàn)象。在城市的中心區(qū)域呈上升氣流，周邊污染物進(jìn)入城區(qū)嚴(yán)重影響空氣質(zhì)量，影響居民健康。同時(shí)，城市熱島還會(huì)影響生態(tài)環(huán)境的多個(gè)要素，包括土壤理化性質(zhì)、物質(zhì)循環(huán)、氣候、水文、能量代謝等[1]。熱島效應(yīng)強(qiáng)度作為國(guó)家園林城市和生態(tài)園林城市申報(bào)的評(píng)估指標(biāo)之一，一直受到相關(guān)部門(mén)的重視，尤其在我國(guó)碳達(dá)峰碳中和目標(biāo)實(shí)現(xiàn)過(guò)程中，城市熱島監(jiān)測(cè)是一種有效的監(jiān)測(cè)手段。衛(wèi)星遙感技術(shù)具有大面積的同步觀測(cè)、覆蓋范圍廣、可靠性高且成本低的特點(diǎn)，被廣泛應(yīng)用到城市地表溫度反演和城市熱環(huán)境的時(shí)空變化研究中[2]。劉施含[3]等在研究北京城市熱島效應(yīng)的過(guò)程中發(fā)現(xiàn)，建成區(qū)的擴(kuò)大和人流量和交通流量的增加會(huì)導(dǎo)致城市熱島的嚴(yán)重程度增加。Wang[4]等運(yùn)用單窗法，使用Landsat-8 TIRS傳感器的第10波段，反演了地表溫度，反演模型的平均誤差在 0.67 K。覃志豪[5]等利用MODIS數(shù)據(jù)，使用分裂窗算法對(duì)大區(qū)域的地表溫度反演問(wèn)題做了研究，結(jié)果表明分裂窗算法不僅效率高而且反演結(jié)果精確，但是MODIS數(shù)據(jù)的地面分辨率較低，只能反映宏觀數(shù)據(jù)。蘇偉忠[6]等利用Landsat數(shù)據(jù)反演地表溫度，結(jié)果表明城市熱島的空間分布受地表覆蓋類型的空間格局所影響，尤其是工業(yè)的分布。

2 研究區(qū)概況

長(zhǎng)沙作為我國(guó)中部地區(qū)迅速發(fā)展的中心城市之一，建成區(qū)不斷擴(kuò)大，城市化水平快速提升，伴隨著城市下墊面性質(zhì)也在發(fā)生快速轉(zhuǎn)變，城市生態(tài)環(huán)境受到極大考驗(yàn)。關(guān)于城市熱島效應(yīng)的研究也有不少成果[7，8]，主要集中在長(zhǎng)時(shí)間序列的熱島效應(yīng)時(shí)空演變特征分析，或者利用氣象資料分析熱島的變化特征，定量分析和地表結(jié)構(gòu)與熱環(huán)境的關(guān)系研究還比較少，且有待進(jìn)一步深入。本研究采用Landsat8衛(wèi)星的多光譜數(shù)據(jù)和熱紅外數(shù)據(jù)，基于單窗法反演城市地表溫度，不僅研究熱島效應(yīng)的時(shí)空變化特征，還對(duì)長(zhǎng)沙市內(nèi)不同土地利用類別與地表溫度的關(guān)系，下墊面性質(zhì)與熱島的分布和強(qiáng)度之間的定量關(guān)系開(kāi)展研究。

長(zhǎng)沙市地理位置介于東經(jīng)111°53′～114°15′，北緯27°51′～28°41′之間，屬亞熱帶季風(fēng)氣候?！笆濉逼陂g，長(zhǎng)沙城市建設(shè)速度明顯加快；“十四五”期間，長(zhǎng)沙將聚焦建設(shè) 1 000萬(wàn)常住人口、建成區(qū)面積近 1 000 km2的特大城市目標(biāo)，進(jìn)一步做大長(zhǎng)沙城市空間容量和人口規(guī)模。因此為分析長(zhǎng)沙城市熱島效應(yīng)特征，本文把市域主要建成區(qū)作為研究重點(diǎn)(圖1)。該區(qū)域人口密度大，交通發(fā)達(dá)，建筑物密集，呈現(xiàn)較大的城市規(guī)模，局部氣候特征明顯。長(zhǎng)沙作為國(guó)家的園林城市，同時(shí)也正在申報(bào)國(guó)家生態(tài)園林城市，對(duì)熱島效應(yīng)指標(biāo)有明確的約束條件，因此研究長(zhǎng)沙市熱環(huán)境效應(yīng)具有一定的現(xiàn)實(shí)意義和地理研究效益。

圖1 研究區(qū)域影像圖

3 溫度反演主要方法和參數(shù)選擇

3.1 衛(wèi)星傳感器和反演方法選擇

高質(zhì)量的衛(wèi)星遙感圖像決定了反演地表溫度的精度。由于熱紅外波段接收的信號(hào)受到大氣中的水汽含量的影響很大，所以無(wú)云或少云的長(zhǎng)沙衛(wèi)星影像是最優(yōu)選擇。另外，在選定衛(wèi)星遙感數(shù)據(jù)時(shí)，要考慮的因素是空間分辨率。衛(wèi)星影像較高的空間分辨率可以使反演的溫區(qū)空間范圍更加精準(zhǔn)，也為鎖定城市熱島熱點(diǎn)以及研究其形成原因提供了有利的條件。

通過(guò)對(duì)相關(guān)衛(wèi)星傳感器和波段數(shù)據(jù)的分析統(tǒng)計(jì)得出下表1，分析對(duì)比后選取了Landsat-8 TIRS傳感器的衛(wèi)星數(shù)據(jù)作為反演地表溫度的數(shù)據(jù)。未選取MODIS傳感器和AVHRR傳感器的原因是空間分辨率過(guò)低，未選取ASTER傳感器的原因是數(shù)據(jù)資源獲取較為麻煩，未選取高分五號(hào)全光譜成像儀(VIMS)是由于衛(wèi)星數(shù)據(jù)出現(xiàn)問(wèn)題，不適合研究使用。本研究選取了Landsat-8衛(wèi)星的6景影像，時(shí)間為2018年7月29日、2019年8月17日和2020年8月3日，每個(gè)時(shí)間兩景影像數(shù)據(jù)，可以覆蓋長(zhǎng)沙市主城區(qū)范圍。目前，利用熱紅外遙感反演地表溫度比較成熟的方法有單窗法[9，10]和劈窗法[11，12]兩種算法，衛(wèi)星數(shù)據(jù)僅包含一個(gè)熱紅外波段的(如Landsat-TM/ETM)用單窗算法，而具有多個(gè)熱紅外波段的遙感數(shù)據(jù)(如Terra-MODIS、NOAA-AVHRR等)則兩種方法都可以使用。

通過(guò)對(duì)反演算法和波段穩(wěn)定性的評(píng)估，本次研究選擇Landsat-8的第10波段數(shù)據(jù)通過(guò)單窗法反演地表溫度。根據(jù)表1可知，Landsat-8衛(wèi)星可用于反演地表溫度的熱紅外波段分別為第10波段(10.6 μm～11.2 μm)和第11波段(11.5 μm～12.5 μm)，而目前Landsat-8的第11波段經(jīng)常處于不穩(wěn)定狀態(tài)。因此選用第10波段來(lái)反演是目前較為主流的做法。選取單窗法的原因是Landsat系列單窗法反演研究多，精度評(píng)定高，反演模型的總體誤差在 1.4 K左右。

表1 具備地表溫度反演熱紅外波段的衛(wèi)星統(tǒng)計(jì)表

3.2 地表溫度反演主要參數(shù)

地表比輻射率是指相同溫度條件下，地表發(fā)射的輻射量與黑體發(fā)射的輻射量的比值。目前比較成熟的地表比輻射率計(jì)算方法有基于NDVI的比輻射率方法、基于分類的比輻射率方法以及獨(dú)立溫度光譜指數(shù)法等三種。本次研究使用基于NDVI的比輻射率方法來(lái)確定地表比輻射率。

大氣透過(guò)率也是衛(wèi)星熱紅外數(shù)據(jù)反演地表溫度的一個(gè)關(guān)鍵參數(shù)，是指電磁波經(jīng)大氣衰減后的電磁輻射通量與入射時(shí)電磁輻射通量的比值。本次研究通過(guò)NASA網(wǎng)站資源獲取大氣透過(guò)率相關(guān)參數(shù)。

4 研究方法

4.1 技術(shù)流程

本研究中城市熱島效應(yīng)的溫度反演流程圖如如圖2所示。

圖2 城市熱島效應(yīng)溫度反演流程圖

4.2 地表比輻射率計(jì)算

地表比輻射率作為區(qū)分星上亮溫和地表溫度的主要參數(shù)之一，被應(yīng)用到單窗法和劈窗法兩種地表溫度反演的方法中。本次研究我們將采用NDVI閾值法來(lái)計(jì)算地表比輻射率。這個(gè)方法將土地覆蓋大致分為三類：植被、土壤和介于植被與土壤之間的材料，依據(jù)它們不同的參數(shù)來(lái)計(jì)算地表比輻射率。我們將Landsat-8 OLI L1T產(chǎn)品的第4波段(Band4：0.64 μm-0.67 μm)和第5波段(Band5：0.85 μm-0.88 μm)提取出來(lái)，利用公式計(jì)算出長(zhǎng)沙城區(qū)的地表比輻射率(如下式)。ε為地表比輻射率，Ln為NDVI值。

4.3 大氣參數(shù)獲取

利用NASA官方網(wǎng)站，根據(jù)成像時(shí)間(GMT)、長(zhǎng)沙市中心點(diǎn)經(jīng)緯度(28.228N，112.939E)和其他的相關(guān)信息可以得到大氣在熱紅外波段的透過(guò)率τ、大氣向上輻射亮度L↑和大氣向上輻射亮度L↓。

表2 大氣透過(guò)率和大氣向上向下輻射值

4.4 溫度反演

根據(jù)普朗克公式計(jì)算地表溫度。如下：

其中，TS是為各波段對(duì)應(yīng)的亮度溫度(K)；B(TS)為黑體輻射亮度；K1、K2為L(zhǎng)andsat-8 TIRS的系數(shù)，從衛(wèi)星影像的頭文件中獲取。

最后得到反演結(jié)果如下(以2019年為例)：

圖3 長(zhǎng)沙市2019年8月17日地表溫度反演圖

4.5 溫區(qū)劃分

考慮到遙感數(shù)據(jù)獲取時(shí)間不同，太陽(yáng)輻射、云量等因素對(duì)反演的影響，本次所選三期遙感數(shù)據(jù)將使用平均溫度及標(biāo)準(zhǔn)差對(duì)反演的長(zhǎng)沙市地表溫度進(jìn)行分割。通過(guò)正態(tài)分布和標(biāo)準(zhǔn)差劃分溫區(qū)的方法，我們將研究范圍劃分成了低溫區(qū)、中低溫區(qū)、中溫區(qū)、中高溫區(qū)、高溫區(qū)，因此每年度熱島效應(yīng)溫區(qū)的劃分由不同日期的不同溫度分布決定。計(jì)算公式為：

T=T0±X×SD

T為地表溫度劃分閾值，T0為地表平均溫度值，SD為地表溫度標(biāo)準(zhǔn)差。T0和SD可以通過(guò)影像的溫度反演統(tǒng)計(jì)數(shù)據(jù)獲取(如表3所示)。根據(jù)地表溫度分布實(shí)際情況，將X取值定為0.5和1.5[10]，將地表溫度劃分為5個(gè)等級(jí)(如表4所示)

表3 各年份地表平均溫度值及地表溫度標(biāo)準(zhǔn)差(℃)

表4 各年份溫區(qū)劃分表(℃)

根據(jù)各年份溫區(qū)劃分表制作2018年～2020年熱島效應(yīng)溫區(qū)圖，各年份溫區(qū)圖如圖4～圖6所示(白色部分為云掩膜去除的影像數(shù)據(jù)，不參與整體統(tǒng)計(jì)分析)：

圖4 2018年長(zhǎng)沙市熱島效應(yīng)溫區(qū)圖

圖5 2019年長(zhǎng)沙市熱島效應(yīng)溫區(qū)圖

圖6 2020年長(zhǎng)沙市熱島效應(yīng)溫區(qū)圖

將每年度按劃分后的溫區(qū)計(jì)算面積，統(tǒng)計(jì)得到下表5：

表5 各年份溫區(qū)面積統(tǒng)計(jì)表

5 數(shù)據(jù)分析

5.1 熱島強(qiáng)度變化分析

通過(guò)分析逐年溫區(qū)的比例變化和面積變化，可以研究熱島強(qiáng)度。根據(jù)各年份溫區(qū)面積表制作2018年～2020年溫區(qū)比例變化圖，從圖7可以看出，近三年長(zhǎng)沙市的熱島效應(yīng)有所緩解。

圖7 2018年～2020年溫區(qū)比例變化

5.2 地類抽樣分析

在長(zhǎng)沙市內(nèi)選取各種地類中具有代表性的地點(diǎn)記錄他們的平均溫度，分析地類與地溫的聯(lián)系，如圖8～圖10所示。本研究選取的地類是：居民居住用地、商業(yè)用地、農(nóng)業(yè)用地，林地，草地，政府用地，交通用地，水體。

圖8 2018年采樣點(diǎn)平均溫度

圖9 2019年采樣點(diǎn)平均溫度

圖10 2020年采樣點(diǎn)平均溫度

從以上3年的8種地類的地表溫度可以得出：水域溫度為地類中最低，交通用地的溫度最高，植被覆蓋率高的林地、農(nóng)業(yè)、公園用地溫度較低，住宅、政府、商業(yè)等建筑指數(shù)高的地類溫度較高。

5.3 下墊面性質(zhì)分析

城市下墊面性質(zhì)會(huì)影響熱島的分布和強(qiáng)度。本次研究選取了約 1 200個(gè)樣本點(diǎn)，根據(jù)樣本點(diǎn)位的正規(guī)化指數(shù)值和地表溫度值進(jìn)行趨勢(shì)擬合，得到各指數(shù)與地表溫度的線性回歸方程，如圖11～圖13所示。

圖11 地表含水散點(diǎn)趨勢(shì)圖

圖12 植被生長(zhǎng)狀況散點(diǎn)趨勢(shì)圖

圖13 地表建設(shè)程度散點(diǎn)趨勢(shì)圖

從以上三圖可以看出，地表水體、植被、建筑與地表溫度有明顯的相關(guān)關(guān)系。地表水體和植被覆蓋與地表溫度之間具有負(fù)相關(guān)關(guān)系，城市建筑與地表溫度之間具有正相關(guān)關(guān)系。簡(jiǎn)單來(lái)說(shuō)，地表含水量越高，地表溫度越低；地表植被生長(zhǎng)狀況越好，越密集，地表溫度越低；地表建設(shè)程度越高，地表溫度越高。因此，改變城市下墊面性質(zhì)可以有效治理城市熱島問(wèn)題。

6 結(jié) 論

城市熱島效應(yīng)受到多個(gè)影響因素的共同作用，進(jìn)行其形成機(jī)制的多因素綜合研究非常必要。

隨著城市的規(guī)?；l(fā)展，建成區(qū)面積的不斷擴(kuò)張，改變了原有的地表覆蓋特征和下墊面性質(zhì)，這是熱島效應(yīng)形成的主要因素。NDBI(歸一化建筑指數(shù))和交通用地與地表溫度升高關(guān)系顯著，是長(zhǎng)沙市城市熱島形成的主要影響因子。從下墊面的性質(zhì)分析得出的回歸方程可以表明，長(zhǎng)沙主城區(qū)的建筑物指數(shù)增加0.1，將會(huì)導(dǎo)致地表溫度增加2.85℃；植被指數(shù)增加0.1，將會(huì)導(dǎo)致地表溫度降低大約1.72℃；水體指數(shù)增加0.1，將會(huì)導(dǎo)致地表溫度降低大約1.56℃。因此保護(hù)和增加城市濕地和綠地面積是降低地表溫度、緩解城市熱島的有效途徑。尤其是在熱島效應(yīng)明顯的道路旁建設(shè)一定規(guī)模的綠帶，在建構(gòu)筑物密集的城區(qū)建設(shè)一定規(guī)模的綜合公園、海綿公園、帶狀公園等綠地，通過(guò)科學(xué)的城市綠地設(shè)計(jì)能力最大化地發(fā)揮植被、水體的降溫作用，能夠有效降解城市的熱污染，解決熱島效應(yīng)治理問(wèn)題。

7 進(jìn)一步思考

城市熱島效應(yīng)研究還可以根據(jù)其應(yīng)用發(fā)展來(lái)進(jìn)一步拓展研究深度和廣度。

7.1 綠色建筑識(shí)別

從2019年8月17日的多色溫區(qū)圖(圖14)中，我們發(fā)現(xiàn)長(zhǎng)沙地標(biāo)建筑IFS(長(zhǎng)沙國(guó)際金融中心)被劃為中溫區(qū)/中高溫區(qū)，平均溫度為39.79℃，顯著低于其周邊建成區(qū)溫度。經(jīng)過(guò)分析發(fā)現(xiàn)，IFS為L(zhǎng)EED(能源與環(huán)境設(shè)計(jì)先鋒)認(rèn)證的鉑金級(jí)(最高級(jí))綠色建筑，碳排放量少，因此在溫區(qū)圖中它的溫度相較其他市中心傳統(tǒng)商業(yè)建筑的溫度較低。

圖14 長(zhǎng)沙IFS局部溫區(qū)示意圖

7.2 高層建筑和低層建筑地表溫度差值分析

基于每年度熱島效應(yīng)溫區(qū)圖，可以對(duì)比分析高層建筑與低層建筑之間的地表溫度差異，對(duì)城市規(guī)劃建設(shè)提供參考決策依據(jù)。

7.3 長(zhǎng)沙市熱島效應(yīng)強(qiáng)度

按照國(guó)家園林城市標(biāo)準(zhǔn)對(duì)城市熱島效應(yīng)的考核要求，可以進(jìn)一步開(kāi)展長(zhǎng)沙市熱島效應(yīng)強(qiáng)度指標(biāo)計(jì)算。