全 文 ：
摇 摇 摇 摇 摇 生 态 学 报
摇 摇 摇 摇 摇 摇 摇 渊杂匀耘晕郧栽粤陨 载哉耘月粤韵冤
摇 摇 第 猿猿卷 第 圆员期摇 摇 圆园员猿年 员员月摇 渊半月刊冤
目摇 摇 次
前沿理论与学科综述
生态系统生产总值核算院概念尧核算方法与案例研究 欧阳志云袁朱春全袁杨广斌袁等 渊远苑源苑冤噎噎噎噎噎噎噎噎
气候变化对传染病爆发流行的影响研究进展 李国栋袁张俊华袁焦耿军袁等 渊远苑远圆冤噎噎噎噎噎噎噎噎噎噎噎噎
好氧甲烷氧化菌生态学研究进展 贠娟莉袁王艳芬袁张洪勋 渊远苑苑源冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
氮沉降强度和频率对羊草叶绿素含量的影响 张云海袁何念鹏袁张光明袁等 渊远苑愿远冤噎噎噎噎噎噎噎噎噎噎噎噎
世界蜘蛛的分布格局及其多元相似性聚类分析 申效诚袁张保石袁张摇 锋袁等 渊远苑怨缘冤噎噎噎噎噎噎噎噎噎噎噎
风向因素对转基因抗虫棉花基因漂移效率的影响 朱家林袁贺摇 娟袁牛建群袁等 渊远愿园猿冤噎噎噎噎噎噎噎噎噎噎
个体与基础生态
长江口及东海春季底栖硅藻尧原生动物和小型底栖生物的生态特点 孟昭翠袁徐奎栋 渊远愿员猿冤噎噎噎噎噎噎噎
长江口横沙东滩围垦潮滩内外大型底栖动物功能群研究 吕巍巍袁马长安袁余摇 骥袁等 渊远愿圆缘冤噎噎噎噎噎噎噎
沣河沿岸土壤和优势植物重金属富集特征和潜在生态风险 杨摇 阳袁周正朝袁王欢欢袁等 渊远愿猿源冤噎噎噎噎噎噎
盐分和底物对黄河三角洲区土壤有机碳分解与转化的影响 李摇 玲袁仇少君袁檀菲菲袁等 渊远愿源源冤噎噎噎噎噎噎
短期夜间低温胁迫对秋茄幼苗碳氮代谢及其相关酶活性的影响 郑春芳袁刘伟成袁陈少波袁等 渊远愿缘猿冤噎噎噎噎
猿圆个切花菊品种的耐低磷特性 刘摇 鹏袁陈素梅袁房伟民袁等 渊远愿远猿冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
年龄和环境条件对泥蚶富集重金属镉和铜的影响 王召根袁吴洪喜袁陈肖肖袁等 渊远愿远怨冤噎噎噎噎噎噎噎噎噎噎
角倍蚜虫瘿对盐肤木光合特性和总氮含量的影响 李摇 杨袁杨子祥袁陈晓鸣袁等 渊远愿苑远冤噎噎噎噎噎噎噎噎噎噎
多噬伯克霍尔德氏菌 宰杂鄄云允怨对草甘膦的降解特性 李冠喜袁吴小芹袁叶建仁 渊远愿愿缘冤噎噎噎噎噎噎噎噎噎噎
金龟甲对蓖麻叶挥发物的触角电位和行为反应 李为争袁杨摇 雷袁申小卫袁等 渊远愿怨缘冤噎噎噎噎噎噎噎噎噎噎噎
种群尧群落和生态系统
白洋淀生态系统健康评价 徐摇 菲袁赵彦伟袁杨志峰袁等 渊远怨园源冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
珠海鹤洲水道沿岸红树林湿地大型底栖动物群落特征 王摇 卉袁钟摇 山袁方展强 渊远怨员猿冤噎噎噎噎噎噎噎噎噎
典型森林和草地生态系统呼吸各组分间的相互关系 朱先进袁于贵瑞袁王秋凤袁等 渊远怨圆缘冤噎噎噎噎噎噎噎噎噎
抚育间伐对油松人工林下大型真菌的影响 陈摇 晓袁白淑兰袁刘摇 勇袁等 渊远怨猿缘冤噎噎噎噎噎噎噎噎噎噎噎噎噎
百山祖自然保护区植物群落 遭藻贼葬多样性 谭珊珊袁叶珍林袁袁留斌袁等 渊远怨源源冤噎噎噎噎噎噎噎噎噎噎噎噎噎
土霉素对堆肥过程中酶活性和微生物群落代谢的影响 陈智学袁谷摇 洁袁高摇 华袁等 渊远怨缘苑冤噎噎噎噎噎噎噎噎
景观尧区域和全球生态
兴安落叶松针叶解剖结构变化及其光合能力对气候变化的适应性 季子敬袁全先奎袁王传宽 渊远怨远苑冤噎噎噎噎
盐城海滨湿地景观演变关键土壤生态因子与阈值研究 张华兵袁刘红玉袁李玉凤袁等 渊远怨苑缘冤噎噎噎噎噎噎噎噎
半干旱区沙地芦苇对浅水位变化的生理生态响应 马赟花袁张铜会袁刘新平 渊远怨愿源冤噎噎噎噎噎噎噎噎噎噎噎
杂宰粤栽模型融雪模块的改进 余文君袁南卓铜袁赵彦博袁等 渊远怨怨圆冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
科尔沁沙地湖泊消涨对气候变化的响应 常学礼袁赵学勇袁王摇 玮袁等 渊苑园园圆冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎
贝壳堤岛 猿种植被类型的土壤颗粒分形及水分生态特征 夏江宝袁张淑勇袁王荣荣袁等 渊苑园员猿冤噎噎噎噎噎噎噎
三峡库区古夫河着生藻类叶绿素 葬的时空分布特征及其影响因素 吴述园袁葛继稳袁苗文杰袁等 渊苑园圆猿冤噎噎噎
资源与产业生态
煤炭开发对矿区植被扰动时空效应的图谱分析要要要以大同矿区为例 黄摇 翌袁汪云甲袁李效顺袁等 渊苑园猿缘冤噎噎
学术信息与动态
叶中国当代生态学研究曳新书推介 刘某承 渊苑园源源冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
期刊基本参数院悦晕 员员鄄圆园猿员 辕 匝鄢员怨愿员鄢皂鄢员远鄢猿园园鄢扎澡鄢孕鄢 预 怨园郾 园园鄢员缘员园鄢猿园鄢圆园员猿鄄员员
室室室室室室室室室室室室室室
封面图说院 百山祖保护区森林植物群落要要要百山祖国家级自然保护区位于浙西南闽浙交界处袁由福建武夷山向东北伸展而成袁
主峰海拔 员愿缘远援苑皂袁为浙江省第二高峰遥 其独特的地形和水文地理环境形成了中亚热带气候区中一个特殊的区域袁
保存着十分丰富的植物种质资源以及国家重点保护野生动植物种袁尤其是 员怨愿苑年由国际物种保护委员会列为世界
最濒危的 员圆种植物之一的百山祖冷杉袁是第四纪冰川的孑遗植物袁素有野活化石冶之称遥 随着海拔的升高袁其植被为
常绿阔叶林尧常绿鄄落叶阔叶混交林尧针阔混交林尧针叶林尧山地矮林和山地灌草丛遥
彩图及图说提供院 陈建伟教授摇 北京林业大学摇 耘鄄皂葬蚤造院 糟蚤贼藻泽援糟澡藻灶躁憎岳 员远猿援糟燥皂
第 33 卷第 21 期
2013年 11月
生 态 学 报
ACTA ECOLOGICA SINICA
Vol.33,No.21
Nov.,2013
http: / / www.ecologica.cn
基金项目:国家自然科学基金资助项目(51174287);国家自然科学基金青年资助项目(41101428,71103182);江苏省普通高校研究生科研创新计
划资助项目(CXZZ12_0942);江苏高校优势学科建设工程资助项目(SZBF2011鄄6鄄B35);山西省科技重大专项资助项目(20121101008)
收稿日期:2012鄄07鄄08; 摇 摇 修订日期:2012鄄10鄄26
*通讯作者 Corresponding author.E鄄mail: wyj4139@ cumt.edu.cn
DOI: 10.5846 / stxb201207080952
黄翌,汪云甲,李效顺.胡召玲,刘国平.煤炭开发对矿区植被扰动时空效应的图谱分析 ———以大同矿区为例.生态学报,2013,33(21):7035鄄7043.
Huang Y, Wang Y J, Li X S, Hu Z L, Liu G P.Graphic analysis of spatio鄄temporal effect for vegetation disturbance caused by coal mining: a case of
Datong Coal Mine Area.Acta Ecologica Sinica,2013,33(21):7035鄄7043.
煤炭开发对矿区植被扰动时空效应的图谱分析
———以大同矿区为例
黄摇 翌1,2,汪云甲1,2,*,李效顺2,胡召玲3,刘国平3
(1. 矿山生态修复教育部工程研究中心 中国矿业大学,徐州摇 221116;
2. 江苏省资源环境信息工程重点实验室 中国矿业大学,徐州摇 221116; 3. 江苏师范大学城市与环境学院,徐州摇 221116)
摘要:植被是综合体现矿区生态环境状况的关键因子。 以 NDVI值为参数,基于像元二分的遥感估算方法测度了 1999—2010年
大同矿区的植被覆盖情况。 利用 Global Moran忆s I、Getis鄄Ord General G、Getis鄄Ord Gi*、Anselin Local Moran I等空间关联指数分别
从全局演变和局部效应的视角揭示了矿区植被受煤炭开发等活动扰动的时空演化和内在作用机制。 研究结果表明:(1)大同
矿区植被盖度变化在空间上呈集聚分布状态,表现出很强的空间自相关性,这一现象在植被盖度最低时期最明显;(2)在植被
盖度上升期,其变化的热点区不断减少或向次热点地区分散,矿区植被变化集聚效应越发不显著,总体上呈离散的态势。 (3)
煤炭开发对矿区植被的扰动呈现整体性特征,与煤炭产量成正比。 (4)不同植被类型受煤炭开发扰动呈现不同的变化特征。
(5)采用循环工艺的高效高产矿井大规模煤炭开发对植被的扰动作用尚不明显,但长期效应有待进一步研究。
关键词: 煤矿区; 植被盖度; 时空效应; 图谱分析
Graphic analysis of spatio鄄temporal effect for vegetation disturbance caused by
coal mining: a case of Datong Coal Mine Area
HUANG Yi1,2, WANG Yunjia1,2,*, LI Xiaoshun2, HU Zhaoling3, LIU Guoping3
1 Engineering Research Center of Mine Ecological Construction, Ministry of Education, China University of Mining & Technology, Xuzhou 221116, China
2 Jiangsu Key Laboratory of Resources and Environmental Information Engineering, China University of Mining & Technology, Xuzhou 221116, China
3 School of Urban and Environmental Science, Jiangsu Normal University, Xuzhou 221116, China
Abstract: Vegetation is a key factor which comprehensively reflects condition of ecology and environment in mining area.
With NDVI index as the parameter, vegetation coverage in Datong coal mining area from 1999 to 2010 were measured based
on dimidiate pixel model by remote sensing estimation method. Spatial correlation indexes such as Global Moran忆s I, Getis鄄
Ord General G, Getis鄄Ord Gi* and Anselin Local Moran I were respectively used to reveal the temporal evolution and
intrinsic mechanism of vegetation disturbed by coal mining from the global evolution and local effect perspective. Research
results showed that: (1) Changing of vegetation coverage in Datong coal mining area presented cluster distribution and
strong spatial autocorrelation in space, meanwhile, the phenomenon was absolutely obvious when vegetation coverage was in
the lowest period. ( 2) In the rising period of vegetation coverage, hotspots of vegetation coverage gradually reduced or
separated to the secondary hotspots area, clustering effect of vegetation changing became less significantly, which presented
discrete state in general. (3) The disturbance of coal mining brought to vegetation showed overall characteristics and in
proportion with coal production. ( 4) Different vegetation types presented different changing characters caused by coal
http: / / www.ecologica.cn
mining.(5)The disturbance effect which large鄄scale coal mining of high efficiency and large production coal mines adopted
circulation techniques brought to vegetation was still inconspicuous, however, the long鄄term effects need to be further
studied.
Key Words: coal mine area; vegetation coverage; spatio鄄temporal effect; graphic analysis
煤炭开采在为国家建设提供大量优质煤炭的同时也严重破坏了矿区的生态环境。 根据相关研究,煤炭开
发活动对矿区生态环境的主要影响有:地表沉陷、土地沙漠化、土壤质量下降、地表及地下水污染、植被破坏、
生态系统退化、生物多样性丧失、景观破坏、农作物减产等[1]。 这些直接或间接影响一方面破坏了矿区植被
的生长环境,使其失去适合的生存空间;另一方面,逐渐裸露的地表又进一步加剧了矿区生态脆弱度。 因此,
植被长势与矿区生态环境之间存在着高度的双向制衡关系。 同时,以生态学视角出发,诸如水土保持、生态系
统调节、生物多样性保护、土壤质量改良、景观修复等矿区资源环境诸多因子都不同程度依赖于植被的长势和
繁茂程度。 因此,植被是考量矿区生态环境状况的关键之一。 在衡量植被繁茂程度的要素中,植被覆盖度又
称为投影盖度,是指植被植株冠层或叶面在地面的垂直投影面积占植被区总面积的比例[2鄄4],其范围分布在
0—1之间,反映了植被在水平方向上的密度情况,数值越大表明植被覆盖度越高。 国内外研究表明,作为重
要的气候、生态水文影响因子,植被覆盖度影响着区域大气圈、水圈、生物圈层间的各种物质转化和能量转移
过程[5]。 因此,研究其演化过程是检验矿区生态环境变化趋势的重要指标之一,对于评价煤炭开发对矿区环
境影响机理、实施切实有效的生态修复工程等具有重要的研究意义。 目前,国内外关于区域植被覆盖度的研
究主要集中在植被覆盖度与气候因子的关系[6鄄7]、植被覆盖度变化对区域环境影响[8鄄9]、以城市化等为代表的
人类活动对植被覆盖度的影响[10鄄11]等方面。 在煤矿区,煤炭开发及相关生产活动显然是影响植被的最主要
因子,将植被覆盖度应用到煤矿区环境问题研究主要包括:基于景观视角的矿区土地利用 \覆被变化及植被景
观破碎化程度研究[12];矿区植被系统恢复过程中的多样性变化[13];基于遥感过程的 CASA 模型和生态环境
状况指数测算矿区植被净初级生产力及生态环境状况变化趋势[14];矿区地下水与植被的相互关系探讨[15];
采用 NDVI与一元线性回归方法,分析植被覆盖和土地沙化的动态变化[16]等,以上研究为本文提供了思路上
的启迪。 就已有的研究而言,植被覆盖度的计算方法已较为成熟,多数文献对植被覆盖度进行分级,对比分析
多时相分级结果,得出区域植被覆盖度的全局趋势和变化规律并加以解释,这些基于分级的对比分析能够在
一定程度上表现区域植被演化特征,但是很难进一步提高定量表达的精度。 同时,全局趋势分析在研究局部
空间问题中也受到限制。 因此,本文综合运用以局部性和空间相关性为主的 GIS空间分析及数理统计方法研
究矿区植被覆盖度受煤炭开放影响的时空效应,以期更加深入地挖掘煤炭开采活动对矿区植被的影响机制。
1摇 数据来源与处理
研究选择山西大同矿区作为研究区域(图 1),主要基于以下几点考虑:其一,煤炭开发历史悠久,产量大,
环境累积效应和扰动量显著。 作为国家“十二五冶规划的 14 个大型煤炭基地之一———晋北煤炭基地主干企
业和国内特大型煤矿集团,同煤集团 2010、2011 年原煤产量连续突破亿吨,采煤塌陷区范围广,15 个主要井
田范围沉陷区面积达到 345.26 km2,平均每个井田周围 20 km2以上的土地都是沉陷区,区域环境受煤炭开发
扰动显著。 其二,地处生态环境脆弱的中西部地区,在全国具有典型性。 由于地处黄土高原边缘,地表湿陷性
黄土受采动影响,保水能力差,加之沉陷区地形复杂,地表沟谷纵横,起伏陡峭,土地养分退化,耕作条件破坏,
亟待复垦修复,是典型的生态脆弱型矿区。 因此,大同矿区区位独特,具有很强的典型性,对此进行深入研究,
对于全国其他矿区有普遍的参考和借鉴意义。 本文以空间关联测度的图谱视角来表达大同矿区植被覆盖度
演变的时空格局,探索煤炭开发对矿区生态环境的影响机制。
研究区自然环境、植被和煤炭开采的总体特征为:
大同地区属大陆性季风气候,干旱、半干旱地区。 矿区地貌以中低土石山群和黄土丘陵构成,相对高差
6307 摇 生摇 态摇 学摇 报摇 摇 摇 33卷摇
http: / / www.ecologica.cn
图 1摇 研究区区位和地形
Fig.1摇 Location and terrain of research area
450m,山体陡峭,树枝状冲沟极为发育,地形支离破碎,沟壑纵横。 全区土壤类型主要有山地栗钙土,淡栗钙
土和少量草甸土及盐潮土。 土质疏松,肥力贫乏,有机质含量少,抗冲力低。 矿区植被总体稀疏,多样性差,种
类贫乏,旱化特征明显,具有雁北干草原过渡地带特征,表现出个体生态与群落生态的高度统一。 据初步统
计,自然植被组成以温性落叶阔叶灌丛为主,天然植物共 18 科、35 属、57 种。 草本植物多为旱生、中生植物,
有华北落叶松、油松、山杨和桦;灌丛主要有沙棘、虎榛子、绣线菊等[17]。
大同煤田为双纪煤田,走向 NE—SW,上部为侏罗纪含煤地层,位于大同煤田的东北部,下部为石炭、二叠
纪含煤地层,除最北端局部地层外,几乎遍布整个大同煤田。 矿区自大规模开发以来的数十年,始终以开采埋
藏较浅的侏罗纪煤炭为主,石炭、二叠纪煤炭开发近几年刚刚起步,因此本文的研究范围限定于侏罗纪煤田界
限。 研究区内,共有 17对大型矿井,5对分布在北部云冈沟,主要开采侏罗纪煤炭;12 对位于南部口泉沟,其
中 10对开采侏罗纪煤炭,2 对开采石炭、二叠纪煤炭,井田边界如图 1 所示。 开采侏罗纪煤层的井田年开采
规模 106—5伊106 t,开采年限普遍在 50a以上;开采石炭、二叠纪煤层的井田开采规模 1.5伊107—2伊107 t,大规
模开采年限不足 5a。
以大同矿区 1999年至 2010 年同期(10 月)12 景 TM / ETM 遥感影像数据,使用 ERDAS 软件首先提取归
一化植被指数(NDVI),根据像元二分法原理[18鄄21],利用 ERDAS的建模工具 Spatial Modeler根据 NDVI计算出
1999—2010年 12个年份大同矿区植被覆盖度。
2摇 植被盖度时空效应的测度方法
已有的研究大多基于植被盖度的空间格局演变,从总体上揭示矿区植被覆盖度的全局演化过程,但是空
间效应可能是大尺度的趋势也可能是局部效应,前者一般称为“一阶效应冶 [22],描述的是某个参数的总体变
化性,“一阶效应冶忽略了局部性和空间相关性等因素,在应用中会带来偏差,需要引入“二阶效应冶的方法,
“二阶效应冶是由空间依赖性和空间异质性产生的。 表达的是空间上邻近位置上数值间的相互趋同或背离的
倾向。 受以煤炭开发为主的多种因素的共同作用,矿区植被覆盖度表现出明显的局部依赖性和异质性,因此,
空间数据分析的“二阶效应冶方法在矿区植被覆盖度研究中具有重要意义。
运用空间关联指数 Global Moran忆s I、Getis鄄Ord General G、Getis鄄Ord Gi*、Anselin Local Moran I来测度全局
和局域的空间聚簇特征[21],前两者是用于探测整个研究区的空间关联结构模式;后两者用于识别不同空间位
置上的高值簇与低值簇,即热点区(hot spots)与冷点区(cold spots)的空间分布。 空间关联指数表达式和原理
如下:
Global Moran忆s I I =移
n
i = 1
移
n
j = 1
X i - 軈( )X X j - 軈( )X / S2移
n
i = 1
移
n
j = 1
Wij (1)
7307摇 21期 摇 摇 摇 黄翌摇 等:煤炭开发对矿区植被扰动时空效应的图谱分析 ———以大同矿区为例 摇
http: / / www.ecologica.cn
式中, X i 为区域 i 的观测值, Wij 为空间权重矩阵,空间相邻为 1,不相邻为 0, S2 为移
n
i = 1
X i - 軈( )X 2 / n 。 对
Moran忆s I结果进行统计检验,采用 z 检验: Z( I) = 1 - E( I) / Var( I) ,式中 E( I) 为数学期望值, E( I) =
- 1 / n - 1 , Var( I) 为变异数, Var( I) = E( I2) - E( I) 2。
在给定显著性水平时,若 Moran忆s I显著为正,则表示覆盖水平较高(或较低)的区域在空间上显著集聚。
反之,若 Moran忆s I显著为负,则表明该区域与其周边地区的覆盖水平具有显著的空间差异。 仅当 Moran忆s I
接近期望值-1 / (n-1)时,观测值之间才相互独立,在空间上随机分布,此时满足传统区域差异度量方法所要
求的独立条件。
Getis鄄Ord General G G(d) =移
n
i = 1
移
n
j = 1
Wij(d)X iX j /移
n
i = 1
移
n
j = 1
X iX j (2)
式中, Wij(d) 为以距离规则定义的空间权重,同样空间范围相邻为 1,不相邻为 0, X i和 X j是是 i区域和 j区域
的观测值。 在空间不集聚的假设下,G(d)的期望值(G)为: E(G) = W / n(n - 1) ,在正态分布的条件下,G(d)
的统计检验值为: Z( I) = G - E(G) / Var(G) ;当 G(d)值高于 E(d)值,且 Z值显著时,检测区出现高值簇;
当 G(d)值低于 E(d)值,且 Z值显著时,检测区出现低值簇,当 G(d)趋近于 E(d)时,检测区变量呈现出随机
分布的特征。
Getis鄄Ord Gi* G*(d) =移
n
j = 1
Wij(d)X j /移
n
j = 1
X j (3)
为了便于解释和比较,对 G*(d) 进行标准化处理: Z(G*) = G*i - E(G*i ) / Var(G*i )
式中, E(G*i ) 和 Var(G*i ) 分别是 G*的数学期望和方差,Wij(d) 为空间权重,权重的计算方法如同 Getis鄄Ord
General G。 如果 Z(G*) 为正,且显著,表明位置 i 周围的值相对较高(高于均值),属高值空间集聚(热点
区);反之,如果 Z(G*) 为负且显著,则表明位置 i 周围的值相对较低(低于均值),属低值空间集聚(冷点
区)。
Anselin Local Moran I Ii = (xi - X)
_
/ S2i 移
n
j = 1,j屹i
w ij(xi - X
_
) (4)
式中, X i 为区域 i的观测值, X为相关观测值的均值, Wij 为空间权重矩阵,空间相邻为 1,不相邻为 0, S2i =
移
n
j = 1,j屹i
w ij / n - 1 - X2
_
,式中 n为观测值的数量。 对 Anselin Local Moran忆s I 结果进行统计检验,采用 Z 检验:
Z( I) = 1 - E( I) / Var( I) ,式中 E( I) 为数学期望值, Var( I) 为变异数。 Anselin Local Moran忆s I可区分具
有统计学上的显著性(0.05 水平)的高值 (HH) 聚类、低值 (LL) 聚类、由低值围绕的高值异常值 (HL) 聚类
以及由高值围绕的低值异常值 (LH)聚类。
3摇 结果分析与讨论
3.1摇 全局演变特征分析
利用 1999—2010年大同矿区植被覆盖度计算 Global Moran忆s I 和 General G 的估计值及其相关指标
(表 1)。
可以看出,1999—2010年间,Global Moran忆s I估计值全部为正,检验结果显著,且数值的大小总体稳定在
0.7—0.8之间,但呈现一定的波动性,2003年和 2008年尤为明显,表明 1999年以来植被覆盖度空间上分布的
趋势,即植被覆盖度相似(高或低)的地区在空间上呈现集中分布的特点,且总体上较稳定;在此期间,全局 G
统计指标的观测值和期望值都十分接近,相差不大,且都大于 0,当设定总体显著性水平 a 为 0.05 时,可知全
局 G统计量除 2000年,2003年,2009年 3个年份不显著外,其余各年份都显著,这说明检测区域高(或低)值
的集聚现象显著,植被盖度的变化可能围绕着热点区域来展开,但在不同的阶段呈现不同的集聚特点:
(1)1999—2002年间,Moran忆s I和 G(d)值分别由 0.765和 0.007782迅速扩大到 0.783和 0.007891,Z( I)
8307 摇 生摇 态摇 学摇 报摇 摇 摇 33卷摇
http: / / www.ecologica.cn
值也由 3.448扩大到 6.629,其中全局 G统计指标达到全期最大值。 表明矿区内植被盖度差异显著扩大,变异
系数由 0.245增加到 0.251。
表 1摇 植被盖度的变异系数、Moran忆s I和 General G估计值(1999—2010年)
Table 1摇 Estimated value of variation coefficient, Moran忆s I and General G of vegetation coverage
年份 植被盖度 变异系数 Moran忆s I G(d) Z( I) Z(d)
1999 0.4259 0.245 0.765 0.007782 3.448 3.433
2000 0.3852 0.189 0.816 0.007924 5.760 1.953
2001 0.3455 0.231 0.770 0.007980 4.509 3.610
2002 0.3307 0.251 0.783 0.007891 6.629 3.703
2003 0.3064 0.211 0.909 0.007879 4.920 1.457
2004 0.3742 0.131 0.707 0.007866 5.310 2.422
2005 0.3592 0.237 0.745 0.007814 8.029 4.272
2006 0.4146 0.149 0.759 0.007759 7.325 4.057
2007 0.4321 0.122 0.763 0.007813 6.585 3.134
2008 0.4200 0.088 0.634 0.007845 7.853 4.454
2009 0.465 0.117 0.756 0.007785 2.651 0.756
2010 0.5173 0.147 0.765 0.007748 7.027 4.307
(2)2003年,Moran忆s I值由 0.783扩大到 0.909,Z值分别由 6.629和 3.703缩小到 4.920和 1.457。 同时,
2003年也是全阶段矿区植被覆盖度最低的时期和转折点, Moran忆s I 值却达到全阶段最大值,矿区的植被受
破坏程度呈现高度正相关,植被覆盖度低的地区在空间上呈现集中分布。
(3)2004—2007年间是矿区植被覆盖度振荡上升的阶段,Moran忆s I和 G(d)值变化不大,Z( I)和 Z(d)值
却有不同程度的提升,均大于 0.05 显著水平上的标准阈值 1.96,表明 Moran忆s I 和 G(d)值在统计上显著性
增强。
(4)2008年启动的全矿区范围内的采煤棚户区搬迁工程对于矿区植被的扰动较明显。 采煤棚户区搬迁
过程中,涉及的人流、物流等要素的大规模流动现象,对矿区植被亦产生一定程度的破坏。 2008 年 Moran忆s I
值降至全阶段最低值,表明植被覆盖度在矿区范围内各处均有不同程度的下降。
(5)2009—2010年间,矿区植被覆盖度迅速提升,Moran忆s I 趋于稳定,并且缓慢增长,G(d)却缓慢下降,
并且 Z( I)和 Z(d)值呈现震荡的特点。 同时,棚户区搬迁后,由于无人耕种,矿区耕地大量减少,多数转化为
草地,在 10月份遥感图像中,农作物已收割,大同地区耕地表现为裸地类型,反映的植被盖度值较低,这也是
2009年以后矿区 10月植被盖度上升的原因之一。
3.2摇 局部效应分析与热点区域的演变
虽然 Global Moran忆s I和 Getis鄄Ord General G统计指标在一定程度上揭示了矿区植被盖度全局演变特征,
但当需要进一步考虑是否存在局部空间集聚时,全局空间自相关往往会掩盖局部状况或小范围的局部不稳定
性及不同位置上的空间变异程度,因此当全局关联特征 Global Moran忆s I和 Getis鄄Ord General G指标不能充分
揭示空间依赖性和异质性时,采用局域空间关联指数 Getis鄄Ord Gi*来探测局域空间的集聚程度,识别不同空
间位置上的高值簇与低值簇,即热点区与冷点区的空间分布[22]。 选取 2001—2010 年相对于 2000 年的植被
覆盖度变化情况作为局域统计指标 Getis鄄Ord Gi*来进一步进行空间关联分析,生成大同矿区植被覆盖度的空
间热点演变图(图 2)。
分析图 2得出:
(1)从整体来看,2000—2010年热点区域的整体格局基本保持稳定,但各种类型区的比例随着时间的推
移而有所变化(图 3),热点区、较热点区和冷点区、较冷点区的数量逐渐减少,次热点区和次冷点区的数量有
所增加,但个别年份变化很大,这说明,处于两极(冷点区和热点区)状态的类型区向中间状态(次热点区和次
9307摇 21期 摇 摇 摇 黄翌摇 等:煤炭开发对矿区植被扰动时空效应的图谱分析 ———以大同矿区为例 摇
http: / / www.ecologica.cn
图 2摇 植被覆盖度的演变热点图(相对于 2000年)
Fig.2摇 Hot spots mapping of vegetation coverage (related to year 2000)
由于缺乏同时相的 2007年、2009年大同矿区东部遥感数据,图中相应年份该处结果缺失
冷点区)分化,两极状态(热点区和冷点区)的类型区集聚效应越来越弱,极化现象越发不显著,总体上呈离散
的态势。
(2)从空间结构上看,2001—2004年,热点区向云冈沟周边地区不断集中,而冷点区逐渐形成了以东南部
地区为核心并在其周边地区集聚的势态,2005年以后云冈沟周边的热点区开始减退,又回到以两条沟为轴的
热点聚集区,并且分别为次热点区和次冷点区所包围,由此,形成了以两条沟为核心的阶梯状“核心鄄外围鄄边
缘冶结构,并且这种结构显得越来越稳固。 热点区集中在固定区域的现象较明显,而这两条沟又是矿区煤炭
生产、运输的主轴和居民生活的轴线,各种资源物质相互交织,植被覆盖度变化十分剧烈,这其中又以 2008 年
采煤棚户区搬迁年最为明显,大同矿区绝大部分的棚户区沿口泉沟和云冈沟布局,搬迁活动对沿线植被的扰
动量显著,2010年,其空间态势又开始减弱,热点区和冷点区范围较往年明显减小。
0407 摇 生摇 态摇 学摇 报摇 摇 摇 33卷摇
http: / / www.ecologica.cn
图 3摇 各种集聚类型区比例演化
Fig.3摇 Revolution of percentage for each centralize type
(3)从各类型区空间分布变化的格局来看,2001 年
热点区沿云冈沟分布,主要集中在燕子山矿附近,而冷
点区主要分布在云冈沟两侧,口泉沟及周边极化现象不
显著;2002—2006年,热点区从条带状演化为聚集在云
冈沟及两侧的巨大组团;值得注意的是,2008 年全国单
体规模最大的矿井———塔山高产矿井(年产量 2伊107 t)
大规模投产后,其所在的口泉沟南侧植被盖度与之前的
2004、2006 年相比变化不大,始终处于冷点区,但 2010
年开始有向热点区演化的趋势。
(4) 2001—2010 年,热点区域总体上聚集在矿区
的两条发展主轴口泉沟和云冈沟周围,而云冈沟的热点
区域范围和热度大于口泉沟,从 2001 年到 2008 年,云
冈沟内矿井的煤炭产量始终高于口泉沟,其植被受扰动
的剧烈程度也相应较高;2009 年,随着塔山、同忻等高
产矿井逐步稳产,口泉沟内煤炭产量超过了云冈沟,但
热点区并未产生明显的转移。
植被盖度热点区域的演变探测了各年份矿区植被
状况空间结构的局域变化,根据图谱分析,可将大同矿
区植被盖度演化分为 3 个典型的阶段:云冈沟两侧的热点集聚阶段(2001—2004 年)、沿轴线演化阶段
(2005—2008年)、热点消退阶段(2009—2010年),为了研究煤炭开发及相关活动对植被扰动过程中是否存
在着空间上的不一致性,分别选取 3 个阶段的典型年份 2001—2002 年、2005—2006 年、2009—2010 年,运用
Anselin Local Moran I对植被盖度的变化进行聚类分析(图 4)。
图 4摇 典型年份植被覆盖度变化聚类分析
Fig.4摇 Clustering analysis of vegetation coverage in typical years
3个时段各增长单元的集聚类型均发生了变化,2001—2002 年、2005—2006 年呈现明显的聚集效应,而
2009—2010年聚类效应总体减弱,但在口泉沟南侧却集聚明显。 全阶段聚类类型中,基本表现为“高鄄高冶聚
类和“低鄄低冶聚类,即植被盖度变化大的地区呈集中分布的特点,几乎未出现异常值,即“高鄄低冶聚类的情况,
表明在高强度煤炭开采扰动下,植被遭受的破坏是整体性的。
1407摇 21期 摇 摇 摇 黄翌摇 等:煤炭开发对矿区植被扰动时空效应的图谱分析 ———以大同矿区为例 摇
http: / / www.ecologica.cn
4摇 结论
本文综合运用以局部性和空间相关性为主的 GIS空间分析方法,研究植被这一综合体现矿区生态环境状
况的关键因子受以煤炭开发为主的矿区活动扰动的时空效应,对比以往大多单纯基于 NDVI 值的矿区植被演
化研究,Global Moran忆s I、Getis鄄Ord General G、Getis鄄Ord Gi*、Anselin Local Moran I等空间关联指数在趋势分析
等一阶整体效应的基础上,分别从全局演变和局部效应的视角揭示了矿区植被受扰动的时空演化和内在作用
机制。
(1)时序上总体把握及影响因子分析,1999 到 2010 年,大同矿区植被状况以 2003 年和 2006 年为分界
点,将植被盖度分为下降阶段、震荡上升阶段和上升阶段,在煤炭开采规模不断加大的过程中,矿区植被状况
总体好转,主要存在以下原因:淤2003年煤炭行业经济形势好转,煤矿企业开始重视矿区绿化,大量种植各类
树木,积极开展植被修复工程,这是矿区植被盖度总体上升的原因之一。 于近年来,煤矿绿色开采技术不断投
入实践,新工艺的采用对环境的扰动呈减小趋势。 盂同时,棚户区搬迁后大量耕地转化为野草地和低矮灌丛,
加之人类生活破坏的减少,客观上对植被盖度上升起到一定作用。 2003 年以前,煤炭开发和矿区人类活动是
植被受扰动较大且破坏明显的主要因素;2003—2008 年,矿区植被受多种因素共同影响,变化波动性较大;
2008年棚户区搬迁以后,种植活动和煤炭开发分别成为影响植被的正负因子。
(2)全局演变分析,从 Global Moran忆s I和 G(d)变化看,当矿区植被受破坏最严重时,在空间格局上,植被
盖度变化相似(高或低)的地区在空间上呈集聚分布状态,表现出很强的空间自相关性,全阶段矿区植被覆盖
度最低时,Moran忆s I值却达到全阶段最大值,矿区植被受破坏程度呈现高度正相关,而这种相关性在植被恢
复期呈减弱的趋势。
(3)局部效应分析,植被盖度变化的热点区不断减少或向次热点地区分散,矿区植被变化集聚效应越来
越弱,总体上呈离散的态势。 热点区集中在固定区域的现象较明显。 基本表现为“高鄄高冶聚类和“低鄄低冶聚
类,即植被盖度变化大的地区呈集中分布的特点,几乎未出现异常值,即“高鄄低冶聚类的情况,表明在高强度煤
炭开采扰动下,植被遭受的破坏是整体性的。
(4)从煤炭开发与植被受扰动相关性看,2001—2010 年,热点区域总体上聚集在矿区的两条发展主轴口
泉沟和云冈沟周围,而云冈沟的热点区域范围和热度大于口泉沟,从 2001年到 2008年,云冈沟内矿井的煤炭
产量始终高于口泉沟,其植被受扰动的剧烈程度也相应较高,表明煤炭开发对植被的破坏与煤炭产量呈正
相关。
(5)从不同植被类型受扰动效应看,大同矿区植被类型分异明显,山谷内种植有山杨和桦等木本植物,随
着海拔的升高,逐渐过渡到以草本和灌丛为主的植被类型,将植被类型与植被覆盖度结合分析,在煤炭开发活
动遍及研究区的前提下,对比人类活动均较少的相对高海拔区(矿区西南部地区)和相对低海拔区(云冈沟北
部地区),高海拔地区植被覆盖度变化较平缓,草本植物受扰动的程度应低于木本植物。 但是,矿区植被受煤
炭开发、土壤质量变化、区域气候等多种因素共同影响,势必存在一定的演替规律,在共同作用下的演替机制
是下一步研究的重点。
此外,口泉沟南侧一直是扰动效应较轻的区域,近年来投产的塔山和同忻矿井虽然是国内单产最大的井
工煤矿,但开采的是埋藏较深的石炭、二叠纪煤炭,对地面环境的影响具有较强的滞后效应,且采用了高效循
环经济工艺,力图使煤炭开发和产能对区域环境造成扰动减小到最低,从投产初期看,其对生态环境的扰动量
不明显,但煤炭开发对矿区环境的破坏存在累积性和滞后性,其长远扰动效应有待深入研究。
References:
[ 1 ]摇 Fan Y H, Lu Z H, Cheng J L,Zhou Z X,Wu G. Major ecological and environmental problems and the ecological reconstruction technologies of the
coal mining areas in China. Acta Ecologica Sinica, 2003, 23(10): 2144鄄2152.
[ 2 ] 摇 Gitelson A A, Kaufman Y J, Stark R, Rundquist D. Novel algorithms for remote estimation of vegetation fraction. Remote Sensing of Environment,
2002, 80(1): 76鄄87.
[ 3 ] 摇 Xing Z R, Feng Y G, Yang G J, Wang P, Huang W J. Method of estimating vegetation coverage based on remote sensing. Remote Sensing
2407 摇 生摇 态摇 学摇 报摇 摇 摇 33卷摇
http: / / www.ecologica.cn
Technology and Application, 2009, 24(6): 849鄄854.
[ 4 ] 摇 Garcia S, Saura S. Estimating landscape fragmentation indices from satellite images: the effect of sensor spatial resolution. Proc of the SPIE, 2004,
52(3): 668鄄675.
[ 5 ] 摇 Wimberly M C, Ohmann J L. A multi鄄scale assessment of human and environmental constraints on forest land cover change on the Oregon (USA)
coast range. Landscape Ecology, 2004, 19(6): 631鄄646.
[ 6 ] 摇 Liu L L, Xiao F J. Spatial鄄temporal correlations of NDVI with precipitation and temperature in Yellow River Basin. Chinese Journal of Ecology,
2006, 25(5): 477鄄481.
[ 7 ] 摇 Yang J P, Ding Y J, Chen R S. NDVI reflection of alpine vegetation changes in the south Yangtze and yellow rivers. Acta Geographica Sinica,
2005, 60(3): 467鄄478.
[ 8 ] 摇 Niu B R, Liu J R, Wang Z W. Remote sensing information extraction based on vegetation fraction in Drought and Half鄄Drought Area. Geomatics and
Information Science of Wuhan University, 2005, 30(1): 27鄄30.
[ 9 ] 摇 Gan C Y, Wang X Z, Li B S. Changes of vegetation coverage during recent 18 years in Lianjiang River Watershed. Scientia Geographica Sinica,
2011, 31(8): 1019鄄1024.
[10] 摇 L L, Ma A Q, Ma Q M. Spatial and temporal variations of vegetation coverage in Coastal Peri鄄urban Area: a case study of Laoshan District,
Qingdao. Environmental Science and Technology, 2002, 35(1): 178鄄185.
[11] 摇 Ji X M, Cheng T, Zheng F. Progress in ecological environment research in peri鄄urban areas. Ecology and Environmental Sciences, 2009, 18(4):
1579鄄1586.
[12] 摇 Hu Z Q, Xie H Q. Study on land use / cover change of coal mining area based on remote sensing images. Journal of China Coal Society, 2005, 30
(1): 44鄄48.
[13] 摇 Guo X Y, Zang J T, Gong H L, Zhang G L, Dong Z. Analysis of changes of the species diversity in the process of vegetation restoration in Antaibao
Mining Field, China. Acta Ecologica Sinica, 2005, 25(4): 763鄄770.
[14] 摇 Hou H P, Zhang S L, Ding Z Y, Gong Y L, Ma C Z. Study on the measurement of ecological loss based on the net primary productivity in coal
mines. Journal of China Coal Society, 2012, 37(3): 445鄄451.
[15] 摇 Wang L, Wei S P, Wang Q J. Effect of coal exploitation on groundwater and vegetation in the Yushenfu Coal Mine. Journal of China Coal Society,
2008, 33(12): 1408鄄1414.
[16] 摇 Wu L X, Ma B D, Liu S J. Analysis to vegetation coverage change in Shendong mining area with SPOT NDVI data. Journal of China Coal Society,
2009, 34(9): 1217鄄1222.
[17] 摇 Zhang L F.Main methods discussion of vegetation recovery in mining area of Datong city. Forestry of Shanxi, 2009, (3): 29鄄30.
[18] 摇 Zhao Y S. Principle and Methods of Remote Sensing Application and Analysis. Beijing: Science Press, 2003.
[19] 摇 Duncan J, Stow D, Franklin J, Hope A. Assessing the relation鄄ship between spectral vegetation indices and shrub cover in the Jornada Basin, New
Mexico. International Journal of Remote Sensing, 1993, 14(18): 3395鄄3416.
[20] 摇 Carlson T N, Ripley D A. On the Relation between NDVI, fractional vegetation cover, and leaf area index. Remote Sensing of Environment, 1997,
62(3): 241鄄252.
[21] 摇 Jiang M, Zang X F, S Q, Tong Q X. Vegetation coverage retrieval scale effect analysis using multi鄄sensor data. Geomatics and Information Science
of Wuhan University, 2011, 36(3): 311鄄315.
[22] 摇 Wang Y F, He H L. Methods of Spatial Data Analysis. Beijing: Science Press, 2007.
[23] 摇 Ma X D, Zhu C G, Ma R H, Pu Y X. Urban spatial growth pattern and its evolution in Suzhou, Jiangsu Province, China. Acta Geographica
Sinica, 2008, 63(4): 405鄄416.
参考文献:
[ 1 ]摇 范英宏, 陆兆华, 程建龙, 周忠轩, 吴钢. 中国煤矿区主要生态环境问题及生态重建技术. 生态学报, 2003, 23(10): 2144鄄2152.
[ 3 ] 摇 邢著荣, 冯幼贵, 杨贵军, 王萍, 黄文江. 基于遥感的植被覆盖度估算方法述评. 遥感技术与应用, 2009, 24(6): 849鄄854.
[ 6 ] 摇 刘绿柳, 肖风劲. 黄河流域植被 NDVI与温度、降水关系的时空变化. 生态学杂志, 2006, 25(5): 477鄄481.
[ 7 ] 摇 杨建平, 丁永健, 陈仁升. 长江黄河源区高寒植被变化的 NDVI记录. 地理学报, 2005, 60(3): 467鄄478.
[ 8 ] 摇 牛宝茹, 刘俊蓉, 王政伟. 干旱半干旱地区植被覆盖度遥感信息提取研究. 武汉大学学报: 信息科学版, 2005, 30(1): 27鄄30.
[ 9 ] 摇 甘春英, 王兮之, 李保生. 连江流域近 18年来植被覆盖度变化分析. 地理科学, 2011, 31(8): 1019鄄1024.
[10] 摇 刘林, 马安青, 马启敏. 滨海半城市化地区植被覆盖度的时空变化———以青岛市崂山区为例. 环境科学与技术, 2002, 35(1): 178鄄185.
[11] 摇 季小妹, 陈田, 郑芳. 半城市化地区生态环境研究进展. 生态环境学报, 2009, 18(4): 1579鄄1586.
[12] 摇 胡振琪, 谢宏全. 基于遥感图像的煤矿区土地利用 /覆盖变化. 煤炭学报, 2005, 30(1): 44鄄48.
[13] 摇 郭逍宇, 张金屯, 宫辉力, 张桂莲, 董志. 安太堡矿区复垦地植被恢复过程多样性变化. 生态学报, 2005, 25(4): 763鄄770.
[14] 摇 侯湖平, 张绍良, 丁忠义, 公云龙, 马昌忠. 基于植被净初级生产力的煤矿区生态损失测度研究. 煤炭学报, 2012, 37(3): 445鄄451.
[15] 摇 王力, 卫三平, 王全九. 榆神府煤田开采对地下水和植被的影响. 煤炭学报, 2008, 33(12): 1408鄄1414.
[16] 摇 吴立新, 马保东, 刘善军. 基于 SPOT卫星 NDVI数据的神东矿区植被覆盖动态变化分析. 煤炭学报, 2009, 34(9): 1217鄄1222.
[17] 摇 张利芳. 大同市矿区植被恢复主要途径探讨. 山西林业, 2009, (3): 29鄄30.
[18] 摇 赵英时. 遥感应用分析原理与方法. 北京: 科学出版社, 2003.
[21] 摇 江淼, 张显峰, 孙权, 童庆禧. 不同分辨率影像反演植被覆盖度的参数确定与尺度效应分析. 武汉大学学报: 信息科学版, 2011, 36(3):
311鄄315.
[22] 摇 王远飞, 何洪林. 空间数据分析方法. 北京: 科学出版社, 2007.
[23] 摇 马晓冬, 朱传耿, 马荣华, 蒲英霞. 苏州地区城镇扩展的空间格局及其演化分析. 地理学报, 2008, 63(4): 405鄄416.
3407摇 21期 摇 摇 摇 黄翌摇 等:煤炭开发对矿区植被扰动时空效应的图谱分析 ———以大同矿区为例 摇
粤悦栽粤 耘悦韵蕴韵郧陨悦粤 杂陨晕陨悦粤 灾燥造援猿猿袁晕燥援圆员 晕燥增援袁圆园员猿渊杂藻皂蚤皂燥灶贼澡造赠冤
悦韵晕栽耘晕栽杂
云则燥灶贼蚤藻则泽 葬灶凿 悦燥皂责则藻澡藻灶泽蚤增藻 砸藻增蚤藻憎
郧则燥泽泽 藻糟燥泽赠泽贼藻皂 责则燥凿怎糟贼院 贼澡藻燥则赠 枣则葬皂藻憎燥则噪 葬灶凿 糟葬泽藻 泽贼怎凿赠 韵哉再粤晕郧 在澡蚤赠怎灶袁在匀哉 悦澡怎灶择怎葬灶袁再粤晕郧 郧怎葬灶早遭蚤灶袁藻贼 葬造 渊远苑源苑冤噎噎噎
粤凿增葬灶糟藻泽 蚤灶 蚤皂责葬糟贼泽 燥枣 糟造蚤皂葬贼藻 糟澡葬灶早藻 燥灶 蚤灶枣藻糟贼蚤燥怎泽 凿蚤泽藻葬泽藻泽 燥怎贼遭则藻葬噪 蕴陨 郧怎燥凿燥灶早袁在匀粤晕郧 允怎灶澡怎葬袁允陨粤韵 郧藻灶早躁怎灶袁藻贼 葬造 渊远苑远圆冤噎噎
耘糟燥造燥早赠 燥枣 葬藻则燥遭蚤糟 皂藻贼澡葬灶藻 燥曾蚤凿蚤扎蚤灶早 遭葬糟贼藻则蚤葬 渊皂藻贼澡葬灶燥贼则燥责澡泽冤 再哉晕 允怎葬灶造蚤袁宰粤晕郧 再葬灶枣藻灶袁 在匀粤晕郧 匀燥灶早曾怎灶 渊远苑苑源冤噎噎噎噎噎
晕蚤贼则燥早藻灶 凿藻责燥泽蚤贼蚤燥灶 葬灶凿 蕴藻赠皂怎泽 糟澡蚤灶藻灶泽蚤泽 造藻葬枣 糟澡造燥则燥责澡赠造造 糟燥灶贼藻灶贼 蚤灶 陨灶灶藻则 酝燥灶早燥造蚤葬灶 早则葬泽泽造葬灶凿
在匀粤晕郧 再怎灶澡葬蚤袁 匀耘 晕蚤葬灶责藻灶早袁 在匀粤晕郧 郧怎葬灶早皂蚤灶早袁藻贼 葬造 渊远苑愿远冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
宰燥则造凿憎蚤凿藻 凿蚤泽贼则蚤遭怎贼蚤燥灶 葬灶凿 皂怎造贼蚤增葬则蚤葬贼藻 泽蚤皂蚤造葬则蚤贼赠 糟造怎泽贼藻则蚤灶早 葬灶葬造赠泽蚤泽 燥枣 泽责蚤凿藻则泽
杂匀耘晕 载蚤葬燥糟澡藻灶早袁在匀粤晕郧 月葬燥泽澡蚤袁 在匀粤晕郧 云藻灶早袁 藻贼 葬造 渊远苑怨缘冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
栽澡藻 蚤灶枣造怎藻灶糟藻 燥枣 憎蚤灶凿 凿蚤则藻糟贼蚤燥灶 燥灶 责燥造造藻灶鄄皂藻凿蚤葬贼藻凿 早藻灶藻 枣造燥憎 蚤灶 贼则葬灶泽早藻灶蚤糟 蚤灶泽藻糟贼鄄则藻泽蚤泽贼葬灶贼 糟燥贼贼燥灶
在匀哉 允蚤葬造蚤灶袁 匀耘 允怎葬灶袁 晕陨哉 允蚤葬灶择怎灶袁 藻贼 葬造 渊远愿园猿冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
粤怎贼藻糟燥造燥早赠 驭 云怎灶凿葬皂藻灶贼葬造泽
耘糟燥造燥早蚤糟葬造 糟澡葬则葬糟贼藻则蚤泽贼蚤糟泽 燥枣 遭藻灶贼澡蚤糟 凿蚤葬贼燥皂泽袁 责则燥贼燥扎燥葬 葬灶凿 皂藻蚤燥遭藻灶贼澡燥泽 蚤灶 贼澡藻 泽藻凿蚤皂藻灶贼泽 燥枣 贼澡藻 悦澡葬灶早躁蚤葬灶早 耘泽贼怎葬则赠 葬灶凿 耘葬泽贼 悦澡蚤灶葬
杂藻葬 蚤灶 泽责则蚤灶早 酝耘晕郧 在澡葬燥糟怎蚤袁 载哉 运怎蚤凿燥灶早 渊远愿员猿冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
酝葬糟则燥遭藻灶贼澡蚤糟 枣怎灶糟贼蚤燥灶葬造 早则燥怎责泽 葬贼 贼澡藻 则藻糟造葬皂葬贼蚤燥灶 葬灶凿 灶葬贼怎则葬造 贼蚤凿葬造 枣造葬贼泽 燥枣 匀藻灶早泽澡葬 耘葬泽贼 杂澡燥葬造袁 贼澡藻 耘泽贼怎葬则赠 燥枣 悦澡葬灶早躁蚤葬灶早 砸蚤增藻则
蕴灾 宰藻蚤憎藻蚤袁 酝粤 悦澡葬灶早忆葬灶袁 再哉 允蚤袁 藻贼 葬造 渊远愿圆缘冤
噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
耘灶则蚤糟澡皂藻灶贼 葬灶凿 藻糟燥造燥早蚤糟葬造 则蚤泽噪 燥枣 澡藻葬增赠 皂藻贼葬造 蚤灶 泽燥蚤造泽 葬灶凿 凿燥皂蚤灶葬灶贼 责造葬灶贼泽 蚤灶 贼澡藻 则蚤责葬则蚤葬灶 燥枣 贼澡藻 云藻灶早澡藻 砸蚤增藻则
再粤晕郧 再葬灶早袁 在匀韵哉 在澡藻灶早糟澡葬燥袁 宰粤晕郧 匀怎葬灶澡怎葬灶袁 藻贼 葬造 渊远愿猿源冤
噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
耘枣枣藻糟贼泽 燥枣 泽葬造蚤灶蚤贼赠 葬灶凿 藻曾燥早藻灶燥怎泽 泽怎遭泽贼则葬贼藻泽 燥灶 贼澡藻 凿藻糟燥皂责燥泽蚤贼蚤燥灶 葬灶凿 贼则葬灶泽枣燥则皂葬贼蚤燥灶 燥枣 泽燥蚤造 燥则早葬灶蚤糟 糟葬则遭燥灶 蚤灶 贼澡藻 再藻造造燥憎 砸蚤增藻则
阅藻造贼葬 蕴陨 蕴蚤灶早袁 匝陨哉 杂澡葬燥躁怎灶袁 栽粤晕 云藻蚤枣藻蚤袁 藻贼 葬造 渊远愿源源冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
耘枣枣藻糟贼泽 燥枣 泽澡燥则贼鄄贼藻则皂 凿葬则噪 糟澡蚤造造蚤灶早 燥灶 造藻葬增藻泽 糟葬则遭燥灶 葬灶凿 灶蚤贼则燥早藻灶 皂藻贼葬遭燥造蚤泽皂 葬灶凿 蚤灶增燥造增藻凿 葬糟贼蚤增蚤贼蚤藻泽 燥枣 藻灶扎赠皂藻泽 蚤灶 皂葬灶早则燥增藻 运葬灶凿藻造蚤葬
燥遭燥增葬贼葬 泽藻藻凿造蚤灶早 在匀耘晕郧 悦澡怎灶枣葬灶早袁 蕴陨哉 宰藻蚤糟澡藻灶早袁 悦匀耘晕 杂澡葬燥遭燥袁 藻贼 葬造 渊远愿缘猿冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
孕则藻造蚤皂蚤灶葬则赠 藻增葬造怎葬贼蚤燥灶 燥灶 贼燥造藻则葬灶糟藻 贼燥 责澡燥泽责澡燥则燥怎泽 凿藻枣蚤糟蚤藻灶糟赠 燥枣 猿圆 糟怎造贼蚤增葬则泽 燥枣 糟怎贼 糟澡则赠泽葬灶贼澡藻皂怎皂
蕴陨哉 孕藻灶早袁 悦匀耘晕 杂怎皂藻蚤袁 云粤晕郧 宰藻蚤皂蚤灶袁 藻贼 葬造 渊远愿远猿冤
噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
耘枣枣藻糟贼泽 燥枣 葬早藻 葬灶凿 藻灶增蚤则燥灶皂藻灶贼葬造 糟燥灶凿蚤贼蚤燥灶泽 燥灶 葬糟糟怎皂怎造葬贼蚤燥灶 燥枣 澡藻葬增赠鄄皂藻贼葬造泽 悦凿 葬灶凿 悦怎 蚤灶 栽藻早蚤造造葬则糟葬 早则葬灶燥泽葬
宰粤晕郧 在澡葬燥早藻灶袁宰哉 匀燥灶早曾蚤袁悦匀耘晕 载蚤葬燥曾蚤葬燥袁藻贼 葬造 渊远愿远怨冤
噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
耘枣枣藻糟贼泽 燥枣 悦澡蚤灶藻泽藻 早葬造造灶怎贼 燥灶 责澡燥贼燥泽赠灶贼澡藻贼蚤糟 糟澡葬则葬糟贼藻则蚤泽贼蚤糟泽 葬灶凿 贼燥贼葬造 灶蚤贼则燥早藻灶 糟燥灶贼藻灶贼 燥枣 砸澡怎泽 糟澡蚤灶藻灶泽蚤泽
蕴陨 再葬灶早袁 再粤晕郧 在蚤曾蚤葬灶早袁 悦匀耘晕 载蚤葬燥皂蚤灶早袁藻贼 葬造 渊远愿苑远冤
噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
栽澡藻 糟澡葬则葬糟贼藻则蚤扎葬贼蚤燥灶 燥枣 早造赠责澡燥泽葬贼藻 凿藻早则葬凿葬贼蚤燥灶 遭赠 月怎则噪澡燥造凿藻则蚤葬 皂怎造贼蚤增燥则葬灶泽 宰杂鄄云允怨 蕴陨 郧怎葬灶曾蚤袁 宰哉 载蚤葬燥择蚤灶袁再耘 允蚤葬灶则藻灶 渊远愿愿缘冤噎噎
耘造藻糟贼则燥葬灶贼藻灶灶燥早则葬责澡蚤糟 葬灶凿 遭藻澡葬增蚤燥怎则葬造 则藻泽责燥灶泽藻泽 燥枣 泽糟葬则葬遭 遭藻藻贼造藻泽 贼燥 砸蚤糟蚤灶怎泽 糟燥皂皂怎灶蚤泽 造藻葬枣 增燥造葬贼蚤造藻泽
蕴陨 宰藻蚤扎澡藻灶早袁 再粤晕郧 蕴藻蚤袁 杂匀耘晕 载蚤葬燥憎藻蚤袁 藻贼 葬造 渊远愿怨缘冤
噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
孕燥责怎造葬贼蚤燥灶袁 悦燥皂皂怎灶蚤贼赠 葬灶凿 耘糟燥泽赠泽贼藻皂
耘糟燥泽赠泽贼藻皂 澡藻葬造贼澡 葬泽泽藻泽泽皂藻灶贼 蚤灶 月葬蚤赠葬灶早凿蚤葬灶 蕴葬噪藻 载哉 云藻蚤袁 在匀粤韵 再葬灶憎藻蚤袁 再粤晕郧 在澡蚤枣藻灶早袁 藻贼 葬造 渊远怨园源冤噎噎噎噎噎噎噎噎噎噎噎
悦澡葬则葬糟贼藻则蚤泽贼蚤糟泽 燥枣 皂葬糟则燥遭藻灶贼澡蚤糟 糟燥皂皂怎灶蚤贼蚤藻泽 蚤灶 皂葬灶早则燥增藻 憎藻贼造葬灶凿泽 葬造燥灶早 贼澡藻 憎葬贼藻则憎葬赠泽 燥枣 晕燥则贼澡 匀藻扎澡燥怎袁 在澡怎澡葬蚤袁 杂燥怎贼澡 悦澡蚤灶葬
宰粤晕郧 匀怎蚤袁 在匀韵晕郧 杂澡葬灶袁 云粤晕郧 在澡葬灶择蚤葬灶早 渊远怨员猿冤
噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
栽澡藻 蚤灶贼藻则葬糟贼蚤燥灶 遭藻贼憎藻藻灶 糟燥皂责燥灶藻灶贼泽 燥枣 藻糟燥泽赠泽贼藻皂 则藻泽责蚤则葬贼蚤燥灶 蚤灶 贼赠责蚤糟葬造 枣燥则藻泽贼 葬灶凿 早则葬泽泽造葬灶凿 藻糟燥泽赠泽贼藻皂泽
在匀哉 载蚤葬灶躁蚤灶袁 再哉 郧怎蚤则怎蚤袁 宰粤晕郧 匝蚤怎枣藻灶早袁 藻贼 葬造 渊远怨圆缘冤
噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
耘枣枣藻糟贼泽 燥枣 贼澡蚤灶灶蚤灶早 燥灶 皂葬糟则燥 枣怎灶早蚤 葬灶凿 贼澡藻蚤则 则藻造葬贼蚤燥灶泽澡蚤责 憎蚤贼澡 造蚤贼贼藻则 凿藻糟燥皂责燥泽蚤贼蚤燥灶 蚤灶 孕蚤灶怎泽 贼葬遭怎造葬藻枣燥则皂蚤泽 责造葬灶贼葬贼蚤燥灶泽
悦匀耘晕 载蚤葬燥袁 月粤陨 杂澡怎造葬灶袁 蕴陨哉 再燥灶早袁 藻贼 葬造 渊远怨猿缘冤
噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
月藻贼葬 凿蚤增藻则泽蚤贼赠 燥枣 责造葬灶贼 糟燥皂皂怎灶蚤贼蚤藻泽 蚤灶 月葬蚤泽澡葬灶扎怎 晕葬贼怎则藻 砸藻泽藻则增藻 栽粤晕 杂澡葬灶泽澡葬灶袁 再耘 在澡藻灶造蚤灶袁 再哉粤晕 蕴蚤怎遭蚤灶袁 藻贼 葬造 渊远怨源源冤噎噎噎噎
耘枣枣藻糟贼 燥枣 韵曾赠贼藻贼则葬藻赠藻造蚤灶藻 渊韵栽悦冤 燥灶 贼澡藻 葬糟贼蚤增蚤贼蚤藻泽 燥枣 藻灶扎赠皂藻 葬灶凿 皂蚤糟则燥遭蚤葬造 糟燥皂皂怎灶蚤贼赠 皂藻贼葬遭燥造蚤糟 责则燥枣蚤造藻泽 蚤灶 糟燥皂责燥泽贼蚤灶早
悦匀耘晕 在澡蚤曾怎藻袁 郧哉 允蚤藻袁 郧粤韵 匀怎葬袁 藻贼 葬造 渊远怨缘苑冤
噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
蕴葬灶凿泽糟葬责藻袁 砸藻早蚤燥灶葬造 葬灶凿 郧造燥遭葬造 耘糟燥造燥早赠
灾葬则蚤葬贼蚤燥灶泽 蚤灶 造藻葬枣 葬灶葬贼燥皂赠 燥枣 蕴葬则蚤曾 早皂藻造蚤灶蚤蚤 则藻枣造藻糟贼 葬凿葬责贼葬贼蚤燥灶 燥枣 蚤贼泽 责澡燥贼燥泽赠灶贼澡藻贼蚤糟 糟葬责葬糟蚤贼赠 贼燥 糟造蚤皂葬贼藻 糟澡葬灶早藻泽
允陨 在蚤躁蚤灶早袁 匝哉粤晕 载蚤葬灶噪怎蚤袁 宰粤晕郧 悦澡怎葬灶噪怎葬灶 渊远怨远苑冤
噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
栽澡藻 泽贼怎凿赠蚤灶早 燥枣 噪藻赠 藻糟燥造燥早蚤糟葬造 枣葬糟贼燥则泽 葬灶凿 贼澡则藻泽澡燥造凿 燥枣 造葬灶凿泽糟葬责藻 藻增燥造怎贼蚤燥灶 蚤灶 再葬灶糟澡藻灶早 悦燥葬泽贼葬造 憎藻贼造葬灶凿
在匀粤晕郧 匀怎葬遭蚤灶早袁 蕴陨哉 匀燥灶早赠怎袁 蕴陨 再怎枣藻灶早袁藻贼 葬造 渊远怨苑缘冤
噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
耘糟燥鄄责澡赠泽蚤燥造燥早蚤糟葬造 则藻泽责燥灶泽藻 燥枣 孕澡则葬早皂蚤贼藻泽 糟燥皂皂怎灶蚤泽 贼燥 憎葬贼藻则 贼葬遭造藻 糟澡葬灶早藻泽 蚤灶 贼澡藻 匀燥则择蚤灶 杂葬灶凿 蕴葬灶凿
酝粤 再怎灶澡怎葬袁在匀粤晕郧 栽燥灶早澡怎蚤袁 蕴陨哉 载蚤灶责蚤灶早 渊远怨愿源冤
噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
陨皂责则燥增藻皂藻灶贼 燥枣 泽灶燥憎皂藻造贼 蚤皂责造藻皂藻灶贼葬贼蚤燥灶 蚤灶 贼澡藻 杂宰粤栽 澡赠凿则燥造燥早蚤糟 皂燥凿藻造 再哉 宰藻灶躁怎灶袁晕粤晕 在澡怎燥贼燥灶早袁在匀粤韵 再葬灶遭燥袁藻贼 葬造 渊远怨怨圆冤噎
砸藻泽责燥灶泽藻泽 燥枣 造葬噪藻 枣造怎糟贼怎葬贼蚤燥灶 贼燥 糟造蚤皂葬贼藻 糟澡葬灶早藻 蚤灶 匀燥则择蚤灶 杂葬灶凿赠 蕴葬灶凿
悦匀粤晕郧 载怎藻造蚤袁 在匀粤韵 载怎藻赠燥灶早袁 宰粤晕郧 宰藻蚤袁 藻贼 葬造 渊苑园园圆冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
宰葬贼藻则 藻糟燥造燥早赠 葬灶凿 枣则葬糟贼葬造 糟澡葬则葬糟贼藻则蚤泽贼蚤糟泽 燥枣 泽燥蚤造 责葬则贼蚤糟造藻 泽蚤扎藻 凿蚤泽贼则蚤遭怎贼蚤燥灶 燥枣 贼澡则藻藻 贼赠责蚤糟葬造 增藻早藻贼葬贼蚤燥灶泽 蚤灶 杂澡藻造造 陨泽造葬灶凿
载陨粤 允蚤葬灶早遭葬燥袁 在匀粤晕郧 杂澡怎赠燥灶早袁 宰粤晕郧 砸燥灶早则燥灶早袁 藻贼 葬造 渊苑园员猿冤
噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
杂责葬贼蚤燥鄄贼藻皂责燥则葬造 凿蚤泽贼则蚤遭怎贼蚤燥灶 燥枣 藻责蚤造蚤贼澡蚤糟 葬造早葬造 糟澡造燥则燥责澡赠造造 葬 蚤灶 则藻造葬贼蚤燥灶 贼燥 贼澡藻 责澡赠泽蚤糟燥鄄糟澡藻皂蚤糟葬造 枣葬糟贼燥则泽 燥枣 郧怎枣怎 砸蚤增藻则 蚤灶 栽澡则藻藻 郧燥则鄄
早藻泽 砸藻泽藻则增燥蚤则 宰哉 杂澡怎赠怎葬灶袁 郧耘 允蚤憎藻灶袁 酝陨粤韵 宰藻灶躁蚤藻袁 藻贼 葬造 渊苑园圆猿冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
砸藻泽燥怎则糟藻 葬灶凿 陨灶凿怎泽贼则蚤葬造 耘糟燥造燥早赠
郧则葬责澡蚤糟 葬灶葬造赠泽蚤泽 燥枣 泽责葬贼蚤燥鄄贼藻皂责燥则葬造 藻枣枣藻糟贼 枣燥则 增藻早藻贼葬贼蚤燥灶 凿蚤泽贼怎则遭葬灶糟藻 糟葬怎泽藻凿 遭赠 糟燥葬造 皂蚤灶蚤灶早院 葬 糟葬泽藻 燥枣 阅葬贼燥灶早 悦燥葬造 酝蚤灶藻 粤则藻葬
匀哉粤晕郧 再蚤袁 宰粤晕郧 再怎灶躁蚤葬袁 蕴陨 载蚤葬燥泽澡怎灶袁 藻贼 葬造 渊苑园猿缘冤
噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
远源园苑 摇 生摇 态摇 学摇 报摇 摇 摇 猿猿卷摇
叶生态学报曳圆园员源年征订启事
叶生态学报曳是由中国科学技术协会主管袁中国生态学学会尧中国科学院生态环境研究中心主办的生态学
高级专业学术期刊袁创刊于 员怨愿员年袁报道生态学领域前沿理论和原始创新性研究成果遥 坚持野百花齐放袁百家
争鸣冶的方针袁依靠和团结广大生态学科研工作者袁探索生态学奥秘袁为生态学基础理论研究搭建交流平台袁
促进生态学研究深入发展袁为我国培养和造就生态学科研人才和知识创新服务尧为国民经济建设和发展服务遥
叶生态学报曳主要报道生态学及各分支学科的重要基础理论和应用研究的原始创新性科研成果遥 特别欢
迎能反映现代生态学发展方向的优秀综述性文章曰研究简报曰生态学新理论尧新方法尧新技术介绍曰新书评价和
学术尧科研动态及开放实验室介绍等遥
叶生态学报曳为半月刊袁大 员远开本袁圆愿园页袁国内定价 怨园元 辕册袁全年定价 圆员远园元遥
国内邮发代号院愿圆鄄苑袁国外邮发代号院酝远苑园
标准刊号院陨杂杂晕 员园园园鄄园怨猿猿摇 摇 悦晕 员员鄄圆园猿员 辕 匝
全国各地邮局均可订阅袁也可直接与编辑部联系购买遥 欢迎广大科技工作者尧科研单位尧高等院校尧图书
馆等订阅遥
通讯地址院 员园园园愿缘 北京海淀区双清路 员愿号摇 电摇 摇 话院 渊园员园冤远圆怨源员园怨怨曰 远圆愿源猿猿远圆
耘鄄皂葬蚤造院 泽澡藻灶早贼葬蚤曾怎藻遭葬燥岳 则糟藻藻泽援葬糟援糟灶摇 网摇 摇 址院 憎憎憎援藻糟燥造燥早蚤糟葬援糟灶
本期责任副主编摇 余新晓摇 摇 摇 编辑部主任摇 孔红梅摇 摇 摇 执行编辑摇 刘天星摇 段摇 靖
生摇 态摇 学摇 报渊杂匀耘晕郧栽粤陨摇 载哉耘月粤韵冤渊半月刊摇 员怨愿员年 猿月创刊冤
第 猿猿卷摇 第 圆员期摇 渊圆园员猿年 员员月冤
粤悦栽粤 耘悦韵蕴韵郧陨悦粤 杂陨晕陨悦粤摇渊杂藻皂蚤皂燥灶贼澡造赠袁杂贼葬则贼藻凿 蚤灶 员怨愿员冤摇灾燥造郾 猿猿摇 晕燥郾 圆员 渊晕燥增藻皂遭藻则袁 圆园员猿冤
编摇 摇 辑摇 叶生态学报曳编辑部
地址院北京海淀区双清路 员愿号
邮政编码院员园园园愿缘
电话院渊园员园冤远圆怨源员园怨怨憎憎憎援藻糟燥造燥早蚤糟葬援糟灶泽澡藻灶早贼葬蚤曾怎藻遭葬燥岳 则糟藻藻泽援葬糟援糟灶
主摇 摇 编摇 王如松
主摇 摇 管摇 中国科学技术协会
主摇 摇 办摇 中国生态学学会
中国科学院生态环境研究中心
地址院北京海淀区双清路 员愿号
邮政编码院员园园园愿缘
出摇 摇 版摇
摇 摇 摇 摇 摇 地址院北京东黄城根北街 员远号
邮政编码院员园园苑员苑
印摇 摇 刷摇 北京北林印刷厂
发 行摇
地址院东黄城根北街 员远号
邮政编码院员园园苑员苑
电话院渊园员园冤远源园猿源缘远猿耘鄄皂葬蚤造院躁燥怎则灶葬造岳 糟泽责早援灶藻贼
订摇 摇 购摇 全国各地邮局
国外发行摇 中国国际图书贸易总公司
地址院北京 猿怨怨信箱
邮政编码院员园园园源源
广告经营
许 可 证摇 京海工商广字第 愿园员猿号
耘凿蚤贼藻凿 遭赠摇 耘凿蚤贼燥则蚤葬造 遭燥葬则凿 燥枣
粤悦栽粤 耘悦韵蕴韵郧陨悦粤 杂陨晕陨悦粤
粤凿凿院员愿袁杂澡怎葬灶早择蚤灶早 杂贼则藻藻贼袁匀葬蚤凿蚤葬灶袁月藻蚤躁蚤灶早 员园园园愿缘袁悦澡蚤灶葬
栽藻造院渊园员园冤远圆怨源员园怨怨
憎憎憎援藻糟燥造燥早蚤糟葬援糟灶
泽澡藻灶早贼葬蚤曾怎藻遭葬燥岳 则糟藻藻泽援葬糟援糟灶
耘凿蚤贼燥则鄄蚤灶鄄糟澡蚤藻枣摇 宰粤晕郧 砸怎泽燥灶早
杂怎责藻则增蚤泽藻凿 遭赠摇 悦澡蚤灶葬 粤泽泽燥糟蚤葬贼蚤燥灶 枣燥则 杂糟蚤藻灶糟藻 葬灶凿 栽藻糟澡灶燥造燥早赠
杂责燥灶泽燥则藻凿 遭赠摇 耘糟燥造燥早蚤糟葬造 杂燥糟蚤藻贼赠 燥枣 悦澡蚤灶葬
砸藻泽藻葬则糟澡 悦藻灶贼藻则 枣燥则 耘糟燥鄄藻灶增蚤则燥灶皂藻灶贼葬造 杂糟蚤藻灶糟藻泽袁 悦粤杂
粤凿凿院员愿袁杂澡怎葬灶早择蚤灶早 杂贼则藻藻贼袁匀葬蚤凿蚤葬灶袁月藻蚤躁蚤灶早 员园园园愿缘袁悦澡蚤灶葬
孕怎遭造蚤泽澡藻凿 遭赠摇 杂糟蚤藻灶糟藻 孕则藻泽泽
粤凿凿院员远 阅燥灶早澡怎葬灶早糟澡藻灶早早藻灶 晕燥则贼澡 杂贼则藻藻贼袁
月藻蚤躁蚤灶早摇 员园园苑员苑袁悦澡蚤灶葬
孕则蚤灶贼藻凿 遭赠摇 月藻蚤躁蚤灶早 月藻蚤 蕴蚤灶 孕则蚤灶贼蚤灶早 匀燥怎泽藻袁
月藻蚤躁蚤灶早 员园园园愿猿袁悦澡蚤灶葬
阅蚤泽贼则蚤遭怎贼藻凿 遭赠摇 杂糟蚤藻灶糟藻 孕则藻泽泽
粤凿凿院员远 阅燥灶早澡怎葬灶早糟澡藻灶早早藻灶 晕燥则贼澡
杂贼则藻藻贼袁月藻蚤躁蚤灶早 员园园苑员苑袁悦澡蚤灶葬
栽藻造院渊园员园冤远源园猿源缘远猿
耘鄄皂葬蚤造院躁燥怎则灶葬造岳 糟泽责早援灶藻贼
阅燥皂藻泽贼蚤糟 摇 摇 粤造造 蕴燥糟葬造 孕燥泽贼 韵枣枣蚤糟藻泽 蚤灶 悦澡蚤灶葬
云燥则藻蚤早灶 摇 摇 悦澡蚤灶葬 陨灶贼藻则灶葬贼蚤燥灶葬造 月燥燥噪 栽则葬凿蚤灶早
悦燥则责燥则葬贼蚤燥灶
粤凿凿院孕援韵援月燥曾 猿怨怨 月藻蚤躁蚤灶早 员园园园源源袁悦澡蚤灶葬
摇 陨杂杂晕 员园园园鄄园怨猿猿悦晕 员员鄄圆园猿员 辕 匝 国内外公开发行 国内邮发代号 愿圆鄄苑 国外发行代号 酝远苑园 定价 怨园郾 园园元摇