0511 马雯思1,3,4,马 超1,2,刘玮玮1
(1.河南理工大学 测绘与国土信息工程学院,河南 焦作 454000; 2.河南理工大学 矿山空间信息国家测绘与地理信息局重点实验室，河南 焦作 454000; 3.中国科学院西北生态环境资源研究院 冰冻圈科学国家重点实验室，甘肃 兰州 730000; 4.中国科学院大学，北京 100049)
摘 要:过去30 a，全球归一化差值植被指数(NDVI)呈增加趋势，植被返青期(SOS)提前，枯黄期(EOS)滞后，生长期(LOS)延长。作为人类经济活动最密集的工矿区，其生态环境问题是环境与发展的焦点问题之一，在全球变化的背景下，研究矿区及周边地区NDVI趋势，能很好地反映区域尺度上非自然生态区植被变化及对全球变化的响应。基于GIMMS AVHRR NDVI3g长时间序列遥感影像，利用IDL编程、线性回归、趋势拟合提取研究区1982—2013年植被覆盖的物候信息，分析矿区、缓冲区(10 km,20 km)及生态校验区月度、季度、年际变化情况，推算出植被SOS,EOS与LOS长度变化趋势，并结合气象观测记录分析气候驱动因子(降水、气温)与NDVI变化相关性。时序相关性分析表明:32 a彬长矿区年均NDVI总量随时间序列缓慢升高，32 a增长百分比达13.31%;生态校验区NDVI总量随时间序列升高趋势明显，32 a增长百分比达19.45%;受人类活动影响，矿区NDVI阶段增长百分比明显低于生态校验区。空间相关性分析表明:彬长矿区NDVI原始基底值高于校验区，植被SOS提前6±3 d，EOS滞后5±3 d，LOS延长11±3 d;生态校验区植被SOS提前3±3 d，EOS滞后3±3 d，LOS延长6±3 d;矿区植被LOS延长天数高于生态校验区5 d。气候相关性分析表明:彬长矿区处于半湿润气候区，植被对气温变化的响应高于降水，植被生长对降水有2～3 a的滞后性，NDVI的增加与该区气温升高、降水减少的共同作用有关。研究认为，采矿活动使得矿区及其周边地区植被活动放缓，年均NDVI增长率明显低于生态校验区;在全球变化作用下，植被LOS延长的基础上，人类活动使得矿区植被LOS被再度延长。
关键词:GIMMS AVHRR NDVI3g;彬长矿区;生态校验区;采矿活动;气候变化;生长期
中图分类号:TD8
文献标志码:A
文章编号:0253-9993(2019)04-1197-10
收稿日期:20180511
修回日期:20181122
责任编辑:韩晋平
基金项目:国家自然科学基金委员会与神华集团有限责任公司联合基金资助项目(U1261206，U1261106)
作者简介:马雯思(1993—)，女，河南开封人，博士研究生。E-mail:969259520@qq.com
通讯作者:马 超(1967—)，男，内蒙古克什克腾旗人，教授，博士。E-mail:d_insar@126.com
MA Wensi1,3,4,MA Chao1,2,LIU Weiwei1
(1.College of Surveying & Land Information Engineering,Henan Polytechnic University,Jiaozuo 454000,China; 2.Key Laboratory of State Bureau of Surveying and Mapping of Mine Spatial Information Technology,Henan Polytechnic University,Jiaozuo 454000,China; 3.State Key Laboratory of Cryospheric Science，Northwest Institute of Eco-Environment and Resources，Chinese Academy of Sciences，Lanzhou 730000，China; 4.University of Chinese Academy of Sciences，Beijing 100049,China)
Abstract:In the past 30 years,the global normalized difference vegetation index (NDVI) showed an increasing trend,the start of the vegetation growing season (SOS) was advanced,the end of the growing season (EOS) was delayed and the length of growing season (LOS) was prolonged.As the most intensive industrial and mining area of human economic activities,its ecological environment is one of the focal issues of environment and development.Under the background of global climate change,the NDVI trend in the mining area and around can well reflect vegetation change in non-natural ecological zones and its response to global change at regional scale.Based on the long-term sequence of GIMMS AVHRR NDVI3g remote sensing images,the phenological information of vegetation cover was extracted by using IDL programming,linear regression and trend fitting in the study area from 1982 to 2013.Analyzing the monthly,quarterly and annually variations of the mining area,buffer zone (10 km,20 km) and its calibration area,the variation trends in SOS,EOS and LOS length were obtained,and finally the correlation between the climatic driving factors (i.e.precipitation,temperature) and the change of NDVI was analyzed combined with meteorological records.The temporal-correlation analysis demonstrates that the average annual NDVI of Binchang mining area increased slowly with an annual growth rate of 13.31% over the past 32 years.But the total NDVI in natural ecological check area (CK) showed apparent growth trend with an annual growth rate up to 19.45%.Therefore,it can be concluded that the annual growth rate of NDVI in mining areas is significantly lower than that in the CK under the influence of human activities.Spatial-correlation analysis indicates that the original NDVI value of the Binchang mining area is higher than that of the CK.In the mining area,the vegetation SOS is 6±3 days ahead of time,the EOS lags behind 5±3 days and the LOS extends 11±3 days.The corresponding values in the CK area are 3±3 days,3±3 days and 6±3 days,respectively.The extension of LOS in the mining area is longer than that in CK for 6 days.Climate-correlation analysis shows that Binchang mining area is located in the semi-humid climate zone.As a result,the response of vegetation to temperature change is higher than to precipitation,and vegetation growth has a lag of 2 to 3 years in response to precipitation.In other word,the increase of NDVI depends on the interactions of temperature rise and precipitation decrease.The study concludes that the mining activities slow down the vegetation growth in the mining area and around,and the annual growth rate of NDVI is obviously lower than that in the CK.Under the effects of global change and based on the extension of LOS,the vegetation LOS extends further in the mining area due to the human activities.
Key words:GIMMS AVHRR/NDVI3g;Binchang mining area;checked area (CK);mining activities;climate change;length of growing season(LOS)
移动阅读
马雯思,马超,刘玮玮.陕西彬长矿区NDVI3g(1982—2013)变化趋势及气候响应[J].煤炭学报,2019,44(4):1197-1206.doi:10.13225/j.cnki.jccs.2018.0633
MA Wensi,MA Chao,LIU Weiwei.Variation and climate response of NDVI3g(1982—2013)in Binchang mining area of Shaanxi Province[J].Journal of China Coal Society,2019,44(4):1197-1206.doi:10.13225/j.cnki.jccs.2018.0633
煤炭资源的开发利用造成了地表沉陷、土地挖损、煤矸石堆存占地等无法弥补的损失，进而造成水土流失加剧、耕地退化、土地荒漠化等更为严重的生态环境问题[1]。矿区生态环境的破坏主要体现在对水、土壤、植被大气的影响，绿色植被作为生态环境最敏感和最主要的环境要素，是矿区生态环境恶化的“指标器”。遥感植被指数可以定量分析植被覆盖及其生长状况，其中归一化差值植被指数(NDVI)为最常用的一种。
利用遥感数据动态监测自然地表植被覆盖变化，国内外已有诸多研究，如Myneni,Fang,Tucker,Los等利用卫星遥感数据对全球范围内植被覆盖变化，研究表明:过去20 a北半球高纬度地区陆地植被活动在显著增强[2-7]。方精云等利用AVHRR/NDVI时间序列数据(1982—1999年)研究中国地区NDVI变化，结果显示:全国平均归一化植被指数(NDVI)年增加7.4%，中国西北、陕西地区植被指数增加高于6.5%[8]，NING等运用归一化植被指数及气候相关因素对中国北方黄土高原研究表明:过去几十年里，年均NDVI最大值上升速率达0.09/10 a，荒漠化面积由54.7%(1999年)降至8.9%(2012年)，气温、降水变化是植被生长状况改变的主要诱因[9]。
针对矿区人为影响植被监测的研究相对较少:RASIM等运用两种植被指数(归一化植被指数NDVI和归一化湿润指数NDWI)对加拿大阿萨巴斯卡油砂区进行检测，研究证实:1984—2012年，NDVI平均减少18.6%，归一化水体指数(NDWI)平均减少31.0%，是气候变化和人为因素的干扰结果[10]。PROSPER等利用遥感和地理信息系统技术评估佤邦东部矿区土地覆盖变化:植被出现下降趋势，NDVI像元均值从0.48降至0.11[11]。HUANG等利用时空自适应反射融合模型对山西大同矿区进行NPP和生物量损失估算表明:2001—2010年，矿区年NPP平均减少24.71 gC/m2，年生物量损失率达33.48 gC/m2，气候变化与平均NDVI值呈正相关，采煤对植被的干扰导致生物量的损失及植被吸收能力的减弱[12]。徐占军等利用Landsat5 TM数据对徐州矿区研究表明，采矿活动对NPP 的影响大于气候变化对NPP的影响，NPP变化对采矿活动具有敏感性[13];郝成元等基于EOS/MODIS卫星遥感数据研究分析，研究时段内(2001—2006年)潞安矿区NPP呈显著减少趋势，其时间异质性主要与煤炭开采、农业耕作等类活动联系密切[14];马超等采用SPOT2/4卫星遥感数据研究采矿扰动区NDVI变化趋势，潞安矿区8个试验工作面在开采前NDVI相对稳定，开采后出现下降，2004—2007年期间，年均降幅达14.83%，认为采矿扰动是矿区地表NDVI变化的主要诱因[15]。
鉴于国外对非自然生态区植被遥感监测研究较少，且国内的多数研究时序间断或过短，缺乏空间对比参考，也没有考虑全球变化的影响，难以量化开采扰动对矿区植被的影响程度。笔者采用时间分辨率较高的GIMMS AVHRR NDVI3g(1982—2013年)长时间序列数据，针对缺乏长期监测的煤炭开采区，设置多重对比实验区(10 km缓冲区、20 km缓冲区、自然生态校验区)，结合气候驱动因子(气温、降水)，分别从时间相关性、空间相关性、气候相关性3方面揭示矿区开采前、后地表植被覆盖的变化特征及植被生长期的变化规律，通过对比分析开采区相对于非开采区的植被变化差异，甄别自然变化和人类活动的影响，旨在为采煤区植被自然恢复、人工修复及其周边地区生态综合治理提供借鉴。
1.1.1 直接影响区
彬长矿区位于黄陇侏罗纪煤田中部、鄂尔多斯盆地南端、太峪背斜以北、彬县和长武的交界处。坐标为东经107.65°～108.33°、北纬34.98°～35.32°，南北距离约36.5 km，东西距离约60.6 km，总面积978 km2。下辖大佛寺井田、小庄井田、文家坡井田等13座井工矿。研究区属渭北黄土高原沟壑地貌，地貌类型以黄土梁塬、黄土沟谷为主。其中黄土沟谷(49.79%)面积最大，其次为黄土梁(36.59%)，黄土塬(6.02%)面积较小，剩余为河流阶地、漫滩及水体。研究区处于暖温带半干旱大陆性季风气候区，年平均气温 11.1 ℃，年平均降水量为 561.3 mm。区内植被以中覆盖度为主，植被类型较为单调，主要包括刺槐林、白羊草、长芒草和农业植被，研究区土地利用类型以草地(44.63%)和耕地(33.06%)为主，林地(11.66%)次之。
1.1.2 间接影响区(缓冲区)和自然生态校验区
煤炭开采不仅对直接影响区域造成一定生态影响，也会影响到周边区域生态系统，为更好的定性分析开采扰动对矿区植被覆盖的影响，笔者增设了10 km缓冲区、20 km缓冲区、自然生态校验区[16] (图1)。其中，在选择自然生态校验区时，参照了遥感“伪不变特征区”(Pseudo Invariant Feature areas,PIFs)的设置原则[17]，即:① 必须是植被覆盖区;② 与研究区均处于同一纬度，尽可能远离扰动区，受人为因素影响较少;③ 与研究区高程相近，太阳光照条件相近。两区经差0.64°，纬差0.94°，几何中心相距150.5 km，鉴于两区原始植被覆盖类型并不一致，因此通过NDVI的相对变化程度(即增长率)来进行对比分析。
图1 研究区与生态校验区的地理位置
Fig.1 Geographical location of the study area and the CK
1.2.1 遥感数据
研究所用的遥感数据为GIMMS NOAA/AVHRRR 的第1,2通道生成的NDVI3g数字影像，该NDVI数据最初是由美国航天局(NASA)全球监测与模型研究组发布的最大合成(Maximum Value Composites,MVC)数据(ftp://ftp.glcf.uniacs.umd.edu/glcf/GIMMS/)。图像分辨率为8 km×8 km，时间分辨率为15 d，时间跨度为1982-01—2013-12，空间范围为70°～140°E，15°～55°N，32 a涵盖了768幅影像。经过辐射校正和几何粗校正的NOAA/AVHRR，再进一步消除云、大气、太阳高度角等的部分干扰，运用国际通用的MVC(最大合成)方法获得半月合成数据。相对于其他NDVI数据源来说，GIMMS NDVI3g时间序列数据适合于长时间、大尺度植被活动变化的分析，用于分析我国的植被活动状况是合适的[8]。
1.2.2 气象数据
气候数据为1982—2013年各省、市、自治区气候资料处理部门逐月上报的“地面气象记录月报表”的年均气温和降水资料，由于研究区域内的气象站点较稀疏，为消除极端数据所带来的影响，笔者结合当地气候地貌环境，以矿区几何中心点为圆心，300 km为半径作一包含34个气象站点的圆形插值区进行插值，土地利用图取自国家基础地理信息中心的30 m地表覆盖数据(GlobeLand30)，矢量文件为1∶25万县界数字线划图。
1.3.1 遥感数据预处理
NDVI半月合成数据再处理方法主要有累加处理、平均值处理与最大值处理。累加处理是将分析某时间段内的NDVI值求和，可反映植被在某时间段内的生物量当量累加值;平均值处理是将某时间段内的NDVI值取平均值，该方法可消除气候异常对NDVI值的影响;最大值处理消除了云、雾等自然条件对NDVI值的影响，是将某时间段内的NDVI值取最大值。矿区的地理位置、所在区域的降水量、气温、以及作物生长状况等条件共同决定合成方法的选择。低植被覆盖度的半干旱地区(年均降雨量<400 mm,如神东矿区)，为突出植被覆盖状况，可取全年NDVI的最大值[18];中等植被覆盖度的半湿润地区(400 mm<年均降雨量<800 mm,如抚顺矿区)，可取植被生长期内的累加值[19-20];高植被覆盖度湿润地区(800 mm<年均降雨量,如淮南矿区)，可取全年NDVI的平均值[21]。
(1)累加值处理。研究区属中度植被覆盖半湿润地区，宜采用累加值方法，运用累加值也可估算区域NDVI总量年际变化，为开采扰动对矿区带来的植被变化提供一定依据，计算数据模型为
(1)
式中，ai为第i个像元NDVI半月合成值;m为累加次数，分别为2(月总值)、6(季总值)、24(年总值);bi为第i个像元累加值;n为像元个数。
(2)线性回归方程。STOW等用采用一元线性回归来计算植被的绿度变化率(Greenness Rate of Change,GRC)，GRC被定义为某时间段内NDVI一元线性回归方程的斜率S[22]。该斜率揭示了在一定时间内每个像元所在地区植被指数的变化趋势，公式为
(2)
其中，xi为自变量，
为因变量。当S>0时，表明该像元所在区域植被状况变好，反之则是变差。
(3)变化率。在单位时间内的变化量就是变化率。方精云[8]提出用变化率判断区域植被指数增加或减少的程度大小，并约定NDVI的变化率为
(3)
式中，S为32 a 的NDVI累积值线性回归直线的斜率;M为32 a的NDVI累加值均值，该参数用以对比分析开采扰动区与非开采区植被活动的趋势。
(4)植被生长期阈值。利用32 a的NDVI3g月均值数据，计算反映植被物候的3个参数:返青期(SOS)，枯黄期(EOS)和生长期(LOS)。月度曲线突然升高和降低处即为生长期开始或结束临界点，生长期长度为枯黄期开始时间减去返青期开始时间[23-24]。因原始数据月NDVI值时间分辨率为30 d，为进一步精确计算植被生长期，利用Matlab对32 a的月NDVI值以3 d为单位进行插值，可以计算出0.1个月的时间变化，即有效计算精度为±3 d。
1.3.2 气象数据预处理
运用反距离加权法(Inverse Distance Weighted Method，IDW),对插值区34个气象站点年平均降水、气温进行插值，得到研究区32 a年均降水及气温(图2)。
(1)矿区年均降水分析。如图2(a)所示，彬长矿区32 a降水量波动较大，1982—2002年降水量整体呈下降趋势，与全国降水量变化趋势一致[25]，2003年后，降水量有上升之势。研究区平均年降水量在1997年降到极小值418.6 mm，1983年出现最高值818.3 mm，1981—2012年均降水量为592.56 mm，降水下降速率为21.55 mm/(10 a)。
(2)矿区年均气温分析。32 a研究区年均气温最低为9.81 ℃(1984年)，最高为12.12 ℃(2006年)和12.05 ℃(1998年)，32 a平均气温为11.11 ℃(图2(b))。矿区升温速率为0.47 ℃/(10 a)，高于陕北地区增温速率0.29 ℃/(10 a)[26]，更高于全国近50 a(1951—2004年)的平均增温速率0.25 ℃/(10 a)[27]，远高于同期全球气温变化速率0.12 ℃/(10 a)[28-29]。
图2 彬长矿区32 a气候变化趋势
Fig.2 Trend of average temperature & precipitation for 32 years of Binchang coalfield
2.1.1 研究区月度NDVI值比较
(1)高斯拟合
将每个月的NDVI求平均，获得植被NDVI的年内特征曲线，图3为4个区域NDVI月均值拟合曲线，可见月度植被指数动态变化曲线符合高斯分布:
其中，y0为基线偏移;A为曲线下方的积分面积;x0为中央峰值;w=2σ近似于峰值半高宽的0.849，4个区域的拟合度R2分别为0.910 0(矿区)，0.921 8(对比区，也称校验区)，0.950 7(10 km缓冲区)，0.927 5(20 km缓冲区)。
图3 研究区32 a月度NDVI趋势
Fig.3 Trend of monthly NDVI for 32 years in the study area
彬长矿区地处彬县与长武县境内，当地以落叶阔叶林和农业植物为主，年内变化曲线整体呈单峰型，1—3月份NDVI值处于低值期，4月份后增加迅速(NDVI值陡升)，5—8月达到并保持峰值状态，9月底至10月初NDVI值大幅度下降，直至2月份到达谷底。此变化特征与陕西省植被年内变化趋势相一致[30]。
(2)生长期分析
由图4可以判断，矿区和缓冲区(10 km,20 km)NDVI基数较高，有些年份12月时仍保持0.25的NDVI值，而且0.3的曲线极不稳定，其返青期和枯黄期宜采用相同阈值0.35;校验区受人为扰动较少，原始植被覆盖类型为草地，其返青期和枯黄期宜采用相同阈值0.20。
图4 研究区时间序列月NDVI晕渲图
Fig.4 Shading map of monthly average NDVI in the study area
① 图4的“倒梯形”柱状结构表明，4个研究区返青期均有提前趋势，枯黄期有滞后趋势，生长期有延长趋势。② 对1982—2013年矿区、10 km缓冲区、20 km缓冲区和校验区植被返青期、枯黄期、生长期进行一元线性回归，得到4个区琙植被生长期变化趋势(图5)。32 a四区生长期均有延长趋势，生长期分别为182 d(矿区)，185 d(10 km缓冲区)，188 d(20 km缓冲区)和180 d(校验区)，非开采区生长期延长13±3 d(10 km缓冲区)，8±3 d(20 km缓冲区)，6±3 d(校验区)，矿区生长期延长11±3 d，线性增长率为0.34 d/a，1984年为矿区生长期最小年(156 d)，生长期最大年是2004年(175 d)。
图5 4类区域植被生长期变化
Fig.5 LOS change in four types of regions
2.1.2 年NDVI累加值比较
对1982—2013年裁剪、去除DN值后的遥感图像运用IDL进行累加值处理，得到32 a累积值遥感图像，进而提取出每年各个区域的NDVI累积值，可反映不同研究区植被生长状况，为矿区植被自然恢复、人工修复提供参考。
图6 32 a NDVI累加值变化趋势
Fig.6 Slope of NDVI during 32 years
根据图6，研究区年NDVI总值对比趋势显示:
(1)通过0.01的显著性水平检验，矿区与10 km缓冲区、20 km缓冲区、校验区NDVI的相关系数分别为0.979 1,0.961 2,0.651 3;说明矿区与缓冲区植被响应高度一致，而矿区植被活动变化趋势与校验区相关性不大，这样更加适合相对变化的比较。32 aNDVI总量从高到低依次为20 km缓冲区(358.33)、10 km缓冲区(237.42)、直接影响区(158.43)、生态校验区(121.93)，4区植被覆盖均在波动中呈增加之势，其中矿区植被活动变化幅度最弱。2006年开始，校验区植被累积值开始逐渐缩小与矿区的差距，至2013年两区年际植被总量基本一致，在自然生态环境条件下，校验区植被自然生长，矿区植被受到人类活动影响，导致两区植被总量差距缩小。
(2)对研究区32 a连续时间序列NDVI累积值做一元线性回归，由回归方程(式(2))得到的绿度变化率(S)分别为0.659 2(矿区)、0.809 0(10 km缓冲区)、1.174 8(20 km缓冲区)、0.741 1(校验区)。32 a NDVI总量均值分别为158.43，237.42，358.33，121.93，结合年变化率式(式(3))，可得研究区32 a NDVI增长百分比为10.49%(20 km缓冲区)<10.9%(10 km缓冲区)<13.31%(矿区)<19.45%(生态校验区)，矿区因受采煤影响植被增长率显著低于无人工开采干扰的自然生态区。若以校验区NDVI上升速率为自然增长率，可初步认为亏欠的部分是采矿活动的影响，则开采活动对矿区、10 km缓冲区和20 km缓冲区NDVI增长率的贡献分别为-6.14%，-8.55%，-8.96%，对覆盖度高的区域影响似更为显著。
降水和气温是植被响应气候变化最主要的2个驱动因子，笔者选取32 a气象数据与研究区NDVI做相关性统计分析(图7)。结果表明:矿区和缓冲区植被与气温呈显著性正相关，相关系数分别为0.443 0,0.412 1,0.501 0(p<0.05)，均通过显著性检验。校验区对气温响应敏感性较差。4区与降水呈不显著负相关:-0.094 9,-0.085 9,-0.152 8,-0.085 1，植被对气温的响应高于降水，这与以往研究结论有一致性[27-30]，气候因子对研究区植被做“超前-滑动增量”表明，降水对植被生长有2～3 a的滞后性(相关系数达-0.469,0.587)，气温对于植被生长无明显滞后性，秦超等对陕西省植被响应气候变化研究中也得出相同结论[31]。
图7 研究区AVHRR NDVI与降水(上)/气温(下)相关性曲线
Fig.7 Correlation curves of AVHRR NDVI and precipitation (above)/temperature (below) in study area
(1)确定NDVI生长期阈值的方法很多(如滑动平均值法[32-33]，NDVI阈值法[34-36]，温度界限法[37]等)，研究方法不同，得出结论也会有所偏差，笔者所采用的NDVI阈值法适合于我国北方植被生长期变化规律研究。由于植被类型不同，不同植被区域的NDVI阈值也有所不同。气候驱动因子只是植被活动的影响因素的一部分，且本文只研究了气候影响因子里最主要的降水和气温，并未研究蒸散量，湿度，日照强度等因素对植被的作用效果。研究结果表明，32 a来研究区植被活动显著增强，这与MYNENI,FANG,TUCKER,LOS等研究成果一致[4-6，8-9];针对RASIM，PROSPER，HUANG，徐占军、马超等对矿区植被研究表明[10-15]，植被在受到开采扰动后生长受到较大影响，生物量呈减少之势，而本文研究矿区植被NDVI累积量呈增加趋势，但低于自然环境校验区NDVI增长率，仍能说明煤炭开采对植被造成了一定影响。
(2)矿区植被生长期延长趋势高于校验区，一方面，说明植被在全球气候变化条件下，生长期有延长倾向;植被受气候变化影响返青期提前，生长期延长的特点，与农牧交错带、陕西及我国北方等地植被生长期研究结论一致[38-41]。而彬长矿区属粮矿复合区，工矿区工农业生产均会导致植被生长周期发生变化，有如下方面:① 工矿区人类生产生活实践，可能形成矿区热岛，改变局部水热条件[42-45]，使区域植被返青期提前，枯黄期滞后;② 矿区绿化、植树造林、退耕还林和土地复垦实践，都可能引入人工植被、改良生态结构[46]，使区域植被返青期提前，枯黄期滞后;③ 矿区周边为传统农业种植区，粮矿复合区的植被物候受农业生产影响较大，越冬作物或大秋作物都可能使区域植被返青期提前，枯黄期滞后;④ 煤矿开采造成土壤理化性质的变化导致植被生理周期被延长。如土壤氮、磷等养分元素损失，开采塌陷造成土壤表面裂缝生成，土壤持水保肥能力弱，植物根系受损，光合作用能力变差，不利于植物生长及植被恢复，植被生理周期被延长[47]。这一点可以从“矿区植被生长期延长，但NDVI增长率却下降”的结论得到一定的支持。
综上所述，矿区及其缓冲区的植被生长期由于人工干扰，不同程度地被二度延长，凸显了人类活动对生态环境的重大影响。
(1)32 a(1982—2013)彬长矿区月度NDVI变化表明，研究区年内植被变化呈双峰型，5—8月份植被指数较高，8月份达到峰值，9月底植被指数骤降，2月份植被指数达到谷值。年际上，4区平均NDVI总量为20 km缓冲区(358.33)、10 km缓冲区(237.42)、直接影响区(158.43)、生态校验区(121.93)，32 a间NDVI增长百分比为生态校验区(19.45%)>矿区(13.31%)>10 km缓冲区(10.9%)>20 km缓冲区(10.49%)，自然环境生长下的植被活动比煤矿开采区活跃，采煤活动对矿区植被造成了一定影响。
(2)1982—2013年研究区植被生长期研究表明:矿区及缓冲区在4月初进入返青期，9月底进入枯黄期，校验区进入返青期日期为4月上旬，枯黄期日期为10月上旬，四区生长期长度在185 d左右。时间序列拟合曲线表明，四区返青期提前，枯黄期推迟，生长期延长，矿区延长天数高于生态校验区5 d。
(3)植被响应气候因子相关性分析:研究区受全球气候变化影响，气温升高，降水量减少，处于半湿润地区的矿区与缓冲区植被对气温的响应更为敏感;校验区NDVI受其气温影响较小，植被与同期降水量相关性不大，但降水对其有明显滞后性。植被受气候变化影响返青期提前，生长期延长的特点，与农牧交错带、陕西及我国北方等地植被生长期研究结论一致。
参考文献:
[1] 李小静.彬长煤矿地表沉陷区植被变化遥感监测研究[D].西安:西安科技大学,2013:1-50.
LI Xiaojing.Research of vegetation change trend with remote sensing technology in coal mine surface subsidence area of Binchang[D].Xi’an:Xi’an University of Science and Technology,2013:1-50.
[2] MYNENI R B,KEELING C D,TUCKER C J,et al.Increased plant growth in the northern high latitudes from 1981-1991[J].Nature,1997,386:698-702.
[3] FANG J Y,PIAO S L,FIELD C,et al.Increasing net primary production in China from 1982 to 1999[J].Frontiers in Ecology and the Environment,2003(1):293-297.
[4] MYNENI R B,DONG J,TUCKER C J,et al.A large carbon sink in the woody biomass of Northern forests[J].Proc.Natl.Acad.Sci.USA,2001,98:14784-14789.
[5] LOS S O,COLLATZ G J,BOUNOUA L,et al.Global inter annual variations in sea surface temperature and land surface vegetation,air temperature,and precipitation[J].J.Climate,2001,14:1535-1549.
[6] TUCKER C J,SLAYBAEK D A,PINZON J E,et al.Higher northern latitude NDVI and growing season trends from 1982 to 1999[J].Int.J.Biometeorol,2001,45:184-190.
[7] ZHOU L M,TUCKER C J,KAUFMANN R K,et al.Variations in northern vegetation activity inferred from satellite data of Vegetation index during,1981 to 1999[J].J.Geophys.Res.,2001,106(D17):20069-20083.
[8] 方精云,朴世龙,贺金生,等.近20年来中国植被活动在增强[J].中国科学:C辑,2003,33(6):554-565.
FANG Jingyun,PIAO Shilong,HE Jinsheng,et al.Vegetation of China invigorated in last 20 years[J].Science in China:Series C,2003,33(6):554-565.
[9] NING Tingting,LIU Wenzhao,LIN Wen,et al.NDVI variation and its responses to climate change on the Northern Loess Plateau of China from 1998 to 2012[J].Advances in Meteorology,2015:1-10.
[10] RASIM Latifovic,DARREN Pouliot.Monitoring cumulative long-term vegetation changes over the athabasca oil sands region[J].Applied Earth Observations and Remote Sensing,2014,7(8):3380-3392.
[11] PROSPER LAARI Basommi,GUAN Qingfeng,CHENG Dandan.Exploring Land use and Land cover change in the mining areas of Wa East District,Ghana using Satellite Imagery[J].Meteorology & Atmospheric Sciences,2015(1):618-626.
[12] HUANG Yi,TIAN Feng,WANG Yunjia,et al.Effect of coal mining on vegetation disturbance and associated carbon loss[J].Environ.Earth Sci.,2015,73:2329-2342.
[13] 徐占军,侯湖平,张绍良,等.采矿活动和气候变化对煤矿区生态环境损失的影响[J].农业工程学报,2012,28(5):232-239.
XU Zhanjun,Hou Huping,ZHANG Shaoliang,et al.Effects of mining activity and climatic change on ecological losses in coal mining areas[J].Transactions of the Chinese Society of Agricultural Engineering,2012,28(5):232-239.
[14] 郝成元,杨志茹.基于MODIS数据的潞安矿NPP时空格局[J].煤炭学报,2011,36(11):1840-1844.
HAO Chengyuan,YANG Zhiru.Net primary production and its spatial-temporal in Lu’an mining area based on MODIS data[J].Journal of China Coal Society,2011,36(11):1840-1844.
[15] 马超,张晓克,郭增长,等.半干旱山区采矿扰动植被指数时空变化规律[J].环境科学研究,2013,26(7):750-758.
MA Chao,ZHANG Xiaoke,Guo Zengzhang,et al.Spatial-temporal variation of vegetation index caused by mining subsidence in semi-arid mountain regions[J].Research of Environmental Sciences,2013,26(7):750-758.
[16] 康萨如拉,牛建明,张庆,等.草原区矿产开发对景观格局和初级生产力的影响——以黑岱沟露天煤矿为例[J].生态学报,2014，34(11):2855-2867.
KANG Sarula,NIU Jianming,ZHANG Qing,et al.Impacts of mining on landscape pattern and primary productivity in the grassland of Inner Mongolia:A case study of Heidaigou open pit coal mining[J].Acta Ecologica Sinica,2014,34(11):2855-2867.
[17] JENSEN J R.遥感数字影像处理导论[M].陈晓玲,龚威,李平湘,等译.北京:机械工业出版社,2007:215-219.
[18] LASAPONARA R.On the use of principal component analysis(PCA) for evaluating inter annual vegetation anomalies from SPOT/VEGETATION NDVI temporal series[J].Ecological Modeling,2006,194(4):429-434.
[19] STOW D,DAESCHNER S,HOPE A,et al.Variability of the seasonally integrated normalized difference vegetation index across the north slope of Alaska in the,1990s[J].Remote Sensing,2003,24(5):1111-1117.
[20] DU PLESSIS W P.Linear regression relationships between NDVI,vegetation and rainfall in Etosha National Park,Namibia[J].Arid Environments,1999(42):235-260.
[21] 吴立新,马保东,刘善军.基于SPOT卫星NDVI数据的神东矿区植被覆盖动态变化分析[J].煤炭学报,2009,34(9):1217-1222.
WU Lixin,MA Baodong,LIU Shanjun.Analysis to vegetation coverage change in Shendong mining area with SPOT NDVI Data[J].Journal of China Coal Society,2009,34(9):1217-1222.
[22] DOUGLAS A Stow,ALLEN Hope,DAVID Verbyla,et al.Remote sensing of vegetation and land-cover change in Arctic Tundra Ecosystems[J].Remote Sensing of Environment,2004,89(3):281-308.
[23] PIAO Shilong,FANG Jingyun,ZHOU Liming,et al.Variations in satellite-derived phenology in China’s temperate vegetation[J].Global Change Biology,2006,12(4):672-685.
[24] 刘宪锋,任志远.西北地区植被覆盖变化及其与气候因子的关系[J].中国农业科学,2012,45(10):1954-1963.
LIU Xianfeng,REN Zhiyuan.Vegetation coverage change and its relationship with climate factors in Northwest China[J].Scientia Agricultura Sinica,2012,45(10):1954-1963.
[25] 马超,马雯思,王孜健,等.中国大陆1951—2012 年400 mm等降水量线的迁移及诱因[J].河南理工大学学报(自然科学版),2016,35(4):520-525.
MA Chao,MA Wensi,WANG Zijian,et al.Migration and its inducements of 400 mm precipitation contour in the mainland China from 1951 to 2012 year[J].Journal of Henan Polytechnic University(Natural Science),2016,35(4):520-525.
[26] 李琦,宋令勇,张文静,等.陕西省气温及日照时间变化特征分析[J].北京师范大学学报(自然科学版),2010,46(3):395-400.
LI Qi,SONG Lingyong,ZHANG Wenjing,et al.Tendency of air temperature and sunshine duration in Shaanxi Province[J].Journal of Beijing Normal University(Natural Science),2010,46(3):395-400.
[27] 任国玉,徐铭志,初子莹,等.近54年中国地面气温变化[J].气候与环境研究,2005,10(4):717-727.
REN Guoyu,XU Mingzhi,CHU Ziying,et al.Changes of Surface Air Temperature in China during,1951—2004[J].Climatic and Environmental Research,2005,10(4):717-727.
[28] 秦大河.气候变化科学与人类可持续发展[J].地理科学进展,2014,33(7):873-883.
QIN Dahe.Climate change science and sustainable development[J].Progress in Geography,2014,33(7):873-883.
[29] IPCC.Climate change,2001:The scientific basis[M].HOUGHTON J T,DING Y,GRIGGS D J,Eds.Cambridge:The Press Syndicate of Cambridge University,2001:365.
[30] 秦超.陕西省植被覆盖时空演变及其驱动因子研究[D].西安:陕西师范大学,2015:1-79.
[31] 李本纲,陶澍.AVHRR NDVI与气候因子的相关分析[J].生态学报,2000,20(5):898-902.
LI Bengang,TAO Shu.Correlation between AVHRR NDVI and climate factors[J].Acta Ecologica Sinica,2000,20(5):898-902.
[32] 罗隆诚,王俊.半干旱地区NDVI对气温和降水响应的多时间尺度分析——以甘肃省榆中县为例[J].干旱区资源与环境,2011(9):167-171.
LUO Longcheng,WANG Jun.The multi-time-scale response of NDVI to temperature and precipitation in semiarid areas-A case study for Yuzhong County of Gansu Province[J].Journal of Arid Land Resources and Environment,2011(9):167-171.
[33] 侯学会,牛铮,高帅,等.基于 SPOT-VGT NDVI 时间序列的农牧交错带植被物候监测[J].农业工程学报,2013,29(1):142-150.
HOU Xuehui,NIU Zheng,GAO Shuai,et al.Monitoring vegetation phenology in farming-pastoral zone using SPOT-VGT NDVI data[J].Transactions of the Chinese Society of Agricultural Engineering,2013,29(1):142-150.
[34] 王宏,李晓兵,韩瑞波,等.利用 NOAA NDVI和MSAVI遥感监测中国北方不同纬度带植被生长期变化[J].应用生态学报,2006,17(12):2236-2240.
WANG Hong,LI Xiaobing,HAN Ruibo,et al.Variability of vegetation growth season in different latitudinal zones of North China:A monitoring by NOAA NDVI and MSAVI[J].Chinese Journal of Applied Ecology,2006,17(12):2236-2240.
[35] REED B C,BROWN J F,VANDER ZEE D,et al.Measuring phonological variability from satellite imagery[J].Journal of Vegetation Science,1994,5:703-714.
[36] LI X B,CHEN Y H,FAN Y D,et al.Detecting inter-annual variations of vegetation growth based on satellite-sensed vegetation index data from 1983 to 1999[J].Proceedings of IGARSS,2003(5):3263-3265.
[37] FISCHER A.A model for the seasonal variations of vegetation indices in coarse resolution data and its inversion to extract crop parameters[J].Remote Sensing of Environment,1994,48:220-230.
[38] MARKON C J,FLEMING M D,BINNIAN E F.Characteristics of vegetation phonology over the Alaskan landscape using AVHRR time-series data[J].Polar Record,1995,31:179-190.
[39] JUSTICE C O,TOWN SHEND J R G,HOLBEN B N,et al.Analysis of the phonology of global vegetation using meteorological satellite data[J].International Journal of Remote Sensing,1985(6):1271-1318.
[40] 郭灵辉,吴绍洪,赵东升,等.近50 a内蒙古不同植被类型区生长期变化[J].干旱区地理,2014,37(3):532-538.
GUO Linghui,WU Shaohong,ZHAO Dongsheng,et al.Variation sand trends of climatic growing season in different vegetation zones,Inner Mongolia over the past 50 years[J].Arid Land Geography,2014,37(3):532-538.
[41] SKOUSEN Jeff,ZIPPER CARL E.Post-mining policies and practices in the Eastern USA coal region[J].International Journal of Coal Science & Technology,2014,1(2):135-151.
[42] ZHOU Jinhua,WANG Lejie.Comprehensive study on ecological restoration and land exploitation of mining subsidence in suburbs of Chinese mining cities[J].International Journal of Coal Science & Technology,2014,1(2):248-252.
[43] BI Yinli,ZOU Hui,ZHU Chenwei.Dynamic monitoring of soil bulk density and infiltration rate during coal mining in sandy land with different vegetation[J].International Journal of Coal Science & Technology,2014,1(2):198-206.
[44] XIAO Wu,HU Zhenqi,FU Yanhua.Zoning of land reclamation in coal mining area and new progresses for the past 10 years[J].International Journal of Coal Science & Technology,2014,1(2):177-183.
[45] 胡振琪,多玲花,王晓彤.采煤沉陷地夹层式充填复垦原理与方法[J].煤炭学报,2018,43(1):198-206.
HU Zhenqi,DUO Linghua,WANG Xiaotong.Pinciple and method of reclaiming subsidence land with inter-layers of filling materals[J].Journal of China Coal Society,2018,43(1):198-206.
[46] 卞正富,雷少刚,金丹,等.矿区土地修复的几个基本问题[J].煤炭学报,2018,43(1):190-197.
BIAN Zhengfu,LEI Shaogang,JIN Dan,et al.Several basic scientific issues related to mined land remediation[J].Journal of China Coal Society,2018,43(1):190-197.
[47] 胡振琪,龙精华,王新静.论煤矿区生态环境的自我修复、自然修复和人工修复[J].煤炭学报,2014,39(8):1751-1757.
HU Zhenqi,LONG Jinghua,WANG Xinjing.Selfhealing,natural restoration and artificial restoration of ecological environment for coal mining[J].Journal of China Coal Society,2014,39(8):1751-1757.