成渝地区双城经济圈，位于“一带一路”和长江经济带交汇处，以成都和重庆为中心，包含四川15个市和重庆29个区县（见图1），是西南地区经济水平最高，发展潜力最大的城镇化区域^[22]。地形与四川盆地基本一致，中部多为平原，周围由低山、丘陵等环绕，气候适宜，自然资源丰富，生态禀赋优良^[23]。

图  1  研究区地理位置示意图

Figure 1.  Schematic diagram of geographical location of the study area

植被NPP数据来源于美国NASA EOS/MODIS的2002—2021年（实际计算植被NPP年份为2001—2020年）的MOD17A3HGF6.0版本的全球NPP数据，空间分辨率为500 m。基于GEE平台对MODIS数据进行预处理，获取统一坐标的成渝地区双城经济圈20年的植被NPP数据。对预处理后的各年份NPP数据乘以转换因子并剔除多年份NPP平均值为0的像元。

数字高程数据来源于地理空间数据云平台，分辨率为30 m。气象数据来自中国气象局数据信息中心的“中国地面气候资料日值数据集(V3.0)”，其数据时间尺度为2001年1月—2020年12月，包括日照时数、风速、降水、湿度、水面蒸发、日平均温度等6个气象要素逐日数据。基于成渝地区双城经济圈内的气象站点日观测数据，采用ANUSPLIN软件，以DME为协变量插值生成空间分辨率为1 km的气象栅格数据，再经过重采样生成所用的气象数据，空间分辨率为500 m。对比其他插值方法，ANUSPLIN插值精确度得到显著提高^[15]。

采用一元线性回归模型来分析2001—2020年成渝地区双城经济圈植被NPP整体变化趋势及空间异质性。

式中，$ \mathrm{\theta } $是NPP年际变化趋势的斜率，n为统计时段的累计年数，$ {NPP}_{i} $为第i年植被npp值，若$ \mathrm{\theta } $>0则表明植被NPP呈上升趋势，若$ \mathrm{\theta } $<0则表明植被NPP呈下降趋势，若$ \mathrm{\theta } $=0则表明植被NPP几乎无变化。

采用F检验法对植被NPP的变化趋势进行显著性分析，其公式如下：

式中，U表示为NPP的方差：$ U={\sum }_{i=1}^{n}{(N\widehat{P}{P}_{i}-\overline{NPP})}^{2} $，Q表示为NPP平方和：$ Q={\sum }_{i=1}^{n}{({NPP}_{i}-N\widehat{P}{P}_{i})}^{2} $，$ N\widehat{P}{P}_{i} $是$ {NPP}_{i} $的回归值，$ \overline{NPP} $是20年内NPP的平均值。

变异系数($ {C}_{v} $)表示为标准差和平均值的比值，是用来描述数据离散程度的指标^[24]，采用$ {C}_{v} $来反映近20年植被NPP的波动特征，其公式如下：

式中，$ {NPP}_{i} $为第i年植被npp值，$ \overline{NPP} $是20年内NPP的平均值；$ {C}_{v} $值越小，植被NPP变化越小，反之亦然^[25]。参考ArcGIS软件中的自然分类法^[26]，将$ {C}_{v} $分为以下五个等级：稳定（$ {C}_{v} $<=0.10）,较稳定（0.10<$ {C}_{v} $<=0.20），不稳定（0.20 <$ {C}_{v} $<=0.30），极不稳定（$ {C}_{v} > 0.30 $）。

基于像元偏相关分析方法计算植被NPP与气温、降水影响因子的偏相关系数，其公式如下：

式中，$ {R}_{xy,z} $表示为控制变量z之后x与y的偏相关系数，x表示NPP值，y表示降水量，z表示气温，$ {R}_{xy} $、$ {R}_{xz} $、$ {R}_{yz} $分别表示NPP与降水量、NPP与气温以及气温与降水量的线性相关系数。采用t检验法对偏相关系数进行显著性检验。

复相关分析法是研究几个要素与某一要素的相关关系^[27]。采用复相关分析法探究气温与降水对植被NPP的共同影响作用，其公式如下：

式中，$ {R}_{xy} $表示NPP与降水量的线性相关系数，$ {R}_{xz,y} $表示NPP与气温的偏相关系数。采用F检验法对复相关系数进行显著性检验。

根据陈云浩等学者研究的植被变化驱动因子分区准则，将引起成渝地区双城经济圈NPP变化的驱动力进行分区，根据偏相关系数与复相关系数的显著性检验结果，将驱动力分为下表五类（见表1）^[27]。当复相关显著时：气温与NPP偏相关的像元定义为气温驱动；降水与NPP偏相关的像元定义为降水驱动，气温、降水分别与NPP偏相关的像元定义为气温降水强驱动；气温、降水与NPP偏相关不显著的像元定义为气温降水共同驱动。不满足复相关显著的像元定义为非气候因子驱动。

研究区内2001—2020年植被NPP均值分布如图2所示，成渝地区双城经济圈植被NPP平均值为563.14 g C·m⁻²·a⁻¹，年均值分布在<400 g C·m⁻²·a⁻¹、400—500 g C·m⁻²·a⁻¹、500—600 g C·m⁻²·a⁻¹、600—700 g C·m⁻²·a⁻¹、700—800 g C·m⁻²·a⁻¹、>800 g C·m⁻²·a⁻¹的区域分别占总面积的5.49%、28.79%、34.31%、14.33%、8.95%、8.13%。从空间分布来看，研究区内植被NPP呈现四周高中间低的分布格局，其与四川盆地地形特征有着密切联系。结合不同海拔下植被NPP均值分布图来看，随着海拔的升高，植被NPP均值也在增加，在海拔>3000 m的情况下，均值开始显著降低，但通过整体趋势线可以看出植被NPP呈增加趋势。在德阳、成都和雅安一带海拔较高的山脉植被NPP均值高，在海拔较低的盆地内部的各城市中心植被NPP均值最低，在重庆黔江区、彭水县等低山丘陵处植被NPP均值较高。

图  2  成渝地区双城经济圈2001—2020年植被NPP均值分布图及不同海拔下植被NPP均值分布图

Figure 2.  Mean NPP distribution of vegetation in Chengdu-Chongqing economic circlein 2001-2020 and mean NPP distribution of vegetation at different altitudes

利用一元线性回归法及变异系数分析研究区植被NPP的时空变化特征（见图3a和3b），并对变化趋势做显著性检验，将结果分为：显著减少、减少、稳定、增加、显著增加五类（见表2和图3c）。

图  3  成渝地区双城经济圈2001—2020植被NPP变化趋势

Figure 3.  Variation trend of vegetation NPP in Chengdu-Chongqing economic circle from 2001 to 2020

2001—2020年植被NPP变异系数如图3a所示，可以看出成渝地区双城经济圈近20年来植被NPP较为稳定，大部分区域变异较弱，变异较强的地方存在于成都市和重庆市主城区周边以及沿长江一带。从斜率分布图(见图3b)及近20年植被NPP变化分类图（见图3c）来看，其变异较强的区域都呈现显著减少趋势；在植被NPP均值较高的德阳、成都和雅安一带的山脉以及重庆市的黔江区、彭水县等低山丘陵处变化甚微，趋于稳定；在城市群内部盆地及丘陵地带植被NPP呈现增加趋势；结合表2统计结果来看，其植被NPP呈现显著增加趋势的面积最大，占研究区总面积的70.30%；呈现减少趋势的面积最小，其面积占比为0.43%。从2001—2020年植被NPP年均值变化（见图3b和3d）来看，整体呈现增加趋势，其中植被NPP呈上升趋势的面积占比达95.08%，呈下降趋势的面积占比为4.92%；2005年研究区植被NPP年均值最小为471.93 g C·m⁻²·a⁻¹，植被NPP年均值最大出现在2020年为644.99 g C·m⁻²·a⁻¹。

将成渝地区双城经济圈海拔划分成不同梯度等级，统计各海拔梯度下植被NPP的变化情况，从图4可以看出，在研究区内，随着海拔的增加，植被NPP的年变化趋势率不断下降，在海拔小于500 m的区域植被NPP增长趋势最显著，年变化趋势率达到8.72 g C·m⁻²·a⁻¹；在海拔大于3000 m的区域植被NPP最稳定，年变化趋势率为0.95 g C·m⁻²·a⁻¹；在海拔500—1000 m、1000—2000 m、2000—2500 m及2500—3000 m高程区间内，年变化趋势率分别为7.24 g C·m⁻²·a⁻¹、4.62 g C·m⁻²·a⁻¹、2.99 g C·m⁻²·a⁻¹、1.94 g C·m⁻²·a⁻¹。从整体来看，成渝地区双城经济圈植被NPP整体处于增长趋势，在海拔大于3000 m的区域人类活动少、植被覆盖度高其植被NPP均值较高且变化小，在海拔小于500 m的区域，由于城市扩张、人类活动干扰其植被NPP均值相对较低，近年来天然林保护、长江流域生态屏障构建等措施致使植被NPP均值得到显著增加。

图  4  成渝地区双城经济圈各海拔梯度植被NPP年际变化图

Figure 4.  Interannual variation of vegetation NPP at different elevations in Chengdu-Chongqing economic circle

通过一元线性趋势分析法，可知研究区20年里气温和降水都呈现增长趋势，年变化趋势率分别为0.0067℃a⁻¹、8.14 mma⁻¹。

研究发现成渝地区双城经济圈植被NPP与气温、降水的偏相关系数范围都在−0.83~0.87之间，且相关性空间格局存在异质性（见图5）。分析表明研究区植被NPP与气温呈显著负相关、负相关的区域面积占比分别为0.39%、25.96%，呈显著正相关、正相关的区域面积占比分别为6.22%、67.43%；与降水呈显著负相关、负相关的区域面积占比分别为0.21%、7.33%，呈显著正相关、正相关的区域面积占比分别为23.40%、69.06%。

图  5  成渝地区双城经济圈植被NPP与气温、降水的偏相关系数（a、b）及显著性（c、d）

Figure 5.  Partial correlation coefficients (a, b) and significance (c, d) of vegetation NPP with air temperatureand precipitation in Chengdu-Chongqing economic circle

从空间上来看，研究区植被NPP与气温正负关系并存，在其西南部即眉山市、乐山市和雅安市等地植被NPP与气温呈负相关关系，在重庆市黔江区及彭水县部分地区存在负相关关系且零星出现显著负相关关系；在研究区中部及东北部植被NPP与气温正相关关系明显，尤其是在重庆城市群及东部城市群的部分地区存在显著正相关关系。研究区大部分区域植被NPP与降水之间都呈现正相关关系，并且在成都城市群与东部城市群的区域正相关关系显著；只有在东部城市群零星存在负相关关系。其结果表明了气温在不同区域对植被NPP变化的影响不同，而降水在研究区内几乎都起着积极影响。

植被NPP的变化除了受气温、降水等单方面因子的作用，同时受气温与降水等因子的共同驱动作用，通过植被NPP与气候因子的相关性及其显著性分析，将驱动因子进行分区，从图6可知，植被NPP与气温、降水驱动因子的复相关系数范围为0.002~0.867，其高值区主要分布在东部城市群及成都市与德阳市的西部，在东部城市群的东北部、雅安市与宜宾市的屏山区系数值较低。研究发现成渝地区双城经济圈植被NPP变化主要受非气候因子驱动，其面积占比达86.87%；受气候因子影响的区域面积占研究区总面积的13.13%，其中受气温驱动、降水驱动、气温降水共同驱动、气温降水强驱动的面积占比分别为1.58%、8.40%、0.16%、2.99%。从空间上来看，受气温驱动的区域主要分布在东部城市群的平行岭谷地带及重庆城市群的江津区、綦江区；受降水驱动的区域分布最广，包括了成都城市群的西北部及南部部分区域、南部城市群的泸州市部分区域、重庆市黔江区部分区域；受气温、降水共同驱动的区域十分有限，分布稀少；在南充仪陇、蓬安县一带受气温、降水强驱动的区域最为密集。

图  6  复相关系数及驱动力分区图

Figure 6.  Partition diagram of complex correlation coefficient and driving forces

基于MOD17A产品NPP数据，利用一元线性趋势法分析2001—2020年成渝地区双城经济圈植被NPP变化特征，并结合气象数据，运用相关分析法，研究植被NPP变化的气候因子驱动力。

研究结果表明成渝地区双城经济圈植被近20年间NPP在471.93—644.99 g C·m⁻²·a⁻¹范围之间，均值为563.14g C·m⁻²·a⁻¹，比全国平均值514.48 g C·m⁻²·a⁻¹高48.66g C·m⁻²·a⁻¹，其总体呈现增长趋势，与冯婉^[12]、张凤英^[28]等学者的关于成渝城市群地区的研究结果一致。从年际变化上来看，研究区内植被NPP整体呈现缓慢增长的趋势，增长率为7.53 g C·m⁻²·a⁻¹。在植被NPP均值较高的德阳、成都和雅安一带的山脉以及重庆市的黔江区、彭水县等低山丘陵处趋于稳定；在城市群内部盆地及丘陵地带植被NPP呈现增加趋势，这一现象得力于“退耕还林、还草政策”、“长防林工程”、“长江保护修复攻坚战行动计划”^[29]。从空间上来看，研究区内植被NPP呈现四周高中间低的分布格局，随着海拔的升高，植被NPP值不断增加。海拔在2500—3000 m的高程区间内植被NPP均值为805.21 g C·m⁻²·a⁻¹，海拔小于500 m的高程区间内植被NPP均值为491.06 g C·m⁻²·a⁻¹。这与四川盆地地形息息相关，这可能是因为研究区靠近川西高原的部分高海拔区域，由于气候影响，生态环境脆弱，植被覆盖率低，NPP低且较为稳定，在盆地城市内部由于地势平坦，水热条件好，其城市不断扩张，NPP值较低，在低山丘陵地区，因气候适宜而少受人类活动影响，植被覆盖广，NPP值最高，该结果与胥晓^[30]、张锐^[31]等学者的研究结论相符合。

此外，基于气象站点数据，利用ANUSPLIN软件对气象数据进行插值，分析气候因子对植被NPP变化的影响。结果表明研究区植被NPP与气温、降水的相关性存在空间异质性，气温在眉山市、乐山市和雅安市等地以及重庆市黔江区、彭水县部分地区对植被NPP起着抑制作用，在重庆城市群及东部城市群的部分地区促进了植被NPP的增长，而降水与植被NPP在整个研究区内都呈正相关的关系，这可能是因为水热条件决定了植被的生长，而降水的升高使植被获得更多的水分，促进植被生长，使得NPP不断增加，但研究表明气温对植被NPP的影响是个复杂的过程，气温升高能够为植被提供更多的热量促进植被光合作用从而导致NPP增加，但过高的温带增长会对植被产生高温胁迫^[10,32]。此外，研究表明研究区内受到非气候因子影响的区域面积占比达到86.87%，植被NPP分异较大程度受到了人为因素干扰，但仅针对气候因子来说，气温、降水仍然是影响植被NPP变化的主要因素。

以往大多数研究多选择长江流域、川西高原或横断山区这样的生态区内^[4,12,33]，但对成渝这样的城市群体研究甚少，成渝城市群作为长江上游最大的城市群体，其植被NPP变化趋势对长江上游生态屏障构建及生态修复有着举足轻重的影响。