黄土高原植被物候变化及其对季节性气候变化的响应

受气候变化影响,全球范围内植被物候发生了显著变化,而目前针对不同植被分区类型下（荒漠草原区、典型草原区、森林草原区、落叶栎林区、落叶栎林亚区）植被物候变化及其对季节性气候变化响应的研究尚少。因此基于MODIS遥感归一化差值植被指数（MODIS NDVI：MOD13Q1）数据、中国植被区划数据及135个气象站点插值数据,利用Sen''s斜率估计、Hurst指数和高阶偏相关分析等方法,研究黄土高原2001-2018年植被物侯变化及其对季节性气候变化的响应。结果表明：（1）黄土高原植被生长季始期（SOS,Start of Growing Season）主要集中在第96-144天,子植被分区由西北向东南方向,逐渐呈现提前趋势,71.0%的像元植被SOS整体提前0-2 d/10a （α=0.05）,且在未来一段时间66%的像元植被SOS继续呈现提前趋势;植被生长季末期（EOS,End of Growing Season）主要集中在第288-304天,各子植被分区植被EOS变化基本保持一致,87.6％的像元植被EOS整体延迟0-3 d/10a （α=0.05）,且在未来一段时间有80%的像元植被EOS继续呈现推迟趋势。（2）黄土高原植被SOS主要受各季节温度的影响;当年春季降水导致植被SOS提前,主要分布在黄土高原中部;上年夏季和上年秋季降水增加会导致植被SOS推迟;当年春季、上年秋季和年初冬季的温度升高均会导致植被SOS提前;各子植被分区植被SOS对不同季节降水的响应存在差异,而对不同季节温度的响应具有一致性。（3）黄土高原植被EOS主要受各季节降水和秋季温度的影响;不同季节降水增加均会导致大部分植被EOS推迟;当年秋季温度导致整体区域植被EOS推迟,且各子植被区植被EOS对当年秋季温度响应具有一致性。该研究可为大尺度植被物候影响因素提供新的认识,也为植被适应未来气候变化提供借鉴。     

黄土高原植被覆盖时空变化及其对气候因子的响应

为研究黄土高原地区退耕还林(草)后,植被覆盖变化及其对水热条件的响应,利用1999—2013年SPOT VGT NDVI 1km/10d分辨率数据,采用最大合成法、一元线性回归法和偏相关分析法,系统分析了黄土高原地区NDVI(归一化植被指数)的时空分布及变化趋势,及其与气候因子的关系。结果表明:黄土高原1999—2013年年最大NDVI的平均值为0.31,NDVI较高的区域位于黄土高原南部,而西北部植被覆盖度较低;自1999年开始,黄土高原地区NDVI呈极显著(P0.01)增加趋势,年最大NDVI的变化斜率为0.0099;不同季节(春、夏、秋、冬)和生长季的植被状况均呈现良性发展趋势;1998—2013年间,黄土高原地区气候呈现不显著的\"冷湿化\"特征;NDVI年际(及生长季和季节)变化与降雨和温度的相关性不显著,而在月时间尺度上,呈显著的相关性,并且月NDVI与当月降雨量的相关性要强于与当月温度的相关性;植被生长对温度的响应存在一个月的滞后期,而对降雨的响应无滞后效应。     

1982-2015年黄土高原植被变化特征及归因

黄土高原植被对于黄河中下游泥沙量减少以及区域生态安全维持发挥着重要作用。气候变化和人类活动是驱动黄土高原植被覆盖度变化的两大因素,然而先前研究多基于统计方法建立植被覆盖度与气候变化的线性关系,忽略了植被覆盖度与气候变化之间复杂的非线性关系,限制了对黄土高原植被变化驱动机制的理解。基于1982-2015年黄土高原植被覆盖度(GIMMS NDVI)数据,以1999年退耕还林工程开始实施为界限,采用随机森林和残差阈值法等手段量化并分析了1982-2015年植被变化特征及气候变化和人类活动对植被变化的影响。结果表明:(1)1982-2015年黄土高原植被覆盖度整体上呈现增长趋势,1999年之后植被覆盖度增加的面积大于1999年之前, 1999年以前显著增长的区域只占总面积的54.68%,1999年以后显著增长的区域面积占比增长至85.11%;(2)在考虑时滞效应的基础上,1982-1998年黄土高原超过80%的区域的逐月温度、降水、日照时数与植被覆盖度存在显著的正相关关系,表明1999年以前植被覆盖度变化主要与气候变化有关。随机森林算法能较好地模拟植被覆盖度与气候因子之间的关系,99.71%的区域拟合的决定系数R²达到0.5以上;(3)2000-2015年,气候因素和人类活动对植被覆盖度的影响具有空间异质性。人类活动对植被覆盖度产生积极影响的面积增加了55.94%,气候消极影响的面积减少了42.58%,气候积极影响的面积减少了12.86%。研究提出的量化气候变化与人类活动对植被覆盖度影响的方法能够更好的区分黄土高原植被覆盖度变化的驱动因素,为黄土高原生态建设以及黄河流域的生态安全和高质量发展提供数据支撑。     

祁连山不同植被类型的物候变化及其对气候的响应

基于1982—2006年GIMMS NDVI和2000—2014年MODIS NDVI遥感数据,利用double logistic拟合方法提取了1982—2014年祁连山区不同植被的生长季始期、生长季末期和生长季长度3个重要的物候参数,分析了不同植被物候期的时间变化趋势、空间分异特征及对气候因子的响应。结果表明:(1)祁连山区不同植被的生长季始期和生长季末期随年际变化表现出波动提前或推迟,其中沼泽植被的变化波动最大;草甸植被、灌丛植被、阔叶林植被和栽培植被生长季长度出现延长趋势;(2)祁连山区植被生长季始期集中在5月初,其中阔叶林植被生长季开始最早,荒漠植被生长季开始最晚,植被生长季末期集中在9月,栽培植被生长季结束较早,荒漠植被、沼泽植被生长季结束较晚,植被生长季长度集中在110—140 d,其中阔叶林植被、针叶林植被生长季长度较长,而荒漠植被、高山植被生长季长度较短;(3)祁连山植被物候期变化趋势的空间分布表明植被生长季始期、生长季末期主要表现为提前不明显和推迟不明显,生长季长度主要表现为缩短不明显和延长不明显;(4)物候要素与气候要素相关性表明前期温度的积累有利于植被的开始生长,但当年3月的降水量对植被生长季始期同样有重要作用,不同植被生长季末期与8月、9月温度相关性较大,而与10月、11月降水的相关性较大。     

黄土高原植被景观多尺度变化及其与地形的响应关系

基于RS和GIS技术,利用黄土高原近30年不同分辨率的归一化植被指数(NDVI),采用基于小波分析的多尺度空间统计学方法,研究不同时期黄土高原植被景观的多尺度变化特征及其与地形的响应关系。结果表明:(1)黄土高原植被景观1982—2011年期间发生了周期性变化。其中,1990年以前为植被恢复阶段,1990—2001年为植被退化阶段,2001年以后为植被恢复阶段。(2)植被景观的空间异质性与数据获取的时间、数据分辨率、空间位置和地形均有关。一方面,植被指数和数据分辨率越高,植被景观的空间异质性就越大,而且沿经度方向植被景观的空间差异性大于纬度方向;另一方面,NDVI的多尺度变化特征与地形因子有不同程度的相关性,表现为:高程地形湿度坡度坡向,这对于今后黄土高原植被景观的合理布局,提高生物多样性,控制水土流失,增强景观的连续性具有特别重要的意义,也可为指导生态环境建设提供基础资料。     

2000-2011年黄土高原植被景观格局变化

在RS和GIS技术的支持下,利用遥感数据MODIS NDVI提取植被覆盖度,运用景观生态学的分析方法,对黄土高原近12年不同时段(2000-2003年、2004-2007年和2008-2011年)以及夏季不同时期(初夏、盛夏和暮夏)的植被景观格局变化进行研究.结果表明:(1)近12年黄土高原的植被覆盖度在景观水平上,破碎度和复杂度在初夏时期呈先增后减的趋势;而在盛夏和暮夏时期整体趋于复杂化,景观优势度、最大斑块的聚集度降低,多样性先增后减,景观连通性无明显变化.(2)在类型水平上,年际变化主要表现为中低度、中度覆盖向中高度、高度覆盖类型的转化,而低度覆盖呈先增后减的变化趋势,变化幅度较小,约为3％～5％;季节内变化主要表现为低度覆盖持续减小,中度、高度覆盖类型在盛夏明显增大,最大增幅为9.76％,中低度覆盖类型在盛夏减小幅度最高,达13.77％,其他类型无明显变化趋势.研究结果揭示了黄土高原植被覆盖情况和景观格局变化规律,可为当地植被建设和生态环境保护提供参考.     

退耕还林后陕西省植被覆盖度变化及其对气候的响应

1999年起,陕西省实施了大规模的退耕还林、封山禁牧等生态建设和保护工程,使植被覆盖得到迅速恢复和增加.为了进一步跟踪评估植被覆盖变化,为生态建设和规划提供依据,本文基于2000-2017年MOD13Q1数据、气象数据,利用像元二分法估算陕西省18年间的植被覆盖度,通过空间插值方法、空间相关分析、统计学方法分区对其时空...     

2001–2020年黄土高原光合植被时空变化及其驱动机制

为揭示实施退耕还林(草)政策20年后黄土高原植被盖度的最新演变趋势及区域差异,定量分析气候和人类活动对该区植被盖度变化的贡献率及空间分布。该研究以光合植被(PV)盖度为植被生长状况指标,基于2001–2020年PV数据及同期气象数据,采用Mann-Kendall检验、Sen分析和残差分析等方法,分析了黄土高原2001–2020年植被覆盖的时空演变特征及其驱动要素。主要结果:20年中黄土高原植被盖度呈显著增加趋势,增速为每年0.8%。全区植被盖度呈增加趋势的区域面积比例为90%,呈显著增加的区域面积占比为71%;对全区植被盖度增加的贡献,主要是黄土丘陵区(约2/5),其次为风沙丘陵区(约1/4)和石质山区(约1/5);不同地貌分区内,黄土丘陵区中陕西榆林和延安两市区境内植被盖度增加迅速,风沙丘陵区中内蒙古鄂尔多斯市植被盖度变化最快;研究时段内人类活动和气候变化对黄土高原植被增加的贡献率分别为76%和24%;人类活动对植被盖度贡献较大区域主要分布在陕西延安以北、山西太原以南、宁夏同心以南和甘肃平凉和庆阳等丘陵、台塬和风沙丘陵等政府生态工程实施较好的地区。     

基于地理探测的黄土高原植被生长对气候的响应

为探讨黄土高原不同植被类型对气候变化的响应机制,以2002-2019年黄土高原归一化植被指数(NDVI)数据为基础,利用趋势分析、Hurst指数、地理探测器等方法分析不同植被类型NDVI变化趋势及其与气象因子的关系.结果 表明:2002-2019年,黄土高原不同植被类型NDVI以增长趋势和同向中持续性为主,仅栽培植被在...     

植被变化对气候的反馈机制及调节效应

植被通过光合作用固定大气中的CO₂来减缓温室效应,同时植被也通过改变地表能量收支影响温室效应。在过去的气候-植被研究中,大多关注气候变化对植被的影响,而植被对气候反馈的研究相对较少。植被通过调节地表能量收支、水通量等重要地气过程影响局地、区域乃至全球气候,在气候变化中的作用十分重要。因此,需要厘清植被对气候的反馈效应机制及其结果,并识别其地域差异。从生物地球物理和生物地球化学过程两方面分析植被与气候之间的作用机制,对全球及关键区域内植被变化对局地、区域乃至全球的气候反馈效应进行了系统总结：（1）生物地球物理反馈的区域特征明显,生物地球化学反馈则表现在全球尺度上,二者相互作用但难以统一;（2）植被破坏带来的气候影响在气温效应方面与生态系统的类型及地理分布相关：热带森林破坏带来增温效应,北方森林破坏带来降温效应,温带森林破坏则会通过增加森林反照率抵消丢失的固碳降温效应,气温效应表现不明显;（3）当前研究对关键过程机制考虑不够完善,不同研究方法的结果差异较大,且缺乏高质量观测数据的验证;同时考虑生物地球物理和生物地球化学的净气候反馈研究尚无法支撑植树造林对气候变化单一减缓作用的常规理解。本文可为科学评估植树造林对气候变化作用的方向与强度提供理论依据。     

鄱阳湖区域植被覆盖变化与降水相互响应关系

利用1982—1999年归一化植被指数(NDVI)数据集和鄱阳湖区域87个气象站点降水数据,采用年代差异变化、经验正交函数分析及相关分析方法,对鄱阳湖区域植被和降水的相互响应关系进行分析,从年内和年际不同尺度分别阐明植被和降水的相互影响时间和程度。结果表明:降水和植被的变化具有很好的空间一致性;多年平均表明,植被滞后于降水的变化;NDVI的年际变化对后期降水存在一定影响;相关分析显示,植被的生长变化滞后于降水变化6、7、8旬的相关性最高,相关系数分别为0.564、0.571和0.549,并且季节的相关分析发现,滞后一个季节的相关性最显著,且春秋季的降水对植被生长尤为重要。     

中国北方草原植被对气象因子的时滞响应

利用1982~1997年的气温、降水和1983~1997年生长季的NOAA/AVHRR的归一化植被指数(Normalized differential vegetation index, NDVI)遥感数据,分析了中国北方温带草原植被生长对气象因子的时滞响应。根据4个时间尺度(1~4个月)和4个时滞期(前0~3个月)将降水数据进行16种组合方式,计算了植被的NDVI与同期及前期(前1~6个月)降水之间的相关系数。同时,计算了植被的NDVI与同期和前一个月气温之间的相关系数。结果表明:1)中国北方温带草原植被的NDVI与同期降水和气温的显著相关。2)植被的NDVI对前一个月降水的时滞响应最强烈,植被的NDVI与当月降水和前两个月降水的累积量相关性最强。3)在生长季的起始阶段,去冬、今春的降水总量对草甸草原植被的生长有重要的作用。在生长季的中期和后期,当月和前一、二个月的降水对典型草原和荒漠草原的植被有显著影响。4)在草甸草原、典型草原区,生长季早期的气温均对植被生长的影响较为显著。在荒漠草原区,气温不仅在生长季初期与植被的NDVI呈现正相关,而且在生长季的中后期,气温与植被的NDVI呈现负相关性。     

基于NOAA-AVHRR数据的中国东部地区植被遥感分类研究

该文采用 19幅 (时间跨 8个月 ) 时间序列的NOAA_AVHRR的归一化植被指数 (NDVI) 最大值合成影像遥感数据, 经过主分量分析 (Principlecomponentanalysis, PCA) 处理后, 用非监督分类方法的ISODATA算法, 对中国东部地区的 (五省一市 ) 植被进行分类, 结果可以分出 2 8种土地覆盖类型, 除了两种类型为水体和城市或裸地外, 其余 2 6种类型均为植被类型, 根据中国植被分类系统, 这 2 6类可以归并为 6大植被类型 :1) 常绿阔叶林 ;2 ) 针叶林 ;3) 竹林 ;4 ) 灌草丛 ;5 ) 水生植被 ;6 ) 农业植被。用 1∶10 0 0 0 0 0数字化《中国植被图集》的植被类型检验遥感分类结果表明, 针叶林、灌草丛、常绿阔叶林和农业植被的分类具有较高的位置精度和面积精度, 位置精度分别为 79.2 %、91.3%、6 8.2 %和 95.9%, 面积精度分别达到 92.1%、95.9%、6 3.8%和 90.5 %。这 6大植被类型在地理空间上的分布规律与中国东部常绿阔叶林区植被的地带性分布基本一致。     

高原鼠兔洞穴数量与其栖息地植被分布格局的GIS分析

高原鼠兔(Ochotona curzoniae)是青藏高原高寒草甸生态系统的关键种,其生境利用特征对高寒草甸植被保护的影响一直是重要的研究内容.我们于2008年9月在四川省石渠县云波沟面积为44 km2的研究区域内进行了高原鼠兔有效洞穴数量样线调查.该地区一直以来放牧压力较小,2006年起执行国家西部退牧还草政策后.放牧完全停止.调查期间共完成样线7条,总长度19 km,获得取样单元2,423个.根据中巴资源2号卫星2008年9月18日拍摄的遥感图像计算并提取研究区域归一化植被指数(normalized difference vegetation index,NDVI)数据.利用石渠县地形数字地图建立数字高程模型计算坡度.用Bonferroni z检验比较取样单元内鼠兔洞穴数量与NDVI以及洞穴数量与坡度间的相互关系发现:NDVI值较低区间(0.1-0.3)内鼠兔洞穴数量显著低于理论预期值,中等NDVI值区间(0.3-0.4)内鼠兔洞穴的数量与理论值没有显著差异,而NDVI值较高区间(0.4-0.6)内,鼠兔洞穴的数量显著高于理论预期.回归分析结果表明:NDVI与洞穴数有显著正相关的线性回归关系y=16.50x+0.87(R2=0.78,P<0.001);鼠兔有效洞穴数量和坡度没有显著的线性关系(R2=0.04,P=0.441).因此,我们认为在无家畜竞争的自然环境中,栖息地植被的质量是影响高原鼠兔空间分布的重要因素.高原鼠兔有选择地利用植被较好的生境而避免使用植被过度退化的生境,这种利用方式在客观上防止了栖息地植被的过度破坏.     

基于NOAA/AVHRR NDVI 监测中国北方典型草原的生长季及变化

 该研究基于Savitzky_Golay滤波算法平滑了1982～1999年NOAA/AVHRR NDVI 时间序列影像,然后利用曲线拟合了锡林郭勒典型草原1982～1999年的每年物候期（返青期、黄枯期）及18年的平均物候期和生长季长度,并对1982～1999年的物候期进行了线性拟合,从而分析了物候期的变化趋势。结果表明：1）1982、1986、1992年的返青期处于正常水平,1985、1988、1989、1991年大部分地区的牧草返青期比正常年份有所提前。1984、1990、1993 年的黄枯期处于正常水平,大部分年份的黄枯期主要处于9月下旬至10月上旬（290～310 d）。2）在整个典型草原,返青期有较大的变异性,而黄枯期变化表现出了锡林郭勒典型草原的西南部较早、中部及东北部较晚的格局,生长季长度的变化格局为西南地区最短,中部地区最长。3）从1982～1999年,不同的地区表现出物候期及生长季长度提前或延迟的变化趋势,返青期大多数地区延迟时间集中在10～20 d,提前日期主要集中在10 d之内。锡林郭勒盟西南地区的黄枯期提前趋势最大。大部分地区的生长季长度变化呈缩短趋势,缩短日期小于 10 d,少部分地区的生长季延长,延长日期主要集中在0～10 d。4）对锡林浩特的物候期研究表明,牧草返青期提前日期小于10 d,黄枯期延迟大约14±5 d,生长季长度延迟大约1 5±5 d。最后利用野外观测数据对锡林浩特牧草返青期的拟合精度作出了评价。     

近20年藏北地区AVHRR NDVI与气候因子的关系

利用藏北那曲地区1981~2001年NOAA/AVHRR的旬合成NDVI资料和6个站的逐日气象资料,分析了(归一化指数NDVI)的年内和年际变化规律以及NDVI与8个气候要素的相关关系,主要结论:影响NDVI年变化最显著的气候因子是温度,其中水汽压与NDVI的相关程度明显大于降水量;日照时数与NDVI呈负相关;NDVI与日照时数的滞后时间约0~10d,与风速没有滞后现象,与潜在蒸散、温度、水汽压和降水约20~40d;影响NDVI年际变化最显著的气候因子是潜在蒸散量。     

基于时序植被指数的华北地区作物物候期/种植制度的时空格局特征

随着气候变化对农业系统的影响不断加剧,为保障粮食安全,须掌握变化中的自然与人文因素,而农业生态系统中变化最为明显直观的是农作物物候特征,如何提取大区域尺度上农作物物候期以及种植制度的时空格局特征,是评价区域粮食安全的重要因素.基于多时相遥感信息可以有效反映年内/年际农作物物候特征变化的原理,首先利用近10a来的SPOT/VGT-NDVI时间序列数据,在进行数据序列平滑重构处理基础上,提取了华北地区农作物典型物候期的数量分布与空间格局特征;然后,基于上述物候期的分异特性建立了一年一熟和一年二熟等种植制度类型的遥感识别标志;最后,重点分析了上述种植制度的空间格局及其时间波动特征,并利用农业统计资料对提取结果进行了简单验证.分析结果表明,作物物候期特征的数量分布和空间分布在不同生长季均具有显著差异,直接体现了与外界环境条件(诸如区域温度、降水和光照等)的匹配程度以及作物类型自身的生长特征;从主要种植制度空间分布来看,华北地区南部地区农作物类型以夏收作物和二熟秋收作物为主,与之对应的农田种植制度以一年两熟为主;华北地区北部主要为一熟制区域,作物类型以一年一熟秋收作物为主,作物种植制度空间分布随着纬度递减呈现出简单到复杂的总体趋势;从近10a的种植情况来看,一年一熟作物种植面积最大,年际变化幅度亦较大,一年二熟的夏收作物种植比例次之,而年际变幅最小,二熟秋收作物比例最低,其年际变幅居中.研究中亦提出,在进一步加强多时相遥感技术监测大区域农业生态系统动态变化的同时,亦需深入探讨作物物候特征及种植制度变化的驱动因素及其对国家粮食安全的影响.     

中国草地植被生物量及其空间分布格局

 草地生态系统是陆地生态系统分布最广的生态系统类型之一,它在全球变化中的作用越来越受到重视。利用中国草地资源清查资料,并结合同期的遥感影像,建立了基于最新修正的归一化植被指数(NDVI)的我国草地植被生物量估测模型,并利用该模型研究了我国草地植被生物量及其空间分布特征。结果表明：草地植被地上生物量与当年最大NDVI值具有很好的相关关系,两者可以用幂函数很好地拟合(R2=0.71, p<0.001)。我国草地植被总地上生物量为146.16 TgC(1 Tg=1012 g),主要集中在北方干旱、半干旱地区和青藏高原;总地下生物量为898.60 TgC,是地上生物量的6.15倍;而总生物量是1 044.76 TgC,占世界草地植被的2.1%～3.7%,其平均密度约等于315.24 gC·m-2,低于世界平均水平。我国草地植被单位面积地上生物量水平分布趋势为：东南地区高,西北地区低,与水热条件的分布趋势一致;从垂直分布看,在海拔1 350 m和3 750 m处分别出现了波谷和波峰,与我国特有的三级阶梯地势有着密切的关系。此外,我国草地植被生物量为森林的1/4左右,显著大于世界平均水平,说明我国草地在碳平衡中的贡献相对较大。     

基于近20年遥感数据的藏北草地分类及其动态变化

利用1982—2000年NOAA/AVHRR的旬合成归一化植被指数(NDVI)资料,采用主成分分析和非监督分类方法对藏北那曲地区植被进行分类,分析不同草地类型代表像元的NDVI年内和年际变化特征;定义那曲地区牧草主要生长期平均NDVI≥0.1的地区为植被区,NDVI<0.1的地区为植被稀少区,进一步分析植被区每个像元NDVI的时空变化特征.结果表明:该地区草地类型可分为高寒草甸、高寒草甸草原、高寒草原和高寒荒漠,分类结果与实际情况相符.4种草地类型的NDVI年内变化均呈单峰型,年最大值出现在8月;近20年来那曲地区植被区7—8月平均NDVI由东南向西北逐渐减少,约在0.1~0.6之间变化,变异系数在0.05~0.40之间.NDVI高的地区,变异系数相对较小;NDVI低的地区,变异系数相对较大,年变率的范围在-0.005~0.008之间.近20年来那曲地区植被变化不明显,约20%地区(主要分布在西部的尼玛和东部的嘉黎、比如、索县和巴青等县)的植被活动在减弱.     

近20年来中国植被活动在增强

为阐明近20年来中国植被覆盖变化的整体状况, 利用归一化植被指数(NDVI)作为植被活动的指标, 使用第3代NOAA-AVHRR/NDVI时间序列数据, 研究了1982~1999年间中国地区NDVI的变化. 为消除地表非植被因素的影响, 参考国际惯例, 定义年NDVI≥0.1的地区为有植被覆盖地区(简称植被地区), NDVI < 0.1的地区为植被稀少地区. 结果表明, 18年来, 我国大多数地区的NDVI都呈现不同程度的增加趋势, 表明我国的植被活动在增强. 与80年代初相比, 90年代末植被地区的面积增加3.5%, 植被稀少地区的面积下降了18.1%. 全国平均年NDVI增加了7.4%. 生长季节的延长和生长加速是我国NDVI增加的主要原因, 而温度上升和夏季降水量的增加以及农业活动的加强可能是其主要的驱动因子. 我国NDVI变化趋势显示了较大的空间异质性: 东部沿海地区呈下降趋势或变化不明显; 农业产区和西部地区增加显著. 这种空间异质性是由于城市化过程、农业生产活动、区域气候特征以及植被对气候变化的区域响应等综合因素作用的结果.