王朝陽，宋辛辛，李海懿

(1.重慶師范大學 地理與旅游學院，重慶 401331；2.鄭州師范學院 地理與旅游學院，河南 鄭州 450044)

1 引言

近年來，我國新型城鎮(zhèn)化建設進程不斷加速，大量建筑和道路集中構(gòu)建，自然地表被破壞，廢氣廢熱高強度排放等引起下墊面性質(zhì)及結(jié)構(gòu)發(fā)生改變，城市出現(xiàn)“高溫孤島”現(xiàn)象[1～3]。相關(guān)學者[4,5]已經(jīng)證實了城市熱島效應是城市化對氣候影響最顯著的特征。

國內(nèi)大多數(shù)學者采用氣象學和遙感學方法對城市熱島效應進行研究。黃群芳[6]通過研究表明城市空間形態(tài)可以顯著改變城市熱島強度。張曉東等[7]利用景觀格局和緩沖區(qū)分析方法，發(fā)現(xiàn)適當增加綠地斑塊形狀復雜度，可以很好地改善城市熱環(huán)境。李玉杰等[8]基于“源-匯”理論，發(fā)現(xiàn)不同匯源景觀的平均熱場變異指數(shù)不同。劉麗香等[9]利用輻射傳輸法反演地表溫度，發(fā)現(xiàn)城市不同功能區(qū)對地表溫度作用程度不同。錢敏蕾等[10]揭示了城市熱環(huán)境的變化與城市用地擴張或縮減的規(guī)律。潘瑩等[11]采用了定量和定性分析的手段，揭示了地表溫度與下墊面的關(guān)系，并表明城市熱島效應與環(huán)境要素關(guān)系密切。Streutker[12]發(fā)現(xiàn)城市熱島效應的增強與人口的增長息息相關(guān)。因此，有必要對城市熱島的形成機制及時空分布進行研究，以緩解地區(qū)的城市熱島效應。

鄭州城市化進程不斷加快，城市規(guī)模迅速擴大，城市結(jié)構(gòu)也隨之發(fā)生變化，這必然會引起熱島效應變化。學界對于鄭州市城市熱島效應的研究開始較早，程炳巖等[13]分析了鄭州的熱島特征及其隨著季節(jié)和高度變化的規(guī)律，并討論了關(guān)于鄭州城市熱島效應的成因機制。李芳芳等[14]基于ETM數(shù)據(jù)分析了鄭州市城市熱島效應及其影響因素。張凡等[15]利用2007～2016年MODIS地表溫度產(chǎn)品，分析了鄭州市城市熱島的時空演變特征。耿麗文[16]在風景園林視角下研究鄭州中心城區(qū)的城市熱島效應時空演變。目前利用長時間序列的遙感資料對鄭州城市熱島效應晝間和夜間年際變化的研究不多，本文以鄭州為實證對象，采用2011～2020年MODIS溫度產(chǎn)品，利用遙感技術(shù)研究鄭州城市熱島效應特點及其演變，分析熱島效應的晝間和夜間的年際變化規(guī)律，以期為鄭州市未來實現(xiàn)綠色發(fā)展提供參考。

2 研究區(qū)概況

鄭州市是河南省會城市，位于河南省中部偏北，屬于我國中部的中心地區(qū)，地處東經(jīng)112°42′～114°14′，北緯34°16′～34°58′。鄭州市氣溫年較差大，屬于溫帶大陸性氣候。鄭州市處于黃河下游，約有124條河流，流域面積11.4 km2。地勢西南高東北低。2020年鄭州市GDP產(chǎn)值已達12003億元。近10年來，鄭州中心城區(qū)建成區(qū)面積迅速擴張，由354.7 km2擴張至709.69 km2，城市的擴張必然會引起熱島效應的變化，因此有必要對鄭州市開展城市熱島效應的研究。

3 數(shù)據(jù)來源和研究方法

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

使用的地表溫度產(chǎn)品(LST，Land Surface Temperature)是通過美國國家航空航天局(NASA)的數(shù)據(jù)下載網(wǎng)站選取8d合成空間分辨率為1km的MOD11A2地表溫度產(chǎn)品，以2011年、2015年和2020年為代表年份，選擇覆蓋研究區(qū)云量較少、圖像質(zhì)量較高的數(shù)據(jù)。

土地利用數(shù)據(jù)來源于中國科學院資源環(huán)境科學數(shù)據(jù)中心提供的2015年Landsat TM影像(空間分辨率30m)解譯數(shù)據(jù)，社會經(jīng)濟數(shù)據(jù)來源于《鄭州統(tǒng)計年鑒(2011-2020)》、鄭州市統(tǒng)計局。

3.2 數(shù)據(jù)處理

本研究利用ENVI和ArcGIS軟件對2011～2020年共計138幅影像數(shù)據(jù)進行分析，先對數(shù)據(jù)進行投影矯正等預處理。利用ENVI中的拓展工具MCTK(The Modis Coversion Toolkit)實現(xiàn)數(shù)據(jù)坐標的轉(zhuǎn)換和重投影，提取LST_Day_1km和LST_Night_1km對城市熱島效應進行分析，利用公式(1)將MODIS溫度產(chǎn)品的單位轉(zhuǎn)為℃。

LST=DN×0.02-273.15

(1)

式(1)中:DN為原始像素值，LST為地表溫度。采用規(guī)則裁剪的方法，利用鄭州市的行政界限在ArcGIS中對處理后的圖形進行掩膜裁剪，去除殘缺數(shù)據(jù)之后求均值，得到不同年份各個像元溫度的均值，利用密度分割法對所得影像數(shù)據(jù)劃分5個不同的地表溫度等級，并利用RS和GIS計算出不同溫度區(qū)所占面積的比例.

3.3 熱島指標

本研究利用UHIER指數(shù)[17]對熱島等級進行劃分，UHIER指數(shù)的計算方法如下：

UHIERn=(Tn-Tmin)/(Tmax-Tmin)

(2)

式(1)中，UHIERn代表第n個像元的城市熱島強度，Tn表示第n個像元的地表溫度值，Tmax則代表研究區(qū)地表溫度最大值，Tmin代表研究區(qū)地表溫度的最小值。UHIER指數(shù)可以反映熱島強度，利用密度分割法[17]將計算結(jié)果劃分為5個等級。

4 結(jié)果與分析

4.1 鄭州熱島平均場空間分布

采用2011～2020年每個像素的平均值表示熱島平均場(圖1)，以高溫區(qū)和較高溫區(qū)面積之和表示熱島范圍，熱島中心為高溫區(qū)最集中區(qū)域，熱島效應的明顯程度反應熱島強度，以高溫區(qū)面積進行定量描述。對2011～2020年熱島平均場空間分布進行研究，可以反映出鄭州市熱島范圍的大體分布。

圖1 2011～2020年研究區(qū)熱島平均場空間分布

從空間上來看，熱島范圍由點狀中心向外擴散，熱島中心的分布范圍相對來說更為集中，大部分位于鄭州市市轄區(qū)，而這些區(qū)域正是鄭州市經(jīng)濟發(fā)展以及城市化水平較高的地區(qū)。高溫區(qū)所占面積的比例為4.84%，較高溫區(qū)在市轄區(qū)的周邊地區(qū)、鞏義市、登封市、新鄭市、新密市及中牟縣均有分布。在范圍上，較高溫區(qū)分布范圍比高溫區(qū)大，所占面積為26.93%。鄭州市多數(shù)區(qū)域還是以中溫區(qū)為主，所占面積為44.95%。低溫區(qū)和較低溫區(qū)分布于鄭州的北部、東北部、登封北部及南部，所占面積為23.25%。

4.2 熱島效應時空變化分析

利用RS和GIS技術(shù)對研究區(qū)不同溫度等級所占面積進行計算，獲得鄭州2011～2020年晝夜各溫度區(qū)所占比例(表1)。近10年間，熱島范圍所占面積基本穩(wěn)定在25%～26%，沒有預期的擴張現(xiàn)象，也沒有明顯的萎縮趨勢，總體趨勢較為平穩(wěn)。

表1 2011～2020年晝夜各溫度區(qū)所占的比例 %

從空間分布來看，2011年高溫區(qū)主要集中分布在鄭州市市轄區(qū)，其他地區(qū)分布較少，所占面積為4.89%。較高溫區(qū)的分布比較分散，研究區(qū)內(nèi)均有分布(圖2a)，所占面積為21.56%；2015年較高溫區(qū)分布相對集中，主要分布在市轄區(qū)及其周邊地區(qū)(圖2b)，所占面積為19.82%；2020年較高溫區(qū)分布區(qū)域與2015年相似，分布范圍有所擴張(圖2c)，所占面積為4.55%.低溫區(qū)和較低溫區(qū)分布范圍保持相對穩(wěn)定，主要分布于鄭州西部、北部以及東北部地區(qū)，其中東北部地區(qū)的較低溫區(qū)在2011年后有縮減趨勢。從熱島中心的空間分布來看，自2011～2020年，市轄區(qū)的熱島中心位置從未發(fā)生變動(圖2)。

圖2 研究區(qū)熱島平均場空間分布

從時間上來看，2011～2020年高溫區(qū)分布范圍先擴大后縮減，較高溫區(qū)分布范圍先縮減后擴大，新密市的熱島范圍在2011～2015年有所縮減，但2015～2020年又被高溫區(qū)覆蓋。中溫區(qū)分布范圍持續(xù)擴大，這與王燕強等[19]對鄭州市城市熱島的研究結(jié)論一致。低溫區(qū)和較低溫區(qū)分布范圍持續(xù)縮減，中牟縣大量較低溫區(qū)轉(zhuǎn)為中溫區(qū)。計算每個像元的溫度求平均值，得到鄭州市2011年平均地表溫度為17.03 ℃，2015年平均地表溫度為18.24 ℃，2020年平均地表溫度為17.38 ℃，2011～2020年平均地表溫度基本維持在17～19 ℃之間(圖3)。

圖3 研究區(qū)2011～2020年晝夜平均地表溫度變化

自2011～2020年市轄區(qū)一直處于熱島范圍內(nèi)。相關(guān)研究[7,19,20]表明，城市熱島效應與城市的植被覆蓋率呈現(xiàn)出負相關(guān)關(guān)系。鄭州城市擴建造成植被覆蓋度減小[19]。汪洋等[20]學者認為在都市區(qū)內(nèi)部的高密度建成區(qū)，植被覆蓋度具有主導性約束力，因此緩解鄭州市的城市熱效應，可以通過發(fā)展城市綠化來實現(xiàn)。

在城市化推進過程中，鄭州市城區(qū)面積不斷擴大。劉宇峰等[21]研究表明：不同因素對于熱島強度變化的影響程度也存在著差異，其中城市常住人口數(shù)、城市常住人口密度以及城市建成區(qū)面積，對城市熱島效應的影響最為顯著。在市轄區(qū)人口較多，極易形成一個強熱島，而周邊地區(qū)較難形成強熱島區(qū)。鄭州城區(qū)建成時間較早。相關(guān)研究[22,23]表明，鄭州城區(qū)面積在2001～2010年間建設用地迅速向四周擴張，城市下墊面性質(zhì)改變，致使城市熱島效應顯著增強，因而可以通過優(yōu)化鄭州市國土空間格局來緩解城市熱島效應。

4.2.1 晝間熱島時空變化分析

晝間熱島范圍(表1)和中心(圖4)波動較大，2011年熱島范圍所占面積達到42.13%， 2015年后，熱島范圍基本上穩(wěn)定在34%～35%，晝間熱島范圍有微弱縮減趨勢，高溫區(qū)和較高溫區(qū)的變化趨勢明顯(表1)。

從空間分布來看，2011年較高溫區(qū)在鄭州市轄區(qū)、新鄭市、中牟縣、新密市和登封市分布廣泛(圖4a)，所占面積達到了39.63%。高溫區(qū)只有零星分布，所占面積為2.50%。低溫區(qū)和較低溫區(qū)主要分布在鄭州北部、東北部以及登封北部，所占面積為24.37%；2015年高溫區(qū)集中分布在鄭州市市轄區(qū)、新鄭市以及中牟縣(圖4b)，所占面積為13.41%。低溫區(qū)與較低溫區(qū)分布與2011年基本保持相似；2020年高溫區(qū)分布區(qū)域極小(圖4c)，所占面積僅為0.90%。較高溫區(qū)主要分布在鄭州市市轄區(qū)、登封市以及中牟縣，所占面積為33.26%。

圖4 研究區(qū)晝間熱島平均場空間分布

從時間上來看，2011～2020年高溫區(qū)分布范圍先擴大后縮小，較高溫區(qū)分布范圍先縮小后擴大，中溫區(qū)分布范圍一直保持擴大的趨勢，較低溫區(qū)分布范圍先擴大后縮小，低溫區(qū)分布范圍先減小后擴大。 2011年鄭州市的高溫區(qū)只有零星分布，較高溫區(qū)分布范圍廣,2015年較高溫區(qū)面積縮減，高溫區(qū)面積擴散，使得市中心熱島強度增加，熱島中心向新鄭市擴散，2020年鄭州市市轄區(qū)以及東南地區(qū)的高溫區(qū)大幅衰減，轉(zhuǎn)而成為較高溫區(qū)。從2011～2020年，鄭州北部、西部以及西南邊緣都為較低溫區(qū)和低溫區(qū)，結(jié)合鄭州土地利用類型分布，鄭州北鄰黃河，西部以及西南邊緣都分布有大量林地。由圖3得到，2011～2020年晝間平均地表溫度先上升后下降，與熱島強度的變化趨勢保持一致。

綜上所述，晝間熱島范圍先縮小后擴大。結(jié)合數(shù)據(jù)[24,25]，鄭州市2011年生產(chǎn)總值為4952.1億元，到2020年增加至12003.0億元；2010年鄭州市中心城區(qū)建成區(qū)面積為316 km2，到2020年增加至709.69 km；2010年鄭州市區(qū)常住人口為862.7萬人，到2020年增加至1262萬人。經(jīng)濟高速發(fā)展以及城鎮(zhèn)化建設是鄭州市中溫區(qū)面積增加的主要原因。熱島效應的變化與城市化關(guān)系密切，伴隨城市化向新鄭市等地區(qū)推進，熱島中心有向新鄭市逐漸轉(zhuǎn)移的趨勢。

4.2.2 夜間熱島時空變化分析

夜間的城市熱島范圍相對于晝間縮小程度更加明顯(表1)。熱島范圍先縮小后擴大，熱島范圍所占面積維持在25%～28%之內(nèi)，各溫度區(qū)所占面積變化較小，市轄區(qū)的熱島中心位置從未發(fā)生變動。

從空間分布來看，2011年高溫區(qū)主要分布在鄭州市市轄區(qū)，其他地區(qū)只是零星分布，高溫區(qū)所占面積為5.82%(表1)。低溫中心多分布于登封北部，低溫區(qū)和較低溫區(qū)主要分布于鄭州的北部、西南以及東南地區(qū)，分布范圍較廣(圖5a)，所占面積達到了33.58%；2015年鄭州市市轄區(qū)仍以高溫區(qū)為主導，鞏義市部分地區(qū)的較高溫區(qū)轉(zhuǎn)為高溫區(qū)(圖5b)，高溫區(qū)所占面積為7.08%，低溫區(qū)和較低溫區(qū)的分布范圍與2011年基本相似，所占面積為35.35%；2020年高溫區(qū)分布仍以市轄區(qū)為主(圖5c)，所占面積為6.47%。低溫區(qū)和較低溫區(qū)仍主要分布于鄭州的北部、西南和東南地區(qū)，所占面積為33.92%；2011～2020年新鄭市、登封市以及中牟縣的南部都保持著較大面積的低溫區(qū)，其他地區(qū)的低溫區(qū)只是零星分布，未形成較大規(guī)模(圖5)。

圖5 研究區(qū)夜間熱島平均場空間分布

從時間上來看，2011～2020年高溫區(qū)分布范圍先擴大后縮小，較高溫區(qū)則出現(xiàn)與高溫區(qū)相反趨勢，低溫區(qū)分布范圍先縮小后擴大，較低溫區(qū)分布范圍先擴大后縮小。夜間的高溫中心不同于晝間，自2011年以來，夜間的高溫中心一直存在于市轄區(qū)。中溫區(qū)基本保持穩(wěn)定趨勢。在鞏義市，2011年有大量較高溫區(qū)分布，2015年部分較高溫區(qū)轉(zhuǎn)變?yōu)楦邷貐^(qū)，2020年強熱島中心開始縮減；鄭州西南、東南地區(qū)始終大面積處于較低溫區(qū)；新鄭市自2011～2020年一直存在著一個小面積的高溫中心。鞏義市的高溫區(qū)在2015～2020年間有縮減的趨勢。經(jīng)過計算得到鄭州市2011、2015和2020年夜間平均地表溫度分別為11.58 ℃、12.89 ℃和13.53 ℃，夜間平均地表溫度有逐漸上升的趨勢(圖3)。

綜上所述，夜間熱島變化范圍比晝間小。結(jié)合數(shù)據(jù)及圖像可得:夜間熱島范圍先縮小再擴大，趨勢與晝間保持一致。夜間城市熱島效應比晝間弱，說明太陽輻射減弱，人口流動減少，交通壓力減小等在一定程度上可緩解或減弱城市熱島效應。

4.3 熱島效應空間分布與用地關(guān)系

為了探究熱島效應的空間分布與城市用地的關(guān)系，選取研究期內(nèi)中間年份2015年為代表，分析該時期的土地利用類型(圖6)與2015年熱島平均場的空間分布(圖2b)關(guān)系。

從空間上看，2015年熱島范圍集中分布與市轄區(qū)及周邊地區(qū)，滎陽市、鞏義市、新密市以及登封市均只分布小面積高溫中心，高溫中心呈現(xiàn)零星分布的特點，結(jié)合圖6這些高溫中心都為城鄉(xiāng)居民用地。鄭州市多數(shù)區(qū)域以中溫區(qū)為主，結(jié)合其土地利用類型，中溫區(qū)以水田旱地為主，可知耕地、草地、林地以及水域在一定程度上可緩解或減弱城市熱島效應。

低溫區(qū)分布在登封的北部、南部以及鄭州的北部，此區(qū)域林地是其主要土地利用類型，鄭州北鄰黃河，所以此區(qū)域為低溫區(qū)，土地利用類型為水域(圖6)。通過對比平均場的空間分布和土地利用類型，鄭州市城市熱島效應的空間分布范圍和形式與鄭州市的下墊面性質(zhì)基本保持一致，下墊面性質(zhì)可以影響城市溫度。

圖6 2015年研究區(qū)土地利用類型

5 結(jié)論

本文利用遙感和GIS技術(shù)，根據(jù)2011年、2015年和2020年的MODIS數(shù)據(jù)、社會經(jīng)濟數(shù)據(jù)以及土地利用數(shù)據(jù)等，經(jīng)過分析和綜合處理等一系列過程，對鄭州市長時序熱島效應變化進行了一系列的時間和空間分析，總結(jié)出鄭州市熱島效應的變化特點，得出以下結(jié)論。

(1)鄭州市的城市熱力環(huán)境主要以市轄區(qū)為主，此地區(qū)為經(jīng)濟發(fā)達、人口密度大、建筑密集地區(qū)；2011～2020年高溫區(qū)有相對縮減的趨勢，2020年高溫區(qū)所占面積較2011年減少了0.34%，(較)低溫區(qū)逐漸減少轉(zhuǎn)為中溫區(qū)。

(2)鄭州市以中溫區(qū)為主，熱島區(qū)的擴展趨勢與鄭州城市化息息相關(guān)，隨著城市化向周邊地區(qū)擴張，熱島效應集中區(qū)也開始向建筑地密集區(qū)及新鄭市轉(zhuǎn)移，城市下墊面性質(zhì)可影響城市溫度，城鄉(xiāng)居民用地多為高溫區(qū)，水域和林地多為(較)低溫區(qū)。

(3)從晝夜熱島平均場分布情況來看，夜間熱島的中心表面呈現(xiàn)出較大的團體狀分布，熱島中心位置未發(fā)生明顯移動，2011～2020年夜間平均熱島范圍所占面積為26.68%，晝間平均熱島范圍所占面積為36.80%。