1. 引言
城市热环境作为城市气候的显著特征,被认为是主导城市生态环境的重要因素之一 [1] [2] 。热岛效应是指由于城市化所引起的城市地表及大气温度高于周边郊区的现象,是城市气候最显著的特征之一 [3] [4] 。引起热岛效应的原因有很多,如地表的热辐射特性差异,不透水层下垫面的蒸腾作用,生产生活中人为的热排放等等 [5] 。如何定量监测城市热岛的动态变化及其特征分析已经成为城市生态环境研究及合理城市规划的重要内容。
城市热岛研究由观测资料的不同可分为基于台站气温资料和基于遥感地表温度数据 [6] [7] [8] [9] 。台站气温资料具有时效长、准确性高的优势,一直被作为研究城市热岛时空演变规律的重要手段。除了常规的台站气温资料,更加精细化的自动站气温资料也被用于研究城市热岛的更多细节特征。目前国内外已经开展许多长时间大范围的城市热环境研究 [10] [11] ,但针对城市内部精细化的热岛日变化研究还较少。王艳姣(2009) [10] 利用TERRAPMODIS遥感反演的地表温度资料,对2000~2006年北京城市热岛季节变化特征进行了研究表明:热岛强度的季节性差异主要与太阳辐射强度、地表植被覆盖状况和城市人为热释放等的季节性变化密切相关。张硕(2017) [11] 通过改进乡村背景选择算法建立一种在区域或全球尺度上相对简单合理的SUHI估算方法,并引入热岛面积、热岛比例指数等指标开展珠三角城市群热岛时空分布特征分析及定量评估。
本研究利用加密自动站气温资料分析厦门精细化城市热岛的日变化特征,以期加深对厦门城市热环境的认识,为城市规划、城市气候承载力分析提供一定的科学依据。
2. 资料与方法
2.1. 数据资料
1) 采用由国家气象信息中心提供的逐小时加密台站气温资料,厦门市内共36个台站数据,以2016年9月1日为例分析00时至23时逐小时气温变化及城市热岛日变化特征。
2) 利用2010年土地利用/土地覆盖资料来划分厦门郊区范围,该资料来源于欧空局CCI (Climate Change Initiative)项目,空间分辨率为300 m。本研究结合厦门市行政区划范围裁切出厦门市土地覆盖数据(如图1),由图1可知,厦门市土地利用/土地覆盖类型主要有旱作农田、灌溉农田、常绿阔叶林、落叶阔叶林、常绿针叶林、灌木丛带、草地、城市、水体。
2.2. 研究方法
2.2.1. 气温插值方法
本研究在对加密台站气温资料进行空间插值时,采用普通克里金(Ordinary Kriging)插值方法,将台站气温插值成空间分辨率为1 km的格点数据。克里金法由Matheron (1971) [12] 提出,是从变量相关性和变异性出发,对区域化变量的取值进行无偏、最优估计的一种统计方法,该方法的最佳适用条件是变量存在着空间相关性。普通克里金(OK)是应用最广、最普遍的克里金方法之一,其计算公式为:
其中,
为待估计点的值,
为待估值点周围的点(xj, yj)处的变量值,
是其对应的权重系数,该系数是通过对样本变异函数计算值进行模拟得到的,样本变异函数值计算式为:
其中h为分离距离,Nh是在(xi + h, xi)之间用来计算样本变异函数值的样本点对数,变异函数的理论模型有球状模型、指数模型、高斯模型等,本研究采用球状模型。
2.2.2. 城市热岛指数计算
为了消除背景场对分析厦门城市热岛24小时变化的影响,本研究对空间插值后的格点气温数值作如下热岛指数化处理。按照多数文献研究城市热岛的方法,城市热岛强度定义为城市和乡村气温之差 [13] ,因此定义热岛指数:
其中Mi为热岛指数,Ti为气温有效像元值,Tj为郊区气温有效像元值,n为研究区有郊像元个数。参考叶彩华等(2011) [13] 提出的热岛强度等级划分标准,并依据热岛指数值的大小进行如下划分:>3,强热岛;2~3,次强热岛;1~2,弱热岛;−1~1,过渡区;−2~−1,弱低温;−3~−2,次强低温;<−3强低温。由土地利用/土地覆盖数据划分厦门郊区范围,将地表类型为常绿阔叶林、常绿针叶林的区域作为郊区。经过处理后,研究区每个格点的气温值就转化成一个热岛指数,有利于不同时刻气温之间直接比较,可较客观地分析热岛日变化特征。
3. 结果与分析
3.1. 气温日变化特征
经过空间插值后的逐小时00时至23时格点气温时空分布如图2,其中,由于14时台站气温缺测数值较多,故将其剔除不参与计算分析。由图2可知,厦门市24小时气温呈现明显日变化特征,从00时气温逐渐升高,在13时气温值达到最大,00时至02时气温比21时至23时气温低。从厦门市气温空间分布上看,整体呈现东南向西北方向逐渐降低,温度较高的区域是集美区、同安区。
Figure 1. Distribution map of land cover type in Xiamen
图1. 厦门市土地覆盖类型分布图
Figure 2. Hourly temperature distribution of Xiamen from 00:00 to 23:00 on September 1st
图2. 厦门地区9月1日00时~23时逐小时气温分布图
3.2. 城市热岛时空演变
厦门城市热岛时空演变如图,整体呈现东南至西北方向城市热岛效应逐渐递减的空间格局。如图3可知,白天气温较高的时候即09时至16时,城市热岛效应较弱,且没有出现强热岛效应。而厦门热岛主要在夜间,热岛强度最大出现在夜间00时,该时刻强热岛影响范围主要在厦门岛内北部、海沧区、集美区、翔安区的环岛区域。
3.3. 热岛日变化曲线
对逐小时城市热岛日变化进行统计分析,将强低温、次强低温、弱低温合并为低温区,过渡区、强热岛与次强热岛、低温区的面积百分比日变化曲线如图4所示。由图可知,强热岛与次强热岛的面积百分比在10时至16时较低,而00时面积百分比最大为46.43%,01时至08时、17时至22时该两个时段内,强热岛与次强热岛的面积百分比在20%至30%之间。低温区的面积百分比变化不大,平均值约为8%。过渡区的百积百分比在15时至16时较大,最大值为75.82%,而其他时段内变化不大,平均值约为30.97%。
Figure 3. Spatial distribution of hourly urban heat island index in Xiamen area from 00:00 to 23:00 on September 1st
图3. 厦门地区9月1日00时~23时逐小时城市热岛指数空间分布图
Figure 4. Daily percentage change of urban heat island area
图4. 城市热岛面积百分比日变化曲线图
4. 结论与讨论
由于本研究以城市热岛日变化的具体个例进行分析,结论的代表性有所欠缺,需要开展不同季节的平均态的城市热岛日变化分析;另外本研究只采用台站观测资料研究气温热岛,需要开展基于静止卫星的逐小时地表热岛的分析研究。综上,本研究探讨了建立精细化城市热岛网格气候信息的方法,以期为城市生态环境、城市气象服务体系建设提供初步分析研究基础。
基金项目
国家自然科学基金项目(41705045)、厦门市科学技术局科技惠民项目(3502Z20174051)、厦门市科学技术局科技惠民项目(3502Z20164083)、中国气象局预报员专项(CMAYBY2017-036)、福建省强降水机理研究创新团队、厦门市气象局海洋精细化预报创新团队共同资助。
参考文献