全球城市人口将在2000年30亿的基础上, 在2050年增加到64亿^[1], 全球一半以上的人口将受到城市热岛效应的影响^[2]。1833年英国气候学家Lake Howard对伦敦城郊温度进行比较, 发现城市中心的温度明显高于郊区^[3]; 1958年, Manley将这种现象形象称之为“城市热岛效应”^[4]。城市热岛效应与许多因素有关：一方面人类活动消耗的各种形式的能量释放到大气中, 这部分热量作为显热直接导致城市温度上升; 另一方面, 城市化带来的土地覆盖从自然景观(如自然植被和农业用地)向不透水表面(如人工材料, 主要是混凝土、沥青、瓷砖、金属)的转变, 导致城市表面的辐射率、热、湿度、粗糙度和发射率的变化^[5]。随着显热通量与潜热通量之比的增加, 导致城市地表及其上层的温度升高^[6]。城市环境中气溶胶或霾污染的集中增强了影响大气和地表之间能量交换的长波辐射强迫^[7]。所有这些能量收支的变化都会影响城市热岛效应的表现。

各位学者利用气象数据和遥感数据从全球尺度^[8-10]、国家尺度^[11-17]、区域尺度^{[16, 18-19]}、城市尺度^[20-27]对城市热岛的时空特征以及影响因素^{[5-7, 10, 16-17, 20, 28]}、缓解措施^[28-30]等方面开展了大量工作。基于气象监测数据研究的大气城市热岛效应, 它具有较高的时间分辨率, 能较好的研究城市热岛的时间演变规律; 根据热红外卫星遥感数据研究的地表城市热岛效应, 相对于气象监测数据, 它具有更好的温度空间分辨率。Jones等人收集了中国多个气象站点多年的气温数据, 发现城区平均气温比乡村的平均气温高0.39℃^[12]。王媛媛等利用地面观测数据作为辅助手段研究了我国城市热岛和空气质量的时空变化及影响因素^[11]。在1972年Rao就利用遥感数据来反演地表温度并应用于城市热岛的研究^[31], 城市热岛效应具有明显的时空变化特征。Santamuris等调查了亚洲、大洋洲的城市热岛效应均发现具有明显的季节特征^[9]。Rui等研究结果显示, 中国大部分城市的城市热岛强度和热岛面积在夏季白天明显增加, 而在夏夜和冬夜分别有大约一半和一半以上的城市显著增加^[13]。Penge等对全球419个大城市的地表城市热岛效应研究发现, 白天地表城市热岛效应高于夜间, 两者没有相关性^[8]。对于我国城市热岛效应规律, 从昼夜尺度来看, 整体大部分城市热岛强度夜晚高于白天, 并且具有一定空间分布特征。如Zhou等研究了我国32个主要城市2003—2011年间的昼夜地表城市热岛效应的时空分布格局, 结果表明夜间地表城市热岛强度高于白天, 且具有较大的空间异质性^[14]。从季节尺度来看, 全国整体夜间冬季城市热岛强度最强, 但因研究方法的不同, 存在各城市热岛强度的昼夜及空间差异。如王媛媛等对全国不同区域的四季昼夜UHI进行统计, 我国白天UHI夏季最强, 冬季最弱, 与白天相反, 夜间UHI冬季最强, 夏季最弱^[11]; 不同于先前关于美国大陆的城市热岛研究, 在美国, 夏季白天和夜晚的UHI均最强^[15]。从空间尺度来看, 我国城市热岛效应的昼夜及季节特征因不同地区气候背景有所差别; 曹畅等按湿润区划分气候区研究中国城市热岛效应的时空分布规律发现各季节、昼夜热岛强度存在明显的气候区差异^[16]。影响城市热岛的因素, 首先是受到地理位置决定的气候条件为主导, 其次是受到来自城市化与工业化进程的影响。自然因素包括地理位置、海拔、植被覆盖以及与地表覆盖有关的因素。孙然好等研究了我国245个不同规模的城市局地气候背景对城市热岛效应日变化的影响, 其解释率分别为春季30%、夏季19%、秋季29%、冬季25%^[17]。金珊合等人运用局部气候带分类, 研究大连市局部气候带分区类型和特征, 进而分析不同局部气候带对城市地表温度的影响^[20]。

综上研究发现, 同一城市热岛效应存在季节性、昼夜性的差异, 不同城市也存在因自然地理背景以及城市化等因素导致的热岛强度差异; 而我国地域辽阔, 地形复杂, 气候类型多样, 南至热带/亚热带, 北至亚北极^[32], 地处东亚季风区, 具有由南至北逐渐减小的温度梯度和由东南至西北逐渐减小的降水梯度^[33], 因此, 对于全国尺度的城市热岛效应研究有必要在不同气候带进行量化和分析, 有助于加深对城市热岛的演变认识。

1 数据与方法 1.1 数据与方法

利用1991—2019年全国各城市逐年气象数据、以柯本气候分类法为标准, 将我国划分为五个气候带, 研究我国各个气候带的空间分布特征以及时间演变特征(图 1), 所选取的城市通过在Google Earth中对城市郊区气象站点相对位置、站点下垫面的土地利用/覆盖类型进行验证, 筛选出各城市站点背景值最佳的1个市区站点和1个郊区站点, 最终确定了全国适合该研究的69个城市。选取的市区郊区站点是以2020年全国气象站点位置和Google Earth实时地图对照选取, 由于城市化进程一般是不可逆的, 所以对于2020年以前的郊区站点也不会有城市化的影响; 另外, 对于市区站点位置选取的均是城市中心区域, 在研究期间不存在市区站点由郊区转化而来的情况, 仅是一个城市化不断发展的过程, 这正好与本研究中城市热岛时间演变过程相契合。选取时间分辨率较高的气象数据作为热岛效应的量化依据, 数据均来自于国家气象科学数据中心(http://data.cma.cn/)。

1.2 柯本气候带分类

根据表 1的分类方法^[34], 通过最冷月平均气温、最热月平均气温、年降水量、年平均气温的计算, 将中国划分为热带、干带、冷温带、暖温带、极地带5个气候带(图 1)。本研究中所涉及的69个城市是以各个城市城郊气象站点1991—2019年累年平均气温、降水数据作为依据进行的气候带分类, 中国其余地区的气候带类型以Google Earth全球累年1986—2010年柯本气候带分布区划图(空间分辨率为5m, 更新时间为2017年)作为对照进行划分。

1.3 大气城市热岛的量化 1.3.1 城市热岛强度

数据为1991—2019年各站点逐年月均温, 考虑到市区、郊区站点海拔差异的影响, 在进行热岛强度量化时, 同时做了海拔温度校正, 即海拔每上升100m气温下降0.6℃; 各城市的城市热岛强度以市区、郊区海拔温度校正温差作为量化标准。其计算方法如公式(1)。

(1)

式中, UHII(Urban Heat Island Intensity)表示城市热岛强度(℃); T_u 表示市区站点温度(℃); T_r表示郊区站点温度(℃); E_u表示市区站点气压传感器海拔高度(m); E_r表示郊区站点气压传感器海拔高度(m); 其中, 若UHII为负值, 视为城市冷岛效应UCI(Urban Cold Island), 其强度以城市冷岛强度表示UHII(Urban Cold Island Intensity)。

1.3.2 城市热岛效应的季节差异

分别计算各城市春季(3月、4月、5月)、夏季(6月、7月、8月)、秋季(9月、10月、11月)、冬季(12月、1月、2月)累年平均热岛强度, 并将热岛强度划分为：强热岛(>2℃)、中等热岛(1—2℃)、弱热岛(0—1℃)、强冷岛(< -2℃)、中等冷岛(-2—-1℃)、弱冷岛(-1—0℃)六个热岛强度等级, 最终空间分布结果在ArcGIS10.2中呈现。另外计算不同季节中各气候带城市热岛强度的平均值, 分别讨论各季节中城市热岛效应的气候带差异以及各气候带中城市热岛效应的季节性差异。

1.3.3 城市热岛强度的时间趋势

在Origin 2018中绘制各气候带平均城市热岛效应的逐年变化趋势图。首先进行各气候带、各季节逐年平均城市热岛强度散点图的绘制, 再利用SG滤波法(Savitzky Golay Filter)进行平滑信号处理^[35]。SG滤波法的核心思想是对窗口内的数据基于局部多项式最小二乘法进行加权滤波拟合, 是由Savitzky和Golay在1964年首次提出^[36]。它的优点在于, 在滤波平滑的同时, 能够更有效地保留信号的变化信息, 在对城市热岛逐年变化趋势的研究中能在最优情况下拟合出变化的连续性。SG滤波法的模拟过程可用公式(2)表示：

(2)

式中, U_j^*为SG滤波法平滑处理的UHII数据; U_j+i为各年UHII原始数据; W_i为滤波系数, 即第i年UHII原始数据U_j在平滑窗口中的权重, 由给定的高阶多项式进行最小二乘拟合得到; N为滑动窗口的数据个数(2m+1), 其中m为窗口宽度。SG滤波函数的窗口宽度m值越大, 滤波拟合越平滑; 多项式阶数越低, 滤波越平滑, 但也会带来更高的误差, 阶数较高可降低此误差, 但“过拟合”会带来噪声影响平滑处理结果^[37]。为了确定SG滤波在不同气候带各季节逐年变化趋势拟合中的最优参数, 本文按照宽度3、5、10、20, 及阶数2、3、4组成12组参数组合进行SG滤波拟合, 经实验最终选取平滑窗口宽度m为10、阶数2为拟合效果最佳参数进行平滑滤波处理。

2 结果与分析 2.1 城市热岛的区域差异

根据前述各城市、各季节累年城市热岛强度的量化方法, 经ArcGIS10.2软件绘图得到全国各城市1991—2019年城市热岛效应季节性空间特征分布图(图 2)。我国城市热岛效应具有明显的气候带差异性, 干带地区四季均较强, 且明显高于其他气候带; 冷温带京津冀地区各季节均较高, 温暖带西南地区四季均较强, 热带、极地带四季均较弱, 且极地带有较强的冷岛效应出现。统计了各季节全国城市热岛效应及冷到效应最强的城市及气候带(表 2)。从气候带来看, 干带平均城市热岛效应在四季均最强, 极地带平均城市热岛效应四季均为最弱(冷岛效应最强)。从全国各个城市来看, 春、夏、冬、年平均热岛强度最高的城市均处于干带地区; 春、夏、秋、冬、年平均热岛效应最弱(冷岛效应最强)的城市均处于冷温带地区。以往研究表明干带地区大部分城市由于郊区植被稀少, 多为反照率更低、吸热能力更强的裸土或荒漠, 城市热岛效应一般较弱、或表现为冷岛效应^{[16, 23, 38-41]}, 这与本文研究结果差异较大。主要原因是由于先前研究结果反映的是地表城市热岛效应, 而本研究反映的是大气城市热岛效应。如干带地区的哈密、库尔勒等城市的市区下垫面和建筑表面反照率高于郊区荒漠, 郊区地表白天更易吸收热量, 从而表现为郊区地表温度高于市区。但对于大气温度而言, 市区人口聚集、建筑密集, 大量的热排放至大气中却无法散失; 而郊区多为地表开阔且吸热能力强的裸土或荒漠, 大气温度更易通过气流散失或被地表吸收, 从而表现为市区大气温度高于郊区温度。冷温带地区冷岛效应最强的城市为西宁和漠河, 其主要原因为：①城郊气流交换阻碍小, 西宁处于冷温带高原地区、漠河处于亚北极冷温带地区, 气温低且降水少导致适宜生长的植被低矮稀少、建筑物也较为低矮稀疏, 因而对城郊气流交换的阻碍也较小; ②市区人为热排放少、散热快, 由于城市规模小、建筑密度低、市区人为热排放也就较低, 加之市区站点海拔较高、风速较大也进一步提高了市区的散热效果; ③市区地表及建筑材料吸收了大量的大气温度, 地表及建筑材料的导热系数高、比热容小, 具有在高温中升温快、在低温中降温快的特性, 因此在寒冷地区反而能吸收大气中的热量。因此城市化规模较小、寒冷、降水较为缺乏的冷温带城市更易表现为大气冷岛效应。

对各季节不同气候带热岛强度进行量化, 讨论各季节中城市热岛效应的气候带差异(图 3、表 3)。城市热岛效应的气候带差异性全年整体表现为干带最强; 极地带最弱, 且表现为较强的冷岛效应。各季节中气候带热岛效应强弱依次为, 春季：干带>热带>冷温带>暖温带>极地带, 其中, 极地带有较强的冷岛效应; 夏季：干带>热带>极地带>冷温带>暖温带, 各气候带均出现热岛效应; 秋季：干带>热带>冷温带>暖温带>极地带, 其中, 极地带有较强的冷岛效应; 冬季：干带>冷温带>暖温带>热带>极地带, 其中, 极地带表现为较强的冷岛效应, 热带表现为较弱的冷岛效应。

2.2 城市热岛的季节性差异

在不同区域气候背景下量化了各季节城市热岛效应(图 4)。其强弱依次为, 热带：夏季>秋季>春季>冬季, 其中冬季出现较弱的冷岛效应; 暖温带：夏季>秋季>春季>冬季, 全年各季节均出现热岛效应; 冷温带：春季>冬季>夏季>秋季, 全年均出现城市热岛效应; 极地带：仅夏季出现较弱的热岛效应, 其他季节均出现强冷岛效应, 其中冷岛效应强度依次为, 冬季>春季>秋季; 干带：春季>冬季>夏季>秋季, 各季节热岛强度均为各气候带最高。

对各个气候带不同季节之间的城市热岛强度的差异进行了显著性分析(表 4), 发现各个气候带中不同季节之间城市热岛强度具有显著性差异。如, 热带夏热冬暖、四季降水充沛, 春、夏、秋季气温差异不大, 仅冬季气温相对较低时城市热岛强度才与其他季节具有显著差异性。暖温带四季均较温暖、降水充沛, 各季节气温、降水等气候差异较小, 因此气候带之内的城市热岛效应的季节差异性并不如其他气候带显著。冷温带夏热冬冷, 气候带之内的季节差异性显著, 但其季节差异性并不如其他气候带显著; 其主要原因可能是由于冷温带地区冬季大量的供暖增加了城市能耗, 加剧了冬季原本较低的城市热岛强度, 从而减小了气候带内的季节差异性。极地带四季气温均较低, 一般仅夏季气温在0℃以上, 夏季表现为较弱的热岛效应, 其他季节均表现为较强的冷岛效应; 又因冬季气温极低且缺乏供暖系统、人为热排放较少, 冬季表现为极强的冷岛效应; 从而冬、夏季与其他季节形成了显著的差异性。干带夏热冬冷、城市热岛效应四季均较强; 由于四季降水缺乏、气候常年干燥, 城市热岛强度受降水调节的影响较为直接, 当夏、秋季节降水较多时, 市区降温效果明显, 但在郊区少量的植被也具有一定的恒温调节能力, 导致夏秋季城市热岛明显下降, 从而夏、秋季与其他季形成了显著的差异性。

统计了全国各季节六个热岛强度等级的城市数量特征(表 5), 发现全国城市春、夏、秋、冬四季普遍存在城市热岛效应, 其发生率(热岛城市数量)分别为：70.6% (48个)、76.5% (52个)、63.23% (43个)、60.3% (41个); 城市热岛效应主要处于低热岛强度(0—1℃)范围; 极少数城市在春、夏、冬三季出现强热岛效应(>2℃), 其中冬季出现强热岛效应的城市最多(4个)。另外, 在各个季节也存在部分城市冷岛效应的现象, 主要集中于弱冷岛强度(-1—0℃)范围; 仅冬季有少量城市(4个)出现强冷岛效应(< -2℃)的现象。总体来说, 秋季城市热岛效应较为稳定, 热岛和冷岛效应均主要处于弱强度范围; 春、冬季城市热岛效应较为波动, 易出现强热岛、强冷岛效应的极端现象; 夏季热岛效应发生率最高, 但热岛强度相对较弱。

2.3 城市热岛的时间演变规律

利用各城市不同季节1991—2019逐年热岛强度进行气候带平均化, 并利用SG滤波法进行平滑处理, 得出各气候带不同季节1991—2019年逐年平均热岛强度变化曲线(图 5)。(1)干带：各年份四季几乎均出现城市热岛效应, 各季节热岛效应逐年变化趋于同一趋势, 呈现先升高后降低的起伏状态; 在2010—2015年四季热岛效应均达到较高的状态, 2010年后逐渐降低。热岛强度季节特征2010年前表现为：冬季>春季>夏季>秋季; 2010年后表现为春季>夏季>冬季>秋季。(2)极地带：各年份各季节中仅夏季出现较弱的热岛效应, 其余均表现为冷岛效应, 且冬季冷岛效应最强。各季节热岛效应(冷岛效应)逐年变化趋于同一趋势, 呈下降-上升-下降-上升的波动状态, 到2010年后呈现明显的上升趋势, 且各个季节开始出现较弱的热岛效应。整体热岛强度季节性特征表现为：夏季>秋季>春季>冬季, (冷岛强度季节性特征表现为：冬季>春季>秋季>夏季)。(3)冷温带：各年份四季节均出现城市热岛效应, 各季节热岛效应逐年变化趋于同一趋势, 2005年前较为平稳, 2005年后各季节城市热岛强度逐渐降低。整体热岛强度季节性特征表现为：春季>冬季>夏季>秋季。(4)暖温带：各年份四季节均出现城市热岛效应, 且强度较为平稳; 各季节热岛效应逐年变化趋于同一趋势, 2010年后各季节城市热岛效应逐渐降低。整体热岛强度季节性特征表现为：夏季>春季≥秋季>冬季。(5)热带：除了冬季少数年份出现城市冷岛效应外, 各年份其他季节均出现城市热岛效应; 各季节热岛效应逐年变化趋于同一趋势, 2010年左右以前较为平稳, 2010年以后有上升趋势。整体热岛强度季节性特征表现为：夏季>秋季>春季>冬季。总体来说, 各季节城市热岛效应演变规律在同一气候带都趋于同一趋势, 不同气候带又具有较大差异。2010年前后是各气候带各季节城市热岛效应的变化拐点; 2010年后, 干带、冷温带、暖温带热岛效应均有所下降, 热带、极地带有所上升。

3 讨论 3.1 本研究结果与以往研究的差异

本研究基于全国各城市气象数据揭示了不同气候带各季节大气城市热岛的时空规律, 目的在于针对不同气候带、不同季节的城市热岛效应的缓解措施给予参考。本研究结论与以往研究存在一定的差异, 如本研究中热岛效应各季节最强的城市主要集中于干带地区, 区域主要位于我国的西北地区; 而Zhou等研究了我国32个主要城市热岛效应的空间格局, 发现东南部地区白天较高、北部地区夜间较高, 尤其是气候湿冷的东北地区, 城市热岛效应最强^[14]。其原因主要是两个方面：其一是数据来源的差异, 以往对我国全国尺度或区域尺度的城市热岛效应多是基于遥感数据反演的地表城市热岛效应, 热岛强度受城市下垫面、建筑材料等更直接的影响; 而本文所利用的是气象数据量化的大气城市热岛效应, 热岛强度受气流、空气质量、城市形态等更直接的影响; 影响因素权重不同, 自然结果也就有所差异。另外, 气象数据和遥感数据时间、空间分辨率差异较大, 气象数据有较连续的时间序列, 能更好的研究城市热岛的时间演变规律, 遥感数据具有较好的空间分辨率, 能更准确的量化各城市整体热岛效应, 二者本身量化热岛效应的时空尺度差异较大, 这在以后的研究中需要尤其注意的问题。

3.2 存在问题

以往研究表明, 城市热岛效应在时间尺度上具有季节性、昼夜性差异。而本文所研究的时间尺度仅仅考虑到季节性差异, 量化的是各季节全天平均城市热岛效应, 缺乏进一步对各季节昼夜差异的探讨, 因此在后期研究中有必要对城市热岛效应的昼夜差异做出细致划分。此外, 本研究采用气象站点数据作为研究对象, 由于各气候带中适合作为城市热岛效应研究的城郊站点数量有限, 无法保证各个气候带都有相同且足够统计量的城市数目进行研究。因此, 对于各个气候带的差异分析可能存在一定的不足或偏差。但对于各气候带中的季节差异性研究, 因为各气候带中各城市的不同季节的数据是一一对应的, 且满足统计量分析, 因此对于各个季节的差异分析具有可信性。

3.3 研究启示

对于城市热岛效应的治理, 我们考虑从解决其两个方面的危害出发, 首先是城市热岛效应加剧了高温天气、降低了城市生活舒适度; 另外就是影响了局部小气候, 使气压环流向城市中心汇聚, 导致了如市区内部极端降水事件的增加增强、有毒有害气体向市区中心汇聚、加速光化学污染速率等现象。对于两种危害而言, 不同气候带、不同季节所呈现的作用效果不同。对于不同季节, 如一些城市热岛效应进一步升高了夏季高温天气气温, 降低了城市生活舒适度、增加了制冷能耗, 但在冬季却在一定程度上增加了舒适度、降低了供暖能耗。与此同时, 对于一些如冬季雾霾天气较多的城市, 城市热岛效应也加剧了污染物的汇聚。再者, 对于不同气候带, 如干带地区夏热冬冷、降水较少、植被较为低矮稀疏、自然风速较大, 植被对热岛效应的缓解效果较低, 但市区与郊区的大气热量和物质交换较快, 一定程度上降低了城市热岛效应; 而对于暖温带地区夏热冬暖、降水充沛、植被高大茂盛, 却也在一定程度上降低了城郊的大气交换效率, 加剧了城市热岛效应。因此针对城市热岛效应的缓解必须首先考虑到各个城市的季节特征以及气候地理适宜性, 从而选择各个气候带适宜的植被、建筑材料、建筑物类型等; 其次再是充分考虑城市发展状况及需求、环境污染特征等, 从而做出适合的城市景观格局规划和工程措施。

城市热岛效应时空特征规律是一个受多维因素综合影响的结果, 如气候带背景、季节、昼夜、地形、城市化程度、城市形态等尺度划分差异都有可能导致最终结果的差异。因此, 在今后的城市热岛的研究中首先有必要对城市热岛效应的时空特征规律进行更精细的量化。充分利用气象数据较高时间分辨率的特点, 以及遥感数据较高的空间分辨率特点, 二者相结合, 增强长时间序列的连续性、细化不同时间阶段的分类, 精确大尺度研究的空间特征及其延续性; 其次, 对多种影响因素进行综合考量, 进一步细化研究城市对象的分类, 为城市热岛效应治理应用提供更为准确可靠的依据。

4 结论

基于我国69个城市的市区-郊区气象站点1991—2019年逐年月均温及月均降水数据, 对我国不同气候带城市热岛效应的时空特征规律进行研究, 主要结论：

(1) 区域差异：我国城市热岛效应干带明显高于其他气候带, 极地带最弱, 且表现为较强的冷岛效应, 热带、暖温带、冷温带较弱; 其中, 冷温带京津冀地区各季节热岛效应较强, 温暖带南部地区热岛效应较为明显。全国各季节城市热岛效应最强的城市主要集中于干带地区, 热岛效应最弱(冷岛效应最强)的城市主要集中于冷温带和极地带地区。

(2) 季节性差异：我国各季节城市热岛效应的强弱关系在不同气候带有所不同; 但从全国范围来看, 秋季城市热岛效应较为稳定, 热岛和冷岛效应均主要处于弱强度范围; 春、冬季城市热岛效应较为波动, 易出现强热岛、强冷岛效应的极端现象; 夏季热岛效应发生率最高, 但热岛强度相对较弱。各气候带四季热岛强度表现为, 热带：夏季>秋季>春季>冬季, 其中冬季出现较弱的冷岛效应; 暖温带：夏季>秋季>春季>冬季, 全年各季节均出现热岛效应; 冷温带：春季>冬季>夏季>秋季, 全年均出现城市热岛效应; 极地带：仅夏季出现较弱的热岛效应, 其他季节均出现强冷岛效应, 其中冷岛效应强度依次为, 冬季>春季>秋季; 干带：春季>冬季>夏季>秋季, 各季节热岛强度均为各气候带最高。

(3) 时间演变规律：各季节城市热岛效应逐年演变规律在同一气候带趋于同一趋势, 不同气候带又具有较大差异。2010年前后是各气候带各季节城市热岛效应的变化拐点, 2010年后, 干带、冷温带、暖温带热岛效应均有所下降, 热带、极地带有所上升。