1982—2020年安徽省净生态系统生产力时空格局变化及其成因

华朗钦,张方敏,*,翁升恒,2,卢燕宇

1 南京信息工程大学应用气象学院气象灾害预报预警与评估协同创新中心/江苏省农业气象重点实验室,南京 210044 2 福建省气象科学研究所,福州 350007 3 安徽省气象局气象科学研究所大气科学与卫星遥感安徽省重点实验室,合肥 230031

研究表明,目前大气CO2浓度较工业化之前升高了约50%[1]。不断增加的大气CO2浓度改变了全球碳循环模式[2],可能加剧了气候变暖[3]。陆地生态系统与大气之间的碳交换是全球碳循环的重要过程,是大气CO2重要的汇[4],也是预估未来气候变化的不确定性来源之一[5]。净生态系统生产力(NEP)为净初级生产力(NPP)减去异养呼吸(Rh),是定量描述陆地生态系统碳循环的重要指标,NEP>0表明陆地生态系统为碳汇,反之为碳源[6]。因此,准确评估区域NEP的变化特征及对气候环境的响应有助于明确区域陆地生态系统碳收支状况,对深入理解区域碳循环变化机制有重要指导意义。

陆地生态系统NEP受到气候变化、大气CO2浓度与植被结构等因素影响[7],表现出较强的年际变化[8],在局地区域尺度上,NEP变化规律及主导因子并不一致。在我国不同区域开展的NEP时空变化及归因研究结果具有明显的区域差异。比如,李洁等[9]结合CEVSA(Carbon Exchange between Vegetation, Soil and Atmosphere)模型研究认为东北地区NEP没有显著变化趋势并且主要影响因素是降水量,刘春雨[10]基于CASA(Carnegie-Ames-Stanford Approach)模型针对甘肃省的研究发现NEP呈下降趋势的主导因素亦为降水量,但是刘凤等[11]基于CASA模型发现青藏高原NEP总体呈上升趋势,而气温对NEP的抑制作用且影响范围大于降水量;庞瑞等[12]利用CEVSA模型研究发现西南高山地区的NEP的主要因素也是气温,且为负贡献,导致NEP有显著下降趋势。综上所述,我国不同区域NEP的变化规律不同,主导因子也不一致,且研究大多集中于气温与降水量对NEP的相关分析,对我国东部不同区域NEP变化及驱动因子有待深入研究。

安徽省是我国东部重要的农业基地和重点林区,省内植被覆盖度高,生态区位极为重要[13],是研究我国东部植被生态系统响应气候与环境变化的典型区域。目前,已有学者研究了安徽省生态系统NPP的时空变化[14],分析了安徽省寿县农田NEP的时间变化[15],然而,安徽省生态系统NEP时空格局及驱动因子尚不明确。在全球尺度上,Chen等[7]指出CO2和氮沉降分别对全球NEP变化贡献47.0%和1.1%,叶面积指数(LAI)和气候分别贡献12.4%和-28.6%。另外,Tian等[16]认为我国东南部出现氮饱和导致其驱动的碳吸收持平的现象。因此,本文基于BEPS模型模拟结果,分析安徽省NEP的时空格局变化,并在前人研究基础上,分析NEP对气候因子、植被因子和大气CO2浓度的时空敏感性,借助通径分析和贡献率方法研究安徽省NEP时空变化归因,以期进一步加强对安徽省区域碳循环影响机制的科学认知。

1 资料与方法

1.1 研究区概况

安徽省(29°23′—34°39′N,114°52′—119°38′E)地处我国长江中下游华东地区,长江三角洲腹地,属于暖温带半湿润区和亚热带湿润季风区。1982—2020年间,年平均气温为15.8℃,年总降水量为1252.4 mm。安徽省以平原、丘陵为主,生态环境良好。39年来,土地利用类型发生改变的面积占全省面积的2.57%。据2010年土地利用资料统计(图1),农田为安徽省主要土地利用类型,占57.46%,其次为森林,占22.74%,主要分布在安徽省南部,草地占5.94%,水域、建设用地及未利用土地等占13.86%。

图1 2010年的安徽省土地利用类型分布图Fig.1 Land use map of 2010 in Anhui Province

1.2 资料

本文使用的气象资料来源于中国国家气候中心(http://data.cma.cn)提供的中国地面气候资料日值数据集,包括1982—2020年全国2420个气象站的平均气温(Air Temperature,Ta)、最高气温、最低气温、相对湿度、日照时数、降水量(Pre)等数据,太阳总辐射(Rad)由日照时数计算得出[17],饱和水汽压差(VPD)作为综合反映大气温湿状态的气候参数,由Ta与相对湿度计算得出[18]。在完成异常值剔除和缺测值线性插补的基础上,利用克里金插值方法[19]得到空间分辨率为1 km的栅格数据。1982—2020年的LAI资料来源于中国科学院地理科学与资源研究所的1 km分辨率的 GLOBMAP LAI V3产品,时间分辨率为8 d[20],产品已经得到了很好的验证。土地利用类型数据来源于中国科学院地理科学与资源研究所资源环境科学与数据中心(http://www.resdc.cn),将地表类型分为耕地、森林、草地、水域、建设用地和未利用地6种类型,在本研究中不考虑水域、建设用地与未利用地,空间分辨率为1 km,并使用1980年、1990年、1995年、2000年、2010年、2015年的土地利用类型作为模型输入的土地利用类型数据。以上插值后的空间数据,借助ArcGIS工具进一步裁剪出安徽省的栅格数据用于本研究的计算分析。1982—2020年的大气CO2浓度资料来源于夏威夷Mauna Loa全球本底站观测数据(https://www.ncei.noaa.gov)。为方便表述,本文分别将研究时段1982—1990年、1991—2000年、2001—2010年与2011—2020年依次表述为1980s、1990s、2000s与2010s。

1.3 NEP计算

BEPS(Boreal Ecosystem Productivity Simulator)模型是在FOREST-BGC(Forest Biogeochemical Cycles)模型基础上发展的基于过程的生物地球化学模型[21]。BEPS应用两叶模型将冠层分为阳叶和阴叶[22],并利用空间与时间的升尺度方案将瞬时的Farquhar光合模型转换为日尺度上的冠层光合模拟,进一步融合CENTURY模型的土壤生物化学模块,实现对陆地生态系统NEP的模拟[23]。BEPS模型在中国区域的碳循环研究中已有广泛的应用基础[24—26],具体计算过程详见Chen等[27]和Ju等[28]。本研究采用BEPS模型计算1982—2020年的1 km NEP,并借助ArcGIS工具裁剪出安徽省的1 km NEP栅格数据用于本研究的结果分析。

1.4 统计方法

1.4.1重心模型

通过分析要素的重心变化,可以突显出区域内要素的变化趋势与空间特征[29],公式如下:

(1)

(2)

式中,Gx、Gy分别为某一栅格数据重心的经向、纬向坐标;n、m为栅格数据的行、列数;Latij、Lonij分别为第i行、第j列格点的中心纬度、经度;NEPij为栅格数据第i行、第j列格点上NEP的值。

1.4.2通径分析

1.4.3敏感性系数与相对贡献率

敏感性系数可以较好地反映因子对某个要素变化的响应程度[32]。本文选取敏感性系数作为衡量各个因子对NEP影响程度的指标,计算方法如下:

(3)

式中,fi代表第i个因子;Sfi为NEP对该因子fi的敏感性系数;Pfi为该因子对NEP的偏导数,Norfi为该因子的归一化参数。敏感性系数的绝对值越大,表明NEP对该因子的变化越敏感。敏感性系数为正,表明该因子的变化会导致NEP往相同趋势变化;敏感性系数为负,表明该因子的变化会导致NEP往相反趋势变化。

各个因子的相对贡献率由该因子的相对变化率(Rfi)与敏感性系数(Sfi)计算得出[33]:

(4)

(5)

式中,Cfi为该因子对NEP的相对贡献率;Rfi为该因子在研究时段内的相对变化率;n为因子个数,本文取6;R2为决定系数;year为研究时长,本文取39;CRfi为该因子多年变化率;Meanfi为该因子研究时段内的平均值。

将因子对NEP的偏导数与某年因子距平相乘,即可计算出该因子贡献的NEP理论变化:

(6)

式中,DNEP表示NEP理论变化;Dfi表示该因子贡献的NEP理论变化。

2 结果分析

2.1 模型验证

与实测数据对比是验证模型的重要方式[34],本文选取已发表文献中不同地区的实测数据来进行模拟结果的验证,这些站点均分布于安徽省内或周边地区。根据不同地区实测数据的研究时段与植被类型,提取出相同地区相同时段的BEPS模型模拟结果进行对比,结果如表1所示。BEPS模型模拟结果与实测数据相近,平均误差为7.65%,可认为BEPS模型适用于模拟安徽省及周边区域生态系统的NEP。

表1 BEPS模拟的NEP和其他来源的NEP的对比Table 1 Comparison of BEPS-modeled NEP with NEP from other sources

2.2 环境植被因子变化特征

如图2所示,各环境植被因子在1982—2020年间表现出了不同的变化趋势。其中,Ta、VPD、LAI与CO2总体呈现显著的增加趋势(P<0.01),线性趋势变化率分别为0.04 °C/a、0.002 kPa/a、0.01 m2m-2a-1与1.87×10-6/a。Rad呈现显著的减少趋势(P<0.01),线性趋势变化率为-0.14 W m-2a-1。Pre变化趋势不显著(P>0.05),线性趋势变化率为1.27 mm/a。由相对变化率可知,LAI的变化最剧烈,相对变化率为24.28%,Pre的变化最缓慢,相对变化率为3.94%。其余因子的相对变化率排序依次为CO2(19.53%),VPD(19.11%),Ta(8.72%)和Rad(-4.28%)。

图2 1982—2020年安徽省环境植被因子的变化趋势Fig.2 Interannual variation of environment and vegetation factors in Anhui Province from 1982 to 2020Ta:气温 Air temperature;Pre:降水 Precipitation;Rad:太阳总辐射 Radiation;VPD:饱和水汽压差 Vapor pressure deficit;LAI:叶面积指数 Leaf area index;CO2:二氧化碳 Carbon dioxide

2.3 NEP时空格局

2.3.1NEP时间变化特征

如图3可知,1982—2020年安徽省总体表现为碳汇,多年年均NEP为651.14 gC/m2,最小值出现在2000年(603.29 gC/m2),最大值出现在2015年(712.73 gC/m2)。研究期间虽然2000年左右NEP出现小幅下降趋势,但是总体表现为显著的增强趋势(P<0.01),39年间的线性趋势变化率为1.10 gC m-2a-1。

图3 1982—2020年安徽省净生态系统生产力(NEP)年际变化 Fig.3 Interannual variation of NEP in Anhui Province from 1982 to 2020

2.3.2NEP空间变化特征

1982—2020年安徽省多年年均NEP的空间分布如图4左部所示,整体呈现“南北部较高、中部较低”的分布格局。NEP小于600 gC/m2的低值区主要分布安徽中东部地区,包括合肥市北部、马鞍山市、滁州市南部、宣城市北部和铜陵市西部。NEP为600—800 gC/m2的面积占62.98%,为安徽省NEP主要分布区间,集中分布在安徽省北部和西部。NEP大于800 gC/m2的高值区占13.16%,集中分布在安徽省南部的黄山市、池州市与宣城市南部。

图4 1982—2020年安徽省年均NEP及变化趋势的空间分布Fig.4 Spatial distribution of average annual NEP and change trend in Anhui Province from 1982 to 2020*表示通过0.05显著性水平

图4右部给出了1982—2020年安徽省NEP的空间变化趋势。由图可知,NEP无显著变化(P>0.05)的区域占40.12%,主要分布在安徽省西部与南部。NEP显著减少(P<0.05)的区域占7.11%,平均线性趋势变化率为-1.44 gC m-2a-1,主要位于安徽省西部与东南部。NEP显著增加(P<0.05)的区域占52.77%,平均线性趋势变化率为2.12 gC m-2a-1,主要位于安徽省北部、中部与东南部。

2.3.3安徽省各地级市NEP格局

图5反映了1982—2020年安徽省各地级市的NEP情况。整体而言,安徽省各地级市的NEP有增强趋势,但绝大部分地级市NEP都呈现“增减增”趋势,即在1980s至1990s有较小的增加趋势,随后出现下降趋势并在2000s以后又开始呈增加趋势。各地级市的多年均NEP最高为黄山市(874.14 gC/m2),最低为马鞍山市(547.49 gC/m2),而各地级市的NEP年均总量最高为六安市(7.56 TgC),最低为铜陵市(1.17 TgC)。

图5 安徽省各地级市的NEP大小及变化趋势Fig.5 Magnitudes and change trends of NEP in all cities of Anhui Province

2.3.4NEP重心分析

1982—2020年安徽省NEP重心移动的年际动态变化如图6。NEP纬度重心有显著的增加趋势(P<0.01),即安徽省NEP重心有向北转移的趋势,线性趋势变化率约0.0012°/a。NEP经度重心的减小趋势不显著(P>0.05),且相比于前2000年前,在2000年以后,经向重心波动较大。整体看来,安徽省显著北移的NEP重心变化与图4右部中NEP变化趋势一致。

图6 1982—2020年安徽省NEP重心经纬度的年际变化Fig.6 Interannual variation of longitude and latitude of NEP barycenter in Anhui Province from 1982 to 2020

2.4 NEP变化的归因分析

2.4.1敏感性分析

安徽省NEP对各环境植被因子的敏感性不尽相同,且在研究时段内表现出了不同的变化趋势(图7)。可以看出,安徽省NEP对CO2的敏感性系数的绝对值最大(0.41),其次分别为Ta(-0.40)、Rad(0.33)、LAI(0.32)和VPD(-0.21),NEP对Pre的敏感性系数绝对值最小(-0.014),即NEP对CO2的变化响应最为敏感,对Pre的变化响应最不敏感。

图7 1982—2020年安徽省NEP对影响因子的敏感性系数的时间变化趋势Fig.7 Time trends of the sensitivity coefficients of NEP to influencing factors in Anhui Province from 1982 to 2020

从敏感性系数的变化趋势可知,NEP对LAI与CO2的敏感性呈显著增强趋势(P<0.01),对Rad的敏感性呈显著减弱趋势(P<0.01),对VPD的敏感性呈显著的减小趋势(P<0.05),对Ta和Pre的敏感性线性变化趋势不显著(P>0.05),且对Ta的敏感性系数呈现“V”型变化。整体而言,安徽省NEP对LAI、CO2的正响应愈发敏感,对VPD变化的负响应也愈发敏感,对Rad变化的正敏感性降低,对Ta和Pre敏感性的变化不显著。

图8给出了安徽省NEP对各环境植被因子的敏感性系数变化趋势空间分布。NEP对Ta的敏感性显著增强的区域占12.30%(P<0.05),主要位于蚌埠市和黄山市,而显著减弱的地区占32.06%(P<0.05),主要位于阜阳市、六安市、芜湖市和宣城市。NEP对Pre和Rad和VPD的敏感性在大部分地区均表现为减小趋势,其中对Pre仅在1.61%的地区表现为显著减小(P<0.05),主要分布在六安市、安庆市和宣城市,对Rad在65.19%的地区表现为显著减小(P<0.05),主要分布在合肥市和马鞍山市以北的地区以及零散的南部地区,对VPD在48.76%的地区表现为显著减小(P<0.05),主要位于安徽省西北部、中部、东南部和西南部。NEP对LAI的敏感性在72.70%的地区表现为显著增强(P<0.05),主要分布在安徽省中部、北部和东南部,减弱的地区主要分布在南部,但大部分未通过显著性检验(P>0.05)。NEP对CO2的敏感性显著增强的地区占52.35%,主要位于安徽省西北部及东南部,显著减弱的地区占30.06%,主要分布在中西部(P<0.05)。

图8 1982—2020年安徽省NEP对影响因子的敏感性的空间变化分布Fig.8 Spatial changes of the sensitivity coefficients of NEP to influencing factors in Anhui Province from 1982 to 2020

2.4.2通径分析

以NEP为因变量,以各环境植被因子(Ta、Pre、Rad、VPD、LAI、CO2)为自变量进行通径分析。首先在进行逐步回归的过程中剔除了Pre,说明Pre对安徽省NEP的影响可以忽略,结果如表2和图9。直接通径系数绝对值排序依次为LAI >CO2> VPD >Ta>Rad,表明LAI对NEP的直接影响最大。间接通径系数之和的绝对值排序依次为VPD >Rad >Ta>LAI >CO2,其中LAI通过CO2路径产生的正效应和通过Ta、Rad与VPD路径的负效应相互补偿,导致LAI对NEP的间接影响较小,CO2间接通径系数较小亦是如此。VPD对NEP的直接影响和间接影响强度相当且作用方向相反,综合补偿导致VPD对NEP的相关系数较小;而Ta与Rad的间接影响强于直接影响,且对NEP的间接影响均主要通过LAI与CO2路径。

表2 安徽省NEP影响因子的通径分析表Table 2 Path analysis statistics of influencing factors of NEP in Anhui Province

图9 安徽省NEP变化通径分析图Fig.9 Path analysis of NEP changes in Anhui Province因子间的箭头表示因子的相关系数,因子与NEP的箭头表示直接通径系数;ε表示剩余因子

通径分析的剩余因子(ε)为0.45,计算得R2为0.79,说明所选的环境植被因子可以综合解释79%的NEP变化。NEP和各因子的相关系数绝对值排序依次为LAI >CO2> Rad >Ta>VPD,决策系数绝对值排序依次为LAI >CO2> VPD >Rad >Ta,其中Ta、Rad与VPD的决策系数均较小且为负数,相关系数未通过0.05显著性检验,表明这些因子对NEP的变化的综合影响能力较小,而LAI和CO2的决策系数大于0.47,且相关系数均通过0.01显著性检验,说明LAI和CO2为NEP变化的决定因子。

2.4.3贡献分析

各环境植被因子对NEP的相对贡献率计算结果如表3。根据R2可得,各环境植被因子综合解释了NEP占比79%的变化,与通径分析结果一致。结合图2和图7可知,LAI与CO2的显著增加(P<0.01)对NEP的增长呈正贡献,而Ta和VPD的显著增加(P<0.01)与Rad的显著减小(P<0.01)对NEP的增长呈负贡献。尽管NEP对Ta与Rad的敏感性高于LAI(图7),但Ta与Rad的相对变化率仅为8.27%与-4.28%,低于LAI(24.28%)与CO2(19.53%),导致其相对贡献率低于LAI与CO2的相对贡献率。由此可知,LAI与CO2的增加对NEP变化的贡献最大,相对贡献率均为25.23%,是NEP变化的主导因子,和通径分析结果一致。

表3 影响因子对安徽省NEP变化的贡献率Table 3 Contribution rate of influencing factors to NEP change in Anhui Province

图10左部给出了环境植被因子的综合决定系数(R2)空间分布。安徽省大部分地区的R2均大于0.6,面积占79.09%,其他20.91%的区域R2小于0.6,主要位于安徽省西北部、中部及南部,其中R2小于0.45的区域占3.76%。由此可知,环境植被因子能够较好的解释安徽省大部分地区的NEP空间变化。图10右部给出了安徽省NEP空间变化的主导因子分布图。其中,5.81%地区的NEP由气候因子主导,分布较破碎,并且大部分由Ta主导。LAI和CO2分别主导49.65%和44.54%地区NEP的变化,其中LAI为主导因子的地区主要分布在安徽省北部、中西部的大部分地区,CO2为主导因子的地区主要分布在安徽省西北部与东南部的大部分地区。

图10 NEP的决定系数(R2)及主导因子的空间分布图Fig.10 Spatial distribution of determinant coefficient (R2) and dominant factor of NEP in Anhui Province

为了进一步研究环境植被因子变化对安徽省NEP的影响,本文选取了NEP最低年(2000年)与最高年(2015年)作为典型年份,将当年各环境植被因子与多年平均值进行对比,并计算当年环境植被因子对NEP变化的贡献(表4)。结果表明,在NEP为最小值的2000年,除Ta偏高0.3℃外,其他环境植被因子均较多年平均偏低。仅VPD对NEP变化有正向贡献,其余环境植被因子对NEP均为负向贡献,尤其LAI偏低0.24,导致了NEP的大幅度减少。在NEP为最大值的2015年,Pre、LAI与CO2较多年平均偏高,Ta和Rad偏低。Pre与Rad对NEP变化有负向贡献,但远小于其他因子对NEP变化的正向贡献,特别是NEP对LAI、CO2与VPD变化的强烈正向响应使得2015年安徽省NEP取到最大值。

表4 典型年份的影响因子及NEP的变化统计Table 4 Statistics of changes of influencing factors and NEP in typical years

3 讨论

NEP受到诸多因子包括气候因子、植被结构、土壤条件和土地利用等影响,表现出较大的分异[40—41]。以往针对区域NEP的研究多关注气候因子的变化(如Ta、Pre等)[9—12],本文将植被因子与CO2浓度变化一并考虑,得出LAI与CO2为影响安徽省NEP的主导因子,这和Chen等[7]的全球研究结论一致。快速增长的LAI(图2)表明植被可以进行光合作用的部分增多,能够更加充分的接收和利用物质与能量,有利于生态系统的碳吸收。CO2对NEP的影响主要表现为CO2的施肥效应,其浓度增加既能加快羧化反应速率来促进光合作用[42],也能提高植被的水分利用效率[43],对植被的碳吸收有很强的促进作用[44]。但值得关注的是,随着CO2浓度的不断上升,这种促进作用可能会逐渐变小[45]。

另一方面,与其他区域如东北[9]和西北[10]相比,安徽省NEP与Pre无显著相关关系(P>0.05),因为安徽省处于亚热带湿润区,研究期间年均降水量基本超过900 mm,水分条件足够,因此,NEP对Pre的变化不敏感。有研究指出,当土壤水分充分供应时植被会出现由VPD驱动的气孔关闭[46],气孔导度的下降会阻碍CO2进入叶片,显著降低光合作用,限制系统的碳吸收[47]。本研究发现NEP对VPD的敏感性远大于对Pre的敏感性,且VPD是对安徽省NEP变化贡献最大的气候因子,与Konings等[48]对美国草地研究的结论一致。这表明,作为综合反映大气温湿状态的VPD,在湿润半湿润地区对生态系统NEP的影响可能大于降水的影响,其作用在未来研究中值得加强关注。

安徽省北部土地利用类型基本均为农田,根据安徽省统计年鉴[49],安徽省农田产量从4324 kg/hm2(1998年)增加至5926 kg/hm2(2020年),间接说明农田植被状况趋向更好。Chen等[41]指出,相比于其他土地利用类型,农田碳交换往往与植被因子更为密切相关。同时,本研究发现安徽省的LAI在西部和北部的大多数区域表现为显著增加趋势(P<0.05),在东南部小区域为显著减少趋势(P<0.05)(图11),这和NEP对LAI敏感性的变化趋势(图8)在空间上高度一致,这种敏感性变化与因子变化较为一致的配置,对安徽省西部和北部NEP的增加十分有利,导致安徽省北部NEP增长速率大于南部,故此,安徽省NEP重心明显北移。

图11 1982—2020年安徽省LAI变化趋势空间分布图 Fig.11 Spatial distribution of LAI change trend in Anhui Province from 1982 to 2020

本研究考虑的环境植被因子对NEP变化的综合解释能力为79%,表明21%的NEP变化未被解释,即仍有其他影响安徽省NEP的因子未被考虑。研究表明,人工对农田与森林的经营管理、树龄、氮沉降等因素也会影响生态系统的NEP变化[50—52],如氮沉降与CO2的协同作用有利于生态系统的碳吸收[53—54],也有学者关注到了物候,尤其是秋季物候对年际NEP变化的影响[55—56]。因此,在今后的研究中可以加以考虑,进一步明确NEP对其他环境因子的响应规律,以提高安徽省碳汇的解释能力。

4 结论

本文在采用经验证的BEPS模型模拟结果,分析了安徽省NEP时空变化及NEP对环境植被因子敏感性的时空变化,进一步量化了环境植被因子对NEP时空变化的相对贡献率,探究了安徽省NEP时空变化的驱动因子,得出以下主要结论:

(1)1982—2020年,安徽省多年年均NEP为651.14 gC/m2,线性趋势变化率为1.10 gC m-2a-1,总体呈现显著增加趋势(P<0.01)。在空间上,NEP表现为“南北部较高,中部较低”的分布,NEP显著减少(P<0.05)的区域占7.11%,主要位于安徽省西部与东南部。NEP显著增加(P<0.05)的区域占52.77%,主要位于安徽省北部、中部与东南部。NEP重心有显著的北移趋势(P<0.01)。

(2)敏感性分析结果表明,安徽省NEP对CO2的变化最为敏感,其次分别为Ta、Rad、LAI、VPD和Pre。LAI与CO2的增加对安徽省NEP有正向贡献,Ta与VPD的增加和Rad的减小对安徽省NEP有负向贡献,Pre无显著变化且对安徽省NEP的影响较小。空间分布上,NEP对气候因子(Ta、Pre、Rad和VPD)的敏感性总体为减少趋势为主,NEP对LAI与CO2的敏感性总体以增加趋势为主。

(3)通径分析结果表明,LAI与CO2对NEP主要表现为直接影响。其他环境植被因子中Ta与Rad对NEP的间接影响较大,且主要通过LAI与CO2路径实现,而VPD对NEP的直接与间接影响强度相当且方向相反,相互补偿后导致VPD与NEP相关系数较小。

(4)环境植被因子综合解释了安徽省NEP 79%的时空变化。其中,LAI与CO2是NEP变化的主导因子,相对贡献率均为25.23%。其余因子相对贡献率依次为VPD(-13.07%)、Ta(-10.85%)和Rad(-5.01%)。空间上,LAI和CO2分别主导了面积占比为49.65%和44.54%的NEP变化,其余因子的主导面积占比依次为Ta(4.55%),VPD(0.69%)和Rad(0.57%)。