基于主成分和聚類分析的德宏橡膠林土壤肥力評(píng)價(jià)

陳玉芹 胡永亮 張麗萍 尹紅星 黃家衛(wèi) 李守嶺

摘 ?要 ?基于德宏橡膠林的78份土壤樣本數(shù)據(jù)，采用主成分分析法和聚類分析法對(duì)橡膠林土壤肥力進(jìn)行綜合評(píng)價(jià)。結(jié)果表明：德宏橡膠林土壤pH均值為4.83;土壤有機(jī)質(zhì)、全鉀、全氮和堿解氮養(yǎng)分較豐富，變異系數(shù)在16.53%～46.90%之間，呈中等變異性;但全磷和有效磷含量普遍低，處于缺磷狀態(tài)。土壤肥力綜合評(píng)價(jià)結(jié)果：橡膠林土壤IFI值范圍10.81～25.29，均值為17.24，屬于中等以上肥力的土壤占到83.34%。采樣區(qū)域肥力水平從高到低依次為盈江農(nóng)場>瑞麗農(nóng)場>畹町農(nóng)場>遮放農(nóng)場。綜上，德宏橡膠林土壤肥力總體水平處于良好，應(yīng)重視測土配方施肥，重點(diǎn)增施磷肥。

關(guān)鍵詞 ?德宏;橡膠;土壤肥力;主成分分析;聚類分析

中圖分類號(hào) ?S794.1; S158.2 ?????文獻(xiàn)標(biāo)識(shí)碼 ?A

Abstract ?The soil fertility of rubber plantation by principal component and cluster analysis based on 78 soil samples of rubber plantation in Dehong was evaluated. The results showed that the pH value was 4.83， the content of soil organic matter， total potassium， total nitrogen and alkaline nitrogen was sufficient and moderate variability， the CV was between 16.53% to 46.90%. But the content of total phosphorus and available phosphorus was very low. Comprehensive assessment of soil fertility revealed that the integrated fertility index IFI was 17.24， ranged in 10.8125.29. There were 83.34% soils belonged to moderate fertility. The soil fertility level was in the order downward Yingjiang Farm， Ruili Farm， Wanding Farm， Zhefang Farm. The soil fertility of rubber plantation is moderate and phosphate should be applied in Dehong.

Keywords ?Dehong; rubber; soil fertility; principal component analysis; cluster analysis

DOI ?10.3969/j.issn.1000-2561.2019.08.001

德宏1904年開始種植橡膠樹，是云南邊境重要的天然橡膠種植地區(qū)之一，種植面積從2010年的1.6萬hm2[1]下降到2018年的9840 hm2[2]。多年來，由于橡膠樹經(jīng)濟(jì)壽命長，隨著膠乳的流出，橡膠樹均要從土壤中吸收大量的養(yǎng)分，而膠農(nóng)對(duì)橡膠林地土壤的養(yǎng)分狀況缺乏認(rèn)識(shí)，對(duì)科學(xué)施肥的意識(shí)及其技術(shù)也不夠，普遍存在重化肥輕有機(jī)肥，重氮肥輕磷肥、鉀肥，以致施肥不足和管理不善，導(dǎo)致土壤肥力下降、橡膠產(chǎn)量和品質(zhì)難以提升。

肥力是土壤的基本屬性和本質(zhì)特征，是衡量土壤能夠提供作物生長所需的各種養(yǎng)分的能力，是土壤物理、化學(xué)和生物學(xué)性質(zhì)的綜合反應(yīng)。而土壤養(yǎng)分是體現(xiàn)土壤肥力的重要部分，在生產(chǎn)上常用有機(jī)質(zhì)、氮、磷、鉀等養(yǎng)分含量來衡量土壤肥力的狀況，土壤酸堿度是土壤肥力的重要指標(biāo)之一，影響土壤肥力形成和質(zhì)量的演變[3];有機(jī)質(zhì)能增強(qiáng)土壤的系水系肥能力，促進(jìn)根系對(duì)礦物質(zhì)營養(yǎng)的吸收[4];堿解氮體現(xiàn)氮素的供應(yīng)情況。以往土壤養(yǎng)分評(píng)價(jià)側(cè)重定性評(píng)價(jià)和單因素評(píng)價(jià)，評(píng)價(jià)標(biāo)準(zhǔn)不一，存在人為的主觀性;隨著數(shù)值化評(píng)價(jià)方法的提出，模糊數(shù)學(xué)方法和多元統(tǒng)計(jì)分析方法等現(xiàn)代研究方法得到廣泛應(yīng)用，土壤養(yǎng)分評(píng)價(jià)逐漸側(cè)重定量評(píng)價(jià)和多因素的綜合評(píng)價(jià)[5]。近年來，植膠區(qū)土壤養(yǎng)分肥力的研究一直是橡膠研究的熱點(diǎn)[6-9]，越來越多的學(xué)者應(yīng)用指數(shù)和法、主成分分析法、聚類分析法和模糊綜合評(píng)價(jià)法等進(jìn)行土壤肥力綜合評(píng)價(jià)。戴余波等[10]學(xué)者采用主成分分析法評(píng)價(jià)了熱帶作物耕地土壤養(yǎng)分肥力情況分為4個(gè)等級(jí)，提出了需要適當(dāng)使用復(fù)合肥、有機(jī)肥來提高土壤的肥力水平;劉少春等[11]學(xué)者用主成分綜合分析將土壤肥力劃分為5個(gè)等級(jí)，評(píng)價(jià)出蔗田土壤肥力總體水平偏低;李世瑤等[12]研究表明主成分分析具有較好的代表性與客觀性，可用于畦灌灌水質(zhì)量的評(píng)價(jià);安云娜等[13]采用主成分分析將土壤肥力分為5個(gè)等級(jí)，整體水平較差。本研究旨在了解德宏橡膠林地土壤養(yǎng)分狀況并基于主成分和聚類分析對(duì)其作出評(píng)價(jià)，為德宏橡膠林地土壤的改良、培肥提供基礎(chǔ)數(shù)據(jù)資料，對(duì)實(shí)現(xiàn)德宏橡膠林地土壤可持續(xù)利用具有重要意義。

1 ?材料與方法

1.1 ?研究區(qū)域概況

德宏州在云南省西部中緬邊境，位于東經(jīng)97°31′～98°43′、北緯23°50′～25°20′之間，年降雨量在1400～1700 mm，年平均氣溫在18.4～20.0?℃，年日照2281～2453 h，年積溫6400～7300?℃，年陸地蒸發(fā)量在1400～1900 mm，干旱指數(shù)在0.4～1.2。屬于南亞熱帶半濕潤季風(fēng)氣候，地形呈山、丘、壩分布，土壤分為6個(gè)土類、11個(gè)亞類、19個(gè)土層、34個(gè)土種，1904年開始種植橡膠樹，主要分布在芒市、瑞麗、盈江3個(gè)縣（市）區(qū)，涉及5個(gè)國營農(nóng)場和18個(gè)鄉(xiāng)鎮(zhèn)[14]，橡膠林土壤主要有黃壤、磚紅壤、紅壤3類。

1.2 ?方法

1.2.1 ?土壤樣品采集 ?以德宏盈江農(nóng)場、遮放農(nóng)場、瑞麗農(nóng)場、畹町農(nóng)場4個(gè)農(nóng)場為采樣區(qū)域，采用GPS定位，在每個(gè)區(qū)域選取有代表性的橡膠林地采用混合采樣法取20 cm以下耕層土壤樣品，每塊橡膠林地根據(jù)面積大小采5～8個(gè)點(diǎn)，混合成一個(gè)樣，用“四分法”取約1 kg鮮樣，經(jīng)室內(nèi)風(fēng)干后剔除異物，過60目和100目尼龍篩，備用（土壤樣品采集參照NY/T 1121.1-2006）。

1.2.2 ?分析項(xiàng)目與方法 ?土壤pH、有機(jī)質(zhì)、全磷、全氮、堿解氮、有效磷的測定參照《土壤農(nóng)業(yè)化學(xué)分析方法》[15]，全鉀采用原子吸收分光光度法測定，測定時(shí)加入空白樣品及標(biāo)準(zhǔn)土壤樣品進(jìn)行分析質(zhì)量控制。

采用因子分析中提取主成分分析法（PCA）和系統(tǒng)聚類分析，對(duì)橡膠土壤養(yǎng)分進(jìn)行綜合評(píng)價(jià)。PCA法是利用降維的思想，把多指標(biāo)轉(zhuǎn)化為少數(shù)幾個(gè)因子作分析的一種多元統(tǒng)計(jì)方法[5]。采用SPSS軟件通過因子分析[16]確定參評(píng)土壤指標(biāo)主成分特征值和特征向量，根據(jù)特征值>1選為關(guān)鍵主成分，計(jì)算各主成分得分，再利用綜合得分公式求出各樣點(diǎn)土壤養(yǎng)分綜合分值（IFI）。再根據(jù)土壤養(yǎng)分綜合分值，采用類平均法對(duì)土壤樣品進(jìn)行系統(tǒng)聚類分析[17]。

1.3 ?數(shù)據(jù)處理

采用SPSS 13.0統(tǒng)計(jì)軟件對(duì)土壤養(yǎng)分進(jìn)行描述性統(tǒng)計(jì)和K-S檢驗(yàn)，檢驗(yàn)土壤數(shù)據(jù)是否服從正態(tài)分布，對(duì)不服從正態(tài)分布的數(shù)據(jù)，剔除異常值數(shù)據(jù)，直至服從正態(tài)分布。

2 ?結(jié)果與分析

2.1 ?德宏橡膠園土壤指標(biāo)描述性統(tǒng)計(jì)

德宏橡膠林地土壤統(tǒng)計(jì)分析結(jié)果顯示（表1，圖1），采樣區(qū)域的樣點(diǎn)土壤pH分布頻數(shù)較集中介于4.4～5.0，平均pH為4.83呈酸性土壤，變異系數(shù)僅有6.32%，屬于弱變異。由表1可知，表層土中有機(jī)質(zhì)分布頻數(shù)較散，主要含量集中在20～40 g/kg，平均含量為30.52 g/kg;全磷含量分布頻數(shù)集中在0.40～0.60 g/kg，平均含量為 ???0.50 g/kg;全鉀含量分布頻數(shù)集中在10～25 g/kg，平均含量17.69 g/kg;全氮含量分布頻數(shù)集中在0.8～1.8 g/kg，平均含量1.42 g/kg;堿解氮含量分布頻數(shù)集中分布在120～170 mg/kg，平均含量145.27 mg/kg;通過土壤樣品的變異系數(shù)可以得到，有機(jī)質(zhì)、全磷、全鉀、全氮和堿解氮等養(yǎng)分變異系數(shù)為16.53%～46.9%呈中等變異性，表明在采樣區(qū)域內(nèi)有一定差異，但變化不大。由圖1可知，能被作物吸收的有效磷含量分布不集中，較散且平均含量為1.35 mg/kg，變異系數(shù)達(dá)到81.55%（表1），呈高變異性，說明不同土樣中有效磷含量差異巨大，可能受人類活動(dòng)和土壤酸堿性影響較大。另外，除堿解氮，其有機(jī)質(zhì)、全氮、全鉀、有效磷等6個(gè)土壤養(yǎng)分指標(biāo)的中位數(shù)均低于其均值，從K-S檢驗(yàn)和分布類型來看（表1，圖1），有機(jī)質(zhì)、全鉀、全氮服從正態(tài)分布;全磷、堿解氮、有效磷和pH剔除異常值后服從正態(tài)分布，其中，全磷和堿解氮偏度為負(fù)值，呈左偏態(tài)分布，有機(jī)質(zhì)、全鉀、全氮、有效磷和pH偏度為正值，呈右偏態(tài)分布。

2.2 ?橡膠林土壤養(yǎng)分的分級(jí)比例

由表2可知，德宏橡膠林土壤中有機(jī)質(zhì)、全鉀、全氮、堿解氮等養(yǎng)分含量較高。約有69.3%的樣點(diǎn)含有25 g/kg的有機(jī)質(zhì)，僅有12.8%的樣點(diǎn)有機(jī)質(zhì)含量低;84.6%的樣點(diǎn)全鉀含量達(dá)到正常水平，52.5%樣點(diǎn)全鉀含量處于高水平;全氮有93.4%超過0.8 g/kg，僅有6.6%樣點(diǎn)處于低水平;堿解氮有78.2%的樣點(diǎn)超過120 mg/kg，呈豐富水平。說明采樣區(qū)域樣點(diǎn)土壤具有較好的保水保肥能力和相對(duì)較高的氮鉀肥力。但全磷和有效磷含量普遍偏低，其中全磷低于0.6 g/kg的樣點(diǎn)占總樣本的71.4%，有效磷含量低于5 mg/kg的樣點(diǎn)占總樣本的100%，說明采樣區(qū)域的橡膠林土壤中全磷和有效磷極度缺乏。另外，從pH來看，85.1%土樣pH介于4.5～5.5，屬于酸性土壤，有9.5%土樣呈強(qiáng)酸性，5.4%土樣呈弱酸性，根據(jù)pH變異系數(shù)較穩(wěn)定來看（表1），德宏橡膠林采樣區(qū)域主要呈酸性土壤。

2.3 ?土壤養(yǎng)分綜合評(píng)價(jià)

2.3.1 ?指標(biāo)間相關(guān)性檢驗(yàn)、KMO和Bartlett球形度檢驗(yàn) ?因子分析的必要條件是各指標(biāo)間存在相關(guān)性，而SPSS軟件中的主成分分析是采用因子分析中的主成分提取，即需對(duì)各土壤養(yǎng)分指標(biāo)進(jìn)行相關(guān)性檢驗(yàn)，將土壤全鉀、有機(jī)質(zhì)、堿解氮、有效磷、全磷、全氮等養(yǎng)分的原始數(shù)據(jù)采用描述性統(tǒng)計(jì)得到標(biāo)準(zhǔn)化數(shù)據(jù)，計(jì)算其相關(guān)系數(shù)，從相關(guān)系數(shù)矩陣（表3）可以看出，有機(jī)質(zhì)、堿解氮、全氮等指標(biāo)間存在不同程度的顯著相關(guān)性，說明有部分信息重疊。經(jīng)KMO和Bartlett球形度檢驗(yàn)，KMO值為0.550，Bartlett球形度檢驗(yàn)的相伴概率

P<0.05（顯著水平）進(jìn)一步說明本研究中土壤養(yǎng)分各指標(biāo)間存在相關(guān)性。因此，本研究采用PCA分析法評(píng)價(jià)土壤肥力較合適。

2.3.2 ?土壤養(yǎng)分主成分分析 ?由表4中的共同度數(shù)據(jù)可以看出，主成分包含了各個(gè)原始變量至少65%以上的信息，其中全氮指標(biāo)的信息量損失最大，有機(jī)質(zhì)指標(biāo)的信息量損失最少。根據(jù)提取特征值>1的為主成分，可提取3個(gè)主成分并分別解釋了34.891%、19.248%、18.326%的方差，累積貢獻(xiàn)率達(dá)到72.464%（表5），這表明3個(gè)主成分集中了原始土壤養(yǎng)分6個(gè)變量信息的72.464%。根據(jù)因子載荷系數(shù)與主成分系數(shù)之間的關(guān)系（）[18]，可看出土壤養(yǎng)分指標(biāo)在主成分上的權(quán)重比例不同，主成分1主要受有機(jī)質(zhì)、堿解氮、全氮的權(quán)重影響，主成分2主要反映全鉀和有效磷的權(quán)重，主成分3主要體現(xiàn)全磷的權(quán)重影響（表6）。

根據(jù)主成分得分系數(shù)矩陣（表6）可得出主成分的綜合得分線性表達(dá)式：

2.3.3 ?德宏橡膠林土壤養(yǎng)分綜合評(píng)價(jià) ?按綜合評(píng)價(jià)模型計(jì)算出不同采樣區(qū)域橡膠地土壤養(yǎng)分綜合得分進(jìn)行統(tǒng)計(jì)描述。從表7可知，德宏橡膠林土壤綜合肥力總體得分范圍10.81～25.29，平均分值17.24，變異系數(shù)為18.94。從各采樣區(qū)域來看，盈江農(nóng)場、畹町農(nóng)場、瑞麗農(nóng)場和遮放農(nóng)場的橡膠林土壤綜合肥力得分范圍分別為11.51～25.29，14.73～18.59，10.81～22.13，11.85～ 22.06。從高到

低依次為盈江農(nóng)場>瑞麗農(nóng)場>畹町農(nóng)場>遮放農(nóng)場;另外，4個(gè)農(nóng)場采樣區(qū)域土壤肥力變異系數(shù)來看畹町農(nóng)場處于低變異性，盈江農(nóng)場、瑞麗農(nóng)場、遮放農(nóng)場3個(gè)農(nóng)場處于中等變異范圍，說明各個(gè)區(qū)域內(nèi)部土壤肥力狀況差異不盡相同。盈江農(nóng)場土壤肥力水平處于中等水平以上，土壤營養(yǎng)狀況較好;畹町農(nóng)場和瑞麗農(nóng)場土壤肥力處于中等水平;遮放農(nóng)場土壤肥力主要處于缺乏到中等水平，需重點(diǎn)開展測土配方施肥。

采用類平均法[19]對(duì)土壤綜合分值（IFI）進(jìn)行系統(tǒng)聚類將78個(gè)樣本分為3類，即第Ⅰ類IFI≥20.79，第Ⅱ類14.19≤IFI≤20.13，第Ⅲ類10.81≤ IFI≤13.68，分別代表土壤肥力高、中、低3個(gè)等級(jí)（表8）。德宏橡膠林土壤采樣點(diǎn)有15.38%和67.95%的土壤肥力為Ⅰ和Ⅱ級(jí)，說明土壤養(yǎng)分總體水平良好，但有16.67%的土壤肥力處于低水平，肥力較差。從采樣區(qū)域肥力土壤所占比例來看，盈江農(nóng)場有55.88%土壤肥力為Ⅱ級(jí)水平，僅有14.71%土壤肥力為Ⅲ級(jí)水平;畹町農(nóng)場整體土壤肥力為Ⅱ級(jí)水平;瑞麗農(nóng)場85.71%土壤肥力為Ⅱ級(jí)水平;遮放農(nóng)場土壤肥力主要為Ⅱ級(jí)和Ⅲ級(jí)水平。

3 ?討論

土壤肥力對(duì)橡膠樹養(yǎng)分吸收、生長代謝以及乳膠的形成等方面具有重要作用。從統(tǒng)計(jì)描述來看，土壤養(yǎng)分頻數(shù)分布不均勻且具有不同程度的變異性，這可能是因?yàn)橄鹉z產(chǎn)量和品質(zhì)下降，膠農(nóng)沒有見到實(shí)質(zhì)性的經(jīng)濟(jì)效益，小部分橡膠林地長期處于無人管理狀態(tài)，橡膠林雜草叢生，大部分橡膠林套種了咖啡、茶葉等經(jīng)濟(jì)作物得到較好的管理。本研究可看出，德宏橡膠林土壤pH頻數(shù)分布集中在4.4～5.0，呈酸性土壤，少部分土壤呈強(qiáng)酸性和弱酸性，且變異性小，說明土壤的酸度受長期的耕作管理的影響不大;有機(jī)質(zhì)、全磷、全鉀、全氮和堿解氮等養(yǎng)分呈中等變異性與黎小清等[9]學(xué)者研究東風(fēng)農(nóng)場橡膠園土壤養(yǎng)分呈現(xiàn)的變異性有相似規(guī)律。土壤中有機(jī)質(zhì)、全鉀、全氮、堿解氮等養(yǎng)分含量較豐富，但此次研究土壤中全磷有71.4%的樣點(diǎn)和有效磷全部樣點(diǎn)含量極低，處于缺磷狀態(tài)，驗(yàn)證了何鵬等[20]研究中提到的植膠區(qū)土壤普遍缺磷的觀點(diǎn)。這可能與磷特有的化學(xué)性質(zhì)，容易被土壤固定而降低有效性以及傳統(tǒng)的施肥方式重氮肥輕磷肥鉀的施肥習(xí)慣，造成土壤中有效磷含量大幅降低。另外，有效磷是反應(yīng)土壤磷素供應(yīng)水平的重要指標(biāo)，橡膠林土壤中能被作物吸收的有效磷含量頻數(shù)分布不集中，且呈高變異性，說明不同土樣中有效磷含量差異巨大，低于橡膠林土壤養(yǎng)分正常5.0～8.0 mg/kg之間[1]，這可能是由于土壤的酸堿度能通過影響礦質(zhì)鹽分的溶解度影響?zhàn)B分的有效性，在土壤酸性條件下，鐵、鋁活性相對(duì)較高，容易與磷發(fā)生作用形成難以溶解的鐵磷和鋁磷，使土壤及人工施入的磷轉(zhuǎn)化為固定態(tài)磷，這是有效磷含量較低且變異系數(shù)大的原因之一。

選取土壤全鉀、有機(jī)質(zhì)、堿解氮、有效磷、全磷、全氮等6個(gè)指標(biāo)進(jìn)行主成分分析并計(jì)算IFI進(jìn)行系統(tǒng)聚類分析。結(jié)果顯示，可提取3個(gè)主成分且集中體現(xiàn)72.464%原始信息，可見采用主成分分析效果是較好的。綜合區(qū)域評(píng)價(jià)來看，橡膠林土壤綜合肥力得分從高到低依次為盈江農(nóng)場>瑞麗農(nóng)場>畹町農(nóng)場>遮放農(nóng)場，其中，盈江農(nóng)場橡膠林地得到較好的管理并充分利用林地空間與其他作物套種，土壤肥力狀況處于中等偏豐富;其次是瑞麗農(nóng)場和畹町農(nóng)場土壤肥力狀況處于中等水平，遮放農(nóng)場土壤肥力狀況中等水平和缺乏水平比例各占到46.15%，要加強(qiáng)水肥管理。通過本研究可以看出，土壤肥力水平分布不均勻，可能受各區(qū)域施肥管理以及農(nóng)業(yè)活動(dòng)等行為有關(guān)。土壤肥力受人為管理影響較大，肥力水平差的應(yīng)加強(qiáng)科學(xué)施肥和測土配方施肥，改善土壤的肥力。

本研究利用主成分分析和聚類分析獲得的土壤肥力綜合評(píng)價(jià)可以較好地反映實(shí)際土壤質(zhì)量，與張子龍等[21]研究主成分分析在三七連作土壤質(zhì)量綜合評(píng)價(jià)中的應(yīng)用有相似的觀點(diǎn)，但是否能推廣應(yīng)用于不同土壤類型還需進(jìn)一步驗(yàn)證。因?yàn)橥寥蕾|(zhì)量受到生態(tài)系統(tǒng)的類型、土壤的功能、土壤利用方式、不同作物等的影響。此外，采用主成分分析和聚類分析方法評(píng)價(jià)土壤肥力，選定指標(biāo)是土壤肥力評(píng)價(jià)的核心，直接關(guān)系到評(píng)價(jià)結(jié)果的客觀性，實(shí)際研究中部分指標(biāo)獲取難度較大，實(shí)用性差，影響了肥力評(píng)價(jià)結(jié)果在實(shí)際中的應(yīng)用推廣。本研究選取的6個(gè)評(píng)價(jià)指標(biāo)較為有限，需進(jìn)一步增加土壤質(zhì)地、微生物、微量元素等指標(biāo)，更加全面的探討橡膠園土壤肥力情況。

4 ?結(jié)論

（1）德宏橡膠林土壤pH頻數(shù)分布集中在4.4～5.0，呈酸性土壤，少部分土壤呈強(qiáng)酸性和弱酸性，且變異性小;有機(jī)質(zhì)、全磷、全鉀、全氮和堿解氮等養(yǎng)分變異系數(shù)在16.53%～46.90%之間，呈中等變異性;有效磷含量頻數(shù)分布不集中，較散且平均含量為1.35 mg/kg，變異系數(shù)達(dá)到81.55%，呈高變異性。

（2）根據(jù)橡膠樹健康生長所需土壤養(yǎng)分含量來看，目前德宏橡膠林的土壤除了全磷和有效磷處于低水平，其他養(yǎng)分均在正常范圍及以上且具有較高的氮鉀肥力。

（3）利用PCA分析法計(jì)算土壤養(yǎng)分綜合分值（IFI），分值范圍10.81～25.29，平均分值17.24，標(biāo)準(zhǔn)差3.27，變異系數(shù)為18.94，屬于中等變異程度。綜合肥力得分從高到低依次為盈江農(nóng)場>瑞麗農(nóng)場>畹町農(nóng)場>遮放農(nóng)場。

（4）利用系統(tǒng)聚類分析法，德宏橡膠林土壤采樣點(diǎn)有83.34%的土壤肥力介于Ⅰ級(jí)～Ⅱ級(jí)水平之間，說明土壤肥力總體水平良好。

參考文獻(xiàn)

王應(yīng)清， 李守嶺， 張麗萍， 等. 云南德宏橡膠園土壤養(yǎng)分分析[J]. 湖南農(nóng)業(yè)科學(xué)， 2015（11）： 55-59.

劉小瓊， 李守嶺， 張麗萍， 等. 膠園間種咖啡對(duì)土壤養(yǎng)分的影響[J]. 熱帶農(nóng)業(yè)科學(xué)， 2018， 38（1）： 1-5.

焉 ?莉， 王 ?寅， 馮國忠， 等. 吉林省農(nóng)田土壤肥力現(xiàn)狀及變化特征[J]. 中國農(nóng)業(yè)科學(xué)， 2015， 48（23）： 4800-4810.

張小琴， 陳 ?娟， 高秀兵， 等. 貴州重點(diǎn)茶區(qū)茶園土壤pH值和主要養(yǎng)分分析[J]. 西南農(nóng)業(yè)學(xué)報(bào)， 2015（1）： 286-291.

呂蘇丹， 汪光宇， 鄔亞浪， 等. 東陽萬畝園區(qū)土壤養(yǎng)分綜合評(píng)價(jià)研究[J]. 浙江大學(xué)學(xué)報(bào)， 2002， 28（3）： 272-276.

杜中杰， 林 ?電， 許 ?杰， 等. 海南省橡膠園土壤養(yǎng)分狀況研究[J]. 廣東農(nóng)業(yè)科學(xué)， 2011， 38（11）： 73-77.

丘秀靈， 劉正富， 阮云澤， 等. 海南橡膠土壤肥力評(píng)價(jià)及影響因素分析[J]. 廣東農(nóng)業(yè)科學(xué)， 2009， 38（9）： 81-85.

楊 ?熙， 林 ?電， 余 ?愛， 等. 海南省主要橡膠園土壤養(yǎng)分限制因子[J]. 熱帶作物學(xué)報(bào)， 2009， 30（12）： 1746-1751.

黎小清， 余凌翔， 李春麗， 等. 云南東風(fēng)農(nóng)場橡膠園土壤養(yǎng)分空間分布特征[J]. 西南農(nóng)業(yè)學(xué)報(bào)， 2015， 28（1）： 292-307.

戴余波， 張麗萍， 李國明， 等. 熱帶作物耕地土壤養(yǎng)分分析及肥力評(píng)價(jià)[J]. 現(xiàn)代農(nóng)業(yè)科技， 2017（18）： 155-157.

劉少春， 張躍彬， 郭家文， 等. 基于養(yǎng)分豐缺分級(jí)的蔗田土壤肥力主成分綜合分析[J]. 西南農(nóng)業(yè)學(xué)報(bào)， 2016， 29（3）： 611-617.

李世瑤， 蔡煥杰， 陳新明. 基于主成分分析的畦灌質(zhì)量評(píng)價(jià)[J]. 農(nóng)業(yè)工程學(xué)報(bào)， 2013， 29（24）： 86-93.

安云娜， 黃義雄， 官紫玲. 福建東山島土壤養(yǎng)分綜合評(píng)價(jià)[J]. 安徽農(nóng)業(yè)科學(xué)， 2007， 35（13）： 3926-3927.

釧相仙， 白燕冰， 李金濤， 等. 云南德宏、臨滄天然橡膠產(chǎn)業(yè)現(xiàn)狀與發(fā)展意見[J]. 熱帶農(nóng)業(yè)科技， 2015， 35（2）： 9-15.

魯如坤. 土壤農(nóng)業(yè)化學(xué)分析方法[M]. 北京： 中國農(nóng)業(yè)出版社， 2000.

黃 ?安， 楊聯(lián)安， 杜 ?挺， 等. 基于主成分分析的土壤養(yǎng)分綜合評(píng)價(jià)[J]. 干旱區(qū)研究， 2014， 31（5）： 819-825.

趙華富， 周國蘭， 劉曉霞， 等. 貴州茶區(qū)土壤養(yǎng)分狀況綜合評(píng)價(jià)[J]. 中國土壤與肥料， 2012（3）： 30-34.

李小勝， 陳珍珍. 如何正確應(yīng)用spss軟件做主成分分析[J]. 統(tǒng)計(jì)研究， 2010， 27（8）： 105-108.

宋勤飛， 牛素貞， 陳正武， 等. 基于主成分分析的花溪古茶樹立地土壤養(yǎng)分評(píng)價(jià)[J]. 浙江農(nóng)業(yè)學(xué)報(bào)， 2017， 29（11）： 1844-1853.

何 ?鵬， 吳 ?敏， 韋家少， 等. 不同磷水平對(duì)橡膠樹磷吸收利用特性的影響[J]. 中國農(nóng)學(xué)通報(bào)， 2011， 27（19）： 14-18.

張子龍， 王文全， 廖作清， 等. 主成分分析在三七連作土壤質(zhì)量綜合評(píng)價(jià)中的應(yīng)用[J]. 生態(tài)學(xué)雜志， 2013， 32（6）： 1636-1644.