馮徽徽,劉慧平,朱壽東,白 穆

(北京師范大學(xué)地理學(xué)與遙感科學(xué)學(xué)院遙感科學(xué)國家重點實驗室,北京 100875）

城市擴張對東莞市陸表溫度的影響分析

馮徽徽,劉慧平＊,朱壽東,白 穆

(北京師范大學(xué)地理學(xué)與遙感科學(xué)學(xué)院遙感科學(xué)國家重點實驗室,北京 100875）

利用1995年12月30日、2003年1月10日兩期Landsat TM/ETM遙感影像反演東莞市陸表溫度(LST）,分析東莞市陸表溫度的總體情況及時空特征,并通過引入溫度植被指數(shù)(Temperature-Vegetation Index,TVX）空間方法,研究城市化對陸表溫度的影響。結(jié)果表明:城市化進程是導(dǎo)致東莞市陸表溫度上升的重要原因,這種作用在TVX空間中主要體現(xiàn)在聚類點由高植被覆蓋-低溫度向低植被覆蓋-高溫度區(qū)域的運動過程。

城市擴張;陸表溫度;時空特征;溫度植被指數(shù)(TVX）;東莞

0 引言

改革開放以來,隨著經(jīng)濟的迅速發(fā)展,東莞市經(jīng)歷了由“三來一補”帶動的“自下而上”的城市化過程[1-3],城市用地在短時期內(nèi)迅速擴張?？焖俪鞘谢M程不僅深刻影響東莞市社會經(jīng)濟變革,同時也對其資源環(huán)境產(chǎn)生了深刻影響。作為地球資源環(huán)境動態(tài)分析的重要指標,陸表溫度(Land Surface Temperature,LST）對城市規(guī)劃、環(huán)境保護等具有重要作用,它受到許多因素的影響,如季節(jié)因素、下墊面特征、土地利用/覆蓋類型等[4-7]。

近年來,土地利用變化尤其是城市化進程對陸表溫度的影響已成為學(xué)者關(guān)注的熱點問題,如Xiao等[6]通過分析中國喀斯特地區(qū)陸表溫度的分布情況,研究了城鎮(zhèn)建筑材料與地表植被對陸表溫度的影響;宮阿都等[8,9]通過陸地衛(wèi)星熱紅外單波段分別反演了北京與珠三角地區(qū)的陸表溫度,研究了土地覆蓋變化特別是城市增長對地表溫度的影響;劉宇等[10]通過分析不同時期東莞市陸表溫度的分布情況,構(gòu)建了熱效應(yīng)貢獻度指數(shù)等3個評價指標,定量分析了不同土地利用類型對城市陸表溫度的貢獻度。這些研究從不同的角度探討了城市化過程對地表溫度的影響,但在定量分析及可視化表達方法的研究上有待進一步深化。

本文在總結(jié)前人研究成果的基礎(chǔ)上,利用多時相遙感數(shù)據(jù)反演了東莞市陸表溫度,并引入溫度植被指數(shù)(TVX）空間方法,通過觀察 TVX空間中聚類點的運動軌跡,分析了城市化進程對陸表溫度動態(tài)變化的影響。

1 研究區(qū)與數(shù)據(jù)處理

本文選取廣東省東莞市(東經(jīng)113°31′～114°15′,北緯22°39′～23°09′）為研究區(qū),區(qū)域東西長約70.45 km,南北寬約46.8 km,陸地面積2 465 km2[11]。改革開放以來,由于國家政策的引導(dǎo),東莞市經(jīng)濟迅速發(fā)展,城市化水平迅速提高,城市在短期內(nèi)占用了大量的農(nóng)業(yè)用地,得到迅速擴張(圖1）。

圖1 1995-2003年東莞市城市擴張Fig.1 Urban growth in Dongguan City from 1995 to 2003

主要數(shù)據(jù)包括1995年12月30日、2003年1月10日兩期Landsat TM/ETM影像,軌道號是122/ 44,成像時間為北京時間上午10時到11時,輔助數(shù)據(jù)包括東莞市行政區(qū)劃圖、2000年東莞市土地利用分類矢量圖(作為參照圖輔助進行其他年份的土地利用分類）。先對兩期影像進行幾何精校正,匹配精度控制在0.5個像元內(nèi),并將所有數(shù)據(jù)統(tǒng)一采樣成30 m的分辨率;根據(jù)《土地利用現(xiàn)狀分類》國家標準(2007年8月10日頒布）及《城鎮(zhèn)地籍調(diào)查規(guī)程》(1989年9月6日頒布）,結(jié)合研究區(qū)的實際情況,通過在圖像上選取一系列樣本點,運用CART算法獲取土地利用分類規(guī)則,在ENV I中構(gòu)建決策樹,將研究區(qū)土地利用類型分成農(nóng)用地、林地、水域、城鎮(zhèn)用地及其它用地五大類;最后以東莞行政區(qū)劃圖生成掩膜,對分類圖像進行裁剪,得到只包含東莞市的土地利用分類圖。

2 研究方法

2.1 陸表溫度LST計算方法

利用遙感數(shù)據(jù)進行陸表溫度反演已有許多成熟的算法,如分裂窗算法[12,13]、單窗算法等[14-16]。對于只有一個熱紅外波段的 TM數(shù)據(jù),一般采用單窗算法反演陸表溫度。本文由于未能獲取到近地表氣溫和大氣水分含量的大氣參數(shù),故無法使用單窗算法反演地表溫度;另一方面,華南地區(qū)的大氣水汽含量高,使用覃志豪等[16]提出的方法估算大氣透射率會存在較大的誤差。故本文采用A rtis及 Sobrino等提出的溫度反演方法[17-19],在天空晴朗的條件下,大氣的影響程度在空間上可認為一致,其不會影響陸表溫度的空間相對分布。具體計算過程如下:

(1）輻射亮度Lb的計算:

式中:Lmax、Lmin分別表示最大和最小熱輻射能量,可從原始影像的頭文件中獲取。

式中:K1、K2為常數(shù)。

式中:λ=11.5μm,α=0.01438 m K,ε表示地表輻射率。ε可參考Sobrino等[18,19]提出的方法進行估計:當NDV I≥0.7時,ε=0.99;NDV I≤0.05時,ε=0.973; 0.05＜NDV I＜0.7時,ε=0.986+0.004×Pv(Pv表示植被覆蓋度,Pv=[(NDV I-0.05）/0.65]2）。

2.2 溫度植被指數(shù)(TVX）空間軌跡分析方法

自1972年 Rao[5]首次提出利用衛(wèi)星影像測定城區(qū)的熱運動軌跡以來,越來越多的學(xué)者致力于研究土地利用變化(尤其是城市化）、植被指數(shù)等對陸表溫度的影響[20-23]。本文引入一種研究陸表溫度時空變化的新方法——溫度植被指數(shù)(TVX）空間方法[5,22],在集成陸表溫度與植被指數(shù)的散點圖空間中研究兩者的相互關(guān)系:即以陸表溫度與植被指數(shù)分別作為TVX空間的橫、縱坐標,構(gòu)建 TVX空間;以某類土地利用類型為聚類點,計算其在空間中的位置。通過分析TVX空間中聚類點的運動軌跡,可以研究陸表溫度的時空變化特征,尤其是土地利用變化(城市化）對陸表溫度變化的影響。

為使不同時相的數(shù)據(jù)有更好的可比性,需對陸表溫度進行歸一化處理,同時以植被覆蓋度作為植被指數(shù)的分析指標,具體計算公式為:

式中:T表示像元平均溫度,Tmin、Tmax分別表示研究區(qū)中溫度最低值和溫度最高值。

式中:NDV I表示像元的平均NDV I值,NDV Imin表示裸土的 NDV I值,NDV Imax表示全植被覆蓋下的NDV I值。本研究采用 Sobrino等[18,19]的估計方法,NDV Imin取值0.05,NDV Imax取值0.7。

3 結(jié)果分析

3.1 陸表溫度LST反演結(jié)果

根據(jù)式(1）～式(3）,在 IDL編程環(huán)境下,對Landsat TM/ETM影像數(shù)據(jù)進行計算、處理,反演出多時相不同季節(jié)的東莞市陸表溫度,并根據(jù)式(4）求出各時期的歸一化溫度,使各時相數(shù)據(jù)間具有更好的可比性。最終在A rcGIS9.3中根據(jù)自然斷裂法對陸表溫度進行等級分類,將其分為低溫、亞低溫、中溫、亞高溫和高溫5個類別。

(1）1995—2003年研究區(qū)的陸表平均溫度上升幅度較大(表1）,圖2(見封2）顯示了兩個時期的陸表溫度分布情況。對比東莞市土地利用圖與圖2的顯示結(jié)果可以發(fā)現(xiàn),高溫區(qū)域主要集中在城鎮(zhèn)地區(qū),在空間上呈現(xiàn)以城鎮(zhèn)為中心、沿道路不斷向兩旁延展的線狀格局,低溫區(qū)域主要集中在林地。表明城市化過程確實導(dǎo)致研究區(qū)陸表溫度上升,主要原因是由于城市建筑材料(如水泥、磚瓦、瀝青、建筑等）熱容量和熱慣量小,而熱傳導(dǎo)率和熱擴散率卻較大,使得接收的太陽輻射很快輻射到周圍大氣中,從而導(dǎo)致城區(qū)的陸表溫度高于周邊地區(qū),引發(fā)城市熱島效應(yīng),提高研究區(qū)陸表總體平均溫度。

表1 1995-2003年東莞市歸一化陸表溫度 T＊S總體情況Table 1 The overall condition of Normal LST(T＊S） of Dongguan from 1995 to 2003

(2）以土地利用分類圖為掩膜,計算不同土地利用類型的陸表溫度分布狀況(表2）,結(jié)果表明,城鎮(zhèn)用地、其它用地(主要為未利用地）的陸表溫度一般較高,而農(nóng)用地、林地、水體相對較低。值得注意的是,1995年的水域平均溫度均大于農(nóng)用地和林地,之后水域溫度逐漸降低,成為各土地利用類型中溫度最低區(qū)域,主要原因可能是1990-1995年東莞市水污染較嚴重,制造業(yè)污水直接排放到河流中[6],導(dǎo)致水體熱容量降低,溫度升高較快。而到后期,隨著東莞市加大對環(huán)境污染的治理力度,水污染情況得到極大改善,水域在溫度序列中逐漸成為最低的區(qū)域。另需關(guān)注的是,各土地利用類型的溫度標準差均降低較大,表明隨著東莞市加大對土地資源的管理力度,土地利用方式逐漸趨于合理。早期土地利用類型的溫度標準差均較大,主要是由于1990年以后的開發(fā)熱導(dǎo)致大量農(nóng)用地和林地被占用,1995年由于投資鏈斷裂而出現(xiàn)了大批土地撂荒現(xiàn)象,導(dǎo)致研究區(qū)內(nèi)混合像元數(shù)量增大,從而使得其地表輻射溫度具有較寬的變化范圍。水域的溫度標準差逐年遞減,從另一方面表明因水體污染造成水域的不均一性得到改善,水域溫度的變化范圍逐漸減小。

表2 不同土地利用類型的歸一化陸表溫度 T＊S分布狀況Table 2 The distribution of Normal LST(T＊S）of Dongguan with the various land use type

3.2 東莞市陸表溫度LST時空變化特征分析

根據(jù)各個時期的土地利用圖,提取在下一階段轉(zhuǎn)變?yōu)槌鞘杏玫氐耐恋乩妙愋?同時根據(jù)式(4）和式(5）,計算其變化前后的植被覆蓋度與歸一化陸表溫度,進而構(gòu)建TVX空間進行分析。由于東莞市城市化過程中占用最多的是農(nóng)用地和林地,故本文主要針對這兩種土地利用類型的變化對陸表溫度的影響進行分析,得到TVX空間變化軌跡(圖3）。

圖3 TVX空間變化軌跡Fig.3 The change trajectory of TVX space

(1）TVX空間中聚類點的運動軌跡揭示了土地利用變化(城市化）對陸表溫度變化的影響。在城市化早期,林地和農(nóng)用地處于 TVX空間的左上角,該區(qū)域?qū)?yīng)著較高的植被覆蓋度和較低的溫度。隨著城市用地不斷占用林地與農(nóng)用地,被占用的土地類型均沿其斜線軌跡向TVX空間的右下角運動,該區(qū)域?qū)?yīng)較低的植被覆蓋度和較高的溫度。聚類點在TVX空間中的初始位置決定了其運動變化的方向與幅度。圖3表明,在農(nóng)用地與林地轉(zhuǎn)變?yōu)槌擎?zhèn)用地的過程中,陸表溫度與植被覆蓋度的變化情況是,農(nóng)用地:ΔT′farmland=0.107,ΔP′vfarmland= -0.025,林地: ΔT′woodland=0.105,ΔP′vwoodland= -0.064。表明城市化過程中,隨著植被覆蓋度不斷減小,陸表溫度迅速上升。

(2）為從空間上直觀地研究城市化進程對陸表溫度的影響,本文通過疊加處理、圖像相減等相關(guān)操作,得到1995-2003年土地利用變化及陸表溫度變化對比情況(圖4,見封2）。從圖4可以看出,東莞市陸表溫度變化情況與土地利用變化情況相符:增溫幅度最大的區(qū)域主要是農(nóng)用地向城鎮(zhèn)用地轉(zhuǎn)化的部分(圖4a中的藍色部分）及林地轉(zhuǎn)變?yōu)槌擎?zhèn)用地的部分(圖4a中的綠色部分）,增溫區(qū)域在空間上一般以現(xiàn)有的城鎮(zhèn)為中心向外蔓延。

4 結(jié)論

(1）城市化過程造成研究區(qū)下墊面性質(zhì)的劇烈變化是引起陸表平均溫度上升的主要原因,其高溫區(qū)域主要以城鎮(zhèn)為中心,沿道路不斷向兩旁延展。此外,1995年與2003年土地利用類型的溫度序列存在細微的變化,主要是由于環(huán)境污染及盲目開發(fā)造成的。需引起注意的是,1995年各土地利用類型的溫度標準差均大于2003年,說明陸表溫度的大小與分布對環(huán)境變化的影響較敏感。

(2）引入TVX空間方法,分析陸表溫度隨土地利用變化(城市化）的變化情況,進一步論證城市化過程與陸表溫度之間存在著較強的正相關(guān),即城市化是造成陸表溫度上升的重要原因。這種利用TVX空間中聚類點的運動軌跡分析土地利用變化(城市化）對陸表溫度的影響的方法在國內(nèi)鮮見報道。

(3）由于缺乏相關(guān)輔助數(shù)據(jù)的支持,在土地利用分類上只將研究區(qū)分成農(nóng)用地等5類,但對東莞市而言,農(nóng)用地還可細分為園地、耕地等類型,其對陸表溫度的影響各不相同。后期工作中希望可以獲取更多的輔助數(shù)據(jù),進行更精細的土地利用分類,進而更準確地分析不同土地利用類型在城市化進程中對陸表溫度的影響。此外,利用多時相的影像數(shù)據(jù)分析不同季節(jié)城市化過程對陸表溫度的影響也是后續(xù)研究的重點內(nèi)容。

[1] 馮徽徽,夏斌,吳曉青,等.基于SLEU TH模型的東莞市區(qū)城市增長模擬研究[J].地理與地理信息科學(xué),2008,24(6）:76-79.

[2] 崔功豪,馬潤潮.中國自下而上城市化的發(fā)展及其機制[J].地理學(xué)報,1999,54(2）:106-115.

[3] 楊青生,黎夏.珠三角中心鎮(zhèn)城市化對區(qū)域城市空間結(jié)構(gòu)的影響——基于CA的模擬和分析[J].人文地理,2007,94(2）:87-91.

[4] HOWARD L.Climate of London Deduced from Meteorological Observation[M].London:Harvey and Darton,1833.

[5] RAO P K.Remote sensing of urban heat islands from an environmental satellite[J].Bulletin of the American M eteo rological Society,1972,53:647-648.

[6] XIAO H L,WENG Q H.The impact of land use and land cover changeson land surface temperature in a karst area of China[J]. Journal of Environmental Management,2007,85:245-257.

[7] WENG Q H.Estimation of land surface-vegetation abundance relationship for urban heat island studies[J].Remote Sensing of Environment,2004,89(4）:467-483.

[8] 宮阿都,李京,王曉娣,等.北京城市熱島環(huán)境時空變化規(guī)律研究[J].地理與地理信息科學(xué),2005,21(6）:16-18.

[9] 錢樂祥,丁圣彥.珠江三角洲土地覆蓋變化對地表溫度的影響[J].地理學(xué)報,2005,60(5）:761-769.

[10] 劉宇,匡耀求,吳志峰,等.不同土地利用類型對城市地表溫度的影響——以廣東東莞為例[J].地理科學(xué),2006,26(5）:597 -601.

[11] 東莞政府網(wǎng).http://www1.dg.gov.cn/.2009-12-30.

[12] 姚云軍,南鵬,張照錄,等.分裂窗算法在熱紅外遙感陸面溫度反演中的應(yīng)用[J].蘭州理工大學(xué)學(xué)報,2007,33(6）:89-92.

[13] PRICE J C.Land surface temperature measurements from the split window channels of NOAA 7 advanced very high resolution radiometer[J].Journal of Geophysical Research,1984,89 (D5）:231-237.

[14] 覃志豪,ZHANG M H,KARN IEL I A,等.用陸地衛(wèi)星TM 6數(shù)據(jù)演算地表溫度的單窗算法[J].地理學(xué)報,2001,56(4）: 456-466.

[15] Q IN Z H,KARN IEL IA,BERL INER P.A mono-w indow algorithm for retrieving land surface temperature from Landsat TM data and its application to the Israel-Egypt border region [J].International Journal of Remote Sensing,2001,22(18）: 3719-3746.

[16] 覃志豪,L IW J,ZHANGM H,等.單窗算法的大氣參數(shù)估計方法[J].國土資源遙感,2003(2）:37-43.

[17] ARTISD A,CAMAHAN W H.Survey of emissivity variability in thermography of urban areas[J].Remote Sensing of Environment,1982,12:313-329.

[18] SOBRINO J A,JIM ENEZ-MUNOZ J C,PAOL IN IL.Land surface temperature retrieval from Landsat TM 5[J].Remote Sensing of Environment,2004,90:434-440.

[19] SOBRINO J A,RA ISSOUN IN,L I Z L.A comparative study of land surface emissivity retrieval from NOAA data[J].Remote Sensing of Environment,2001,75:256-266.

[20] DOUSSET B,GOURMELON F.Satellitemulti-sensor data analysisof urban surface temperature and land cover[J].ISPRSJournal of Photogrammetry and Remote Sensing,2003,58:43-54.

[21] OWEN T W,CARLSON T N,GILL IES R R.A ssessment of satellite remotely sensed land cover parametersin quantitatively describing the climatic effect of urbanization[J].International Journal of Remote sensing,1998,19:1663-1681.

[22] AM IRIR,WENG Q H,AL IMOHAMMADIA,et al.Spatialtemporal dynamics of land surface temperature in relation to fractional vegetation cover and land use/cover in the Tabriz urban area,Iran[J].Remote Sensing of Environment,2009, 113(12）:2606-2617.

[23] QUA TTROCH ID A,LUVALL JC.Thermal Remote Sensing in Land Surface Processes[M].CRC Press,2004.440.

Influence Analysis of Urban ization on Land Surface Temperature in Dongguan City

FENG Hui-hui,L IU Hui-ping,ZHU Shou-dong,BA IM u
(School of Geography and RS,Beijing N orm al University,the State Key Laboratory of Remote Sensing Science,Beijing 100875,China）

In this paper,the influence of urbanization on land surface temperature(LST）was investigated in Dongguan City. Through the temperature retrieval method p roposed by fo rmer researchers,the LSTs of Dongguan were retrieved from two Landsat TM/ETM imagesw hich were taken in 1995 and 2003.Then,the general situation and spatio-temporal characteristicsof LSTs of Dongguan were discussed,the temperature-vegetation index(TVX）space was also introduced to investigate the influence of urbanization over LST.The result indicated that urbanization was an important reason for the LST increase of Dongguan,w hich was confirmed by the trajecto ry from the dense vegetation-low temperature condition to the sparse vegetation-high temperature condition in TVX space.

urbanization;land surface temperature;spatio-tempo ral characteristics;temperature-vegetation index;Dongguan

TP79;P407

A

1672-0504(2010）03-0087-04

2010-02-27;

2010-04-15

國家863計劃項目(2006AA 120102）;國家自然科學(xué)基金項目(40671127）

馮徽徽(1986-）,男,博士研究生,主要研究方向為遙感城市監(jiān)測與城市時空模擬。＊通訊作者E-mail:hpliu@bnu.edu.cn