蘇俊如，呂 芳，楊 俊,2

(1. 遼寧師范大學(xué)自然地理與空間信息科學(xué)遼寧省重點實驗室，遼寧 大連 116029； 2. 中國科學(xué)院地理科學(xué)與資源研究所，北京 100101)

城市的高速發(fā)展致使城市空間急速增長，引起了一系列的城市問題。城市熱島對城市人居環(huán)境、城市微氣候變化、空氣質(zhì)量和公共健康等方面均產(chǎn)生了深遠影響，同時還會對其他城市氣候效應(yīng)產(chǎn)生影響，如大氣污染擴展、降水空間分布和植被物候期變化等[1]。

近年來，城市熱環(huán)境效應(yīng)引起了國內(nèi)學(xué)者的高度關(guān)注，通常用單窗算法和劈窗算法反演地表溫度。覃志豪等據(jù)地表熱輻射傳導(dǎo)方程，推導(dǎo)出一個簡單易行并且精度較高的演算方法——單窗算法。該算法通過熱紅外波段數(shù)據(jù)提取亮度溫度同時結(jié)合地表輻射率、大氣透射率和大氣平均溫度3個參數(shù)反演地表溫度[2-3]。劈窗算法通常利用2個熱紅外波段，最初應(yīng)用于推算海洋表面溫度，自20世紀80年代中期以來，劈窗算法逐漸應(yīng)用于地表溫度的反演研究[4-5]。

在城市熱環(huán)境效應(yīng)的影響研究中，熊淑娣等通過劈窗算法反演地表溫度，分析了綠地和不同下墊面在城區(qū)與郊區(qū)中對地表溫度的作用[6]；龔珍等通過單窗算法對HJ-1B衛(wèi)星反演地表溫度，對水體和地表溫度進行了相關(guān)分析，提出地表平均溫度受到研究區(qū)域內(nèi)的水體面積比重的影響[7]；孫茂存和李俊峰運用亮度溫度計算和二分模型計算等方法，分析了造成烏魯木齊市城市熱島空間分布的原因，提出熱島中心主要分布在城區(qū)人口密度大及工業(yè)區(qū)聚集的地方，植被和水體對減弱熱島效應(yīng)起到一定作用[8]；陳紹杰等通過建筑用地、土地覆蓋圖和地表熱環(huán)境信息，分析了研究區(qū)城市擴展動態(tài)，同時發(fā)現(xiàn)城市熱島現(xiàn)象顯著，建筑用地不斷增加是熱島效應(yīng)加重的主要因素[9]；Feng Yang等以上海市3個高層住宅小區(qū)為例基于建筑布局、密度和綠化等因子運用回歸分析分析了上海市城市熱島強度，發(fā)現(xiàn)布局、密度和綠地對城市熱島和熱島日夜間模式有不同的影響[10]；Stewart等將建筑密度作為局部氣候帶的一個影響因子，對局部氣候帶的熱環(huán)境效應(yīng)進行了研究[11]；Feng和Myint通過研究不同的土地覆蓋特征和建筑，探索在其鄰近地區(qū)的中心建筑對象和土地覆蓋格局的地表溫度之間的關(guān)系[12]；葛亞寧等通過目視解譯對建筑密度進行劃分，結(jié)合地表溫度反演數(shù)據(jù)，對城市建筑密度分布與城市熱島效應(yīng)進行了分析研究[13]。

綜上所述，城市建筑會對城市熱環(huán)境效應(yīng)產(chǎn)生一定影響。本文研究選用Landsat8和Landsat5遙感影像數(shù)據(jù)，對Landsat8的熱紅外Band10波段和Landsat5的熱紅外Band6波段進行地表溫度反演，并結(jié)合建筑密度因子，依據(jù)雙變量相關(guān)模型分析建筑密度對地表溫度的影響，通過分析研究區(qū)內(nèi)建筑水平特征與地表溫度之間的作用關(guān)系，使城市熱環(huán)境效應(yīng)的研究由理論研究轉(zhuǎn)向解決城市發(fā)展實際問題的研究。

1 研究區(qū)域與數(shù)據(jù)來源

1.1 研究區(qū)概況

中山區(qū)為遼寧省大連市的一個市轄區(qū)，位于大連市區(qū)東南部，東、南、北臨海，西與西崗區(qū)接壤。本文以中山區(qū)現(xiàn)有行政陸地界線為研究范圍，如圖1所示。

圖1 研究區(qū)示意圖

1.2 數(shù)據(jù)來源

原始數(shù)據(jù)為Landsat8 OLI、TIRS數(shù)據(jù)和Landsat5 TM遙感數(shù)據(jù)影像(見表1)，OLI、TIRS、TM影像與同期MODIS大氣溫度、濕度廓線產(chǎn)品(MOD07)產(chǎn)品參照空間插值獲取研究區(qū)溫度和濕度數(shù)據(jù)，將同期 MOD07數(shù)據(jù)得到的氣溫、相對濕度數(shù)據(jù)與氣象觀測點數(shù)據(jù)進行空間插值，獲取溫度及相對濕度空間數(shù)據(jù)[14]。

表1 數(shù)據(jù)來源與說明

2 研究方法與實例分析

2.1 單窗算法

本文選取覃志豪提出的單窗算法[3]，通過熱紅外波段數(shù)據(jù)提取亮度溫度，同時結(jié)合地表比輻射率、大氣平均溫度和大氣透射率進行地表溫度的演算，公式如下

Ts={a(1-C-D)+[b(1-C-D)+C+D]T-DTa}/C

(1)

C=τε

(2)

D=(1-τ)[(1-ε)τ]

(3)

式中，a和b為根據(jù)熱輻射強度和亮度溫度的關(guān)系擬合出的系數(shù)，a=-67.355 351，b=0.458 606；ε為地表比輻射率；τ為整層大氣透射率；T為傳感器的亮度溫度，單位為K；Ta為大氣平均作用溫度，單位為K。本文將同期 MOD07數(shù)據(jù)得到的氣溫、相對濕度數(shù)據(jù)與氣象觀測點數(shù)據(jù)進行空間插值獲取溫度及相對濕度空間數(shù)據(jù)。

通過對地表溫度進行網(wǎng)格化，研究區(qū)共分為54 291個網(wǎng)格，如圖2所示。

圖2 地表溫度空間分布

從地表溫度空間分布的演變過程來看(如圖2所示)，2002—2014年研究區(qū)主體溫度由294～306 K向306～314 K轉(zhuǎn)移，研究區(qū)溫度整體都有所升高。中山區(qū)東部由于植被覆蓋較多建筑較少，溫度相對較低。中山區(qū)的西北部(人民路街道、青泥洼橋街道、昆明街道、桂林街道及海軍廣場街道)及中部(解放路和中南路的沿線)多為建筑區(qū)域，溫度相對較高，而且建筑密度越大，溫度相對越高。2002—2014年，中山區(qū)東北部(東港街道)隨著填海造陸的進行，溫度明顯升高。

2.2 建筑密度

建筑密度是指在一定空間范圍內(nèi)所有建筑物基底面積的總和與占用地面積的百分比[11，15]。建筑密度是衡量一個區(qū)域中建筑空間分布的重要指標，對建筑密度進行分類，見表2。

表2 建筑密度分類

從研究區(qū)建筑密度分布的演變過程(如圖3 所示)來看，高和較高密度的地塊主要分布在青泥洼天津街商業(yè)區(qū)(包括青泥和人民路街道與昆明街道的東北地區(qū))。大多數(shù)低、較低、中密度地塊分布沿解放路和中南路沿線。在2002—2014年，高和相對高密度的地塊在桂林和海軍廣場街道變成了中等和較低密度的地塊；同時，在研究區(qū)中部建筑密度明顯增加。

圖3 建筑密度分布

2.3 雙變量相關(guān)分析

雙變量相關(guān)分析主要用于分析2個或多個變量間的相關(guān)性。本文研究選擇Pearson 簡單相關(guān)系數(shù)，其計算公式為

(4)

Pearson簡單相關(guān)系數(shù)的對應(yīng)檢驗統(tǒng)計量為t統(tǒng)計量，其計算公式為

(5)

SPSS將根據(jù)t統(tǒng)計量和自由度，依照t分布表自動給出t統(tǒng)計量所對應(yīng)的相伴概率。

通過對2002、2014年的地表溫度數(shù)據(jù)和建筑密度進行相關(guān)分析，得到以下結(jié)果，見表3。

表3 2002—2014年地表溫度與建筑密度相關(guān)性

注：**.表示在0.01水平(雙側(cè))上顯著相關(guān)。

總體上，地表溫度和建筑密度呈現(xiàn)中等相關(guān)。2002、2014年，地表溫度與建筑密度的相關(guān)系數(shù)分別為0.514、0.537，呈中等相關(guān)。

3 結(jié)論與討論

本文研究以大連市中山區(qū)為例，通過單窗算法進行地表溫度反演，將地表溫度與建筑密度進行雙變量相關(guān)分析，在一定程度上探討了地表溫度與建筑密度間的相關(guān)性，研究結(jié)果如下：

(1) 2002—2014年研究區(qū)主體溫度由294～306 K向306～314 K轉(zhuǎn)移，研究區(qū)溫度整體都有所升高。中山區(qū)的西北部(人民路街道、青泥洼橋街道、昆明街道、桂林街道及海軍廣場街道)及中部(解放路和中南路的沿線)溫度相對較高；東北部(東港街道)由于填海造陸，溫度明顯升高。

(2) 高密度和較高密度的地塊主要分布在青泥洼天津街商業(yè)區(qū)(包括青泥和人民路街道和昆明街道的東北地區(qū))；大部分低、較低、中密度地塊分布在解放路和中南路沿線。

(3) 2002、2014年，地表溫度與建筑密度呈現(xiàn)中等相關(guān)性，相關(guān)系數(shù)分別為0.514、0.537。

城市熱效應(yīng)受建筑、道路、綠地分布等多方面因素影響。而本文對于城市熱效應(yīng)的描述只側(cè)重于建筑密度，未考慮城市其他因素影響，同時兩期研究數(shù)據(jù)時間有一定偏差，反演溫度可能存在一定誤差，城市熱環(huán)境效應(yīng)機制還需要進一步研究和總結(jié)。