王文舒，武志博，楊衛(wèi)超，劉世英

（1.阿拉善盟林業(yè)草原研究所，內(nèi)蒙古 巴彥浩特 750306;2.阿拉善盟航空護林站，內(nèi)蒙古 巴彥浩特 750306;3.阿拉善盟林工站，內(nèi)蒙古 巴彥浩特 750306）

阿拉善現(xiàn)有天然梭梭荒漠灌叢植被96.67萬hm2，梭梭人工造林面積已達52.38萬hm2，其中：成林面積11.80萬hm2。這149.20萬hm2梭梭林是阿拉善三大生態(tài)屏障之一。國內(nèi)外市場對肉蓯蓉的需求量不斷增加，肉蓯蓉資源因被大量采挖而日益減少，乃至瀕臨滅絕[1]。為保護肉蓯蓉這一珍貴的藥用植物資源，劉銘庭、李天然等在新疆、內(nèi)蒙古等荒漠地區(qū)開展人工接種肉蓯蓉試驗，并取得成功[2,3]。阿拉善盟人工接種肉蓯蓉技術(shù)也不斷提高，截止到2021年底，人工接種面積已達到9.8萬hm2，年產(chǎn)值達到3億元以上。但由于各地區(qū)經(jīng)營分散，組織化程度低，非公有制林業(yè)發(fā)展迅速等原因，使得目前全盟肉蓯蓉產(chǎn)量數(shù)據(jù)不明確，嚴重影響了肉蓯蓉產(chǎn)業(yè)的發(fā)展。國內(nèi)對肉蓯蓉產(chǎn)量的研究主要集中在接種深度、物質(zhì)生產(chǎn)和分配規(guī)律等方面對肉蓯蓉產(chǎn)量的影響，如孫永強等研究表明，接種深度對接種成功率有一定影響[4]；杜友等探討了檉柳管花肉蓯蓉的源庫關(guān)系和非結(jié)構(gòu)性碳水化合物的生產(chǎn)和分配規(guī)律，通過相應的栽培措施合理調(diào)控檉柳和管花肉蓯蓉間的物質(zhì)分配，為管花肉蓯蓉的高產(chǎn)和規(guī)范化栽培提供了理論依據(jù)[5]。但研究未從宏觀角度來評價和衡量某一地區(qū)的肉蓯蓉單穴產(chǎn)量、公頃產(chǎn)量以及總產(chǎn)量，不能宏觀指導肉蓯蓉產(chǎn)業(yè)發(fā)展。因此，亟須利用較為系統(tǒng)完善的統(tǒng)計和分析方法來對肉蓯蓉產(chǎn)量進行測定。

2018—2020年阿拉善盟林業(yè)草原研究所技術(shù)人員從每年春秋兩季在全盟范圍內(nèi)開展肉蓯蓉產(chǎn)量測定工作，結(jié)合不同起源梭梭林分布區(qū)域，對肉蓯蓉的產(chǎn)量分別進行了測定，總結(jié)形成了符合當?shù)貙嶋H、操作性強的產(chǎn)量測定方法。

1 樣地布設方法

依據(jù)全盟天然梭梭林分布區(qū)域和人工梭梭林大面積種植區(qū)域確定了3條調(diào)查樣線，科泊爾（人工）—吉蘭泰（天然）—銀根（天然）為一條樣線，蘇海圖（人工）—塔木素（天然）—板灘井（人工）為一條樣線，古日乃湖（天然）—拐子湖（天然）為一條樣線。每條樣線上的天然林代表地點各設置3個樣地，共15個樣地，每個樣地大小為100 m×100 m（每公頃株數(shù)≥210株）或者200 m×200 m（每公頃株數(shù)＜210株）；人工林代表地點各設置3個樣地，共9個樣地，樣地大小為5行×100株，其中“行”為隨機選取的集中5行，“株”為每行連續(xù)選擇100株梭梭。

2 產(chǎn)量測定方法

最終產(chǎn)量按照天然梭梭林中肉蓯蓉產(chǎn)量和人工梭梭林中肉蓯蓉產(chǎn)量2種方式分別測定。同時，還結(jié)合入戶調(diào)查法確定最終的產(chǎn)量計算模型。入戶調(diào)查法是由當?shù)亓止ふ竟ぷ魅藛T選取肉蓯蓉采挖大戶，在春、秋兩季肉蓯蓉采收結(jié)束后，統(tǒng)計每戶肉蓯蓉采收量，為推測肉蓯蓉總產(chǎn)量提供依據(jù)，使得肉蓯蓉產(chǎn)量測定的準確性得到了進一步提高。

通過對天然梭梭林和人工梭梭林樣地內(nèi)的肉蓯蓉產(chǎn)量進行調(diào)查和研究，分析研究梭梭單株產(chǎn)量和公頃產(chǎn)量，從而通過回歸分析擬合出肉蓯蓉產(chǎn)量計算模型，推算肉蓯蓉總產(chǎn)量，解決阿拉善地區(qū)肉蓯蓉產(chǎn)量的科學估算方法，為資源的合理開發(fā)利用以及沙產(chǎn)業(yè)的發(fā)展提出合理化建議。

3 數(shù)據(jù)分析方法

利用SPSS軟件對數(shù)據(jù)進行統(tǒng)計分析。以種植密度與谷子產(chǎn)量為技術(shù)指標，通過SPSS曲線估計，以決定系數(shù)R2（擬合優(yōu)度）、F值最大為原則選擇最優(yōu)預測模型，同時兼顧運算簡單、便于操作的實用性原則。

4 結(jié)果與分析

根據(jù)收回的調(diào)查表，我所對數(shù)據(jù)進行了整理和分析，分析結(jié)果如下：

4.1 天然梭梭林中肉蓯蓉產(chǎn)量

4.1.1 天然梭梭林樣地內(nèi)數(shù)據(jù)統(tǒng)計結(jié)果。經(jīng)過初步統(tǒng)計，阿左旗郭新軍樣地內(nèi)，平均單株梭梭肉蓯蓉產(chǎn)量和肉蓯蓉平均公頃產(chǎn)量最高，分別為124.54g和59777.55g；阿右旗塔木素阿拉騰圖雅嘎查樣地內(nèi)，平均單株梭梭肉蓯蓉產(chǎn)量和肉蓯蓉平均公頃產(chǎn)量最低，分別為22.37g和2550g。具體見表1。

表1 天然梭梭林樣地內(nèi)肉蓯蓉產(chǎn)量統(tǒng)計表

4.1.2 天然梭梭林樹體結(jié)構(gòu)因子與單株梭梭肉蓯蓉產(chǎn)量的回歸分析。梭梭的樹體結(jié)構(gòu)會影響到肉蓯蓉的產(chǎn)量，所以，合理的樹體結(jié)構(gòu)可以保持肉蓯蓉的豐產(chǎn)穩(wěn)定性。利用SPSS統(tǒng)計軟件試驗分析了天然梭梭林中3個樹體結(jié)構(gòu)因子（冠幅、樹高、地徑）與單株梭梭肉蓯蓉產(chǎn)量的相關(guān)性，結(jié)果見表2~4。

從表2中可以看出，調(diào)整的R2=0.370，因變量的變化有37.0%可由自變量來解釋。樣本數(shù)有158個，所以，只要有大于15.8%的因變量可以由自變量來解釋就說明模型有價值。因此，可以對該模型進行擬合；

表2 模型匯總

從方差分析表（表3）看，F(xiàn)=31.686，sig=0.000，說明回歸方程是顯著的。因此，說明該回歸方程可用，且效果良好。

表3 方差分析表

從系數(shù)表格（表4）中可以擬合出單株梭梭肉蓯蓉產(chǎn)量回歸方程。即：

單株梭梭肉蓯蓉產(chǎn)量（g）=-571.117+5.913×冠幅（cm）-0.537×樹高（cm）+76.712×地徑（cm）。

通過回歸方程可以看出，冠幅和地徑與單株梭梭肉蓯蓉產(chǎn)量呈正相關(guān)關(guān)系，樹高和常量與單株梭梭肉蓯蓉產(chǎn)量呈負相關(guān)關(guān)系；由表4還可以看出，冠幅和地徑因子與單株梭梭肉蓯蓉產(chǎn)量之間有顯著的相關(guān)性，而樹高和常量與單株梭梭肉蓯蓉產(chǎn)量之間的相關(guān)性不顯著。4.1.3 天然梭梭林中肉蓯蓉畝產(chǎn)量回歸分析。利用SPSS統(tǒng)計軟件對天然梭梭林中樣地內(nèi)的梭梭密度和

表4 系數(shù)a

平均畝產(chǎn)量進行回歸分析，從表5中可以看出，調(diào)整的R2=0.919，因變量的變化有91.9%可由自變量來解釋，擬合優(yōu)度很好；

表5 模型匯總

從方差分析表（表6）看，F(xiàn)=1777.890，sig=0.000，說明回歸方程是顯著的。因此，說明該回歸方程可用，且效果良好。

表6 方差分析表

從系數(shù)表格（表7）中可以擬合出畝產(chǎn)量的回歸方程。即：畝產(chǎn)量（g）=-14 362.7+153.504*密度（株/hm2）。

表7 系數(shù)a

由回歸方程可以看出，公頃產(chǎn)量跟密度與常量呈正相關(guān)關(guān)系，且相關(guān)性極顯著。理論產(chǎn)量與實際產(chǎn)量之間的誤差平均為19.89%。

4.1.4 阿拉善盟天然梭梭林中肉蓯蓉總產(chǎn)量。根據(jù)回歸方程，阿拉善盟天然梭梭林密度平均按照150株/hm2，面積按照96.67萬hm2計算，全盟天然梭梭林中的肉蓯蓉總產(chǎn)量=[-14 362.7+153.504×1590（株/hm2）]×966 667 hm2=837 414 150 0 g=8 374.141 5 t?？紤]誤差值后，阿拉善盟天然梭梭林中肉蓯蓉總產(chǎn)量應在6 708.524 8~10 039.7582 t。

4.2 人工梭梭林中肉蓯蓉產(chǎn)量

4.2.1 人工梭梭林樣地內(nèi)數(shù)據(jù)統(tǒng)計結(jié)果。經(jīng)過初步統(tǒng)計（表8），阿拉善右旗板灘井人工梭梭林平均公頃產(chǎn)量較高，最大為5 544.50 g。

表8 人工梭梭林樣地內(nèi)肉蓯蓉產(chǎn)量統(tǒng)計表

4.2.2 人工梭梭林樹體結(jié)構(gòu)因子與單株梭梭肉蓯蓉產(chǎn)量的回歸分析。利用SPSS統(tǒng)計軟件對人工梭梭林中3個樹體結(jié)構(gòu)因子（冠幅、樹高、地徑）與單株梭梭肉蓯蓉產(chǎn)量的相關(guān)性進行線性回歸分析，結(jié)果見表9~11。

從表9中可以看出，調(diào)整的R2=0.569，因變量的變化有56.9%可由自變量來解釋。因此，可以對該模型進行擬合；

表9 模型匯總

從方差分析表（表10）看，F(xiàn)=27.708，sig=0.000，說明回歸方程是顯著的。因此，說明該回歸方程可用，且效果良好。

表10 方差分析表

從系數(shù)表格（表11）中可以擬合出單株梭梭的肉蓯蓉產(chǎn)量回歸方程。即：

表11 系數(shù)a

單株梭梭肉蓯蓉產(chǎn)量（g）=-476.975+5.056×冠幅（cm）-0.288×樹高（cm）+234.995×地徑（cm）。

通過回歸方程可以看出，冠幅和地徑與單株梭梭肉蓯蓉產(chǎn)量呈正相關(guān)關(guān)系，樹高和常量與單株梭梭肉蓯蓉產(chǎn)量呈負相關(guān)關(guān)系；另外，冠幅和地徑因子與單株梭梭肉蓯蓉產(chǎn)量之間有顯著的相關(guān)性，而樹高和常量與單株梭梭肉蓯蓉產(chǎn)量之間的相關(guān)性不顯著。這與天然梭梭林中各因素與單株梭梭肉蓯蓉產(chǎn)量的相關(guān)性一致。

4.2.3 人工梭梭林中肉蓯蓉畝產(chǎn)量回歸分析。利用SPSS統(tǒng)計軟件對人工梭梭林中樣地內(nèi)的梭梭密度和平均畝產(chǎn)量進行回歸分析，從表12中可以看出，調(diào)整的R2=0.525，因變量的變化有52.5%可由自變量來解釋，可以對該模型進行擬合；

表12 模型匯總b

從方差分析表（表13）看，F(xiàn)=71.822，sig=0.000，說明回歸方程是顯著的。因此，說明該回歸方程可用，且效果良好。

表13 方差分析表

從系數(shù)表格（表14）中可以擬合出公頃產(chǎn)量的回歸方程。即：公頃產(chǎn)量（g）=-34 670.1+126.935*密度（株/hm2）。

表14 系數(shù)a

由回歸方程可以看出，公頃產(chǎn)量跟密度呈正相關(guān)關(guān)系，且相關(guān)性極顯著；與常量呈負相關(guān)關(guān)系，相關(guān)性顯著。理論產(chǎn)量與實際產(chǎn)量之間的誤差平均為13.58%。

4.2.4 阿拉善盟人工梭梭林中肉蓯蓉總產(chǎn)量。根據(jù)回歸方程，阿拉善盟人工梭梭林密度平均按照660株/hm2，可接種肉蓯蓉的梭梭林面積按照11.8萬hm2計算，全盟人工梭梭林中的肉蓯蓉總產(chǎn)量=[-34 670.1+126.935×660（株/hm2）]×118 000 hm2=5 794 626 000 g=5 794.626 t?？紤] 誤 差 值 后，阿拉善盟人工梭梭林中肉蓯蓉總產(chǎn)量應在5 007.72~6 581.54 t。

5 結(jié)論與討論

5.1 結(jié)論

全盟肉蓯蓉總產(chǎn)量在11 7 16.2 44 8~16 621.298 2 t，按照1∶7的干肉蓯蓉產(chǎn)出比來計算，全盟干肉蓯蓉年產(chǎn)量在1 673.749 2~2 374.471 2 t之間。

5.2 討論

肉蓯蓉作為一種全寄生植物，其自身不能合成養(yǎng)分來保證生長的需要。梭梭作為其寄主，提供了肉蓯蓉寄生后所有的養(yǎng)分。因此在梭梭-肉蓯蓉復合體內(nèi)，寄主梭梭的長勢對肉蓯蓉的產(chǎn)量影響極為顯著。關(guān)于通過提高肉蓯蓉人工栽培技術(shù)來增加肉蓯蓉產(chǎn)量的研究不斷有新的報道，如郭東峰的研究表明，施肥能夠促進寄主植物的光合作用，從而提高肉蓯蓉產(chǎn)量[6]。楊國濤等研究認為，當土壤含水量為70%~80%時，管花肉蓯蓉產(chǎn)量和質(zhì)量較高[7]。因此，如何通過適當?shù)霓r(nóng)藝措施來提高肉蓯蓉的產(chǎn)量，明確其相關(guān)機制就顯得尤為重要。本試驗雖然達到了預期的試驗目的，初步完成了天然梭梭林和人工梭梭林的產(chǎn)量測定模型，并對影響產(chǎn)量的各因素進行了分析，但在研究中依然存在很多不足之處，如對于人工梭梭林，梭梭樹齡及接種肉蓯蓉時的接種深度還有灌水量等都會對肉蓯蓉產(chǎn)量產(chǎn)生很大影響，但在本次試驗中還未考慮這些因素。