基于元胞自动机民勤绿洲湖区荒漠化演化预测

文章编号:1000 694X (2009)05 0802 06

基于元胞自动机民勤绿洲湖区荒漠化演化预测

收稿日期:2008 12 10;改回日期:2009 01 16

基金项目:国家自然科学基金项目(40501073,40501068)资助

作者简介:宋冬梅(1973 ),女(汉族),博士,副教授,硕士生导师,现从事GIS 与遥感的应用研究。E mail:songdongmei1973@

宋冬梅1

,吴远龙1

,张志诚2

,王建3

,李小玉

4

(1.中国石油大学地球资源与信息学院,山东青岛266555;2.中国动物卫生与流行病学中心,山东青岛266032;3.中国科学院寒区旱区环境与工程研究所,甘肃兰州730000; 4.中国科学院沈阳应用生态研究所,辽宁沈阳110016)

摘 要:民勤湖区是民勤绿洲中生态环境最为恶劣的地区,土地荒漠化问题十分突出。以民勤绿洲湖区为例,解译1992年、1998年、2002年和2006年T M 卫星影像,分析其荒漠化动态变化情况,利用A r cO bjects 模块结合地理元胞自动机理论构造荒漠化动态模拟模型,通过对比2006年的预测数据与实际数据,对模型进行参数调整和预测检验。预测结果表明,模型预测的准确性达到90%。最后对2012年该区土地利用状况做出预测,进而对荒漠化的发展趋势进行预测分析。

关键词:荒漠化;元胞自动机;蒙特卡罗方法;A rcObjects 模块;民勤绿洲中图分类号:X144

文献标识码:A

民勤县地处河西走廊东北部,位于石羊河流域下游,现有人口30 7万人,东西长203km,南北宽150km ,总面积1 6万km 2

,平均海拔为1350m [1]

,东、北、西三面被腾格里沙漠和巴丹吉林沙漠包围,它镶嵌在这两大沙漠之中,阻止两大沙漠的汇合,在西北地理变化和稳定中起着重要的作用。自20世纪60、70年代,随着石羊河流域土地大规模开发,该地区地下水水位的急剧下降[2],水资源匮乏,土地沙化严重,已成为中国沙尘暴的4大策源地之一[3-5],而民勤湖区是民勤绿洲中生态环境最为恶劣的区域,地下水矿化度持续升高,地下水位的下降导致了湖区大量弃耕地及生态灾民的产生,北部沙漠的威胁更加剧了该区环境恶化程度[6-7]。

近几年来,遥感(RS)、地理信息系统(GIS)和全球定位系统(GPS)及其集成技术迅速发展,大量应用于灾害监测、资源监测等方面,在多元数据定量分析与综合研究方面取得了很好的效果,但动态模拟演化研究还显不足[8]。笔者基于地理元胞自动机理论(Geo CA)将此地区荒漠化动态演化规律与其空间分布特征相结合,利用ArcObjects 开发技术,结合数学模型探讨该区土地荒漠化时空动态演变规律,并预测其未来发展趋势。

1 研究区域概况

研究区选择民勤绿洲生态环境最为恶劣的湖区(38 42!45∀ 39 05!51∀N,103 24!21∀ 103 48!54∀

E)。该区位于民勤绿洲的东北部,包括西渠乡、中渠乡、东湖镇和收成乡,绝大部分粮食生产受到地下水质的影响与制约[9]。

研究区属典型的温带大陆性干旱气候,20世纪50年代以来多年平均降水量110m m 左右,蒸发量2644mm ,风沙大,盛行西北风,年平均风速2 8m #s -1,风沙及沙尘暴频繁,年均风沙日数达139d,最多时达150d;8级以上大风日数70多天,强沙尘暴日数29d 。灌区内作物以春小麦、夏玉米、青稞等粮食作物为主,兼种棉花、茴香、籽瓜等经济作物,土壤类型为灌淤土。

2 研究方法

2.1 Geo CA 模型模拟民勤绿洲湖区荒漠化可行

性分析 CA 模型是离散动态系统概念和应用建模的一种方法。标准的CA 系统是一个由元胞(cells)、元胞状态(Cellular States)、邻域(neighbors)和规则(T ransitio n Rules/Functions)构成的四元组[10-11]。从数学模型的角度看,CA 模型是基于元胞的动态模拟系统,它能够与基于栅格数据的GIS 组件很好地集成,这样,一方面增强了当前GIS 软件所缺乏的动态建模能力,提高了操作性能,并为处理时态维提供了一个很好的方法,另一方面,GIS 强大的空间数据处理能力也可为CA 模型准备数据和定义有效

第29卷 第5期

2009年9月

中 国 沙 漠

JO U RN A L O F DESERT R ESEA RCH

V ol.29 No.5S ep.2009

的转换规则,还可以对结果进行直观的显示。因此, CA与GIS组件的集成,可以克服各自的缺点,形成一个全新的优势互补的动态系统,用来对复杂时空现象、行为和过程进行动态建模分析[12]。

荒漠化指包括气候变异和人类活动在内的种种因素造成的干旱、半干旱、亚湿润干旱区的土地退化,它是土地景观变化的一种,也是一种行为很复杂的地学现象[13-14]。通常景观变化的模拟模型不仅要了解一种景观现状变化到另一种景观现状,更重要的是要清楚景观发生变化的原因。

使用ArcObjects模块建立模型,在确定模型中的转换规则时,将能够量化采集的各种因素考虑进去,充分体现荒漠化的变化原因。民勤绿洲荒漠化变化集中表现为原有耕地和新开荒的弃耕地以及有林地、沙化林地、沙化草地、盐渍化土地被无序的开荒而转变为沙化耕地以及沙地[1]。运用CA方法来研究这种变化过程,不仅能兼顾荒漠化的空间维,还能顾及其时间维,为荒漠化动态模拟模型的研究提供了一个很好的方向。

2.2 Geo CA中土地类型参数的转换规则

本研究主要针对民勤绿洲荒漠化问题,即耕地与完全荒漠化土地之间的动态演变关系,因此,为便于模型的应用,将土地类型分为绿地,沙化地及沙地,其中绿地包括耕地、草地以及林地,沙化地包括沙化林地,沙化草地和沙化耕地。其中沙地指完全荒漠化的土地,植被覆盖度很低,沙化地指介于绿地和沙地之间具有一定植被覆盖的土地;盐碱地依据地表植被覆盖程度的大小,分别归入沙化地和沙地。归一化植被指数(NDVI)能够表征地表植被的覆盖状况,在一定范围内,NDVI值与地表植被覆盖度呈现正相关,据此,将归一化植被指数(NDVI)取值的不同作为各地类的遥感分类依据之一,如表1。

表1 土地类型NDVI分布

Table1 The NDVI of diff erent land use types

土地类别NDVI

绿地0.10~0.80

沙化地0.00~0.10

沙地-0.10~0.00

研究中分类和信息提取是用监督分类方法结合人工解译来实现的。首先对研究区域进行实地考测,记录不同地物类型的空间分布,然后结合遥感影像数据和NDVI值进行监督分类。根据遥感影像解译的结果,判定1992年、1998年、2002年和2006年民勤绿洲湖区土地利用类型,计算土地面积和土地类型转移矩阵,如图1和表2、表3

所示。

图1 1992 2002土地类型转移图

F ig.1 M a p show ing tr ansfer of land

use types f rom1992to2002

表2 1992 2002土地类型转移矩阵

Table2 Transfer matrix of land use types from1992to2002

类型

绿地(2002年)沙化地(2002年)沙地(2002年)转移面积/hm2转换率/%转移面积/h m2转换率/%转移面积/hm2转换率/%绿地(1992年)15361.8462.695155.1121.043986.9216.27沙化地(1992年)3157.2619.5210082.2062.342934.5018.14沙地(1992年)1738.258.064192.1419.4415631.9672.50

表3 1992 2002土地面积变化

Table3 C hanges of land area of different land

use types from1992to2002

类型1992年1998年2002年绿地/hm224107.4927571.3624177.75沙化地/hm216928.1514219.1915799.89沙地/hm221613.6220451.5022747.20 1992 2002年,整个湖区绿地面积先增加后减少;沙化地和沙地面积先减少后增加,但沙地面积在后期快速增加。1992 2002年经历一次大规模的开荒过程,到1998年耕地面积达到最大值,之后由于缺乏管理以及水资源的减少,使耕地逐渐退化,这种现象在湖区的西渠,中渠以及收成乡较为明显。

803

第5期宋冬梅等:基于元胞自动机民勤绿洲湖区荒漠化演化预测