全 文 ：
摇 摇 摇 摇 摇 生 态 学 报
摇 摇 摇 摇 摇 摇 摇 渊杂匀耘晕郧栽粤陨 载哉耘月粤韵冤
摇 摇 第 猿猿卷 第 员怨期摇 摇 圆园员猿年 员园月摇 渊半月刊冤
目摇 摇 次
中国生态学学会 圆园员猿年学术年会专辑摇 卷首语
生态系统服务研究文献现状及不同研究方向评述 马凤娇袁刘金铜袁粤援 耘早则蚤灶赠葬 耘灶藻躁蚤 渊缘怨远猿冤噎噎噎噎噎噎噎
非人灵长类性打搅行为研究进展 杨摇 斌袁王程亮袁纪维红袁等 渊缘怨苑猿冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
密度制约效应对啮齿动物繁殖的影响 韩群花袁郭摇 聪袁张美文 渊缘怨愿员冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
食物链长度远因与近因研究进展综述 王玉玉袁徐摇 军袁雷光春 渊缘怨怨园冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
粤酝真菌在植物病虫害生物防治中的作用机制 罗巧玉袁王晓娟袁李媛媛袁等 渊缘怨怨苑冤噎噎噎噎噎噎噎噎噎噎噎
保护性耕作对农田碳尧氮效应的影响研究进展 薛建福袁赵摇 鑫袁杂澡葬凿则葬糟噪 月葬贼泽蚤造藻 阅蚤噪早憎葬贼造澡藻袁等 渊远园园远冤噎噎噎
圈养大熊猫野化培训期的生境选择特征 张明春袁黄摇 炎袁李德生袁等 渊远园员源冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎
利用红外照相技术分析野生白冠长尾雉活动节律及时间分配 赵玉泽袁王志臣袁徐基良 袁等 渊远园圆员冤噎噎噎噎
风速和持续时间对树麻雀能量收支的影响 杨志宏袁吴庆明袁董海燕袁等 渊远园圆愿冤噎噎噎噎噎噎噎噎噎噎噎噎噎
白马雪山自然保护区灰头小鼯鼠的巢址特征 李艳红袁关进科袁黎大勇袁等 渊远园猿缘冤噎噎噎噎噎噎噎噎噎噎噎噎
生境片段化对千岛湖岛屿上黄足厚结猛蚁遗传多样性的影响 罗媛媛袁刘金亮袁黄杰灵袁等 渊远园源员冤噎噎噎噎噎
基于 圆愿杂袁 悦韵陨和 悦赠贼遭基因序列的薜荔和爱玉子传粉小蜂分子遗传关系研究
吴文珊袁陈友铃袁孙伶俐袁等 渊远园源怨冤
噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
高榕榕果内 耘怎责则蚤泽贼蚤灶葬属两种榕小蜂的遗传进化关系 陈友铃袁孙伶俐袁武蕾蕾袁等 渊远园缘愿冤噎噎噎噎噎噎噎噎
镉胁迫下杞柳对金属元素的吸收及其根系形态构型特征 王树凤袁施翔袁孙海菁袁等 渊远园远缘冤噎噎噎噎噎噎噎噎
邻苯二甲酸对萝卜种子萌发尧幼苗叶片膜脂过氧化及渗透调节物质的影响
杨延杰袁王晓伟袁赵摇 康袁等 渊远园苑源冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
极端干旱区多枝柽柳幼苗对人工水分干扰的形态及生理响应 马晓东袁王明慧袁李卫红袁等 渊远园愿员冤噎噎噎噎噎
贝壳砂生境酸枣叶片光合生理参数的水分响应特征 王荣荣袁夏江宝袁杨吉华袁等 渊远园愿愿冤噎噎噎噎噎噎噎噎噎
陶粒覆盖对土壤水分尧植物光合作用及生长状况的影响 谭雪红袁郭小平袁赵廷宁 渊远园怨苑冤噎噎噎噎噎噎噎噎噎
不同林龄短枝木麻黄小枝单宁含量及养分再吸收动态 叶功富袁张尚炬袁张立华袁等 渊远员园苑冤噎噎噎噎噎噎噎噎
珠江三角洲不同污染梯度下森林优势种叶片和枝条 杂含量比较 裴男才袁陈步峰袁邹志谨袁等 渊远员员源冤噎噎噎噎
粤酝真菌和磷对小马安羊蹄甲幼苗生长的影响 宋成军袁曲来叶袁马克明袁等 渊远员圆员冤噎噎噎噎噎噎噎噎噎噎噎
盐氮处理下盐地碱蓬种子成熟过程中的离子积累和种子萌发特性 周家超袁付婷婷袁赵维维袁等 渊远员圆怨冤噎噎噎
悦韵圆浓度升高条件下内生真菌感染对宿主植物的生理生态影响 师志冰袁周摇 勇袁李摇 夏袁等 渊远员猿缘冤噎噎噎噎
预处理方式对香蒲和芦苇种子萌发的影响 孟摇 焕袁王雪宏袁佟守正袁等 渊远员源圆冤噎噎噎噎噎噎噎噎噎噎噎噎噎
镉在土壤鄄金丝垂柳系统中的迁移特征 张摇 雯袁魏摇 虹袁孙晓灿袁等 渊远员源苑冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎
马尾松人工林近自然化改造对植物自然更新及物种多样性的影响 罗应华袁孙冬婧袁林建勇袁等 渊远员缘源冤噎噎噎
濒危海草贝克喜盐草的种群动态及土壤种子库要要要以广西珍珠湾为例
邱广龙袁范航清袁李宗善袁等 渊远员远猿冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
毛乌素沙地南缘沙丘生物结皮对凝结水形成和蒸发的影响 尹瑞平袁吴永胜袁张摇 欣袁等 渊远员苑猿冤噎噎噎噎噎噎
塔里木河上游灰胡杨种群生活史特征与空间分布格局 韩摇 路袁席琳乔袁王家强袁等 渊远员愿员冤噎噎噎噎噎噎噎噎
短期氮素添加和模拟放牧对青藏高原高寒草甸生态系统呼吸的影响 宗摇 宁袁石培礼袁蔣摇 婧袁等 渊远员怨员冤噎噎
松嫩平原微地形下土壤水盐与植物群落分布的关系 杨摇 帆袁王志春袁王云贺袁等 渊远圆园圆冤噎噎噎噎噎噎噎噎噎
广州大夫山雨季林内外空气 栽杂孕 和 孕酝圆援缘浓度及水溶性离子特征 肖以华袁李摇 炯袁旷远文袁等 渊远圆园怨冤噎噎噎
马鞍列岛岩礁生境鱼类群落结构时空格局 汪振华袁赵摇 静袁王摇 凯袁等 渊远圆员愿冤噎噎噎噎噎噎噎噎噎噎噎噎噎
黄海细纹狮子鱼种群特征的年际变化 陈云龙袁单秀娟袁周志鹏袁等 渊远圆圆苑冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
三种温带森林大型土壤动物群落结构的时空动态 李摇 娜袁张雪萍袁张利敏 渊远圆猿远冤噎噎噎噎噎噎噎噎噎噎噎
笔管榕榕小蜂的群落结构与物种多样性 陈友铃袁陈晓倩袁吴文珊袁等 渊远圆源远冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎
海洋生态资本理论框架下的生态系统服务评估 陈摇 尚袁任大川袁夏摇 涛袁等 渊远圆缘源冤噎噎噎噎噎噎噎噎噎噎噎
中国地貌区划系统要要要以自然保护区体系建设为目标 郭子良袁崔国发 渊远圆远源冤噎噎噎噎噎噎噎噎噎噎噎噎噎
生态植被建设对黄土高原农林复合流域景观格局的影响 易摇 扬袁信忠保袁覃云斌袁等 渊远圆苑苑冤噎噎噎噎噎噎噎
华北农牧交错带农田鄄草地景观镶嵌体土壤水分空间异质性 王红梅袁王仲良袁王摇 堃袁等 渊远圆愿苑冤噎噎噎噎噎
中国北方春小麦生育期变化的区域差异性与气候适应性 俄有浩袁霍治国袁马玉平袁等 渊远圆怨缘冤噎噎噎噎噎噎噎
中国南方喀斯特石漠化演替过程中土壤理化性质的响应 盛茂银袁刘摇 洋袁熊康宁 渊远猿园猿冤噎噎噎噎噎噎噎噎
气候变化对东北沼泽湿地潜在分布的影响 贺摇 伟袁布仁仓袁刘宏娟袁等 渊远猿员源冤噎噎噎噎噎噎噎噎噎噎噎噎噎
内蒙古不同类型草地土壤氮矿化及其温度敏感性 朱剑兴袁王秋凤袁何念鹏袁等 渊远猿圆园冤噎噎噎噎噎噎噎噎噎噎
黑河中游荒漠绿洲区土地利用的土壤养分效应 马志敏袁吕一河袁孙飞翔袁等 渊远猿圆愿冤噎噎噎噎噎噎噎噎噎噎噎
成都平原北部水稻土重金属含量状况及其潜在生态风险评价 秦鱼生袁喻摇 华袁冯文强袁等 渊远猿猿缘冤噎噎噎噎噎
大西洋中部延绳钓黄鳍金枪鱼渔场时空分布与温跃层的关系 杨胜龙袁马军杰袁张摇 禹袁等 渊远猿源缘冤噎噎噎噎噎
夏季台湾海峡南部海域上层水体的生物固氮作用 林摇 峰袁陈摇 敏袁杨伟锋袁等 渊远猿缘源冤噎噎噎噎噎噎噎噎噎噎
北长山岛森林乔木层碳储量及其影响因子 石洪华袁王晓丽袁王摇 嫒袁等 渊远猿远猿冤噎噎噎噎噎噎噎噎噎噎噎噎噎
植被类型变化对长白山森林土壤碳矿化及其温度敏感性的影响 王摇 丹袁吕瑜良袁徐摇 丽袁等 渊远猿苑猿冤噎噎噎噎
油松遗传结构与地理阻隔因素的相关性 孟翔翔袁狄晓艳袁王孟本袁等 渊远猿愿圆冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎
基于辅助环境变量的土壤有机碳空间插值要要要以黄土丘陵区小流域为例
文摇 雯袁周宝同袁汪亚峰袁等 渊远猿愿怨冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
基于生命周期视角的产业资源生态管理效益分析要要要以虚拟共生网络系统为例
施晓清袁李笑诺袁杨建新 渊远猿怨愿冤
噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
生态脆弱区贫困与生态环境的博弈分析 祁新华袁叶士琳袁程摇 煜袁等 渊远源员员冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎
野世博冶背景下上海经济与环境的耦合演化 倪摇 尧袁岳文泽袁张云堂袁等 渊远源员愿冤噎噎噎噎噎噎噎噎噎噎噎噎
期刊基本参数院悦晕 员员鄄圆园猿员 辕 匝鄢员怨愿员鄢皂鄢员远鄢源远源鄢扎澡鄢孕鄢 预 怨园郾 园园鄢员缘员园鄢缘缘鄢圆园员猿鄄员园
室室室室室室室室室室室室室室
封面图说院 毛乌素沙地南缘沙丘的生物结皮要要要生物土壤结皮广泛分布于干旱和半干旱区袁它的形成和发育对荒漠生态系统
生态修复过程产生重要的影响遥 组成生物结皮的藻类尧苔藓和地衣是常见的先锋植物袁它们不仅能在严重干旱缺
水尧营养贫瘠恶劣的环境中生长尧繁殖袁并且能通过其代谢方式影响并改变环境遥 其中一个重要的特点是袁生物结皮
表面的凝结水显著大于裸沙遥 研究表明袁凝结水是除降雨之外最重要的水分来源之一袁在水分极度匮乏的荒漠生态
系统袁它对荒漠生态系统结构尧功能和过程的维持产生着重要的影响遥
彩图及图说提供院 陈建伟教授摇 北京林业大学摇 耘鄄皂葬蚤造院 糟蚤贼藻泽援糟澡藻灶躁憎岳 员远猿援糟燥皂
第 33 卷第 19 期
2013年 10月
生 态 学 报
ACTA ECOLOGICA SINICA
Vol.33,No.19
Oct.,2013
http: / / www.ecologica.cn
基金项目:国家自然科学基金资助项目(40925003)
收稿日期:2013鄄05鄄03; 摇 摇 修订日期:2013鄄07鄄25
*通讯作者 Corresponding author.E鄄mail: yfwang@ rcees.ac.cn
DOI: 10.5846 / stxb201305030916
文雯,周宝同,汪亚峰,黄勇.基于辅助环境变量的土壤有机碳空间插值———以黄土丘陵区小流域为例.生态学报,2013,33(19):6389鄄6397.
Wen W, Zhou B T, Wang Y F, Huang Y.Soil organic carbon interpolation based on auxiliary environmental covariates:a case study at small watershed scale
in Loess Hilly region.Acta Ecologica Sinica,2013,33(19):6389鄄6397.
基于辅助环境变量的土壤有机碳空间插值
———以黄土丘陵区小流域为例
文摇 雯1,2,周宝同1,汪亚峰2,*,黄摇 勇2,3
(1. 西南大学地理科学学院,重庆摇 400715; 2. 中国科学院生态环境研究中心城市与区域生态国家重点实验室,北京摇 100085;
3. 东北林业大学林学院,哈尔滨摇 150040 )
摘要:利用普通克里格法(OK)、反距离加权法(IDW)、径向基函数法(RBF)、基于土地利用类型修正的普通克里格法(OK_LU)
4种插值方法,对黄土丘陵羊圈沟小流域的土壤有机碳含量进行空间插值。 预测结果的准确性通过 Pearson相关系数(R),平均
绝对误差(MAE),均方根误差(RMSE),准确度(AC)来评价。 研究结果表明:(1)在前 3 种常规空间插值方法中,OK 对刻画区
域土壤有机碳的空间分布趋势效果最佳,其预测 MAE值和 RMSE值均为最小,Pearson相关系数(R)和准确度(AC)最大,说明
其预测结果的准确性最好、预测的极端误差也最小;其次为 RBF;IDW预测的效果最差。 (2)OK_LU 在空间特征表达方面能够
更好地反映复杂地形区的局部变异,其插值结果的精度相比 OK有一定程度的提高,其平均绝对误差(MAE)从 0.900%降到了
0.567%,均方根误差(RMSE)从 1.101%降到了 0.777%,Pearson相关系数(R)从 0.4026提高到 0.5589,准确度(AC)从 0郾 9081提
高到 0.9505。 综合比较,在黄土丘陵地区,OK_LU能使插值结果的精度有较大提高,是土壤有机碳空间制图的有效途径。
关键词:土壤有机碳;空间插值;土地利用类型;黄土丘陵地区
Soil organic carbon interpolation based on auxiliary environmental covariates:
a case study at small watershed scale in Loess Hilly region
WEN Wen1,2, ZHOU Baotong1, WANG Yafeng2,*, HUANG Yong2,3
1 School of Geographic Sciences, Southwest University, Chongqing 400715,China
2 State Key Laboratory of Urban and Regional Ecology, Research Center for Eco鄄Environmental Sciences, Chinese Academy of Sciences, Beijing
100085, China
3 School of Forestry, Northeast Forestry University, Harbin 150040, China
Abstract: Maps of regional spatial distribution of soil organic carbon (SOC) in detail are needed to guide sustainable soil
uses and management decisions. Spatial interpolation methods include deterministic methods and geo鄄statistical methods.
Deterministic method is based on mathematical functions, whichsearch for the spatial similarities within a certain
geographical area, such as inverse distance weighting method( IDW), or create prediction surface by the known samples
based on smoothness such as radial鄄based function method(RBF); Geo鄄statistical method is one of the important means for
soil mapping supported by geographical and ecological parameters. It can achieve unbiased optimal estimation of regionalized
variables by using the statistical property of the known samples, and on the basis of theoretical analysis of semi鄄variogram,
such as ordinary kriging method(OK). Many of these methods perform well in areas with gentle terrains. However it is
uncertain how these methods perform to capture SOC variations in complex terrains, especially in those areas where land
uses are highly influenced by human activities, such as in the Loess Plateau of China. Therefore, four methods were applied
to predict SOC content spatial explicitly at typical small watershed in the Loess Plateau of China including Ordinary Kriging
method(OK), Inverse Distance Weighted ( IDW) method, Radial鄄Based Function method(RBF), and Ordinary kriging
http: / / www.ecologica.cn
modified with land鄄use type method(OK_LU) in the present research to evaluate their usability.
The purpose of this study is to find appropriate methods which are suitable to the complex terrain in Loess Plateau
region of China. The study area was a typical watershed in Loess Plateau with complex hilly gully terrain and various land鄄
use types. A field sampling dataset of 188 points was divided into two parts randomly:75% for model building and 25% for
accuracy validation. The accuracy validation points were used to compare the predicted value with measured value to check
the similarity at each point. Prediction results were validated using Pearson忆s Correlation Coefficient(R), Mean Absolute
Error (MAE), Root Mean Square Error ( RMSE) and Accuracy ( AC). The model with higher R and lower MAE and
RMSE, while the AC value is closer to 1 was regarded as more relevant, less biased and more accurate.
Results showed that: (1) The values of MAE and RMSE produced by OK method were the lowest, while the values of
R and AC were the largest, which indicated OK忆s highest effectiveness in SOC interpolation to yield the most accurate
results among the three ordinary methods. IDW method ranked behind OK and followed by RBF method in interpolation
effectiveness. (2) OK_LU was the best method for predicting SOC in this study. The prediction by OK_LU had more details
of depicting the local variation under land鄄use types. MAE by OK_LU decreased from 0.900% to 0.567%, RMSE decreased
from 1.101% to 0.777%, R increased from 0.4026 to 0.5589, and AC increased from 0.9081 to 0.9505 comparing with OK.
In conclusion, OK_LU gave more precise evaluation which was critical for obtaining accurate spatial distribution of SOC in
Loess Hilly regions. This will be of interests for related research in similar environments.
Key Words: soil organic carbon;spatial interpolation;land鄄use type;Loess Plateau
土壤碳库是陆地生态系统中最大的碳库,准确估算土壤碳库储量对正确评价土壤在陆地生态系统碳循环、全球碳循环以及
全球环境变化中的作用有重要意义[1鄄3] 。 土壤有机碳是土壤中较为活跃的部分,其含量和动态在土壤质量演变和全球碳循环
中起着十分重要的作用,也是估算土壤碳库储量的关键因素[4鄄5] 。 对土壤有机碳含量的空间差异进行研究,获得区域土壤有机
碳含量的连续空间分布,是精准农业技术体系中管理和分析土壤背景数据的一项基础性工作,然而大规模的土壤有机碳空间分
布信息的采集需要花费大量的人力、物力和财力。 因此,研究如何利用有限的样本数据来获得更为详尽的土壤有机碳空间分布
信息的方法具有重要意义。
空间插值模型是实现土壤有机碳含量从离散的点状信息向面状连续信息转换的有力工具,是表征土壤有机碳空间分布特
征的重要手段,可以将有限的采样点用来预测整个研究区域的土壤属性空间信息[6鄄9] 。 空间插值方法包括确定性方法和地理
统计方法[10] :确定性方法是使用数学函数进行插值,以研究区域内部的相似性(如反距离加权法),或以平滑度为基础(如径向
基函数法),由已知样点来创建预测表面的插值方法;地理统计方法是利用已知样点的统计属性,在半变异函数理论分析基础
上,对区域化变量的取值实现无偏最优估计的一种方法,如普通克里格法。 从学者们的研究中不难发现,普通克里格插值法是
使用频率最高的空间预测方法,但在地形复杂的区域,普通克里格方法的应用有一定的局限性[11] ,因而一些结合地貌类型、土
壤类型、土地利用类型和海拔高度等信息的克里格插值模型相继提出,通过辅助信息可以提高目标变量的预测精度[10鄄12] 。
土地利用和管理水平是影响土壤有机碳含量最普遍、最直接的因素之一,很大程度上影响着土壤有机碳变化的程度和方
向[13] ,国内外学者已对不同土地利用方式下土壤碳库特征进行了对比分析[14鄄16] 。 黄土丘陵地区,地形复杂,土地利用类型多
样,土地利用类型是土壤有机碳含量空间分布差异的重要影响因素之一,有必要在土壤有机碳的空间插值中考虑土地利用因
素。 对土壤属性空间预测的插值方法的对比研究很多,然而对于适宜黄土丘陵地形复杂地区插值方法的对比研究还较少,黄土
丘陵地区基于土地利用类型修正的土壤有机碳含量空间插值研究更不多见。 本研究针对黄土丘陵地区现有的土壤有机碳研究
中较少考虑的不同空间插值方法对预测精度的影响、土地利用类型修正对插值方法精度的影响等问题,以羊圈沟流域作为研究
区域,分析比较 3种常用的普通克里格法(OK)、反距离加权法(IDW)、径向基函数法(RBF),及经过土地利用类型修正的普通
克里格法(OK_LU)进行精度对比分析,揭示了最适宜地形复杂的黄土丘陵沟壑区土壤有机碳空间插值研究的方法,以期为类
似区域的土壤属性空间分布信息研究提供参考依据。
1摇 材料与方法
1.1摇 流域自然概况
羊圈沟流域位于延安市宝塔区东北方向 14 km处,属宝塔区李渠镇辖区,为延河左岸的二级支沟(36毅42忆N,109毅31忆E),该
流域属于碾庄沟流域的一级支流。 流域内存在一系列早期建造的淤地坝,自 1999年流域内实施了退耕还林政策。 该区域属黄
土丘陵沟壑区第二副区,年太阳总辐射量 5800 kJ / cm2,年日照时数 2528.8 h,一般年份能完全满足小麦、玉米、谷类作物的生长
要求。 流域植被在区划上属于森林草原过渡带。 流域内由于人为活动的干扰,自然植被破坏殆尽,多为人工种植而形成的次生
0936 摇 生摇 态摇 学摇 报摇 摇 摇 33卷摇
http: / / www.ecologica.cn
图 1摇 研究区域样点分布图
Fig.1摇 Map of sample distribution
植被,植物种类主要有刺槐 (Robiniapseudoacacia)、柳树 ( Salix
spp.)和杨树(Poplar spp.)等[17] 。
2011年羊圈沟流域的八种主要土地利用类型为林地、灌木
林地、草地、果园、梯田、坝地、水域和建设用地(表 1)。 主要的分
布格局为:林地、灌木林地、草地遍布在整个流域;果园分布在流
域中部;水域、居住用地和梯田主要分布在流域南部;坝地沿沟
分布(图 1)。
1.2摇 样品采集与分析
2011年分别在延安羊圈沟流域进行了两次野外全流域采
样。 本研究采用地形剖面线法进行坡面采样,采样时按土地利
用类型和土壤断面布设采样点,同时考虑微地貌,每个土壤断面
取 3至 6个点,取样器为荷兰 Eijkekamp公司产直径为 6 cm的半
圆凿型土钻。 采样深度为 0—20 cm,深度间隔为 0—5 cm,5—10
cm和 10—20 cm。 为体现采样点分布的均匀性,根据各地类所
占比例的大小,分地类进行样点的采集,以保证样点在不同土地
利用类型中的合理分布,如图 1、表 1 所示。 将野外采集的土壤
样品带回实验室后,先经过风干、研磨、过筛、称重和装盒等过程
后,再进入实验室分析阶段。 有机碳含量的测定使用 K2Cr2O4外
加热法。
1.3摇 数据处理
在 Mapgis软件中通过 GPS 将野外均匀布点采集样点的经
纬度坐标与羊圈沟流域的空间分布坐标相连接,将测得的数据
导入羊圈沟流域的点位表中,添加到土地利用类型图中,得到研
究区域采样点的分布图(图 1)。 将研究区域内采样点 0—5 cm,
5—10 cm,10—20 cm测得的土壤有机碳含量值综合,得到 0—20
cm土层厚度的土壤有机碳含量值。 本研究主要采用 0—20 cm
土层的土壤有机碳含量值作为研究对象。 通过对比分析 3 种常
用的插值方法与基于土地利用类型修正的方法的预测精度,分
析本研究区域的最优插值方法。
1.3.1摇 插值方法的选择
本研究中土壤有机碳的空间插值的 4 种方法均基于
ArcGIS10.0中的 ArcMap平台实现。
表 1摇 黄土丘陵区羊圈沟流域土地利用面积、采样数及其所占比例
Table 1摇 Land鄄use type area,sample size and their proportion in Yangjuangou watershed of hilly region,Loess Plateau
土地利用类型
Land use type
面积 / m2
Area
面积比例 / %
Area proportion
样本数
Sample number
样本比例 / %
Sample proportion
草地 Grass land 527745.31 25.69 43 22.87
灌木林地 Shrub land 517662.40 25.20 40 21.28
林地 Forest land 655663.03 31.92 66 35.11
梯田 Terrace land 218743.42 10.65 12 6.38
果园 Orchard land 10263.60 0.50 4 2.13
坝地 Dam land 88437.12 4.31 23 12.23
水域 Water 2677.85 0.13 0 0.00
居住和建筑用地 Construction land 32934.44 1.60 0 0.00
合计 Summation 2054127.16 100.00 188 100.00
(1)普通克里格法(Ordinary Kriging,简称 OK)
普通克里格法是利用区域化变量的原始数据和变异函数的结构特点,对未采样点的区域化变量的取值进行线性无偏最优
1936摇 19期 摇 摇 摇 文雯摇 等:基于辅助环境变量的土壤有机碳空间插值———以黄土丘陵区小流域为例 摇
http: / / www.ecologica.cn
估计的一种方法。 使用公式表示为:
Z(x0) = 移
n
i = 1
姿iz(xi) (1)
式中,Z为待估计的土壤有机碳含量栅格值,n为用于插值的采样点数目,Z(xi)为采样点土壤有机碳含量值,姿i为赋予土壤有机
碳含量数据的一组权重系数,并且其和等于 1,即:移
n
i = 1
姿i = 1,选取的 姿i使 Z(x0)的估计无偏,并且使方差小于任意观测值线性组
合的方差[18] 。
(2)逆距离权重法(Inverse Distance Weighting,简称 IDW)
逆距离权重法[19]是对采样点进行线性的加权来决定输出的栅格值,加权与距离成反比,输入点离输出栅格越远,它对输出
栅格的影响越小。 使用公式表示为:
Z =
移
n
i = 1
1
(di)p
z(xi)
移
n
i = 1
1
(di)p
(2)
式中,Z为待估计的土壤有机碳含量值栅格值, z(xi)为第 i( i= 1,2,…,n)个采样点的土壤有机碳含量值,n 为用于土壤有机碳
含量值插值的采样点个数,di为插值点到第 i个采样点的距离,p为距离的幂。 在本文中,p取为 2,即反距离平方插值。
(3)径向基函数插值
径向基函数插值[18]是假设有一个实值函数 F = F(x),要逼近其在某些点 x处的函数值。 用径向基函数来完成逼近的方
法是构造一个具有以下形式的逼近函数 S(x):
S(x) = 移
N
i = 1
ai渍( | | x - xi | | ),x 沂 Rd (3)
其中, 渍( t) ,t逸0,是一个确定的实值函数,即径向基函数, | |· | |表示欧几里得距离,ai,i = 1,…,N是待定系数。 点 xi,i= 1,
…,N为径向基函数逼近的节点。 相比于其他需要复杂处理的方法,径向基函数曲面的简单结构使得径向基函数逼近方法更加
具有吸引力。
(2)结合土地利用类型修正的克里格法(OK_LU)
以采样点分布图为底图,导出样点土壤有机碳含量与土地利用类型的关系数据表,并利用 SPSS18.0 软件对数据进行统计
分析,分别计算平均值、标准差、变异系数和相关系数,并利用单因素方差分析中的 LSD 进行多重比较。 明确土地利用类型对
土壤有机碳含量分布的影响后,在普通克里格插值方法的基础上,结合土地利用类型修正的克里格法即结合土地利用类型因子
进行修正,具体分两步进行:
第一步,先将所有预测点按照土地利用类型分为 6类,求出各土地利用类型的平均土壤有机碳含量值。 将每一个样本的土
壤有机碳含量值 Z(Xkj)分为相同类型均值 滋(Tk) 和离差 R(Xkj)之和。 用公式表示为:
Z X( )kj = 滋 T( )k + R X( )kj (4)
式中, Z(Xkj) 为样品的土壤有机碳含量值; 滋(Tk) 为相同土地利用类型样本的土壤有机碳含量均值; R(Xkj) 为样本土壤有机
碳含量与其相同土地利用类型样品的均值之差,称为“离差冶。 研究将离差作为一个新的区域变量 R(Xkj) 进行普通克里格插
值,空间插值利用 ArcGIS10.0软件中的地统计分析模块完成。
第二步,将插值过后的离差值图层 A按土地利用类型的不同,分别加上其地类 SOC均值,得到 SOC初步预测值图层 B。 各
待估点的土壤有机碳含量预测值为各土地利用类型均值与离差估计值之和。 由此得到研究区域土壤有机碳含量空间插值
结果。
1.3.2摇 检验方法
为检验不同插值方法在研究区域的预测精度,将 188 个实测样本在 Arcgis10.0 地统计工具下的 Subset 中随机构建两个子
集,75%的样本(N= 141)做预测子集,25%(N= 47)的样本做验证子集(图 1)。 本研究对验证样点的预测值与相同位置的实测
值进行相关性分析,对 SOC的空间插值结果进行精度检验。 用预测子集中的样点进行插值预测,得到验证子集中样点的 SOC
值,然后通过计算所有验证样点的预测值和观测值的 Pearson相关系数(R),平均绝对误差(MAE),均方根误差(RMSE),准确
度(AC)作为检验插值精度的标准,具体公式如下[20] :
R = 1
n - 1移
n
i = 1
Xi - X
-
S
æ
è
ç
ö
ø
÷
x
Yi - Y
-
S
æ
è
ç
ö
ø
÷
y
(5)
式中,
Xi-軈X
Sx
,軈X,Sx 分别为标准化变量,样本均值和样本标准差。
2936 摇 生摇 态摇 学摇 报摇 摇 摇 33卷摇
http: / / www.ecologica.cn
MAE = 1
n 移
n
i = 1
ABS(Voi - Vpi) (6)
RMSE = 1
n 移
n
i = 1
(Voi - Vpi) 2 (7)
式中,Voi是验证点土壤有机碳含量的的实测值,Vpi是验证点土壤有机碳含量的值,n是验证点数量。
AC = 1 - nRMSE
2
PE
(8)
其中 PE = 移
n
j = 1
Pi - O + Oi -[ ]O
2
(9)
式中,n为样本数;PE 为可能的误差变化( Potential error variance),Pi和 Oi分别为观测值和预测值;O 为实测值的平均值。 AC
取值范围为 0到 1,1代表预测值和实测值完全一致,而 0代表预测值和实测值完全不同。
2摇 研究结果
2.1摇 描述性统计分析
由表 2可知,羊圈沟流域各土地利用类型中土壤有机碳含量平均值的顺序为草地>灌木林地>林地>坝地>果园>梯田。 不
同的土地利用方式下,土壤有机碳含量值存在显著性差异,在探究土壤有机碳含量信息时,有必要考虑土地利用类型为辅助因
子,分地类进行插值。
表 2摇 不同土地利用方式下土壤有机碳含量状况
Table 2摇 Soil organic carbon content in different land鄄use type
土地利用类型
Land鄄use type
土壤有机碳含量 Soil organic carbon content / (g / kg)
平均值
Average value
标准差
Standard deviation
最小值
Minimum value
最大值
Maximum value
草地 4.90 A 1.27 2.46 9.19
灌木林地 4.89 AC 1.48 2.27 8.76
林地 4.87 AC 1.38 2.92 9.12
坝地 4.32 B 1.10 2.34 6.39
果园 4.16 C 1.01 3.31 5.56
梯田 3.34 C 1.16 2.48 4.48
总和 4.73 1.33 2.27 9.19
摇 摇 右上角大写字母表示平均值差异的显著性水平为 0.05(显著相关)
2.2摇 SOC空间结构分析
空间自相关理论是地理统计插值法进行预测的理论基础[21] ,所谓空间自相关是指在一定区域单元上,某种地理现象或某
一属性值与临近区域单元上同一现象或属性值的相关性,它用于衡量空间变量的变异对邻近区域相同变量的变异的依赖程
度[22鄄23] ,半变异函数是地统计学中研究土壤变异性的关键函数,是用来描述土壤性质的空间连续变异的一个连续函数,反映土
壤性质的不同距离观测值之间的变化。 存在空间依赖性的数据点对,距离越近,其值越相近;反之越远,其值越不同[24] 。 半变
异函数中,N表示其块金值,S表示基台值,本文应用 N / S和 S / (N+S)量化原始的 SOC 数据(应用于 OK)及 SOC 离差数据(应
用于 OK_LU)的空间结构强度。 N / S比值按空间结构强度的大小分为三类(N / S> 0.6,0.3示强烈的空间相关性;0.3 0.6,表示较弱空间相关性[25] 。 S / (N+S)比率越接近于 1,显
示其具有越强的空间相关性,反之,比率越接近于 0显示其空间相关性越弱[26] 。 OK和 OK_LU的半变异函数模型如图 2所示。
2.3摇 插值结果的验证与分析
图 3中别为运用 OK、IDW、RBF和 OK_LU插值法得到的羊圈沟流域 SOC空间分布图。 从整体上看,四图均反映出在流域
的西北部 SOC含量较高,流域南部含量低。 然而,各插值方法的结果图在局部上存在很大差异。 由图 3 可以看出,OK 的预测
值较连续,有机碳含量值分布很均匀,但是没有体现出不同土地利用类型的差异;图 3中,IDW和 RBF插值的结果图出现了“牛
眼冶现象,这是由于样点空间分布不均匀导致的[10] ,反距离权重法的插值表面不平滑,“牛眼冶现象较为严重;样条函数法的插
值牛眼现象较反距离权重法有所减轻,但表面仍不光滑;结合图 3和图 2的 2011年羊圈沟流域的土地利用现状图可以发现,OK
_LU插值表面比较连续、光滑,在林地、灌木林地和草地地区体现了较高的土壤有机碳含量值,而在梯田、坝地、裸地表现出的土
壤有机碳含量值较低,更详细具体地体现了不同土地利用方式对于土壤有机碳空间分布的影响,对于流域尺度的土壤有机碳的
空间变化,是一种比较理想的插值方法。
如表 3所示,SOC的 N / S比值(0.666)大于 0.6,OK_LU离差的 N / S比值(0.485)在 0.3—0.6之间,SOC离差的 S / (N+S)值
3936摇 19期 摇 摇 摇 文雯摇 等:基于辅助环境变量的土壤有机碳空间插值———以黄土丘陵区小流域为例 摇
http: / / www.ecologica.cn
图 2摇 黄土丘陵区羊圈沟流域土壤有机碳空间插值半变异函数图
Fig.2摇 Semi鄄variagram map of spatial interpolation of SOC in YangJuanGou watershed of hilly region,Loess Plateau
(0.8175)大于 SOC的 S / (N+S)值(0.7997),即 SOC和 SOC离差均具有较强的空间相关性,其中 SOC离差的空间相关性更强,
表明 OK和 OK_LU两种插值方法在此研究区域内较为适宜。
表 3摇 半变异函数模型
Table 3摇 Semi鄄variagram models
参数 Parameter 土壤有机碳含量 Soil organic carbon OK_LU离差 Residue of OK_LU
模型 Model 指数模型 Exponential 指数模型 Exponential
步长值 Range / km 100 100
块金值 N / (km / m2) 2 0.1551 0.3104
基台值 S / (km / m2) 2 0.6191 1.3901
N / S 0.2505 0.2233
S / (N+S) 0.7997 0.8175
从表 4可以看出,4种方法的实测值与预测值均具有一定的相关性,IDW 的 R = 0.3146(P = 0.0397)和 RBF 的 R = 0.3012
(P= 0.0385)是显著相关;OK的 R= 0.4026(P= 0.0050)和 OK_LU的 R= 0.5589(P= 0.0000)是极显著相关。 同时,通过均方根误
差(MAE)、平均绝对误差(RSME)和准确度(AC)来衡量 4种预测方法的系统误差及预测结果的准确性。 OK_LU平均绝对误差
和均方根误差(MAE= 0.7930,RSME= 0.2136)均小于 OK(MAE= 0.8665,RSME= 0.2927)、IDW(MAE= 0.8665,RSME= 0郾 2927)、
RBF(MAE= 0.8665,RSME= 0.2927);而 OK_LU(R= 0.5589,AC= 0.9505)相关性和准确度均高于 OK(R = 0.4026,AC= 0.9081)、
IDW(R= 0.3146,AC= 0.0397)、RBF(R= 0.3012,AC= 0.0385)。 可以看出,在对土壤有机碳进行空间插值时(特别是在黄土丘陵
沟壑地区) 结合土地利用类型,能使插值结果的精度得到很大的提高,是获得精确的土壤有机碳空间分布图是关键因素。
表 4摇 验证样点实测值与估算值相关性检验
Table 4摇 The correlation test of the measured value and the estimated value
普通克里格法
OK
反距离插值法
IDW
径向基函数法
RBF
基于土地利用类型修正的
普通克里格法
OK_LU
相关性 Correlation 0.4026 0.3146 0.3012 0.5589
显著性 Significance 0.0050** 0.0397* 0.0385* 0.0000**
平均绝对误差 Mean Absolute Error (MAE) 0.8665 0.8778 0.8707 0.7930
均方根误差 Root Mean Square Error(RMSE) 0.2927 0.3265 0.3109 0.2136
准确度 Accuracy(AC) 0.9081 0.8577 0.8647 0.9505
摇 摇 *表示相关显著性水平为 0.05;**表示相关显著性水平为 0.01
3摇 讨论
3.1摇 OK_LU的优越性
地统计法在地形变化不十分剧烈、比较平坦的平原地区可以得到较好的探测和模拟,其不仅考虑了待预测点与邻近样点数
据间的空间距离关系,还考虑了各参与预测的样点间的位置关系,充分利用了各样点数据的空间分布结构特征,地理统计方法
的这种优越性对于常常出现样本不充分、采点未经优化等问题的土壤属性数据的空间插值来说尤为重要[22] 。 本研究区的土壤
4936 摇 生摇 态摇 学摇 报摇 摇 摇 33卷摇
http: / / www.ecologica.cn
图 3摇 黄土丘陵区羊圈沟流域土壤有机碳空间分布图
Fig.3摇 Distribution of SOC in Yangjuangou watershed of hilly region,Loess Plateau
有机碳数据来自随机采样,OK使其插值结果比 IDW、RBF等传统方法更精确,更有效避免了系统误差的出现。
然而在黄土丘陵沟壑区,OK的应用却有一定的局限性[10] 。 由本文分析可知,在研究区内土地利用类型是影响 SOC 空间
分布差异的重要因素,因此结合土地利用类型对 SOC空间插值修正是非常必要的。 该区地形破碎、土地利用类型多样,基于半
变异函数的克里格插值难以将复杂的地形参数用于研究 SOC空间变化。 本文通过对普通克里格插值方法进行土地利用类型
修正,在黄土丘陵沟壑区对 SOC空间格局的模拟得到了较好效果。 范胜龙[11]等以福建省龙海市为例,设计了结合地貌类型、土
壤类型和土地利用类型等信息的克里格插值模型,结果表明设计的插值模型预测精度均高于普通克里格法,进一步验证了 OK_
LU的优越性。
3.2摇 土壤属性数据插值方法的选择
由于受到采样点密度、地表的复杂程度等影响,不同的插值方法生成的结果可能大相径庭[22] ,没有任何一种插值方法适合
所有地区,应依据对数据的定性分析、对研究区的地形分析,及对数据的空间探索分析,经过反复实验,选择最优的空间插值
5936摇 19期 摇 摇 摇 文雯摇 等:基于辅助环境变量的土壤有机碳空间插值———以黄土丘陵区小流域为例 摇
http: / / www.ecologica.cn
方法[26鄄27] 。
若研究目的仅在于对研究区域内几个粗略划分的亚区进行比较,如对黄土丘陵区最高值和最低值的区域进行大致定位,则
此时较小尺度的空间变异可以被忽略,只需把握住整体特征和全局趋势,以实施一系列插值优化措施为前提,本研究选用的 4
种方法均可符合要求。 若 SOC数据样点密度足够大、研究区域的地形变化不大时,可以直接采用 OK 进行插值,插值方法简便
快捷,预测结果准确性较高;但使用 OK插值相对其他地统计方法,步骤繁多,参数设置要求较高,耗时耗力,在研究目的许可的
情况下,简单、灵活的 IDW、RBF也是合适的选择。 如需利用插值结果进行更细致,更高级的空间分析,如图形叠加、构建生态
系统模型等,而数据具有某种变化趋势,则必须选择插值精度最高、在小空间尺度表现最好的插值方法,可以考虑将某种因素对
克里格插值进行修正,剔除趋势再进行预测[28] ,如本研究中的 OK_LU,因为误差会发生上行传递和放大,造成更深远的影响。
4摇 结论
在不同研究区域应根据不同研究目的选择不同的插值方法,不同插值方法各有利弊。 当常规方法不能满足预测精度时,可
通过辅助信息提高目标变量的预测精度。 在本研究区内,普通克里格法是 3种常规方法中插值效果最好的方法,径向基函数法
一般,反距离权重法插值效果最差。 黄土丘陵沟壑区生态环境脆弱,地形较复杂,土地利用类型多样,不同土地利用类型的土壤
有机碳含量平均值差异明显,基于土地利用类型修正的普通克里格插值法对 SOC 含量空间分布的预测精度有一定程度的提
高,是地形复杂的黄土丘陵沟壑区土壤有机碳空间插值较适宜的方法。
References:
[ 1 ]摇 Post W M, Peng T H, Emanuel W R, King A W, Dale V H, De Angelis D L. The global carbon cycle. American Scientist, 1990, 78: 310鄄326.
[ 2 ] 摇 Batjes N H. Total carbon and nitrogen in the soils of the world. European Journal of Soil Science, 1996, 47(2): 151鄄163.
[ 3 ] 摇 Lal R. Soil carbon sequestration impacts on global climate change and food security. Science, 2004, 304(5677): 1623鄄1627.
[ 4 ] 摇 Post W M, Kwon K C. Soil carbon sequestration and land鄄use change: processes and potential. Global Change Biology, 2000, 6(3): 317鄄327.
[ 5 ] 摇 Lal R, Kimble J, Levine E, Whitman C. World soils and greenhouse effect: an overview / / Lal R, Kimble J, Levine E, Stewart B A, eds. Soils
and Global Changes. Boca Raton: CRC Lewis Publishers, 1995: 1鄄25.
[ 6 ] 摇 Hu K L, Li B G, Lin Q M, Li G T, Chen D L. Spatial variability of soil nutrient in wheat field. Transactions of the Chinese Society of Agricultural
Engineering, 1999, 15(3): 33鄄38.
[ 7 ] 摇 Wang J, Qiu Y. Review on spatial variability and scale effects of land quality. Progress in Geography, 2005, 24(4): 28鄄35.
[ 8 ] 摇 Zhang S J, Fang H, He Y. Sampling strategies of field information on precision agriculture. Transactions of the Chinese Society for Agricultural
Machinery, 2004, 35(4): 88鄄92.
[ 9 ] 摇 Sauer T J, Meek D W. The spatial variation of plant鄄available phosphorus in pastures with contrasting management. Soil Science Society of America
Journal, 2003, 67(3): 826鄄836.
[10] 摇 Yao X L, Fu B J, L俟 Y H, Wang S, Liu M. Comparison of four spatial interpolation methods for estimating soil moisture in a complex terrain
catchment. PLoS One, 2013, 1(8): 1鄄13.
[11] 摇 Zhang T C, Chang Q R, Liu J. Comparison of spatial interpolation methods for soil nutrient elements A case study of Yuyang County, Shaanxi
Province. Agricultural Research in the Arid Areas, 2010, 28(2): 177鄄182.
[12] 摇 Liu X N, Huang F, Wang P. Principles and Methods of GIS Spatial Analysis. Beijing: Science Press, 2008: 185鄄223.
[13] 摇 Liu J, Cai G X, Liu H B. Research of the interpolation methods for determining soil organic matter content in the hilly regions of southwestern
China. Journal of Southwest University: Natural Science Edition, 2008, 30(3): 107鄄112.
[14] 摇 Fan S L, Huang Y H, Lin J S. The optimized interpolation models and its relationship with soil sampling density on detecting spatial variability of
farmland soil organic carbon A case study in Longhai city, Fujian province. Research of Soil and Water Conservation, 2011, 18(6): 1鄄5.
[15] 摇 Allmaras R R, Schomberg H H, Douglas C L, Dao T H Jr. Soil organic carbon sequestration potential of adopting conservation tillage in U. S.
croplands. Journal of Soil and Water Conservation, 2000, 55(3): 365鄄373.
[16] 摇 Chen J Y, Sun B, Li Z P, Li H Y, Pan J J. Pool size of soil organic carbon and dynamics under different land use. Journal of Soil and Water
Conservation, 2008, 22(1): 91鄄95, 179鄄179.
[17] 摇 Wang Y F, Fu B J, Chen L D, L俟 Y H, Luo C Y. Effects of land use change on soil erosion intensity in small watershed of Loess Hilly Region: A
quantitative evaluation with 137鄄Cesium tracer. Chinese Journal of Applied Ecology, 2009, 20(7): 1571鄄1576.
[18] 摇 Fang X, Xu G L, Hong Y, Jing W F. Study on soil organic carbon density under different land use types in central hilly area of Hu忆nan province.
Journal of Soil and Water Conservation, 2009, 23(1): 54鄄58, 63鄄63.
[19] 摇 Zhu H Y, Liu S L, Jia S F. Problems of the spatial interpolation of physical geographical elements. Geographical Research, 2004, 23( 4):
425鄄432.
[20] 摇 Yang L, Zhu A X, Qin C Z, Li B L, Pei T, Liu B Y. Soil property mapping using fuzzy membership A case study of a study area in Heshan
farm of Heilongjiang province. Acta Pedologica Sinica, 2009, 46(1): 9鄄15.
[21] 摇 Lin L, Li C H, Dai M, Cai W G, Lin Q. Optimization of the spatial interpolation for marine phytoplankton abundance. Acta Ecologica Sinica,
6936 摇 生摇 态摇 学摇 报摇 摇 摇 33卷摇
http: / / www.ecologica.cn
2007, 27(7): 2880鄄2888.
[22] 摇 Gringarten E, Deutsch C V. Teacher忆s aide variogram interpretation and modeling. Mathematical Geology, 2001, 33(4): 507鄄534.
[23] 摇 Li X, Cheng G D, Lu L. Comparison study of spatial interpolation methods of air temperature over Qinghai鄄Xizang plateau. Plateau Meteorology,
2003, 22(6): 565鄄573.
[24] 摇 Leenaers H, Okx J P, Burrough P A. Comparison of spatial prediction methods for mapping floodplain soil pollution. Catena, 1990, 17(6):
535鄄550.
[25] 摇 Cambardella C A, Karlen D L. Spatial analysis of soil fertility parameters. Precision Agriculture, 1999, 1(1): 5鄄14.
[26] 摇 Lam N S N. Spatial interpolation methods: A review. American Cartographer, 1983, 10(2): 129鄄149.
[27] 摇 Liu J S, Chen H, Yang B Y, Wang W, Xiang Y, Zhao C. Comparison of interpolation methods on annual mean precipitation in Hebei Province.
Acta Ecologica Sinica, 2009, 29(7): 3493鄄3500.
[28] 摇 Lado L R, Polya D, Winkel L, Berg M, Hegan A. Modelling arsenic hazard in Cambodia: A geostatistical approach using ancillary data. Applied
Geochemistry, 2008, 23(11): 3010鄄3018.
参考文献:
[ 6 ]摇 胡克林, 李保国, 林启美, 李贵桐, 陈德立. 农田土壤养分的空间变异性特征. 农业工程学报, 1999, 15(3): 33鄄38.
[ 7 ] 摇 王军, 邱扬. 土地质量的空间变异与尺度效应研究进展. 地理科学进展, 2005, 24(4): 28鄄35.
[ 8 ] 摇 张淑娟, 方慧, 何勇. 精细农业田间信息采样策略. 农业机械学报, 2004, 35(4): 88鄄92.
[11] 摇 张铁婵, 常庆瑞, 刘京. 土壤养分元素空间分布不同插值方法研究———以榆林市榆阳区为例. 干旱地区农业研究, 2010, 28( 2):
177鄄182.
[12] 摇 刘湘南, 黄方, 王平. GIS 空间分析原理与方法. 北京: 科学出版社, 2008: 185鄄223.
[13] 摇 刘静, 蔡国学, 刘洪斌. 西南丘陵地区土壤有机质含量的空间插值法研究. 西南大学学报: 自然科学版, 2008, 30(3): 107鄄112.
[14] 摇 范胜龙, 黄炎和, 林金石. 表征土壤有机碳区域分布的优化空间插值模型研究———以福建省龙海市为例. 水土保持研究, 2011, 18(6):
1鄄5.
[16] 摇 陈锦盈, 孙波, 李忠佩, 李海鹰, 潘剑君. 不同土地利用类型土壤有机碳各库大小及分解动态. 水土保持学报, 2008, 22(1): 91鄄 95,
179鄄179.
[17] 摇 汪亚峰, 傅伯杰, 陈利顶, 吕一河, 罗春燕. 黄土丘陵小流域土地利用变化的土壤侵蚀效应: 基于 137Cs示踪的定量评价. 应用生态学
报, 2009, 20(7): 1571鄄1576.
[18] 摇 方晰, 徐桂林, 洪瑜, 金文芬. 湘中丘陵区不同土地利用方式下土壤有机碳密度. 水土保持学报, 2009, 23(1): 54鄄58, 63鄄63.
[19] 摇 朱会义, 刘述林, 贾绍凤. 自然地理要素空间插值的几个问题. 地理研究, 2004, 23(4): 425鄄432.
[20] 摇 杨琳, 朱阿兴, 秦承志, 李宝林, 裴韬, 刘宝元. 运用模糊隶属度进行土壤属性制图的研究———以黑龙江鹤山农场研究区为例. 土壤学
报, 2009, 46(1): 9鄄15.
[21] 摇 林琳, 李纯厚, 戴明, 蔡文贵, 林钦. 海洋浮游植物丰度的空间插值优化. 生态学报, 2007, 27(7): 2880鄄2888.
[23] 摇 李新, 程国栋, 卢玲. 青藏高原气温分布的空间插值方法比较. 高原气象, 2003, 22(6): 565鄄573.
[27] 摇 刘劲松, 陈辉, 杨彬云, 王卫, 相云, 赵超. 河北省年均降水量插值方法比较. 生态学报, 2009, 29(7): 3493鄄3500.
7936摇 19期 摇 摇 摇 文雯摇 等:基于辅助环境变量的土壤有机碳空间插值———以黄土丘陵区小流域为例 摇
粤悦栽粤 耘悦韵蕴韵郧陨悦粤 杂陨晕陨悦粤 灾燥造援猿猿袁晕燥援员怨 韵糟贼援袁圆园员猿渊杂藻皂蚤皂燥灶贼澡造赠冤
悦韵晕栽耘晕栽杂
粤 则藻增蚤藻憎 燥枣 藻糟燥泽赠泽贼藻皂 泽藻则增蚤糟藻泽 葬灶凿 则藻泽藻葬则糟澡 责藻则泽责藻糟贼蚤增藻泽 酝粤 云藻灶早躁蚤葬燥袁蕴陨哉 允蚤灶贼燥灶早袁 粤援 耘早则蚤灶赠葬 耘灶藻躁蚤 渊缘怨远猿冤噎噎噎噎噎噎噎噎噎杂藻曾怎葬造 蚤灶贼藻则枣藻则藻灶糟藻 蚤灶 灶燥灶鄄澡怎皂葬灶 责则蚤皂葬贼藻泽 再粤晕郧 月蚤灶袁 宰粤晕郧 悦澡藻灶早造蚤葬灶早袁 允陨 宰藻蚤澡燥灶早袁 藻贼 葬造 渊缘怨苑猿冤噎噎噎噎噎噎噎噎噎噎噎噎阅藻灶泽蚤贼赠鄄凿藻责藻灶凿藻灶贼 藻枣枣藻糟贼 燥灶 则藻责则燥凿怎糟贼蚤燥灶 燥枣 则燥凿藻灶贼泽院 葬 则藻增蚤藻憎 匀粤晕 匝怎灶澡怎葬袁 郧哉韵 悦燥灶早袁在匀粤晕郧 酝藻蚤憎藻灶 渊缘怨愿员冤噎噎噎噎噎噎噎孕则燥曾蚤皂葬贼藻 葬灶凿 怎造贼蚤皂葬贼藻 凿藻贼藻则皂蚤灶葬灶贼泽 燥枣 枣燥燥凿 糟澡葬蚤灶 造藻灶早贼澡 宰粤晕郧 再怎赠怎袁 载哉 允怎灶袁 蕴耘陨 郧怎葬灶早糟澡怎灶 渊缘怨怨园冤噎噎噎噎噎噎噎噎噎噎噎酝藻糟澡葬灶蚤泽皂 燥枣 遭蚤燥造燥早蚤糟葬造 糟燥灶贼则燥造 贼燥 责造葬灶贼 凿蚤泽藻葬泽藻泽 怎泽蚤灶早 葬则遭怎泽糟怎造葬则 皂赠糟燥则则澡蚤扎葬造 枣怎灶早蚤蕴哉韵 匝蚤葬燥赠怎袁 宰粤晕郧 载蚤葬燥躁怎葬灶袁 蕴陨 再怎葬灶赠怎葬灶袁 藻贼 葬造 渊缘怨怨苑冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎粤凿增葬灶糟藻泽 蚤灶 藻枣枣藻糟贼泽 燥枣 糟燥灶泽藻则增葬贼蚤燥灶 贼蚤造造葬早藻 燥灶 泽燥蚤造 燥则早葬灶蚤糟 糟葬则遭燥灶 葬灶凿 灶蚤贼则燥早藻灶载哉耘 允蚤葬灶枣怎袁 在匀粤韵 载蚤灶袁 杂澡葬凿则葬糟噪 月葬贼泽蚤造藻 阅蚤噪早憎葬贼造澡藻袁 藻贼 葬造 渊远园园远冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎匀葬遭蚤贼葬贼 泽藻造藻糟贼蚤燥灶 燥枣 贼澡藻 责则藻鄄则藻造藻葬泽藻凿 早蚤葬灶贼 责葬灶凿葬 蚤灶 宰燥造燥灶早 晕葬贼怎则藻 砸藻泽藻则增藻在匀粤晕郧 酝蚤灶早糟澡怎灶袁 匀哉粤晕郧 再葬灶袁 蕴陨 阅藻泽澡藻灶早袁 藻贼 葬造 渊远园员源冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎粤糟贼蚤增蚤贼赠 则澡赠贼澡皂 葬灶凿 遭藻澡葬增蚤燥则葬造 贼蚤皂藻 遭怎凿早藻贼泽 燥枣 憎蚤造凿 砸藻藻增藻泽忆泽 孕澡藻葬泽葬灶贼 渊杂赠则皂葬贼蚤糟怎泽 则藻藻增藻泽蚤蚤冤 怎泽蚤灶早 蚤灶枣则葬则藻凿 糟葬皂藻则葬在匀粤韵 再怎扎藻袁 宰粤晕郧 在澡蚤糟澡藻灶袁 载哉 允蚤造蚤葬灶早袁 藻贼 葬造 渊远园圆员冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎栽澡藻 藻灶藻则早赠 遭怎凿早藻贼 燥枣 贼则藻藻 泽责葬则则燥憎泽 孕葬泽泽藻则 皂燥灶贼葬灶怎泽 蚤灶 憎蚤灶凿 凿蚤枣枣藻则藻灶贼 泽责藻藻凿 葬灶凿 凿怎则葬贼蚤燥灶再粤晕郧 在澡蚤澡燥灶早袁 宰哉 匝蚤灶早皂蚤灶早袁 阅韵晕郧 匀葬蚤赠葬灶袁藻贼 葬造 渊远园圆愿冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎晕藻泽贼 泽蚤贼藻 糟澡葬则葬糟贼藻则蚤泽贼蚤糟泽 燥枣 孕藻贼葬怎则蚤泽贼葬 糟葬灶蚤糟藻责泽 蚤灶 月葬蚤皂葬 杂灶燥憎 酝燥怎灶贼葬蚤灶 晕葬贼怎则藻 砸藻泽藻则增藻蕴陨 再葬灶澡燥灶早袁 郧哉粤晕 允蚤灶噪藻袁 蕴陨 阅葬赠燥灶早袁 匀哉 允蚤藻 渊远园猿缘冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎耘枣枣藻糟贼泽 燥枣 澡葬遭蚤贼葬贼 枣则葬早皂藻灶贼葬贼蚤燥灶 燥灶 贼澡藻 早藻灶藻贼蚤糟 凿蚤增藻则泽蚤贼赠 燥枣 孕葬糟澡赠糟燥灶凿赠造葬 造怎贼藻蚤责藻泽 燥灶 蚤泽造葬灶凿泽 蚤灶 贼澡藻 栽澡燥怎泽葬灶凿 陨泽造葬灶凿 蕴葬噪藻袁 耘葬泽贼悦澡蚤灶葬 蕴哉韵 再怎葬灶赠怎葬灶袁 蕴陨哉 允蚤灶造蚤葬灶早袁 匀哉粤晕郧 允蚤藻造蚤灶早袁 藻贼 葬造 渊远园源员冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎栽澡藻 皂燥造藻糟怎造葬则 早藻灶藻贼蚤糟 则藻造葬贼蚤燥灶泽澡蚤责 遭藻贼憎藻藻灶 贼澡藻 责燥造造蚤灶葬贼燥则泽 燥枣 云蚤糟怎泽 责怎皂蚤造葬 增葬则援 责怎皂蚤造葬 葬灶凿 云蚤糟怎泽 责怎皂蚤造葬 增葬则援 葬憎噪藻燥贼泽葬灶早宰哉 宰藻灶泽澡葬灶袁 悦匀耘晕 再燥怎造蚤灶早袁 杂哉晕 蕴蚤灶早造蚤袁 藻贼 葬造 渊远园源怨冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎栽澡藻 早藻灶藻贼蚤糟 藻增燥造怎贼蚤燥灶葬则赠 则藻造葬贼蚤燥灶泽澡蚤责泽 燥枣 贼憎燥 耘怎责则蚤泽贼蚤灶葬 泽责藻糟蚤藻泽 燥灶 云蚤糟怎泽 葬造贼蚤泽泽蚤皂葬悦匀耘晕 再燥怎造蚤灶早袁 杂哉晕 蕴蚤灶早造蚤袁 宰哉 蕴藻蚤造藻蚤袁 藻贼 葬造 渊远园缘愿冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎酝藻贼葬造 怎责贼葬噪藻 葬灶凿 则燥燥贼 皂燥则责澡燥造燥早蚤糟葬造 糟澡葬灶早藻泽 枣燥则 贼憎燥 增葬则蚤藻贼蚤藻泽 燥枣 杂葬造蚤曾 蚤灶贼藻早则葬 怎灶凿藻则 糟葬凿皂蚤怎皂 泽贼则藻泽泽宰粤晕郧 杂澡怎枣藻灶早袁 杂匀陨 载蚤葬灶早袁 杂哉晕 匀葬蚤躁蚤灶早袁 藻贼 葬造 渊远园远缘冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎耘枣枣藻糟贼泽 燥枣 责澡贼澡葬造蚤糟 葬糟蚤凿 燥灶 泽藻藻凿 早藻则皂蚤灶葬贼蚤燥灶袁 皂藻皂遭则葬灶藻 造蚤责蚤凿 责藻则燥曾蚤凿葬贼蚤燥灶 葬灶凿 燥泽皂燥则藻早怎造葬贼蚤燥灶 泽怎遭泽贼葬灶糟藻 燥枣 则葬凿蚤泽澡 泽藻藻凿造蚤灶早泽再粤晕郧 再葬灶躁蚤藻袁 宰粤晕郧 载蚤葬燥憎藻蚤袁 在匀粤韵 运葬灶早袁 藻贼 葬造 渊远园苑源冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎栽澡藻 皂燥则责澡燥造燥早蚤糟葬造 葬灶凿 责澡赠泽蚤燥造燥早蚤糟葬造 则藻泽责燥灶泽藻泽 燥枣 栽葬皂葬则蚤曾 则葬皂燥泽蚤泽泽蚤皂葬 泽藻藻凿造蚤灶早 贼燥 凿蚤枣枣藻则藻灶贼 蚤则则蚤早葬贼蚤燥灶 皂藻贼澡燥凿泽 蚤灶 贼澡藻 藻曾贼则藻皂藻造赠葬则蚤凿 葬则藻葬 酝粤 载蚤葬燥凿燥灶早袁 宰粤晕郧 酝蚤灶早澡怎蚤袁 蕴陨 宰藻蚤澡燥灶早袁 藻贼 葬造 渊远园愿员冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎砸藻泽责燥灶泽藻 糟澡葬则葬糟贼藻则蚤泽贼蚤糟泽 燥枣 责澡燥贼燥泽赠灶贼澡藻贼蚤糟 葬灶凿 责澡赠泽蚤燥造燥早蚤糟葬造 责葬则葬皂藻贼藻则泽 蚤灶 在蚤扎蚤责澡怎泽 躁怎躁怎遭葬 增葬则援 泽责蚤灶燥泽怎泽 泽藻藻凿造蚤灶早 造藻葬增藻泽 贼燥 泽燥蚤造憎葬贼藻则 蚤灶 泽葬灶凿 澡葬遭蚤贼葬贼 枣燥则皂藻凿 枣则燥皂 泽藻葬泽澡藻造造泽 宰粤晕郧 砸燥灶早则燥灶早袁 载陨粤 允蚤葬灶早遭葬燥袁 再粤晕郧 允蚤澡怎葬袁 藻贼 葬造 渊远园愿愿冤噎噎噎噎噎噎噎噎噎耘枣枣藻糟贼泽 燥枣 糟藻则葬皂泽蚤贼藻 皂怎造糟澡蚤灶早 燥灶 泽燥蚤造 憎葬贼藻则 糟燥灶贼藻灶贼袁 责澡燥贼燥泽赠灶贼澡藻贼蚤糟 责澡赠泽蚤燥造燥早蚤糟葬造 糟澡葬则葬糟贼藻则蚤泽贼蚤糟泽 葬灶凿 早则燥憎贼澡 燥枣 责造葬灶贼泽栽粤晕 载怎藻澡燥灶早袁 郧哉韵 载蚤葬燥责蚤灶早袁 在匀粤韵 栽蚤灶早灶蚤灶早 渊远园怨苑冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎阅赠灶葬皂蚤糟泽 燥枣 贼葬灶灶蚤灶 糟燥灶糟藻灶贼则葬贼蚤燥灶 葬灶凿 灶怎贼则蚤藻灶贼 则藻泽燥则责贼蚤燥灶 枣燥则 遭则葬灶糟澡造藻贼泽 燥枣 悦葬泽怎葬则蚤灶葬 藻择怎蚤泽藻贼蚤枣燥造蚤葬 责造葬灶贼葬贼蚤燥灶泽 葬贼 凿蚤枣枣藻则藻灶贼 葬早藻泽再耘 郧燥灶早枣怎袁 在匀粤晕郧 杂澡葬灶早躁怎袁 在匀粤晕郧 蕴蚤澡怎葬袁 藻贼 葬造 渊远员园苑冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎杂怎造枣怎则 糟燥灶贼藻灶贼泽 蚤灶 造藻葬增藻泽 葬灶凿 遭则葬灶糟澡藻泽 燥枣 凿燥皂蚤灶葬灶贼 泽责藻糟蚤藻泽 葬皂燥灶早 贼澡藻 贼澡则藻藻 枣燥则藻泽贼 贼赠责藻泽 蚤灶 贼澡藻 孕藻葬则造 砸蚤增藻则 阅藻造贼葬孕耘陨 晕葬灶糟葬蚤袁 悦匀耘晕 月怎枣藻灶早袁 在韵哉 在澡蚤躁蚤灶袁 藻贼 葬造 渊远员员源冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎陨皂责葬糟贼泽 燥枣 葬则遭怎泽糟怎造葬则 皂赠糟燥则则澡蚤扎葬造 枣怎灶早蚤 葬灶凿 责澡燥泽责澡燥则怎泽 燥灶 早则燥憎贼澡 凿赠灶葬皂蚤糟泽 燥枣 月葬怎澡蚤灶蚤葬 枣葬遭藻则蚤 泽藻藻凿造蚤灶早泽杂韵晕郧 悦澡藻灶早躁怎灶袁 匝哉 蕴葬蚤赠藻袁 酝粤 运藻皂蚤灶早袁 藻贼 葬造 渊远员圆员冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎悦澡葬则葬糟贼藻则蚤泽贼蚤糟泽 燥枣 蚤燥灶 葬糟糟怎皂怎造葬贼蚤燥灶 葬灶凿 泽藻藻凿 早藻则皂蚤灶葬贼蚤燥灶 枣燥则 泽藻藻凿泽 枣则燥皂 责造葬灶贼泽 糟怎造贼怎则藻凿 葬贼 凿蚤枣枣藻则藻灶贼 糟燥灶糟藻灶贼则葬贼蚤燥灶泽 燥枣 灶蚤贼则葬贼藻灶蚤贼则燥早藻灶 葬灶凿 泽葬造蚤灶蚤贼赠 在匀韵哉 允蚤葬糟澡葬燥袁 云哉 栽蚤灶早贼蚤灶早袁 在匀粤韵 宰藻蚤憎藻蚤袁 藻贼 葬造 渊远员圆怨冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎孕澡赠泽蚤燥鄄藻糟燥造燥早蚤糟葬造 藻枣枣藻糟贼泽 燥枣 藻灶凿燥责澡赠贼藻 蚤灶枣藻糟贼蚤燥灶 燥灶 贼澡藻 澡燥泽贼 早则葬泽泽 憎蚤贼澡 藻造藻增葬贼藻凿 悦韵圆 杂匀陨 在澡蚤遭蚤灶早袁 在匀韵哉 再燥灶早袁 蕴陨 载蚤葬袁 藻贼 葬造 渊远员猿缘冤
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎耘枣枣藻糟贼泽 燥枣 责则藻贼则藻葬贼皂藻灶贼 燥灶 早藻则皂蚤灶葬贼蚤燥灶 燥枣 栽赠责澡葬 凿燥皂蚤灶早藻灶泽蚤泽 葬灶凿 孕澡则葬早皂蚤贼藻泽 葬怎泽贼则葬造蚤泽酝耘晕郧 匀怎葬灶袁 宰粤晕郧 载怎藻澡燥灶早袁 栽韵晕郧 杂澡燥怎扎澡藻灶早袁 藻贼 葬造 渊远员源圆冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎栽则葬灶泽枣藻则 糟澡葬则葬糟贼藻则蚤泽贼蚤糟泽 燥枣 糟葬凿皂蚤怎皂 枣则燥皂 泽燥蚤造 贼燥 杂葬造蚤曾 伊 葬怎则藻燥鄄责藻灶凿怎造葬 在匀粤晕郧 宰藻灶袁 宰耘陨 匀燥灶早袁 杂哉晕 载蚤葬燥糟葬灶袁 藻贼 葬造 渊远员源苑冤噎噎耘枣枣藻糟贼 燥枣 悦造燥泽藻鄄贼燥鄄晕葬贼怎则藻 皂葬灶葬早藻皂藻灶贼 燥灶 贼澡藻 灶葬贼怎则葬造 则藻早藻灶藻则葬贼蚤燥灶 葬灶凿 泽责藻糟蚤藻泽 凿蚤增藻则泽蚤贼赠 蚤灶 葬 皂葬泽泽燥灶 责蚤灶藻 责造葬灶贼葬贼蚤燥灶蕴哉韵 再蚤灶早澡怎葬袁 杂哉晕 阅燥灶早躁蚤灶早袁 蕴陨晕 允蚤葬灶赠燥灶早袁 藻贼 葬造 渊远员缘源冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎孕燥责怎造葬贼蚤燥灶 凿赠灶葬皂蚤糟泽 葬灶凿 泽藻藻凿 遭葬灶噪泽 燥枣 贼澡藻 贼澡则藻葬贼藻灶藻凿 泽藻葬早则葬泽泽 匀葬造燥责澡蚤造葬 遭藻糟糟葬则蚤蚤 蚤灶 孕藻葬则造 月葬赠袁 郧怎葬灶早曾蚤匝陨哉 郧怎葬灶早造燥灶早袁 云粤晕 匀葬灶早择蚤灶早袁 蕴陨 在燥灶早泽澡葬灶袁 藻贼 葬造 渊远员远猿冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎耘枣枣藻糟贼泽 燥枣 遭蚤燥造燥早蚤糟葬造 糟则怎泽贼泽 燥灶 凿藻憎 凿藻责燥泽蚤贼蚤燥灶 葬灶凿 藻增葬责燥则葬贼蚤燥灶 蚤灶 贼澡藻 杂燥怎贼澡藻则灶 耘凿早藻 燥枣 贼澡藻 酝怎 哉泽 杂葬灶凿赠 蕴葬灶凿袁 晕燥则贼澡藻则灶 悦澡蚤灶葬再陨晕 砸怎蚤责蚤灶早袁宰哉 再燥灶早泽澡藻灶早袁 在匀粤晕郧 载蚤灶袁 藻贼 葬造 渊远员苑猿冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎蕴蚤枣藻 澡蚤泽贼燥则赠 糟澡葬则葬糟贼藻则蚤泽贼蚤糟泽 葬灶凿 泽责葬贼蚤葬造 凿蚤泽贼则蚤遭怎贼蚤燥灶 燥枣 孕燥责怎造怎泽 责则怎蚤灶燥泽葬 责燥责怎造葬贼蚤燥灶 葬贼 贼澡藻 怎责责藻则 则藻葬糟澡藻泽 燥枣 栽葬则蚤皂 砸蚤增藻则匀粤晕 蕴怎袁 载陨 蕴蚤灶择蚤葬燥袁 宰粤晕郧 允蚤葬择蚤葬灶早袁 藻贼 葬造 渊远员愿员冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎陨灶贼藻则葬糟贼蚤增藻 藻枣枣藻糟贼泽 燥枣 泽澡燥则贼鄄贼藻则皂 灶蚤贼则燥早藻灶 藻灶则蚤糟澡皂藻灶贼 葬灶凿 泽蚤皂怎造葬贼藻凿 早则葬扎蚤灶早 燥灶 藻糟燥泽赠泽贼藻皂 则藻泽责蚤则葬贼蚤燥灶 蚤灶 葬灶 葬造责蚤灶藻 皂藻葬凿燥憎 燥灶 贼澡藻栽蚤遭藻贼葬灶 孕造葬贼藻葬怎 在韵晕郧 晕蚤灶早袁 杂匀陨 孕藻蚤造蚤袁 允陨粤晕郧 允蚤灶早袁 藻贼 葬造 渊远员怨员冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
栽澡藻 糟燥则则藻造葬贼蚤燥灶 遭藻贼憎藻藻灶 泽燥蚤造 憎葬贼藻则 泽葬造蚤灶蚤贼赠 葬灶凿 责造葬灶贼 糟燥皂皂怎灶蚤贼赠 凿蚤泽贼则蚤遭怎贼蚤燥灶 怎灶凿藻则 皂蚤糟则燥鄄贼燥责燥早则葬责澡赠 蚤灶 杂燥灶早灶藻灶 孕造葬蚤灶再粤晕郧 云葬灶袁 宰粤晕郧 在澡蚤糟澡怎灶袁 宰粤晕郧 再怎灶澡藻袁 藻贼 葬造 渊远圆园圆冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎悦燥皂责葬则蚤泽燥灶 燥枣 栽杂孕袁 孕酝圆援缘葬灶凿 贼澡藻蚤则 憎葬贼藻则鄄泽燥造怎遭造藻 蚤燥灶泽 枣则燥皂 遭燥贼澡 蚤灶泽蚤凿藻 葬灶凿 燥怎贼泽蚤凿藻 燥枣 阅葬枣怎泽澡葬灶 枣燥则藻泽贼 责葬则噪 蚤灶 郧怎葬灶早扎澡燥怎凿怎则蚤灶早 则葬蚤灶赠 泽藻葬泽燥灶 载陨粤韵 再蚤澡怎葬袁蕴陨 允蚤燥灶早袁运哉粤晕郧 再怎葬灶憎藻灶袁藻贼 葬造 渊远圆园怨冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎云蚤泽澡 糟燥皂皂怎灶蚤贼赠 藻糟燥造燥早赠 蚤灶 则燥糟噪赠 则藻藻枣 澡葬遭蚤贼葬贼 燥枣 酝葬忆葬灶 粤则糟澡蚤责藻造葬早燥 域援 杂责葬贼蚤燥鄄贼藻皂责燥则葬造 责葬贼贼藻则灶泽 燥枣 糟燥皂皂怎灶蚤贼赠 泽贼则怎糟贼怎则藻宰粤晕郧 在澡藻灶澡怎葬袁 在匀粤韵 允蚤灶早袁 宰粤晕郧 运葬蚤袁 藻贼 葬造 渊远圆员愿冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎陨灶贼藻则葬灶灶怎葬造 增葬则蚤葬贼蚤燥灶 蚤灶 贼澡藻 责燥责怎造葬贼蚤燥灶 凿赠灶葬皂蚤糟泽 燥枣 泽灶葬蚤造枣蚤泽澡 蕴蚤责葬则蚤泽 贼葬灶葬噪葬藻 蚤灶 贼澡藻 再藻造造燥憎 杂藻葬悦匀耘晕 再怎灶造燥灶早袁 杂匀粤晕 载蚤怎躁怎葬灶袁 在匀韵哉 在澡蚤责藻灶早袁 藻贼 葬造 渊远圆圆苑冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎杂责葬贼蚤葬造 葬灶凿 贼藻皂责燥则葬造 增葬则蚤葬贼蚤燥灶 燥枣 泽燥蚤造 皂葬糟则燥鄄枣葬怎灶葬 糟燥皂皂怎灶蚤贼赠 泽贼则怎糟贼怎则藻 蚤灶 贼澡则藻藻 贼藻皂责藻则葬贼藻 枣燥则藻泽贼泽蕴陨 晕葬袁 在匀粤晕郧 载怎藻责蚤灶早袁 在匀粤晕郧 蕴蚤皂蚤灶 渊远圆猿远冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎悦燥皂皂怎灶蚤贼赠 泽贼则怎糟贼怎则藻 葬灶凿 泽责藻糟蚤藻泽 遭蚤燥凿蚤增藻则泽蚤贼赠 燥枣 枣蚤早 憎葬泽责泽 蚤灶 泽赠糟燥灶蚤葬 燥枣 云蚤糟怎泽 泽怎责藻则遭葬 酝蚤择援 增葬则援 躁葬责燥灶蚤糟葬 酝蚤择援 蚤灶 云怎扎澡燥怎悦匀耘晕 再燥怎造蚤灶早袁 悦匀耘晕 载蚤葬燥择蚤葬灶袁宰哉 宰藻灶泽澡葬灶袁藻贼 葬造 渊远圆源远冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎酝葬则蚤灶藻 藻糟燥造燥早蚤糟葬造 糟葬责蚤贼葬造院 增葬造怎葬贼蚤燥灶 皂藻贼澡燥凿泽 燥枣 皂葬则蚤灶藻 藻糟燥泽赠泽贼藻皂 泽藻则增蚤糟藻泽 悦匀耘晕 杂澡葬灶早袁砸耘晕 阅葬糟澡怎葬灶袁载陨粤 栽葬燥袁藻贼 葬造 渊远圆缘源冤噎噎郧藻燥皂燥则责澡燥造燥早蚤糟 则藻早蚤燥灶葬造蚤扎葬贼蚤燥灶 燥枣 悦澡蚤灶葬 葬蚤皂藻凿 葬贼 糟燥灶泽贼则怎糟贼蚤燥灶 燥枣 灶葬贼怎则藻 则藻泽藻则增藻 泽赠泽贼藻皂 郧哉韵 在蚤造蚤葬灶早袁 悦哉陨 郧怎燥枣葬 渊远圆远源冤噎噎噎噎噎陨皂责葬糟贼 燥枣 藻糟燥造燥早蚤糟葬造 增藻早藻贼葬贼蚤燥灶 糟燥灶泽贼则怎糟贼蚤燥灶 燥灶 贼澡藻 造葬灶凿泽糟葬责藻 责葬贼贼藻则灶 燥枣 葬 蕴燥藻泽泽 孕造葬贼藻葬怎 宰葬贼藻则泽澡藻凿再陨 再葬灶早袁 载陨晕 在澡燥灶早遭葬燥袁 匝陨晕 再怎灶遭蚤灶袁 藻贼 葬造 渊远圆苑苑冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎杂责葬贼蚤葬造 澡藻贼藻则燥早藻灶藻蚤贼赠 燥枣 泽燥蚤造 皂燥蚤泽贼怎则藻 葬糟则燥泽泽 葬 糟则燥责造葬灶凿鄄早则葬泽泽造葬灶凿 皂燥泽葬蚤糟院 葬 糟葬泽藻 泽贼怎凿赠 枣燥则 葬早则燥鄄责葬泽贼怎则葬造 贼则葬灶泽蚤贼蚤燥灶 蚤灶 灶燥则贼澡 燥枣悦澡蚤灶葬 宰粤晕郧 匀燥灶早皂藻蚤袁 宰粤晕郧 在澡燥灶早造蚤葬灶早袁 宰粤晕郧 运怎灶袁 藻贼 葬造 渊远圆愿苑冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎栽澡藻 则藻早蚤燥灶葬造 凿蚤增藻则泽蚤贼赠 燥枣 糟澡葬灶早藻泽 蚤灶 早则燥憎蚤灶早 凿怎则葬贼蚤燥灶 燥枣 泽责则蚤灶早 憎澡藻葬贼 葬灶凿 蚤贼泽 糟燥则则藻造葬贼蚤燥灶 憎蚤贼澡 糟造蚤皂葬贼蚤糟 葬凿葬责贼葬贼蚤燥灶 蚤灶 晕燥则贼澡藻则灶悦澡蚤灶葬 耘 再燥怎澡葬燥袁 匀哉韵 在澡蚤早怎燥袁酝粤 再怎责蚤灶早袁藻贼 葬造 渊远圆怨缘冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎砸藻泽责燥灶泽藻 燥枣 泽燥蚤造 责澡赠泽蚤糟葬造鄄糟澡藻皂蚤糟葬造 责则燥责藻则贼蚤藻泽 贼燥 则燥糟噪赠 凿藻泽藻则贼蚤枣蚤糟葬贼蚤燥灶 泽怎糟糟藻泽泽蚤燥灶 蚤灶 杂燥怎贼澡 悦澡蚤灶葬 运葬则泽贼杂匀耘晕郧 酝葬燥赠蚤灶袁 蕴陨哉 再葬灶早袁 载陨韵晕郧 运葬灶早灶蚤灶早 渊远猿园猿冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎孕则藻凿蚤糟贼蚤燥灶 燥枣 贼澡藻 藻枣枣藻糟贼泽 燥枣 糟造蚤皂葬贼藻 糟澡葬灶早藻 燥灶 贼澡藻 责燥贼藻灶贼蚤葬造 凿蚤泽贼则蚤遭怎贼蚤燥灶 燥枣 皂蚤则藻 蚤灶 晕燥则贼澡藻葬泽贼藻则灶 悦澡蚤灶葬匀耘 宰藻蚤袁月哉 砸藻灶糟葬灶早袁蕴陨哉 匀燥灶早躁怎葬灶袁藻贼 葬造 渊远猿员源冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎杂燥蚤造 灶蚤贼则燥早藻灶 皂蚤灶藻则葬造蚤扎葬贼蚤燥灶 葬灶凿 葬泽泽燥糟蚤葬贼藻凿 贼藻皂责藻则葬贼怎则藻 泽藻灶泽蚤贼蚤增蚤贼赠 燥枣 凿蚤枣枣藻则藻灶贼 陨灶灶藻则 酝燥灶早燥造蚤葬灶 早则葬泽泽造葬灶凿泽在匀哉 允蚤葬灶曾蚤灶早袁 宰粤晕郧 匝蚤怎枣藻灶早袁 匀耘 晕蚤葬灶责藻灶早袁 藻贼 葬造 渊远猿圆园冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎耘枣枣藻糟贼泽 燥枣 造葬灶凿 怎泽藻 燥灶 泽燥蚤造 灶怎贼则蚤藻灶贼 蚤灶 燥葬泽蚤泽鄄凿藻泽藻则贼 藻糟燥贼燥灶藻 蚤灶 贼澡藻 皂蚤凿凿造藻 则藻葬糟澡 燥枣 贼澡藻 匀藻蚤澡藻 砸蚤增藻则酝粤 在澡蚤皂蚤灶袁 蕴譈 再蚤澡藻袁 杂哉晕 云藻蚤曾蚤葬灶早袁 藻贼 葬造 渊远猿圆愿冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎粤泽泽藻泽泽皂藻灶贼 燥灶 澡藻葬增赠 皂藻贼葬造 责燥造造怎贼蚤燥灶 泽贼葬贼怎泽 蚤灶 责葬凿凿赠 泽燥蚤造泽 蚤灶 贼澡藻 灶燥则贼澡藻则灶 悦澡藻灶早凿怎 孕造葬蚤灶 葬灶凿 贼澡藻蚤则 责燥贼藻灶贼蚤葬造 藻糟燥造燥早蚤糟葬造 则蚤泽噪匝陨晕 再怎泽澡藻灶早袁 再哉 匀怎葬袁 云耘晕郧 宰藻灶择蚤葬灶早袁 藻贼 葬造 渊远猿猿缘冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎砸藻造葬贼蚤燥灶泽澡蚤责 遭藻贼憎藻藻灶 贼澡藻 贼藻皂责燥则葬造鄄泽责葬贼蚤葬造 凿蚤泽贼则蚤遭怎贼蚤燥灶 燥枣 造燥灶早造蚤灶藻 枣蚤泽澡蚤灶早 早则燥怎灶凿泽 燥枣 赠藻造造燥憎枣蚤灶 贼怎灶葬 渊栽澡怎灶灶怎泽 葬造遭葬糟葬则藻泽冤 葬灶凿 贼澡藻贼澡藻则皂燥糟造蚤灶藻 糟澡葬则葬糟贼藻则蚤泽贼蚤糟泽 蚤灶 贼澡藻 悦藻灶贼则葬造 粤贼造葬灶贼蚤糟 韵糟藻葬灶 再粤晕郧 杂澡藻灶早造燥灶早袁 酝粤 允怎灶躁蚤藻袁在匀粤晕郧 再怎袁 藻贼 葬造 渊远猿源缘冤噎噎噎噎噎月蚤燥造燥早蚤糟葬造 灶蚤贼则燥早藻灶 枣蚤曾葬贼蚤燥灶 蚤灶 贼澡藻 怎责责藻则 憎葬贼藻则 糟燥造怎皂灶 蚤灶 贼澡藻 泽燥怎贼澡 栽葬蚤憎葬灶 杂贼则葬蚤贼 凿怎则蚤灶早 泽怎皂皂藻则 圆园员员蕴陨晕 云藻灶早袁 悦匀耘晕 酝蚤灶袁 再粤晕郧 宰藻蚤枣藻灶早袁 藻贼 葬造 渊远猿缘源冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎杂贼燥则葬早藻 葬灶凿 凿则蚤增藻则泽 燥枣 枣燥则藻泽贼泽 糟葬则遭燥灶 燥灶 贼澡藻 月藻蚤糟澡葬灶早泽澡葬灶 陨泽造葬灶凿 燥枣 酝蚤葬燥凿葬燥 粤则糟澡蚤责藻造葬早燥杂匀陨 匀燥灶早澡怎葬袁宰粤晕郧 载蚤葬燥造蚤袁宰粤晕郧 粤蚤袁藻贼 葬造 渊远猿远猿冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎陨皂责葬糟贼 燥枣 糟澡葬灶早藻泽 蚤灶 增藻早藻贼葬贼蚤燥灶 贼赠责藻泽 燥灶 泽燥蚤造 悦 皂蚤灶藻则葬造蚤扎葬贼蚤燥灶 葬灶凿 葬泽泽燥糟蚤葬贼藻凿 贼藻皂责藻则葬贼怎则藻 泽藻灶泽蚤贼蚤增蚤贼赠 蚤灶 贼澡藻 悦澡葬灶早遭葬蚤 酝燥怎灶贼葬蚤灶枣燥则藻泽贼泽 燥枣 悦澡蚤灶葬 宰粤晕郧 阅葬灶袁 蕴譈 再怎造蚤葬灶早袁 载哉 蕴蚤袁 藻贼 葬造 渊远猿苑猿冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎粤灶葬造赠泽蚤泽 燥枣 则藻造葬贼蚤燥灶泽澡蚤责 遭藻贼憎藻藻灶 早藻灶藻贼蚤糟 泽贼则怎糟贼怎则藻 燥枣 悦澡蚤灶藻泽藻 孕蚤灶藻 葬灶凿 皂燥怎灶贼葬蚤灶 遭葬则则蚤藻则泽酝耘晕郧 载蚤葬灶早曾蚤葬灶早袁 阅陨 载蚤葬燥赠葬灶袁 宰粤晕郧 酝藻灶早遭藻灶袁 藻贼 葬造 渊远猿愿圆冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎杂燥蚤造 燥则早葬灶蚤糟 糟葬则遭燥灶 蚤灶贼藻则责燥造葬贼蚤燥灶 遭葬泽藻凿 燥灶 葬怎曾蚤造蚤葬则赠 藻灶增蚤则燥灶皂藻灶贼葬造 糟燥增葬则蚤葬贼藻泽院葬 糟葬泽藻 泽贼怎凿赠 葬贼 泽皂葬造造 憎葬贼藻则泽澡藻凿 泽糟葬造藻 蚤灶 蕴燥藻泽泽 匀蚤造造赠则藻早蚤燥灶 宰耘晕 宰藻灶袁 在匀韵哉 月葬燥贼燥灶早袁 宰粤晕郧 再葬枣藻灶早袁 藻贼 葬造 渊远猿愿怨冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎耘糟燥鄄皂葬灶葬早藻皂藻灶贼 遭藻灶藻枣蚤贼 葬灶葬造赠泽蚤泽 燥枣 蚤灶凿怎泽贼则蚤葬造 则藻泽燥怎则糟藻泽 枣则燥皂 造蚤枣藻 糟赠糟造藻 责藻则泽责藻糟贼蚤增藻院葬 糟葬泽藻 泽贼怎凿赠 燥枣 葬 增蚤则贼怎葬造 泽赠皂遭蚤燥泽蚤泽 灶藻贼憎燥则噪杂匀陨 载蚤葬燥择蚤灶早袁蕴陨 载蚤葬燥灶怎燥袁再粤晕郧 允蚤葬灶曾蚤灶 渊远猿怨愿冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎栽澡藻 早葬皂藻 葬灶葬造赠泽蚤泽 遭藻贼憎藻藻灶 责燥增藻则贼赠 葬灶凿 藻灶增蚤则燥灶皂藻灶贼 蚤灶 藻糟燥造燥早蚤糟葬造造赠 枣则葬早蚤造藻 扎燥灶藻泽 匝陨 载蚤灶澡怎葬袁再耘 杂澡蚤造蚤灶袁悦匀耘晕郧 再怎袁 藻贼 葬造 渊远源员员冤噎栽澡藻 糟燥怎责造蚤灶早 凿藻增藻造燥责皂藻灶贼 燥枣 藻糟燥灶燥皂赠 葬灶凿 藻灶增蚤则燥灶皂藻灶贼 怎灶凿藻则 贼澡藻 遭葬糟噪早则燥怎灶凿 燥枣 宰燥则造凿 耘曾责燥 蚤灶 杂澡葬灶早澡葬蚤晕陨 再葬燥袁 再哉耘 宰藻灶扎藻袁 在匀粤晕郧 再怎灶贼葬灶早袁 藻贼 葬造 渊远源员愿冤噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎噎
远圆源远 摇 生摇 态摇 学摇 报摇 摇 摇 猿猿卷摇
叶生态学报曳圆园员猿年征订启事
叶生态学报曳是由中国科学技术协会主管袁中国生态学学会尧中国科学院生态环境研究中心主办的生态学
高级专业学术期刊袁创刊于 员怨愿员年袁报道生态学领域前沿理论和原始创新性研究成果遥 坚持野百花齐放袁百家
争鸣冶的方针袁依靠和团结广大生态学科研工作者袁探索生态学奥秘袁为生态学基础理论研究搭建交流平台袁
促进生态学研究深入发展袁为我国培养和造就生态学科研人才和知识创新服务尧为国民经济建设和发展服务遥
叶生态学报曳主要报道生态学及各分支学科的重要基础理论和应用研究的原始创新性科研成果遥 特别欢
迎能反映现代生态学发展方向的优秀综述性文章曰研究简报曰生态学新理论尧新方法尧新技术介绍曰新书评价和
学术尧科研动态及开放实验室介绍等遥
叶生态学报曳为半月刊袁大 员远开本袁猿园园页袁国内定价 怨园元 辕册袁全年定价 圆员远园元遥
国内邮发代号院愿圆鄄苑袁国外邮发代号院酝远苑园
标准刊号院陨杂杂晕 员园园园鄄园怨猿猿摇 摇 悦晕 员员鄄圆园猿员 辕 匝
全国各地邮局均可订阅袁也可直接与编辑部联系购买遥 欢迎广大科技工作者尧科研单位尧高等院校尧图书
馆等订阅遥
通讯地址院 员园园园愿缘 北京海淀区双清路 员愿号摇 电摇 摇 话院 渊园员园冤远圆怨源员园怨怨曰 远圆愿源猿猿远圆
耘鄄皂葬蚤造院 泽澡藻灶早贼葬蚤曾怎藻遭葬燥岳 则糟藻藻泽援葬糟援糟灶摇 网摇 摇 址院 憎憎憎援藻糟燥造燥早蚤糟葬援糟灶
本期责任副主编摇 陈利顶摇 摇 摇 编辑部主任摇 孔红梅摇 摇 摇 执行编辑摇 刘天星摇 段摇 靖
生摇 态摇 学摇 报渊杂匀耘晕郧栽粤陨摇 载哉耘月粤韵冤渊半月刊摇 员怨愿员年 猿月创刊冤
第 猿猿卷摇 第 员怨期摇 渊圆园员猿年 员园月冤
粤悦栽粤 耘悦韵蕴韵郧陨悦粤 杂陨晕陨悦粤摇渊杂藻皂蚤皂燥灶贼澡造赠袁杂贼葬则贼藻凿 蚤灶 员怨愿员冤摇灾燥造郾 猿猿摇 晕燥郾 员怨 渊韵糟贼燥遭藻则袁 圆园员猿冤
编摇 摇 辑摇 叶生态学报曳编辑部
地址院北京海淀区双清路 员愿号
邮政编码院员园园园愿缘
电话院渊园员园冤远圆怨源员园怨怨憎憎憎援藻糟燥造燥早蚤糟葬援糟灶泽澡藻灶早贼葬蚤曾怎藻遭葬燥岳 则糟藻藻泽援葬糟援糟灶
主摇 摇 编摇 王如松
主摇 摇 管摇 中国科学技术协会
主摇 摇 办摇 中国生态学学会
中国科学院生态环境研究中心
地址院北京海淀区双清路 员愿号
邮政编码院员园园园愿缘
出摇 摇 版摇
摇 摇 摇 摇 摇 地址院北京东黄城根北街 员远号
邮政编码院员园园苑员苑
印摇 摇 刷摇 北京北林印刷厂
发 行摇
地址院东黄城根北街 员远号
邮政编码院员园园苑员苑
电话院渊园员园冤远源园猿源缘远猿耘鄄皂葬蚤造院躁燥怎则灶葬造岳 糟泽责早援灶藻贼
订摇 摇 购摇 全国各地邮局
国外发行摇 中国国际图书贸易总公司
地址院北京 猿怨怨信箱
邮政编码院员园园园源源
广告经营
许 可 证摇 京海工商广字第 愿园员猿号
耘凿蚤贼藻凿 遭赠摇 耘凿蚤贼燥则蚤葬造 遭燥葬则凿 燥枣
粤悦栽粤 耘悦韵蕴韵郧陨悦粤 杂陨晕陨悦粤
粤凿凿院员愿袁杂澡怎葬灶早择蚤灶早 杂贼则藻藻贼袁匀葬蚤凿蚤葬灶袁月藻蚤躁蚤灶早 员园园园愿缘袁悦澡蚤灶葬
栽藻造院渊园员园冤远圆怨源员园怨怨
憎憎憎援藻糟燥造燥早蚤糟葬援糟灶
泽澡藻灶早贼葬蚤曾怎藻遭葬燥岳 则糟藻藻泽援葬糟援糟灶
耘凿蚤贼燥则鄄蚤灶鄄糟澡蚤藻枣摇 宰粤晕郧 砸怎泽燥灶早
杂怎责藻则增蚤泽藻凿 遭赠摇 悦澡蚤灶葬 粤泽泽燥糟蚤葬贼蚤燥灶 枣燥则 杂糟蚤藻灶糟藻 葬灶凿 栽藻糟澡灶燥造燥早赠
杂责燥灶泽燥则藻凿 遭赠摇 耘糟燥造燥早蚤糟葬造 杂燥糟蚤藻贼赠 燥枣 悦澡蚤灶葬
砸藻泽藻葬则糟澡 悦藻灶贼藻则 枣燥则 耘糟燥鄄藻灶增蚤则燥灶皂藻灶贼葬造 杂糟蚤藻灶糟藻泽袁 悦粤杂
粤凿凿院员愿袁杂澡怎葬灶早择蚤灶早 杂贼则藻藻贼袁匀葬蚤凿蚤葬灶袁月藻蚤躁蚤灶早 员园园园愿缘袁悦澡蚤灶葬
孕怎遭造蚤泽澡藻凿 遭赠摇 杂糟蚤藻灶糟藻 孕则藻泽泽
粤凿凿院员远 阅燥灶早澡怎葬灶早糟澡藻灶早早藻灶 晕燥则贼澡 杂贼则藻藻贼袁
月藻蚤躁蚤灶早摇 员园园苑员苑袁悦澡蚤灶葬
孕则蚤灶贼藻凿 遭赠摇 月藻蚤躁蚤灶早 月藻蚤 蕴蚤灶 孕则蚤灶贼蚤灶早 匀燥怎泽藻袁
月藻蚤躁蚤灶早 员园园园愿猿袁悦澡蚤灶葬
阅蚤泽贼则蚤遭怎贼藻凿 遭赠摇 杂糟蚤藻灶糟藻 孕则藻泽泽
粤凿凿院员远 阅燥灶早澡怎葬灶早糟澡藻灶早早藻灶 晕燥则贼澡
杂贼则藻藻贼袁月藻蚤躁蚤灶早 员园园苑员苑袁悦澡蚤灶葬
栽藻造院渊园员园冤远源园猿源缘远猿
耘鄄皂葬蚤造院躁燥怎则灶葬造岳 糟泽责早援灶藻贼
阅燥皂藻泽贼蚤糟 摇 摇 粤造造 蕴燥糟葬造 孕燥泽贼 韵枣枣蚤糟藻泽 蚤灶 悦澡蚤灶葬
云燥则藻蚤早灶 摇 摇 悦澡蚤灶葬 陨灶贼藻则灶葬贼蚤燥灶葬造 月燥燥噪 栽则葬凿蚤灶早
悦燥则责燥则葬贼蚤燥灶
粤凿凿院孕援韵援月燥曾 猿怨怨 月藻蚤躁蚤灶早 员园园园源源袁悦澡蚤灶葬
摇 陨杂杂晕 员园园园鄄园怨猿猿悦晕 员员鄄圆园猿员 辕 匝 国内外公开发行 国内邮发代号 愿圆鄄苑 国外发行代号 酝远苑园 定价 怨园郾 园园元摇