周媛媛，周向睿，周志宇，金茜，李金輝，宋鑫

（1.草地農(nóng)業(yè)生態(tài)系統(tǒng)國家重點實驗室，蘭州大學草地農(nóng)業(yè)科技學院，甘肅 蘭州730020；2.甘肅農(nóng)業(yè)大學草業(yè)學院，甘肅 蘭州730070）

毛烏素沙漠南緣地區(qū)屬于多層次過渡帶，生態(tài)環(huán)境脆弱，是沙丘移動的活躍地帶。在干旱氣候與不合理的人類活動下［1］，成為我國沙漠化較為嚴重的地區(qū)之一。雖然國家早就認識到要對土壤沙漠化進行預(yù)防和治理，并采取行政手段大面積封育、飛播和人工造林等手段，經(jīng)過遙感影像分析發(fā)現(xiàn)：在過去30年里，毛烏素沙漠及周邊地區(qū)的沙漠化得到全面逆轉(zhuǎn)［2］，但2000－2005年期間，僅有2.33%的沙漠化土地發(fā)生了動態(tài)變化。目前，沙塵暴依舊影響著北方地區(qū)，并有南侵的趨勢。因此防止沙漠化已成為一項長期而艱巨的工程。

榆林地區(qū)蒙古巖黃芪在沙地生長良好，可以有效固定流動沙丘，改善沙丘土壤狀況［3－4］，因此對其他地區(qū)沙漠化防治和治理有一定的借鑒和指導(dǎo)意義。蒙古巖黃芪（Hedysarummongolicum），又名踏郎、羊柴、沙欠、山花子等，是一種生長快、耐貧瘠、寒冷、干旱及沙埋等具有很強根蘗能力的多年生豆科灌木［5］。具體來說，蒙古巖黃芪的根狀莖具有極強的營養(yǎng)繁殖能力，使它能在半固定、流動沙地上生存，并表現(xiàn)出極強的生命力［5］。因其幼苗生長速率大于土壤干沙層厚度的變化，故是目前飛播效果最好的樹種之一。此外，它是良好的蜜源，可以改良土壤，又有較高的飼用價值，被喻為沙漠里的“紅花苜蓿”。學者對蒙古巖黃芪的研究主要集中在生物學特性［6］、生長環(huán)境與群落特征［7］、抗逆性［8］、育苗造林、撫育管理及病蟲害防治技術(shù)等方面。

變異是生物界基本特征之一，對形態(tài)變異的度量、表述和分析是研究進化問題的基礎(chǔ)［9］。時間變異和空間變異是變異的兩種類型，目前多數(shù)學者［10－11］研究的是植物形態(tài)隨空間變化而產(chǎn)生的變異，梁坤倫等［12］對青藏高原的3個紫穗槐（Amorphafruticosa）居群的主要形態(tài)特征變異進行了分析，楊東華和趙雨森［4］研究了不同生境對踏郎（Hedysrummongolicum）生長及土壤養(yǎng)分空間差異的影響，研究發(fā)現(xiàn)：由于環(huán)境不同，不同居群內(nèi)及居群間的形態(tài)特征存在一定程度的變異。為了更好、更全面地了解這種西部荒漠地區(qū)適宜發(fā)展的優(yōu)良灌木，本文主要應(yīng)用主成分分析法（principal component analysis，簡稱PCA）［13］從時間軸這個方向來研究蒙古巖黃芪的表型變異特征，通過本文的研究，可以了解蒙古巖黃芪生長的過程中表型性狀變異的趨勢及影響這些變異的主要因子，初步判定了這些變異與土壤環(huán)境之間的關(guān)系，為沙區(qū)種植蒙古巖黃芪和沙漠化可持續(xù)治理提供理論依據(jù)，為生物進化研究奠定理論基礎(chǔ)。

1 材料與方法

1.1 研究區(qū)概況

試驗區(qū)位于107°28′－111°15′E，36°57′－39°34′N，海拔980～1534m，絕對氣溫－27.0～37.6℃，平均氣溫7.9℃，無霜期150d，年平均降水438.4mm（主要集中在7－9月），屬于半干旱內(nèi)陸性季風氣候的陜西省榆林市靖邊縣毛烏素沙漠南緣。地貌主要為流動、半固定或固定沙丘；土壤為各類風沙土；植被類群主要為沙生植物：樟子松（Pinussytvestrisvar.mongolica）、河北楊（Populushopeiensis）、花棒（Hedysarumscoparium）、紫穗槐、檸條（Caraganakorshinskii）、酸刺（acidthorn）、蒙古韭（Alliummongolicum）、苦豆子（Sophoraalopecuroides）、沙柳（Salixpsammophila）、沙蒿（Artemisiaordosica）、沙棗（Elaeagnusangustifolia）、沙打旺（Astragalusadsurgens）、沙蓬（Agriophyllumsquarrosum）等。

1.2 樣品的采集

2012年8月在榆林市靖邊縣毛烏素沙漠南緣，選取1957年（55齡）、1984年（28齡）、2004年（8齡）種植的蒙古巖黃芪各3個樣地（每個樣地面積為100m×100m），在各樣地內(nèi)隨機選取3株中等大小樣株，選擇莖干相對粗壯、冠幅大的植株（母株）作為實驗材料，與由種子掉落而生長出來的或是根狀莖產(chǎn)生的弱小幼齡新株相互區(qū)別，供采樣分析。單株植物測定冠幅、自然高度、葉厚、葉長、葉寬（20片葉子）、主枝（莖干較粗較高）直徑、主枝重量、一級側(cè)枝長度、一級側(cè)枝直徑、一級側(cè)枝重量、一級側(cè)枝總長度、二級側(cè)枝長度、二級側(cè)枝直徑、二級側(cè)枝重量、二級側(cè)枝總長度、主根直徑、根重量等形態(tài)指標。其中葉長、葉寬、葉厚、主枝直徑、一級側(cè)枝直徑、二級側(cè)枝直徑及根（基部）直徑使用游標卡尺測定，枝序的測定采用離心法。

1.3 土樣采集與理化分析

去除落葉層，用土壤刀從植株基部開始逐段、逐層挖去上層覆土，追蹤根系的伸展方向，沿側(cè)根找到須根部分，剪下分枝，然后輕輕抖動，落下的為非根際土，仍粘在根上的為根際土，收集保存，供分析用［14－15］。將野外帶回的樣品放置在干燥通風的室內(nèi)，使自然風干（防止污染）。樣品風干后，揀去雜質(zhì)，研細并全部通過1mm孔徑的篩子，裝、封袋后儲藏備用。

土壤全氮全磷用凱氏消化法測定，硝態(tài)氮和氨態(tài)氮用2mol／L KCl浸提，并用FIA star 5000全自動流動注射儀（瑞典FOSS公司生產(chǎn)）測定；土壤有機碳采用重鉻酸鉀氧化－外加熱法測定；速效磷用0.5mol／L NaHCO3（pH值8.5）鉬銻抗比色法測定；全鉀用微波消解法消化，之后利用火焰分光光度計測定。

1.4 數(shù)據(jù)分析

采用SPSS 17.0和Microsoft Excel軟件進行統(tǒng)計分析，應(yīng)用單因素方差分析來比較不同株齡蒙古巖黃芪各形態(tài)指標的差異性；用主成分分析法（PAC）先對各指標的數(shù)據(jù)進行標準化，即將原始數(shù)據(jù)歸一化處理成均值為0、方差為1的分析數(shù)據(jù)，通過對相似矩陣和協(xié)方差進行分析，提取特征值和特征向量，確定主要因子［16］，并對前兩個主成分進行排序（先算出每個主成分的得分再排序）；最后用Person相關(guān)分析確定前面主成分分析得出的主要因子與根際土壤的全氮、銨態(tài)氮、硝態(tài)氮、全磷、速效磷、全鉀及有機碳的關(guān)系。

表1 不同株齡蒙古巖黃芪17個表型性狀的差異Table 1 The difference about 17indicators of H.mongolicumin different ages

2 結(jié)果與分析

2.1 不同株齡蒙古巖黃芪的表型特征

如表1，不同株齡的蒙古巖黃芪，除了葉長、葉寬、一級側(cè)枝直徑、二級側(cè)枝直徑及根直徑之外，其余指標都是55齡最高，且均是8齡＜28齡＜55齡，其中自然高度，主枝重量，一級側(cè)枝及二級側(cè)枝重量在各年份中的值差異顯著，8，28，55齡的自然高度分別為103.6，150.6和197.7cm，三者的一級側(cè)枝重量分別為11.2，56.4和99.8 g。差異幅度最大的是主枝重量和二級側(cè)枝重量，55齡的（分別為231.6與556.1g）比8齡的（14.2與32.6g）分別增加了15和16倍。其次，最高值比最低值增加倍數(shù)的大小順序是：一級側(cè)枝重量（7.9倍）＞根重（3.0倍）（55齡59.0g，8齡14.7g）＞二級側(cè)枝總長（2.8倍）（55齡100.6m，8齡26.5m）＞二級側(cè)枝直徑（2.2倍）（28齡3.0mm，8齡6.6mm）＞二級側(cè)枝長（1.1倍）（55齡68.9cm，8齡32.4cm）＞一級側(cè)枝總長（2.1倍）（55齡18.8m，8齡9.0m）。剩余指標最大值比最小值增加百分數(shù)的排列順序是：一級側(cè)枝長（97.6%）（55齡89.3 cm，8齡45.2cm）＞自然高度（91.0%）＞冠幅（77.0%）（55齡236.7cm，8齡133.7cm）＞主枝直徑（62.6%）（55齡17.4mm，8齡10.7mm）＞一級側(cè)枝直徑（51.9%）（28齡5.0mm，8齡7.6mm）＞根直徑（17.1%）（8齡13.7mm，55齡11.7mm）＞葉厚（1.8%）（55齡0.56mm，8齡0.55mm）。而比較特殊的是葉長（8，28和55齡分別為26.8，21.1和20.8mm）和葉寬（7.5，7.2和6.9mm），株齡小的反而大（圖中平均值為絕對平均值）。

2.2 不同株齡蒙古巖黃芪形態(tài)變異主成分分析

通過對不同株齡蒙古巖黃芪各指標主成分分析，找出了4個主成分（特征值大于1）（表2），其中前兩個主成分累積貢獻率達78.4%，基本可以概括生長情況的大部分信息，所以只取前兩個主成分，表3列出了各指標對前兩個主成分的載荷系數(shù)，即特征向量。第一主成分的特征值為9.428，貢獻率達55.5%，自然高度（0.946）、根重（0.910）、二級側(cè)枝重量（0.892）、二級側(cè)枝總長（0.891）、冠幅（0.850）、主枝重量（0.845）、一級側(cè)枝長（0.822）這7個指標所占的負荷量絕對值較大。說明這些指標之間有很高的相關(guān)性，存在信息上的重疊，它們代表了蒙古巖黃芪枝、根重量和長度方面的信息。第二主成分占總效益權(quán)重的22.907%（特征值為3.894），特征向量較大的是一級側(cè)枝直徑（0.819）和根直徑（0.810），代表了蒙古巖黃芪根、枝直徑方面的信息。這里選取了負荷量絕對值大于0.8的指標作為主要因子，因此，自然高度、根重、二級側(cè)枝重量、二級側(cè)枝總長、冠幅、主枝重量、一級側(cè)枝長、一級側(cè)枝直徑及根直徑9個指標是造成不同株齡蒙古巖黃芪表型差異的主要因素。

表2 主成分的特征值、貢獻率和累積貢獻率Table 2 The character value，contribution proportion and cumulative contribution proportion

表3 17個性狀對前2個主成分的負荷量Table 3 Loading dose of two principal components for seventeen characters

2.3 不同株齡蒙古巖黃芪前兩個主成分排序

通過對第一、二主成分排序，可以看出在第一主成分上不同株齡的蒙古巖黃芪隨年齡的增加而劃分為3個較為明顯的區(qū)域，其中1，2，3指55齡的3個重復(fù)，4，5，6指28齡的3個重復(fù)，依此，1，2，3是8齡的3個重復(fù)。且55與28齡蒙古巖黃芪居群之間的距離小于28與8齡之間的距離。除第一主成分外，其余各主成分說明的是樣本內(nèi)部的各方面特征，這些重復(fù)在第二主成分上排序較為混亂（圖1）。

2.4 不同株齡蒙古巖黃芪形態(tài)差異的主要因子與根際土壤生境因子的相關(guān)性分析

由表4可以看出，在P＜0.05的水平上時，一級側(cè)枝長，一級側(cè)枝重量及二級側(cè)枝總長與氨態(tài)氮含量呈顯著負相關(guān)（相關(guān)系數(shù)為分別為－0.767，－0.763，－0.721），一級側(cè)枝直徑與速效磷和全鉀含量呈顯著負相關(guān)關(guān)系（相關(guān)系數(shù)分別為－0.794和－0.792），其他指標與全氮、全磷、有機碳、氨態(tài)氮、硝態(tài)氮、速效磷和全鉀有相關(guān)性，但相關(guān)系數(shù)不大，即相關(guān)關(guān)系不顯著。

綜上，由表2可以看出，3個不同株齡中，55齡蒙古巖黃芪長勢最好，方差分析很好地反映了不同株齡蒙古巖黃芪間各形態(tài)指標的變異情況，但由于指標多，為了用盡可能少的維數(shù)去反映樣本之間的聯(lián)系，進行降維，分析結(jié)果顯示蒙古巖黃芪17個表型性狀最后可以歸為4個主要成分，由于前兩個主成分占總效益權(quán)重的78.365%，因此這兩個主要成分就可以反饋原來變量所要反映的主要信息。而自然高度、冠幅等9個指標所占的負荷量大，故是造成這種變異的主要因子。相關(guān)性分析表明氨態(tài)氮，速效磷和全鉀的含量與這些主要因子中的某些因子之間存在顯著的相關(guān)關(guān)系。

圖1 第一、二主成分排序圖Fig.1 Coordinate diagram of principal components

表4 各主要因子與根際土壤養(yǎng)分的相關(guān)性系數(shù)矩陣Table 4 The correlation coefficient matrix between every main factors and the nutrient of the rhizosphere soil

3 討論

種內(nèi)變異是個體和居群在時空變化發(fā)展中的產(chǎn)物，對變異規(guī)律的深入研究，使我們可以全面準確地了解種內(nèi)個體的形態(tài)差異，對揭示物種的形成與進化機制提供幫助［16］。通過本文研究發(fā)現(xiàn)，隨著株齡的增長，蒙古巖黃芪的部分表型形狀發(fā)生了不同程度的變異，呈不均衡性，表明這些性狀對環(huán)境適應(yīng)的敏感程度不同。

選擇合適的數(shù)學模型對所得數(shù)據(jù)進行分析，有助于揭示變異規(guī)律以及對環(huán)境的適應(yīng)能力［17］。本文先計算出不同株齡蒙古巖黃芪隨時間變異的總趨向及差異性。發(fā)現(xiàn)55齡的蒙古巖黃芪長勢優(yōu)于其他株齡，除葉長、葉寬、一級側(cè)枝直徑、二級側(cè)枝直徑、根直徑之外，其他指標都隨年齡增加呈增大趨勢，并且自然高度和一級側(cè)枝重量這兩個指標在各株齡間差異顯著。主枝重量和二級側(cè)枝重量變異最大。這些現(xiàn)象說明蒙古巖黃芪在干旱沙化地區(qū)生長的過程中，隨著株齡的增加，地上和地下指標為適應(yīng)環(huán)境變化向最優(yōu)、最適宜生存的方向發(fā)展，即在干旱沙地的惡劣氣候環(huán)境下，生長和繁衍所采取的向空間擴展來提高對空間資源利用的生存策略，使蒙古巖黃芪更容易獲得水分、光照和養(yǎng)分等［18－19］。葉片是植物的重要組成部分，葉長、葉寬隨株齡減小而增大，這與主成分分析出來的結(jié)果（載荷量為－0.694，－0.585的負向作用）一致，然而葉厚呈遞增趨勢，引起這一現(xiàn)象的原因可能是由于株齡大的蒙古巖黃芪面臨較多的逆境而形成的適應(yīng)缺水，強光，溫差大等環(huán)境氣候來減少水分損失和保持養(yǎng)分的機制［20－21］。

主成分最大的優(yōu)點就是方差貢獻率大，其中各性狀載荷值的大小體現(xiàn)了各性狀在主成分中的重要程度，于是自然高度、根重、二級側(cè)枝重量、二級側(cè)枝總長、冠幅、主枝重量、一級側(cè)枝長、一級側(cè)枝直徑及根直徑這9個指標所占負荷量較大，故是造成表型差異的主要因素，也就是蒙古巖黃芪的枝與根的重量和直徑、枝的長度能反映出其生長的大體情況。因此在育種時注意這些方面可以來評價蒙古巖黃芪種的優(yōu)劣性。這些主要因素的變異越大說明植株適應(yīng)環(huán)境的能力越強，豐富的變異有助于保持物種和生態(tài)系統(tǒng)的多樣性和穩(wěn)定性［22－23］。趙俊卉等［24］對長白山云杉、冷杉和紅杉的研究發(fā)現(xiàn)，競爭是引起枝長和直徑變異的主要因素。

對各重復(fù)前兩個主成分進行排序，在第一數(shù)軸上，不同株齡的重復(fù)在時間增加的基礎(chǔ)上進行分類，聚為3個不同年份的團體，與實際情況相吻合，即個體形態(tài)特征在一定程度上發(fā)生了分化［25］。這個結(jié)果與第一主成分綜合原始指標信息的能力最強，可以全面反映各指標狀況的綜合指標，能說明原始數(shù)據(jù)變動的總趨勢的功能相呼應(yīng)［26］，很好地反映了蒙古巖黃芪隨年齡增加表型性狀在可遺傳的、有序的、穩(wěn)定的變異軌跡中的變化，符合進化論的基本規(guī)律，表明蒙古巖黃芪植株在形態(tài)學功能上對時間變異的適應(yīng)性，同時也印證了關(guān)系越密切的植株排序越近的理論［27］。類似的實驗結(jié)果有：通過環(huán)境因子對青南地區(qū)植物群落進行排序發(fā)現(xiàn)植物群落在空間地理上分布很有規(guī)律［28］。

影響生物表型變異的因素主要有兩個，即遺傳因素和環(huán)境因素。遺傳因子是受基因控制對植物表型變異影響較為緩慢、穩(wěn)定，而環(huán)境因子能夠在短期內(nèi)影響和改變植物表型變異，二者對種內(nèi)形態(tài)變異相輔相成，共同控制著表型變異的趨向［29－30］。影響生物表型的環(huán)境因子［16］主要包括：土壤環(huán)境、物候環(huán)境，種的生態(tài)分布區(qū)等。本研究從土壤環(huán)境這方面進行研究，對9個主要因子與之對應(yīng)的根際土壤全氮等7個養(yǎng)分因子做了相關(guān)性分析，一級側(cè)枝長度，二級側(cè)枝總長及重量與氨態(tài)氮之間存在顯著的負相關(guān)性（即氨態(tài)氮與第一主成分的三個主要因子的相關(guān)關(guān)系顯著，因此可以把第一主成分的主要因子看作是氨態(tài)氮主控因子），速效磷和全鉀與一級側(cè)枝直徑之間也呈顯著的負相關(guān)性。表明蒙古巖黃芪的這9個主要因子與土壤養(yǎng)分存在一定的聯(lián)系，這是因為植物在生長的過程中受環(huán)境影響，同時也影響著土壤環(huán)境［31－33］，而作為影響蒙古巖黃芪的主要因子，與土壤養(yǎng)分之間必然有密切的關(guān)系。土壤是蒙古巖黃芪營養(yǎng)元素的主要源泉，而蒙古巖黃芪根系、落葉物殘體為土壤提供腐殖質(zhì)，并隨著株齡的增加對土壤改良具有顯著的效果［34］。氨態(tài)氮之所以與第一主成分的三個主因子成負相關(guān)關(guān)系，可能是因為在蒙古巖黃芪自身生長過程中吸收了大量氨態(tài)氮以使側(cè)枝快速生長。表型變異是變異的重要線索，本文從植物進化年齡，即株齡變化這個方向?qū)γ晒艓r黃芪表型性狀進行了研究，但這種變異具體是因遺傳因素占的份額大還是環(huán)境因素大，還有待于從分子方面進行進一步研究。

Reference：

［1］Hao C Y，Wu S H，Yang Q Y.Research on relationship between sandy desertification and land－use in Mu Us Region.Journal of Desert Reseach，2005，25（1）：33－39.

［2］Guo J，Wang T，Han B S，etal.Dynamic process of aeolian desertification land variation in Mu Us sandy land and its surrounding area in recent 30years.Journal of Desert Reseach，2008，28（6）：1017－1021.

［3］Zhang C S，Liu G B，Xue S，etal.Anti－erodibility and microaggregate fractal features of soils under different psammophytes.Bulletin of Soil and Water Conservation，2012，32（2）：1－6.

［4］Yang D H，Zhao Y S.Effects of site condition on growth and spatial variations ofHedysarumlaeve.Scientia Silvae Sinicae，2011，47（2）：180－186.

［5］Lv Y L.The preliminary observation on the shape and anatomical structure on rhizome ofHedysarummongolicum.Journal of Beijing Forestry University，1986，（1）：91－95.

［6］Gu S.The relationship between the growth characteristics and climatic factors ofHedysarumscoparium，HedysarumleaveandKorshinsk peashrubat seedling stage.Arid Zone Research，1994，11（2）：60－63.

［7］Jiang G M.Species and habitat variablity of photosynthesis，transpiration and water use efficiency of different plant speicies in Maowusu sand area.Acta Botanica Sinica，1999，41（10）：1114－1124.

［8］Guo J P，Gao S H，Wang L M，etal.Response ofHedysarummongolicumturcz.to high CO2concentration and soil drought stress.Journal of Plant Resources and Environment，2002，11（1）：14－16.

［9］Ge S，Hong D Y.Biosystematic studies onAdenophorapotaniniikorsh.Complex（campanulaceae）.Acta Phytotaxonomica Sinica，1995，33（5）：433－443.

［10］Zhang S P，Wang R Q，Zhang Z G，etal.Study on morphological variation ofPhragmitesaustralisin the Yellow River downstream wetland.Acta Phytoecologica Sinica，2003，27（1）：78－85.

［11］Zhao X L，Ling Y，Zhang G F，etal.Community characteristics of beach wetland vegetations along a habitat gradient in Dafeng Milu reserve of Jiangsu province.Chinese Journal of Ecology，2010，29（2）：244－249.

［12］Liang K L，Jiang W Q，Zhou Z Y，etal.Variation in main morphological characteristics ofAmorphafruticosaplants in the Qinghai－Tibet Plateau.Acta Ecologica Sinica，2012，32（1）：311－318.

［13］Yu S X.Introduction to mathematical ecology.Beijing：Science and Technology Publishing House，1995：207－213.

［14］Li W H.The research on the difference of nitrogen content and pH between rhizosphere and bulk soil ofQuercusvariabilis，Cunninghamia lanceolataandPinustaeda.Journal of Nanjing Forestry University，1996，20（2）：49－52.

［15］Zeng S C，Su Z Y，Chen B G，etal.A review on the rhizospherenutritionecology research.Journal of Nanjing Forestry University（Nature Science Edition），2003，27（6）：79－83.

［16］Yang J.Infraspecific variation in plant and the exploring methods.Journal of Wuhan Botanical Research，1991，9（2）：185－195.

［17］Sneath P，Sokal R.Mumerical Taxonomy［M］.Beijing：Science Press，1984.

［18］He M Z，Zhang J G，Wang H，etal.Analysis of branching architecture factors of desert plants.Journal of Desert Reseach，2006，26（4）：625－630.

［19］Zhao Z J.The Study of Branch and Root Morphometries ofMedicugofalcate［D］.Changchun：Northeast Normal University，2007.

［20］Lambrecht S C，Dawson T E.Orrelated variation of floral and leaf traits along a moisture availability gradient.Oecologia，2007，151（4）：574－583.

［21］Santiago L S，Kitajima K，Wright S J，etal.Coordinated changes in photosynthesis，water relations and leaf nutritional traits of canopy trees along aprecipitation gradient in lowland tropical forest.Oecologia，2004，139（4）：495－502.

［22］Yu H，Kiang Y T.Genetic variation in South Korean natural populations of wild soybean（Glycinesoja）.Euphytica，1993，68（3）：213－221.

［23］Jiang X L，Yue J，Zhang W G，etal.Biodiversity，ecosystem functioning and spatio－temporal scales.Acta Prataculturae Sinica，2010，19（1）：219－225.

［24］Zhao J H，Kang X G，Zhang H D，etal.Relationships between coefficient of variation of diameter and height and competition index of main coniferous trees in Changbai Mountains.Chinese Journal of Applied Ecology，2009，30（8）：1832－1837.

［25］Ma H J，He P，Deng H P，etal.Analyses of morphological differentiation within six populations ofMalushupehensis（Rosaceae）.Journal of Southwest China Nomal University（Natural Science），2002，27（4）：559－562.

［26］Zhang P.The Principal Component Analysis Based Research in Comprehensive Assessment［D］.Nanjing：Institutes of Technology of Nanjing，2004.

［27］Luan Q S，Sun J，Song S，etal.Canonical correspondence analysis of summer phytoplankton community and its environment in the Yangtze river estuary，China.Journal of Plant Ecology（Chinese Version），2007，31（3）：445－450.

［28］Wang W Y，Wang C Y，Wang Q J，etal.The gradient analysis of plant communities at sounthern Qinghai of Qinhai－Tibet Plateau.Journal of Lanzhou University（Natural Science），2004，40（3）：60－65.

［29］Soltis P S，Soltis D E.Genetic variation in endemic and wide spread plant species：examples from Saxifragaceae and Polystichum（Dryopteridaceae）.Aliso，1991，13：215－223.

［30］Ge S.Genetic Diversity［M］.Hangzhou：Zhejiang Science and Technology Press，1997：11－19.

［31］Huang Y P，Chen Y，Guo F X，etal.Effect of N，P，K combinations fertilizion on the cultivation quality ofGlycyrrhizauralensisseedlings.Acta Prataculturae Sinica，2012，21（2）：233－240.

［32］Zhu Q，Pan Y Z，Zhao L.etal.The effects of N，P，K，and Ca on plant growth and nutrient content of Lily leaves.Acta Prataculturae Sinica，2012，21（5）：274－284.

［33］Zhang L，Wang C T，Liu W，etal.Relationships of dominant species root activity，plant community characteristics and soil micro－environment in artificial grassland over different cultivation periods.Acta Prataculturae Sinica，2012，21（5）：185－194.

［34］Wang Y W，Liao C Y，Sun C Z，etal.Soil chemical properties of sand－fixing forests in Maowusu Sandland.Chinese Journal of Soil Science，2009，40（4）：776－780.

［1］郝成元，吳紹洪，楊勤業(yè).毛烏素地區(qū)沙漠化與土地利用研究.中國沙漠，2005，25（1）：33－39.

［2］郭堅，王濤，韓邦帥，等.近30a來毛烏素沙地及其周邊地區(qū)沙漠化動態(tài)變化過程研究.中國沙漠，2008，28（6）：1017－1021.

［3］張昌勝，劉國彬，薛萐，等.不同沙生植被土壤微團聚體分形特征及抗蝕性.水土保持通報，2012，32（2）：1－6.

［4］楊東華，趙雨森.不同生境對踏郎生長及土壤養(yǎng)分空間差異的影響.林業(yè)科學，2011，47（2）：180－186.

［5］呂悅來.蒙古巖黃茵根狀莖的形態(tài)及解剖構(gòu)造的初步觀察.北京林業(yè)大學學報，1986，（1）：91－95.

［6］古松.花棒、楊柴、毛條苗期生長特性與氣候因子的關(guān)系.干旱區(qū)研究，1994，11（2）：60－63.

［7］蔣高明.毛烏素沙地若干植物光合作用、蒸騰作用和水分利用效率種間及生境間差異.植物學報，1999，41（10）：1114－1124.

［8］郭建平，高素華，王連敏，等.楊柴對高CO2濃度和土壤干旱脅迫的響應(yīng).植物資源和環(huán)境學報，2002，11（1）：14－16.

［9］葛頌，洪德元.泡沙參復(fù)合體（桔?？疲┑奈锓N生物學研究Ⅲ.性狀的遺傳變異及其分類價值.植物分類學報，1995，33（5）：433－443.

［10］張淑萍，王仁卿，張治國，等.黃河下游濕地蘆葦形態(tài)變異研究.植物生態(tài)學報，2003，27（1）：78－85.

［11］趙小雷，凌云，張光富，等.大豐麋鹿保護區(qū)不同生境梯度下灘涂濕地植被的群落特征.生態(tài)學雜志，2010，29（2）：244－249.

［12］梁坤倫，姜文清，周志宇，等.青藏高原紫穗槐主要形態(tài)特征變異分析.生態(tài)學報，2012，32（1）：311－318.

［13］余世孝.數(shù)學生態(tài)學導(dǎo)論［M］.北京：科學技術(shù)文獻出版社，1995：207－213.

［14］厲婉華.栓皮櫟、杉木和火炬松根際與非根際土壤氮素及pH差異的研究.南京林業(yè)大學學報，1996，20（2）：49－52.

［15］曾曙才，蘇志堯，陳北光，等.植物根際營養(yǎng)研究進展.南京林業(yè)大學學報（自然科學版），2003，27（6）：79－83.

［16］楊繼.植物種內(nèi)形態(tài)變異的機制及其研究方法.武漢植物學研究，1991，9（2）：185－195.

［17］史尼斯P，索卡爾R.數(shù)值分類學［M］.趙鐵橋，譯.北京：科學出版社，1984.

［18］何明珠，張景光，王輝.荒漠植物枝系構(gòu)型影響因素分析.中國沙漠，2006，26（4）：625－630.

［19］趙占軍.黃花苜蓿分枝與根系形態(tài)的研究［D］.長春：東北師范大學，2007.

［23］江小雷，岳靜，張衛(wèi)國，等.生物多樣性，生態(tài)系統(tǒng)功能與時空尺度.草業(yè)學報，2010，19（1）：219－225.

［24］趙俊卉，亢新剛，張慧東，等.長白山主要針葉樹種胸徑和樹高變異系數(shù)與競爭因子的關(guān)系.應(yīng)用生態(tài)學報，2009，30（8）：1832－1837.

［25］馬洪菊，何平，鄧洪平，等.湖北海棠不同居群形態(tài)變異的數(shù)量分析.西南師范大學學報（自然科學版），2002，27（4）：559－562.

［26］張鵬.基于主成分分析的綜合評價研究［D］.南京：南京理工大學，2004.

［27］欒青杉，孫軍，宋書，等.長江口夏季浮游植物群落與環(huán)境因子的典范對應(yīng)分析.植物生態(tài)學報，2007，31（3）：445－450.

［28］王文穎，王春燕，王啟基，等.青藏高原青南地區(qū)植物群落生態(tài)梯度分析.蘭州大學學報（自然科學版），2004，40（3）：60－65.

［30］葛頌.遺傳多樣性［M］.杭州：浙江科學技術(shù)出版社，1997：11－19.

［31］黃亞萍，陳垣，郭鳳霞，等.氮磷鉀配施對甘草育苗質(zhì)量的影響.草業(yè)學報，2012，21（2）：233－240.

［32］朱嶠，潘遠智，趙莉.氮、磷、鉀、鈣對香水百合生長及葉片養(yǎng)分含量的影響.草業(yè)學報，2012，21（5）：274－284.

［33］張莉，王長庭，劉偉，等.不同建植期人工草地優(yōu)勢種植物根系活力、群落特征及其土壤環(huán)境的關(guān)系.草業(yè)學報，2012，21（5）：185－194.

［34］王彥武，廖超英，孫長忠，等.毛烏素沙地固沙林土壤化學性質(zhì)差異.土壤通報，2009，40（4）：776－780.