全 文 ：福建林学院学报 2009, 29(2):97-102
JournalofFujianCollegeofForestry
第 29卷 第 2期
2009年 4月
收稿日期:2008-12-24 修回日期:2009-03-02
基金项目:教育部博士点学科专项基金资助项目(200803890010)。
作者简介:侯丹(1982-),女 ,辽宁沈阳人 ,硕士研究生 ,从事森林生态学研究。通讯作者吴承祯(1970-),男 ,江西吉安人 ,教授 ,博士生导
师 ,从事森林生态学研究。 E-mail:fjwcz@ 126.com。
邓恩桉胸径生长对模拟硫 、氮复合沉降的响应
侯 丹 , 吕广林 , 张艳荷 , 包和林 , 洪 伟 , 吴承祯
(福建农林大学桉树研究中心 ,福建 福州 350002)
摘要:应用二次回归旋转设计方法 ,安排硫 、氮复合沉降对邓恩桉幼林影响的模拟实验 , 研究邓恩桉幼龄林胸径增长量对模
拟硫 、氮复合沉降的响应。排除年龄生长的影响 , 建立邓恩桉胸径生长响应的回归方程:y =0.776 3 +0.193 5 x1 -
0.125 6 x2-0.510 0x1x2 -0.245 0 x21 -0.597 5 x22 ,分析硫 、氮复合沉降对胸径增长的影响。 结果表明 , 硫编码值为 0.394 9
(纯硫为28.20 kg· hm-2· a-1)和氮编码值为 -0.105 1(纯氮为 64.78 kg· hm-2· a-1)时 , 邓恩桉胸径年增长量可望达到
最大 0.84 cm。
关键词:邓恩桉人工林;硫 、氮复合沉降;二次回归旋转设计;胸径生长
中图分类号:S718.51+2 文献标识码:A 文章编号:1001-389X(2009)02-0125-06
EfectofsimulatedNandScombineddepositiononDBHofEucalyptusdunni
HOUDan, L Guang-lin, ZHANGYan-he, BAOHe-lin, HONGWei, WUCheng-zhen
(ResearchCenterforEucalyptus, FujianAgricultureandForestryUniversity, Fuzhou, Fujian350002, China)
Abstract:Theexperimentabouttheefectofsimilatedsulphur(S)andnitrogen(N)combineddepositiononthegrowthof
diameteratbreastheight(DBH)ofyoungEucalyptusdunniiplantationwasarrangedbythemethodofquadraticregressionrevolution
design.TheresponseofDBHinE.dunniilayertosulphur(S)andnitrogen(N)combineddepositionwasstudied.Eliminatedthe
growthofage, themathematicalmodels:y =0.776 3+0.193 5x
1
-0.125 6x
2
-0.510 0x
1
x
2
-0.245 0x2
1
-0.597 5x2
2
describing
theaboverelationswereestablished.AsforthegrowthofDBH, theirsingleefect, marginaleffectandcombinedefectbetween
sulphur(S)andnitrogen(N)wereanalyzed.TheresultsshowedthatthebestcombinationwasS28.20 kg· hm-2· a-1 andN
64.78kg· hm-2· a-1 depositionrespectively, whichindicatedthatthegrowthofdiameteratbreastheightwouldreach0.84cm· a-1.
Keywords:Eucalyptusdunniplantation;sulphurandnitrogencombineddeposition;quardraticregressionrevolutiondesign;
growthofdiameteratbreastheight
由于工农业生产 、交通运输 、城市化进程加快 ,导致全球生态环境的逐步恶化 ,在北美 ,科学界和公众
对酸雨的广泛认识始于 20世纪 70年代初 ,美国和加拿大还实施了一系列的酸雨研究计划 ,如美国的 “全
球大气沉降计划 ”(nationalatmosphericdepositionprogram, NADP, 1978)、“国家酸雨评价计划 ”(national
acidprecipitationassessmentprogram, NAPAP, 1982)、加拿大的 “国家酸雨早期预警系统 ”(acidrainnation-
alearlywarningsystems, ARNEWS, 1984)等。自 Nihlgard[ 1]首次提出氮沉降对森林生态系统的危害以来 ,
这一课题在欧美己成为当今酸沉降研究领域的一个热点问题。 20世纪 80年代末 ,欧共体委员会资助了
“氮饱和试验”(nitrogensturationexperiments, NITREX)和 “欧洲森林生态系统实验操作”(experimentalma-
nipulationofforestecosystermsinEurope, EXMAN)两大研究项目 ,分别研究增加或减少大气氮沉降对欧洲
森林生态系统特别是针叶林生态系统的影响和改变周边大气沉降的化学组成或数量从而研究其对森林生
态系统的影响[ 2-3] 。挪威于 1991年启动跨学科研究项目 “氮—从山地到峡湾 ”,该项目旨在了解氮平衡
与氮沉降对挪威南部 Bjerkrein流域和 Auli流域造成的影响[ 4] 。
我国在 20世纪 70年代末 80年代初 ,酸沉降被确定为一个潜在的环境问题 ,但对氮沉降的起步较
晚 [ 5-6] 。鉴于我国氮沉降的严峻形势 , 2002年 10月在广东鼎湖山生物圈保护区选择亚热带代表性的森
林生态系统建立了永久性的试验样地 ,通过人工模拟方法系统地研究氮沉降对南亚热带森林生态系统结
构和功能及其机理[ 7-12] 。樊后保等 [ 13 -16]模拟氮沉降增加对亚热带杉木人工林植被养分状态 、土壤化学
DOI :10.13324/j.cnki.j fcf.2009.02.019
和生物活性 ,以及系统养分循环过程的影响 ,揭示杉木人工林生态系统对氮沉降及氮饱和的响应特征和机
制 ,以及酸化系统自我恢复的状况。然而关于硫沉降对森林生态环境的影响的研究国内外均鲜有报道 ,利
用试验设计方法系统阐述硫 、氮复合沉降对森林生态系统的影响的研究未见报道 。
福建省是酸雨严重的地区 [ 17] ,桉树(Eucalyptus)作为福建省主要的引种树种 ,具有抗寒性且较速生 ,
目前正显示出良好的发展势头 ,其对酸沉降的响应及其程度对进一步扩大种植面积及制订应对全球酸沉
降措施具有重要的理论指导意义。本研究采用野外模拟硫 、氮复合沉降的方法 ,研究邓恩桉(Eucalyptus
dunni)幼龄林林分胸径生长量对硫 、氮复合沉降的响应 ,并进一步探讨邓恩桉幼龄林对于外加硫 、氮的承
载能力 ,预期为我国桉树速生丰产林的可持续经营提供科学依据 ,为酸沉降增加环境下的亚热带阔叶林森
林管理提供理论基础 。
1 材料与方法
1.1 研究区概况
试验地位于福建省建阳市 ,位于武夷山南麓 ,地处北纬 27°06′-27°43′,东经 117°31′-118°38′,属亚
热带季风性气候 ,年平均气温 18℃,无霜期 282 d,年平均降雨量 170-240m,年平均日照数 1 802.7 h,土
壤以山地红壤为主[ 18] 。邓恩桉造林地为杉木 、马尾松采伐迹地 , 2006年 2月份种植邓恩桉 ,造林密度为
1 750株·hm-2 ,后期人工管理进行每年 3次的除草追肥 ,每次施肥 500g·株 -1 ,现存密度 1 650株· hm-2 ,
平均胸径 8.57 cm,平均树高 9.1 m。
1.2 试验设计
应用二因素正交旋转组合设计[ 19] ,参照 NITREX项目 [ 2]和段雷等 [ 20]对酸沉降临界负荷的研究 ,确定
模拟硫 、氮沉降的梯度和沉降量 ,其中二水平全因子试验点个数 mc=4,变量为零水平的重复试验点个数
m0 =8,与因子个数有关的参数 mr=4,空白样地 3个 ,共 19次试验。于 2007年 7月 23日在 3年生邓恩桉
人工纯林中随机设置立地条件基本相似的 20 m×20 m的试验样地 19块 ,各样地之间设置保护带 ,按硫
(x1)、氮(x2)沉降试验因素与水平的无量线性编码值(表 1)安排具体模拟试验。
表 1 试验因素与编码
Table1 Testfactorsandtheircoding (kg· hm-2· a-1)
因素 变化区间 编码值(r=1.414)+r 1 0 -1 -r
x1 16 44 38 22 6 0
x2 50 140 120 70 20 0
1.3 研究方法
2007年 7月建立样地后进行本底调查 ,对林木进行每木检尺 ,实测林木单株胸径(D)。于 8月开始每
月月初以溶液的形式对样地进行模拟硫 、氮复合沉降处理 ,试验持续 1 a。以硫酸钠(Na2SO4)作为外加硫
源 ,尿素 [ CO(NH2)2 ]作为外加氮源 ,将不同处理组合的硫酸钠及尿素溶解于 20 L水中 ,以背式喷雾器在
样地内人工均匀喷洒 ,对照喷施等量的清水。每 3个月进行每木检尺 ,记录林木单株胸径和树高 。
1.4 数据处理
采用 DPS2000数据处理软件进行数据处理 。对试验处理前各处理的胸径进行方差分析 ,结果表明各
处理间在林分平均胸径上不存在显著差异 ,可判断各处理林分生长状况基本一致 ,可利用处理 1a后的平
均胸径 ,计算每一试验处理下邓恩桉胸径生长量的年增长量 ,并运用二次正交旋转组合设计对胸径年增长
量进行回归分析 。
2 结果分析
2.1 胸径生长对硫 、氮沉降响应模型的建立与检验
由于方差分析表明试验前各处理之间林分平均胸径不存在显著差异 ,因此可利用邓恩桉幼林模拟硫 、
氮沉降处理 。 1 a后观察获得胸径生长调查数据 ,计算每一处理下的邓恩桉林分平均胸径Di(i=1、…、16)
·126· 福 建 林 学 院 学 报 第 29卷
和对照处理平均胸径D0 ,利用■Di=Di-D0计算不同试验处理下邓恩桉林分胸径生长量年增量■Di(i为
第 i处理)(表 2)。
表 2 试验方案与数据
Table2 Designoftheexperimentanditsdata
试验号 试验因素x1 x2
胸径年增长量
y 试验号
试验因素
x1 x2
胸径年增长量
y
1 1 1 -0.65 9 0 0 0.45
2 1 -1 0.95 10 0 0 0.60
3 -1 1 0.19 11 0 0 0.83
4 -1 -1 -0.25 12 0 0 0.98
5 r 0 0.58 13 0 0 0.75
6 -r 0 -0.26 14 0 0 0.92
7 0 r -0.49 15 0 0 0.99
8 0 -r -0.60 16 0 0 0.69
根据回归设计原理 ,试验所建立的数学模型为
y =f(xj)=b0 +∑2j=1bjxaj+∑i
j=1bjjx2aj (1)
经计算获得硫 、氮沉降下邓恩桉人工林胸径年增长量与硫 、氮沉降之间的二元二次回归方程
y =0.776 3 +0.193 5x1 -0.125 6x2 -0.510 0x1x2 -0.245 0x21 -0.597 5x22 (2)
对回归方程进行显著性分析(表 3), F1 =3.91
硫 、氮沉降下的邓恩桉胸径年增长量变化规律 。
表 3 方差分析
Table3 Varianceanalysis
变异来源 平方和 自由度 均方 F P
x1 0.299 5 1 0.299 5 4.366 6 0.063
x2 0.126 1 1 0.126 1 1.838 6 0.205
x21 0.480 2 1 0.480 2 7.000 8 0.025
x22 2.856 0 1 2.856 0 41.638 0 0.000
x1x2 1.040 4 1 1.040 4 15.167 9 0.003
回归 4.802 3 5 0.960 5 F2 =14.00 0.005
剩余 0.685 9 10 0.068 6
失拟 0.429 5 3 0.143 2 F1 =3.91 0.044
误差 0.256 4 7 0.036 6
总和 5.488 2 15
2.2 胸径生长对硫 、氮沉降响应数学模型的解析
2.2.1 硫 、氮沉降对胸径增长量的独立效应分析 对于式(2),将其中一个因素分别固定在 -1.414、 -1、
0、1、1.414水平上 ,得到邓恩桉胸径增长量与另一个因素的回归方程 ,按方程作图(图 1)。结果表明 ,硫沉
降对胸径增长量的效应因不同氮沉降编码水平而变化 。在较低的氮沉降水平时(编码 1、2、3),刚开始胸
径增长量较低 ,但随着硫沉降编码水平的增加 ,胸径增长量逐渐增加 ,最后趋于平缓。在较高氮沉降水平
时(编码 4、5),刚开始胸径增长量较高 ,但随着硫沉降编码值的水平增加而缓慢增加后逐渐减少。当氮沉
降编码值为零水平时 ,整体效应最好 ,而且变化较平缓 。以上说明低氮沉降情况下 ,硫沉降能够提高胸径
增长量 ,但在高氮沉降情况下 ,效应相反。氮沉降对胸径增长的作用因不同硫沉降水平而变化 ,在较低硫
沉降水平时(编码 1、2、3),氮沉降对胸径增长表现为正向效应 ,随着硫沉降编码水平的提高 ,正向效应逐
渐减弱 。而在较高硫沉降水平时(编码为 4、5),氮沉降对胸径增长量表现:氮编码值(-1.414, -1),表现
为正效应;氮编码值(-1, 0),正效应逐渐减弱转化为负效应;随着氮沉降编码水平继续增加 ,胸径增长量
迅速下降 ,负效应增强 。
·127· 第 2期 侯丹等:邓恩桉胸径生长对模拟硫 、氮复合沉降的响应
图 1 胸径增长量与硫 、氮沉降的回归关系
Figure1 RelationshipbetweenDBHandS, Ncombineddeposition
2.2.2 硫 、氮沉降对胸径增长量的交互效应 将一个因素固定为零水平进行降维分析 ,即可得到二因素
的交互作用回归子模型。
y=0.776 3+0.193 5 x1 -0.245 0 x21 (3)
y=0.776 3-0.125 6 x2 -0.597 5 x22 (4)
图 2 交互作用曲线
Figure2 Interactioncurve
根据以上方程做出每个因子对胸径增长量的交互作用曲
线(图 2)。硫沉降与氮沉降交互作用在低氮沉降(编码值为
-1.414、 -1、0)时为正交互作用 ,当氮沉降编码值超过 0时 ,
硫沉降与氮沉降交互作用转变为负交互作用。
2.2.3 硫 、氮沉降对胸径增长量的边际效应 边际效应分析
在经济学领域应用广泛 ,其实质是分析每增加一个单位的货
币投入所能带来的收益数量。现将其应用于分析酸沉降因子
(硫 、氮),每增加一个单位的量所导致胸径增长量的变化 ,由
此可推导出各种环境背景下环境对于酸沉降的容载力。
对零水平的回归子模型(3)和(4)分别求出 y对 x1、x2的
偏导数 ,可分别得到硫沉降和氮沉降的边际效应模型。
y x1 =0.193 5-0.490 0x1 (5)
y x2 =-0.125 6-1.195 0 x2 (6)
表 4 各因素边际效应值
Table4 Marginalefectvaluesofthefactors
因素 编码值-1.414 -1 0 1 1.414
硫 0.886 0.684 0.194 -0.297 -0.499
氮 1.564 1.069 -0.126 -1.321 -1.815
将硫 、氮各个水平的编码值分别代入相应
模型中 ,即可得到在特定条件下各因素不同水
平的边际效应(表 4、图 3)。从数值的绝对值
大小看 ,在整个浓度范围内 ,两因子的边际效应
值先减后增 ,硫 、氮在低浓度时 ,边际效应值逐
渐减小 ,对胸径增长的促进作用逐渐减弱 ,而在
高浓度时 ,对胸径增长抑制作用增强 。利用边际效应模型 ,令导数等于零 ,当 S编码值为 0.394 9(即纯硫
28.20kg· hm-2·a-1)、N编码值为 -0.105 1(即纯氮 64.78 kg· hm-2·a-1)时 ,邓恩桉胸径增长量可望
达到最大 0.84cm。硫沉降的容载力与前人提出的土壤硫沉降临界负荷值 32kg·hm-2·a-1比较接近 ,而
对于氮容载力较土壤的临界负荷值 25-40 kg· hm-2·a-1大[ 20] 。森林生态系统对于环境的承载能力较
强 ,不同林龄的树种表现出不同的反馈机制 ,研究表明 ,中氮处理(100 kg·hm-2·a-1)对亚热带 3种阔叶
树幼苗基径生长的促进作用最大[ 8] 。可见 ,树木幼龄阶段对于氮肥的需求量大。前人研究表明 ,与阔叶
林地相比针叶林地不容易达到氮饱和 [ 21-22] ,该邓恩桉幼龄林地原本为杉木 、马尾松林采伐迹地 ,邓恩桉作
·128· 福 建 林 学 院 学 报 第 29卷
图 3 边际效应曲线
Figure3 Marginalefectcurve
为速生树种 ,对林地养分需求量大 ,模拟氮沉降具
有施肥效应 。
3 结论与讨论
本研究能消除邓恩桉年龄增长对胸径生长量
的影响 ,较好地反应硫 、氮复合沉降对邓恩桉胸径
增长量的影响。研究表明 ,酸沉降对邓恩桉胸径
增长量有双面影响 ,在适宜范围内 ,硫 、氮沉降的
增加有利于胸径增长 ,但是随着沉降水平的增加 ,
胸径增长量不断减少 ,说明林地内硫与氮的含量
足够林木生长利用 ,继续外加硫源和氮源 ,会引起
负面效应。证明了当生态系统生产力不受硫与氮限制 ,或者说生态系统达到硫 、氮饱和时 ,外加硫 、氮处理
不但不增加生产力 ,反而减少森林的活力最终可能导致森林衰退 。本研究与温带森林进行的外加氮模拟
试验结果不同 ,本研究阔叶林胸径的生长对外加硫 、氮的响应速度较快(1a),这可能与研究林分的具体林
龄及其所处生长阶段有关 。而在丹麦 , 75年生的云杉林在喷加 NH4NO3 5 a后(35 kg·hm-2· a-1),外加
氮没有对树木生长产生显著的影响 , Gundersen[ 24]认为出乎意料的 4 a干旱影响了试验结果。同样 ,在瑞
典一针叶林集水区进行的外加氮模拟试验 ,持续 5 a的氮输入增加(约 40 kg· hm-2·a-1)也没有改变植
物生长的速度[ 25] 。在 Harvard的红松林和硬木林在进行 3 a的氮处理(喷加 NH4NO3和 150 kg·hm-2·
a-1的 N)后初级生产力才开始发生变化 [ 26] 。Gundersenetal[ 2]认为 ,温带森林植物生长对氮沉降的响应
较慢 ,一般需要 5 a以上。邓恩桉具有较强的速生性 [ 18] ,在幼龄阶段对水肥的吸收量较大 [ 23] ,本研究对象
是邓恩桉 3年生幼龄林 ,因此表现出对于外加硫 、氮反应比较敏感 ,响应速度较快 。
经过 1a的研究 ,结果表明 ,硫 、氮复合沉降对邓恩桉胸径的增长影响符合收益递减规律 ,适当增加硫
和氮输入量能够促进胸径的增长 ,但是持续的增加沉降量会造成负面影响。所以邓恩桉幼龄林胸径增长
量达到最大时的酸沉降输入值为纯硫 28.20 kg·hm-2·a-1 ,纯氮 64.78kg· hm-2· a-1。
运用现代回归设计理论 ,利用较小的试验规模 ,建立了硫沉降及氮沉降与邓恩桉幼龄林胸径年增长量
数学模型 ,对于进一步研究南方阔叶林地对输入硫 、氮的响应机制具有重要价值 。本研究既能够为我国桉
树速生丰产林的可持续经营提供科学依据 ,也能够为酸沉降增加环境下的亚热带阔叶林森林管理提供理
论基础 。由于野外模拟试验 ,影响试验因素不能够掌控 ,因此长期的观测以及对于不同地位级的阔叶林地
的试验有待于进一步研究 。
参考文献
[ 1] NihlgardB.Theammoniumhypothesis-anadditionalexplanationtotheforestdiebackinEurope[ J] .Ambio, 1985, 14:2-8.
[ 2] GundersenP, EmmetBA, KjonaasOJ, etal.Impactofnitrogendepositiononnitrogencyclinginforest, asynthesisof
NITREXdata[ J] .ForestEcologyandManagement, 1998, 101:37-55.
[ 3] WrightRF, RasmussenL.IntroductiontotheNITREXandEXMANprojects[ J] .ForestEcologyandManagement, 1998,
101:1-7.
[ 4] 曹军 , SembAK.挪威南部两流域氮硫和氯元素的大气沉降 [ J] .人类环境杂志 , 1997, 25(5):254-260;282.
[ 5] 樊后保 , 苏兵强 ,林德喜 , 等.杉木人工林生态系统的生物生态地球化学循环:氮素沉降动态 [ J] .应用与环境生物学报 ,
2000, 6(2):133-137.
[ 6] FanHB, HongW.EstimationofdrydepositionandexchangeinChinesefirplantations[ J] .ForestEcologyandManagement,
2001, 147:99-107.
[ 7] 方运霆 , 莫江明 , PerG, 等.森林土壤氮素转换及其对氮沉降的响应 [ J] .生态学报 , 2004, 24(7):1 523-1 531.
[ 8] 李德军 , 莫江明 , 方运霆 , 等.模拟氮沉降对三种亚热带树苗生长和光合作用的影响 [ J] .生态学报 , 2004, 24(5):
876-882.
[ 9] 李德军 , 莫江明 ,方运霆 , 等.模拟氮沉降对南亚热带两种乔木幼苗生物良及其分配的影响 [ J] .植物生态学报 , 2005, 29
(4):543-549.
·129· 第 2期 侯丹等:邓恩桉胸径生长对模拟硫 、氮复合沉降的响应
[ 10] 方运霆 ,莫江明 , 江远清 ,等.鼎湖山森林土壤渗透水酸度和无机氮含量对模拟氮沉降增加的早期响应 [ J] .热带亚热
带植物学报 , 2005, 13(2):123-129.
[ 11] 薛瑾花 ,莫江明 , 李炯 ,等.氮沉降对外生菌根真菌的影响 [ J] .生态学报 , 2004, 24(8):1 789-1 794.
[ 12] 徐国良 ,莫江明 , 周国逸 , 等.氮沉降下鼎湖山森林凋落物分解及与土壤动物的关系 [ J] .生态环境 , 2005, 14(6):
901-907.
[ 13] 袁颖红 ,樊后保 , 王强 , 等.模拟氮沉降对杉木人工林土壤有效养分的影响 [ J] .浙江林学院学报 , 2007, 24(4):
437-444.
[ 14] 樊后保 ,刘文飞 , 杨跃霖 , 等.杉木人工林凋落物分解对氮沉降增加的响应 [ J] .北京林业大学学报 , 2008, 30(2):
8-13.
[ 15] 刘文飞 ,樊后保 , 张子文 ,等.杉木人工林针叶养分含量对模拟氮沉降增加的响应 [ J] .应用与环境生物学报 , 2008, 14
(3):319-323.
[ 16] 樊后保 ,刘文飞 , 徐雷 ,等.氮沉降下杉木人工林凋落叶分解过程中 C、N元素动态变化 [ J] .生态学报 , 2008, 28(6):
2 546-2 553.
[ 17] 王文兴 ,丁国安.中国降水酸度和离子浓度的时空分布 [ J] .环境科学研究 , 1997, 10(2):1-6.
[ 18] 陈建忠.邓恩桉人工林幼林生长规律 [ J] .福建林学院学报 , 2007, 27(1):88-91.
[ 19] 洪伟 ,吴承祯.试验设计与分析———原理 、操作 、案例 [ M] .北京:中国林业出版社 , 2004:179-200.
[ 20] 段雷 ,郝吉明 , 谢绍东 ,等.用稳态法确定中国土壤的硫沉降和氮沉降临界负荷 [ J] .环境科学 , 2002, 23(2):7-12.
[ 21] 方运霆 ,莫江明 , 周国逸 ,等.鼎湖山主要森林类型植物胸径生长对氮沉降增加的初期响应 [ J] .热带亚热带植物学报 ,
2005, 13(3):198-204.
[ 22] MoJM, BrownS, PengSL, etal.Nitrogenavailabilityindisturbed, rehabilitatedandmatureforestsoftropicalChina[ J] .
ForestEcologyandManagement, 2003, 175:573-583.
[ 23] 李林.邓恩桉幼龄林施肥效应研究 [ J] .华东森林经理 , 2007, 21(2):9-14.
[ 24] GundersenP.EfectsofenhancednitrogendepositioninaspruceforestatKlosterhede, Demark, examinedbymoderate
NH4NO3 addition[ J] .ForestEcologyandManagement, 1998, 101:251-268.
[ 25] KjonaasOJ, StuanesAO, HuseM.EffectsofweeklynitrogenadditionsonNcyclinginaconiferousforestcatchment, Gard-
sjon, Sweden[ J] .ForestEcologyandManagement, 1998, 101:227-250.
[ 26] MagilAH, AberJD, HendricksJJ, etal.Biogeochemicalresponseofforestecosystemstosimulatedchronicnitrogendepo-
sition[ J] .EcologicalApplications, 1997, 7(2):402-415.
(责任编辑:江 英)
·130· 福 建 林 学 院 学 报 第 29卷