全 文 ：　 　 倡 国家自然科学基金重大项目“荒漠化综合防治优化模式研究”（３９９９０４９０）和地方合作项目“大同市数字农业与数字生态建立与应用 ”资助
　 　 倡倡 通讯作者
收稿日期 ：２００６唱０６唱０９ 　 改回日期 ：２００７唱０１唱０９
大同市南郊区土地利用生态服务价值空间分布及变化 倡
郭青霞１ ，２ 　陈焕伟１ 倡倡 周 　欣３
（１ ．中国农业大学资源环境学院 　 北京 　 １０００９４ ； ２ ．山西农业大学资源环境学院 　 太谷 　 ０３０８０１ ；
３畅 大同市南郊区农业局 　 大同 　 ０３７００１）
摘 　 要 　 本研究以中国生态服务价值当量因子表为基础计算了大同市南郊区不同土地利用生态服务价值单价 ，结
合该区 １９８７ 年和 ２００２ 年 Landsat TM／ ETM ＋ 影像分类数据 ，借助 F ragstats 和 ENVI 软件支持进行了矩形单元土
地利用生态服务价值的空间模拟 ，然后采用 ENVI 的波段运算功能 、分类功能进行大同市南郊区土地利用生态服
务价值空间分布及变化研究 。结果显示 ：大同市南郊区 １９８７ 年和 ２００２ 年土地利用正生态服务区面积大于负生态
服务区面积 ，但生态服务质量等级都不高 。 １９８７ 年平原区生态服务质量好于山丘区 ，２００２ 年山丘区生态服务质量
好于平原区 。 大同市南郊区从 １９８７年至 ２００２ 年土地利用生态服务增加区面积大于减小区 ，且以轻度增加为主 。
平原区生态服务增加区面积小于山丘区 ，生态服务减小区面积大于山丘区 。 平原区生态服务增加区面积小于减小
区面积 ；山丘区生态服务增加区面积大于减小区面积 。
关键词 　 土地利用 　 生态服务价值 　 矩形单元 　 空间模拟 　 空间分布 　 空间变化
Spatial distribution and changes of ecological service value of land use in the southern suburb of Datong City ．GUO Qing唱
Xia１ ，２ ，CHEN Huan唱Wei１ ，ZHOU Xin３ （１畅College of Resources & Environmental Sciences ，China Agricultural Universit y ，
Beijing １０００９４ ，China ；２畅College of Resources & Environmental Sciences ，Shanxi Agricultural University ，Taigu ０３０８０１ ，
China ；３畅Agricultural Bureau of Southern Suburb of Datong City ，Datong ０３７００１ ，China） ，CJEA ，２００７ ，１５（５） ：１８８ ～ １９２
Abstract 　 Based on established equivalen t w eight factor of ecosystem service value of China ，ecological service value of land
use in the sou thern suburb of Datong City was estimated ．The spatial dist ribu tion and changes of ecological service value of
land use were analyzed using land use classification data f rom Landsat TM／ETM ＋ images fo r １９８７ and ２００２ and spatially
imitating land use ecological service value within rectangular units by Fragstat and ENVI soft ．Results show that positive
ecological service value land use area outweighs that of the negative for both １９８７ and ２００２ ．Ecological service value land
use quality of the plains is higher than that of the massifs in １９８７ ，while it is higher for the massifs in ２００２ ．Decreasing
ecological service value land area is larger than increasing ecological service value land area in the plains from １９８７ to
２００２ ，while the opposite t rend is the case in massifs and the entire southern suburb research area ．On the whole ，decreas唱
ing ecological service value land area in the plains is larger than that of t he massifs ，while increasing ecological service value
land area in the plains is smaller than that of the massifs ．
Key words 　 Land use ，Ecological service value ，Rectangular unit ，Spatial imitation ，Spatial dist ribu tion ，Spatial change
（Received June ９ ，２００６ ；revised Jan ．９ ，２００７）
大同市南郊区位于黄土高原典型生态脆弱区 ，水土流失严重 、水资源缺乏 、生物多样性小 。 因此 ，改善
大同市南郊区生态环境质量是保障当地社会经济可持续发展的重要前提 。 由于土地利用深刻影响着区域
水循环 、环境质量 、生物多样性及陆地生态系统生产力而成为全球环境变化研究的重要内容［１］ 。 因此合理
配置大同市南郊区土地利用方式对当地生态环境可持续发展具有重要意义 。 但是对于给定地块 ，选择合理
的土地利用方式 ，首先应了解土地利用的生态服务价值［１４］ 。
生态系统服务的概念是在 ２０世纪 ６０年代中期和 ７０ 年代早期提出的［１５ ，１６］ 。之后 ，国际上对生态系统服务研
究非常重视 ，开展了从生态系统过程 、生态服务功能及其生态经济价值等多方面的综合研究 ，不断充实与丰富了生
态系统服务的内涵［１７ ，１８］ ，完善了生态系统服务的分类［１８ ，１９］ ，探讨了生态服务评价技术及生态经济价值的评估方
法 。 ２０世纪 ９０ 年代末以来 ，国外有关生态系统服务的概念 、生态效益的价值理论及评估方法等开始引入国内［２］ ，
第 １５卷第 ５ 期 中 国 生 态 农 业 学 报 Vol ．１５ 　 No ．５
２ ０ ０ ７ 年 ９ 月 Chinese Journal of Eco唱Agriculture Sept ．，　 ２００７
国内也开始从区域生态系统服务价值 、单个生态系统服务价值和生态系统单项服务价值 ３ 个层面开展研
究［３］ 。 在区域生态系统服务研究方面 ，主要采用 ４ 种方法 ：一是应用国际上对生态系统服务评价技术和生
态经济价值评估方法进行区域生态系统服务的测算［４ ～ ８］ ；二是利用国际生态服务价值单价进行国内不同等
级 、不同类型区域生态系统服务的测算［２ ，９ ，１０］ ；第三种方法是以谢高地等［１１］建立的中国生态服务价值单价
和当量因子表为基础 ，进行区域生态系统价值测算［１１ ，１２］ 。 另外白晓飞等［１３］则是考虑区域生态系统服务项
目的实际 ，对具有生态服务价值单价的采用已有单价 ，对没有单价的采用生态系统服务评价技术和生态经
济价值评估方法进行生态服务价值单价的计算 ，从而快速地进行了区域内相对比较完整的生态服务项目的
价值估算 。
生态系统服务价值估算在国内外的进展为其应用于土地利用管理决策提供了基础 。 用于土地利用管
理决策的生态环境指标要具有可测性 、易用性 、时空可对比特征 ，基于价值单价基础上的生态服务价值测算
就具有这方面的特点 。 由于谢高地等［１１］建立的中国生态系统服务价值当量因子表更符合我国实际 ，且其采
用的全国标准农田的生态服务价值当量因子为 １ 的提法可以外推到更小的区域 ，通过小区域标准农田生态
服务价值单价的计算就可以得到小区域生态服务价值单价 。 同时宗跃光等［５］将 Costanza 等人单纯自然资
本的测算推广到自然资本与经济资本和社会资本的综合测算上 ，为此结合这两种方法进行大同市南郊区土
地利用生态服务价值单价计算具有可靠的理论基础 ，可为该区土地合理利用方式的选择提供依据 。
1 　 研究数据和方法
1畅1 　 土地利用现状数据的获取
土地利用研究数据来自 １９８７年 ９ 月 ８ 日的 Landsat５ TM 影像和 ２００２ 年 ９ 月 ２５ 日 Landsat７ ETM ＋ 影像 ，
为了保持分类条件的一致性 ，统一采用两期影像 ６ 个波段数据 ，以高分辨率卫星影像 SPOT 多光谱波段与全波
段的融合影像以及相应时期土地利用现状图为辅助数据 ，结合地面调查 ，确定影像训练区 。 为了提高分类精
度 ，采用光谱信息与纹理信息相结合的方法 ，在进行多种分类方法比较后 ，以分类精度较高的最大似然分类法
进行分类 ，两期影像分类结果的 Kappa 系数分别达到 ０畅９５６７ 和 ０畅９６２０ ，可以满足后续分析使用 。
1畅2 　 土地利用生态服务价值单价的确定
大同市南郊区土地利用生态服务价值计算共考虑 ９ 项生态服务功能 ：气体调节 、气候调节 、水源涵养 、土
壤形成与保护 、废物处理 、生物多样性维持 、食物生产 、原材料生产 、休闲娱乐 。 并假设各土地利用类型的各
项生态服务之间具有谢高地［１１］研究确定的中国生态系统服务价值当量因子之间的关系 ，取耕地食物生产生
态服务价值当量为 １ ，其他土地利用类型的各项生态服务价值当量为该土地利用类型的各项生态服务价值
与耕地食物生产生态服务价值的比值 。 耕地食物生产生态服务价值单价为 １hm２ 平均产量的农田每年自然
粮食产量的经济价值［１２］ ，计算公式为 ：
Cc ro p ＝ １７ ×
Tp v
Ta （１）
式中 ，Tp v为当年粮食作物总产值（元） ；Ta 为当年粮食作物播种面积（hm２ ） ；１／７ 是指没有人力投入的自然生
态系统的生态服务价值是现有单位面积农田提供的食物生产经济价值的 １／７ 。
依据大同市南郊区统计局统计资料 ，采用 １９９０ 年不变价分别计算 １９８７ 年和 ２００２ 年大同市南郊区耕地食
物生产生态服务价值单价后 ，采用公式（２）计算得到森林 、草地 、水域和未利用地的生态服务价值单价（表 １） 。
V i ＝ ∑
m
j ＝ １
Cc ro p × f i j 　 （ m ＝ １ ，２ ，… ，９） （２）
式中 ，V i 为第 i种土地利用生态服务的价值单价（元／hm２ ·a） ；f i j代表第 i种土地利用第 j 种生态服务的价值
当量因子 ；m 为土地利用生态服务项目数 。
　 　 由于大同市南郊区缺乏相
关资料 ，居民点及独立工矿用地
生态服务价值单价采用大同市
城区和矿用地有关资料计算 。
主要考虑大气调节 、水循环调节
和废弃物处理 ３ 项生态服务 ，并
假定其他方面的生态服务价值
表 1 　 大同市南郊区 1987和 2002 年土地利用生态系统服务价值单价
Tal畅１ 　 Ecological service value of land use in southern suburb of Datong City in １９８７ and ２００２
年份
Year
土地利用生态服务价值单价 ／元 ·hm － ２·a － １ Ecolo gical service value of land
农田
Cropland
森林
Forest
草地
Grassland
居民点及工矿用地
Build ing and f act ory
水域
Wat er
未利用地
Unused land
１９８７ １５９３畅８９ ５０４０畅０２ １６７０畅０１ － １１１９畅９４ １０６０３畅６４ ９６畅８８
２００２ １８３４畅３０ ５８００畅２２ １９２１畅９０ － １７２５畅６９ １２２０３畅０２ １１１畅４９
第 ５期 郭青霞等 ：大同市南郊区土地利用生态服务价值空间分布及变化 １８９　
单价为零 。 其中大气调节 、废弃物生态服务价值主要采用防治成本法 ，水循环价值采用替代成本法估计 。
城市废弃物处理生态服务价值单价等于城市工业废弃物和生活废弃物年排放量乘以各自处理价格加和后
除以居民点及独立工矿用地面积 ；水循环生态服务价值单价等于生活和工业用水总量乘以各自的价格加和
后除以居民点及独立工矿用地面积 ；气体调节价值等于城市废气排放总量乘以城市单位废气控制成本除以
居民点及独立工矿用地面积 。
1畅3 　 矩形单元土地利用生态服务价值的模拟方法
生态服务价值单价是以土地利用斑块作为单元的 ，但实践中研究各土地利用类型组合价值更有意义 。 因
此对于任一由诸多土地利用类型组成的实体 ，采用矩形单元分割并表达其价值 ，有利于空间对比分析 。 景观类
型水平指数缀块类型总面积 CA或缀块类型面积百分比 PLAND 以及景观的移动窗口分析方法为进行矩形单
元生态服务价值模拟提供了便利 。 这里主要采用移动窗口计算景观缀块类型面积百分比 PLAND 的方法 ，首
先借助 Fragstats软件矩形移动窗口分析法计算景观类型水平 PL AND 指标值 ，然后将其 AscⅡ码格式的结果文
件导入 ENVI软件中 ，采用公式（３） ，则可以模拟得到代表矩形单元的中心像元生态服务价值 ：
Ser v icece ll ＝ Scell × ∑
n
i ＝ １
PLANDi × Vi （３）
式中 ，PLANDi 为窗口内第 i种土地利用所占的面积比重 ，Scell代表每个矩形单元中心像元的面积大小（hm２） 。
进一步借助 ENVI 的波段运算功能得到这 １５ 年期间大同市南郊区土地利用生态服务价值变化空间分
布图 。 根据各土地利用生态服务价值及变化分布图的统计直方图信息 ，确定分类阈值 ，然后采用监督分类
法得到研究区土地利用生态服务价值及变化分区 。
2 　 结果与分析
经测算 ，大同市南郊区 １９８７ 年和 ２００２ 年的土地利用生态服务总价值分别为 １１１２９畅５１ 万元和
１７６９３畅２９ 万元 ，２００２ 年较 １９８７ 年净增加 ６５６３畅７８ 万元 。
2畅1 　 1978 年土地利用生态服务价值空间分布特征
表 2 　大同市南郊区 1987年和 2002年土地利用生态服务价值分区特征
Tab畅２ 　 Characteristics of ecological service value regions of land use
in southern suburb of Datong Cit y in １９８７ and ２００２
生态服务价值分区
Region of ecological
se rvice value
１９８７ ２００２
面积／h m２
Area
百分比 ／ ％
Percentage in
to tal a rea
面积 ／hm ２
Area
百分比 ／ ％
Percentag e in
t otal area
负生态服务区 ８５５６畅１２ ７畅８８ １３５８２畅７１ １２畅５２
零生态服务区 ３６５７０畅２４ ３３畅７０ ２１８８８畅２７ ２０畅１７
低生态服务区 ５６２４３畅６１ ５１畅８２ ５２５０３畅０３ ４８畅３８
中生态服务区 ６３４２畅３０ ５畅８４ １９６７１畅１２ １８畅１３
高生态服务区 ８１６畅８４ ０畅７５ ８８３畅９８ ０畅８１
　 　 根据 １９８７ 年和 ２００２ 年大同市南郊区土
地生态服务价值分布直方图 ，影像像元主要
集聚分布于 － １００ 元 、０ 元 、１５０ 元 、４５０ 元 、
１０００ 元附近 。 采用阈值分割法对影像进行分
类 ，阈值分割标准为 ：－ ２００ ～ － １０ 元为负生
态服务区 ，－ １０ ～ １０ 元为零生态服务区 ，１０ ～
３００ 元为低生态服务区 ，３００ ～ ６００ 元为中生
态服务区 ，６００ ～ １２００ 元为高生态服务区 。 分
类结果见表 ２ 、表 ３ 和表 ４ 。 从表 ２ 可以看出 ，
１９８７ 年大同市南郊区土地利用正生态服务区
所占面积比重大于负生态服务区 ，但生态服
　 　务质量等级不
高 ，以低生态服
务区占比重最
大 ，其次是零生
态服务区和负
生 态 服 务 区 。
从表 ３ 可以看
出 ，在不同地貌
单元下 ，平原区
总体生态服务
优 于 山 丘 区 。
　 　 　 　
表 3 　 1987 年不同地貌单元土地利用生态服务价值分区特征
Tab畅３ 　 Characteristics of ecological service value regions of land use in different relief in １９８７
生态服务
价值分区
Region of
ecological
se rvice value
平原区
面积 ／hm ２
Area of
plain
占平原面积
百分比／ ％
Percent age
in plain
占各生态服务区
面积百分比／ ％
Percent age in the
corresponding
se rvice region
山丘区
面积 ／hm ２
Area of
massif
占山丘
区面积
百分比／ ％
Percent age
in massif
占各生态服务区
面积百分比 ／ ％
Percent age in the
corresponding
se rvice region
负生态服务区 ５６６７畅６６ １０畅８０ ６６畅２４ ２８８８畅４６ ５畅１８ ３３畅７６
零生态服务区 ８８５８畅６７ １６畅８７ ２４畅４０ ２７４４８畅１８ ４９畅２２ ７５畅６０
低生态服务区 ３３９０７畅８２ ６４畅５９ ６０畅２９ ２２３３５畅８０ ４０畅０５ ３９畅７１
中生态服务区 ３３５２畅２３ ６畅３９ ５２畅８６ ２９９０畅０７ ５畅３６ ４７畅１４
高生态服务区 ７１３畅２５ １畅３６ ８７畅３２ １０３畅５９ ０畅１９ １２畅６８
虽然平原区负生态服务价值比重大于山丘区 ，但零生态服务区面积小于山丘区 ，各正生态服务区面积都大
１９０　 中 国 生 态 农 业 学 报 第 １５ 卷
于山丘区 。 从平原区内部看 ，以低生态服务价值区占面积比重最大 ，其次是零生态服务区和负生态服务区 。
而在山丘区内部 ，以零生态服务区所占面积比重最大 ，然后是低生态服务区 。 因此 ，各地貌类型内部土地利
用生态服务等级都不高 。
2畅2 　 2002 年土地利用生态服务价值空间分布特征
　 　 从表
２ 可 以 看
出 ，大同市
南 郊 区
２００２ 年正
生态服务
区面积大
于负生态
服务区 ，以
低生态服
　 　 　 　
表 4 　 2002 年不同地貌单元土地利用生态服务价值分区特征
Tab畅４ 　 Characteristics of ecological service value regions of land use in different relief in ２００２
生态服务
价值分区
Region o f
ecological
se rvice value
平原区
面积 ／hm ２
Area of
p lain
占平原面积
百分比／ ％
Percentag e
in plain
占各生态服务区
面积百分比 ／ ％
Percentage in the
corresponding
service region
山丘区
面积／h m２
Area o f
massif
占山丘
区面积
百分比 ／ ％
Percent age
in massif
占各生态服务区
面积百分比／ ％
Percent age in the
corresponding
se rvice reg ion
负生态服务区 １０６７５畅０８ ２０畅３３ ７８畅５９ ２９０７畅６３ ５畅２１ ２１畅４１
零生态服务区 ９４０７畅１３ １７畅９２ ４３畅５０ １２２１７畅７５ ２１畅９１ ５６畅５０
低生态服务区 ２７０７６畅０５ ５１畅５７ ５１畅５７ ２５４２６畅９８ ４５畅６０ ４８畅４３
中生态服务区 ４６３６畅３１ ８畅８３ ２３畅５７ １５０３４畅８２ ２６畅９６ ７６畅４３
高生态服务区 ７０５畅０６ １畅３４ ７９畅７６ １７８畅９２ ０畅３２ ２０畅２４
务区占面积比重最大 ，其次是零生态服务区和中生态服务区 ，总体生态服务质量不高 。 但与 １９８７ 年相比 ，虽
然负生态服务区面积增加 ，但是零生态服务区和低生态服务区面积都有所下降 ，而中生态服务区面积和高
生态服务区面积增加 。 这与大同市南郊区城镇用地扩展和生态环境建设是密不可分的 。
从表 ４ 可以看出 ，２００２ 年平原区的土地利用生态服务状况较差 ，负生态服务区和零生态服务区占比重
大于山丘区 ，中生态服务区占的比重小于山丘区 。 而山丘区生态服务状况较平原区好 。 在平原区内部 ，
２００２ 年以低生态服务区占面积比重最大 ，其次是负生态服务区和零生态服务区 。 在山丘区 ，以低生态服务
区占比重最大 ，其次是中生态服务区和零生态服务区 。 各地貌单元内部生态服务都处于较低水平 。
2畅3 　 15 年间土地利用生态服务价值变化的空间分布特征
表 5 　 大同市南郊区 1987 ～ 2002 年土地利用生态服务价值变化分区特征
Tab畅５ 　 Characteristic of ecological service value change regions of land
use in southern suburb of Datong Cit y betw een １９８７ and ２００２
生态服务价值变化分区
Region of eco lo gical service value change
面积 ／hm ２
Area
百分比／ ％
Percentage
高度增加区 １８０畅９９ ０畅１７
中度增加区 ９７８４畅９８ ９畅０１
轻度增加区 ３２６５５畅６９ ３０畅０９
不变化区 ４０６４４畅８１ ３７畅４５
轻度减小区 ２４１８４畅５３ ２２畅２８
中度减小区 １０９１畅４３ １畅０１
　 　 将大同市南郊区土地利用生态服务价值变化
进行等间距分割 ，进行生态服务价值变化分区 ：
＜ － １０５０ 元为极高减小区 ，－ １０５０ ～ － ７５０元为高
度减小区 ，－ ７５０ ～ － ４５０ 元为中度减小区 ，－ ４５０ ～
－ １５０ 元为轻度减小区 ，－ １５０ ～ １５０ 元为不变化
区 ，１５０ ～ ４５０ 元为轻度增加区 ，４５０ ～ ７５０ 元为中
度增加区 ，７５０ ～ １０５０ 元为高度增加区 ，＞ １０５０ 元
为极高增加区 ，分区结果见表 ５ 。 从表 ５ 可以看
出 ，大同市南郊区在 １９８７ 年至 ２００２ 年期间 ，土地
利用生态服务价值增加区的面积为 ４２６２１畅６６hm２ ，
　 　 　 　大于生态服务价值减小区的面积 ２５２７５畅９６hm２ ，且在土地生态服务增加区中 ，主要以轻度生态服务增加为
主 。 不变化区在各变化等级分区中面积最大 ，为 ４０６４４畅８１hm２ 。
　 　 在 不
同 地 貌 单
元间 ，平原
区 的 生 态
服 务 价 值
不 变 化 区
面 积 大 于
山丘区 ，生
态 服 务 价
值 增 加 区
面 积 小 于
山丘区 ，生
表 6 　 1987 年 ～ 2002 年不同地貌单元土地利用生态服务价值变化分区特征
Tab畅６ 　 Characteristic of ecological service value change regions of land use in different relief between １９８７ and ２００２
生态服务价
值变化分区
Region o f
ecological se rvice
value change
平原区
面积 ／hm ２
Area of
p lain
占平原面积
百分比／ ％
Percentag e
in plain
占各生态服务变化
区面积百分比／ ％
Percentage in the
corresponding
se rvice chang e region
山丘区
面积／h m２
Area o f
massif
占山丘
区面积
百分比 ／ ％
Percent age
in massif
占各生态服务变化
区面积百分比 ／ ％
Percent age in the
corresponding
service change region
高度增加区 　 １３６畅０８ ０畅２６ ７５畅１９ 　 ４４畅９１ ０畅０８ ２４畅８１
中度增加区 　 １０４０畅８５ １畅９８ １０畅６４ 　 ８７４４畅１３ １５畅６８ ８９畅３６
轻度增加区 　 １００９８畅９９ １９畅２３ ３０畅９３ 　 ２２５５６畅７０ ４０畅４５ ６９畅０７
不变化区 　 ２１２６４畅３２ ４０畅４９ ５２畅６６ 　 １９１１７畅１０ ３４畅２８ ４７畅３４
轻度减小区 　 １９１５５畅８７ ３６畅４８ ７９畅２１ 　 ５０２８畅６６ ９畅０２ ２０畅７９
中度减小区 　 ８１６畅６６ １畅５６ ７４畅８２ 　 ２７４畅７７ ０畅４９ ２５畅１８
第 ５期 郭青霞等 ：大同市南郊区土地利用生态服务价值空间分布及变化 １９１　
　 　 　 　态服务价值减小区的面积大于山丘区 。 在平原区内部 ，生态服务价值不变化区占比重最大 ，轻度减小区面
积为第二 ，生态服务价值增加区面积小于生态服务价值减小区面积 ；在山丘区 ，生态服务价值轻度增加区面
积最大 ，不变化区面积居第二 ，生态服务价值增加区面积大于生态服务价值减小区面积（表 ６） 。
3 　 结论与讨论
矩形单元的土地利用生态服务价值模拟方法 ，可以突破土地利用生态服务价值计量只能局限于土地利
用单元的限制 ，进行组合土地利用单元生态服务价值的准确计算 。 有利于实现组合土地利用单元生态服务
价值的空间对比及生态服务价值的空间分区研究 ，为土地利用生态服务价值分区管理提供依据 。
大同市南郊区土地利用生态服务价值研究显示 ：大同市南郊区 １９８７ 年和 ２００２ 年的土地利用生态服务
总价值分别为 １１１２９畅５１ 万元和 １７６９３畅２９ 万元 ，２００２ 年较 １９８７ 年净增加了 ６５６３畅７８ 万元 。 生态服务价值
增加区面积大于生态服务价值减小区面积 。 分区研究显示 ：１９８７ 年和 ２００２ 年正生态服务区面积都大于负
生态服务区 ，但服务等级不高 ，以低生态服务区占面积比重最大 ；但 ２００２ 年较 １９８７ 年负生态服务区面积增
加 ，零生态服务区和低生态服务区面积下降 ，中生态服务区面积和高生态服务区面积增加 。 这与大同市南
郊区城镇用地扩展和生态环境建设有关 。 不同地貌间土地利用生态服务价值研究显示 ：１９８７ 年平原区总体
生态服务优于山丘区 ，２００２ 年山丘区生态服务状况较平原区好 。 大同市南郊区 １９８７ 年至 ２００２ 年期间的土
地利用变化分区研究显示 ，土地利用生态服务价值增加区面积为 ４２６２１畅６６hm２ ，大于生态服务价值减小区面
积 ２５２７５畅９６hm２ 。 平原区生态服务价值增加区面积小于山丘区 ，生态服务价值减小区面积大于山丘区 。 平
原区内部生态服务价值增加区面积小于生态服务价值减小区面积 ；而山丘区则是生态服务价值增加区面积
大于生态服务价值减小区面积 。 平原区面临生态退化风险 ，而山丘区总体生态状况有所改善 ，但平原区和
山丘区在 ２００２ 年生态服务等级都相对较低 。 因此今后在土地利用决策中 ，不仅要注意山丘区生态服务的提
高 ，也应该改善平原区土地利用 。 尤其是城市用地的增加 ，必然加大生态服务负价值 ，应通过林草地等高生
态服务用地的补充来减缓平原区生态服务价值减小的趋势 。 在平原区耕地压力比较大的情况下 ，应注重农
林牧综合发展 ，加大城区绿化用地的数量 。
参 　 考 　 文 　 献
１ 　 李秀彬 ．全球环境变化研究的核心领域 ——— 土地利用 ／土地覆被变化的国际研究动向 ．地理学报 ，１９９６ ，５１（６） ：５５３ ～ ５５８
２ 　 张志强 ，徐中民 ，王 　 建 ，等 ．黑河流域生态系统服务的价值 ．冰川冻土 ，２００１ ，２３（４） ：３６０ ～ ３６７
３ 　 徐中民 ，张志强 ，程国栋 ．生态经济学理论方法与应用 ．郑州 ：黄河水利出版社 ，２００３ ．１２６ ～ １３０
４ 　 欧阳志云 ，王效科 ，苗 　 鸿 ．中国陆地生态系统服务功能及其生态经济价值的初步研究 ．生态学报 ，１９９９ ，１９（５） ：６０７ ～ ６１３
５ 　 宗跃光 ，陈红春 ，郭瑞华 ，等 ．地域生态系统服务功能的价值结构分析 ．地理研究 ，２００２ ，１９（２） ：１４８ ～ １５５
６ 　 高旺盛 ，董孝斌 ．黄土高原丘陵沟壑区脆弱农业生态系统服务评价 ．自然资源学报 ，２００３ ，１８（２） ：１８２ ～ １８８
７ 　 何 　 浩 ，潘耀忠 ，朱文泉 ，等 ．中国陆地生态系统服务价值测量 ．应用生态学报 ，２００５ ，１６（６） ：１１２２ ～ １１２７
８ 　 谢高地 ，肖 　 玉 ，甄 　 霖 ，等 ．我国粮食生产的生态服务价值研究 ．中国生态农业学报 ，２００５ ，１３（３） ：１０ ～ １３
９ 　 陈仲新 ，张新时 ．中国生态系统效益的价值 ．科学通报 ，２０００ ，４５（１） ：１７ ～ ２３
１０ 高清竹 ，何立环 ，黄晓霞 ，等 ．黄河上游农牧交错地区生态系统服务价值的变化 ．自然资源学报 ，２００２ ，１７（６） ：７０６ ～ ７１２
１１ 谢高地 ，鲁春霞 ，冷允法 ，等 ．青藏高原生态资产的价值评估 ．自然资源学报 ，２００３ ，１８（２） ：１８９ ～ １９６
１２ 肖 　 玉 ，谢高地 ，安 　 凯 ．莽措湖流域生态系统服务功能经济价值变化研究 ．应用生态学报 ，２００３ ，１４（５） ：６７６ ～ ６８０
１３ 白晓飞 ，陈焕伟 ．不同土地利用结构生态系统服务功能价值的变化研究 ．中国生态农业学报 ，２００４ ，１２（１） ：１８０ ～ １８２
１４ Goulder L ．H ．，Kennedy D ．Valuing ecosyst em services ：Philosophical bases and empirical methods ．Na ture’s Services ：Societal Dependence on
Natural Ecosyst em s ．Washingt on ，D ．C ．：Island Press ，１９９７
１５ King R ．T ．Wildli fe and man ．New York Conservat ionist ，１９６６ ，２０（６） ：８ ～ １１
１６ Helli well D ．R ．Valuat ion o f w ildlife resources ．Regional S tudies ，１９６９ ，３ ：４１ ～ ４９
１７ De Groot R ．S ．F unct ions o f nature ：evaluation of nature in environment al planning ，management and decision making ．Nature’s Services ：Soci唱
et al D ependence on Natural Ecosystems ．Washing ton ，D ．C ．：I sland Press ，１９９７
１８ C ostnaza R ．，D’Arge R ．，De Groo t R ．S ．，et a l ．T he value o f the w orld’s ecosystem services and natural capital ．Nature ，１９９７ ，３８７（６６３０） ：
２５３ ～ ２６０
１９ De Groot R ．S ．，Wilson M ．A ．，Boumans R ．M ．J ．A typolog y fo r the classificat ion ，descrip tion and valua tion o f ecosystem funct ions ，g oods and
services ．Ecological Economics ，２００２ ，４１（３） ：３９３ ～ ４０８
１９２　 中 国 生 态 农 业 学 报 第 １５ 卷