城市热环境是指城市区域内与热有关的、影响居民生存与发展的各种外部要素组成的物理环境系统^[1], 它是城市生态系统最重要的大气环境要素之一, 也是城市环境状况的综合表征^[2—4]。联合国发布的最新报告显示, 到2030年全球城市人口将增加至60%, 而人口超过1000万的特大城市也将从33个增加到43个。大量研究表明, 全球气候变暖和城市扩张联合作用引起的热环境恶化对城市人居环境、生态系统服务乃至人类福祉产生严重影响^[5—8]。随我国城市化进程的加速推进, 城市热环境的时空变化过程及其驱动机制备受学界关注。

在时间与空间维度获得足量的热环境数据是进行城市热环境研究的重要前提。近地面气象监测一直是城市热环境研究最直接、最传统的技术手段^[9—11], 由政府部门管理的城市气象站能提供高时间频率的多种热环境参数, 但是站点较为分散, 无法表达城市内部复杂的微气候状况。为了弥补城市气象站空间分辨率的不足, 学者们一直开展自主布设监测点收集包括气温^[12—15]、风速^{[14, 16—17]}、地面温度^[18—19]和黑球温度^[17]等热环境参数。通过实地监测能够获得丰富的城市热环境信息, 但是设备的采购和维护管理带来的人力和物力成本会大大增加热环境实地监测的难度, 在一定程度上阻碍了热环境研究工作的推进。

近年来物联网和传感器技术的进步为生态环境监测与管理提供了新的理论和方法, 使得生态环境动态监测更具高效性和快捷性^[20]。目前, 国内外学者已在生态环境监测领域逐步使用物联网技术, 如常立侠^[21]提出了我国海岛生态系统物联网监测的概念和模型, 李新^[22]紧密围绕国家生态监测需求, 研发生态系统关键参量监测设备和生态物联网关键技术, Matasov^[23]利用物联网技术开展了城市绿色基础设施水分含量和光谱特征的在线监测工作, Martín-Baos^[24]基于物联网技术开发了低成本交通流量和空气质量指数在线式监测设备, 这些研究成果都为城市热环境的高效监测工作提供了新思路。本研究尝试基于物联网技术构建在线式热环境监测设备, 对城市校园空间热环境状况进行连续高频监测, 深入挖掘微尺度下城市热环境的时空变化特征。

1 研究区与研究方法 1.1 研究区选取和监测点设置

研究区位于广州市番禺区的广州大学校园内, 实验位置地处校园西南角, 面积约0.16km²(图 1)。该位置绿化程度较高, 有大型人工湖, 校园道路主要由沥青和透水铺砖构成, 区域内基本包含了校园典型景观类型, 可作为热环境实验的示范场所。本研究选择了10个热环境监测点, 同时针对不同的科学问题, 在10个监测点上布设了不同的监测设备, 详细描述如表 1所示。

1.2 热环境相关参数获取方法

本研究关注的热环境参数包括太阳辐射、风速、气温和地面温度, 所述地面温度是指沥青和透水铺砖表面的温度平均值, 典型热环境测点的设备布设示意如图 2所示。实验涉及的设备均使用了无线传输技术, 即获取的所有热环境参数会通过公共网络传输至云平台, 用户可以直接在线地从云平台客户端随时实时查看或下载数据。根据不同热环境传感器的供电要求, 本研究设计了相应的在线传输和供电方式, 不同的监测设备数据采集同步进行。

气温和风速：气温传感器测量范围为-40—80℃, 精度±0.3℃, 风速传感器测量范围0—30m/s, 精度0.1m/s。这两种参数使用了窄带物联网(NB-IoT)通讯, 能够实现设备的超低功耗运行, 设备使用了3.6V 20000mAh电池, 数据采集间隔设置为15min, 上传设置为720min, 该设计在理论上能使设备稳定运行4年。

太阳辐射和地面温度：太阳辐射使用光电式太阳辐射计进行测量, 其测量范围为0—1500W/m², 精度为±10W/m²。地面温度监测使用K型热电偶, 其温度测量范围为-200—200℃, 精度为±0.2℃, 地面温度监测是在每个测点上分别测量沥青和透水铺砖两种典型材质, 然后用两种材质的平均温度作为该测点的地面温度。两种传感器使用了12V供电和4G通讯, 为了尽量降低设备功耗和延长电池使用时间, 本研究使用12V 20000mAH锂电池, 同时研发了一种定时开关(图 3), 实现电源开关间隔和通电时长的控制。对太阳辐射和地面温度的采集间隔和上传间隔分别设置均为15min, 即定时开关通电间隔为15min, 通电时间为1min。对比发现, 在没有定时开关, 即两个传感器常开的情况下, 设备只能运行5d时间, 而加入定时开关后同样的采集间隔设置设备可持续运行3个月, 大大提升了热环境监测效率。

1.3 数据处理方法

所有设备在2022年10月15日安装完毕, 并在16日0时开始收集数据。同年10月25日, 从云平台网页端下载数据, 通过01测点的降雨量数据可知, 10月16日至24日内研究区无降雨。对所有获得的数据进行小时平均处理, 随后的分析均在小时平均数据上进行。分析不同测点的热岛效应时, 考虑到测点06的乔木覆盖程度较高, 能代表自然/半自然状态下的热环境状况, 因此选定其作为热岛效应的比较基准, 即将热岛强度定义为其它测点与测点06之间的气温差。此外, 在分析时将7:00—18:00定义为日间, 19:00—06:00(次日)定义为夜间。

2 结果与分析 2.1 太阳辐射和风速变化状况

图 4结果表明, 研究时期内多为晴朗天气, 仅在10月18日、19日出现多云的情况。10月18日研究区太阳辐射量较低, 正午时楼顶的自动气象站01和测点03仅有400W/m², 被树木遮蔽的02测点则更低。其它日期太阳辐射总量逐时变化较为平稳, 具体表现为：正午时刻太阳辐射总量均达到1100W/m²左右, 位于楼顶的01测点正午太阳辐射总量最高；因03测点周边较为开敞, 无任何遮挡物, 其太阳辐射的变化特征与楼顶测点相似；在东南侧被树木遮蔽的测点02也表现出较强的规律性, 即正午之前太阳辐射总量较低, 平均值仅约180W/m²左右, 而在午后, 测点开始受到太阳辐射的直接影响, 辐射总量较高。

图 5为各测点的风速变化情况, 分析结果表明10月16日至18日研究区内出现大风天气, 楼顶自动气象站测得三天内小时平均风速为2m/s, 最高风速5m/s出现在10月17日正午。其它9个测点在这三天也测得较高的风速, 其中测点07和测点09的风速较其它测点高, 原因是它们分别位于相对开敞的人行道上和湖面旁, 通风条件相对较好。10月18日至24日期间各测点风速呈现周期性波动变化, 变化规律相似。

为了进一步分析不同测点之间的风速变化是否存在关联, 对各测点逐时风速值进行皮尔森相关性分析, 结果如图 6所示。结果表明各测点的风速值之间存在较强关联, 相关性系数大部分都在0.5以上。位于楼顶的01测点与各测点均表现出非常强的关联, 相关性系数都在0.7以上(仅与测点06除外, 相关性仅为0.49), 说明楼顶测点能代表研究区的边界条件, 并主导各测点的风环境。同时, 分析也发现, 位置较近或者景观特征相似的测点风速相关性会比较强, 如测点07和测点08的相关性高达0.88, 它们位置仅相隔20m, 因此风速变化的同步性较强, 但又由于测点08旁有树木遮挡, 所以风速值会较测点07低。此外, 测点07与测点09、05之间的相关性也较高, 表明在同一气象条件下它们的风环境变化具有一定的同步性。部分测点之间的相关性相对较低, 如测点10和06与其它测点的相关性大部分在0.5左右, 显示这些测点受研究区边界条件的影响较低, 形成相对独立的小气候。

2.2 日间和夜间热岛强度和风速变化特征

将所有日期的热岛强度和风速数据划分为日间和夜间, 绘制各测点的小提琴变化图, 结果如图 7、图 8所示。从图 7可以看出, 不同测点的风速变化特征差异较大, 07和09测点的风速变化最大, 与楼顶自动气象站的风速变化较为一致, 大部分时间风速集中在0—2m/s左右, 最大风速达到5m/s, 其它测点风速变化较小。对比日间和夜间风速分布特征发现, 日间的风速相对较高, 01测点日间平均风速为1.05m/s, 且有36%的时间风速在1m/s以上, 相比之下夜间平均风速为0.8m/s, 风速在1m/s以上时间只占17%。受到外围边界条件的影响, 其它测点表现出类似的特征。与图 6的结果相似, 测点10和06日间和夜间风速分布类似, 与其它测点的变化不同, 进一步说明了这两个测点的热环境特的独特性。

以自然/半自然状态的测点06为基准点, 统计其它测点日间和夜间热岛效应变化特征, 结果如图 8所示。结果发现, 楼顶01测点所有时间均呈现热岛效应, 日间平均热岛强度为1.73℃, 日间最高热岛强度达4.72℃, 夜间平均热岛强度为0.89℃, 相比之下, 其它测点都出现热岛和冷岛效应。测点05、07和03的热岛效应变化幅度较大, 日间热岛效应都十分强烈, 测点07日间平均热岛强度为0.46℃, 最高热岛强度达2.85℃, 夜间大部分时间会出现冷岛, 最强冷岛达到-1.27℃, 其它几个测点均有类似的变化特征。这些测点在正午左右均受太阳辐射直接影响, 使得日间的热岛效应尤为明显, 到了晚上, 良好的通风条件也促进了近地面的热量散失, 形成冷岛效应。测点10和测点04位置相近, 且周边被建筑物和树木包围, 从图 7看出它们的风速在研究时间内都较小, 因此热岛强度变化较小。测点04受建筑物的阻挡日间没有受到太阳辐射的直接影响, 因此日间有一半时间都表现为冷岛, 到了晚上围合式的建筑也阻碍了通风, 使得日间吸收的热量在夜间没有及时扩散到大气当中, 夜间的热岛效应特别明显。

2.3 地面温度与气温相关系分析

地面温度和气温的关系一直是热环境研究的重要内容, 本研究尝试考虑风速的情况, 分析02和03测点地面温度与气温的关系。同样地将数据划分为日间和夜间, 同时为了更清晰地对比二者关系的差异, 将散点图横坐标和纵坐标设置为同一个边界范围, 结果如图 9所示。结果表明, 日间和夜间地面温度与气温的关系存在明显的差异, 02和03测点的日间温差都比较大, 回归决定系数都在0.8以上, 表示地面温度能表达气温约80%左右的变异程度。02测点决定系数都较高, 日间和夜间分别为0.8082和0.8769, 03测点日间的决定系数较低, 仅为0.7774, 并且数据全部集中在 1:1 直线的下方, 表明日间地面温度普遍高于气温, 这是因为该测点比较开阔, 太阳辐射全天直接加热地面, 使得气温快速升高。上述结果可知, 在有树木遮挡的情况下, 地面由于缺少太阳辐射的直接加热, 气温与地面温度更为接近。夜间两个测点地面温度和气温的关系特征非常相似, 由于夜间没有太阳辐射的加热, 地面温度与气温的关系更为密切, 数据都分布在 1:1 直线附近。此外, 也发现不同大小的风速分布特征不明显, 表明风速对地面温度和气温的关系影响不大, 周边景观特征对二者的关系影响会更大。

3 讨论

基于10月16日至24日的校园微尺度热环境监测结果, 本研究发现, 局部的热环境特征与天气边界条件有直接关系, 但由于景观特征的局部差异, 比如植被条件、建筑物分布、地面状况(邻近人行道或者水体)等, 会直接影响局部热环境特征, 使得小尺度的热环境状况存在一定的时空差异, 表现为热环境的时空异质性。处于楼顶的01测点, 可代表研究区热环境特征的平均状态, 其它测点受所在位置景观特征的影响, 热环境特征会表现不同, 如02测点因树木的遮挡, 下午才能获得更多直接的太阳辐射, 03测点因为周边都没遮挡物, 表现出与01测点相似的太阳辐射变化特征。边界条件的风速状况基本决定了局部位置的风速变化特征, 但具体风速大小差异与树木、建筑等可形成遮挡的景观密切相关。如07测点地处周边无树木遮挡的人行道上, 日间与夜间的风速都比较高, 而距离仅20m远的08测点因为旁边树木的遮挡风速较低。除了风速的差异, 不同树木条件提供的遮阴、蒸散发差异直接影响了07和08两个邻接测点的热岛强度差异。本研究的发现与Tong^[14]的研究相似, 他指出50m范围内绿色植被的比例对局部热环境产生直接影响。刘琳^[17]的研究结果也表明建筑物及树木等遮挡物对区域内部的遮蔽作用会对包括风速、黑球温度在内的热环境参数产生重要影响。因此, 日后的热环境实测工作需要进一步关注微尺度的时空变化特征, 同时也需要在大范围的实测中选取更有代表性的测点。

本研究发现城市热环境在日间和夜间具有不同的时空差异, 研究区日间的风速较高, 夜间风速较低, 热岛强度也是日间高于夜间, 并且日间风速和热岛效应的波动都比夜间大, 因此日间能更加凸显景观类型间的热环境差异。日间在太阳辐射的作用下, 地面直接被太阳的辐射能量加热, 因此地面温度要比近地面空气温度高的多, 建筑和树木的遮挡形成的阴影能极大地降低近地面温度^[25—26]。过往不少热环境实验都只在日间进行, 一方面是实验设备安全管理需要, 二是日间城市热环境的空间异质性会更加明显^[19], 日间也是人们在室外活动的主要时间段, 热环境对人体舒适度的影响更为明显^[27—29]。夜间一般风速都比较低, 大气相对稳定, 在没有太阳辐射影响的情况下, 不同景观表现出来的热环境特征与自身热性质直接相关。有学者研究发现, 日间植被的蒸散发降温作用会被强烈的空气流动所抑制^[30], 因此夜间植被的降温作用可能会更大。因此, 需要开展长时间持续的监测, 在站点尺度上获取更为全面和细粒度的热环境时间变化数据, 并根据日间和夜间不同的热环境形成机理进行有差别的分类探讨, 以全面挖掘城市热环境时空变化过程^{[12, 31—32]}。

在时空维度上获取足量的温度数据是热环境研究的前提, 遥感影像能快速地获取城市大范围全覆盖的热环境数据, 具有较高的空间分辨率, 但时间分辨率较低, 不足以在精细的时间尺度上开展城市热环境研究工作^[33]。而城市气象站只能提供非常稀疏的站点数据, 无法实现细粒度的气温数据获取, 因此使用相关设备开展实地观测是热环境研究的主要手段^{[14, 31]}。许多学者都认为热环境实验需要花费很多的人力和物力^{[16, 19]}, 其核心问题在于使用的设备大部分都是离线的, 即所获取的数据都存放在设备的存储卡内, 科研人员需要花额外的时间精力在线下读取数据才能进行下一步数据分析研究工作。相比于传统的线下数据读取, 本研究利用物联网技术实现了包括气温、风速、太阳辐射、地面温度四种典型热环境参数的在线传输, 设备只需一次安装, 科研人员就能随时随地通过网络查看和下载数据, 在大大降低了热环境实验管理成本的同时, 也提高了数据获取与分析的效率。已有研究表明, 物联网和信息技术的发展可以为生态环境监测和管理提供新的理论和方法, 使生态环境动态监测具备更好的实用性和便捷性^{[20, 22—23, 34]}。因此, 亟需强化物联网技术在热环境实测的应用工作, 进一步推进热环境研究向精细化的时空尺度深入发展。

4 结论

本研究利用物联网技术构建了在线式校园热环境监测设备, 获取了2022年10月16日至24日校园温度、风速、太阳辐射、地面温度四种参数, 探讨了小尺度下城市热环境的时空异质性特征, 主要结论如下：1)小尺度上校园热环境状况具有高度的时空异质性, 这种异质性主要源于周边建筑、植被等微尺度景观特征的影响。小尺度下大部分测点的太阳辐射和风速特征具有一定的同步性, 位于楼顶的01测点能够代表研究区的热环境总体特征, 表现为太阳辐射和风速与其它测点的变化规律十分相似。

2)不同测点之间的风速具有较好的相关性, 其中测点07和测点08的位置相隔仅有20m, 它们的相关性系数为0.88, 而测点10与06、07、08和09测点的相关性都比较低, 个别测点受局部景观特征的影响会形成相对独立的小气候。

3)各测点日间和夜间的风速和热岛强度都表现出明显的时空差异, 07和09测点的日、夜间风速变化幅度较大, 夜间最大风速达5m/s, 其它测点风速变化较小。不同测点的热岛强度也十分明显, 日间热岛动态变化都较大, 而夜间变化较小, 各测点之间的热环境特征差异与所在位置的植被、建筑特征直接相关。

4)地面温度与气温的关系较为复杂, 这种关系受风速的影响不大, 周边植被和建筑遮挡对二者的关系影响会更大。二者相关性在夜间会更为密切, 这是由于夜间没有太阳辐射的直接加热, 使得地面温度变化更为平缓, 直接影响了近地面气温的状况。