生态地理建模中的多尺度问题_岳天祥

合集下载
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

第23卷 第3期2003年5月
第 四 纪 研 究
QUATERNARY SCIENCES
Vol.23,No.3
May,2003
生态地理建模中的多尺度问题*
岳天祥 刘纪远
(中国科学院地理科学与资源研究所,100101 北京)
摘要 本文在分析生态地理建模内涵的基础上,讨论了生态地理建模中的尺度转换问题、跨尺度相互作用问题、空间尺度与时间尺度的关联问题和多尺度数据处理问题。

由于生态地
理问题的非线性、生态环境的异质性和随机事件,简单的线性尺度转换方法远不能满足生态
地理建模的要求。

为了从根本上解决生态地理建模中的时空尺度问题,除需要运用微分几何
学和等级理论等经典方法外,还需要引入格点生成法和网格计算等现代理论和技术手段。

主题词 生态地理建模 时空尺度 格点生成 网格计算 系统模拟
1 引言
生态地理建模(ecogeographical modeling)是对各种时空尺度的生态系统结构、空间格局和过程的定量描述。

它的直接目标是定量描述和分析分类单元(ta xa)的空间格局和过程,它的最终目标是模拟生态系统及其服务功能的空间演化规律[1]。

生态地理分类是生态地理建模的首要工作,已有很长的研究历史[2],例如,Humboldt (1907)关于气候与物种分布关系的描述;Schouw(1823)对植物区系和植被类型分布与各种气候因素相互关系的探讨;Griesbac h(1872)关于植物全球分布的描述;Merriam(1892)提出的基于温度的北美植物区系和动物区系的地带分布;Koeppen和Geiger(1930)提出的气候分类系统;在Koeppen方法的基础上,Thornthwaite(1931)提出的侧重植物、土壤和流域特征的气候分类;为了改进前人分类的缺陷,Holdridge(1947)提出了Holdridge生命地带模型,根据年平均生物温度、年平均总降雨量和潜在土壤水分蒸发蒸腾损失率,将全球划分为100多个生命地带。

我国在生态地理分类方面也做了许多工作。

例如,中国综合自然区划根据温度、干湿度和地形将中国划分为12个温度带和4类干湿区[3,4];中国综合农业区划将我国划分为10个一级农业区和38个二级农业区[5];为了使Holdridge生命地带模型更适合中国的实际情况,张新时[6]对其进行了修正,并将其应用于中国潜在净第一性生产力的估算; 中国生态区划方案 将我国划分为3个生态大区、13个生态地区和57个生态区[7];中国生态环境胁迫过程区划根据人类活动对自然资源和生态环境构成的压力,将中国划分为两个一级区、10个二级区和29个三级区[8];中国土地覆盖分类首先根据气候数据将中国划分
第一作者简介:岳天祥 男 39岁 研究员 资源环境模型与系统模拟专业 E mail:yue@
*国家 973 项目(批准号:2002CB412500)和国家自然科学重点基金项目(批准号:90202002)资助2003-01-26收稿,2003-02-12收修改稿
为9个生物气候区,并在集成AVHRR 卫星遥感影像与自然地理数据集的基础上对每个生物气候区进行了土地覆盖分类1)。

在完成生态地理分类之后,生态系统边界的变化探测(detection of boundary change)是生态地理建模的最重要步骤。

边界变化探测包括遥感、统计分析和数值模拟3种方法[9]。

遥感方法包括Sobel 算子、Laplace 算子、Hough 变换和边界跟踪法。

Sobel 算子通过计算相邻像元属性值的线性差或一阶导数找出变化率最高的地方[10];Laplace 算子通过计算临近像元属性值的二阶导数来分析边界的宽度[11];Hough 变换是为了减小搜寻计算的复杂性对几何形状的参数描述[12]。

边界跟踪法是指以某一初值为起点对满足边界阈值条件要素的连续探测[13]。

根据数据的具体情况,边界变化的统计分析包括4种途径:如果数据是沿横断面的采样值,则边界是相邻样点值方差最高的地方[14];如果样点规则地分布于研究区域,4个相邻样点可形成一个正方形,则用每个正方形的顶点样值可计算其一阶偏导数,一阶偏导数高的地方即为边界[15];如果样点随机地或不规则地分布于研究区域,则
相邻的3个样点可形成一个三角形,一阶偏导数高的三角形即为边界[16];如果有效数据
不仅包括定量变量,而且也包括半定量和定性变量,则可同时计算所有变量在相邻样点的异点(dismililarity),总体上异点高的位置即为边界[15]。

地理信息系统为集成遥感数据处理和采样数据统计分析提供了有效的工具,为数值模拟提供了便捷的平台。

然而,边界变化的数值模拟需要数据的足够积累和高性能的计算机硬件和软件。

近年来,在生态地理建模研究中,Holdridge 生命地带模型得到了较广泛的应用。

例如,20世纪50年代以来在气候变化和土地利用变化双重驱动力作用下,中国Holdridge 生命地带变化规律及其对生态环境稳定性影响的模拟分析[2];Saint Lucia 岛生态系统分布评估[17];气候变化对Costa Rica 热带森林地理分布和多样性的区域影响预测等[18]。

目前,生态地理建模主要使用的是统计归纳的方法,局限于特定的空间尺度和分类层次,妨碍了对其他空间尺度和分类层次问题的分析,很难用于分析跨尺度问题。

也就是说,生态地理建模目前需要研究的理论核心是多尺度问题。

2 多尺度问题
20世纪60年代,生态学家和地理学家就注意到了尺度问题的重要性[19,20]。

20世纪90年代以来,多尺度问题被称为生态学和地理学的新前缘,受到生态地理学界的高度重视[1]。

例如,为了认识生态格局、过程和尺度之间的关系及解决有关科学问题,美国环境
保护局建立了多尺度实验生态系统研究中心(MEERC)
[21];为了确定地质变化和植被动态之间的相互作用,J.D.Phillips [22]提出了4种尺度指标;为了解决全球变化影响的跨尺度
问题,G.D.Peternson [23]引入了等级理论; A.Stein 等[24]用地统计方法确定环境变量的最恰当空间和时间尺度;A.Veldkamp 等[25]提出了农业经济研究的多尺度系统方法[17];R.H.Gardner 等[26]提出了实验生态学的多尺度分析理论;R.Schulze [27]分析了气候变化农业水文响应的尺度问题、尺度类型和尺度转换的关键问题; ne 和W.B.Cohen [28]根
1)Liu J Y,Zhuang D F,Luo D et nd cover cl as sification of China:Integrated analysis of A VHRR i magery and geophysical data.International J ournal of Remote Sensing,2003(i n printi ng)
2573期 岳天祥等:生态地理建模中的多尺度问题
258第 四 纪 研 究 2003年
据分形的自相似性建立了针对MODIS数据的尺度转换方法;K.M.Konarska等[29]通过比较NOAA-AVHRR和Landsat TM数据集的分析结果,提出了空间尺度对生态系统服务功能评价影响的分析方法;Q.Gao等[30]通过对草地模拟模型中的非线性函数进行泰勒展开,减小了由于空间尺度变化产生的误差;吕一河和傅伯杰[31]对生态学中常见的尺度转换方法进行了分析和总结;岳天祥和刘纪远等[2]建立了多尺度生态多样性模型、提出了多尺度信息融合模型的框架[32,33]。

通过分析多尺度问题的研究进展可以发现,大多数研究侧重于某一个主题或某一主题的一个方面,在生态地理的大背景下对多尺度问题进行全面系统分析和研究的成果很少,远不能满足生态地理建模的需要。

生态地理是生态学和地理学内容的结合,生态地理建模必须涵盖生态学和地理学的所有尺度问题。

尺度(scale)是观测(observation)、过程(process)或过程模型(model of that process)在空间或时间方面的基准尺寸(a c haracteristic dimension or size),它囊括过程的离散状态或过程状态间的临界点[27]。

尺度转换、跨尺度(cross scale)相互作用、空间尺度与时间尺度相互关联和多空间尺度数据处理问题是生态地理建模需要研究的重要内容。

(1)尺度转换问题 尺度转换可区分为自下而上的尺度转换(upscaling)和自上而下的尺度转换(downscaling)。

自下而上的尺度转换是指高分辨率研究结果向低分辨率的转换;自上而下的尺度转换是指低分辨率研究结果向高分辨率的转换[20]。

生态地理建模需要研究的主要尺度转换问题包括:如何进行观测过程和数学模型的尺度转换,尺度转换如何影响变量灵敏性、空间异质性和系统可预测性,非线性响应被放大或减小的环境条件等。

(2)跨空间尺度相互作用问题 当某一空间尺度的事件或现象影响其他空间尺度的事件或现象时,就产生了跨尺度的联系和相互作用。

然而,以往的大多数生态地理研究是在特定的空间尺度进行的,对跨尺度相互作用的分析非常有限。

生态地理建模需要研究的主要问题包括:在不同空间尺度同时发生作用的驱动力分析,跨尺度相互作用的识别,生态地理事件在不同空间尺度的相互作用机制,跨尺度相互作用产生的非线性问题,跨尺度相互作用如何影响环境管理以及在政策制定中如何考虑跨尺度相互作用问题。

(3)空间尺度和时间尺度相互关联问题 生态地理过程的空间尺度和时间尺度往往是密切相关的。

例如,食物生产可以在一年的时间尺度和局部的空间尺度进行模拟分析;生态系统的水调节功能可以在多年的时间尺度和区域的空间尺度进行模拟分析;生态系统的气候调节功能必须在至少几十年的时间尺度和全球的空间尺度进行模拟分析。

也就是说,大空间尺度的变化对应大时间尺度的生态地理过程。

空间尺度和时间尺度的这种关联是否在生态地理建模中可以作为一种普遍规则运用是需要研究的重要内容之一。

(4)多空间尺度数据处理问题 生态地理建模必须处理各种不同空间尺度的数据。

目前,在生态地理问题的研究中,一般使用低(粗)分辨率数据,即使有高(细)分辨率数据,由于计算机计算容量的制约,大多数模拟分析也是将其进行平均处理后转换为低分辨率数据。

这种处理几乎损失了所有局部格局信息和非线性特征信息。

如何将高分辨率数据和低分辨率数据结合起来分析生态系统的结构、空间格局和过程是生态地理建模面临的首要挑战。

3 解决尺度问题的重要理论、方法和技术
有关研究表明[32,33],等级理论(hierarchy theory)、格点生成法(grid generation)和网格计
算(grid computing)是解决生态地理建模尺度问题的最重要理论基础和现代技术手段。

(1)等级理论 等级理论是生态系统尺度效应的传统分析方法,它是系统理论的一种延伸,用于分析尺度对复杂系统组织的影响。

等级理论把系统看作一个相互作用要素的集合,这个集合又形成一个更大系统的要素,它并不假定系统具有必要的等级结构,而是把系统分成不同的等级层次,以便于跨尺度相互作用的分析[23]。

(2)格点生成法 19世纪60年代出现了格点生成法[34]。

数值化格点生成方法的发展形成了一个独立的数学分支。

20世纪90年代,格点生成技术的发展达到了一个新的阶段。

格点生成可以被看作是寻找地图上有用的参数化信息的过程,格点生成技术应用的独特性在于它不局限于任何特定的公式和需要任何特定的理论基础。

格点生成技术最重要的步骤是找到合适的从自然现象到数值化信息的转换方法。

基于内插的格点生成方法有两大优势,即相对于局部微分方程的快速计算能力和对格点的直接控制。

(3)网格计算 网格是借鉴电力网的概念提出来的。

在使用电力时,不需要知道它是从哪个地方的发电站输送出来的,也不需要知道该电力是通过什么样的发电机产生的,不管是水力、火力还是核力。

网格也希望给用户提供的是与地理位置无关、与具体的计算机设备无关的通用的计算能力。

网格计算是利用互联网把分散在不同地理位置的计算机组织成一个 虚拟的超级计算机 ,其中每一台参与计算的计算机就是一个 节点 ,而整个计算是由成千上万个 节点 组成的 一张网格 ,所以这种计算方式叫网格计算[35]。

在网格计算中,首先要查清网格里所有可用资源,比如,哪些主机可供访问、还空置多少处理能力、数据库里可供使用的数据是什么、共享的应用程序是否已准备好以及共享主机采用何种文件系统等。

用户提交的任务要由系统来分配资源并控制其运行,包括要将其分配到哪些主机上运行、调用哪些数据、启动何种应用程序以及何时开始运行等。

4 讨论
在以往的大多数研究中,地理学者和生态学者往往在特定的地理和分类尺度开展研究。

这些地理和分类尺度包括生物个体(organism),种群(population),群落(community),生态系统(ecosystem),生物群落区(biome),生物圈(biosphere),岩石圈(lithosphere),大气圈(atmosphere)和水圈(hydrosphere)等。

生态地理建模的主要任务是定量化地描述生态地理类型在各种时空尺度的空间格局和过程。

不同尺度过程不是独立的,小尺度过程受大尺度过程的制约,大尺度过程是许多小尺度过程相互作用的结果。

在许多情况下,某一尺度过程和格局的相互作用会在较大尺度形成自组织系统。

在这些自组织系统中,起主导作用的变量只有不多几个,也就是说,对这些主导变量及其相互作用的识别是解决跨尺度问题的关键。

目前,大多数研究使用简单的线性尺度转换方法,例如,加权平均、内插和跨尺度数据转换等。

然而,由于生态地理问题的非线性、生态环境的异质性、跨尺度相互作用和随机事件,这些简单的线性尺度转换方法远不能满足生态地理建模的要求。

初步研究结果表明[34,33,36]
,运用等级理论,通过确定恰当的搜索半径,建立跨尺度分析模型,可以解决跨尺度问题中的多层次问题;在空间分类系统背景下的尺度转换方法可归纳为:运用格点生成法,通过建立规范化的格点生成规则、构建多尺度数据融合模型,实现生态地理建模中的多尺度转换;高分辨率和低分辨率遥感数据的尺度转换方法可归纳2593期 岳天祥等:生态地理建模中的多尺度问题
260第 四 纪 研 究 2003年
为:运用微分几何学曲面论基本定理,建立尺度转换模型的基本方程,运用遥感数据,反演尺度转换模型的一次近似表达形式,根据点上的采样数据或观测数据,修正一次近似表达形式,实现尺度转换模型的第二次近似表达,如果有更新的遥感数据和点上的采样数据或观测数据,重复以上过程,直至理论模型与实际需求完全相符;运用网格计算技术,解决由于跨尺度分析模型和尺度转换模型采用新理论方法带来的超计算机负荷问题。

参 考 文 献
1 Allen H D.M edi terranean Ecogeography.London:Prentice Hall,2001.2~10
2 Yue T X,Liu J Y,J rgensen S E e t al.Changes of Holdri dge life z one di versity in all of China over a half century.Ec olo gical
Modelling,2001,144:153~162
3 中国科学院自然区划工作委员会.中国综合自然区划(初稿).北京:科学出版社,1959.2~3
4 黄秉维.中国综合自然区划纲要(1986年).中国科学院地理研究所编.地理集刊(第21号).北京:科学出版社,
1990.10~20
5 全国农业区划委员会.中国农业自然资源和农业区划.北京:农业出版社,1991.71~99
6 张新时.植被的PE(可能蒸散)指标与植被-气候分类(三).植物生态学与地植物学学报,1993,17(2):97~109
7 傅伯杰,刘国华,陈利顶等.中国生态区划方案.生态学报,2001,21(1):1~6
8 苗 鸿,王效科,欧阳志云.中国生态环境胁迫过程区划研究.生态学报,2001,21(1):7~13
9 Fortin M J,Ols on R J,Ferson S et al.Is sues related to the de tection of ndsc ape Ecology,2000,15:453~466
10 Pitas I,Venetsanopolos A N.Edge de tectors based on order stati stics.IEEE Transactions on Pattern Anal ysis and M achine
Inte lligence,1986,8(4):538~550
11 Pi tas I.Di gi tal Image Proces sing Algori thms.Ne w York:Prentice Hall.1993.109~127
12 Chatzis V,Pitas I.Fuzzy cell Hough trans form for curve detecti on.Patte rn Rec ognition,1997,30:2031~2042
13 Pi tas I.Di gi tal Image Proces sing Algori thms and Applications.New York:John Wiley&Sons,Inc.2000.241~272
14 Ludwig J A,Corneli us J M.Locating discontinuiti es along ec ological gradients.Ecolo gy,1987,68:448~450
15 Fortin M J,Drapeau P.Deli neati on of ecological boundaries:Comparis on of approaches and significance tes ts.Oikos,1995,
72:323~332
16 Fortin M J.Edge de tection algorithms for two di mensi onal ecological data.Ec ology,1994,75:956~965
17 Isaac C,Bourque C.Ecological life zones of Sai nt Lucia.Global Ecology&Biogeography,2001,10:549~566
18 Enquis t C A F.Predic ted regional impacts of climate c hange on the geographical distri bution and diversity of tropical forests in
Costa Rica.Journal o f Bio geography,2002,29:519~534
19 Stommel H.Varieties of oceanographic experience.Sc ie nce,1963,139:572~576
20 Schumm S A,Lichty R W.Ti me,Space,and causality in Geomorphology.American Journal o f Scie nce,1965,263:110~119
21 Peternson D L,Parker V T.Ecological Scale.New York:Columbia University Press,1998.203~244
22 Phillips J D.Biogeo morphology and landscape evol ution:The problem of scale.Ge omorphology,1995,13:337~347
23 Peternson G D.Scaling ecological dynamics:Sel f organization,hierarchical structure,and ecological resilience.Climate
Change,2000,44:291~309
24 Stein A,Riley J,Halberg N.Iss ues of scale for environmental i ndicators.Agric ulture,Ecosyste ms and Environme nt,2001,82:
215~232
25 Veldkamp A,Kok K,de Koning G H J e t al.Multi scale s ys tem approaches in agronomic research at the landscape level.Soil
&Tillage Research,2001,58:129~140
26 Gardner R H,Kemp W M,Kennedy V S et al.Scaling Relations in Experi mental Ecology.Ne w York:Columbia University
Press,2001.9~19
27 Sc hulze R.Transcending scales of s pace and time in impact s tudies of cli mate and cli mate change on agrohydrological res ponses.
Agriculture,Ecosyste ms and Environme nt,2000,82:185~212
28 M i lne B T,Cohen W B.M ultiscale assess ment of binary and continuous landcover variables for MODIS vali dation,mapping,and
modeli ng applications.Re mote Se nsing o f Environme nt ,1999,70:82~98
29 Konarska K M ,Sutton P C,Castellon M.Evaluating scale dependence of ecosystem service valuation:A compari son of NOAA
AV HRR and Landsat TM datasets.Ecological Economics ,2002,41:491~507
30 Gao Q,Yu M,Yang X et al .Scaling si mulati on models for spati all y heterogeneous ecos ys tems with diffusi ve transportation.
L andsc ape Ecology ,2001,16:289~300
31 吕一河,傅伯杰.生态学中的尺度与尺度转换方法.生态学报,2001,21(12):2096~2105
32 岳天祥,刘纪远.多源信息融合数字模型.世界科技研究与发展,2001,23(5):19~25
33 岳天祥,刘纪远.第四代地理信息系统研究中的尺度转换数字模型.中国图像图形学报,2001,6(9):907~91134 Lis ei kin V D.Grid Generation Methods.Berlin:Springer Verlag,1999.3~8
35 Fos ter I,Kess el man C,Tuecke S.The anatomy of the grid:Enabling s calable virtual organiz ati ons.Inte rnational Journal
Supe rcompute r Applications ,2001,15(3):6~7
36 刘纪远,岳天祥,王英安等.中国人口密度数字模拟.地理学报,2003,58(1):17~24
ISSUES ON MULTI SCALES IN ECOGEOGRAPHIC AL MODELING
Yue Tianxiang Liu Jiyuan
(Institute o f Geographic al Sciences and Natural Resources Rese arch ,Chine se Academy o f Sc ienc es,Bei jing 100101)
Abstract
Ec ogeographical modeling is quantitatively describing ecosyste m structures,spatial patterns and processes on various temporal and spatial scales.Its direct aim is to describe and analyze spatial patterns and processes of ta xa.Its final aim is to simulate the spatial evolution la w of ecosyste ms and ecosystem services.Ecogeography is an integration of ec ology and geography so that ecogeographical modeling must deal with all scale issues involved in ecology and geography.Scale transformation,cross scale interaction,linkage between spatial scale and temporal scale,and multiscale data fusion are essential works of ecogeographical modeling.
The primary results show that hierarchy theory might be employed to solve problems of the cross scale interac tions by choosing appropriate search radius and constructing cross scale model.For spatial hierarchical systems,grid generation method is useful for solving problems of nonlinear scale transformation by establishing standardized grid generation re gulation and constructing model for multi scale information fusion.For scale transformation between different resolution pixels of remotely sensed data,the first step is to establish basic equation of digital model by means of surface theory;the second step is to retrieve first approximate formulation using available informa tion whic h would be repeated until required acc uracy is reached.Grid computing is an useful technique to solve problems of computing capacity overload.
Key words ec ogeographical modeling,temporal and spatial scale,grid generation,grid computing,syste m simulation 2613期 岳天祥等:生态地理建模中的多尺度问题。

相关文档
最新文档