第35卷第10期2019年10月赤峰学院学报(自然科学版)JournalofChifengUniversity(NaturalScienceEdition)Vol.35No.10
Oct.2019
内蒙古四子王旗干旱变化特征及与植被覆盖关系研究韩美娟1袁金额尔德木吐1,2,3袁郭恩亮1,2袁苏志刚1渊1.内蒙古师范大学地理科学学院曰2.
内蒙古自治区蒙古高原灾害与生态安全重点实验室曰
呼和浩特
010022冤
3.内蒙古自治区遥感与地理信息系统重点实验室袁内蒙古
摘
要院本文基于1961年-2018年四子王旗地面降水和气温日值0.5毅伊0.5毅格点数据集(V2.0)袁计算该
区域23个网格点标准降水蒸散指数渊SPEI冤袁并且运用森式斜率法分析了研究区干旱空间分布情况袁运用2000年-2018年ModisNDVI数据袁利用像元二分法模型计算了植被覆盖度袁分析研究区生长季植被覆盖度的变化情况.最后利用相关性分析方法进一步研究干旱和植被覆盖之间的关系.研究表明院渊1冤四子王旗地区1961年-2018年整体上SPEI为下降趋势袁但2000年-2018年在SPEI年际变化上呈现上升趋势袁干旱情况有所改善.渊2冤2000年-2018年生长季植被覆盖度整体上呈现上升趋势.渊3冤2000年-2018年来SPEI与DNVI呈现较好的正相关关系袁63.9%区域SPEI与NDVI呈正相关关系.本文的研究结果表明干旱对植被覆盖情况有较大影响.
关键词院标准降水蒸散指数(SPEI)曰归一化植被指数渊NDVI冤曰像元二分法曰相关性曰四子王旗中图分类号:Q948
文献标识码:A
文章编号:1673-260X(2019)10-0118-05
间均趋于干旱情况[6].徐建文尧局辉等发现黄淮海平
原在春夏两个季节中随着温度的升高有趋向于干旱化的现象袁同时该地区蒸散以及降水由于在空间上的不平衡性袁导致了该地区的干旱特征出现了明显的南北差异.严晓瑜学者研究若尔盖湿地植被年际的生长状况与降水尧气温等气候因子有明显的相关性袁其中影响植被最重要的因素就是气温袁同时植被生长季的NDVI与气候因子又存在一定的空间差异等[8].包刚尧包玉海等建立了植被覆盖度与高光谱反射率数据的三种转换方式的关系模型袁对不同转换方式在植被覆盖度估算中的应用优势做了对比分析袁得出植被冠层反射率数据与植被覆盖度呈现负相关的关系[9].
地处中温带性季风气候区袁春季多风干旱袁夏季炎热短暂袁秋季相对凉爽袁年蒸发量是降水量近8.1倍袁是导致作物袁牧草受旱的重要原因之一.地势高寒袁昼夜温差较大袁风沙不断袁寒暑变化比较剧烈.因此本文通过利用NDVI数据与SPEI指标袁来研究我国内蒙古自治区四子王旗干旱变化特征与植被覆盖的关系.2研究区概况
四子王旗隶属我国内蒙古自治区乌兰察布市袁地理位置在大青山北麓袁乌兰察布市西北部袁北接
1前言
近年来干旱已经成为我们国家甚至整个地球都面对的自然灾害问题袁它所涉及的范围广袁时间也相对比较长袁频率也是越来越高袁导致负面影响越来越大.在社会经济尧生态系统尧气候等方面都造成了一定的影响.而对生态系统的影响更为显著袁植被的减少是一部分是因为干旱导致植被减少可以在一定程度上造成土地退化袁生物的多样性减少袁水土流失等袁进一步影响生态环境恶化[1].因此袁对于研究一个地区的干旱变化特征与植被覆盖关系具有很重要的意义.
近年一部分学者在某些区域或者全球范围内袁在干旱变化特征及植被覆盖方面做出了很多研究袁如国外学者利用SPEI分析了1982年要2006年地中海和1961年要2012年捷克地区的干旱变化时空特征[3-4].胡子瑛尧周俊菊等学者关于中国北方气候干湿变化及干旱演变特征方面得出1960s是最为湿润年代袁2000s是最为干旱年代.不同的干湿区干湿变化趋势也是差异显著袁不同干湿区在不同季节干湿变化也是差别较大等结论[5].李虹雨学者在
关于内蒙古地区气候因子变化及干旱特征研究中得出内蒙古地区除西部年际尧春尧秋尧冬四季袁东部春尧冬两季尧中部春尧秋两季气候趋湿润外袁其他时
收稿日期院2019-07-19-118-
蒙古国袁东部与乌兰察布市察右中旗尧察右后旗及锡林郭勒盟苏尼特右旗毗邻袁南与乌兰察布市卓资县尧西连包头市达尔汗茂名安联合旗袁地理坐标为北纬41毅13忆耀43毅24忆袁东经110毅23忆耀113毅03袁总面积241.33万hm2[10]袁国境104.12公里渊图1冤.
2000年要2018年研究区植被覆盖度遥感图像袁共228幅.运用ENVI5.0对遥感影像进行裁剪处理后袁通过IDL编程进行植被覆盖度计算并进行最大合成等处理袁最后研究分析四子王旗地区干旱变化特征及与植被覆盖关系.
3.2研究方法
3.2.1标准化降水蒸散指数SPEI
四子王旗地区23个网格点1961年要2018年地面降水日值0.5毅伊0.5毅格点数据集(V2.0)和地面气温日值0.5毅伊0.5毅格点数据集(V2.0)来源于中国气象
图1研究区位置示意图
王旗研究区及周边共23个网格点地面降水日值0.5毅伊0.5毅格点数据集(V2.0)和地面气温日值0.5毅伊0.5毅格点数据集(V2.0).NDVI数据来源于美国国家航天航空局NASA的MODIS数据产品MOD13A1袁空间分辨率为500m伊500m袁时间分辨率为16日.选取
渊http://data.cma.cn/冤,运用了1961年要2018年四子
3数据来源及研究方法3.1数据来源
本文的SPEI数据来源于中国气象数据网
表1
SPEI干旱情况
臆-2.0极端干旱
-2.0约SPEI臆-1.5严重干旱
科学数据共享服务网.标准化降水蒸散指数SPEI是指标准化潜在蒸散与降水的差值袁它用来衡量一个地区的干湿度与该地区平均状态的偏离程度[2].其中SPEI干旱等级分类标准[11]如下表1:,
-1.5约SPEI臆-1.0中度干旱
-1.0约SPEI臆-0.5轻度干旱
-0.5约SPEI臆0.5正常
SPEI干旱等级划分表格
并且袁为了更好地了解SPEI的性能袁采用了滑动累计方法统计了包括1尧3尧6和12月尺度的SPEI的年变化趋势袁以及春夏秋冬四个季节的SPEI变化趋势袁利用以上表格划分标准作2000年要2018年年际SPEI干旱等级空间分布图.3.2.2像元二分模型
像元二分模型近来常用于表征植被覆盖度袁它又是混合像元分解模型中较为简单的模型袁模型中输入的参数相对容易获取并且计算出的结果可靠袁所以也是它被广泛运用的有利因素.植被覆盖度其实是指一些草尧灌尧乔和农作物等植被的枝叶尧冠层在一定区域的垂直投影面积占该区域总面积的百分比.假设像元由全裸土和全植被组成袁其中S1为纯植被所覆盖的像元袁S2为纯裸土所覆盖的像元袁S为研究区域的整体覆盖的像元[12].则植被覆盖度F1的计算公式如下院
F1=渊S-S2冤/渊S1-S2冤渊6冤
由于植被覆盖度与NDVI有显著的正相关关系袁因此可以将植被覆盖度与像元二分模型相结合估算出植被覆盖度袁式子如下院
F1=渊NDVI-NDVImin冤/渊NDVImax-NDVImin冤渊7冤
其中院区域内NDVI的最大值为NDVImax袁ND鄄
的噪声影响袁根据图像的实际情况来决定NDVI鄄
max和NDVImin取置信度范围以内的最大值和最小值[13].提取研究区统计结果中NDVI值得积累率NDVImax为95%袁NDVImin为5%.3.2.3Mann-Kendall突变分析法
本文利用M-K渊Mann-Kendall冤突变分析法来分析研究四子王旗地区干湿变化的突变点.M-K趋势分析法相对于线性回归方法来说袁它非参数检验袁不需要方差正态性假设袁可以检验时间序列的变化趋势袁同时还允许缺值等特点[11].UF尧UB是M-K统计量袁它们的取值范围是负无穷到正无穷的袁并且当UF掖0时袁说明这时的时间序列为下降趋势袁反之为上升趋势袁当UF超于临界值袁就说明趋势比较明显.当UF尧UB它们交于临界值范围之内时袁该交点被称为突变点袁此时为突变时间开始点[2].3.2.4
VImin为区域内NDVI的最小值.由于存在不可避免
森式斜率估计法
Sen式斜率估计法表示在一个线性趋势单位时间内趋势减少或者增加的量[14].本文采用Sen式斜率计算变化趋势袁避免了数据分布形态对分析结
果以及时间序列缺失的影响袁同时也剔除了一些异常值对时间序列的影响.计算公式如下院
-119-
Q=median(Xa-Xb)1式中袁表示为Sena-b式斜率袁Xa和Xb分别为a
和3.2.5b时间通相序列值[15].常关研究系数法
分析两个变量之间的相互关系通过相关分析法袁这种方法早被广泛利用.通过计算四子王旗2000年要2018年年际SPEI值与NDVI之间的相关系数R袁来直观有效的说明两者之间的密切联系[16]袁其公式如n
下院
R=
式中袁姨n表移n移i=1i=1
示(xi-x1)(xi)-x1(yi-y1)2
年份数19姨移ni=1
渊2000-2018(yi-y1)2
渊9冤
冤袁xi和yi表示相关分析的两个变量曰xl和yl分别表示变量样本值的平均值[16]4.4.1结果由内蒙古四子王旗干旱变与分析
图2a我们可以看到袁近化58特征
年来四子王旗地区年SPEI指数以0.105伊10a-1的速度减少袁整体呈现1965出干旱他年尧1966的趋势年.1961尧1997年年出要2000现了年中近度干旱40年现内象只袁在2014大多年份都属于正常范围内的袁但是在2000其状态袁年在份内2005袁除年了达2003到了年极袁端其他年干旱袁份2001都处年于尧干旱年-2007年2015尧2009趋势图年年干旱也是比较严重袁达到了重度干旱.到可要看2018出变年化有波很动比大的较改大观.恢图复三了渊b正冤M-K常袁但趋由势分UB析地区值图检年袁均上验下结SPEI两果虚袁指线红数为色波临实动界线较线为为.频分UF繁析值袁显袁蓝色虚线为由示袁UF四子王旗曲线可见袁从20世纪70年代中期以后整体呈现下降趋势袁以2006年要2015年下降趋势最为明显渊超过临界线冤袁在0.05显著水平下袁19年为突变时间的开始.
图2四子王旗1961年-2018年年均SPEI指数年际变化和M-K突变检验
为了更好地说明四子王旗地区干旱的变化特征袁绘制了图2中1月尧3月尧6月袁12月不同尺度的随时间变化的SPEI值图.通过对四幅图的比较袁四个时间尺度的SPEI均呈现下降趋势.但是在不同时间段SPEI值也是有所不同的袁在1961年要
-120-
2000范围袁年但袁是在不同尺2000度下年的要SPEI20多年期时间间袁四都个处于时间正常
尺度SPEI波动较大袁出现不同程度的干旱现象.
图3四子王旗多时间尺度标准降水蒸散指数1961年要2018年动态特征
图3四幅图以1961年要2018年四子王旗地区3月尺度的SPEI指数来表示干旱的季节特征袁可SPEI以清晰地看到袁研究区在春尧夏尧秋三个季而冬季指的数呈SPEI现下指降数趋是势呈袁现表上明升趋干旱势程袁度也有就所是增节说加的近袁年来冬季呈现湿润化趋势.在春季和秋季SPEI变化的趋势不是很明显袁冬夏两季节的变化相对明显些袁秋季到冬季明显的表现由干旱转湿润袁而整体上看夏季加剧了干旱的程度.
图4四子王旗地区1961年要2018年SPEI指数四季变化趋势
通过四子王旗2000年要2018年SPEI年际变化统计图袁可以得出在2000年要2018年SPEI年际变化线性趋势整体上呈现为上升袁说明一些区域的干旱情况有了一定程度的改善.通过显著性95%的四子王旗干旱sen斜率变化趋势分析图袁彩色区域为通过显著性95%的区域袁在这部分区域内袁偏绿色部分为干旱减少变化显著的区域袁偏红色部分为干旱减少变化不显著的区域袁这两幅图反映的情况具有一致性袁整体上干湿度有所改善恢复.
图5
四子王旗2000年要2018年SPEI年际变化统计图与通过显著性95%的四子王旗干旱森式斜率4.2
内蒙古四子王旗植变化趋势分析图
被覆盖变化特征
2000利用了四子王旗及周边23个网格点提供的228长季幅年月袁-2018份利渊用文Arcgis年每个月中所取的4-10像份元的月统计数植被覆盖遥感图像共份冤的据植将被覆盖每年植度物生作最大合成袁根据19年的遥感数据计算的19幅植被覆盖度分类结果图袁本文选取其中六个年份的图袁进行象征性分析19年的变化趋势.从图6反应的情况看袁在2000年要2018年植被覆盖情况整体上是恢复的状态袁植被覆盖度处于上升趋势.从不同年份的变化趋势上看袁从2000年到2001年袁植被覆盖度出现下降趋势袁低植被覆盖度区域大面积增加2003袁2002年到2004年植被覆盖情况有所增加袁度整体年上的又变出化现情了况下最降为的明趋显.势在袁在20052006年年植到被覆盖2013年每年的植被覆盖情况都发生不同程度的上升现象袁直到2014年植被覆盖再次出现明显的下降袁低植被覆盖度和较低植被覆盖度的区域大面积增加.基于20182014年的植被覆盖情况袁从盖度和年高植植被覆盖被覆盖上度升区速域度相大面对积较增大2015加袁.
较年高起植直被覆到图6四子王旗植被覆盖度分类图
4.3内蒙古四子王旗干旱及植被覆盖关系研究
2000利用四子王旗及周边23个网格点提供的并做出年了如要2018图7年的的植年被覆盖均SPEI度和与NDVI干旱的数相据关袁计算性分布图袁统计发现研究区63.9%区域SPEI与NDVI呈正相关关系袁表明整体变化上看大部分干旱程度
减弱的情况下袁植被覆盖度也有一定程度的增加.同时在其他外界等各种原因的影响下袁也有27.6%的区域SPEI与NDVI呈现负相关关系袁说明研究区19年以来SPEI与NDVI存在较好的正相关.再利用以上通过显著性95%的四子王旗干旱与植被覆盖度森式斜率变化趋势分析结果袁做出如图7所示的植被覆盖度与干旱的显著性图袁同以上分析相似95%.图显著的片的区彩区域色域袁区袁偏域偏绿色为植被覆盖度与干旱通过显著性红色部为分为变化植相被覆盖对不显著度与的干旱变区域袁化整体上看研究区大面积上出现不同趋势程度的变化状态.
图7四子王旗植被覆盖度与干旱的相关性渊a冤和植被5结论覆盖与展望
度与干旱的显著性渊b冤
利用四子王旗23个网格点提供的气象数据袁分析四子王旗干旱变化特征及与植被覆盖关系影响的研究SPEI在SPEI为上渊1冤结季升趋四子王旗果如下院
节变势化袁上在不地区袁春同1961尧夏尺年要2018年整体上尧度秋季上也呈是为下降下趋降势趋袁冬季势袁为上升趋势.但2000年要2018年在SPEI年际变化上呈现上呈现渊上2上冤升趋2000势升趋年袁势要干旱情.相2018况对于年有2000生所长改善.
年季袁植到被覆盖2018年度整低体及较低植被覆盖度的区域较少2.5%袁较高及高植被覆盖度区域增加了2.2%.
好的渊正3冤相2000关关年系要袁63.9%2018年区来域SPEISPEI与与DNVINDVI呈呈现正较相关关系袁说明研究区63.9%的区域干旱程度在减弱
-121-
的情况下袁植被在生长季的覆盖度也有一定的提高.
干旱变化特征与植被覆盖之间的关系还存在其他元素的影响袁同时干旱变化特征与植被覆盖的研究是一个长时间系列研究.本文对四子王旗干旱变化特征与植被覆盖关系研究的时间序列较短袁只用23个网格点所获取的数据还不够全面袁在今后的研究当中需要进一步考虑到数据的精确性.文中对植被覆盖度的分类还不够细袁使用的标准化降水指数以降水量为基础袁所以分析出的干旱变化特征未考虑到气温尧土壤尧以及连续无降水日等方面对干旱变化的影响.对于以上存在问题袁在今后的研究中要要要
参要需考要要进文献要要一要要步完善袁提高准确性.院
要要要要要要要要要也1页孔冬冬,张强,顾西辉袁王月袁李华贞.植被对不同尺度干旱事件的响应特征及成因分析[J].生态学报袁2016,36(24):7908-7918.
也2页万红莲袁王静.多尺度下宝鸡地区干旱动态格局演变及其与植被覆盖的关系[J].生态学报袁2018袁38渊19冤院6941-6950.也3页POTOPOVAV,PETRSTEPANEKC
P,
MOZNYM,etal.Performanceofthestandard鄄izedprecipitationevapotranspirationindexatvariouslagsforagriculturaldroughtriskassess鄄mentintheCzechRepublic[J].AgriculturalandForestMeteorology,2015,202(202):26-38.也4页GOUVEIAC,TRIGORM,BEGUERIAS,etal.DroughtimpactsonvegetationactivityintheMediterraneanregion:anassessmentusingremotesensingdataandmulti-scaledroughtindicatore[J].GlobalandPlanetaryChange,2016,151(6):15-27.
也5页胡子瑛袁周俊菊袁张利利袁魏伟袁曹建军.中国北
方气候干湿变化及干旱演变特征.生态学报,2018,38(6):1908-1919.
也6页李虹雨.内蒙古地区气候因子变化及干旱特征研究[D].内蒙古院内蒙古农业大学袁2017.
也7页徐建文袁局辉袁刘勤袁杨建莹.黄淮海地区干旱变化特征及其对气候变化的响应[J].生态学报袁2014袁34渊2冤院460-470.
也8页严晓瑜.不同时间尺度若尔盖湿地植被变化及
其与气候的关系[D].北京院中国气象科学研究
-122-
院,2008.
也9页包刚袁包玉海袁覃志豪袁周义袁黄明祥袁张宏斌.高
光谱植被覆盖度遥感估算研究[J].自然资源学报袁2013袁28渊7冤.
也10页郭强袁米福贵袁殷国梅袁张富贵.气象因子对内
蒙古四子王旗草原退化的影响[J].畜牧与饲料科学.2008袁29渊6冤院28-30.
也11页马志婷袁武志涛袁卫洁.京津风沙源区干旱时空
特征及对植被变化影响[J].山地学报,2014,34(4):536-546.
也12页季浩宇袁用环境小卫星遥感数据分析攀枝花市苏保区植被变化情况[J].四川环境袁2016袁35(4冤.也13页佟斯琴袁包玉海袁张巧凤袁都瓦拉.基于像元二
分法和强度分析方法的内蒙古植被覆盖度时空变化规律分析[J].生态环境学报袁2016袁25渊5冤院737-743.也14页陈淑芳.湘江径流泥沙变异规律及其对气候变
化与人类活动的响应研究[D].湖南师范大学,2016.
也15页李运刚,胡金明,何大明,柳江.1960-2007年红
河流域强降水事件频次和强度变化及其影响[J].地理研究,2013,32(01):-72.
也16页李跃鹏袁刘海艳袁周维博.陕西省1982-2015NDVI时空分布特征与气候因子相关性[J].生态科学袁2017袁36渊6冤院153-160.
也17页曲学斌袁孙小龙袁冯建英袁范雪松.呼伦贝尔草
原NDVI时空变化及其对气候变化的响应[J].干旱气象,2018,36(1),97-103.
也18页王宇航袁赵鸣飞,康慕谊袁左婉怡.黄土高原地区
NDVI与气候因子空间尺度依存性及非平稳性研究[J].地理研究袁2016袁35渊3冤.
也19页刘世梁袁田韫钰袁尹艺洁袁安南南袁董世魁.云南省植被NDVI时间变化特征及其对干旱的响应[J].生态学报袁2016袁36渊15冤院4699-4707渊7冤.也20页赵平伟袁李宏波袁郭萍袁杨文.滇西南地区NDVI
变化及其对不同时间尺度干湿变化的响应[J].山地学报袁2018袁36渊2冤院229-238.也21页何鑫袁吴吉东袁李颖袁冯玉.基于SPEI的辽东地
区气象干旱时空分布特征[J].干旱区地理袁2017袁40渊2冤.
也22页邓少福.祁连山气候变化对植被的影响研究
[D].甘肃院兰州大学袁2013.