随着城市化的发展，城市的规模不断扩大和人口的急剧增加，自然地表覆盖正逐步被人工建筑所替代，改变了当地的景观和气候条件，影响了生态系统功能和过程^{[1, 2]}，产生了诸多的生态环境问题，其中最为突出的现象之一为城市热岛效应日益显著。城市热岛是指城市地区整体或局部温度高于周围地区，温度较高的城市地区被温度较低的郊区所包围或部分包围的现象^[3]。它是由于城市建筑及人类活动导致热量在城区空间范围内集聚^[4]，从而导致城市和郊区的温差明显。城市区域范围的大气温度分布与土地利用强度紧密相关，因而形成城市经典的钟型温度剖面层^[5]，研究表明不透水地面的增加和城市能源的消耗是城市热岛形成的两个主要方面^{[6, 7]}。城市热岛效应已经严重影响了城市生态环境。热岛效应会促使形成复杂的大气污染物环流，加重城区的大气污染^[8]。还会导致地表臭氧浓度升高^[9]，甚至人口死亡率增加^[10]。由于夏季热岛效应形成城市热浪，城市居民更容易受到高温疾病与流行病的困扰^[11]，2003年，法国14800人死于城市热浪^[12]。此外，城市内建筑内部空调使用的能源以及交通的能源需求不断增加，增加了温室气体的排放，加速了全球变暖而直接导致热岛效应越来越严重^[8]。

采用气象台站观测城市地表温度和利用遥感影像反演地表温度的方法是目前研究城市热岛主要的两种研究方法。一些学者利用气象台站的数据分析了北京城区的城市热岛^{[3, 13, 14, 15]}。但受到城市地区观测点数量的限制、空间上的不连续性，通过气象台站的数据拟合城市地区的地表温度具有一定的局限性^{[16, 17]}。由于城市是一个复杂的系统，下垫面性质不断变化，不同类型下垫面的热惯性、热容量、热传导和热辐射的不同导致温度的空间差异十分明显。与气象台站实测数据相比，卫星遥感数据具有空间分辨率高、适用范围大的特点，目前遥感数据已成为研究城市热岛的主要数据来源，遥感信息技术也成为其研究的主要技术手段^[17]。

北京是我国快速城市化发展的典型区域，城市总体规划(2004—2020年)以来，城市规模不断扩大，强烈的人类活动改变了自然地表，大量自然植被地表被水泥、沥青等硬质地表所替代，由于地表不透水面积与地表温度呈正相关^[18]，加上城市能源消费的不断增加，使得城市的热岛格局和热环境特征发生了很大变化。本研究基于遥感反演地表温度的方法，分析上一版城市总体规划以来的北京城市六环内热景观及热岛变化趋势和特征，旨在为城市管理部门指定减缓城市热岛措施提供依据。

研究范围为北京城市六环内城市区域，该范围是城市快速发展的区域，土地利用强度大，社会经济发展的自然生态资产占用迅速，对于该地区的热岛变化特征研究具有典型意义。本文使用了具有30 m空间分辨率的美国陆地资源卫星Landsat TM影像，2003年5月25日、2007年5月28日、2011年6月8日(北京地区)三景数据，轨道号为123/32，三景数据质量较好，云量均低于2%。采用Erads Imagine9.0遥感影像处理软件对TM影像进行大气校正、几何校正、图像增强、合成等过程，将遥感影像与北京地形图(1 ∶ 5万)进行地理坐标配准后，在Arcgis9.2中进行计算并制作专题图进行分析。

城市下垫面的辐射亮温是未经大气校正以像元为单位的平均地表温度，而并非实际意义的地表实际温度，由于该研究主要针对城市热岛的变化趋势进行分析，地表实际温度和地面亮温的趋势基本表现为一致，所以直接采用地面亮温进行分析。

TM遥感影像第6波段是以数据像元值来表示的，数值在0—255之间，数值越大，亮度越大。地面亮温反演方法参考NASA官方的Landsat用户手册提供的算法进行计算^[19]。

地面亮温是地面每一个像元的温度值，为了研究城市区域范围内的热岛效应强度，本文采用热场强度指数这一监测指标。热岛效应的强弱是衡量城市与乡村或郊区之间的温差所形成的城市小气候状况的优劣，根据马蔼乃^[20]遥感信息模型研究中所提出的地理相识准则，类似地引入相对亮温来表示城市热岛效应的强度。即公式：

式中,Ⅱ(T)为各地区亮温的相似准则，T_i为研究范围内i点的亮温，T_a为城市与城郊的平均温度。

为了明确地表示城市的热岛强度，对相似准则进行等级划分，表示不同的热岛等级，采用了孙飒梅和卢昌义的划分方法(表 1)^[21]。

为定量分析城市热景观格局的变化特征，本研究采用景观生态学的指数分析方法来分析城市热岛景观的空间格局变化特征，分别从景观类型水平和景观水平选取相关的指数进行分析，分析工具采用Fragstats 3.3。

表征植被状况的指数较多，但NDVI指数是最为广泛应用的指数之一，在植被遥感中，对于landsat5而言，其计算方法已经成熟，此处不再具体表述。NDVI指数也用于计算植被覆盖度，本文利用像元二分模型估算研究区的植被覆盖度，计算公式如下：

植被覆盖度的计算公式：

式中,NDVI_min和NDVI_max为0.5%—99.5%置信区间的NDVI最小值和最大值。分别表示无植被覆盖区域和纯植被像元的NDVI值。

由于本研究没有获得植被NDVI实测数据，王浩等^[22]认为可以采用一定的置信区间对NDVI取值，可以消除遥感影像噪声所带来的误差。置信区间的设定取决于遥感影像的质量，本文采用0.5%和99.5%的置信区间，计算NDVI最大值与最小值。NDVI指数计算和植被覆盖度的计算分别在Erdas和ArcGIS上完成。

对植被覆盖度的计算结果采用我国土壤侵蚀强度等级划分方法，将植被覆盖度分为5级：f＜10%、10%≤f＞30%、30%≤f＞50%、f≥70%，从小到大依次分别记为Ⅰ、Ⅱ、Ⅲ、Ⅳ、Ⅴ级，生成植被覆盖度等级图进行分析。

从图 1可以看出，2003—2011年城市LST分布格局发生了明显的变化。2003年夏季北京六环内平均LST为292.63 K，最高LST为302.95 K，最低LST为284.87 K。城市中高的LST主要分布在城市三环外南部的丰台和大兴地区，城市北部五环之间六环之间的昌平和顺义地区，其中西南部干涸的永定河河道下游地区高温分布比较明显，而二环内的高LST出现在天安门广场附近。2007年夏季北京六环内的平均LST为303.45 K，最高LST为315.84 K，最低LST为289.03 K。该时段的LST总体呈现星型的分布特征，中心城二环外围地区的热环境变化迅速，丰台、大兴、朝阳、海淀、昌平、顺义地区的高LST分布范围比2003年增长明显，二环内的东西城地区的LST也较2003年有所增加。2011年夏季北京六环内的平均LST为301.75K，最高和最低LST分别为315.84 K和286.72 K。城市总体的LST分布特征相比2007年变化也十分明显，高LST分布呈现围绕中心城集聚的特点，而城市外围的高LST分布面积较之2007年减少也十分明显，其中永定河地区的高温地区基本消失，这可能与2010年北京市实施的永定河绿色生态发展带建设有很大关系。

图 1 2003—2011年城市六环LST变化图 Fig. 1 Land surface temperature changes in study area during 2003—2011

图选项

从2003—2011年夏季城市LST的分布特征分析来看，城市高温地区和低温地区的温度差异在逐渐增大，2003年高低温差18 K，2007年达到26.8 K，而2011年相差达到29.1 K。同时期的地面亮温的高温也从2003年的302.95 K上升至2011年的315.84 K。但六环地区的城市平均LST从2003—2011年经历了先增高后降低的过程，2007年是一个转折点，这与北京为举办奥运会的城市环境建设有很大关系，奥运会举办前期北京市政府相继关停了诸多大耗能工厂，典型的如首钢和焦化厂，奥运后，市政府也积极实施了改善生态环境的建设，如永定河绿色生态带建设、绿化隔离地区建设以及郊野公园建设，城市绿化覆盖率大大提高，从而使得2011年六环内城市整体LST相比2007年有所下降。

由图 2可以看出，城市的绿化覆盖逐年增长，2003年城市五环内Ⅰ等级的植被覆盖地区分布范围较广，同时五环外地区由于保留了成片的农田，Ⅲ、Ⅳ、Ⅴ等级的植被覆盖地区也较多，并且高植被覆盖度分布地区能基本连成一片，绿地之间的连通性较好。2007年五环内的植被覆盖度相比2003年有所提高，北五环外地区的植被覆盖度也大大提高，但南部地区的植被覆盖则相对减少。至2011年，五环内Ⅰ等级的植被覆盖地区基本消失，但五环外围由于城市化的建设，大量农田被占用，外围的高等级的植被覆盖呈现出破碎化趋势，比较典型的是永定河下游沿岸的高等级植被覆盖逐渐向Ⅱ类地区过渡。

图 2 2003—2011年城市六环内植被覆盖度图 Fig. 2 Urban plants cover changes in study area during the year of 2003—2011

图选项

为分析植被覆盖与地面亮温之间的关系，本研究提取了城市不同环线内的植被覆盖度和地面亮温值(LST)做对比分析，见表 2。

北京市六环范围内平均植被覆盖度从2003—2011年呈现逐年上升趋势，六环内的植被覆盖度分别为31.6%、35%和38.9%，城市范围的绿化建设有了明显提高。2003年二环至六环的城市平均LST变化差异不大，热环境格局整体保持在同一水平。2007年五环内平均LST最高，而到2011年，城市热岛已经比较突出，六环内的平均温度相比四环内的平均温度相差1.66 K。从3a的LST对比分析来看，2007年城市的各环线平均温差尽管不十分明显，但六环内的整体温度要比2003年和2011年要高。

为比较城市绿化建设和地面亮温之间的关系，采取随机采样的方法，在六环内提取了2200个样本(不包含水域)的植被覆盖和地面亮温值，分析两者之间的相关关系，对植被覆盖度和地面亮温进行线性回归。分析软件采用spss 11。结果如表 3，其中y为地面亮温，x为植被覆盖度。

表 3通过相关性分析和线性回归分析发现，城市植被覆盖度和城市的热场存在极显著的负相关关系，回归系数揭示了不同年份相同时期城市热场受植被覆盖的影响程度。对于城市土地景观来说，植被覆盖度越高，将导致更高的蒸腾率，同时加速陆地表面和大气层之间潜热与显热的交换，从而降低地表温度^[23]。

图 3是夏季北京城市的热岛等级变化图。2003年，城市的热岛现象尚不十分明显，中等热岛以上等级的分布范围很小。至2007年，二环内中等热岛等级增加明显，同时永定河沿岸地区和东南三环至四环之间的地区热岛等级也增加，由中等热岛向强热岛、较强热岛以及极强热岛过渡，而且，城市西南三环至五环地区的热岛强度也增加明显，五环外大兴亦庄地区热岛相比2003年更加突出，北部顺义、昌平地区的热岛等级均变强。2011年，城市的热岛现象已经十分突出，从图中可以明显看到，五环内的城市热岛范围基本连成一片，相对集中于四环以内，五环至六环之间的西南部房山良乡地区和顺义首都机场地区热岛强度提高明显。在城市五环内弱热岛等级仅分布在部分大型绿地和水面地区。永定河绿色生态带建设实施后，永定河河道下游地区的热岛等级则相比2007年下降十分明显。

图 3 2003—2011年北京城市六环内热岛强度分布图 Fig. 3 The urban heat island intensity in study area during the year of 2003—2011

图选项

表 4是统计的各环线内强热岛(4)等级以上分布面积所占的比例，可以看出2011年城市六环内的强热岛等级以上面积增加很快，其中四环内的热岛等级所占百分比最大，五环内其次。四环内2003年同期还尚未出现热岛，从五环至六环开始出现热岛分布区域。至2007年时，各环线的强热岛以上等级的分布面积都在增加。2011年时出现跳跃式的增长，五环内热岛面积增加明显，但六环内热岛面积相比2007年增加幅度相对较小。

北京市城市总体规划(2004—2020年)提出了城市建设“两轴两带多中心”的方向，城市沿东西长安街、南北中轴线不断扩展，而且强化了城市中心城区的核心功能，同时促进了城市村镇的整合，城市建设的高度和幅度都不断提高。2003—2011年，城市六环内建成区比例由36.3%上升至46.3%，同时城市能源消耗迅速增加，全市小汽车数量突破500万辆，能源消耗的增加和建成区面积的增加导致了城市热岛效应的增强。

2003—2011年城市的热岛景观发生了很大变化，为定量分析城市热岛格局的变化，采用景观指数来表征城市热岛景观格局的变化特征。在景观类型上的景观指数分别采用斑块数(NP)、斑块密度(PD)、斑块类型比例(PLAND)、平均斑块面积(AREA_MN)、景观形状指数(LSI)、面积加权分形指数(FRAC_AM)和聚类指数(CLUMPY)来研究热岛景观类型的固有特性，由于对城市热岛起主导作用的为强热岛级别以上的类型，所以在热岛景观类型上着重选择了强热岛、较强热岛和极强热岛(4、5、6)3个级别；选取面积加权平均形状指数(SHAPE_AM)、香农多样性指数(SHDI)、香农均匀度指数(SHEI)、蔓延度指数(CONTAG)以及景观分离度指数(DIVISION)来研究景观斑块在景观水平上的特征，表征的是各类景观在整体水平上的特征。

由表 5可看出，2003年到2007年，各等级的热岛景观斑块数(NP)、斑块密度(PD)、斑块类型比例(PLAND)和平均斑块面积(AREA_MN)都增加迅速，说明在此时段内，城市强热岛等级以上的斑块在不断的增加，城市的热岛面积和热岛强度相比2003年要大。2011年强热岛等级斑块数(NP)、斑块密度(PD)相比2007年有所下降，但是较强热岛、极强热岛等级的斑块数(NP)、斑块密度(PD)比2007年要大，尽管较强热岛和极强热岛平均斑块面积(AREA_MN)有所减少，但各级斑块类型比例在增加，导致2011年比2007年城市热岛强度和热岛范围更大。2003年至2011年，各热岛等级的景观形状指数一直呈现增长趋势，说明城市的热岛等级斑块的形状更加复杂，呈现多样化。面积加权平均分形指数(FRAC_AM)也有所增加，该指标能反映人类活动对于景观格局的影响，值越大，人类活动影响越大，说明由于城市化过程，人类活动对于热岛景观格局产生了很大的影响。聚类指数是用来指示斑块类型的分布状态的，当聚类指数(CLUMPY)接近于1时，斑块类型最大程度的聚集。从表 4中可以看出，2007年各等级的聚类程度达到最大值，说明该时期的热岛斑块集聚程度高，2011年该指数略有下降。

面积加权平均形状指数(SHAPE_AM)反映景观的规则性，值越大，景观形状越不规则。从城市热岛在景观水平上的比较(表 6)可以看出，2003年，热景观斑块的形状最为丰富，到2011年热景观斑块的形状指数又相比2007年有所增加，说明热岛景观形状又呈现多样化。2007年城市热岛斑块的多样性指数(SHDI)和均匀度指数(SHDI)最高，反映了热岛景观斑块类型多，而且在城市范围内分布相比2003年和2011年要更为均匀。蔓延度指数(CONTAG)是描述景观内各斑块类型的蔓延趋势，当某种优势斑块类型形成了好的连接性，则该值较高，2003—2011年，六环内热岛景观的蔓延度下降后又有升高，总体上呈现下降趋势，低等级的热岛景观的比重下降，优势度降低。2007年景观分离度(DIVISION)达到最大值，说明热岛斑块的破碎化程度高。

本文以2003年、2007年、2011年夏季相同时段的Landsat TM遥感影像为基础，运用RS和GIS工具以及景观分析软件，分析了北京城市六环内的城市热景观变化特征，得出了以下结论：

(1)2003—2011年，城市化的迅速发展，导致了北京城市六环内的热环境分布格局发生了明显变化，城市高温的范围出现了由分散到集中的转变，高温呈现围绕城市中心城集聚的特点，而且城市与城郊的温差逐步增大。

(2)城市的植被覆盖度与地表温度之间的相关性分析表明，二者之间存在显著的负相关关系。城市的绿化建设有了明显提高，植被覆盖度逐年增长，绿化建设和城市化的发展影响着地表温度的分布。相比2003年，2007年的植被覆盖度也有明显提高，但是由于城市化发展，城市中心外围大型绿地农田被城市硬化地表占用，导致各环线内的平均温度要比2003年高。

(3)城市的热岛等级和分布范围变化明显。至2011年，高等级热岛斑块(4、5、6)等级的分布范围和强度都明显提高。而且城市五环内的高热岛等级范围基本连成一片，城市中心的热岛效应更加明显。

(4)随着城市建设规模的不断扩大，人类活动强度增大，城市中高等级热岛斑块(4、5、6)的数量、密度和比例都增加迅速，而且形状更加复杂多样化，同等级类型的斑块集聚程度高。2007年城市的热景观斑块类型最多，而且分布均匀度高，同时热景观斑块的破碎度高。热景观斑块的蔓延度呈现下降趋势，低等级热岛(1、2、3)斑块的优势度逐渐降低。

土地利用/覆盖影响着城市热岛的分布^{[24, 25]}。遥感影像反演的地表温度是地物所释放的辐射能量，包括建筑屋顶，地面铺装，植被，裸露地表和水体，城市中高的地表温度往往出现在高比例的不透水地表地区，两者之间存在显著的线性相关关系，这些不透水地表主要包括道路和屋顶，因此，城市土地覆盖类型会成为影响地表温度的重要因素^[26]。除此以外，城市土地利用类型的在景观尺度上的空间分布也是影响城市热岛的重要因素，因为不同的城市用地类型对城市的热岛贡献存在差异，土地利用的变化会导致城市中的物质和能量的流动，这会导致城市热岛发生变化。Li^[27]等的研究表明，住宅用地对城市中热岛形成的贡献最大，其次是工业用地，而且，中低密度或低植被覆盖的住宅用地比高密度或髙植被覆盖的住宅用地的温度要高，所以，建筑的空间形态以建筑密度和容积率对城市的热环境形成发挥着主导的作用^[28]。

自北京市城市总体规划(2004—2020年)实施以来，城市化快速发展，城市城镇人口规模由2003年的1151万人增长到2011年的1740万人，城市建设用地规模不断扩大，与此同时，城市化过程中能源消费急剧增长，2003年全市能源消费总量4648万吨标煤，2011年消费6995万吨标煤^[29]。尽管研究区内城市的绿化覆盖率相比2003年有了大幅的提升，特别是四环内绿化覆盖率提高最为迅速，提高了近16%，但城市热岛分布面积比例却也相应的增加最为迅速。由于植被覆盖率与热岛存在显著的负相关，而城市热岛与地表覆盖类型和能源消耗直接相关，所以城市热岛比例的增加在很大程度上是由城市能源消耗导致的温室气体排放所造成的。

在景观斑块类型水平的指数分析表明，各热岛等级的景观形状指数呈现增长趋势，热岛等级斑块形状复杂化，面积加权平均分形指数和各斑块类型的比例都有增加。景观水平的分析表明，热岛斑块的多样性和均匀度指数在2007年达到最大，2011年略有下降，蔓延度指数分析表明城市低等级的热岛景观优势度降低。城市用地类型的空间布局和地表覆盖类型将会影响到这些指标，工业企业以及大面积低密度住宅区或城中村的空间分布会影响到能源排放的空间分布，进而影响到热岛的景观指数。绿地类型(大小、植被类型)及空间分布同样会影响到热岛的景观指数。

为减缓城市的热岛效应，城市规划和管理部门需要进行多方面的考虑：(1)优化城市用地的空间布局，特别是居住用地和工业用地的空间布局；对于中低密度的城中村或老住宅区进行节能改造。加快工业用地和村庄的整合力度，集约用地。(2)疏解城市人口，降低城区人口密度，城市人口增加导致的能源消耗增加对于城市热岛的贡献有间接的关系；(3)积极推广清洁能源的利用；发展低碳交通，降低交通热源排放；(4)继续扩大城市公园绿地的比重，增加城市植被覆盖度，对绿化的植被类型和配置进行合理设计，优化植被群落的生态效应。对于城市中的未利用地、裸地、干涸河道、滩地等进行生态恢复。(5)有条件地实施雨水收集措施，增加城市水面，水体具有“冷岛”效应，对于缓解城市热岛同样具有积极作用。