内蒙古荒漠草原植被对气候响应的滞后性研究

　　摘要：草原植被驱动机制研究一直备受关注。以1982-2015年内蒙古荒漠草原为研究对象，基于GIMMSNDVI3g及气象数据，采用皮尔逊相关系数、滞后系数探讨内蒙古地区荒漠草原NDVI对气候因子的滞后性。结果表明：1)荒漠草原NDVI、降水、气温和太阳总辐射年内变化都是单峰型。8月草原NDVI值达到最大，NDVI值年内变化不明显;降水也是8月达到最大;气温和太阳总辐射7月达到最大。2)从草原NDVI与气候因子的相关关系和显著相关面积占比来看，气候因子对草原NDVI影响从大到小的排列顺序为：降水>太阳辐射>气温。3)荒漠草原生长季植被的NDVI与气候的滞后性分析表明，5-9月NDVI与降水相关关系显著，滞后时间为1个月左右。4月NDVI与1-4月的降水不显著，降水不是草地返青的主导气候因子;6-9月草原NDVI与气温成负相关，滞后时间为0～1个月。4-5月草原NDVI与气温成正相关，该季气温的回升是植被返青的关键因子;6-9月草原NDVI与太阳辐射成负相关，滞后时间为1个月。5月草原NDVI与太阳辐射成正相关，滞后时间为1个月;4月NDVI与1-4月的太阳辐射关系不显著。

　　关键词：荒漠草原;NDVI;滞后时间;气候因子

　　草原是内陆半干旱到半湿润气候条件下特有的一种生态系统类型，不仅为人类提供生产、生活原料，还是支撑与维持地球生命的支持系统，维持生命的物质循环与水循环，维持生物物种与遗传多样性、净化环境、维持大气化学平衡与稳定[1-4]。

　　荒漠草原是草原区向荒漠区过渡的旱生化草原生态系统，在内蒙古约占草原总面积的10.7%[5-6]，该地区干旱少雨，属于典型的生态脆弱带，稳定性差，极易遭到破坏，对全球变化的响应十分的敏感[7-9]。草原植被驱动机制研究一直受到广泛关注[10]。归一化植被指数(NDVI)是目前公认的作为区域或全球植被覆盖度及植被生长状态的最佳指示因子[11-12]。研究结果表明，NDVI与气候因子有着显著的相关关系[13-15]。同时，植被对气候因子响应存在着明显的滞后性，植被生长除了受当时的气候因子影响外，前一段时间的气候条件的累积效应对植被的生长也产生很大的影响，并且不同地区植被对气象因子的响应时段不尽相同[16-19]。

　　辽宁省生长季植被从滞后期的时间长短看，植被生长对气温的变化最为敏感，其次是降水量，而对日照响应相对缓慢[20]。1982-2000年内蒙古地区植被NDVI与降水有很好的相关性，降水对植被的生长影响很大[21-22]。2001-2010年内蒙古草原植被覆盖度在年际水平上主要受降雨影响，月植被覆盖度对降雨量的响应存在时滞效应[23-24]。

　　目前关于气候对草地生长的影响研究，主要集中在降水方面。针对气温，尤其是太阳辐射对草地的生长研究比较少。同时研究的时间跨度比较小且缺少对空间分布的描述。基于以上原因，利用长时间序列(1982-2015年)GIMMSNDVI3g、ERA5降水、气温、太阳辐射数据，对内蒙古荒漠草原NDVI对降水、气温、太阳辐射的滞后性进行研究，探讨荒漠草原植被对不同气候因子的响应和时滞的差异性，以期揭示荒漠草原生长的规律特征及其气候驱动机制。

　　1研究区概况

　　内蒙古荒漠草原位于37°-46°N，105°-115°E，总面积约为11.2万km2，占全区草地总面积的10.68%。主要分布在阴山北麓和西鄂尔多斯，是内蒙古干旱区代表性的植被类型之一，面积辽阔，占据了乌兰察布高原大部及西鄂尔多斯。与荒漠带紧邻，在典型草原和草原化荒漠之间呈带状由东北向西南分布。该区为典型的温带大陆性气候，干旱少雨。

　　年均气温为4.79℃，最低气温在1月份，均温为-14.65℃，最高气温7月份，均温达22.32℃。平均年降水量为217mm，且年内分布不均匀，主要集中在6-8月。草原植被类型主要为戈壁针茅(Stipatianschanicavar.gobica)草原、石生针茅(Stipatianschanicavar.klemenzii)草原、沙生针茅(Stipacaucasicasubsp.glareosa)草原和短花针茅(Stipabreviflora)草原组成的草原群落。

　　2数据与研究方法

　　2.1数据获取及预处理

　　2.1.1NDVI数据

　　植被指数数据为美国宇航局(NASA)最新提供的第3代NDVI数据集(GIMMSNDVI3gV1.0)(https：//ecocast.arc.nasa.gov/data/pub/gimms/3g.v1/)，该数据集时间跨度为1982-2015年，时间分辨率为15d，空间分辨率为8km，经过几何精校正、辐射校正、大气校正、图像增强等预处理。使用Matlab对1982-2015年覆盖内蒙古草地的GIMMSNDVI3g数据进行格式转换、图像镶嵌、图像裁剪、Albers等面积投转换等处理，形成涵盖内蒙古草地的月度DNVI时空数据集。

　　2.1.2气象数据选用的气象再分析资料为EC-MWF公开的ERA5数据(http：//apps.ecmwf.int/datasets/)，ERA5是迄今最强大的全球气候监测数据集，与之前的ERA-Interim产品相比，ERA5具有更高的时空分辨率，时间分辨率为逐小时，空间分辨率高达0.125°。选用覆盖内蒙古草地1982-2015年1-12月0.125°的ECMWF-ERA5气象模式2m气温、降水、太阳总辐射数据。使用Matlab对数据进行格式转换、年度合成、图像裁剪、Albers等面积投转换、重采样等处理。

　　2.1.3荒漠草原数据荒漠草原类型数据是在中科院植物所2000年制作的1∶100万全国植被类型数据基础上，完成对内蒙古地区荒漠草原范围的提取。并对草原图进行投影转换、栅格化等处理，获得投影和分辨率与遥感数据一致数字化产品。

　　3结果与分析

　　3.1草地的NDVI年内生长变化规律

　　给出了1982-2015年荒漠草原类型月平均NDVI的变化趋势。荒漠草原ND-VI变化呈现单峰生长规律。草原植被从4月底到5月快速生长，到8月NDVI值达到最大，9月就开始衰退，到11月彻底枯黄。荒漠草原由于受自然条件和降水的影响，NDVI值一直变化不太明显(0.105～2.210)，但也表现出了一定季节性规律。

　　3.2气候因子年内变化规律

　　给出了1982-2015年荒漠草原月平均降水、气温、太阳总辐射的变化趋势，降水、气温、太阳总辐射年内变化都是单峰型。荒漠草原从1月开始降水量逐渐增加，到8月份达到最大，降水量为53.40mm，9-12月份开始逐渐下降;由中可知荒漠草原气温也是从1月份开始逐渐增加，到7月份气温达到最大值22.75℃，之后开始快速下降;月太阳总辐射量也是前7个 月逐渐递增，到7月份达到最大值860.18MJ/m2，之后开始快速下降。

　　3.3荒漠草原降水滞后时间及空间分布

　　荒漠草原9月份NDVI与1-9月的降水相关性最高的是8月份，其平均相关系数为0.433。由此可得，9月份NDVI与8月份降水有很好的相关性，9月份NDVI对应降水滞后时间为1个月。

　　其中，通过P<0.1显著性检验的面积占整个荒漠草原面积的63.30%，显著相关空间分布主要集中在锡林郭勒盟的苏尼特左旗、苏尼特右旗，乌兰察布盟的四子王旗，包头，巴彦淖尔的乌拉特中旗、乌拉特后旗，鄂尔多斯的杭锦旗、鄂托克旗等地区;8月份NDVI与1-8月的降水相关性最高的月份是7月份，其平均相关系数为0.480。8月份的滞后时间为也是1个月。其中，通过P<0.1显著相关的区域占整个草原面积的81.47%，除了阿拉善左旗、鄂托克前旗，显著相关地区几乎覆盖了整个荒漠草原;7月份NDVI与1-7月的降水相关性最高的月份是6月份，其平均相关系数为0.513。7月份的滞后时间也是1个月左右。

　　其中，通过P<0.1显著相关的区域占整个草原面积的93.00%。相关系数和面积占比是所有月份里面最高的，几乎覆盖荒漠草原整个地区;6月份NDVI与1-6月的降水相关性最高的月份是5月份，其平均相关系数是0.486。6月份的滞后时间为1个月。其中，通过P<0.1显著相关的面积占整个草原面积的88.86%，显著相关的区域除了乌拉特前旗和部分乌拉特后旗，覆盖了剩下整个荒漠草原;5月份NDVI与1-5月的降水相关性最高的月份是4月份，其平均相关系数是0.408，5月份的滞后时间也是1个月。显著相关的区域主要分布在锡林郭勒盟的苏尼特左旗、苏尼特右旗，乌兰察布盟的四子王旗，鄂尔多斯的杭锦旗、鄂托克旗等地区;4月份NDVI与1-4月的降水显著性不明显，且1-3月平均相关系数为负，4月份正好是草地返青季，由此可见降水并不是草地返青的主导气候因子，导致4月份草地的生长可能还受其他气候因子控制。

　　4结论与讨论

　　4.1讨论

　　荒漠草原NDVI对气候因子的响应存在一定的滞后性，且不同学者研究结论差异比较大[25]。张颖等[26]研究发现三江源地区荒漠草原受降水影响较大，在月尺度上草地覆盖度对气温表现出明显的滞后效应，而降水滞后效应不明显;李晓兵等[27]认为内蒙古荒漠草原分布区夏季的降水对当年植被生长的影响最明显;王举凤等[18]认为内蒙古荒漠草原植被覆盖度对降水量存在明显的1～2个月滞后和累积效应，且表现出越干旱的草原类型滞后效应越明显的特征，荒漠草原植被对降水量变化更加敏感。

　　王雨晴等[23]认为内蒙古荒漠草原的月NDVI对降水、温度存在1～2个月滞后效应，且降水对NDVI的影响大于温度;刘成林等[28]通过气象站点资料收集方式研究发现二连浩特、苏尼特左旗、朱日和荒漠草原NDVI对降水滞后时间都为64d以上，那仁宝力格和多伦滞后48d;李霞等[29]也是通过气象站数据研究荒漠草原气候因子时滞效应，结果发现ND-VI对降水滞后时间1个月，气温的滞后时间不明显。

　　草原NDVI与当月的气温相关关系最强，且气温在植被生长季前期与NDVI成正相关，中后期成负相关关系，这个结论与本文的研究基本一致。由此可见，目前绝大多数研究主要集中在水热滞后性方面，而本研究发现太阳辐射对荒漠草原植被ND-VI影响要大于气温的影响，且与气温对NDVI响应比较类似，也是在植被生长季前期与NDVI成正相关，中后期成负相关关系，滞后时间为1个月。ND-VI与降水的相关最强，滞后也是1个月。NDVI与气温的相关关系在三者中最弱，滞后约0～1个月。总的来看，以上结论有的和本研究一致，有的不同，主要原因可能是研究区域不一致，也可能是数据、研究方法不同造成的。

　　4.2结论

　　利用1982-2015年的GIMMS月度NDVI、ECMWF-ERA5月度降水、气温、太阳辐射数据，对内蒙古地区荒漠草原NDVI的气候因子滞后性进行研究。荒漠草原NDVI、降水、气温、太阳总辐射年内变化都是单峰型。草原到8月NDVI值达到最大，NDVI值(0.105～2.210)年内变化不太明显;降水8月份达到最大，降水量为53.40mm;气温7月份达到最大值22.75℃;太阳总辐射也是7月份达到最大值860.18MJ/m2。

　　生态论文投稿刊物：应用生态学报主要报道应用生态学领域的创新性科研成果与科研进展，反映我国应用生态学的学术水平和发展方向，跟踪学科发展前沿，注重理论与应用结合，促进国内外学术交流合作与人才培养。

　　荒漠草原生长季植被的NDVI与降水的滞后性分析表明，5-9月NDVI与降水相关关系比较显著，滞后时间为1个月，显著相关的区域主要分布在锡林郭勒盟的苏尼特左旗、苏尼特右旗，乌兰察布盟的四子王旗，鄂尔多斯的杭锦旗、鄂托克旗等地区。4月份NDVI与1-4月的降水不显著，降水不是草地返青的主导气候因子。

　　参考文献：

　　[1]王伯荪，彭少麟.植被生态学群落与生态系统[M].北京：中国环境科学出版社，1997.

　　[2]ZHAOTQ，OUYANGZY，JIALQ，etal.Ecosystemserv-icesandtheirvaluationofChinagrassland[J].ActaEcologicaSinica，2004，24(6)：1101-1110.

　　[3]ADAMSJM，FAUREL，FAURE-DENARDJM，etal.Wood-wardincreasesinterrestrialcarbonstoragefromthelastglacialmaximumtothepresent[J].Nature，1990，348(6303)：711-714.

　　[4]CURLOCKJM，HALLDO.Theglobalcarbonsink：agrass-landperspective[J].GlobalChangeBiology，2010，4(2)：229-233.

　　[5]吴永胜，马万里，李浩，等.内蒙古退化荒漠草原土壤有机碳和微生物生物量碳含量的季节变化[J].应用生态学报，2010，21(2)：312-316.WUYS，MAWL，LIH，etal.Seasonalvariationsofsoilor-ganiccarbonandmicrobialbiomasscarbonindegradeddesertsteppesofInnerMongolia[J].ChineseJournalofAppliedE-cology，2010，21(2)：312-316.(inChinese)

　　作者：滑永春，高润宏*，敖敦