基金项目:十三五国家重点研发计划子课题(2018YFC0704505); 国家自然科学基金青年基金(51608419); 陕西省自然科学基金面上基金(2018JM5073); 陕西省教育厅专项基金(18JK0453); 西安建筑科技大学自然科学专项基金(ZR19001)
第一作者:冯晓刚(1979-),男,博士,副教授,主要从事城市环境遥感及智慧城市设计.E-mail:fendao_ren@163.com.
(School of Architecture,Xi'an Univ. of Arch. & Tech.,Xi'an 710055,China)
urban heat island; thermal landscape; pattern evolution; Xianyang and Xi'an city
DOI: 10.15986/j.1006-7930.2021.03.013
土地利用类型的急剧变化是快速城镇化最显著的特征[1].原有耕地、植被、水体等自然地表类型快速转变为建设用地,地表热通量增大[2].同时,伴随区域经济一体化快速发展,人为热排放不断增加,两者共同作用导致近地面热辐射通量交换过程显著变化,城市逆温加剧,导致城市热环境持续恶化.因此,城市热环境研究已经成为当下构建低碳、适宜与可持续发展新型城镇化过程中人居热生态安全研究的热点 [3].
城市热环境研究起源于19世纪初期的英国[4].此后,国内外众多学者从多视角、多方法体系对不同样区分别展开了实验研究[5-6],其中城市热岛效应的研究最为集中.梳理已有研究成果发现:分形理论、景观生态学理论等与气候学的交叉与融合为城市热环境在机理研究、时空格局演变分析、成因及对策机制等研究中提供了理论支撑,有力的促进了多学科的交叉与融合.陈云浩等采用景观指数对上海市热力景观演变过程进行了剖析,结果表明景观格局指数可以有效表征城市热环境的时空演化特征[7].冯晓刚等利用景观格局指数对西安市长时序的热环境展开了剖析,结果表明景观格局指数能够深入刻画景观类型水平下的热环境时空演变过程与趋势[8].沈中健等探索了绿地景观格局对城市热环境时空分布的影响,结果表明两者间存在显著的空间自相关性[9].因此,景观生态学理论的融入为城市热环境的研究提供了时空表征的新思路.
在“一带一路”建设的大背景下,西咸一体化进程不断加快,导致区域热环境恶化趋势明显.本文以景观生态学理论为指导,基于分异性原理,以Landsat TM/ETM+与Landsat8为主要数据源,辅以气象与实测数据,综合运用单窗算法[10]深入剖析了西咸一体化背景下咸阳市1988~2019年31年间城市热环境的演变过程与趋势特征,以期为有效改善西咸一体化进程中的城市热环境问题提供科学参考.
咸阳市是中国大地原点所在地,位于陕西省中部(东经107°38'~109°10',北纬34°11'~35°32'),东邻陕西省省会西安市,渭水穿南,西接杨凌高新农业示范区,北与铜川、渭南相接壤.辖2区2市9县,全市总面积10 196 km2,截止2019年底,全市总人口500余万,常住人口435万人; 地势由东南向西北呈阶梯形态势分布.属暖温带大陆性季风气候,南部地区四季分明、雨量充沛; 北部地区,气候稍寒.全境年均气温9~13.7 ℃,最高41.2 ℃,最低-26.7 ℃,全年降水量500~600 mm[11].
以热红外遥感数据Landsat TM/ETM+/TIRS为基础,影像成像时间分别为1988-6-4、1995-5-4、2016-5-5和2019-6-5.轨道号均为p127r36,多光谱波段、热红外波段与全色波段空间分辨率分别为30 m、60/120 m和15 m.辅助数据包括咸阳市2016年1:1万地形图,2016年SPOT卫星数据,30 m分辨率的DEM及1988~2019年咸阳市统计年鉴数据.
本文以2016年ETM+数据为基准与咸阳市2016年1:1万地形图进行空间配准,配准误差控制在0.5个像元内; 在进行辐射定标的基础上,基于MODTRAN5辐射传输模型对四个时相的多光谱数据分别进行大气辐射校正; 然后分别确定多时相遥感数据对应的地表比辐射率、大气平均作用温度和大气透过率等核心反演参数,最后综合运用基于Python自主编写的热环境反演程序计算得到研究区实际地表温度,并根据需要将反演结果依次划分为强绿岛区、绿岛区、常温区、热岛区和强热岛区五种类型[14],在此基础上,结合景观格局指数对咸阳市1988~2019年31年间城市热环境时空割据演化特征展开剖析.
基于Landsat热红外波段数据计算地表温度,首先需利用亮温模型计算得到亮温温度,然后采用单窗算法计算得到实际地表温度[12].
(1)亮温温度计算
Lλ=DN×(Lmax-Lmin)/255+Lmin (1)
T6=K2/Ln(K1/Lλ+1) (2)
式中:Lλ为热辐射强度值,DN为热波段的像元亮度值; Lmax、Lmin分别为传感器可探测到的最大、最小辐射亮度; 对Landsat 5:Lmax= 1.560 0(W·m-2·sr-1·μm-1)、Lmin = 0.123 8(W·m-2·sr-1·μm-1),K1 = 60.776(W·m-2·sr-1·μm-1)、K2 = 1 260.56(K); Landsat 7:Lmax =17.04(W·m-2·sr-1·μm-1)、Lmin=0.0(W·m-2·sr-1·μm-1),K1 = 666.09(W·m-2·sr-1·μm-1)、K2 =1 282.71(K)[10].通过ENVI5.2构建地表亮温反演模型,计算得到研究区不同年份地表亮温温度.
(2)参数计算
由于Landsat影像大都是混合像元,基于混合像元及其分解理论,确定地表比辐射率ε= fεv+(1-f )εi +dε,其中:εv和εi分别为植被比辐射率和裸露地比辐射率,通常取值0.985和0.960[13]; f为植被盖度,f =( NDVI-NDVImin )/(NDVImax-NDVImin ); dε为地表几何分布和内散射效应,dε= 0.06 × f ×(1-f ).由于研究区位于北半球中纬度地区,因此在获取当日近地面气温T0基础上,利用Ta =16.011 0 + 0.926 21 T0计算得到大气平均作用温度[13].由于大气透过率较难获取,本文借鉴覃志豪[10]等提出的大气透过率估算方程进行大气透过率参数的估算(表1).
(3)地表温度计算
以地表比辐射绿、大气平均作用温度和大气透过率为实际地表温度计算的核心参数,基于上述计算结果,利用公式(3)计算得到实际地表温度.
TS=1/C{a(1-C-D)+[b(1-C-D)+C=D]×Tb-D×Ta} (3)
式中:TS为实际地表反演温度; Ta为大气平均作用温度; Tb为亮度温度; C和D为中间过程变量,C = ετ,D = (1-τ)×[1+(1–ε)τ]; a和b为拟合系数,通常取-67.355 351和0.458 606; ε和τ分别为地表比辐射率和大气透过率[8].基于(3)计算得到咸阳市1988、1995、2006和2019年四个时期的地表温度数据.在对计算结果进行归一化处理[14]的基础上得到咸阳市归一化温度数据(图2).
表1 大气平均作用温度与近地面气温关系式(覃志豪,2003)[13]
Tab.1 The relationship between the average temperature of atmosphere and the ground temperature
基于上述得到的咸阳市1988、1995、2006和2019年温度数据.运用Vector Utilities模块,将影像由Raster转换为Vector格式,并在ArcGIS10.5平台下,经归一化处理后分别得到不同时相的热环境数据,结果见表2.
由表2可知:1988年常温区面积为431.31 km2,占咸阳市总面积的79.75%; 热岛区和强热岛区面积为57.86 km2,占总面积的10.66%; 绿岛区和强绿岛区面积为51.97 km2,占总面积的9.87%,由此可见1988年咸阳市常温区占主导地位,其次为热岛区域,最后为绿岛区
表2 不同年份咸阳市归一化温度面积统计表
Tab.2 The Statistical table of normalized temperature area in Xianyang in 1988、1995、2006 and 2019
域.1995年常温区面积为361.24 km2,占咸阳市总面积的66.56%; 热岛区和强热岛区面积为90.58 km2,占总面积的16.69%; 绿岛区和强绿岛区面积为90.92 km2,占总面积的16.75%,表明1995年咸阳市常温区继续占主导地位,其次为热岛区域,最后为绿岛区域.2006年常温区面积为451.55 km2,占咸阳市总面积的86.05%; 热岛区和强热岛区面积为 57.56 km2,占总面积的10.96%; 绿岛区和强绿岛区面积为33.65 km2,占总面积的6.42%,表明2006年咸阳市常温区继续占主导地位,且面积显著较1995年增加,其次为热岛区,最后为绿岛区.2019年常温区面积为326.66 km2,占咸阳市总面积的60.19%; 热岛区和强热岛区面积为175.18 km2,占总面积的32.27%; 绿岛区和强绿岛区面积为39.03 km2,占总面积的7.53%,由此表明尽管2019年咸阳市常温区继续占主导地位,但热区域显著增加,常温区减少面积主要流向了热区域; 同时,绿区域较2006年表现为微小幅度增加的趋势.由此可见,1988~2019年间咸阳市始终常温区占主导地位,且2006年占比最大为86.05%,2019年占比最小为60.19%; 热区域面积次之,且2019年占比最大为32.27%.综合分析:1988~2019年31年间热区域变化幅度不大,绿区域占总面积最小,但整体呈现下降的趋势.
为了进一步剖析1988~2019年咸阳市热力景观的变化情况,以1988、1995、2006和2019年为时间节点分别计算咸阳市不同热力斑块的变化情况(图3).结合图3和表2可知:1988~1995年,咸阳市常温区和强绿岛区面积持续减少,且常温区面积减幅最大,累计减少70.07 km2,年均减少10.01 km2.绿岛区、热岛区和强热岛区面积均增加,其中绿岛区面积累计增加41.1 km2,年均增加5.88 km2; 强热岛区面积增幅最小,累计增加2.98 km2.通过转移矩阵分析发现:常温区减少的面积主要流向绿岛区和热岛区,强热岛区面积的增加主要源自于热岛区域.1995~2006年,常温区和强热岛区面积持续增加,其中常温区增幅最为显著,累计增加90.31 km2,年均增加10.03 km2.同时,热岛区、绿岛区和强绿岛区均表现为减少的趋势; 热岛区和绿岛区面积减幅显著,9年间热岛区和绿岛区累计减少33.2 km2和57.28 km2,强绿岛区面积基本保持不变.通过转移矩阵分析发现,减少的热岛区和绿岛区均流入常温区.因此,2006年咸阳市整体热环境较1995年呈现恶化的趋势.2006~2019年,除了常温区面积急剧减少外,强热岛区、热岛区、绿岛区和强绿岛区均表现为面积增加的趋势,且热岛区增幅最为显著,13年间累计增加102.71 km2,其次为强热岛区域,累计增加14.91 km2,绿岛区和强绿岛区变幅最小,累计增加面积3.88 km2和1.52 km2.综合分析1988~2019年31年间咸阳市不同热力景观面积变化特征,发现31年间热岛区域显著增加,累计增幅99.25 km2,强热岛区增加18.07 km2; 同时,常温区减幅最为显著,31年间累计减少104.65 km2,变化比占咸阳市总面积的19.28%.进一步分析发现,咸阳市31年间热环境整体表现为加剧的趋势,同时,局部热环境恶化显著.
景观格局既是景观异质性的体现,又是人类活动与各种生态过程在不同尺度上综合作用的结果.城镇热环境景观是城镇与周围环境相互作用而形成的,是人工干预城镇局地生态系统导致局地气候变异的热力学表现[15].本文从斑块、类型和景观三个层次择斑块总数(NP)、景观形状指数(LSI)、周长-面积分维数(PAFRAC)和聚集度指数(CONTAG)4个指标[16]用以分析咸阳市1988~2019年31年间热力景观格局的演变特征.基于ENVI5.2平台,采用中误差划分的方法将研究区热力景观类型划分为五大类,利用ARCGIS10.5平台将分类后数据转化为ARCGRID格式,最后利用景观分析软件FRAGSTATS4.2分别计算了不同时相反演结果的景观指数(表3,表4).
图2 1988-2019年咸阳市不同热岛类型数量变化图
Fig.2 The changes of the different urban heat island types of Xianyang City from 1988 to 2019
NP指数变化分析:咸阳市强绿岛区斑块数1988~1995年呈现小幅增加的趋势,1995~2006年呈现较大幅度增加的趋势,11年间斑块数累计增加21个; 2006~2019年,强绿岛区NP个数呈现减少的趋势,13年间累计减少12个.分别提取强绿岛区并利用ArcGIS空间叠加分析功能,对比分析1988年和1995年NP变化发现,该时间段增加的强绿岛斑块主要集中在渭河沿岸地区.绿岛区斑块数1988~1995年呈现减少的趋势,7年累计减少270个.1995年之后呈现较大幅度增加的趋势,截至2006年,绿岛区斑块数较1995年增加459个,2019年绿岛斑块数较2006年呈现小幅度减少的趋势,13年间累计减少32个.分析原因主要是由于城市的快速发展,使得具有高吸热的建设用地取代自然地表覆盖所致.常温区2006~2019年13年间,NP指数呈现先增加后减少的趋势,1988~2006年,该时间段常温区NP指数增加明显,累计增加斑块697个,主要是由于土地利用方式的改变所致,通过空间叠加对比发现,变化最为剧烈的区域主要集中在城市周边地区.2006年之后,常温区斑块数表现为减少的趋势,13年间累计减少147,分析原因在于耕地等较低热熔的土壤覆被类型被水泥、沥青等高吸热材料代替所致.热岛区斑块指数呈现先持续减少,后较大幅度增加的趋势,且变化的区域主要集中在咸阳市的东北方向.同时相比1988年,高温区变化主要呈现基于原热区域的片状扩展与延伸.强热岛斑块呈现先增加后持续减少的趋势,1988~1995年累计增加199个,之后截止到2006年,NP指数累计减少124,之后表现为增加的趋势,13年间累计增加强热岛斑块63个.分析原因在于土地利用类型的转变导致的区域地表热辐射平衡被打破所致.
(2)LSI变化分析:咸阳市强绿岛区景观形状指数1988~1995年7年间持续减少,1995~2006年11年间又急剧增加,之后2006~2019年13年间较小幅度减少.表明强绿岛区斑块形状在城市化的不同阶段整体趋于简单化.绿岛区1988~1995年7年间斑块复杂度总体表现为降低的趋势,累计减少11.8,表明在改革开放不久后的咸阳市该时间段绿岛区景观形状趋于规则化,与当时的农业经济发展模式相一致.1995年后绿岛区景观形状指数呈现持续上升的趋势,至2019年累计增加9.6,表明自1995年后市场经济发展以来,咸阳市绿岛区景观形状伴随经济市场复苏与多元化发展模式,绿岛区斑块逐渐趋于破碎化和复杂化,这与城市在这一时间段开始规模化建设起步相一致.常温区1988~2006年18年间斑块形状指数持续变小,表明常温斑块在这一时期趋于复杂化; 2006年后则呈现规则化特征,前者是因为组成常温地表覆被不断遭遇无序建设与开发的结果,而后者则是城镇化体系及其空间格局基本形成,同时城镇化建设更加注重规划、有序开发的结果.热岛区景观形状指数1988~2019年31年间呈现先小幅增加,后大幅度降低,降幅为13.1,之后又增加的趋势,分析原因在于2002年西咸一体化政策实施后,城市开发强度和速度不断加快,城市地表覆被急剧变化,且已有的城镇空间格局中城镇生态建设中的公园、绿地等低容物质的大量出现,使得城镇热岛区域形状复杂度降低后又增大的趋势.热岛区景观形状复杂度指数呈现先升高,后降低又持续走高的趋势.分析原因一方面是由于城市组团类型建筑在2000年后不断涌现,另一方面城市开发强度、空间和格局更加紧凑的结果.
(3)PAFRAC变化分析: 31年间咸阳市强绿岛区斑块分维数呈现小幅降低的趋势,表明强绿岛区的斑块趋于规则化,主要原因是城市进一步开发建设使得城市水体及高密度绿化区呈现单一化、有序化.绿岛区和强热岛区呈现较为相近的变化规律,即先持平后小幅增加.常温区和热岛区则整体表现为小幅增加的趋势,表明改革开放经济快速增长的同时,热斑块复杂度增加,城市土地利用异质性加剧; 2006年后又表现为小幅减少的趋势,表明西咸一体化进程不断加快过程中,咸阳市通过合理的城市规划、增加绿化及水域面积等相关措施使得常温区域斑块趋于规则化.
(4)CONTAG变化分析:咸阳市景观聚集度指数由1988年的51.224减少到1995年的49.651,表明咸阳市7年间土地利用整体趋于复杂化与破碎化,进而导致许多离散小斑块的产生,究其原因是大面积优势斑块被不断扩张的离散型城镇建设用地不断替代和冲击的结果; 到2006年CONTAG指数增加到54.918,较1995年增加了5.267; 至2019年景观聚集度指数则高达59.254.1988~2019年31年间CONTAG指数整体增加8.03,表明咸阳市热环境景观格局趋于整体化,究其原因是咸阳市大规模建设过程相对完成,地表覆被格局已基本形成所致.
为了进一步剖析不同时间节点,咸阳市强热岛区和热岛区两类热力景观的空间演化过程,基于GIS提供的度量地理分布功能,采用中心要素和利菲弗方向性分布算法,刻画不同时间节点的热环境区域空间位置及伴随时间的推移轨迹,结果见图3所示,用以分析随时间变化和城市不断发展,城市强热岛区和热岛区两类热力景观的空间演变特征.
图3 强热岛区(左)和热岛区(右)1988~2019年重心分布空间移动轨迹图
Fig.3 Spatial movement track of gravity center distribution in strong heat island area(left)and heat island area(right)from 1988 to 2019
由图3可知,1988~1995年7年间强热岛区中心分布呈现由西向东的趋势,之后向咸阳市东北方向偏移,2006~2019年13年间强热岛区中心点位继续向东北方向偏移,且该时间段偏移量最大.结合咸阳市地图可以看出,2019年咸阳市强热岛区的中心点位已经由原咸阳市旧城区开始转向空港新城方向,特别是2010年以后空港新城的高强度开发与建设已经导致该区域成为新的强热岛区域,由于该区域的存在与影响,整体导致咸阳市强热岛区域的空间走向偏向东北方.对于热岛区而言,1988~1995年7年间咸阳市热岛区域中心偏向东南方向,之后偏向北方,2006年以后,热岛区域继续偏向东北方,但增加幅度小于强热岛区,进一步分析发现,强热岛区和热岛区的空间发展方向相一致,但强热岛区比热岛区的发展空间跨度更大,表明在空港附近由原有的低热容的地表覆盖被高热熔的水泥、柏油等下垫面所取代,具有高热容吸热、储热能力强的建筑材料的广泛使用,使得咸阳市强热岛区和热岛区变化显著.以标准差椭圆走势可以看出,无论是强热岛区还是热岛区,它们的宏观走势均由西南方向走向东北方向.且热岛区域在偏向东北方的同时,更多的偏向东方,即西咸一体化所在的西咸新区.由此表明,咸阳市空间热力景观格局随着西咸新区、空港新城等的不断建设与完善,已经形成了新的热力格局.认识和了解咸阳市热力中心的动态演变有利于后期合理规划和制定相应的缓解城市热环境的策略.
本文以1988、1995、2006和2019年四个时相的热红外波段为主要数据源,基于单窗算法和中误差划分方法得到了不同时相地表温度数据,归一化处理后,借助景观生态学理论和地理空间度量方法对咸阳市1988~2019年31年间的热力景观格局变化特征进行分析,得到以下结论:
(1)咸阳市1988~2019年31年间不同热力景观类型面积大小整体表现为:常温区>绿岛区>热岛区>强热岛区>强绿岛区; 2019年常温区和绿岛区占地面积较1988年减少了21.75%,且主要流向热岛区和强热岛区,表明咸阳市热力环境呈现恶化的趋势.
(2)常温区是咸阳市1988~2019年31年间热力景观斑块变化数量最多的类型,累计减少面积104.65 km2; 其次是热岛区,31年间累计增加99.5 km2; 变化幅度最小的是强绿岛区,31年间累计减少了2.24 km2.经过对比分析不同热力景观格局的演变过程,发现热岛区增加的面积主要源自于常温区.
(3)热力景观指数分析表明:1988~2019年31年间咸阳市常温区热力景观类型占主导地位; 景观斑块形状复杂度整体降低,这与周长~面积分维数指数和聚集度指数反映结果相一致.分析其成因最根本的是城市化进程导致城市用地类型占主导地位及城镇化大规模建设过程基本完成所致.
(4)咸阳市热岛景观和强热岛景观整体呈现由西南向东北方向的走势.其中强热岛区中心点位的走势受空港新城版块的建设显著,而热岛区走势受西咸一体化政策影响显著.