全 文 ：
　 　 　 　 　 生 态 学 报
　 　 　 　 　 　 　 （ＳＨＥＮＧＴＡＩ ＸＵＥＢＡＯ）
　 　 第 ３４卷 第 ４期　 　 ２０１４年 ２月　 （半月刊）
目　 　 次
前沿理论与学科综述
富营养化湖泊溶解性有机碳生物可利用性研究进展 叶琳琳，孔繁翔，史小丽，等 （７７９）………………………
黄河下游平原农业景观中非农生境植物多样性 卢训令，梁国付，汤　 茜，等 （７８９）……………………………
个体与基础生态
锰胁迫对杠板归细胞超微结构的影响 王　 钧，邬　 卉，薛生国，等 （７９８）………………………………………
不同渗氧能力水稻品种对砷的耐性和积累 吴　 川，莫竞瑜，薛生国，等 （８０７）…………………………………
弱光下水分胁迫对不同产地披针叶茴香幼苗生理特性的影响 曹永慧，周本智，陈双林 （８１４）…………………
不同分枝数对桑树幼苗生长发育的影响 郇慧慧，胥　 晓，刘　 刚，等 （８２３）……………………………………
斑膜合垫盲蝽若虫在国槐上的空间分布型及抽样技术 朱惠英，沈　 平，吴建华，等 （８３２）……………………
连作苹果园土壤真菌的 Ｔ⁃ＲＦＬＰ 分析 尹承苗，王功帅，李园园，等 （８３７）………………………………………
棉隆对苹果连作土壤微生物及平邑甜茶幼苗生长的影响 刘恩太，李园园，胡艳丽，等 （８４７）…………………
两株具有芘降解功能的植物内生细菌的分离筛选及其特性 孙　 凯，刘　 娟，李　 欣，等 （８５３）………………
种群、群落和生态系统
温度对柑橘始叶螨实验种群生长发育繁殖的影响 李迎洁，王梓英，张国豪，等 （８６２）…………………………
高原鼠兔有效洞穴密度对青藏高原高寒草甸群落植物生态位的影响 贾婷婷，毛　 亮，郭正刚 （８６９）…………
三工河流域琵琶柴群落特征与土壤因子的相关分析 赵学春，来利明，朱林海，等 （８７８）………………………
岷江干旱河谷造林对土壤微生物群落结构的影响 王卫霞，罗　 达，史作民，等 （８９０）…………………………
滩涂围垦和土地利用对土壤微生物群落的影响 林　 黎，崔　 军，陈学萍，等 （８９９）……………………………
福寿螺对稻田水生植物群落结构的影响 赵本良，章家恩，戴晓燕，等 （９０７）……………………………………
４种木本植物在潜流人工湿地环境下的适应性与去污效果 陈永华，吴晓芙，郝　 君，等 （９１６）………………
基于静态箱式法和生物量评估海北金露梅灌丛草甸碳收支 李红琴，李英年，张法伟，等 （９２５）………………
初始 ｐＨ值对碱性和酸性水稻土微生物铁还原过程的影响 吴　 超，曲　 东，刘　 浩 （９３３）……………………
景观、区域和全球生态
库姆塔格柽柳沙包年层稳定碳同位素与气候环境变化 张锦春，姚　 拓，刘长仲，等 （９４３）……………………
资源与产业生态
大棚甜瓜蒸腾规律及其影响因子 张大龙，常毅博，李建明，等 （９５３）……………………………………………
盐胁迫下荒漠共生植物红砂与珍珠的根茎叶中离子吸收与分配特征 赵　 昕，杨小菊，石　 勇，等 （９６３）……
普通鹿蹄草品质与根际和非根际土壤的关系 耿增超，孟令军，刘建军 （９７３）……………………………………
作物种植前后土壤有机质及养分因子的空间变异分析 方　 斌，吴金凤 （９８３）…………………………………
城乡与社会生态
城市河流健康评价指标体系构建及其应用 邓晓军，许有鹏，翟禄新，等 （９９３）…………………………………
西藏生态足迹与承载力动态分析 安宝晟，程国栋 （１００２）…………………………………………………………
研究简报
三峡库区岸坡消落带草地、弃耕地和耕地土壤微生物及酶活性特征 马　 朋，李昌晓，雷　 明，等 （１０１０）……
盐胁迫对 ２种栎树苗期生长和根系生长发育的影响 王树凤，胡韵雪，孙海菁，等 （１０２１）………………………
恒温和变温驯化对大蟾蜍蝌蚪热耐受性的影响 王立志 （１０３０）…………………………………………………
学术信息与动态
国际生物土壤结皮研究发展态势文献计量分析 贺郝钰，侯春梅，迟秀丽，等 （１０３５）……………………………
期刊基本参数：ＣＮ １１⁃２０３１ ／ Ｑ∗１９８１∗ｍ∗１６∗２６４∗ｚｈ∗Ｐ∗ ￥ ９０􀆰 ００∗１５１０∗３０∗２０１４⁃０２
􀥒􀥒􀥒􀥒􀥒􀥒􀥒􀥒􀥒􀥒􀥒􀥒􀥒􀥒
封面图说： 大蟾蜍蝌蚪群———大蟾蜍别名癞蛤蟆，体长达 １０ｃｍ以上，身体肥胖，四肢短，步态及齐足跳的姿势具特征性。 其背
部皮肤厚而干燥，通常有疣，呈黑绿色，常有褐色花斑，趾间具蹼。 毒腺在背部的疣内，受惊后毒腺分泌或射出毒液。
大蟾蜍早春在水中繁殖，可迁移至 １．５ｋｍ外或更远的适合繁殖的池塘，产卵量很大，产卵数天后蝌蚪即可孵出，１—３
个月后发育为蟾。 大蟾蜍常作为实验动物或药用动物，其耳后腺和皮肤腺的白色分泌物可制成“蟾酥”，可治疗多种
疾病。 研究表明，大蟾蜍蝌蚪最高逃避温度和最高致死温度比最适温度产生的影响要大。
彩图及图说提供： 陈建伟教授　 北京林业大学　 Ｅ⁃ｍａｉｌ： ｃｉｔｅｓ．ｃｈｅｎｊｗ＠ １６３．ｃｏｍ
第 34 卷第 4 期
2014年 2月
生 态 学 报
ACTA ECOLOGICA SINICA
Vol.34,No.4
Feb.,2014
http: / / www.ecologica.cn
基金项目:国家自然基金资助项目(40971105, 41271189);江苏省高校哲学社会科学研究重点资助项目(2010ZDIXM049);江苏高校优势学科建
设工程资助项目
收稿日期:2012鄄09鄄30; 摇 摇 修订日期:2013鄄04鄄08
*通讯作者 Corresponding author.E鄄mail: wenyanfang731@ 163.com
DOI: 10.5846 / stxb201209301364
方斌, 吴金凤.作物种植前后土壤有机质及养分因子的空间变异分析.生态学报,2014,34(4):983鄄992.
Fang B, Wu J F.Spatial variation analysis of soil organic matter and nutrient factor for before and after planting crops.Acta Ecologica Sinica,2014,34(4):
983鄄992.
作物种植前后土壤有机质及养分因子的空间变异分析
方摇 斌1,*,吴金凤1,2
(1. 南京师范大学地理科学学院, 南京摇 210023; 2. 中国科学院地理与资源研究所, 北京摇 100101)
摘要:以作物土壤试验数据为基础,结合多年的农户调查数据,采用 Moran忆s I指数分析、Pearson相关系数分析和 Kriging插值 3
种地统计学方法,探讨了浙江省浦江县作物种植前后土壤有机质及氮磷钾等养分因子的空间演化规律。 研究发现 Moran忆s I指
数分析证实随着区域范围的加大,各因子间自相关程度产生较大差异。 这一点在 Pearson 相关系数分析中也得到证实,并在
kriging插值图中直观表达。 结果还表明:有机质和碱解氮的形成与水的关系较为密切,水源区两因子的含量较周边高,且会促
进前者向后者转化。 养分投入后,过多的养分在土壤中降解过程与区域自然环境条件的关系较密切。 同一作物的农户养分投
入在不同区域差异较大,不同作物的农户养分投入在区域内却存在量的相关。
关键词:作物种植前后;土壤养分;空间变异分析
Spatial variation analysis of soil organic matter and nutrient factor for before and
after planting crops
FANG Bin1,*, WU Jinfeng1,2
1 College of Geography, Nanjing Normal University, Nanjing 210023, China
2 Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China
Abstract: The soil organic matter and nutrient factor space evolution for before and after planting crops was discussed,
based on data of tested soil and household survey over the years. Results could provide the theoretical basis for reasonable
input of crop nutrient and planting structure adjustment to divisional guidance in different regional.
The spatial variability of soil factor before and harvest after crops was analyzed by using Moran忆s I index analysis,
Pearson correlation coefficient analysis, and Kriging interpolation diagram. Results suggested that three methods can
interactly examine and complemente. Moran忆s I index was proved to the great different of nutrient factor of all kinds with the
enlarged area coverage,. The correlation degree among all kinds was tested through Pearson correlation. The feature of space
distribution in different area could be directly reacted by Kriging interpolation..
The characteristics of used nutrient was showed as following: 1) The soil organic matter and nitrogen was a very close
relationship with water, and they had a high content in the water resource district. Results resulted from the accelerated the
transformation of nitrogen by water. 2) The a period of time was needed to degradation of a plethora of nutrients in the soil
after input nutrient, and this process was closely related with the regional natural environment condition. 3) The calculated
result showed a significant correlation between the content of P and K Because of mainly using the compound fertilizer.The
content of P and K of soil were higher after harvest crops than before planting crops, which was contributed to the used
compound fertilizer with fixed ratio of P and K by farmers. In addition, another reasons is the slower degradation rate of P
http: / / www.ecologica.cn
and K. For the same crop, nutrient inputs by farmers was differend in different regional, but it was correlations in
quantitative in the area for different crop. Input nutrient in the same area was the relatively small differences between in
conventional rice crop and in scale management grape. However, input nutrient was bigger different for other crops. Spatial
analysis results was basically agreed with the peasant household survey results. Nitrogen input should be properly controlled
in high mountain, suburban vegetable growing areas, and pot source river downstream zone due to high input nitrogen. The
high content of soil P was also properly needed to control in Mountain surrounding. The content of input K is usually not
high. Although the the K content in some areas is higher, it is mainly distributed in the high mountains areas, which may be
related to the trees.
Key Words: before and after planting crops ; soil nutrients; spatial variation analysis
摇 摇 土壤养分的空间变异规律可作为土壤资源合理
利用与改良、精准施肥和环境物质修复的基本依
据[1鄄3]。 土壤的变异来源可分为系统变异和随机变
异, 其中前者可认为是由母质、气候、水文、地形、生
物等相对稳定的因素产生,后者可认为是由不同土
地利用方式和行为,以及取样、分析等的误差产
生[4鄄6]。 大量的研究表明地统计学方法能较好地揭
示土壤异质性机理[7鄄11]。 近几年, 国内外学者对不
同尺度( 乡镇级、村级、田间尺度和微尺度) 的土壤
性质( 土壤水分、盐分、养分、理化性质及生物学性
质等)进行了广泛研究, 如,连纲等通过对黄土高原
县域土壤养分空间变异分析后,认为不同土地利用
类型的养分含量差异显著,水田和川地的有机质和
全氮含量较高,全磷含量以梯田最高。 高义民等通
过对黄土区村级农田土壤养分空间变异特征分析,
认为施肥是导致该村土壤养分空间变异性的主要原
因。 张博等通过对江苏滩涂围垦区土壤养分空间变
异分析得出该区域应多施有机肥和氮肥,少施或不施
钾肥和磷肥的结论。 这些研究都在一定程度上揭示
了土壤有机质、水解氮、速效磷、速效钾等土壤因子在
各自环境区内的空间变异规律[12鄄19]。 在方法论的研
究上,GIS技术被引入进行景观空间异质性分析,并将
景观斑块结构、形态、数量的变化特征,作为景观连通
性、稳定性和破碎化的重要依据[20鄄24]。 但是,运用地
统计学方法探讨作物种植前后土壤养分空间变异性
的文献则鲜有发现。 本文以浦江县作物土壤试验数
据为基础,结合多年的农户调查数据,以地统计学方
法为工具,探讨作物种植前后土壤有机质及其养分因
子的空间演化规律,为不同区域作物养分合理实施与
种植结构的调整进行分区指导提供理论依据。
1摇 研究对象与方法
图 1摇 浦江县的地理特征
Fig.1摇 Geographical Feature of PuJiang county
1.1摇 研究区简介
浙江省浦江县位于浙江省中部偏西,金衢盆地
的北缘,是浦阳、壶源两江之发源地。 地理坐标界于
东经 119毅42忆—120毅 07忆、北纬 29毅21忆—29毅41忆。 东西
宽 39. 25 km,南北长 36. 5 km。 全县土地面积为
907郾 6637km2。 境内地貌分西北山区和浦江盆地两
大类型。 西北部山区属浙中丘陵山地地貌,以低、高
丘为主。 浦江盆地属于金衢丘陵盆地地貌。 以平
畈、岗地、低丘为主。 根据海拔高度和农业利用特点
可将浦江的地形分这平原区(海拔<150 m)、半山区
(150 m<海拔< 500 m)、中低山区(海拔>500 m),符
合“七山二水一分田冶中国东部地形特征。
浦江县属亚热带季风气候区,四季分明,气温适
489 摇 生摇 态摇 学摇 报摇 摇 摇 34卷摇
http: / / www.ecologica.cn
中,光照充足,雨量丰沛,自然资源丰富。 由于农业
用地分别邻近浦阳江、壶源江和大陈江,因而又形成
3种小气候带,分别是浦江盆地温暖少雨区(玉区)、
壶源江谷温和次多雨区(域区)、中低山温凉多雨区
(芋区)(见图 1)。 各气候带特征见表 1。
浦江县区域经济发展水平略高于全国平均水
平,随着产业扩张速度的不断加快,耕地保护的压力
也面临前所未有的挑战如何优化区域资源是亟待解
决的重要问题。 本研究正是在欧盟项目 IRMLA (综
合农业资源管理与土地利用分析)成果的基础上,以
正在执行的国家基金为支助,以连续 10a 的农户土
地利用与管理的调查数据和土壤、作物的试验数据
为基础,探讨县域农地资源的优化利用。
表 1摇 浦江各农业气候带的主要土壤特征
Table 1摇 The main soil characteristics of each agricultural climate in PuJiang
编号
No.
海拔
EL / m
年平均温度
MAT / 益
年降雨量
ARF / mm
逸10 益的积温
AT逸10益
土壤特性 Soil characteristics
种类 T 质地 CT 表层SL / cm pH
有机质含量
OM / %
I 70—300 16—17 1250—1400 4900—5200 水稻土 CL 17—19 6.85 2.61
域 28—500 15—16 1400—1500 4600—4900 水稻土 C 15—17 5.80 2.42
芋 500—1200 13.5—15 1500—1577 3900—4600 黄 壤 CL 14—17 5.23 4.49
摇 摇 C:粘性 Cohesive soil;L:壤性 Loam soil;CL:粘壤性 Cohesive and Loam soil;EL:海拔 Elevation;MAT:年平均温度 Mean annual temperature;ARF:
年降雨量 Annual rainfall;AT: 积温 Accumulated temperature;T:类型 Type;CT:质地 Character;SL:表层 Surface layer,;OM:有机质 Organic matter
1.2摇 主要研究方法
1.2.1摇 土地利用与管理调查
(1)调查方法
在浦江县选取了 107 个代表性农户进行调查,
农户的选取采取点、线、面相结合的方式。 面:尽量
做到覆盖面广。 将全县分 3 片:平原、半山区、中低
山区,并兼顾到全县的每个乡镇。 同时把农作物种
类和种植面积、各地气候差异、典型产品和重点产业
相结合综合考虑,尽量做到平衡。 线:以公路为一
线,池塘、沟渠为二线,山间小路为三线。 点:除考虑
各种不同的作物外,还兼顾好、中、差 3 类生产水平
的农户。 具体调查点分布见图 1。 在调查人员的安
排上,采取乡镇从事农业技术推广专业人员和在校
教师、研究生相结合的方式,根据选取的农户、统一
标准、逐户调查。
摇 摇 (2)调查内容
在调查之前先设计详细的调查表,与本研究相
关的主要内容包括:1)农户家庭土地利用方式与农
业生产管理,包括养分投入,特别是 N 投入;2)农户
耕种田块的基本信息,如:种植面积、规模化程度、地
形条件、地理位置、土壤性质等,为土地利用方式的
空间分析奠定基础。
1.2.2摇 土壤采样与分析方法
(1)采样点的选取及农户投入
根据种植作物的不同从 107 户农户家庭的责任
田中取样 183个样点(即 183块田的采样),其中,水
稻及其它粮油作物共 69块, 葡萄 36块,蔬菜 40块,
桃形李及其它水果类 31 块,其它作物 7 块,各类作
物农户平均投入调查结果见表 2。 采样地块还需考
虑坡度、坡向、相邻性等特征,并做好地理特征的相
关信息采集。 各样点见图 1。
表 2 摇 不同作物 N、P、K农户投入均值
Table 2摇 Farmers mean input of N, P, K in different crops
养分
Nutrient /
(kg / hm2)
粮油作物
Grain and oil crops
水稻
Rice
油菜
Oilseed
高粱
Sorghum
玉米
Maize
小麦
Wheat
葡萄及其他水果
Grapes and other fruits
葡萄
Grap
西瓜 W
Watermelon
甘蔗
Sugarcane
梨
Pear
桃形李
Plum like
peach
蔬菜
Vegetables
蔬菜
Vegetable
N 200 100 45 145 210 430 387 290 350 130 450
P 18 11 60 20 40 105 189 120 180 30 210
K 60 20 25 65 25 200 178 125 154 35 160
摇 摇 N: 碱解氮,P: 速效磷,K: 速效钾
589摇 4期 摇 摇 摇 方斌摇 等:作物种植前后土壤有机质及养分因子的空间变异分析 摇
http: / / www.ecologica.cn
摇 摇 (2) 土壤取样方法
1)采样时间摇 第 1次取样为作物播种前(3 月 3
日—12日,共 10 d),第 2次取样为作物收获后、施肥
前(11月 25日到 12月 4 日共 10 d),这一取样时间
是综合了该区域的季节特征和农事活动的特点选择
的。 两次取样的田块完全相同。
2)采样部位和深度 摇 根据耕层厚度,确定采样
深度,一般取样深度 0—20 cm。
3)采样方法、数量 摇 采用多点混合土样采集方
法,每个混合农化土样由 20 个样点组成。 每个点的
取土深度及重量尽量均匀一致,土样上层和下层的
比例基本相同。 样品处理、储存等过程没接触金属
器具和橡胶制品。 每个混合样品一般取 1kg左右。
(3) 试验分析方法
样品主要参照鲁如坤的《土壤农业化学分析方
法》中的扩散吸收法对土壤有效 N 进行测定[25]。 本
文所指的氮、磷、钾分别是指碱解氮、速效磷和速
效钾。
1.2.3摇 地统计学分析
(1)空间自相关分析可检验土壤养分空间相关
性的存在与否(即是否存在区域的集聚或集中),一
般用全域 Moran忆s I 指数、局域 Moran忆s I 指数表达。
本文利用 ArcGIS空间统计模块中的全域 Moran忆s I,
对县域内两个时点因子空间自相关性进行分析。 计
算模型已在先前发表的论文中作了详细阐述,不再
赘述[26]。
(2)半方差函数是地统计学的基石,其块金值、
基台值和变程等参数能表达区域化变量在一定尺度
上的空间变异和相关程度,是 Kriging 精确插值的基
础,也行预测和模拟的关键[1鄄10]。
模型是否最优,主要看拟合的理论模型与实际
变异函数的残差平方和 ( )RSS 和决定系数 R( )2 的
大小,其中,残差平方和是实际变异函数值 酌 h( )i 和
理论模型变异函数值 酌^ h( )i 之和:
RSS =移
n
i = 1
酌 h( )i - 酌^ h( )[ ]i
2
(3)
式中, RSS是回归精度的重要参数,其取值愈小,说
明实际观测值与回归线靠近,拟合曲线与实际配合
愈好。
但是残差平方和受因变量和自变量绝对值大小
的影响, RSS 值大的拟合模型不一定不及 RSS 值小
的拟合模型,但是决定系数 R2 就不存在这个问题。
决定系数 R2,也称拟合度,其计算公式为:
R2 =
移
n
i = 1
酌^ h( )i - 酌 ( )[ ]h 2
移
n
i = 1
酌 ( )h - 酌^ h( )[ ]i 2
(4)
式中, 酌 h( )i 为实际变异函数值, 酌^ h( )i 为理论模型
计算出的变异函数值, 酌 ( )h 为实际变异函数平
均值。
3摇 研究结果与分析
本文利用单样本 K鄄S检验对各样本因子分别进
行检验,发现通过对数转换后各因子均较好地符合
正态分布。 可以进行 Moran忆s I 指数和 Kriging 插值
分析。
3.1摇 土壤因子空间自关性分析
不同养分在不同时点的空间自相关性与步长的
总变化趋势基本一致,在作物种植前后均呈现波动
性,说明随着尺度的变化,土壤养分之间存在一定的
空间关联,但又各具特点(这一点与高义民、史利江
等学者得出的结果基本一致[7,12,21] ):(1)各因子的
振幅在步长﹤ 1000 m时均达到最高点,并随着步长
的加大而逐渐衰减,说明区域越近,作物种植结构越
相似,养分投入越接近。 而随着步长的加大,范围的
加大,种植多样性,引起投入的多样化,但大范围的
多样化投入又引致相关性程度加大。 (2)作物收获
后,土壤氮和磷的负振幅超过作物种植前,说明不同
农户作物养分投入的差异,导致的相关性增强是由
于随机变异引起的,而经过一段时间土壤表层的氮、
磷含量趋于稳定的原因是由于土壤本身的吸附与消
化能力作用的结果,它属系统变异。 (3)作物收获后
钾的自相关强度明显弱于作物种植前,而且,随着步
长的加大,其相关性增强。 前者与不同区域人们对
钾肥投入的重视程度不一致有关,后者则说明钾的
自相关性与土壤因子系统变异的关系更加密切。 具
体见图 2。
2.2摇 土壤因子交互相关性分析
为分析土壤因子间的交互相关性,论文运用
SPSS软件进行 Pearson 相关系数分析,并在 0.01 的
置信水平下进行 T 检验。 得出:(1)作物种植前,在
置信0 .01水平下,碱解氮 鄄1与有机质 鄄1有较强相
689 摇 生摇 态摇 学摇 报摇 摇 摇 34卷摇
http: / / www.ecologica.cn
图 2摇 作物种植前后 Moran忆s I指数随步长的变化对照
Fig.2摇 Moran忆s I index after and before planting crops with step change
*-1、-2分别为作物种植前和收获后的土壤取样时间
关性,这一点也与史利江等学者得出的结果相一
致[7]。 而碱解氮鄄 2 与有机质鄄 1 之间的相关性较前
两者明显减弱。 说明有机质的含量与碱解氮间的关
系紧密,但前者转化为后者需要一定的时间。 在
0郾 05置信水平下,碱解氮鄄 2与速效磷鄄 2间存在一定
的相关性,与速效磷鄄 1、速效钾之间基本上不相关,
说明氮、磷投入,农户有配施的意识,但对钾的配施
意识不强,这与农户调查结果相符合。 (2)在 0.01
置信水平下,速效磷鄄 1与速效磷鄄 2之间存在极显著
的相关性,说明农户对磷的投入相对稳定且其投入
磷与土壤内磷素间成较为固定的比例。 (3)速效磷
与速效钾之间在种植前后都具很强的相关性,这可
能与调查中 65.4%的农户喜欢使用复合肥有关。 具
体见表 3。
表 3摇 作物种植前后土壤养分交互相关特征
Table 3摇 Correlation characteristic between different soil nutrient in the after and before planting crops
有机质鄄 1
OM鄄1
有机质鄄 2
OM鄄2
碱解氮鄄 1
N鄄1
碱解氮鄄 2
N鄄2
速效磷鄄 1
P鄄1
速效磷鄄 2
P鄄2
速效钾鄄 1
K鄄1
速效钾鄄 2
K鄄2
有机质鄄 1 OM鄄1 1
有机质鄄 2 OM鄄2 0.208** 1
碱解氮鄄 1 N鄄1 0.534** 0.149* 1
碱解氮鄄 2 N鄄2 0.203** 0.081 0.329** 1
速效磷鄄 1 P鄄1 0.178* 0.074 0.040 0.090 1
速效磷鄄 2 P鄄2 0.137 0.055 0.081 0.156* 0.666** 1
速效钾鄄 1 K鄄1 0.108 -.0039 0.042 0.031 0.451** 0.426** 1
速效钾鄄 2 K鄄2 0.115 -.0029 0.005 0.078 0.317** 0.458** 0.513** 1
摇 摇 * 表示在 0.05水平上有较显著的相关系数;** 表示在 0.01水平上有显著的相关系数; 表中 XXX鄄1表明是作物种植前土壤养分,XXX鄄2
表示是作物收获后的土壤养分
789摇 4期 摇 摇 摇 方斌摇 等:作物种植前后土壤有机质及养分因子的空间变异分析 摇
http: / / www.ecologica.cn
2.3摇 土壤因子空间分布特征分析
作物种植前、收获后的数据表明:(1)自相关系
数,种植前碱解氮、有机质、速效磷和速效钾 4 种因
子的自相关系数均< 25%,且均在 11.4% —15.8%之
间,表现出强烈的空间相关。 收获后,它们的自相关
系数均增大,且除碱解氮外,均处于中等空间相关程
度。 (2)变程,种植前,碱解氮和有机质的变程最小,
均为 720.00 m,最大的速效磷为 1830.00 m。 收获
后,各因子的变程普遍加大,最小的速效磷也达到了
6000.00 m。 (3)拟合模型,种植前各土壤因子均符
合指数同性模型,拟合效果较好。 收获后,拟合模型
有差异,根据 R2,有机质、碱解氮、速效磷拟合效果也
好于前者(表 4)。
表 4摇 作物种植前后土壤养分自相关特征
Table 4摇 The autocorrelation characteristic of soil nutrient after and before planting crops
土壤养分
Soil nutrient
C0 (C+C0)
C0 /
(C+C0)
A0 m 模型(Model) R2 RSS CV
有机质鄄 1 OM鄄1 5.10 41.46 0.12 240.00 720.00 Exponential 0.025 124.00000 22.63
有机质鄄 2 OM鄄2 0.08 0.33 0.26 38590.00 16839.84 Gaussian 0.119 0.09970 103.31
碱解氮鄄 1 N鄄1 0.00 0.03 0.12 240.00 720.00 Exponential 0.021 0.00009 18.43
碱解氮鄄 2 N鄄2 0.06 0.07 0.89 16365.53 16365.53 Gaussian 0.100 0.00217 33.24
速效磷鄄 1 P鄄1 1.23 7.81 0.16 610.00 1830.00 Exponential 0.524 3.43000 68.77
速效磷鄄 2 P鄄2 3.99 7.98 0.50 2000.00 6000.00 Exponential 0.765 2.38000 62.12
速效钾鄄 1 K鄄1 0.02 0.20 0.11 590.00 1770.00 Exponential 0.568 0.00171 45.59
速效钾鄄 2 K鄄2 0.23 0.62 0.36 83360.00 13360.00 Exponentia 0.194 0.12400 56.00
摇 摇 *C0 为块金值,C为偏基台值,C+C0 为基台值,C0 / (C+C0)为自相关系数,A0 为相关距离,m为变程,R2为决定系数,RSS为残差平方和,CV
为变异系数
摇 摇 通过对模型的精度进行分析,各土壤因子空间
分布基本符合最优模型的条件,即标准平均值最接
近 0,均方根预测误差最接近平均标准误差,标准均
方根预测误差最接近 1。 根据半方差分析的结果,利
用 ArcGIS 9. 0 中的 Geostatistic analysis 模块采用普
通 Kriging法进行内插,并根据土壤养分的含量对插
值图进行合理的分级,得到各养分含量等级面积的
空间变异(图 3—图 6)。
图 3摇 作物种植前后有机质的空间分布特征
Fig.3摇 Spacial Character of organic matter in the after and before planting crops
摇 摇 有机质含量范围在施肥前后基本保持一致,
90%以上处于 20—40 g / kg范围。 作物种植前,有机
质含量分布较为不均,含量较高的区域主要分布在
园地和林地区。 作物收获后其含量分布呈现总体有
差异、局部有均衡的现象,体现 2点特征:(1)作物收
获后,有机质含量处于 30—40 g / kg范围的区域面积
889 摇 生摇 态摇 学摇 报摇 摇 摇 34卷摇
http: / / www.ecologica.cn
减少,原有林地区的有机质下降较快,主要原因:其
一可能是种植区内作物在生长过程中吸收了土壤原
有的碳,而新的有机质的形成需要一定时间转化,因
此,总体表现为普遍性下降;其二可能是园地和林地
较近,林地的有机质含量通过园地的取样点有机质
含量体现,而水果在生长过程中消耗土壤中原有的
有机质,而短期内又无法及时得到补充。 (2)作物收
获后的有机质含量以水源地周边居高,并增加较快,
原因可能是有机质的形成与水有较密切的关系。 一
般而言同一地点土壤有机质在短时内不会有大的变
化,而两次试验结果却在不同地方有不同表现,这一
现象有必要针对不同地区产生的效应验证其碳收支
现象(图 3)。
作物种植前后碱解氮含量的空间分布格局变化
不大,作物种植前碱解氮含量 75%处于 50—60mg /
kg,作物收获后各区域的氮含量明显增加,范围在
70—90mg / kg占 65%。 主要特点:(1)林地和园地的
土壤碱解氮的含量较高,这可能与农户长期对经济
作物高施肥有较强的相关性。 (2)水源地,特别是水
库和水源下游地带氮的含量更高。 可能过多投入的
氮随着水分的流动而富集在水库、坑塘或流向地势
较低的区域。 (3)作物收获后,土壤碱解氮的含量普
遍高于作物种植前,主要原因在于土壤中的氮降解
和固化需要一定的时间。 这一点从第 2 年作物种植
前的土壤抽样结果得到证明。 具体见较图 4 作物收
获后速效磷的含量较种植前普遍性地有所提高,局
部性增加较大,其中城郊结合部和杭坪镇的磷含量
更高,超过了 100 mg / kg。 原因可能有 3点:(1)目前
农户主要使用的复合肥,导致速效磷含量普遍性有
所增加,但受使用量的影响,增加的幅度不大。 (2)
蔬菜和水果种植区速效磷的浓度增加大。 分析原因
可能是蔬菜的养分投入较大,除了氮投入超过一般
作物的 2—3倍外,磷的投入也有明显增加,这一点
可从农户养分投入调查中得到证实。 从地理位置上
看,速效磷含量最高的地方正是城郊蔬菜和高山蔬
菜所在地(图 5)。
图 4摇 作物种植前后碱解氮的空间分布特征
Fig.4摇 Spacial character of alkali鄄hydrolyzable N in the after and before planting crops
摇 摇 作物收获后速效钾的含量较种植前普遍性提高
不多,但局部性增加却有超过 150 mg / kg。 原因可能
有 3 点:(1)目前农户主要使用的复合肥,导致速效
钾含量普遍性有所增加,但钾含量并不高,因而,增
加的幅度不大。 (2)从全县整体上看,钾的投入没有
得到人们的重视,蔬菜区域钾浓度增大的原因与农
户蔬菜养分投入量大有直接关系。 此外,高山蔬菜
多是直接销售到上海,其养分及其主要成份会被检
测,县内有专职人员定期对养分的施用进行指导,钾
的投入也正逐步引起重视。 因此,对于农村大部分
地区,特别是经济作物种植较为密切的地区,配备养
分指导专家有针对性地进行土壤合理配方是必要的
(图 6)。
989摇 4期 摇 摇 摇 方斌摇 等:作物种植前后土壤有机质及养分因子的空间变异分析 摇
http: / / www.ecologica.cn
图 5摇 作物种植前后速效磷的空间分布特征
Fig.5摇 Spacial character of available P in the after and before planting crops
图 6摇 作物种植前后速效钾的空间分布特征
Fig.6摇 Spacial character of available K in the after and before planting crops
3摇 结论
(1)运用 Moran忆s I 指数分析、Pearson 相关系数
分析,以及 kriging 插值图等方法对作物种植前和收
获后的土壤因子空间变异进行分析,发现 3 种方法
所得出的结果可以相互应证,互为补充。 Moran忆s I
指数分析证实随着区域范围的加大,各因子间自相
关程度存在差异。 Pearson 相关系数分析可检验出
各因子间的相关性程度。 kriging 插值图则能直观反
应各因子在不同区域的分布特征。
(2)从养分利用特征看:1)有机质和碱解氮与水
的关系较为密切,水源区两因子的含量较周边高,可
能是水的作用加速了有机质向氮的转化。 2)养分投
入后,过多的养分在土壤中降解有一段时间过程,这
种过程与区域自然环境条件的关系较密切。 3)由于
主要是通过复合肥被施入到土壤的,测算结果反映
磷、钾之间有极强的相关性,作物收获后两者在土壤
中的含量均高于作物种植前,这可能与农户好用固
定磷、钾比的复合肥有较密切的关系,此外,磷、钾在
土壤中降解较慢也是原因之一。
099 摇 生摇 态摇 学摇 报摇 摇 摇 34卷摇
http: / / www.ecologica.cn
(3)同一作物的农户养分投入在不同区域差异
较大,不同作物的农户养分投入在区域内却存在量
的相关。 水稻等常规作物和规模化经营的葡萄的养
分区域内投入差异相对较小,而其它作物则表现较
大差异,空间分析结果与农户调查结果基本一致。
(4)高山蔬菜、城郊蔬菜种植区和壶源江下游地
带氮的投入过高应适当管控,山体周边土壤磷的含
量高也需要适当管控。 钾的投入普遍不高,虽然部
分区域钾含量呈较高的现象,但主要分布在高山区,
可能与该区域山体林木有关。
References:
[ 1 ]摇 Sheng J D, Xiao H, Wu H Q, Chen B, Wang J, Yang X J.
Primary study of spatial variability of available soil nutrients in
farmland with different sampling intervals. Agricultural Research
in the Arid Areas, 2005, 23(2): 63鄄67.
[ 2 ] 摇 Gao H, Zhang H C, Xing L, Wang Y, Zhu D J, Yan G Z. Study
on spatial variability of soil nutrients in different soil groups and
fertilization strategy. Journal of Yangzhou University: Agricultural
and Life Science Edition, 2011, 32(2): 42鄄47.
[ 3 ] 摇 Gao X Z, Hu K L, Guo Y, Li B G, Ma W D, Du S, Wang Y H.
Spatial variability of soil nutrients and crop yield and site鄄specific
fertilizer management. Scientia Agricultura Sinica, 2002, 35(6):
660鄄666.
[ 4 ] 摇 Liu Y, Sun Z M, Yang J, Liu X, Wang X M, Liu S Q, Wang W
J. Spatial variability of soil N, P, and K in main production area
of Castanea mollissima. Chinese Journal of Applied Ecology,
2010, 21(4): 901鄄907.
[ 5 ] 摇 Yuan X Y, Huang Y F, Gao R T, Cai X R, He Y. Spatial
variability characteristics of farmland soil organic matter in Pinggu
District, Beijing, China. Transactions of the Chinese Society of
Agricultural Engineering, 2008, 24(2): 70鄄76.
[ 6 ] 摇 Yang Q Y, Yang J S, Yao R J, Li D S. GIS鄄based soil suitability
evaluation of cultivated land in saline soil improvement district.
Journal of Natural Resources, 2011, 26(3): 1376鄄1384.
[ 7 ] 摇 Shi L J, Zheng L B, Liu Y L, Yu L Z, Jia Z C. Research on the
characteristics of spatial variability of soil nutrients and general
assessment in farmland in the Yangtze river delta. Resources and
Environment in the Yangtze Basin, 2008, 17(6): 839鄄846.
[ 8 ] 摇 Haefele S M, Wopereis M C S. Spatial variability of indigenous
supplies for N, P and K and its impact on fertilizer strategies for
irrigated rice in West Africa. Plant and Soil, 2005, 270 ( 1):
57鄄72.
[ 9 ] 摇 Sugihara S, Funakawa S, Kilasara M, Kosaki T. Dynamics of
microbial biomass nitrogen in relation to plant nitrogen uptake
during the crop growth period in a dry tropical cropland in
Tanzania. Soil Science and Plant Nutrition, 2010, 56 ( 1 ):
105鄄114.
[10] 摇 Zhao B, Sun L, Hu H Q, Sun Z H. Post鄄fire soil microbial
biomass and nutrient content of Larix gmelinii forest in autumn.
Journal of Natural Resources, 2011, 26(3): 450鄄459.
[11] 摇 Ma Q, Zhao G X. Effects of different land use types on soil
nutrients in intensive agricultural region. Journal of Natural
Resources, 2011, 25(11): 1834鄄1844.
[12] 摇 Gao Y M, Tong Y A, Hu Z Y, Wang C R, Zhang S L, Zhao H
B. Spatial variability of soil nutrients in farmland of Loess plateau.
Chinese Journal of Soil Science, 2006, 37(1): 1鄄6.
[13] 摇 Lian G, Guo X D, Fu B J, Hu C X. Spatial variability and
prediction of soil nutrients on a county scale on the Loess plateau鄄a
case study of Hengshan county, Shaanxi province. Acta Pedologica
Sinica, 2008, 45(4): 577鄄584.
[14] 摇 Hu W, Shao M A, Wang Q J. Scale鄄dependency of spatial
variability of soil moisture on a degraded slope鄄land on the Loess
Plateau. Transactions of the Chinese Society of Agricultural
Engineering, 2005, 21(8): 11鄄16.
[15] 摇 Gan H H, Peng L Y. Spatial variability of nutrients in cultivated
soils of Xinhui District, Jiangmen City. Chinese Journal of Applied
Ecology, 2005, 16(8): 1437鄄1442.
[16] 摇 Li W L, Zhao G X, Xu Z W. Cultivated land soil nutrient
precision management zoning and prescription fertilization at town
scale. Journal of Natural Resources, 2011, 26(5): 746鄄754.
[17] 摇 Wang J M, Wang J A, Hu Y M, Zhao Z Z. Studies on spatial
variability of soil nutrients and sampling methods in a field of
chestnut soil. Agricultural Research in the Arid Areas, 2006, 24
(5): 59鄄63.
[18] 摇 Herai Y, Kouno K, Hashimoto M, Nagaoka T. Relationships
between microbial biomass nitrogen, nitrate leaching and nitrogen
uptake by corn in a compost and chemical fertilizer鄄amended
regosol. Soil Science and Plant Nutrition, 2006, 52 ( 2 ):
186鄄194.
[19] 摇 Cambardella C A, Moorman T B, Parkin T B, Karlen D L,
Novak J M, Turco R F, Konopka A E. Field鄄scale variability of
soil properties in central Iowa soils. Soil Science Society of
America Journal, 1994, 58(5): 1501鄄1511.
[20] 摇 Che Q J, Cao Y H, Yu L, Song Z N, Dong Y W. Landscape
spatial heterogeneity is associated with urbanization: an example
from Yangtze River in Jiangsu Province. Acta Ecologica Sinica,
2011, 31(23): 7261鄄7270.
[21] 摇 Shi L J, Zheng L B, Liu Y L. Research on the characteristics of
spatial variability of soil nutrients in farmland. Journal of Henan
Agricultural University, 2008, 42(1): 51鄄56, 70鄄70.
[22] 摇 Dibari J N. Evaluation of five landscape鄄level metrics for
measuring the effects of urbanization on landscape structure: the
case of Tucson, Arizona, USA. Landscape and Urban Planning,
2007, 79(3 / 4): 308鄄310.
[23] 摇 Deng J S, Wang K, Hong Y, Qi J G. Spatio鄄temporal dynamics
199摇 4期 摇 摇 摇 方斌摇 等:作物种植前后土壤有机质及养分因子的空间变异分析 摇
http: / / www.ecologica.cn
and evolution of land use change and landscape pattern in response
to rapid urbanization. Landscape and Urban Planning, 2009, 92
(3 / 4): 187鄄198.
[24] 摇 Thomas I, Frankhauser P, Biernacki C. The morphology of built鄄
up landscapes in Wallonia ( Belgium): a classification using
fractal indices. Landscape and Urban Planning, 2008, 84(2):
99鄄115.
[25] 摇 Lu Y K. Soil Agricultural Chemical Analysis Method. Beijing:
China忆s Agricultural Science Press,2000.
[26] 摇 Fang B, Wu J f, Ni S X. Correlation analysis of spatial variability
of soil available nitrogen and household nitrogen inputs at Pujiang
County. Acta Ecologica Sinica, 2012, 32(20): 6489鄄6500.
参考文献:
[ 1 ]摇 盛建东, 肖华, 武红旗, 陈冰, 王军, 杨新建. 不同取样尺度
农田土壤速效养分空间变异特征初步研究. 干旱地区农业研
究, 2005, 23(2): 63鄄67.
[ 2 ] 摇 高辉, 张洪程, 邢琳, 王云, 朱德进, 严桂珠. 不同土类土壤
养分空间变异与施肥策略研究. 扬州大学学报: 农业与生命
科学版, 2011, 32(2): 42鄄47.
[ 3 ] 摇 高祥照, 胡克林, 郭焱, 李保国, 马韫韬, 杜森, 王运华. 土
壤养分与作物产量的空间变异特征与精确施肥. 中国农业科
学, 2002, 35(6): 660鄄666.
[ 4 ] 摇 刘杨, 孙志梅, 杨军, 刘霞, 王小敏, 刘树庆, 王文江. 京东
板栗主产区土壤氮磷钾的空间变异. 应用生态学报, 2010, 21
(4): 901鄄907.
[ 5 ] 摇 苑小勇, 黄元仿, 高如泰, 柴旭荣, 贺勇. 北京市平谷区农用
地土壤有机质空间变异特征. 农业工程学报, 2008, 24(2):
70鄄76.
[ 6 ] 摇 杨奇勇, 杨劲松, 李晓明, 李冬顺. 基于 GIS的盐渍土改良区
耕地土壤适宜性评价. 自然资源学报, 2011, 26 ( 3 ):
1376鄄1384.
[ 7 ] 摇 史利江, 郑丽波, 柳云龙, 俞立中, 贾正长. 长三角地区农田
土壤养分空间变异及养分综合评价. 长江流域资源与环境,
2008, 17(6): 839鄄846.
[10] 摇 赵彬, 孙龙, 胡海清, 孙志虎. 兴安落叶松林火后对土壤养分
和土壤微生物生物量的影响. 自然资源学报, 2011, 26(3):
450鄄459.
[11] 摇 马群, 赵庚星. 集约农区不同土地利用方式对土壤养分状况
的影响. 自然资源学报, 2011, 25(11): 1834鄄1844.
[12] 摇 高义民, 同延安, 胡正义, 王彩绒, 张树兰, 赵护兵. 黄土区
村级农田土壤养分空间变异特征研究. 土壤通报, 2006, 37
(1): 1鄄6.
[13] 摇 连纲, 郭旭东, 傅伯杰, 虎陈霞. 黄土高原县域土壤养分空间
变异特征及预测———以陕西省横山县为例. 土壤学报, 2008,
45(4): 577鄄584.
[14] 摇 胡伟,邵明安,王全九. 黄土高原退耕坡地土壤水分空间变异
的尺度性研究. 农业工程学报, 2005, 21(8): 11鄄16.
[15] 摇 甘海华, 彭凌云. 江门市新会区耕地土壤养分空间变异特征.
应用生态学报, 2005, 16(8): 1437鄄1442.
[16] 摇 李文璐,赵庚星, 许志伟. 乡镇级耕地土壤养分精确化管理分
区及配方施肥研究. 自然资源学报, 2011, 26(5): 746鄄754.
[17] 摇 王晋民, 王俊鹏, 胡月明, 赵之重. 栗钙土农田土壤养分空间
变异特性及采样方法研究. 干旱地区农业研究, 2006, 24
(5): 59鄄63.
[20] 摇 车前进, 曹有挥, 于露, 宋正娜, 董雅文. 景观空间异质性及
城市化关联———以江苏省沿江地区为例. 生态学报, 2011, 31
(23): 7261鄄7270.
[21] 摇 史利江, 郑丽波, 柳云龙. 农田土壤养分空间变异特征研究.
河南农业大学学报, 2008, 42(1): 51鄄56, 70鄄70.
[25] 摇 鲁如坤. 土壤农业化学分析方法. 北京: 中国农业科技出版
社, 2000.
[26] 摇 方斌,吴金凤,倪绍祥. 浦江县土壤碱解氮的空间变异与农户
N投入的关联分析. 生态学报, 2012, 32(20): 6489鄄6500.
299 摇 生摇 态摇 学摇 报摇 摇 摇 34卷摇
ＡＣＴＡ ＥＣＯＬＯＧＩＣＡ ＳＩＮＩＣＡ Ｖｏｌ．３４，Ｎｏ．４ Ｆｅｂ．，２０１４（Ｓｅｍｉｍｏｎｔｈｌｙ）
ＣＯＮＴＥＮＴＳ
Ｆｒｏｎｔｉｅｒｓ ａｎｄ Ｃｏｍｐｒｅｈｅｎｓｉｖｅ Ｒｅｖｉｅｗ
Ｔｈｅ ｂｉｏａｖａｉｌａｂｉｌｉｔｙ ｏｆ ｄｉｓｓｏｌｖｅｄ ｏｒｇａｎｉｃ ｃａｒｂｏｎ ｉｎ ｔｈｅ ｅｕｔｒｏｐｈｉｃ ｌａｋｅｓ ＹＥ Ｌｉｎｌｉｎ， ＫＯＮＧ Ｆａｎｘｉａｎｇ， ＳＨＩ Ｘｉａｏｌｉ，ｅｔ ａｌ （７７９）…………
Ｐｌａｎｔ ｓｐｅｃｉｅｓ ｏｆ ｔｈｅ ｎｏｎ⁃ａｇｒｉｃｕｌｔｕｒａｌ ｈａｂｉｔａｔｓ ｉｎ ｔｈｅ ｌｏｗｅｒ ｒｅａｃｈｅｓ ｏｆ ｔｈｅ Ｙｅｌｌｏｗ Ｒｉｖｅｒ ｐｌａｉｎ ａｇｒｏ⁃ｌａｎｄｓｃａｐｅ
ＬＵ Ｘｕｎｌｉｎｇ， ＬＩＡＮＧ Ｇｕｏｆｕ， ＴＡＮＧ Ｑｉａｎ， ｅｔ ａｌ （７８９）
…………………………
………………………………………………………………………………
Ａｕｔｅｃｏｌｏｇｙ ＆ Ｆｕｎｄａｍｅｎｔａｌｓ
Ｍａｎｇａｎｅｓｅ ｓｔｒｅｓｓ ｏｎ ｔｈｅ ｕｌｔｒａｓｔｒｕｃｔｕｒｅｓ ｏｆ ａ ｍａｎｇａｎｅｓｅ ｔｏｌｅｒａｎｔ ｐｌａｎｔ，Ｐｏｌｙｇｏｎｕｍ ｐｅｒｆｏｌｉａｔｕｍ Ｌ．
ＷＡＮＧ Ｊｕｎ， ＷＵ Ｈｕｉ， ＸＵＥ Ｓｈｅｎｇｇｕｏ， ｅｔ ａｌ （７９８）
………………………………………
…………………………………………………………………………………
Ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ａｒｓｅｎｉｃ （Ａｓ） ｔｏｌｅｒａｎｃｅ ａｎｄ ａｃｃｕｍｕｌａｔｉｏｎ ｉｎ ｒｉｃｅ （Ｏｒｙｚａ ｓａｔｉｖａ Ｌ．） ｇｅｎｏｔｙｐｅｓ ｗｉｔｈ ｄｉｆｆｅｒｅｎｔ ｒａｄｉａｌ ｏｘｙｇｅｎ ｌｏｓｓ
ＷＵ Ｃｈｕａｎ， ＭＯ Ｊｉｎｇｙｕ， ＸＵＥ Ｓｈｅｎｇｇｕｏ， ｅｔ ａｌ （８０７）
……
………………………………………………………………………………
Ｅｆｆｅｃｔｓ ｏｆ ｗａｔｅｒ ｓｔｒｅｓｓ ｏｎ ｐｈｙｓｉｏｌｏｇｉｃａｌ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ｄｉｆｆｅｒｅｎｔ Ｉｌｌｉｃｉｕｍ ｌａｎｃｅｏｌａｔｕｍ ｅｃｏｔｙｐｅｓ ｕｎｄｅｒ ｌｏｗ ｌｉｇｈｔ ｉｎｔｅｎｓｉｔｙ
ＣＡＯ Ｙｏｎｇｈｕｉ， ＺＨＯＵ Ｂｅｎｚｈｉ， ＣＨＥＮ Ｓｈｕａｎｇｌｉｎ （８１４）
……………
……………………………………………………………………………
Ｅｆｆｅｃｔ ｏｆ ｂｒａｎｃｈ ｎｕｍｂｅｒ ｏｎ ｔｈｅ ｇｒｏｗｔｈ ａｎｄ ｄｅｖｅｌｏｐｍｅｎｔ ｏｆ Ｍｏｒｕｓ ａｌｂａ ｓａｐｌｉｎｇｓ ＨＵＡＮ Ｈｕｉｈｕｉ， ＸＵ Ｘｉａｏ， ＬＩＵ Ｇａｎｇ， ｅｔ ａｌ （８２３）…
Ｓｐａｔｉａｌ ｄｉｓｔｒｉｂｕｔｉｏｎ ｐａｔｔｅｒｎ ａｎｄ ｓａｍｐｌｉｎｇ ｔｅｃｈｎｉｑｕｅ ｆｏｒ Ｏｒｔｈｏｔｙｌｕｓ （Ｏ．） ｓｏｐｈｏｒａｅ ｎｙｍｐｈｓ ｏｎ Ｓｏｐｈｏｒａ ｊａｐｏｎｉｃａ
ＺＨＵ Ｈｕｉｙｉｎｇ， ＳＨＥＮ Ｐｉｎｇ，ＷＵ Ｊｉａｎｈｕａ， ｅｔ ａｌ （８３２）
…………………………
………………………………………………………………………………
Ａｓｓｅｓｓｍｅｎｔ ｏｆ ｆｕｎｇａｌ ｄｉｖｅｒｓｉｔｙ ｉｎ ａｐｐｌｅ ｒｅｐｌａｎｔｅｄ ｏｒｃｈａｒｄ ｓｏｉｌｓ ｂｙ Ｔ⁃ＲＦＬＰ ａｎａｌｙｓｉｓ
ＹＩＮ Ｃｈｅｎｇｍｉａｏ， ＷＡＮＧ Ｇｏｎｇｓｈｕａｉ， ＬＩ Ｙｕａｎｙｕａｎ， ｅｔ ａｌ （８３７）
……………………………………………………
……………………………………………………………………
Ｅｆｆｅｃｔｓ ｏｆ ｄａｚｏｍｅｔ ｏｎ ｅｄａｐｈｏｎ ａｎｄ ｇｒｏｗｔｈ ｏｆ Ｍａｌｕｓ ｈｕｐｅｈｅｎｓｉｓ ｒｅｈｄ． ｕｎｄｅｒ ｃｏｎｔｉｎｕｏｕｓ ａｐｐｌｅ ｃｒｏｐｐｉｎｇ
ＬＩＵ Ｅｎｔａｉ， ＬＩ Ｙｕａｎｙｕａｎ， ＨＵ Ｙａｎｌｉ， ｅｔ ａｌ （８４７）
………………………………
…………………………………………………………………………………
Ｉｓｏｌａｔｉｏｎ， ｉｄｅｎｔｉｆｉｃａｔｉｏｎ， ａｎｄ ｐｅｒｆｏｒｍａｎｃｅ ｏｆ ｔｗｏ ｐｙｒｅｎｅ⁃ｄｅｇｒａｄｉｎｇ ｅｎｄｏｐｈｙｔｉｃ ｂａｃｔｅｒｉａ ＳＵＮ Ｋａｉ， ＬＩＵ Ｊｕａｎ， ＬＩ Ｘｉｎ， ｅｔ ａｌ （８５３）……
Ｐｏｐｕｌａｔｉｏｎ， Ｃｏｍｍｕｎｉｔｙ ａｎｄ Ｅｃｏｓｙｓｔｅｍ
Ｅｆｆｅｃｔｓ ｏｆ ｄｉｆｆｅｒｅｎｔ ｔｅｍｐｅｒａｔｕｒｅｓ ｏｎ ｔｈｅ ｇｒｏｗｔｈ ａｎｄ ｄｅｖｅｌｏｐｍｅｎｔ ｏｆ Ｅｏｔｅｔｒａｎｙｃｈｕｓ Ｋａｎｋｉｔｕｓ （Ｅｈａｒａ）
ＬＩ Ｙｉｎｇｊｉｅ，ＷＡＮＧ Ｚｉｙｉｎｇ，ＺＨＡＮＧ Ｇｕｏｈａｏ，ｅｔ ａｌ （８６２）
…………………………………
……………………………………………………………………………
Ｅｆｆｅｃｔ ｏｆ ａｖａｉｌａｂｌｅ ｂｕｒｒｏｗ ｄｅｎｓｉｔｉｅｓ ｏｆ ｐｌａｔｅａｕ ｐｉｋａ （Ｏｃｈｏｔｏｎａ ｃｕｒｚｏｎｉａｅ） ｏｎ ｐｌａｎｔ ｎｉｃｈｅ ｏｆ ａｌｐｉｎｅ ｍｅａｄｏｗ ｃｏｍｍｕｎｉｔｉｅｓ ｉｎ ｔｈｅ
Ｑｉｎｇｈａｉ⁃Ｔｉｂｅｔ Ｐｌａｔｅａｕ ＪＩＡ Ｔｉｎｇｔｉｎｇ， ＭＡＯ Ｌｉａｎｇ， ＧＵＯ Ｚｈｅｎｇｇａｎｇ （８６９）…………………………………………………………
Ｃｏｒｒｅｌａｔｉｏｎ ｂｅｔｗｅｅｎ ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ Ｒｅａｕｍｕｒｉａ ｓｏｏｎｇａｒｉｃａ ｃｏｍｍｕｎｉｔｉｅｓ ａｎｄ ｓｏｉｌ ｆａｃｔｏｒｓ ｉｎ ｔｈｅ Ｓａｎｇｏｎｇ Ｒｉｖｅｒ ｂａｓｉｎ
ＺＨＡＯ Ｘｕｅｃｈｕｎ， ＬＡＩ Ｌｉｍｉｎｇ， ＺＨＵ Ｌｉｎｈａｉ， ｅｔ ａｌ （８７８）
…………………
……………………………………………………………………………
Ｅｆｆｅｃｔｓ ｏｆ ａｆｆｏｒｅｓｔａｔｉｏｎ ｏｎ ｓｏｉｌ ｍｉｃｒｏｂｉａｌ ｃｏｍｍｕｎｉｔｙ ｓｔｒｕｃｔｕｒｅ ｉｎ ｔｈｅ ａｒｉｄ ｖａｌｌｅｙ ｏｆ Ｍｉｎｊｉａｎｇ Ｒｉｖｅｒ
ＷＡＮＧ Ｗｅｉｘｉａ， ＬＵＯ Ｄａ， ＳＨＩ Ｚｕｏｍｉｎ， ｅｔ ａｌ （８９０）
……………………………………
………………………………………………………………………………
Ｅｆｆｅｃｔｓ ｏｆ ｒｅｃｌａｍａｔｉｏｎ ｏｎ ｔｉｄａｌ ｆｌａｔ ａｎｄ ｌａｎｄ ｕｓｅ ｏｎ ｓｏｉｌ ｍｉｃｒｏｂｉａｌ ｃｏｍｍｕｎｉｔｙ
ＬＩＮ Ｌｉ， ＣＵＩ Ｊｕｎ， ＣＨＥＮ Ｘｕｅｐｉｎｇ， ＦＡＮＧ Ｃｈａｎｇｍｉｎｇ （８９９）
…………………………………………………………
………………………………………………………………………
Ｅｆｆｅｃｔｓ ｏｆ Ｐｏｍａｃｅａ ｃａｎａｌｉｃｕｌａｔａ ｏｎ ａｑｕａｔｉｃ ｍａｃｒｏｐｈｙｔｅ ｃｏｍｍｕｎｉｔｙ ｓｔｒｕｃｔｕｒｅ ｉｎ ｐａｄｄｙ ｆｉｅｌｄｓ
ＺＨＡＯ Ｂｅｎｌｉａｎｇ， ＺＨＡＮＧ Ｊｉａｅｎ， ＤＡＩ Ｘｉａｏｙａｎ， ｅｔ ａｌ （９０７）
…………………………………………
………………………………………………………………………
Ｔｈｅ ａｄａｐｔａｂｉｌｉｔｙ ａｎｄ ｄｅｃｏｎｔａｍｉｎａｔｉｏｎ ｅｆｆｅｃｔ ｏｆ ｆｏｕｒ ｋｉｎｄｓ ｏｆ ｗｏｏｄｙ ｐｌａｎｔｓ ｉｎ ｃｏｎｓｔｒｕｃｔｅｄ ｗｅｔｌａｎｄ ｅｎｖｉｒｏｎｍｅｎｔ
ＣＨＥＮ Ｙｏｎｇｈｕａ， ＷＵ Ｘｉａｏｆｕ， ＨＡＯ Ｊｕｎ， ｅｔ ａｌ （９１６）
………………………
………………………………………………………………………………
Ｃａｒｂｏｎ ｂｕｄｇｅｔ ｏｆ ａｌｐｉｎｅ Ｐｏｔｅｎｔｉｌｌａ ｆｒｕｔｉｃｏｓａ ｓｈｒｕｂｌａｎｄ ｂａｓｅｄ ｏｎ ｃｏｍｐｒｅｈｅｎｓｉｖｅ ｔｅｃｈｎｉｑｕｅｓ ｏｆ ｓｔａｔｉｃ ｃｈａｍｂｅｒ ａｎｄ ｂｉｏｍａｓｓ ｈａｒｖｅｓｔｉｎｇ
ＬＩ Ｈｏｎｇｑｉｎ， ＬＩ Ｙｉｎｇｎｉａｎ， ＺＨＡＮＧ Ｆａｗｅｉ， ｅｔ ａｌ （９２５）
…
……………………………………………………………………………
Ｅｆｆｅｃｔ ｏｆ ｉｎｉｔｉａｌ ｐＨ ｖａｌｕｅ ｏｎ ｍｉｃｒｏｂｉａｌ Ｆｅ （Ⅲ） ｒｅｄｕｃｔｉｏｎ ｉｎ ａｌｋａｌｉｎｅ ａｎｄ ａｃｉｄｉｃ ｐａｄｄｙ ｓｏｉｌｓ ＷＵ Ｃｈａｏ，ＱＵ Ｄｏｎｇ，ＬＩＵ Ｈａｏ （９３３）…
Ｌａｎｄｓｃａｐｅ， Ｒｅｇｉｏｎａｌ ａｎｄ Ｇｌｏｂａｌ Ｅｃｏｌｏｇｙ
Ｃｌｉｍａｔｅ ｅｎｖｉｒｏｎｍｅｎｔａｌ ｃｈａｎｇｅ ａｎｄ ｓｔａｂｌｅ ｃａｒｂｏｎ ｉｓｏｔｏｐｅｓ ｉｎ ａｇｅ ｌａｙｅｒｓ ｏｆ Ｔａｍａｒｉｘ ｓａｎｄ⁃ｈｉｌｌｏｃｋｓ ｉｎ Ｋｕｍｔａｇ ｄｅｓｅｒｔ
ＺＨＡＮＧ Ｊｉｎｃｈｕｎ，ＹＡＯ Ｔｕｏ， ＬＩＵ Ｃｈａｎｇｚｈｏｎｇ， ｅｔ ａｌ （９４３）
……………………
…………………………………………………………………………
Ｒｅｓｏｕｒｃｅ ａｎｄ Ｉｎｄｕｓｔｒｉａｌ Ｅｃｏｌｏｇｙ
Ｔｈｅ ｃｒｉｔｉｃａｌ ｆａｃｔｏｒｓ ｏｆ ｔｒａｎｓｐｉｒａｔｉｏｎ ｏｎ ｍｕｓｋｍｅｌｏｎ ｉｎ ｐｌａｓｔｉｃ ｇｒｅｅｎｈｏｕｓｅ ＺＨＡＮＧ Ｄａｌｏｎｇ，ＣＨＡＮＧ Ｙｉｂｏ， ＬＩ Ｊｉａｎｍｉｎｇ，ｅｔ ａｌ （９５３）……
Ｉｏｎ ａｂｓｏｒｐｔｉｏｎ ａｎｄ ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｆ ｓｙｍｂｉｏｔｉｃ Ｒｅａｕｍｕｒｉａ ｓｏｏｎｇｏｒｉｃａ ａｎｄ Ｓａｌｓｏｌａ ｐａｓｓｅｒｉｎａ ｓｅｅｄｌｉｎｇｓ ｕｎｄｅｒ ＮａＣｌ ｓｔｒｅｓｓ
ＺＨＡＯ Ｘｉｎ， ＹＡＮＧ Ｘｉａｏｊｕ， ＳＨＩ Ｙｏｎｇ， ｅｔ ａｌ （９６３）
…………………
…………………………………………………………………………………
Ｔｈｅ ｒｅｌａｔｉｏｎｓｈｉｐ ｂｅｔｗｅｅｎ ｓｅｌｅｃｔｅｄ ｒｈｉｚｏｓｐｈｅｒｅ ａｎｄ ｎｏｎ⁃ｒｈｉｚｏｓｐｈｅｒｅ ｓｏｉｌ ｐｒｏｐｅｒｔｉｅｓ ａｎｄ ｔｈｅ ｑｕａｌｉｔｙ ｏｆ Ｐｙｒｏｌａ ｄｅｃｏｒａｔａ
ＧＥＮＧ Ｚｅｎｇｃｈａｏ，ＭＥＮＧ Ｌｉｎｇｊｕｎ，ＬＩＵ Ｊｉａｎｊｕｎ （９７３）
…………………
………………………………………………………………………………
Ｓｐａｔｉａｌ ｖａｒｉａｔｉｏｎ ａｎａｌｙｓｉｓ ｏｆ ｓｏｉｌ ｏｒｇａｎｉｃ ｍａｔｔｅｒ ａｎｄ ｎｕｔｒｉｅｎｔ ｆａｃｔｏｒ ｆｏｒ ｂｅｆｏｒｅ ａｎｄ ａｆｔｅｒ ｐｌａｎｔｉｎｇ ｃｒｏｐｓ
ＦＡＮＧ Ｂｉｎ， ＷＵ Ｊｉｎｆｅｎｇ （９８３）
…………………………………
………………………………………………………………………………………………………
Ｕｒｂａｎ， Ｒｕｒａｌ ａｎｄ Ｓｏｃｉａｌ Ｅｃｏｌｏｇｙ
Ｅｓｔａｂｌｉｓｈｍｅｎｔ ａｎｄ ａｐｐｌｉｃａｔｉｏｎ ｏｆ ｔｈｅ ｉｎｄｅｘ ｓｙｓｔｅｍ ｆｏｒ ｕｒｂａｎ ｒｉｖｅｒ ｈｅａｌｔｈ ａｓｓｅｓｓｍｅｎｔ
ＤＥＮＧ Ｘｉａｏｊｕｎ，ＸＵ Ｙｏｕｐｅｎｇ，ＺＨＡＩ Ｌｕｘｉｎ，ｅｔ ａｌ （９９３）
……………………………………………………
………………………………………………………………………………
Ｄｙｎａｍｉｃ ａｎａｌｙｓｉｓ ｏｆ ｔｈｅ ｅｃｏｌｏｇｉｃａｌ ｆｏｏｔｐｒｉｎｔ ａｎｄ ｃａｒｒｙｉｎｇ ｃａｐａｃｉｔｙ ｏｆ ｔｉｂｅｔ ＡＮ Ｂａｏｓｈｅｎｇ， ＣＨＥＮＧ Ｇｕｏｄｏｎｇ （１００２）…………………
Ｒｅｓｅａｒｃｈ Ｎｏｔｅｓ
Ｒｅｓｐｏｎｓｅｓ ｏｆ ｓｏｉｌ ｍｉｃｒｏｏｒｇａｎｉｓｍｓ ａｎｄ ｓｏｉｌ ｅｎｚｙｍｅ ａｃｔｉｖｉｔｉｅｓ ｔｏ ｄｉｆｆｅｒｅｎｔ ｌａｎｄ ｕｓｅ ｐａｔｔｅｒｎｓ ｉｎ ｔｈｅ ｗａｔｅｒ⁃ｌｅｖｅｌ⁃ｆｌｕｃｔｕａｔｉｎｇ ｚｏｎｅ ｏｆ
ｔｈｅ Ｔｈｒｅｅ Ｇｏｒｇｅｓ Ｒｅｓｅｒｖｏｉｒ ｒｅｇｉｏｎ ＭＡ Ｐｅｎｇ， ＬＩ Ｃｈａｎｇｘｉａｏ， ＬＥＩ Ｍｉｎｇ， ｅｔ ａｌ （１０１０）……………………………………………
Ｅｆｆｅｃｔｓ ｏｆ ｓａｌｔ ｓｔｒｅｓｓ ｏｎ ｇｒｏｗｔｈ ａｎｄ ｒｏｏｔ ｄｅｖｅｌｏｐｍｅｎｔ ｏｆ ｔｗｏ ｏａｋ ｓｅｅｄｌｉｎｇｓ
ＷＡＮＧ Ｓｈｕｆｅｎｇ， ＨＵ Ｙｕｎｘｕｅ， ＳＵＮ Ｈａｉｊｉｎｇ， ｅｔ ａｌ （１０２１）
………………………………………………………………
…………………………………………………………………………
Ｔｈｅ ｅｆｆｅｃｔｓ ｏｆ ｃｏｎｓｔａｎｔ ａｎｄ ｖａｒｉａｂｌｅ ｔｈｅｒｍａｌ ａｃｃｌｉｍａｔｉｏｎ ｏｎ ｔｈｅｒｍａｌ ｔｏｌｅｒａｎｃｅ ｏｆ ｔｈｅ ｃｏｍｍｏｎ ｇｉａｎｔ ｔｏａｄ ｔａｄｐｏｌｅｓ （Ｂｕｆｏ ｇａｒｇａｒｉｚａｎｓ）
ＷＡＮＧ Ｌｉｚｈｉ （１０３０）
…
…………………………………………………………………………………………………………………
２４０１ 　 生　 态　 学　 报　 　 　 ３４卷　
叶生态学报曳圆园员源年征订启事
叶生态学报曳是由中国科学技术协会主管袁中国生态学学会尧中国科学院生态环境研究中心主办的生态学
高级专业学术期刊袁创刊于 员怨愿员年袁报道生态学领域前沿理论和原始创新性研究成果遥 坚持野百花齐放袁百家
争鸣冶的方针袁依靠和团结广大生态学科研工作者袁探索生态学奥秘袁为生态学基础理论研究搭建交流平台袁
促进生态学研究深入发展袁为我国培养和造就生态学科研人才和知识创新服务尧为国民经济建设和发展服务遥
叶生态学报曳主要报道生态学及各分支学科的重要基础理论和应用研究的原始创新性科研成果遥 特别欢
迎能反映现代生态学发展方向的优秀综述性文章曰研究简报曰生态学新理论尧新方法尧新技术介绍曰新书评价和
学术尧科研动态及开放实验室介绍等遥
叶生态学报曳为半月刊袁大 员远开本袁圆愿园页袁国内定价 怨园元 辕册袁全年定价 圆员远园元遥
国内邮发代号院愿圆鄄苑袁国外邮发代号院酝远苑园
标准刊号院陨杂杂晕 员园园园鄄园怨猿猿摇 摇 悦晕 员员鄄圆园猿员 辕 匝
全国各地邮局均可订阅袁也可直接与编辑部联系购买遥 欢迎广大科技工作者尧科研单位尧高等院校尧图书
馆等订阅遥
通讯地址院 员园园园愿缘 北京海淀区双清路 员愿号摇 电摇 摇 话院 渊园员园冤远圆怨源员园怨怨曰 远圆愿源猿猿远圆
耘鄄皂葬蚤造院 泽澡藻灶早贼葬蚤曾怎藻遭葬燥岳 则糟藻藻泽援葬糟援糟灶摇 网摇 摇 址院 憎憎憎援藻糟燥造燥早蚤糟葬援糟灶
本期责任副主编摇 董摇 鸣摇 摇 摇 编辑部主任摇 孔红梅摇 摇 摇 执行编辑摇 刘天星摇 段摇 靖
生摇 态摇 学摇 报
渊杂匀耘晕郧栽粤陨摇 载哉耘月粤韵冤
渊半月刊摇 员怨愿员年 猿月创刊冤
第 猿源卷摇 第 源期摇 渊圆园员源年 圆月冤
粤悦栽粤 耘悦韵蕴韵郧陨悦粤 杂陨晕陨悦粤
摇
渊杂藻皂蚤皂燥灶贼澡造赠袁杂贼葬则贼藻凿 蚤灶 员怨愿员冤
摇
灾燥造郾 猿源摇 晕燥郾 源 渊云藻遭则怎葬则赠袁 圆园员源冤
编摇 摇 辑摇 叶生态学报曳编辑部
地址院北京海淀区双清路 员愿号
邮政编码院员园园园愿缘
电话院渊园员园冤远圆怨源员园怨怨
憎憎憎援藻糟燥造燥早蚤糟葬援糟灶
泽澡藻灶早贼葬蚤曾怎藻遭葬燥岳 则糟藻藻泽援葬糟援糟灶
主摇 摇 编摇 王如松
主摇 摇 管摇 中国科学技术协会
主摇 摇 办摇 中国生态学学会
中国科学院生态环境研究中心
地址院北京海淀区双清路 员愿号
邮政编码院员园园园愿缘
出摇 摇 版摇
摇 摇 摇 摇 摇 地址院北京东黄城根北街 员远号
邮政编码院员园园苑员苑
印摇 摇 刷摇 北京北林印刷厂
发 行摇
地址院东黄城根北街 员远号
邮政编码院员园园苑员苑
电话院渊园员园冤远源园猿源缘远猿
耘鄄皂葬蚤造院躁燥怎则灶葬造岳 糟泽责早援灶藻贼
订摇 摇 购摇 全国各地邮局
国外发行摇 中国国际图书贸易总公司
地址院北京 猿怨怨信箱
邮政编码院员园园园源源
广告经营
许 可 证摇 京海工商广字第 愿园员猿号
耘凿蚤贼藻凿 遭赠摇 耘凿蚤贼燥则蚤葬造 遭燥葬则凿 燥枣
粤悦栽粤 耘悦韵蕴韵郧陨悦粤 杂陨晕陨悦粤
粤凿凿院员愿袁杂澡怎葬灶早择蚤灶早 杂贼则藻藻贼袁匀葬蚤凿蚤葬灶袁月藻蚤躁蚤灶早 员园园园愿缘袁悦澡蚤灶葬
栽藻造院渊园员园冤远圆怨源员园怨怨
憎憎憎援藻糟燥造燥早蚤糟葬援糟灶
泽澡藻灶早贼葬蚤曾怎藻遭葬燥岳 则糟藻藻泽援葬糟援糟灶
耘凿蚤贼燥则鄄蚤灶鄄糟澡蚤藻枣摇 宰粤晕郧 砸怎泽燥灶早
杂怎责藻则增蚤泽藻凿 遭赠摇 悦澡蚤灶葬 粤泽泽燥糟蚤葬贼蚤燥灶 枣燥则 杂糟蚤藻灶糟藻 葬灶凿 栽藻糟澡灶燥造燥早赠
杂责燥灶泽燥则藻凿 遭赠摇 耘糟燥造燥早蚤糟葬造 杂燥糟蚤藻贼赠 燥枣 悦澡蚤灶葬
砸藻泽藻葬则糟澡 悦藻灶贼藻则 枣燥则 耘糟燥鄄藻灶增蚤则燥灶皂藻灶贼葬造 杂糟蚤藻灶糟藻泽袁 悦粤杂
粤凿凿院员愿袁杂澡怎葬灶早择蚤灶早 杂贼则藻藻贼袁匀葬蚤凿蚤葬灶袁月藻蚤躁蚤灶早 员园园园愿缘袁悦澡蚤灶葬
孕怎遭造蚤泽澡藻凿 遭赠摇 杂糟蚤藻灶糟藻 孕则藻泽泽
粤凿凿院员远 阅燥灶早澡怎葬灶早糟澡藻灶早早藻灶 晕燥则贼澡 杂贼则藻藻贼袁
月藻蚤躁蚤灶早摇 员园园苑员苑袁悦澡蚤灶葬
孕则蚤灶贼藻凿 遭赠摇 月藻蚤躁蚤灶早 月藻蚤 蕴蚤灶 孕则蚤灶贼蚤灶早 匀燥怎泽藻袁
月藻蚤躁蚤灶早 员园园园愿猿袁悦澡蚤灶葬
阅蚤泽贼则蚤遭怎贼藻凿 遭赠摇 杂糟蚤藻灶糟藻 孕则藻泽泽
粤凿凿院员远 阅燥灶早澡怎葬灶早糟澡藻灶早早藻灶 晕燥则贼澡
杂贼则藻藻贼袁月藻蚤躁蚤灶早 员园园苑员苑袁悦澡蚤灶葬
栽藻造院渊园员园冤远源园猿源缘远猿
耘鄄皂葬蚤造院躁燥怎则灶葬造岳 糟泽责早援灶藻贼
阅燥皂藻泽贼蚤糟 摇 摇 粤造造 蕴燥糟葬造 孕燥泽贼 韵枣枣蚤糟藻泽 蚤灶 悦澡蚤灶葬
云燥则藻蚤早灶 摇 摇 悦澡蚤灶葬 陨灶贼藻则灶葬贼蚤燥灶葬造 月燥燥噪 栽则葬凿蚤灶早
悦燥则责燥则葬贼蚤燥灶
粤凿凿院孕援韵援月燥曾 猿怨怨 月藻蚤躁蚤灶早 员园园园源源袁悦澡蚤灶葬
摇 陨杂杂晕 员园园园鄄园怨猿猿
悦晕 员员鄄圆园猿员 辕 匝
国内外公开发行 国内邮发代号 愿圆鄄苑 国外发行代号 酝远苑园 定价 怨园郾 园园元摇