亚热带湿润区农业用地热环境效应及其生物物理因子分析

刘显男， 张良侠， 闫章美

(南京信息工程大学应用气象学院， 南京 210044)

土地利用对局地、区域甚至全球气候有着重要影响，成为除化石燃料燃烧外，人类影响气候变化的最大因素[1-5]。农业作为最主要的人为土地利用活动之一，通过调节陆地大气间能量、水分和动量交换在气候系统中发挥着至关重要的作用，特别是在热带和亚热带湿润区[6]。一方面，自然植被转化为农业用地会改变地表粗糙度、反照率、叶片电导和其他特性，从而增加或降低温度；另一方面，农田管理(如施肥、灌溉、松土和种植等)能显著影响当地气候[7-8]。

然而，目前关于农业用地热环境效应的认识大多基于模型模拟的方法[9-10]。由于模型的结构、参数化以及输入数据的差异，导致不同研究结果分歧较大，甚至完全相反[9-12]。同时，数值模型的分辨率通常较低(约100 km)，不能清晰反映温度变化的空间异质性。快速发展的卫星遥感技术，为评估农业用地热环境效应提供了最直接有效的观测证据。然而，已有研究大多基于代表性网格内不同土地利用类型的温度差异来量化毁林或造林的气候效应[10,13-14]，不能全面反映农业用地对区域温度的影响，特别是大面积农业用地集中分布区域。在我国，农业用地面积占比超过50%，其通过各种因素对气候造成的增温影响不可忽视，但目前尚缺乏单独对农业热环境效应长期变化趋势的遥感监测研究。因此，急需发展新的遥感评估方法，全面评估农业用地对区域温度的影响。

本文以江西省为例，基于MODIS地表温度数据，以天然林地为参照，采用空间代替时间的方法，从像元尺度定量评估近20年农业用地的热环境效应及其时空演变规律，并结合地表蒸散、反照率和气象数据，探讨农业用地热环境效应的可能驱动因子。旨在揭示我国亚热带湿润地区农业土地利用对温度的影响特征及其控制因子，为制定缓解区域气候变化的土地利用策略提供参考。

1 数据与方法

1.1 研究区域

江西省位于中国东南部(图1)，长江中下游南岸，属亚热带温暖湿润季风气候，年降雨量为1 341～1 943 mm，年平均气温为16.3～17.5 ℃。江西省是一个农业人口多、农村地域大、农业比重相对较高的省份，全省耕地面积约3 万km2；粮、猪、油、菜、水产等主要农产品产量在全国占有重要地位，是长三角、珠三角和港澳等地重要的农产品供应基地，其中稻谷产量居全国第2位，始终保持全国农业主产省的地位。

图1 江西省地形图Figure 1 Topographic map in Jiangxi Province注:此图基于国家测绘地理信息局标准地图(审图号:GS(2020)4619号)绘制，底图无修改。

1.2 数据来源与预处理

1.2.1 数据来源 本文利用中分辨率成像光谱仪(Moderate Resolution Imaging Spectroradiometer，MODIS)合成的地表温度数据(Land Surface Temperature，LST)来反映研究区的热环境条件。该遥感数据误差在1 ℃之内，已被广泛用于评估人类土地利用活动的热环境效应[15-18]。文中选取了覆盖研究区范围内的2003—2019年最新版本的MODIS/Aqua 8天合成的LST数据产品(MYD16A2，V6)，空间分辨率为1 km，通过Earthdata数据共享平台(https:∥search.earthdata.nasa.gov/)免费获取。Aqua卫星过境时间为地方太阳时13:30和01:30，此时地表温度分别接近日最高值和日最低值，分别用来代表白天、夜晚温度。其中蒸散数据来源于MODIS 8天合成产品(MYD16A2，V6)，空间分辨率为500 m。反照率数据来源于MODIS逐日数据产品(MCD43A1，V6)，空间分辨率为1 km。

土地利用数据来源于清华大学宫鹏教授领衔研发的第二代FROM-GLC土地覆盖产品,该数据参考IGBP(InternationaI Geosphere Biosphere Programme)分类体系，空间分辨率为30 m[15]。

气象数据来自江西省气象局的26个气象站(图1)2003—2019年的气温和降雨量资料，海拔数据来源于ASTER GDEM。

1.2.2 数据预处理 利用MRT(MODIS Reprojection Tool)软件对遥感数据进行批量拼接和投影处理。为了减少误差，本文仅提取平均发射率误差≤0.02(质量控制(QC)标志为“01”)、平均LST误差≤1 K(QC标志为“00”)和天顶角≤65°的LST像元；通过ANUSPLIN软件对气象数据进行插值。

1.3 方法

1.3.1 土地利用分类 本文获取2015年的土地利用类型数据，经拼接、投影、裁剪后得到江西土地利用类型图(图2)，主要包括农田、林地、草地、灌丛、湿地、水体、不透水面和裸地8种土地利用类型。其中林地和农田是土地利用的两大主要类型，分别占研究区总面积的63%和24%。

图2 江西省土地利用图Figure 2 Land use map in Jiangxi province

由于水体白天的冷却效应和夜间增暖效应会影响周围温度，本文首先剔除了湿地、水体及大型水体周围1 km的缓冲区域内的土地利用数据。为了使土地利用数据与MODIS LST数据的空间分辨率一致，本文分别计算了空间分辨率为1 km的农田和林地的面积占比，并把农田占比超过50%的像元定义为农业用地。

1.3.2 农业用地热环境效应估算方法 本文将农业用地的热环境效应(ΔT)定义为：

ΔT=Tagr-Tfor，

(1)

其中，Tagr为农业用地像元的实际地表温度，可通过MODIS地表温度数据直接提取；Tfor为天然林地的地表温度。ΔT>0表示农业用地对区域环境为增温效应，ΔT<0表示农业用地对区域环境为降温效应。由于缺乏对应每个农业用地像元位置的Tfor观测数据，本文基于农田周围天然林地的LST数据，以海拔为协变量，利用协同克里金插值方法获取Tfor。为了准确估算Tfor，降低混合用地对Tfor的影响，本文仅筛选林地面积百分比为100%的纯天然林地像元为参照。为了增加参考林地像元分布范围，同时选取了研究区周围50 km范围内的天然林地像元。

1.3.3 农业用地热环境效应时空变化及其驱动因子分析 本文首先计算2003—2019年的多年平均ΔT，分析农业用地热环境效应空间分布及其昼夜与季节变化特征。进而基于2003—2019年的逐年ΔT数据，分析农业用地热环境效应的年际变化趋势。为了探究ΔT变化与主要地表生物物理因子、气候因子的关系，本文首先采用与估算ΔT相同的方法，分析了农业用地对蒸散(evapotranspiration,ET)和反照率的影响。反照率数据包括白空反照率(White Sky Albedo，WSA)和黑空反照率(Black Sky Albedo，BSA)，由于二者的空间分布高度类似，本文采用WSA表征反照率的大小。然后，分析了ΔT的时空变化与ΔET、ΔWSA、气候因子(气温和降雨量)的相关关系。

2 结果与讨论

2.1 结果

2.1.1 农业用地热环境效应的时空变化特征 由热环境效应的空间分布(图3)可知：(1)在白天，农业用地的整体表现为增温效应，增温幅度最高达8 ℃，且不同季节差异较小，在夏季的增温效应略低于在其他季节的(图3A至图3E)。(2)在夜晚，农业用地降温效应明显，最高降温幅度达-6 ℃，且不同季节差异明显(图3F至图3J)：在冬季的降温幅度最大，秋季次之，在春季的增温和降温效应并存，绝大部分区域在夏季表现为增温效应。(3)南昌、上饶市与九江市的交界地带在夏季呈现微弱的降温效应(图3B)，这可能与周围大型水体(鄱阳湖)的冷岛效应有关。

图3 江西省农业用地热环境效应空间分布Figure 3 Spatial distribution of the thermal environment effects of agricultural lands in Jiangxi Province

从地区分布(表1)可知：(1)在白天，赣州市的增温效应最明显，增幅超过3 ℃；上饶市的最低，增幅低于1.5 ℃。就整个区域平均来看，农业用地导致白天年均增温1.97±1.02 ℃。(2)在夜晚，新余市的降温效果最明显，年均降温-0.51±0.44 ℃；景德镇市次之；赣州市的降温幅度最低，仅-0.07±0.41 ℃。整体上看，农业用地导致夜晚年均温度降低0.31 ℃。

表1 江西省各地级市农业用地平均热环境效应Table 1 The mean thermal environment effects of agricultural lands in all the prefecture-level cities of Jiangxi Province ℃

从逐月变化趋势(图4A)可知：在6月份的增温效应最低(1.3 ℃)。从年际变化趋势(图4B)可知：2003—2019年，江西省农业用地的白天增温效应和夜晚降温效应整体上均呈加剧态势，在白天每10年平均上升0.105 ℃，在夜晚每10年平均下降0.084 ℃。

从空间分布(图5)可知：(1)在白天，约67%的农业用地像元的增温效应呈增加趋势，主要集中在靠近城市和水体的地方。(2)在夜晚，除靠近城市和水体少部分农业用地像元外，约72%的农业用地像元呈降温效应；大部分农业用地的ΔT年际变化趋势的显著性水平高于0.05，这可能与本研究时段较短有关。

图5 农业用地热环境效应年际变化趋势空间分布图Figure 5 Spatial distributions of the inter-annual trends of the thermal environment effects of agricultural lands

2.1.2 农业用地热环境效应与驱动因子的关系 由农业用地对地表蒸散和反照率的影响(图6，图7)可知：(1)农业用地导致绝大部分地区蒸散降低，年均降幅达160 mm，且不同季节差异较大。夏季时的蒸散降低幅度最大，特别是在7月份的降低幅度接近40 mm；冬季时的整体变化幅度较小，大部分农业用地(54%)的蒸散呈现略高于天然林地的态势。(2)农业用地导致大部分地区(87%)的反照率增加，平均增加幅度为0.02，且不同季节影响差异极大。在秋季和冬季的增幅明显高于在春季和夏季的，其中在2月份的增加幅度最大(接近0.03)，在5月份的最低(平均变化几乎接近0)。(3)从年际变化趋势来看，在2003—2019年间，ΔET略有增加趋势，但并不显著；ΔWSA则呈显著下降态势(图7B)。

图6 ΔET和ΔWSA的空间分布Figure 6 Spatial distributions of ΔET and ΔWSA

图7 ΔET与ΔWSA的年内和年际变化趋势Figure 7 Intera-annual and inter-annual trends of ΔET and ΔWSA

农业用地热环境效应的时空变化与地表生物物理因子(蒸散和反照率)以及气候背景变化密切相关。由图8可知：白天的ΔT与ΔET呈显著负相关，与ΔWSA呈显著正相关；夜晚的ΔT与ΔET、ΔWSA均呈显著负相关；白天的ΔT与气温呈明显的正相关，与降雨量的相关关系相对较弱。

图8 年均ΔT与驱动因子的相关关系Figure 8 Correlations between annual ΔT and driving factors

2.2 讨论

针对农业用地生物物理温度效应研究薄弱问题，本文以地处亚热带湿润区的江西省为例，基于卫星观测技术，全面评估了农业用地的热环境效应与驱动因子，结果发现农业用地在白天具有明显的增温效应，而在夜晚的降温作用明显(图3)，这与已有的以国家和全球为范围进行研究得到的结果[10,13]一致。此外，本文发现农业用地的热环境效应随季节变化明显，白天的增温效应在夏季较弱，可能与大面积水稻种植过程中的农田灌溉活动有关。大量研究证实灌溉可导致蒸散增加，进而对环境有明显降温效果[7-8]。夜晚的降温效应在夏季变为微弱的增温效应，这可能与夏季蒸散加强和反照率差异明显降低有关(图7A)。从长期变化趋势来看，本文的研究结果表明农业用地在白天的增温效应逐渐加强(图7B)，但日益加剧的夜晚降温效应在一定程度上可抵消白天增温效应的影响。全球气候变暖、农业活动加强以及快速城市化是导致农业用地白天增温效应加强的可能原因。例如，随着温度的升高，可加剧农业用地与天然林地间的ET差异(图7B和图9B)；快速城市化可通过城市热岛效应影响周围农田的气候效应[19]。

与已有研究[10,13,20]类似，本文发现蒸散是控制农业用地白天的热环境效应的主要因子(图8)。与天然林地相比，农业用地的蒸散降低，进而导致潜热减少，从而增加地表温度。但意外的是，本文发现农业用地在白天的增温效应与反照率呈显著正相关关系(图8)，这与普遍认为的“反照率增加将降低地表净辐射，进而降低温度”[21]相矛盾。导致这种现象的可能原因是亚热带湿润区农田蒸散的作用远远超过反照率的影响[10],掩盖了反照率的可能降温效应。另外，相对于森林而言，农田在夜晚的反照率较高，间接起到了降温效果[22]。 对比之下，本文发现农业用地在夜晚的热环境效应受蒸散和反照率的共同控制。与已有研究结果[23-24]不同的是，本文发现气候因子对农业用地的热环境效应空间分布的影响相对较小(图8)，这可能与研究尺度有关。本研究区内气候背景差异相对较小，其影响被局地生物物理因子所掩盖。但从年内、年际变化趋势(图9)来看，气候因子的作用和生物物理因子作用相当，甚至超过生物物理因子的影响，特别是气温的影响，强调了气候变化在预估未来农业用地气候效应中的重要作用。

图9 ΔT的逐月和年际变化与各驱动因子的相关系数Figure 9 Correlation coefficients between month-to-month and inter annual changes of ΔT and driving factors注：为33 338个像元r值的平均值。

3 结论

本文以江西省为例，分析了亚热带湿润地区农业用地热环境效应的时空变化特征及其与生物物理和气候因子的关系，主要结论如下：

(1)农业用地在白天以增温效应为主，在夜晚主要表现为降温效应，且季节变化大；

(2)2003—2019年，大部分农业用地的白天增温效应和夜晚降温效应均呈加剧态势；

(3)农业用地热环境效应的空间分布在白天主要受地表蒸散的影响，在夜晚受蒸散和反照率共同控制，而且年际变化趋势与气候因子(特别是气温)密切相关。

本文的科学和实践意义如下：首先，文中使用的农业用地热环境效应的定量评估方法可有效避免天然林地地理分布不均的局限性，为区域尺度土地利用活动的气候效应评估提供了一种新的视角。其次，虽然农业用地在白天具有明显增温效应，但在夜晚表现为降温效应，这与文献[25-26]的结论一致。因此，在未来评估土地利用活动对区域气候的影响时，应加强对农业用地的研究。最后，本文发现气候变化可能加剧农业用地的热环境效应，说明在未来区域气候变化评估中，需要综合考虑农业用地气候效应与气候变化的互馈机制。