馬尾松初級種子園土壤養(yǎng)分豐缺診斷與評價(jià)1)

李鵬 楊章旗 陳虎梁遠(yuǎn)毅 梁賢禎 蘇沃榜

(廣西馬尾松工程技術(shù)研究中心(廣西壯族自治區(qū)林業(yè)科學(xué)研究院)，南寧，530002)(廣西藤縣大芒界種子園)

土壤肥力是土壤質(zhì)量的重要組成部分，用于表征土壤供給作物生長各種養(yǎng)分元素的綜合能力，土壤養(yǎng)分是肥力的重要構(gòu)成因素，容易檢測和控制[1-2]。土壤養(yǎng)分豐缺診斷是針對指定區(qū)域內(nèi)土壤養(yǎng)分的限制性因子和養(yǎng)分失衡等做出的判斷[3]，主要目的是為了更好的利用土地資源，為科學(xué)施肥、林分調(diào)整以及制定和防止養(yǎng)分貧瘠化的措施和策略提供科學(xué)依據(jù)[4]。在明晰土壤養(yǎng)分障礙因子的基礎(chǔ)上，進(jìn)行區(qū)域內(nèi)土壤肥力綜合評價(jià)，對于挖掘土壤肥力潛力，提高生產(chǎn)力，實(shí)現(xiàn)土地資源集約化管理和利用具有重要意義。

馬尾松(Pinusmassoniana)是我國重要鄉(xiāng)土樹種，因其具有速生、耐貧瘠、用途廣等優(yōu)良特性，一直是我國南方主要工業(yè)材用和脂用樹種[5]。近年來，我國在馬尾松種子園建設(shè)[6-8]、經(jīng)營管理[9]和育種[10-11]等方面已做了大量研究，并取得了重要成果，為種子園的科學(xué)營建和管理提供了支撐。種子園作為林木遺傳育種的重要基地，具有良種生產(chǎn)、新種質(zhì)創(chuàng)制等重要功能[12]。馬尾松在進(jìn)入盛果期時(shí)，土壤養(yǎng)分貧瘠化是限制馬尾松雌雄球花芽分化、結(jié)實(shí)率和種子產(chǎn)量的重要因子[13]，及時(shí)補(bǔ)充生長所需的養(yǎng)分元素，調(diào)節(jié)生長和各器官養(yǎng)分平衡，對提高其產(chǎn)量具有重要作用。秦雪等[14]對都勻馬尾松種子園的土壤調(diào)查研究表明，土壤養(yǎng)分各元素間存在一定促進(jìn)和抑制作用；潘瓊?cè)氐萚15]對馬尾松1代種子園進(jìn)行母樹和土壤營養(yǎng)診斷，表明磷、鈣、鉬、鎂、硼是影響馬尾松種子園產(chǎn)量的主導(dǎo)因子；陸梅等[16]和陳元鎮(zhèn)等[17]對馬尾松2代種子園土壤營養(yǎng)診斷表明，氮、鎂、鉬等營養(yǎng)元素是影響馬尾松球果產(chǎn)量的主導(dǎo)因子。雖然以上研究豐富了馬尾松種子園的土壤貧瘠化研究體系，但仍有不足，一是研究主要集中于貴州地區(qū)，忽視了不同立地條件馬尾松林地本身土壤理化性質(zhì)存在差異，土壤養(yǎng)分的豐缺狀況和障礙因子并不相同；二是研究主要集中于單因子土壤養(yǎng)分的診斷，而未反映區(qū)域土壤養(yǎng)分的綜合肥力水平，不利于土壤養(yǎng)分的管理和集約化管理。因此，本研究以桂東馬尾松1代種子園為研究對象，開展不同立地類型(坡度、坡向、坡位)土壤養(yǎng)分豐缺診斷及綜合評價(jià)，診斷識別土壤養(yǎng)分的障礙(或主導(dǎo))因子及表征土壤肥力的最小數(shù)據(jù)集指標(biāo)，綜合評價(jià)區(qū)域內(nèi)土壤肥力指數(shù)，為土壤養(yǎng)分貧瘠障礙因子消除、科學(xué)配方施肥，提高種子園產(chǎn)量，實(shí)現(xiàn)土地資源集約化管理和利用提供技術(shù)支持。

1 研究區(qū)概況

研究區(qū)位于廣西藤縣大芒界種子園(111°6′44″～111°10′25″E，23°2′53″～23°5′14″N)，屬亞熱帶季風(fēng)氣候，日照充足，氣候溫暖，雨量充沛，年均降水量1 472 mm，年蒸發(fā)量1 057 mm，年平均氣溫21.3 ℃，極端最低溫-3.0 ℃，極端最高溫39.5 ℃，年均日照時(shí)間1 758 h，無霜期323 d。種子園于1978年定植砧木，1981年采用124個(gè)無性系進(jìn)行嫁接，面積53 hm2，共分為8個(gè)大區(qū)，32個(gè)小區(qū)，種植密度240株/hm2，保存密度204株/hm2，株行距7 m×7 m。2015年進(jìn)行截頂、拉枝定型，之后每年對母樹進(jìn)行適當(dāng)修剪。土壤為紫色砂頁巖發(fā)育而成的赤紅壤。

2 研究方法

2.1 土壤樣品的采集與測定

根據(jù)種子園的基本調(diào)查情況，選擇不同立地類型(坡度、坡向、坡位)的42棵母樹作為目標(biāo)樹，距離母樹周圍1 m處隨機(jī)選擇9個(gè)采樣點(diǎn)，用土鉆法采集0～30 cm的土壤樣品，采集的土樣混合均勻。采用四分法，取其中一份約500 g裝入無菌樣品袋帶回實(shí)驗(yàn)室，并放入專用土樣風(fēng)干盤中，攤成2～3 cm薄層，室溫下自然風(fēng)干，避免陽光直射或暴曬。去除植物殘?bào)w、土壤動(dòng)物、根系、石塊，碾碎、研磨過篩(2 mm)，用于土壤化學(xué)性質(zhì)測定。

根據(jù)林業(yè)行業(yè)標(biāo)準(zhǔn)(LY/T 1210～1275—2016)，pH采用酸度計(jì)法測定(m(土)∶V(水)=1.0 g∶2.5 mL)；采用重鉻酸鉀高溫外熱法測定有機(jī)質(zhì)的質(zhì)量分?jǐn)?shù)；采用堿解擴(kuò)散法測定速效氮的質(zhì)量分?jǐn)?shù)；采用Mehlich 3浸提測定有效磷、速效鉀、鈣、鎂、鐵、鋅的質(zhì)量分?jǐn)?shù)；采用甲亞胺比色法-間斷化學(xué)分析儀測定速效硼的質(zhì)量分?jǐn)?shù)。

2.2 種子產(chǎn)量的測定

所有目標(biāo)樹進(jìn)行統(tǒng)一訂牌、編號，采用“全收獲法”在每株目標(biāo)樹的樹冠東、西、南、北4個(gè)方向收集所有球果稱質(zhì)量；具體每個(gè)方向隨機(jī)收集5粒球果，每株目標(biāo)樹收集20個(gè)球果帶回實(shí)驗(yàn)室分析稱質(zhì)量，經(jīng)洗滌置于烘箱105 ℃殺青2 h后65 ℃烘干至恒質(zhì)量，去殼后測算其含水率，計(jì)算種子產(chǎn)量。

2.3 土壤養(yǎng)分豐缺分級標(biāo)準(zhǔn)

根據(jù)國標(biāo)西南區(qū)耕地質(zhì)量等級劃分標(biāo)準(zhǔn)(GB/T 33469—2016)，參照孫波等[18]和曹繼釗等[19]對南方紅壤區(qū)土壤養(yǎng)分貧瘠劃分標(biāo)準(zhǔn)，結(jié)合廣西馬尾松人工林土壤實(shí)際情況，制定出廣西低山丘陵區(qū)土壤養(yǎng)分豐缺的分級標(biāo)準(zhǔn)(見表1)，對馬尾松種子園土壤養(yǎng)分狀況進(jìn)行豐缺狀況分析。

表1 土壤養(yǎng)分豐缺分級標(biāo)準(zhǔn)

2.4 數(shù)據(jù)分析與處理

2.4.1 土壤養(yǎng)分障礙因子的確定

參考陳海艷[20]和陳元鎮(zhèn)[17]等的方法，采用線性逐步回歸方程(公式1)確定土壤養(yǎng)分單一效應(yīng)作用，二次多項(xiàng)式逐步回歸方程(公式2)確定元素間的交互效應(yīng)，從而進(jìn)行土壤養(yǎng)分障礙因子的確定。其診斷模型如下：

(1)

(2)

式中：Y為種子產(chǎn)量(因變量)，Xi表示第i個(gè)土壤養(yǎng)分指標(biāo)(自變量)，k0表示常數(shù)項(xiàng)系數(shù)，kij表示一次項(xiàng)常數(shù)項(xiàng)系數(shù)，kii表示二次項(xiàng)常數(shù)項(xiàng)系數(shù)，m表示指標(biāo)的個(gè)數(shù)。

2.4.2 土壤肥力最小數(shù)據(jù)集的確定

首先對所選擇的指標(biāo)進(jìn)行主成分分析，根據(jù)分析結(jié)果，提取特征值≥1的主成分，在荷載矩陣中，選擇每個(gè)主成分下的因子荷載>0.6的指標(biāo)，并納入各主成分指標(biāo)；然后對每個(gè)主成分中高荷載指標(biāo)進(jìn)行皮爾森相關(guān)性分析，若相關(guān)性系數(shù)r≤0.6則全部選入，若r>0.6則考慮綜合荷載值的大小，選擇綜合荷載值較大的指標(biāo)和其在10%范圍內(nèi)的指標(biāo)，作為最大影響因子的候選指標(biāo)；最后，對選擇的候選指標(biāo)進(jìn)行綜合對比，去掉重復(fù)，選入最大影響因子。

綜合荷載值的計(jì)算公式為：

(3)

式中：Nij是第i個(gè)變量在特征值≥1的前j個(gè)主成分的綜合荷載值；Rij為第i個(gè)變量在第j個(gè)主成分的荷載；λj為第j個(gè)主成分的特征值。

2.4.3 土壤綜合肥力評價(jià)

(1)土壤養(yǎng)分的標(biāo)準(zhǔn)化，運(yùn)用隸屬度函數(shù)對數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化(歸一化)處理，隸屬度公式為：

(4)

式中:a為函數(shù)的最小值;b為函數(shù)的最大值;x為指標(biāo)的測量值。

(2)權(quán)重的確定，根據(jù)綜合荷載值(公式3)的計(jì)算結(jié)果，確定各項(xiàng)養(yǎng)分指標(biāo)的權(quán)重。

(3)綜合指數(shù)的計(jì)算，采用加權(quán)隸屬度值公式，計(jì)算綜合肥力指數(shù)(I)，其計(jì)算公式為：

(5)

式中：Ij表示第j種類型土壤綜合肥力指數(shù);Wi和Ni分別表示第i個(gè)指標(biāo)權(quán)重和隸屬度值。

采用Excel 2016和IMB SPSS Statistics 20.0軟件對數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析。

3 結(jié)果分析

3.1 土壤養(yǎng)分含量和種子產(chǎn)量及其分布特征

由表2可知，土壤pH平均4.41，變化幅度3.91～5.18，變異系數(shù)最小為7.59%。有機(jī)質(zhì)質(zhì)量分?jǐn)?shù)為2.53～67.05 g·kg-1，平均25.68 g·kg-1，變異系數(shù)55.61%。速效養(yǎng)分元素鋅、鈣、鉀和鐵得質(zhì)量分?jǐn)?shù)變異系數(shù)較大，分別為86.74%、69.82%、57.77%和54.69%，均在50%以上，屬于強(qiáng)變異，表明這些元素隨立地因子變化的敏感性較強(qiáng)，其平均質(zhì)量分?jǐn)?shù)分別為0.87、103.01、35.15和29.7 mg·kg-1；速效養(yǎng)分元素氮、磷、鎂和硼的質(zhì)量分?jǐn)?shù)的變異系分別為41.78%、41.89%、35.97%和39.75%，變異范圍35%～50%，屬中等變異性，其平均質(zhì)量分?jǐn)?shù)分別為120.7、1.56、6.93和0.32 mg·kg-1。單株種子產(chǎn)量平均為150.24 g，變化幅度為80.4～311 g，變異系數(shù)為37.51%，屬中等變異。

由圖1可知，馬尾松種子園土壤養(yǎng)分和種子產(chǎn)量均符合正態(tài)分布特征。

表2 馬尾松種子園土壤養(yǎng)分和種子產(chǎn)量統(tǒng)計(jì)

3.2 土壤養(yǎng)分豐缺分級診斷

由表3可知，根據(jù)土壤養(yǎng)分豐缺分級標(biāo)準(zhǔn)對42個(gè)土樣進(jìn)行土壤養(yǎng)分豐缺分級診斷，土壤pH均小于5.5，土壤pH在酸性4.5～5.5范圍的為33.33%，土壤強(qiáng)酸性環(huán)境(pH<4.5)占66.67%，表明馬尾松種子園土壤主要屬于強(qiáng)酸性土壤且酸化較嚴(yán)重；38.1%的土壤有機(jī)質(zhì)處于較高水平，19.05%的土壤有機(jī)質(zhì)處于中等水平，16.66%的土壤有機(jī)質(zhì)處于適宜水平，19.05%的土壤有機(jī)質(zhì)處于低水平，7.14%的土壤有機(jī)質(zhì)處于極缺水平，表明種子園土壤有機(jī)質(zhì)質(zhì)量分?jǐn)?shù)比較適宜或具有較高水平；土壤速效氮豐富，僅有4.76%的土壤處于低水平；有效磷缺乏嚴(yán)重，16.67%的土壤有效磷處于極缺水平，71.43%的土壤有效磷處于低水平；速效硼缺乏嚴(yán)重，28.57%的土壤速效硼處于極缺水平，61.9%的土壤速效硼處于低水平；速效鉀、速效鈣和速效鎂的豐缺規(guī)律類似，均有21.43%的土壤處于極缺水平，50.00%、47.62%和64.29%處于低水平；速效鐵質(zhì)量分?jǐn)?shù)比較均衡，只有33.33%的土壤處于低水平；45.24%的土壤速效鋅處于極缺水平，23.81%的土壤速效鋅處于低水平。

頻數(shù)表示出現(xiàn)在指標(biāo)數(shù)據(jù)范圍內(nèi)的樣點(diǎn)數(shù)量。

表3 種子園土壤養(yǎng)分豐缺分級結(jié)果

3.3 土壤養(yǎng)分障礙因子的確定

由表4可知，根據(jù)土壤養(yǎng)分的診斷模型，構(gòu)建種子產(chǎn)量與土壤養(yǎng)分指標(biāo)的線性和二次多項(xiàng)式的逐步回歸方程，線性回歸方程中土壤養(yǎng)分單一效應(yīng)的主導(dǎo)因子是N、P、Mg和Zn，二次多項(xiàng)式逐步線性回歸方程中元素交互作用的主導(dǎo)因子是土壤有機(jī)質(zhì)、P、N×Mg、K×Zn和Ca×Mg。

表4 線性和二次多項(xiàng)式的逐步回歸方程

由表5可知，對比土壤養(yǎng)分的各主導(dǎo)因子的t值及其顯著性水平，單一元素效應(yīng)方面，土壤N、P、Mg和Zn是影響種子產(chǎn)量的主導(dǎo)因子，影響程度由大到小的順序?yàn)椋篘、Mg、P、Zn；元素間交互作用方面，N×Mg、K×Zn和Ca×Mg等養(yǎng)分效應(yīng)顯著，能顯著促進(jìn)種子產(chǎn)量，元素間交互效應(yīng)由大到小的順序?yàn)椋篘×Mg、K×Zn和Ca×Mg。

表5 各主導(dǎo)因子的回歸參數(shù)

3.4 最小數(shù)據(jù)集的確定

由表6可知，對10個(gè)土壤養(yǎng)分指標(biāo)進(jìn)行主成分分析，測度(KMO)統(tǒng)計(jì)量=0.732>0.6，巴特利特(Bartlett)球形度檢驗(yàn)F=340.39，自由度(df)=45，p<0.001，表明適合進(jìn)行主成份分析。選擇特征值>1的主成分有2個(gè)，累計(jì)方差貢獻(xiàn)率為72.04%，這2個(gè)主成分解釋了大部分土壤養(yǎng)分信息。第1主成分(PC1)方差貢獻(xiàn)率為44.55%，因子載荷>0.6的有：有機(jī)質(zhì)、速效氮、有效磷、速效鐵、速效鋅、速效硼；第2主成分(PC2)方差貢獻(xiàn)率為27.49%，因子載荷>0.6的有：pH、速效鈣、速效鎂。根據(jù)相關(guān)性和綜合荷載值進(jìn)行指標(biāo)篩選，PC1中有機(jī)質(zhì)的綜合荷載值最大為2.039，10%范圍為1.835 1～2.242 9，其與有效磷的相關(guān)系數(shù)r<0.6，與速效氮、速效鐵、速效鋅和速效硼呈顯著正相關(guān)(r>0.6)，則選擇有機(jī)質(zhì)、速效氮、有效磷納入最佳量化指標(biāo)；PC2中pH的綜合荷載值最大為1.676，10%范圍為1.508～1.843 6，其與速效鈣、速效鎂呈顯著負(fù)相關(guān)(r>0.6)，則選擇pH、速效鎂納入最佳量化指標(biāo)。綜上所述，pH、有機(jī)質(zhì)、速效氮、有效磷、速效鎂共5項(xiàng)指標(biāo)，作為反映土壤肥力的最小數(shù)據(jù)集指標(biāo)。

3.5 土壤綜合肥力指數(shù)及其與種子產(chǎn)量的相關(guān)性

由表7可知，馬尾松單株種子產(chǎn)量與土壤有機(jī)質(zhì)、速效養(yǎng)分質(zhì)量分?jǐn)?shù)均具有顯著正相關(guān)關(guān)系(p<0.05)，而與pH無顯著相關(guān)。

表6 土壤養(yǎng)分指標(biāo)載荷矩陣和綜合荷載值

由表8可知，全部指標(biāo)土壤綜合肥力指數(shù)平均為40.28%，變化幅度為20.93%～70.49%，最小數(shù)據(jù)集指標(biāo)的土壤肥力指數(shù)平均為41.82%，變化幅度為21.62%～71.88%，屬于中等土壤肥力水平。

表7 土壤養(yǎng)分指標(biāo)和單株種子產(chǎn)量的相關(guān)性系數(shù)

表8 全指標(biāo)和最小數(shù)據(jù)集指標(biāo)土壤綜合肥力指數(shù)

由圖2可知，基于10個(gè)全部指標(biāo)的綜合肥力指數(shù)與基于5個(gè)指標(biāo)的綜合肥力指數(shù)間具有極顯著相關(guān)性(R2=0.899 1)。

IT為全部指標(biāo)土壤綜合肥力指數(shù)，IM為最小數(shù)據(jù)集指標(biāo)土壤綜合肥力指數(shù)。

由圖3可知，兩種類型下土壤綜合肥力指數(shù)(IT和IM)與單株種子產(chǎn)量均具有極顯著相關(guān)性(p<0.001)。表明土壤綜合肥力的大小能夠顯著影響種子產(chǎn)量變化，而使用5個(gè)指標(biāo)計(jì)算的綜合肥力指數(shù)能夠代表整個(gè)土壤肥力指數(shù)，減小數(shù)據(jù)的冗余，且擬合效果較好。

◆為IT數(shù)據(jù)，——為IT線性；▲為IM數(shù)據(jù)，----為IM線性。

4 結(jié)論與討論

馬尾松種子園經(jīng)營管理過程中，豐富的土壤養(yǎng)分狀況是保證種子園高產(chǎn)重要因素。本研究中，馬尾松種子園土壤的pH<5.5，其中，樣地66.67%土壤的pH<4.5，處于酸性和強(qiáng)酸性環(huán)境，主要是由于種子園長期施用化肥導(dǎo)致[21-22]，一般認(rèn)為馬尾松適宜土壤pH為5.5～6.5。土壤有機(jī)質(zhì)是土壤微生物代謝活動(dòng)的能量來源，種子園有機(jī)質(zhì)質(zhì)量分?jǐn)?shù)適宜，樣地73.83%的土壤有機(jī)質(zhì)質(zhì)量分?jǐn)?shù)處于中等以上水平(≥15 g·kg-1)，這與曹勝等[23]研究結(jié)果相似。種子園速效氮豐富，樣地95.24%的土壤速效氮質(zhì)量分?jǐn)?shù)處于中等以上水平(≥60 mg·kg-1)，由于氮肥的施用比例偏高，氮肥過量投入是導(dǎo)致土壤酸化的重要原因[24]。土壤有效磷、速效鋅和速效硼缺乏嚴(yán)重，樣地71.73%的土壤有效磷質(zhì)量分?jǐn)?shù)處于低水平(1.0～2.5 mg·kg-1)，樣地16.67%的土壤有效磷質(zhì)量分?jǐn)?shù)處于極缺狀態(tài)(<1.0 mg·kg-1)，樣地61.9%的土壤速效硼質(zhì)量分?jǐn)?shù)處于低水平(0.25～0.5 mg·kg-1)，樣地28.57%的土壤速效硼質(zhì)量分?jǐn)?shù)處于極低水平(<0.25 mg·kg-1)。樣地23.81%和45.24%土壤速效鋅質(zhì)量分?jǐn)?shù)處于低水平(0.5～1.0 mg·kg-1)和極缺狀態(tài)(<0.5 mg·kg-1)。曹繼釗等[19]發(fā)現(xiàn)廣西八角林土壤速效氮質(zhì)量分?jǐn)?shù)處于中上水平，速效磷和速效硼質(zhì)量分?jǐn)?shù)處于極缺狀態(tài)，速效鋅屬中下水平，這樣本研究結(jié)果相一致。土壤速效鉀、速效鈣和速效鎂的豐缺規(guī)律類似均呈現(xiàn)中度缺乏，21.43%土壤處于極缺水平，45%～65%土壤處于低水平，這與杜賢明等[25]研究的結(jié)果相一致。由于土壤酸化嚴(yán)重，加之降水豐富，土壤中可移動(dòng)性K+、Ca2+和Mg2+離子流失造成土壤養(yǎng)分匱乏[26]。南方紅壤是經(jīng)過脫硅富鐵鋁等一系列地球化過程形成，鐵雖然不是葉綠素組成成分，但卻是葉綠素合成的重要物質(zhì)，本研究中，樣地66.67%的速效鐵質(zhì)量分?jǐn)?shù)處于中等以上水平(>20 mg·kg-1)，相對比較適中。

養(yǎng)分障礙因子分析發(fā)現(xiàn)，氮、磷、鎂和鋅是影響種子園產(chǎn)量單一效應(yīng)主導(dǎo)因子，氮-鎂、鉀-鋅和鈣-鎂相互作用效應(yīng)顯著，均對種子產(chǎn)量具有顯著促進(jìn)作用。陳海艷等[20]研究表明，磷、鎂、硼和鉬及氮-硼、鎂-硼、磷-鎂、鎂-鉬的交互作用是影響杉木種子園產(chǎn)量主導(dǎo)因子；陳元鎮(zhèn)等[17]研究表明氮、鎂和鉬及氮-鈣、氮-硼、磷-鉀、磷-鈣、鉀-鎂、鉀-硼的交互作用是影響馬尾松種子園產(chǎn)量主導(dǎo)因子，這與本研究結(jié)果并不完全一致。雖然立地類型以及種子園管理措施不一致，但均呈現(xiàn)出氮、磷和鎂對產(chǎn)量具有較強(qiáng)的作用。此外，本研究選擇pH、有機(jī)質(zhì)、速效氮、有效磷、速效鎂5項(xiàng)指標(biāo)作為反映土壤綜合肥力的最小數(shù)據(jù)集指標(biāo)。土壤pH作為土壤生物化學(xué)環(huán)境反映，顯著影響土壤微生物代謝、土壤速效元素釋放，以及植物養(yǎng)分吸收和果實(shí)產(chǎn)量和質(zhì)量[27]；有機(jī)質(zhì)作為土壤重要組成部分，影響土壤中許多反應(yīng)發(fā)生[28]；土壤氮和磷可用于土壤診斷和肥力評估[29]；鎂是葉綠素的主要成分，與產(chǎn)量密切相關(guān)[30]。因此，pH、有機(jī)質(zhì)、速效氮、有效磷、速效鎂可以作為反映種子園土壤綜合肥力最小數(shù)據(jù)集指標(biāo)。

根據(jù)全部指標(biāo)計(jì)算的土壤綜合肥力指數(shù)平均值為40.28%，最小數(shù)據(jù)集指標(biāo)平均值為41.82%，屬于中等土壤肥力水平，土壤綜合肥力水平不高。兩種方法計(jì)算的土壤綜合肥力指數(shù)之間，以及分別與單株種子產(chǎn)量均具有極顯著相關(guān)性，表明土壤綜合肥力大小能夠顯著影響種子產(chǎn)量的高低；而使用5個(gè)指標(biāo)計(jì)算的綜合肥力指數(shù)能夠代表整個(gè)土壤肥力指數(shù)，減小數(shù)據(jù)的冗余，且擬合效果較好。因此，在種子園的經(jīng)營和管理過程中應(yīng)該注意配方施肥或培肥技術(shù)的應(yīng)用，提高種子園土壤綜合肥力，從而提高種子產(chǎn)量。應(yīng)重視氮、磷、鎂、鋅肥的合理施用，補(bǔ)充磷、鋅、硼肥；采取相應(yīng)措施改善土壤pH，使土壤由強(qiáng)酸性向酸性、弱酸性轉(zhuǎn)變。本研究僅作為種子園土壤養(yǎng)分診斷及土壤綜合肥力的基礎(chǔ)調(diào)查研究，為種子園合理利用土壤養(yǎng)分和精確測土配方施肥提供基礎(chǔ)數(shù)據(jù)。