左青 陶前輝 李博 楊波

摘 要：以常德市規(guī)劃區(qū)為研究區(qū)，利用2003年、2008年、2013年、2017年及2018年5景遙感影像數(shù)據(jù)，提取植被覆蓋度及建筑指數(shù)，反演其地表溫度，并評(píng)價(jià)熱島強(qiáng)度。結(jié)果表明，城市地表溫度與建筑指數(shù)、植被覆蓋度具有相關(guān)性，分別呈正相關(guān)和負(fù)相關(guān);常德市熱島效應(yīng)顯著，2003—2013年熱島強(qiáng)度比例指數(shù)由0.222升高為0.245，熱島面積增加14.035km2;2013—2018年熱島比例指數(shù)下降為0.221，熱島面積減少19.466km2;在2013—2018年期間，海綿城市的試點(diǎn)建設(shè)區(qū)域熱島格局有破碎狀趨勢(shì)，說(shuō)明建設(shè)海綿城市對(duì)熱島效應(yīng)具有一定的緩解作用。

關(guān)鍵詞：植被覆蓋度;歸一化建筑指數(shù);城市熱島強(qiáng)度;海綿城市;常德市

中圖分類(lèi)號(hào) X16 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 1007-7731（2020）05-0136-05

Change and Analysis of Heat Island Effect in Changde City in Recent 15 Years

Zuo Qing et al.

（College of Resources and Environmental Sciences/Key Laboratory of Geospatial Big Data Mining and Application， Hunan Normal University， Changsha 410081， China）

Abstract： Taking Changde city Planning Area as the research area， vegetation coverage and building index were extracted from five remote sensing images in 2003， 2008， 2013， 2017 and 2018， to retrieve the land surface temperature and evaluate the urban heat island intensity. The results show that there is a correlation between urban land surface temperature and building index and vegetation coverage， with positive correlation and negative correlation respectively; Changde City has a significant urban heat island effect. From 2003 to 2013， the urban-heat-island ratio index increased from 0.222 to 0.245， and the area of heat island increased by 14.035 km2; From 2013 to 2018， urban-heat-island ratio index decreased to 0.221， and the area of heat island decreased by 19.466 km2.During the period from 2013 to 2018， the pattern of urban heat island in the pilot construction area of sponge city has a trend of fragmentation， which indicates that the construction of sponge city has a certain mitigation effect on the urban heat island effect.

Key words： Vegetation Coverages; Normalized Difference Build-up Index （NDBI）; Urban Heat Island Intensity （UHII）; Sponge City; Changde City

1 引言

隨著社會(huì)經(jīng)濟(jì)飛速發(fā)展，城市人口迅猛增長(zhǎng)，城市化進(jìn)程的不斷加快，城市下墊面受到人類(lèi)活動(dòng)的干擾增強(qiáng)，其結(jié)構(gòu)發(fā)生改變，城市生態(tài)對(duì)熱環(huán)境的調(diào)節(jié)能力已經(jīng)超過(guò)其閾值，熱島效應(yīng)（Urban Heat Island，UHI）給城市系統(tǒng)帶來(lái)的氣候、生態(tài)環(huán)境及人類(lèi)健康等問(wèn)題也愈發(fā)嚴(yán)峻[1-4]。

目前，國(guó)內(nèi)外學(xué)者從宏觀和微觀角度對(duì)城市熱島的演變格局、影響因子及緩解策略等方面進(jìn)行了大量研究[5-9]。但國(guó)內(nèi)外學(xué)者對(duì)熱島效應(yīng)的研究主要聚焦在北京、上海、廣州、歐洲城市群等大城市[10-13]，較少涉及對(duì)中小城市的熱島效應(yīng)的研究。本研究以中小型城市常德市為對(duì)象，該市近年來(lái)城市化進(jìn)程加快，城市熱島效應(yīng)明顯，研究常德市近15年夏季（6—8月）城市熱島的分布及變化情況，分析城市熱島與植被覆蓋度及建筑指數(shù)的關(guān)系，從而緩解常德市的城市化發(fā)展過(guò)程中引起的熱島效應(yīng)，為其它中小城市的發(fā)展合理規(guī)劃景觀格局及改善生態(tài)環(huán)境提供理論基礎(chǔ)與決策依據(jù)。

2 研究區(qū)與數(shù)據(jù)源

2.1 研究區(qū)概況 常德市（110°27′33″E—112°17′52″E， 28°23′31″N—30°07′53″N）位于湖南省西北部，大體處于西洞庭湖區(qū)，屬亞熱帶季風(fēng)濕潤(rùn)氣候區(qū)，四季分明，春秋短，夏冬長(zhǎng)。常德市作為長(zhǎng)江經(jīng)濟(jì)帶、環(huán)洞庭湖生態(tài)經(jīng)濟(jì)圈的重要城市，其經(jīng)濟(jì)實(shí)力穩(wěn)居湖南省第3，僅次于長(zhǎng)沙市與岳陽(yáng)市。本研究選取常德市城市規(guī)劃區(qū)作為研究區(qū)，包括長(zhǎng)庚街道、白馬湖街道、穿紫河街道、丹陽(yáng)街道等20個(gè)街道（鄉(xiāng)/鎮(zhèn)），總面積613.6km2（圖1）。

2.2 數(shù)據(jù)源及其預(yù)處理 選取近15年5景夏季（6—9月）的Landsat5-TM （2003/07/27、2008/06/06）及Landsat8-OLI（2013/08/07、2017/08/18、2018/07/20）的遙感影像資料作為反演研究區(qū)LST數(shù)據(jù)源。選取的5景影像云層覆蓋較少，在研究區(qū)內(nèi)幾乎無(wú)云層覆蓋，可以確保反演的精度。5景Landsat影像數(shù)據(jù)均從美國(guó)地質(zhì)調(diào)查局（USGS）官網(wǎng)（http：//glovis.usgs.gov）獲取，其條帶號(hào)為124/040。

LST反演之前均在ENVI5.3中進(jìn)行幾何校正、輻射定標(biāo)、大氣校正等圖像預(yù)處理工作，以消除成像過(guò)程中的幾何畸變、噪音等影響，然后對(duì)研究區(qū)進(jìn)行掩膜裁剪。

3 研究方法

本研究對(duì)2003年、2008年、2013年、2017年及2018年5景Landsat TM/OLI 數(shù)據(jù)，采用仿歸一化指數(shù)法、混合像元法、輻射傳輸法分別提取了NDBI、植被覆蓋度以及地表溫度，同時(shí)采用均值—標(biāo)準(zhǔn)差法劃分熱島強(qiáng)度等級(jí)，計(jì)算熱島強(qiáng)度比例指數(shù)，常德城區(qū)溫度變化趨勢(shì)，探討NDBI和植被覆蓋度對(duì)地表溫度和熱島強(qiáng)度的影響。

3.1 植被覆蓋度提取

3.1.1 NDVI計(jì)算 NDVI為歸一化植被指數(shù)，本研究采用非線性歸一化處理得到NDVI：

[NDVI=DNNIR-DNRDNNIR+DNR] （1）

式中，[DNNIR]表示近紅外波段;[DNR]表示熱紅外波段。

3.1.2 植被覆蓋度計(jì)算 植被覆蓋度（[FV]），可用以表征地表植被覆蓋情況。計(jì)算公式如下：

[FV=NDVI-NDVISNDVIV-NDVIS] （2）

式中，[NDVI]為歸一化植被指數(shù);[NDVIV]、[NDVIS]分別為純植被和純土壤的植被指數(shù);取[NDVIV=0.7]和[NDVIS=0]，且[FV]的值域?yàn)閇[0，1]]，當(dāng)影像中某個(gè)像元的[NDVI]大于0.7時(shí)，[FV]取值為1;當(dāng)[NDVI]小于0，[FV]取值為0。

3.2 NDBI提取 歸一化建筑指數(shù)（NDBI）是基于歸一化植被指數(shù)提出：

[NDBI=DNMIR-DNNIRDNMIR+DNNIR] （3）

式中，[DNMIR]表示中紅外波段;[DNNIR]表示近紅外波段。

3.3 LST反演 本研究采用熱輻射傳輸法對(duì)Landsat TM/OLI數(shù)據(jù)進(jìn)行LST反演[14]，基于輻射強(qiáng)度反演地表溫度。

3.3.1 地表比輻射率計(jì)算 比輻射率[ε]是用以研究城市地表熱環(huán)境的一個(gè)重要參數(shù)，把地物類(lèi)型分為水體、城鎮(zhèn)和自然表面3類(lèi)并給水體的[ε]賦值為0.995。 通過(guò)以下公式分別對(duì)城鎮(zhèn)的比輻射率（[εsurface]）和自然表面的比輻射率（[εbuilding]）進(jìn)行計(jì)算：

[εsurface=0.9625+0.0614FV-0.0461F2V] （4）

[εbuilding=0.9589+0.086FV-0.0671F2V] （5）

式中，[FV]表示植被覆蓋度。

3.3.2 輻射亮度值計(jì)算 根據(jù)輻射傳輸方程，溫度為[T]的黑體在熱紅外波段的輻射亮度[B（TS）]為：

[B（TS）=[Lλ-L↑-τ?（1-ε）L↓]/τ?ε] （6）

式中，[ε]為地表輻射率;[TS]為地表真實(shí)溫度;[B（TS）]為黑體在[TS]的熱輻射亮度;[τ]為大氣在熱紅外波段的透過(guò)率;[L↑]為大氣向上輻射亮度，單位[W/（m2?sr?μm）];[L↓]為大氣向下輻射到達(dá)地面后反射的能量，單位[W/（m2?sr?μm）]。Landsat影像的大氣剖面參數(shù)（[τ]、[L↑]、[L↓]）可以通過(guò)NASA提供的網(wǎng)站（http：//atmcorr.gsfc.nasa.gov/）直接獲取。

3.3.3 地表溫度計(jì)算 根據(jù)普朗克公式的反函數(shù)，求得地表真實(shí)溫度[TS]：

[TS=K2/ln（K1/B（TS）+1）] （7）

對(duì)于Landsat5 TM數(shù)據(jù)，[K1][=][607.76 W/（m2][?sr][?μm）]，[K2=1260.56? K];對(duì)于Landsat8 TIRS數(shù)據(jù)，[K1][=][774.89W/（m2][?sr][?μm）]，[K2=1321.08? K]。

3.4 城市熱島強(qiáng)度等級(jí)評(píng)估與計(jì)算

3.4.1 UHII等級(jí)評(píng)估 在實(shí)際工作中，根據(jù)均值-標(biāo)準(zhǔn)差法評(píng)價(jià)熱島強(qiáng)度等級(jí)時(shí)，一般分為4～7級(jí)[15]。本研究參考將熱島強(qiáng)度劃分為5級(jí)：強(qiáng)冷島區(qū)、冷島區(qū)、無(wú)熱島、熱島區(qū)及強(qiáng)熱島區(qū)（見(jiàn)表1）。

3.4.2 URI計(jì)算 城市熱島比例指數(shù)（URI）由徐涵秋[16]提出，計(jì)算過(guò)程如下：

[URI=1100mi=1nWiPi] （8）

式中，m為熱島強(qiáng)度等級(jí)總數(shù)量，在此取值為5;i為熱島等級(jí)序號(hào);n為熱島等級(jí)數(shù)量;Wi為熱島的權(quán)重值，取第i級(jí)的級(jí)數(shù);Pi為第i級(jí)熱島面積百分比。一般認(rèn)為，地表溫度高于平均溫度的區(qū)域?qū)儆跓釐u區(qū)域，指數(shù)值越大，熱島現(xiàn)象越嚴(yán)重。

4 結(jié)果與分析

4.1 地表溫度分布及變化 以常德市中心城區(qū)為研究區(qū)，利用Landsat遙感影像數(shù)據(jù)，通過(guò)熱紅外波段反演得到各年地表溫度（圖2），并統(tǒng)計(jì)各年平均溫度（表2）。從圖2可以看出，在2003—2013年期間研究區(qū)溫度整體呈上升趨勢(shì)，在2013—2018年期間研究區(qū)溫度整體有所下降。2013年受到副熱帶高氣壓帶的影響出現(xiàn)異常高溫現(xiàn)象[17-18]。同時(shí)，沅江北岸的長(zhǎng)庚街道、白馬湖街道、穿紫河街道、丹陽(yáng)街道、府坪街道、啟明街道及灌溪鎮(zhèn)中心高于白鶴鎮(zhèn)、蘆荻山鄉(xiāng)、東江街道、丹洲鄉(xiāng)及灌溪鎮(zhèn)外圍地區(qū)的溫度，沅江南岸的武陵鎮(zhèn)東北部、德山街道及樟木橋街道北部的溫度高于斗姆湖街道、樟木橋街道及石門(mén)橋鎮(zhèn)的溫度;即NDBI對(duì)地表溫度的影響高于植被覆蓋度對(duì)地表溫度的影響。以常德市西北部產(chǎn)業(yè)新城灌溪鎮(zhèn)為例，灌溪鎮(zhèn)以機(jī)械制造業(yè)作為主產(chǎn)業(yè)發(fā)展城市化，土地利用類(lèi)型以工業(yè)用地為主，呈現(xiàn)出1個(gè)高溫中心。而水體在整個(gè)研究區(qū)呈現(xiàn)明顯的低溫狀態(tài)，如柳葉湖和沅江，呈現(xiàn)出低溫中心。根據(jù)研究區(qū)地表溫度反演統(tǒng)計(jì)（表2），2003年、2008年平均溫度在30℃左右，到2013年平均溫度達(dá)45℃，溫度呈上升趨勢(shì)，上升近15℃;2017年和2018年平均溫度在35℃左右，溫度呈下降趨勢(shì)，同2013年相比，下降近10℃。