李可相，肖玖軍，陳 遠(yuǎn)，謝 剛，張藍(lán)月,董艷艷

(貴州科學(xué)院貴州省山地資源研究所，貴州 貴陽(yáng) 550001)

1 引言

地表溫度是城市下墊面能量平衡的核心，也是城市熱島研究的關(guān)鍵。遙感監(jiān)測(cè)法是目前宏觀尺度獲取城市地表溫度，研究城市熱島效應(yīng)的主要研究方法之一，具有監(jiān)測(cè)面積廣、時(shí)間序列長(zhǎng)、可在時(shí)空上進(jìn)行長(zhǎng)期連續(xù)監(jiān)測(cè)等優(yōu)點(diǎn)。許多研究者對(duì)此進(jìn)行了研究，如裴志方[1]等以鄭州市為例研究了城市化下城市熱環(huán)境與下墊面關(guān)系，結(jié)果表明城市化下城市熱島效應(yīng)明顯，地表溫度與植被覆蓋高度負(fù)相關(guān)。魏冠軍[2]基于Landsat TM/ETM +對(duì)蘭州市城市熱環(huán)境格局動(dòng)態(tài)變化進(jìn)行了研究分析。侯浩然[3]等基于遙感方法對(duì)福州市近20年城市熱環(huán)境變化進(jìn)行了研究分析。多是利用Landsat TM/ETM +等高空間分辨率的遙感數(shù)據(jù)，基于高時(shí)間分辨率遙感的長(zhǎng)期監(jiān)測(cè)研究較少見，欠缺針對(duì)地表溫度的季節(jié)變化特征研究。

本文在前人研究的基礎(chǔ)上，選取典型喀斯特山地城市(貴陽(yáng)市)作為研究區(qū)域，采用近20年的MODIS遙感數(shù)據(jù)，對(duì)貴陽(yáng)市地表熱環(huán)境進(jìn)行研究，探討其城市熱島效應(yīng)的分布及強(qiáng)度，分析研究區(qū)域地表熱環(huán)境的長(zhǎng)期變化情況，以期為緩解城市熱島效應(yīng)提供參考。

2 研究區(qū)域與研究方法

2.1 研究區(qū)概況

貴陽(yáng)市地處長(zhǎng)江和珠江水系分水嶺地帶，云貴高原東斜坡，東部平原向西部高原的過(guò)渡地帶，屬巖溶峰叢洼地、峰叢槽谷和峰林谷地地貌，海拔高，緯度低，是一座典型的喀斯特發(fā)育完全、城市生態(tài)脆弱的高原山區(qū)城市[4]。轄6個(gè)區(qū)，3個(gè)縣，1個(gè)縣級(jí)市，為云貴高原上僅有的兩個(gè)省會(huì)城市之一。常年受西風(fēng)帶控制，屬于亞熱帶濕潤(rùn)溫和型氣候，兼有高原性和季風(fēng)性氣候特點(diǎn)。

2.2 數(shù)據(jù)來(lái)源

本文選取的數(shù)據(jù)包括MODIS產(chǎn)品中的土地利用數(shù)據(jù)和MOD11A2地表溫度數(shù)據(jù)。其中，遙感數(shù)據(jù)為2000～2019年MOD11A2數(shù)據(jù)，分辨率1km，是已經(jīng)過(guò)定標(biāo)定位，大氣和氣溶膠校正的8天合成產(chǎn)品，此數(shù)據(jù)消除了部分的云干擾，提高了LST數(shù)據(jù)的有效性；土地利用數(shù)據(jù)為2006年分辨率為500m的土地利用覆被數(shù)據(jù)MCD12Q1。為保障數(shù)據(jù)的連續(xù)性和完整性MOD11A2和MCD12Q1均獲取自NASA官網(wǎng)(https://ladsweb.modaps.eosdis.nasa.gov/)。為便于季節(jié)變化分析，四季按照MOD11A2原始數(shù)據(jù)的時(shí)間編碼(001～361)進(jìn)行劃分，春季(3～5月份)為065～145，夏季(6～8月份)為153～241，秋季(9～11月份)為249～329，冬季(12月～次年2月份)為337～361和次年的001～057。

2.3 研究方法

利用MODIS產(chǎn)品中的土地利用數(shù)據(jù)和地表溫度數(shù)據(jù)，對(duì)貴陽(yáng)市地表熱場(chǎng)的時(shí)空特征進(jìn)行分析。在空間上，對(duì)城市的地表熱場(chǎng)強(qiáng)度等級(jí)進(jìn)行劃分，分析地表熱環(huán)境空間變化趨勢(shì)；時(shí)間上，分析城市地表熱場(chǎng)的年際變化特征。

采用MODIS地表溫度產(chǎn)品(LST)表示熱場(chǎng)分布，LST是反映熱場(chǎng)最直接、最重要的定量指標(biāo)[5]。由于白天貴陽(yáng)市地表溫度數(shù)據(jù)空值較多且不穩(wěn)定，而夜間的熱島強(qiáng)度較穩(wěn)定[6]，因此使用夜間地表溫度進(jìn)行研究。本文采用地表熱島強(qiáng)度指數(shù)[7](Urban Heat Island Intensity Index，UHII)和斜率分析法對(duì)研究區(qū)地表溫度的時(shí)空特征進(jìn)行研究。地表熱島強(qiáng)度指數(shù)是將地表熱島強(qiáng)度劃分為強(qiáng)冷島、較強(qiáng)冷島、弱冷島、無(wú)熱島、弱熱島、較強(qiáng)熱島和強(qiáng)熱島等7級(jí)，以反映城區(qū)溫度高于郊區(qū)鄉(xiāng)村的程度；斜率分析法能對(duì)不同的柵格變化趨勢(shì)進(jìn)行分析，得到地表溫度的年際變化特征。使用的軟件環(huán)境為ArcGIS 10.4、ENVI5.3、SPSS20等軟件。

3 結(jié)果與分析

3.1 貴陽(yáng)市地表熱場(chǎng)時(shí)空變化特征

貴陽(yáng)市作為典型喀斯特山地城市，熱島效應(yīng)最典型的特點(diǎn)是：大部分城區(qū)沒(méi)有熱島，而在沒(méi)有大面積集中建設(shè)和高耗能工礦企業(yè)的城郊區(qū)域卻出現(xiàn)了大面積的熱島，這與蔡宏等人[8]的研究結(jié)果一致。說(shuō)明貴陽(yáng)市存在區(qū)別于常規(guī)熱島的異常熱島。

圖1為基于MODIS數(shù)據(jù)得到的貴陽(yáng)市2005、2010、2015、2019年夜晚熱島強(qiáng)度等級(jí)分布圖。由圖可知，貴陽(yáng)市夜晚熱島強(qiáng)度等級(jí)僅包含7個(gè)等級(jí)中的4個(gè)等級(jí)(即：弱冷島、無(wú)熱島、弱熱島、較強(qiáng)熱島)。2005～2019年，貴陽(yáng)市的熱島效應(yīng)以弱熱島為主，主要分布在云巖區(qū)、南明區(qū)、花溪區(qū)及北部低洼河谷地帶，其次是較強(qiáng)熱島，較集中地分布于云巖區(qū)、南明區(qū)及開陽(yáng)縣東部的低海拔區(qū)域。弱熱島在觀山湖區(qū)、白云區(qū)、烏當(dāng)區(qū)僅有零星分布，而在云巖區(qū)、南明區(qū)及花溪區(qū)則呈現(xiàn)出擴(kuò)散趨勢(shì)。較強(qiáng)熱島未出現(xiàn)明顯擴(kuò)散趨勢(shì)，歷年均較集中地分布于云巖區(qū)、南明區(qū)及花溪區(qū)北部的低海拔區(qū)域。

圖1 2005～2019年貴陽(yáng)市熱島強(qiáng)度變化

統(tǒng)計(jì)2005～2019年貴陽(yáng)市夜晚熱島強(qiáng)度范圍平均地表面積占貴陽(yáng)市總面積的比例情況，發(fā)現(xiàn)貴陽(yáng)市82.27%的區(qū)域無(wú)熱島效應(yīng)，僅有17.73%的區(qū)域存在熱島現(xiàn)象，且以弱熱島為主，占貴陽(yáng)市總面積的17.41%，較強(qiáng)熱島僅占0.32%，說(shuō)明貴陽(yáng)市雖然存在熱島效應(yīng)，但并不強(qiáng)烈。

由圖2可見，貴陽(yáng)市地表溫度變化幅度均超過(guò)1.5 ℃/a，中部高海拔區(qū)域地表溫度變化幅度較四周低海拔區(qū)域低，大部分區(qū)域變化幅度均在1.5～1.8 ℃/a，而北部、西部及西南部城郊的大部分區(qū)域地表溫度變化幅度均超過(guò)1.8 ℃/a，尤其是南部的南明區(qū)和花溪區(qū)，部分區(qū)域變化幅度均超過(guò)2.0 ℃/a。從地表溫度平均變化幅度來(lái)看，周邊修文縣、開陽(yáng)縣、息烽縣及清鎮(zhèn)市平均變化幅度為1.78 ℃/a，云巖區(qū)、南明區(qū)及花溪區(qū)平均變化幅度為1.81 ℃/a，而中部觀山湖區(qū)、白云區(qū)及烏當(dāng)區(qū)平均變化幅度為1.72 ℃/a。

圖2 2000～2019年貴陽(yáng)市地表溫度年際變化率

3.2 貴陽(yáng)市地表熱場(chǎng)的時(shí)間變化特征

2000～2019年貴陽(yáng)市地表溫度年際變化呈增加趨勢(shì)，這與全球氣候變暖趨勢(shì)相一致，但變化率很小，僅為0.0617 ℃/a，增溫幅度不明顯(圖3)。分析不同季節(jié)地表溫度年際變化，發(fā)現(xiàn)貴陽(yáng)市地表溫度四季均有微小的增加趨勢(shì)，其中冬季的變化率最大，達(dá)0.0832 ℃/a，春季、夏季和秋季的變化率分別為0.0678 ℃/a、0.0336 ℃/a和0.0239 ℃/a。

圖3 2000～2019年貴陽(yáng)市地表溫度年際變化

2000～2019年，貴陽(yáng)市地表溫度季節(jié)差異明顯(圖4)，夏季地表溫度最高，冬季最低，夏季平均地表溫度達(dá)20.08 ℃，冬季平均地表溫度達(dá)7.18 ℃。夏季與冬季地表溫度平均溫差為12.90 ℃，變異系數(shù)達(dá)12.09%。春季和秋季的平均地表溫度較為接近，分別為15.90 ℃和15.10 ℃。四季中以冬季地表溫度的變異系數(shù)最大，達(dá)到19.11%，其次是春季和秋季地表溫度變異系數(shù)，分別為10.72%和8.33%，而夏季地表溫度波動(dòng)幅度最小，變異系數(shù)只有3.92%。冬季地表溫度波動(dòng)幅度最大，可能是由于氣候變化，貴陽(yáng)市冬季各年降雪量發(fā)生較大變化，導(dǎo)致下墊面反射和吸收太陽(yáng)輻射的能力出現(xiàn)明顯差異，引起冬季地表溫度波動(dòng)幅度高于其他季節(jié)。

圖4 2000～2019年貴陽(yáng)市四季地表溫度年際變化

4 結(jié)論

2000～2019年貴陽(yáng)市地表溫度年際變化呈增加趨勢(shì)，這與全球氣候變暖趨勢(shì)相一致，但變化率很小，僅為0.0617 ℃/a，增溫幅度不明顯；貴陽(yáng)市雖然存在熱島效應(yīng)，但并不強(qiáng)烈，大部分城區(qū)沒(méi)有熱島，而在沒(méi)有大面積集中建設(shè)和高耗能工礦企業(yè)的城郊區(qū)域卻出現(xiàn)了大面積的熱島；貴陽(yáng)市的熱島效應(yīng)以弱熱島為主，較強(qiáng)熱島僅有少量分布；地表溫度季節(jié)差異明顯，夏季地表溫度最高，冬季最低，平均溫差為12.90 ℃，春季和秋季的平均地表溫度較為接近，平均溫差小于1 ℃。本文使用的MOD11A2數(shù)據(jù)，雖然時(shí)間分辨率較高，但由于其空間分辨率較低，未能針對(duì)不同地類及地表環(huán)境與地表溫度間的關(guān)系進(jìn)行詳細(xì)分析，未來(lái)還需結(jié)合精度更高的數(shù)據(jù)研究下墊面與地表溫度間的相互影響過(guò)程。