齐齐哈尔高程图见 图1 ，齐齐哈尔市位于中国东北地区，地处松嫩平原，总面积达到42,469平方千米。这座城市聚集了汉、满、蒙古、朝鲜等多个民族，总人口为516.5万人。齐齐哈尔市拥有丰富的自然资源，其土地肥沃，适合农业发展。在社会经济方面，齐齐哈尔市是国务院批复确定的中国重要的工业基地和商品粮基地，也是黑龙江西部中心城市 [4] 。近年来，该市的经济持续增长，2023年实现了地区生产总值1327.6亿元，比上年增长1.6%。其中，第一产业、第二产业和第三产业均有所增长，显示出多元化的经济发展趋势，此外，齐齐哈尔市还拥有丰富的历史文化遗产和自然景观，如扎龙自然保护区等，为城市的发展注入了独特的文化内涵和生态价值 [5] 。

本文所需数据由GIS平台进行处理，所用数据来源简要说明见 表1 。

本文以齐齐哈尔市为研究对象，基于环境学、景观生态学及可持续发展等理论，依托于ArcGIS、ENVI

等平台，利用MODIS地表温度产品数据、NDVI等数据，研究齐齐哈尔市蓝绿空间格局时空演化过程中的热岛效应变化，并针对齐齐哈尔市热岛效应提出相应的管理对策。

土地利用数据集采用武汉大学的杨杰和黄昕教授在GEE上制作了第一个源自Landsat的中国年度土地覆盖数据集(CLDC)，其精度、数据整体质量良好。土地利用类型共9类，分别为农田、森林、灌木、草原、水域、冰雪、裸地、不透水面和湿地 [6] 。根据分类基础，利用ArcGIS平台将农田、森林、灌木、草地合并为绿色空间，将湿地、水体合并为蓝色空间。由于裸地和某些不透水下垫面(如水泥)在夏半年日最高地表温度差异不显著，故将不透水表面和裸地合并为不透水下垫面空间。

由于水体在近、中红外波段极高的吸收率以及微波波段极低的发射率，所以利用红外波段对水体提取，具体公式如下 [7] ：

$\text{MNDWI}=\left(\text{Green}-\text{MIR}\right)/\left(\text{Green}+\text{MIR}\right)$(1)

地理探测器是探测空间分异性、揭示其驱动因素的一组统计学方法，能够探测单一驱动因子以及两驱动因子的交互作用对因变量空间分异性的影响，并检验其统计显著性 [8] 。本文运用地理探测器交互作用探测，探究NDVI、MNDWI两因子的相互作用对地表温度的影响。

1) 分异及因子探测：探测因变量Y的空间分异性以及探测某一因子X对因变量Y空间分异的解释程度，用q表示。公式为：

$q=1-\frac{{\displaystyle \underset{h=1}{\overset{L}{\sum }}{N}_h{\sigma }_h^2}}{N{\sigma }^2}$(2)

2) 交互作用探测：通过计算某两个不同因子X1和X2对因变量Y解释力q(X1)和q(X2)、交互作用时对Y的解释力q(X1∩X2)，并比较q(X1∩X2)、q(X1)、q(X2)、min(q(X1)、q(X2))和max(q(X1), q(X2))的大小，探测X1和X2的交互作用对因变量Y空间分异的解释力是否减弱或者增强 [9] 。影响因子的交互作用类型见 表2 。

通过2012年、2015年、2018年2021年的土地利用数据，构建出齐齐哈尔市的蓝绿空间格局，结果见 表3 。如 表3 所示，齐齐哈尔土地利用类型主要以绿色空间为主，其次为不透水面空间，而蓝色空间不足总面积的2%。随着齐齐哈尔的发展，不同的空间类型呈现出各自独特的变化趋势。

齐齐哈尔市2012年、2015年、2018年、2021年蓝绿空间格局图可见 图2 。如 图2 所示，在研究期间内，不透水面空间的持续增长和绿色空间的轻微缩减，共同揭示了齐齐哈尔的城市化进程仍在快速推进中。具体而言，过去十年中不透水面空间的占比持续增加，从2012年占比4.68%至2021年占比6%。这一显著的增长证明了齐齐哈尔的迅速发展，更凸显了其城市化水平的持续提升。与此同时，虽然绿色空间的面积减少了1.42%。另一方面，蓝色空间面积呈现波动上升，这一变化不仅彰显了生态政策在平衡城市化步伐与改善城市风貌方面的积极作用，更意味着未来城市发展中，生态保护与可持续发展将成为不可忽视的重要考量。

2012~2021年齐齐哈尔地表温度可见 图3 。由 图3 可知，地表温度的统计结果表明，齐齐哈尔市温度整体呈现增长趋势。具体而言，2012年地表温度高值区主要分布在西南部；2015年研究区高温趋势向北部、东北部逐渐扩散，并在北部地区逐渐增强；从2018年开始，齐齐哈尔市地表温度最高值虽有所上升，但热分布区域减小，主要集中于西南部地区；2021年研究区最高温度和最低温度均有所上升，高温地区主要集中于西南部地区。总体来看，从2012年到2021年，齐齐哈尔的地表温度高值区经历了由扩散到聚集再到增强的变化过程，热岛效应体现明显。

齐齐哈尔市因子探测结果见 表4 。如 表4 所示，通过地理探测器分析、识别影响齐齐哈尔热环境空间分异的主要驱动因子，包括NDVI和MNDWI。分析结果显示，虽然MNDWI的影响因子略有下降，

但MNDWI、NDVI均对齐齐哈尔市地表温度具有一定的解释力，表明研究区内2012~2021年间齐齐哈尔蓝色空间、绿色空间均能在一定情况下对地表温度产生显著影响。且NDVI∩MNDWI均大于0.5，说明NDVI和MNDWI间相互作用对地表温度解释力较高，表明绿色空间和蓝色空间相互作用也能很好地影响齐齐哈尔的地表温度。