-
土壤养分作为土壤最重要的属性之一,其含量及存在形态直接决定着土壤的肥力状况,进而影响着植被的生长及生理作用[1-2];在植被演替过程中,土壤中C、N、P等元素的循环相互影响[3-4],这必然会引起生态系统循环及其稳定性的变化。C:N和C:P反映了植物生长速度,C:N 在植被演替过程中相对稳定而C:P呈逐渐减小趋势[5],N:P可以定位对植物生长起限制性作用的营养元素[6],经研究,在植被演替过程中,会由受到N素限制转变为P素限制[7]。据统计,黄土高原地区每年因水土流失损失约3×107 t N、P、K等养分,严重制约了农林业的发展[8]。据研究,植树造林在增加植被覆盖度的同时,还可以有效改善土壤内养分状况、提高土壤肥力[9]。目前,国内外学者对植被生长与土壤养分的关系做了大量研究并取得了显著成果[10-13]。杨佳佳[14]等研究了黄土高原地区不同坡向对槐树林土壤中C、N、P含量的影响,发现叶片、枯落物、土壤C:N、C:P和N:P在阴、阳坡均无显著性差异(P < 0.05);张珂[15]等对阿拉善地区不同水分条件下的荒漠植物中C、N、P含量进行测定,探讨了C:N和N:P随水分梯度的变化规律;李玮[16]等通过对比不同年限的茶园表层土壤中养分含量,证明茶园年限越长养分含量越高,且差异显著;阎恩荣[17]等通过N:P值的变化分析了森林不同演替阶段的限制性元素;刘甲午[18]等综述了不同林龄华北落叶松林土壤的C、P、K变化特征,将林龄、坡度、降雨量等多种指标进行主成分分析,证明林龄对土壤养分含量的变化影响最大。但是,前人的研究多为单一林分下土壤的养分含量特征,或者仅仅对表层土壤中养分含量的变化规律进行定性研究,缺乏对林分在不同林龄条件下各土层间土壤养分含量变化的综合定量研究。
因此,本研究选取了黄土高原地区3种主要的人工植被(山杨、油松、侧柏),于2011年和2016年分别对其林内各土层土壤有机碳、全N、全P含量进行了分析,研究了林分在不同林龄以及不同土层间的C、N、P化学计量关系,以期对该区域森林经营及管理进行一定的指导,也可以为树种立地条件的选择提供理论依据。同时,研究林分在不同林龄条件下养分含量的变化,可以更确切地衡量退耕还林工程对改善生态环境的成效如何。
HTML
-
试验地分布在临汾市吉县境内,地处黄河中游东岸,山西省西南部,地理坐标为东经110°30′~110°43′,北纬36°10′~36°19′。地势东高西低,海拔落差大,最高海拔(高天山)为1 820 m,最低海拔(黄河畔)为450 m,试验地海拔位于1 200~1 488 m。县境内山峦起伏、沟壑纵横、地形复杂,试验地土质主要是黄土和褐土。吉县境内气候四季分明,光照充足,日照时数约为2 538 h·a-1,大于10℃的有效积温为3 361.5℃。无霜期年平均172 d,年均气温10.2℃,年均日较差11.5℃,年均降水522.8 mm。属暖温带大陆性季风气候,全年降雨集中在夏季。试验地植物资源比较丰富,且大多为人工林,乔木以油松(Pinus tabuliformis Carr.)、落叶松(Larix gmelinii(Rupr.)Kuzen)、刺槐(Robinia pseudoacacia Linn.)、侧柏(Platycladus orientalis(L.)Franco)、山杨(Populus davidiana Dode.)、白皮松(Pinus bungeana Zucc.)、辽东栎(Quercus liaotungensis Koidz.)、榆树(Ulmus pumila Linn.)为主,灌木有小叶锦鸡儿(Caragana microphylia Lam.)、三裂绣线菊(Spiraea trilobata Linn.)、山楂叶悬钩子(Rubus crataegifolius Bge.)等,草本植物主要包括画眉草(Eragrostis pilosa(L.)Beauv)、铁杆蒿(Artemisia sacrorum Ledeb.)、黑麦草(Lolium perenne Linn.)等,试验地概况如表 1。
林种Species 土壤类型Soil type 密度Density/(株·hm-2) 坡度Declivity/(°) 海拔Altitucle/m 调查时间Time of investigation 林龄Age of stand/a 平均树高Mean treeheight/m 平均胸径Mean DBH/cm 郁闭度Crown density 山杨P.davidiana 褐土 1 125 2~8 1 385~1 452 20112016 1015 9.0613.88 14.5816.59 0.600.70 油松P.tabuliformis 黄土 1 425 6~17 1 321~1 348 20112016 1520 8.099.52 12.9814.68 0.700.75 侧柏P. orientalis 褐土 1 650 10~25 1 225~1 240 20112016 4045 13.3516.71 20.5523.55 0.800.85 Table 1. Basic state of different forest types in the sample plots
-
不同林分类型和林龄对土壤全N含量的影响程度不同,如图 1所示。在0~100 cm土层中,各林地之间的平均土壤全N含量差异显著(P < 0.05),山杨林(10 a)、油松林(15 a)、侧柏林(40 a)分别为2.67、2.71、1.99 g·kg-1;5年后,山杨林(15 a)、油松林(20 a)、侧柏林(45 a)平均土壤全N含量分别为2.92、2.91、2.14 g·kg-1,增长速率分别为0.05、0.04、0.03 g·kg-1·a-1。在0~50 cm土层中,山杨林(10 a)、油松林(15 a)、侧柏林(40 a)平均全N含量为分别为2.97、3.05、2.27 g·kg-1,山杨林(15 a)、油松林(20 a)、侧柏林(45 a)平均全N含量分别为3.23、3.25、2.44 g·kg-1,5年内山杨、油松和侧柏林地内平均全N含量分别增长了9.4%、7.4%、7.5%;显然,林内平均全N含量均随着林龄的增加而增加,且山杨林5年内全N增长速率大于油松林和侧柏林。
由图 1可以看出,林内土壤N素集中在表层土壤:在山杨林(10 a)、山杨林(15 a)、油松林(15 a)、油松林(20 a)、侧柏林(40 a)、侧柏林(45 a)林地内,0~30 cm土层中全N含量分别占0~100 cm土层的75.3%、76.4%、73.9%、73.4%、72.4%、71.6%;而30~100 cm土壤含N量急剧减少,占比不到30%。同时,侧柏林(40 a)各土层土壤全N含量差异不显著,侧柏林(45 a) 0~10 cm土层土壤全N含量仅与50~100 cm差异显著,而山杨和油松林内0~10 cm土层土壤全N含量与10~100 cm土层均有显著性差异(P < 0.05)。
-
不同林种及林龄对土壤含P量的影响,如图 1所示。在0~100 cm土层中,各林地之间土壤的平均全P含量差异显著(P < 0.05):山杨、油松、侧柏林内土壤全P含量平均值为0.40、0.46、0.36 g·kg-1,这与林地之间N素分布规律基本一致。在0~50 cm土层中,山杨林(10 a)、油松林(15 a)、侧柏林(40 a)平均全P含量分别为0.43、0.49、0.42 g·kg-1,5年后,山杨林(15 a)、油松林(20 a)、侧柏林(45 a)平均全P含量分别为0.48、0.55、0.51 g·kg-1,5年内山杨、油松、侧柏林地内0~50 cm土层中土壤全P含量分别增长了11.6%、12.2%、21.4%。即:土壤中全P含量随林龄增加而增加;0~100 cm土层侧柏林内全P含量最低,而0~50 cm土层侧柏林内全P含量增长最快,说明侧柏主要增加中上层土壤中P素、对深层土壤基本没有影响。
林内土壤含P量随土层加深逐渐降低,也呈现出表层土壤集中分布的现象,如图 1所示。山杨林(10 a)、山杨林(15 a)、油松林(15 a)、油松林(20 a)、侧柏林(40 a)、侧柏林(45 a) 0~30 cm土层中全P含量占0~100 cm土层的百分比分别为79.3%、79.5%、71.2%、71.0%、83.6%、83.9%,30 cm以下土壤中P素极少。山杨和侧柏林地内0~10 cm土层全P含量与10~100 cm土层差异显著(P < 0.05),而油松林(20 a)中0~20 cm与20~100 cm间差异显著,油松林(15 a)各土层间无显著影响。就林内土壤全磷的P均值而言,各林地土壤含P量均较低,根据全国土壤普查结果的分级标准[22],该试验区土壤P素处于4~5级水平。
-
由表 2可以看出,各林地平均土壤有机碳含量差异显著(P < 0.05),且同一林分类型有机碳含量随着林龄增加而增加,在100 cm土层内平均有机碳含量为:山杨林(15 a)(17.51 g·kg-1)>山杨林(10 a)(15.93 g·kg-1)>油松林(20 a)(12.07 g·kg-1)>油松林(15 a)(11.28 g·kg-1)>侧柏林(45 a)(10.12 g·kg-1)>侧柏林(40 a)(9.63 g·kg-1)。5年内,山杨、油松、侧柏林地0~100 cm土层土壤平均有机碳含量分别增长了1.58、0.79、0.49 g·kg-1,0~30 cm土层含碳量分别增长了2.12、1.00、0.67 g·kg-1,0~10 cm土层平均碳含量分别增长了3.07、1.37、0.85 g·kg-1,显然林分生长改良表层土壤碳含量效果好,但是对深层土壤影响较小。总体来看,林分对中上层土壤有机碳的改善能力为山杨>油松>侧柏。
土层厚度Soil thickness /cm 有机碳含量Organic carbon content/(g· kg -1) 山杨林(10a)Populus davidiana 山杨林(15a)Populus davidiana 油松林(15a)Pinus tabulaeformis 油松林(20a)Pinus tabulaeformis 侧柏林(40a)Platycladus orientalis 侧柏林(45a)Platycladus orientalis 0~10 37.06±2.35a 40.14±3.44a 21.52±6.21a 22.90±1.82a 16.07±7.22a 16.92±5.40a 10~20 13.89±1.27b 15.82±1.95b 13.49±6.85b 14.49±1.62b 13.35±7.78a 14.18±1.41ac 20~30 12.33±4.19bc 13.69±1.53bc 9.94±3.21bc 10.58±0.48bc 9.85±5.81ac 10.47±2.48bc 30~50 10.27±1.53bc 10.91±0.46bc 6.97±5.35c 7.54±0.23c 6.41±4.34bc 6.64±3.60bd 50~100 6.08±0.72c 6.99±0.48c 4.47±4.04c 4.84±0.08c 2.16±0.05b 2.39±0.06d 平均 15.93±2.01A 17.51±1.57B 11.28±5.13C 12.07±0.85D 9.63±5.04E 10.12±2.59F 注:表中数据为平均值±标准差;同列数据后不同小写字母表示同一林地不同土层间有机碳含量差异显著(P < 0.05),不同大写字母表示不同林地之间平均有机碳含量差异显著(P < 0.05)。 Table 2. Soil organic carbon content in different soil layers of forest lands
林内有机碳的分布也具有表聚性,在0~30 cm土层内,山杨林(10 a)、山杨林(15 a)、油松林(15 a)、油松林(20 a)、侧柏林(40 a)、侧柏林(45 a)有机碳含量分别占0~100 cm土层的79.5%、79.6%、79.7%、74.5%、82.2%、82.3%,说明在30 cm以下有机碳累积速度缓慢,处于非常低的状态,这与土壤氮素和磷素含量的分布规律相一致,都呈现出由表层向深层逐渐迁移的动态过程。在0~100 cm土层内有机碳含量随着深度增加而减小:山杨和油松林内0~10 cm土层土壤有机碳含量显著(P < 0.05)高于10~100 cm土层;侧柏林0~20 cm土层内碳素分布较均匀,但明显高于30 cm以下土层碳含量。
-
如图 2所示,不同林内土壤碳储量均随土层深度逐渐减小,各林分土壤50~100 cm土层碳储量占土壤剖面碳储量比例为山杨林(10 a)(26.9%)、山杨林(15 a)(25.4%)、油松林(15 a)(24.9%)、油松林(20 a)(29.0%)、侧柏林(40 a)(17.1%)、侧柏林(45 a)(18.1%),均不超过0~100 cm土层土壤总碳储量的30%;表层(0~10 cm)土壤碳储量显著高于其他土层,各林分之间表现为山杨林(15 a)(56.42 t·hm-2)>山杨林(10 a)(51.36 t·hm-2)>油松林(20 a)(32.05 t·hm-2)>油松林(15 a)(26.69 t·hm-2)>侧柏林(45 a)(21.70 t·hm-2)>侧柏林(40 a)(21.69 t·hm-2),这与土壤含碳量分布规律相一致。
经过5年的时间,山杨林、油松林、侧柏林0~100 cm土层内土壤碳总储量分别由150.35、101.42、85.63 t·hm-2增加到171.52、125.15、87.84 t·hm-2,相当于山杨、油松、侧柏分别为该区固碳21.17、23.74、2.21 t·hm-2。0~30 cm土层内,山杨、油松、侧柏林地碳储量分别增加了12.58、10.05、0.97 t·hm-2,所以试验区山杨和油松比侧柏更有利于土壤中有机碳含量的提高。
-
不同林地土壤中C:N、N:P、C:P值随着土层深度变化规律如图 3所示。山杨林的C:N值为4.14~6.71,均值为5.72,变异系数为15.5%;油松林的C:N值为3.10~4.85,均值为3.90,变异系数为18.0%;侧柏林的C:N值为2.51~6.14,均值为4.40,变异系数为31.6%,3种林分类型都是中等变异。研究0~10 cm土层土壤C:N值,可知山杨林(15 a)(6.71)>山杨林(10 a)(6.46)>侧柏林(40 a)(6.14)>侧柏林(45 a)(5.67)>油松林(20 a)(4.85)>油松林(15 a)(4.83),表明山杨林和油松林表层土壤C:N值随林龄增大而增大,侧柏林则相反,且差异显著(P < 0.05)。
山杨林的N:P值为4.71~10.71,均值为7.89,变异系数为19.9%;油松林的N:P值为5.41~8.45,均值为6.62,变异系数为14.0%;侧柏林的N:P值为2.96~10.43,均值为8.01,变异系数为34.5%。侧柏林(45 a)土壤N:P值随着土层加深逐渐增大且0~10 cm土层N:P值显著低于10~100 cm土层,山杨和油松林内N:P值在不同土层之间差异不显著。随着林龄的增加,山杨林、油松林、侧柏林土壤表层(0~10 cm)的N:P值分别减少了0.36、0.10、0.12,表明在养分循环的过程中,林分有利于改善晋西黄土区磷素严重缺乏的状况。
山杨林的C:P值为29.34~67.87,油松林的C:P值为17.40~35.02,侧柏林的C:P值为18.03~50.43,总体来看山杨林土壤C:P值显著高于油松林和侧柏林。在0~10 cm土层中,C:P值为山杨林(10 a)(47.03)>山杨林(15 a)(45.49)>油松林(15 a)(35.02)>油松林(20 a)(34.69)>侧柏林(40 a)(18.03)>侧柏林(45 a)(16.99),即随着林龄增加C:P值逐渐减小,这也说明林分生长对改善晋西黄土区养分循环以及缓解其磷素严重缺乏的状况有明显作用。