Response of Vegetation to Extreme Temperature in the Three Gorges Reservoir Area

三峡库区位于长江上游地段，位于105°50′E—111°40′E、28°31′N—31°44′N之间，西起重庆市江津区，东至湖北省夷陵区，面积约5.8万km²。库区属于我国第二级阶梯的东部，地势起伏较大，地貌复杂多样，有山地、丘陵、台地、平原等地形（见图1）。其气候为亚热带季风性湿润气候，该气候季节变化特征明显，年平均气温在16~18℃之间，年平均降水量集中在1000~1300 mm，雨热同期，集中在夏秋两季。

Figure 1.  Location and distribution map of vegetation types in the study area

气象数据来自中国气象局国家气象信息中心（http://data.cma.cn/），共选取22个气象站点，并选取了2000—2019年各个气象站点的日最高气温、日最低气温、日平均气温以及日降水量等数据，并控制数据质量，具体做法为：检验日降水量为异值的数据，一天中最低气温比最高气温大的值，以及对个别缺失的数据采用相邻日期数据进行填补，以EXCEL进行统计。

NDVI数据源于美国国家航天局（NASA，https://adsweb.nascom.nasa.gov/）所提供的2000—2019年中分辨率成像光谱仪（MODIS）光谱反射率数据产品，产品名为MOD13A3月合成植被指数，其空间分辨率为1 km，时间分辨率为30 d，此数据集已经采用几何校正、辐射校正以及大气校正等方法进行了初步的处理。此数据集缺失了MODIS NDVI 2000年1月份的数据，因此使用2001年1月份的数据进行了填补。

采取国家地球系统科学数据中心（http://www.geodata.cn/）发布的2019年1km植被类型空间分布数据，主要采用的植被种类有针叶林、阔叶林、灌丛、草丛。

RClimDex模型是一个专门的极端气候计算模型，它是由加拿大气象局气候研究部所研发的，这种模型在国际上有一定的影响力，通过输入编好格式的逐日的降水量、逐日的最低温以及逐日的最高温数据，便可得出代表极端降水的11个指数以及代表极端气温的16个指数。其中代表昼夜极端温度变化的指数见表1^[20-22]。

趋势分析是指时间序列的气候要素在长时期内呈现出来的持续增加或者减少的变化分析。本研究拟采用一次线性方程表示极端气候的变化趋势^[23~25]，即：

式中，θ_slope表示的是线性趋势，其单位是℃·10a⁻¹、mm·10a⁻¹等；n研究年份；T_i表示某一年份的NDVI或极端气候指数。θ_slope>0表示研究对象随时间上升，θ_slope<0表示研究对象随时间减少，θ_slope的值越大代表研究对象的变化幅度就越大。

在MATLAB平台中利用 Pearson 相关系数来度量极端气候事件与植被覆盖之间的相关性。该方法是英国统计学家 Pearson提出的一种计算直线相关的方法。若随机变量X，Y的联合分布是二维正态分布$\left( {{x_1},{y_1}} \right),\left( {{x_2},{y_2}} \right), \cdots ,\left( {{x_n},{y_n}} \right) $为n次独立观测值，r代表 Pearson 相关系数^[26-27]，其计算公式为：

其中，相关系数的取值介于[−1,1]之间，|r|越大则代表研究变量间的相关程度越高。

据图2，三峡库区白天与夜间的极端气温变化呈不对称性，冷夜日数（TN10p）、暖夜日数（TN90p）、冷昼日数（TX10p）和暖昼日数（TX90p）的年际变化率分别为：−6.59 d·10a⁻¹、 6.53 d·10a⁻¹、−2.67 d·10a⁻¹、5.54 d·10a⁻¹。因全球变化的原因，可以看出TN10p、TX10p两个冷指数呈现减少的趋势，TN90p、TX90p两个暖指数呈现增加的趋势；而以昼夜来看，TN10p、TN90p两个夜间温度指数的变化幅度要比TX10p、TX90p两个白天温度指数的变化幅度大。

Figure 2.  Variation trend of extreme temperature in the Three Gorges reservoir area from 2000 to 2019

NDVI与TX10p在绝大部分地区呈负相关性，同时重庆东北部和东南部的相关性较大，两者的相关性多在−0.6~−0.2区间（见图3），而通过显著性检验的地区主要的分布在渝东北的奉节、巫溪、云阳、开州与渝东南丰都、石柱南部；呈正相关负相关性的地区主要在重庆主城区及湖北巴东的南部区域，但仅有27个像元通过了显著性检验，并且其零散分布。

Figure 3.  Correlation coefficient between vegetation and temperature extremes.

NDVI和TX90p多呈现正相关的关系，多集中在0.2~0.6区间（见图3），主要分布在江津区、巫山县、巫溪县以及湖北段西部等区域，通过显著性检验的区域在渝东北巫山、巫溪、奉节南部，渝东南武隆、石柱北部以及江津；呈负相关的地区主要分布在重庆主城区及其周边的涪陵北部、丰都中部，湖北宜昌的西南部地区，通过显著性检验的地区主要在丰都中部与宜昌西南部。

NDVI和TN10p几乎均为负相关的关系，渝东北部区域和渝东南部区域的相关系数相对而言较高，大多在−0.4~−0.6区间（见图3），通过显著性检验的地区主要为渝东南的丰都石柱、渝东北巫山、巫溪、奉节、云阳等地；而重庆主城区的NDVI与TN10p多呈正相关，同时通过显著性检验的像元也多分布在重庆主城区，此外在涪陵北部、巴东南部也有较为零散的地区呈正相关。

NDVI和TN90p之间多呈现出正相关性，大约有一半的地区相关系数是0.4~0.8（见图3），其中有40%的区域通过了显著性检验，表明库区内近4成地区的NDVI与TN90p之间存在极为显著的相关性，且主要集中在重庆主城周边的江津、长寿和渝东北的忠县至湖北巴东一带；呈负相关的区域为重庆主城、涪陵北部、宜昌西南部等地区，通过了显著性检验的负相关系数较少，在重庆主城及宜昌西南部有较少显著负相关地区。

通过上述分析, 可见三峡库区植被对昼夜极端气温指数的变化具有不同的响应特征, 其中NDVI对夜间极端气温指数的响应更加显著。

通过对三峡库区不同植被类型的NDVI与冷昼日数（TX10p）和暖昼日数（TX90p）、冷夜日数（TN10p）、暖夜日数（TN90p）的相关性分析（见表2），发现不同类型的植被对昼夜极端气温变化的相关性存在着差异。TX10p、TN10p对针叶林和草丛的生长起的阻碍作用最大，TX90p、TN90p的增加对各植被类型的生长均产生了积极的作用，且TN90p的增多所产生的积极影响更多，尤其是对阔叶林和草丛。

三峡库区2000—2019年昼夜温度均呈上涨的趋势与全球气候变暖具有一致性，同时也引起TN10p、TX10p两个冷指数呈减少趋势，TN90p、TX90p两个暖指数呈现了增加的趋势，以昼夜为划分，TX10p、TX90p两个昼指数的变化趋势不如TN10p、TN90p两个夜指数的变化趋势明显，但是植被的光合作用通常发生在白天，夜间进行生长，因此晚上的气温变化对植被的影响作用更大，而这种昼夜指数的不对称增减对植被生长的影响程度必然会加深。

整个库区范围内，渝东北的巫山、巫溪、奉节、云阳等地及重庆主城区的核心对这四个指数的响应更加明显；渝东北地区因自然植被较多以及山地较多，植被生长对气候变化的响应较强烈，也符合植被与冷指数呈负相关、与暖指数呈正相关的一般规律，而重庆主城区的核心则因城市化高，这一区域内的植被多为栽培的观赏性植物，受人为影响的程度较高，因此重庆主城的部分地区植被与冷指数呈正相关，而与暖指数则呈负相关。而库区湖北段、重庆主城区周围的江津、涪陵、长寿等地对这四个指数的响应较弱。

以2000—2019年三峡库区内22个气象站点的逐日最高温、最低温、平均温等气象数据，归一化植被指数数据（NDVI）以及植被类型数据，采用RClimDex极端气候指数模型、趋势分析法以及Pearson相关分析法，分析了三峡库区白天与夜间极端气温的不对称变化及其对植被的动态影响，得出以下结论：

（1）三峡库区近20年来昼夜温度均为上升趋势，但夜间极端气温指数的变化幅度大于白天极端气温指数的变化幅度。

（2）三峡库区植被生长对昼夜极端气温的响应具有一定的差异性，植被对夜极端气温指数的响应比对昼极端气温指数的响应更加突出。

（3）三峡库区植被生长对昼夜极端气温的响应在空间上具有较大的差异，重庆主城区的核心以及渝东北对极端气温的响应更加明显，库区湖北段、重庆主城区周围的江津、涪陵、长寿等地对极端气温的响应较弱。

（4）三峡库区不同植被类型与昼夜极端气温具有不同的相关性，TX10p、TN10p两个冷指数对植被生长多起阻碍作用，尤其对针叶林和草丛的阻碍作用最大；TX90p、TN90p两个暖指数对各植被类型的生长均产生了积极的作用，且TN90p的增多所产生的积极影响更多，尤其是对阔叶林和草丛。

植被的动态变化会受到地形地貌、气候条件及人为因子的影响，本文主要探讨了气温和植被间的关系特征，而未将降水、地形等因素考虑进去；此外使用的归一化植被指数数据分辨率为1km，且植被类型的划分是根据2019年的数据，对结果的精度会形成一定的影响。所以在今后的研究中应对多个要素综合考虑，从而尽可能地减少误差。