1、珍稀濒危植物长苞铁杉种群密度效应模型38 卷 4 期2002li7 月林业科学SCIENFIASIL-,ESINICAE珍稀濒危植物长苞铁杉种群密度效应模型吴承祯洪伟(福建农林大学林学院南平 353001)关键词:长苞铁杉,密度,密度效应ASTUDYoNTHEDENSITYEFFECTMoDELoFRAREANDENDANGEREDE,6ALD6 忸 RACZPoPULATIoNWuChengzhenHongWei(Forestry.CollegeojFujianAgriczdtureandForestryUniversityMmping353001)Abstract:Thedensityef
2、feetofrareandendangeredTsugalongibracteatapopulationwasanalyzedbythetheoreticalmodeldescribingthewholedensitychangeprocessofplantpopulation,i.e.N=exp(1,lnS+binS+C),whereVandSrepresentpopulationdensityandmeanbasalareaofT.1ongibracteataatbreastheightrespectively,and.b,careparanlclers.Then,theauthorspr
3、oposedadensityeffectmodelconsideringtheelfeetofdifferentsiteandassociatedspe(ies0【IdensityregulationofT.1ongibracteata.whichwasdescribedasN=exp(一 0.42426HInS一 5.60238H“InS 一6.75209H“.)+exp(0.204371nS 一 0.095858InS 一 0.14548),whereHandSrepresentedmeanheightofdominatlttreesandmeanbasalareaofassociated
4、speciesatbreastheight.ThesimulatedresultshowedthatthemodelfittedwellobserveddatafromT.1ongibracteatapopulation.Thereforeitisveryusefulinpractice,whichwillprovideascientificbasis)rdensitymanagementofT.1ongibracteatapopulation.Keywords:Tsugalongibracteata,Density,Densityeffect植物种群的调节是一个多因索的复杂过程,通过对种群的
5、出生,死亡和生长等过程的控制,使种群的数量维持在生境可以容纳的最大水平上.尽管影响种群动态的因素可分为密度制约因素和非密度制约因素,但在植物中,尤其对木本植物,主要是密度制约起着种群调节与控制的作用,这种调控的过程反映了种群对群落结构及相应环境条件的适应与利用,主要表现为植物种群密度效应(江洪 ,1992).种群密度效应一直是种群生态学领域研究的热点.前人在研究植物种群密度拥挤效应规律时,提出了大量理论或经验的表达式(Yadaeta1.,1963; 唐守正,1993;方精云,1995;Hozumi,1977;1983;Jobnetnz.,1994;张大勇等.1985;吴承祯等,2000c;20
6、00d;安藤贵,1968;洪伟等,1996; 陈辉等,1992),反映了植物种群数量动态受多种限制因素影响,如有限的资源,空间以及种内,种问植物个体对资源的竞争,导致种群数量发生变化,即为密度效应.美国 Rieneke 首先提出了林分密度竞争效应的乘幂式:N=KD 一(1)其中:J7v 为种群密度;D 为平均直径;n,为因树种而变化的常数.日本的吉良龙夫(Yadaeta1.,1963)发现草木植物密度与生长之问存在:V=KN(2)其中:为平均个体重量;N 为单位面积株数.1982 年 Silvertoan 在研究林分自然稀疏时,也同样提出这种模型.但模型(2)必须在林分疏密度相同时使用才有意义
7、.安藤贵认为林分生长过程中,不断分化出来的生长落后的被压木,因光照和营养不足,最后导致自然枯死(安藤贵,1968). 于是提出:收稿日期:2001.04.09.基金项 H:福建省自然科学基金资助项目(B0110026); 福建省科委重大资助项目(闽财指【2oo1338 号).l58 林业科学 38 卷V:aN一 b(3)立木密度与林分单株材积产量的密度效应关系也可用密度二次效应模型表示:0 一 c|7v(4)前人研究中较多地应用(2)式来研究森林自疏和密度效应规律,0Neill 和DeAngelis 研究娑罗双和柚木,赵学农(1991) 研究木果石栎均采用(2) 式.但是由于模型(2)必须在林
8、分疏密度相同时才有意义,因此,陈辉等(1992) 提出密度效应新模型:口(5)#其中:为单株平均材积;N 为单位面积密度;H 为林分的优势木平均高 ;a,卢,卢为模型参数.笔者(吴承祯等,2000d) 曾提出一个能描述植物种群密度变化全过程的密度效应模型:N=exp(alnS+binS+c)(6)其中:J7v 为单位面积株数(株?hm);S 为林分单株平均胸高断面积(rn2?株);a,b,c 为待定参数.长苞铁杉(Tsugalongibracteata)属松科铁杉属,为国家三级保护植物(国家环境保护局,1987),是珍贵稀有树种之一.其分布范围相对狭窄,散生于我国的贵州,湖南,广东,广西和福建
9、等地(林金星等,1995).目前长苞铁杉的残存个体数量有限,在福建省天宝岩自然保护区保存有一片较完好的长苞铁杉针阔混交林,面积约 72hm2 左右.因此,研究其种群生态学与群落生态学 ,对于保护珍稀濒危植物长苞铁杉具有一定的生态意义和理论价值.但有关长苞铁杉种群生态学与群落生态学研究报道较少(吴承祯等,2000a;2000b;2001),本文探讨长苞铁杉种群密度效应,旨在考察长苞铁杉种群密度调节的过程与机制,揭示其在不同自疏阶段中的密度变化规律.1 研究区概况福建省天宝岩自然保护区位于北纬 2555 一 25o58,东经 117D31 一 117D33.5,面积约 1976.5hm2,是戴云山
10、余脉,为中低山地貌,最高峰天宝岩 1604.8m.保护区年平均温度 23左右,最冷月(1 月)绝对最低温度一 11,最热月(7 月)绝对最高温度 40,全年10的活动积温在45205800左右,年降雨量 2000mm,空气相对温度 80%左右.土壤的垂直分布大致是海拔 800m 以下为红壤,8001350m 为黄红壤,1350m 以上为黄壤.长苞铁杉群落外貌整齐 ,色彩暗绿,林内阴暗湿润,土层较厚,表土质地松软,枯枝落叶层厚 35cm,为中亚热带山地垂直地带性植被类型.长苞铁杉主要伴生树种有甜槠(Castanopsiseyrei),杉木(Cunninghamialanceolata),木荷(S
11、chimasuperba),深山含笑(Micheliamaudiae),薯豆(Elaeocarpusjaponica),小叶青冈 (Cyclobalanopsismyrsinaefolia),冬青(1lexpurpurea)及猴头杜鹃(Rhododen dronsimiarum)等.2 材料与方法2.1 调查方法在保护区内选择生长发育正常的长苞铁杉天然混交林,依不同生境条件,不同发育阶段设置 8 块样地,每块样地面积 600m.对样地内林木进行每木检尺 ,记录样地所有个体的种名,胸径,树高,冠幅(起测径阶4cm), 同时每 100m2 选择 1 株长苞铁杉优势木计算优势木平均高;每一样地内随机
12、设置 6 个 1m1m 小样方调查长苞铁杉群落内灌木和幼苗幼树种类,高度及盖度,分布情况.8 块样地基本情况如表 1.2.2 密度效应模型考虑到各模型本身的特性及其适应性,本文应用模型(6)来研究珍稀濒危植物长苞铁杉种群密度效应规律,一方面动态描述长苞铁杉种群密度调节过程与机制,另一方面进一步检验本文所提出的密度效应模型的适应性.3 结果与分析由(6)式可知, 此方程为非线性方程,非线性方程一般采用麦夸方法拟合.但麦夸方法复杂,需要求偏导,因此计算量大而繁琐.本文采用遗传算法进行最优拟合(吴承祯等,1997;1998).第 4 期吴承祯等:珍稀濒危植物长苞铁杉种群密度效应模型 159(1)Ph
13、yllostachyspuhcens,(2)Cydoba/anopslsmyrsinaefolia,(3),.purpurea,(4)Sch/ma5l,(5)Cunn/ngham/a/anceo/ata,(6)Sorbushems 一,(7)Castanopsiseyrei,(8)M.maud.kle,(9)E.japonica,(10)Prumufok/enens/s,(11)P/nusta/wanens/s,(12)R.5,l,(13)施,0m 口,l,(14)L/tharpuspolystachya;(15)Castanops/sfat,(16)Choerospond/asax/ar/s
14、,(17)Cryptomeriafortunei,(18)s.珊 tsumu.根据遗传算法的基本思想,应用 BASIC 语言编制了长苞铁杉种群密度效应规律模拟的遗传算法计算机运行程序,经计算机运算得到长苞铁杉种群密度效应模型:N=exp(一 1.779595.330301nS 一 0.80921In2S)r:0.9830(7)式中,为长苞铁杉种群密度(株?hmI2);S 为平均单株胸高断面积(IIl2? 株).该式表明随着单位面积株数增加,种群平均单株胸高断面积逐步减小,较符合植物种群密度效应规律.但长苞铁杉种群生长受环境和年龄影响,不同生境,不同年龄间长苞铁杉种群生长必然存在差异,即使种群密
15、度一致,长苞铁杉种群平均胸高断面积也有差异.因此,尽管模型(7)模型相关指数较大,但从理论上说用(6)式研究长苞铁杉林各林分类型显然不适宜.考虑到种群密度随林分生长,立地条件的变化而变化,引入用长苞铁杉优势木平均高日表示的不同立地条件林分类型参数值:林业科学 38 卷0=口 lH,b=02H,c=03H3这样得到:N=exp0lH,InS+2HzInS+3(8)模型中集中反映林地立地质量,S 可以反映不同发育阶段林木的平均大小,因此从理论上讲该模型能客观地描述植物种群随立地,发育阶段而变化的密度效应规律.经遗传算法最优拟合,得到随优势木平均高变化的密度效应规律模型为:N:exp(一 0.416
16、69 月.“lnS5.57788no.25lnS 一 6.65283no25)r:0.9989(9)生境不同,种群密度调节控制的强度也不同,由于模型(9)综合考虑林分类型长苞铁杉优势木平均高,因此较好地反映了生境立地质量对种群调控能力的影响.根据计算,模型残差平方和为 1906.27,仅为模型(7)的 6.64%,说明用长苞铁杉优势木平均高改进长苞铁杉种群密度拥挤效应规律模型效果较好.由于长苞铁杉林由多树种组成,长苞铁杉种群动态不仅受其自身的种内竞争的影响,还受伴生树种对其种间关系的影响,为此,本研究进一步探讨长苞铁杉种群伴生树种对其密度效应规律的影响机制,即在其密度效应规律模型中增加伴生树种
17、为影响因子.研究表明伴生树种平均单株胸高断面积对长苞铁杉种群密度存在显着的影响,而伴生树种优势木平均高对其影响不显着.经遗传算法拟合得到考虑其伴生树种影响的长苞铁杉种群密度效应规律模型:N:exp(一 0.42426no 玛 lnS 一 5.60238no.InS 一 6.75209nof嘴.)+exp(0.204371nS 一 0.0958581nS 一 0.14548)r=0.9995(10)其中:s 为伴生树种平均单株胸高断面积(nl2?株).模型(10)的残差平方和为 930.21,仅为模型(7)的3.24%,模型(9)的 48.80%;利用模型(7),(9),(10)计算长苞铁杉种群
18、密度理论值 (表 2),表明模型(10)的效果最好.可知在探讨长苞铁杉种群密度效应时,考虑伴生树种对其影响是必要的.表 2 密度效应模型模拟结果Tab.2Thesimulatingresultsofdensityeffectmodel4 讨论密度效应反映由各种密度所形成的群体以及组成群体的个体之间相互作用的规律.长苞铁杉群落为多树种组成的异龄林,每个个体在群落内的地位和作用不同,对资源环境的要求也不相同.在相同资源条件下,生活力强,遗传性强的个体尽可能占据林冠层,获取尽可能多的资源,并控制其它个体对资源的获取,成为支配者;生长弱小,遗传性差的个体,因竞争能力远不及较大个体 ,只能占据剩余空间和
19、依靠剩余资源生存,处于林冠下层,生长发育受到抑制,从而被淘汰.而处于林冠层上层的长苞铁杉个体,由于生长发育旺盛,占据了大部分营养空间;加之与伴生树种间的种间竞争 ,导致处于劣势的长苞铁杉第 4 期吴承祯等:珍稀濒危植物长苞铁杉种群密度效应模型 l61种群个体生长受抑并逐渐被淘汰乃至死亡,从而使得林分中长苞铁杉种密度下降,其平均胸高断面积则随之上升.利用收集的长苞铁杉种群密度变化资料建立了其密度效应规律模型,并考虑用长苞铁杉优势木平均高和伴生树种胸高断面积指标改进其密度效应模型,改进后的模型相关指数为 0.9995,残差平方和仅为原模型的 3.24%,拟合效果较好,模型较好地描述了长苞铁杉种群密
20、度效应规律,可应用于指导其密度调节与控制.改进模型模拟表明,长苞铁杉种群密度变化机制既受其自疏作用的影响,也受其伴生树种他疏作用的制约.由于模型(6)中的参数.,b,c 反映了植物种群密度效应受树种特性影响的制约程度,而根据所建立的改进模型(10)式,不难发现长苞铁杉种群所对应的模型制约参数值的负值程度均比其伴生树种对应的模型制约参数值大,因此可以认为长苞铁杉种群密度效应受其自身的制约作用要比受其伴生树种的制约大得多.由此,我们有理由认为长苞铁杉种群密度变化过程 fj 的自疏作用大于其伴生树种的他疏作用,自疏是影响长苞铁杉种群密度减少的主要原因.本文建立的密度效应规律模型(10)式由于既考虑了
21、立地质量(平均优势高),又考虑了种群平均大小(平均胸高断面积), 因此适用于不同立地条件,不同类型林分的各发育阶段的长苞铁杉种群密度效应规律研究;不仅克服了前人模型只适用于某一林分整个生长过程的密度效应规律研究的局限性,而且充分考虑了种间竞争对密度效应规律的影响.但是由于该模型为非线性模型,因此必须选择适当的非线性拟合方法加以拟合,非线性模型的拟合方法有夸麦方法,三次设计法,改进单纯形法等方法,本文提出采用遗传算法实现模型的最优拟合.实例拟合结果表明,该方法具有较强的非线性拟合能力,是生态学,生物学及林学领域大量非线性模型的理想参数估计方法.参考文献安藤贵.同龄单纯林密度管理汇鞠寸为生憨学的研
22、究,林试研报,1968,210:1 152陈辉,何宗明,洪伟.杉木人工林密度效应模型的研究福建林学院,1992,12(3):277282方精云.种描述植物种群自然稀疏过程的经验模型.林业科学,1995,31(3):247253国家环境保护局,中国科学院植物研究所.中国珍稀濒危植物名录,第一册.北京:科学出版社,1987,68洪伟,吴承祯.闽北杉木人工林密度效应新模型.浙江林学院,1996,13(1):l520林金星,胡玉熹,王献溥等.中国特有植物长苞铁杉的生物学特性及其保护.生物多样性,1995,3(3):147152江洪.云杉种群生态学.北京:中国林业出版社,1992,5192唐守正同龄纯林
23、自然稀疏规律的研究.林业科学,1993,29(3):234241吴承祯,洪伟,用遗传算法改进约束条件下造林规划设计研究.林业科学,1997,34(2):133141吴承祯,洪伟.人工林经营过程密度最优控制研究.自然资源,1998,13(2):175180吴承祯,洪伟,吴继林等.珍稀濒危植物长苞铁杉种群生命表分析.应用生态,2000a,I1(3):333336吴承祯,洪伟,吴继林等珍稀濒危植物长苞铁杉的分布格局.植物资源与环境,2000b,9(1):3134吴承祯,洪伟.杉木林自疏过程密度调节规律的研究.热带亚热带植物,2000c,8(1):2834吴承祯,洪伟.杉木人工林自疏规律研究.林业科学
24、.2000d,36(4):97101吴承祯,洪伟,吴继林等.长苞铁杉群落种间竞争的研究.西北植物,2001,21(1):154158赵学农.哀牢山果石栎林种群调节与竞争的初步研究.植物生态学与地植物学,1991,15(2):183189张大勇,赵松龄.森林自疏过程中密度变化规律的研究.林业科学,1985,21(4):369374HozumiK?Eeologicalandmathematicalconsiderationsonself-thinninginever-agedpurestandsI,Meanplantweight-densitytrajectoryduringthecourseof
25、self-thinning.Bot,Mag.Tokyo,1977,90:165179lozumiK.Ecolcalandmathematicalconsiderationsonself-thinninginever-agedpurestandsIll,Effectofthelineargrowthfactor0nself-thinning.Bot.Mag.Tokyo,1983,96:171191JohnPM,WilliamFS,Standandtreedynamicsofuneven.agedPonderosapinesFor.Sci.,1994,40(2):289302YadaK,KiraT,OgawaHeta1.Self-thinninginovercrowdedpurestandsundercuhivatedandnaruralc0nditi0ns.J.Bi01107129
