李成義，強(qiáng)正澤，王燕，王明偉，李碩

(甘肅中醫(yī)學(xué)院藥學(xué)院，甘肅 蘭州 730000)

·中藥農(nóng)業(yè)·

基于因子分析探討影響甘肅紅芪地理分布的環(huán)境因子△

李成義，強(qiáng)正澤，王燕，王明偉，李碩*

(甘肅中醫(yī)學(xué)院藥學(xué)院，甘肅 蘭州 730000)

目的：采用因子分析探討影響甘肅紅芪分布的關(guān)鍵環(huán)境因子。方法：應(yīng)用SPSS 21.0對(duì)紅芪生長(zhǎng)的海拔、經(jīng)度、緯度、年平均氣溫、最熱月平均氣溫、極端最高氣溫、極端最低氣溫、平均無(wú)霜期、≥10 ℃積溫日、≥0 ℃積溫日、≥10 ℃積溫、≥0 ℃積溫、年平均降水量、年均日照數(shù)、年輻射量、年平均干燥度等環(huán)境因子進(jìn)行因子分析。結(jié)果：提取了5個(gè)公因子，分別為熱量因子(F1)”、“地理因子(F2)”、“降水因子(F3)”、“大氣因子(F4)”以及“光照因子(F5)，其中熱量因子是影響紅芪分布的關(guān)鍵性因子。影響紅芪生長(zhǎng)的主要環(huán)境因素為≥10 ℃積溫、緯度、平均降水量、平均無(wú)霜期、年均日照數(shù)。結(jié)論：熱量因子是影響紅芪分布的關(guān)鍵環(huán)境因子，因子得分及排名高的地區(qū)紅芪分布的適應(yīng)性較強(qiáng)。

紅芪；地理分布；因子分析；環(huán)境因子；關(guān)鍵環(huán)境因子

紅芪是甘肅地產(chǎn)大宗藥材，始載于《神農(nóng)本草經(jīng)》黃芪項(xiàng)下，藥用歷史悠久，具有補(bǔ)氣固表，利尿托毒，排膿，斂瘡生肌的功效。近年來(lái)，人們對(duì)自身保健問(wèn)題的高度重視，掀起了“回歸自然”的浪潮，作為藥食同源中藥材紅芪受到人們的親睞，但由于紅芪產(chǎn)區(qū)分布分散，各地區(qū)環(huán)境因子的差異性較大，加上嚴(yán)重的病蟲(chóng)害，導(dǎo)致各地流通的紅芪在品質(zhì)上參差不齊，不僅影響了紅芪產(chǎn)業(yè)的發(fā)展，而且嚴(yán)重阻礙了地方經(jīng)濟(jì)，探尋紅芪的適宜生長(zhǎng)區(qū)域，對(duì)獲得高品質(zhì)紅芪具有重大意義。同時(shí)，隨著科學(xué)技術(shù)的飛躍式發(fā)展，關(guān)于紅芪的質(zhì)量標(biāo)準(zhǔn)、化學(xué)成分、藥理作用、分子特點(diǎn)等方面研究已相對(duì)深入[1-9]，但對(duì)影響紅芪生長(zhǎng)的環(huán)境因子研究較少，不利于紅芪野生資源的合理開(kāi)發(fā)利用和GAP種植基地的建設(shè)。

本文通過(guò)收集31個(gè)產(chǎn)地紅芪生長(zhǎng)的16個(gè)環(huán)境因子，采用因子分析的方法，研究了影響紅芪分布的關(guān)鍵性環(huán)境因子，以期為紅芪的適宜分布區(qū)域研究、資源科學(xué)利用及種植區(qū)域的合理選擇提供參考依據(jù)。

1 材料與方法

1.1 材料

以甘肅地產(chǎn)紅芪為研究對(duì)象，經(jīng)由甘肅中醫(yī)學(xué)院藥學(xué)院李成義教授鑒定，為豆科植物多序巖黃芪HedysarumpolybotrysHand.Mazz.

紅芪分布的環(huán)境因子主要有海拔(E1)、經(jīng)度(E2)、緯度(E3)、年平均氣溫(E4)、最熱月平均氣溫(E5)、極端最高氣溫(E6)、極端最低氣溫(E7)、平均無(wú)霜期(E8)、≥10 ℃積溫日(E9)、≥0 ℃積溫日(E10)、≥10 ℃積溫(E11)、≥0 ℃積溫(E12)、年平均降水量(E13)、年均日照數(shù)(E14)、年輻射量(E15)、年平均干燥度(E16)，當(dāng)?shù)貧庀缶痔峁┘t芪產(chǎn)地氣象數(shù)據(jù)(本文略)。

甘肅紅芪分布信息見(jiàn)表 1。

1.2 方法

由于紅芪的環(huán)境因子變量有16個(gè)，故采用因子分析法，達(dá)到降維的效果，使降維后的主成分能夠代表原來(lái)16個(gè)環(huán)境因子的變異，使分析過(guò)程變得簡(jiǎn)單、準(zhǔn)確。第一對(duì)環(huán)境因子進(jìn)行相關(guān)性分析，判斷各環(huán)境因子之間相關(guān)性程度，第二對(duì)環(huán)境因子數(shù)據(jù)進(jìn)行KMO抽樣適度和Bartlett球形檢驗(yàn)確定是否采用因子分析的方法，第三，在此基礎(chǔ)上進(jìn)行主成分提取和最大方差法旋轉(zhuǎn)，得到主成分的特征值和貢獻(xiàn)率以及各主成分載荷量，第四，以主成分載荷量的大小對(duì)主成分進(jìn)行命名，據(jù)特征值和貢獻(xiàn)率確定影響紅芪分布的關(guān)鍵生態(tài)因子，計(jì)算紅芪分布各產(chǎn)區(qū)的得分和排名，以得分及排名較高的地區(qū)作為紅芪分布的有利區(qū)域。

表1 甘肅紅芪分布信息

2 結(jié)果

2.1 環(huán)境因子相關(guān)性分析

采用Spearman相關(guān)性分析方法，研究紅芪生長(zhǎng)區(qū)域環(huán)境因子之間的相關(guān)性。見(jiàn)表2。由表2可以看出，大部分環(huán)境因子之間存在相關(guān)性，故此可以對(duì)環(huán)境因子進(jìn)行降維分析。KMO抽樣適度檢驗(yàn)值為0.698，大于0.6，Bartlett 球形檢驗(yàn)值為428.32，P<0.05，進(jìn)一步說(shuō)明紅芪生長(zhǎng)區(qū)域環(huán)境因子數(shù)據(jù)可用于因子分析。見(jiàn)表2。

表2 紅芪分布區(qū)域16個(gè)環(huán)境因子的相關(guān)性

2.2 提取公因子

采用主成分提取方法提取能夠代表原始變量的公因子，公因子的數(shù)量由初始特征值和方差貢獻(xiàn)率(%)確定。主成分方法提取的16個(gè)環(huán)境因子的初始特征值和方差貢獻(xiàn)率見(jiàn)表3，成分與特征值的變化情況見(jiàn)圖1。由圖1可知，前5個(gè)成分的變化急劇，累計(jì)方差貢獻(xiàn)率達(dá)到85.90%，包含了初始環(huán)境因子信息的85.90%，僅僅損失了14.10%信息量，可以代表各環(huán)境因子對(duì)紅芪分布區(qū)域的影響情況，故提取前5個(gè)成分作為公因子，分析影響紅芪分布的環(huán)境因子。

表3 公因子的特征值及貢獻(xiàn)率

圖1 16個(gè)環(huán)境因子的碎石圖

2.3 公因子載荷量及解釋公因子

采用主成分分析方法提取公因子后，應(yīng)用具有Kaiser標(biāo)準(zhǔn)化的正交旋轉(zhuǎn)法，旋轉(zhuǎn)在7次迭代后收斂，得到5個(gè)公因子的載荷量，見(jiàn)表4。因子載荷量的值越大，表明因子與變量的相關(guān)性越強(qiáng)，因子對(duì)變量的代表性也越大[10]。據(jù)表4中的數(shù)據(jù)提取每個(gè)因子中因子載荷量較大一個(gè)或幾個(gè)變量，依據(jù)其性質(zhì)給予合理的解釋。

第一公因子F1中，因子載荷量較大的環(huán)境因子有E11、E12，即≥10 ℃積溫、≥0 ℃積溫，這兩個(gè)變量與紅芪生長(zhǎng)區(qū)域的熱量有關(guān)，而且≥10 ℃積溫與≥0 ℃積溫呈正相關(guān)，表明其與F1的變化一致，因此F1可稱之為“熱量因子”。

第二公因子F2中，因子載荷量較大的環(huán)境因子有E1、E3，即海拔、緯度，其反應(yīng)了紅芪生長(zhǎng)區(qū)域的地理環(huán)境，故F2可稱之為“地理因子”。

第三公因子F3中，因子載荷量較大的環(huán)境因子是E13，即年平均降水量，代表了紅芪生長(zhǎng)區(qū)域的降水情況，故F3可稱之為“降水因子”。

第四公因子F4中，因子載荷量較大的環(huán)境因子有E4、E8，即年平均氣溫、平均無(wú)霜期，是紅芪生長(zhǎng)區(qū)域內(nèi)大氣情況的代表性因素，故F4可稱之為“大氣因子”。

第五公因子F5中，因子載荷量較大的環(huán)境因子有E14、E15，即年均日照數(shù)、年輻射量，反應(yīng)的是紅芪生長(zhǎng)區(qū)域的光照情況，故F5可稱之為“光照因子”。

綜上所述，影響紅芪生長(zhǎng)的公因子有“熱量因子(F1)”、“地理因子(F2)”、“降水因子(F3)”、“大氣因子(F4)”、“光照因子(F5)”，5個(gè)公因子代表了影響紅芪生長(zhǎng)的16個(gè)環(huán)境因子的變異。由表3數(shù)據(jù)可以看出，熱量因子方差貢獻(xiàn)率為42.440%，大于其他4個(gè)公因子的方差貢獻(xiàn)率，說(shuō)明了熱量因子對(duì)紅芪分布的影響最顯著，地理因子、降水因子、大氣因子、光照因子對(duì)紅芪對(duì)紅芪的分布影響較熱量因子弱。

表4 公因子載荷量

2.4 確定紅芪生長(zhǎng)適宜區(qū)域

以各紅芪產(chǎn)區(qū)綜合得分及排名為標(biāo)準(zhǔn)，確定紅芪生長(zhǎng)的適宜區(qū)域。紅芪產(chǎn)區(qū)得分的高低，表示紅芪對(duì)生長(zhǎng)區(qū)域的適應(yīng)性的強(qiáng)弱，得分越高，表明此區(qū)域越適合紅芪的生長(zhǎng)。以各公因子的方差貢獻(xiàn)率為權(quán)重系數(shù)，計(jì)算紅芪個(gè)產(chǎn)地的綜合得分，Y=0.424F1+0.173F2+0.132F3+0.073F4+0.057F5，各因子得分計(jì)算公式Fi=∑bijXj(bij為因子得分矩陣第i行、第j列的系數(shù)，文中未列出；Xj為生態(tài)因子值)[11]。由表5可知，總得分及排名較高的紅芪產(chǎn)區(qū)為天水武山、宕昌縣南陽(yáng)鎮(zhèn)、宕昌縣車(chē)?yán)l(xiāng)、隴西縣首陽(yáng)鎮(zhèn)、宕昌縣賈河鄉(xiāng)、岷縣梅川鎮(zhèn)、宕昌、武都區(qū)馬街鄉(xiāng)、宕昌縣城關(guān)鎮(zhèn)馬鞍山、武都區(qū)甘泉鄉(xiāng)、武都安化鎮(zhèn)米倉(cāng)山及武都區(qū)魚(yú)龍鄉(xiāng)，上述區(qū)域的生態(tài)環(huán)境紅芪的生長(zhǎng)適應(yīng)性較強(qiáng)，這與張鵬等[12]學(xué)者對(duì)紅芪栽培地的選擇一致。

表5 紅芪各產(chǎn)地總得分及排名

3 討論

本文收集了甘肅地區(qū)紅芪生長(zhǎng)區(qū)域的16個(gè)環(huán)境因子，通過(guò)因子分析，確定了影響紅芪分布的5個(gè)顯著性環(huán)境因子，分別是“熱量因子(F1)”、“地理因子(F2)”、“降水因子(F3)”、“大氣因子(F4)”以及“光照因子(F5)”，其中熱量因子對(duì)紅芪生長(zhǎng)區(qū)域的貢獻(xiàn)率是最高的，據(jù)此確定“熱量因子”為影響紅芪分布的關(guān)鍵性因子，說(shuō)明紅芪在生長(zhǎng)過(guò)程中對(duì)自然環(huán)境中熱量的變化是最敏感的；地理因子對(duì)紅芪的影響僅次于熱量因子，表明紅芪生長(zhǎng)時(shí)對(duì)地理位置的要求較高。

通過(guò)因子得分及排名，可以對(duì)紅芪分布的各鄉(xiāng)鎮(zhèn)進(jìn)行相互比較。因子得分及排名較高的地區(qū)，綜合了各項(xiàng)生態(tài)因子，相對(duì)而言是最適合紅芪分布的地區(qū)。本研究統(tǒng)計(jì)的31個(gè)紅芪分布區(qū)中，排名及得分較高的區(qū)域有天水武山、宕昌縣南陽(yáng)鎮(zhèn)、宕昌縣車(chē)?yán)l(xiāng)、隴西縣首陽(yáng)鎮(zhèn)、宕昌縣賈河鄉(xiāng)、岷縣梅川鎮(zhèn)、武都區(qū)馬街鄉(xiāng)、宕昌縣城關(guān)鎮(zhèn)馬鞍山、武都區(qū)甘泉鄉(xiāng)、武都安化鎮(zhèn)米倉(cāng)山及武都區(qū)魚(yú)龍鄉(xiāng)，這些地區(qū)應(yīng)該是紅芪分布的最有利區(qū)域。

本研究研究結(jié)果顯示，熱量因子包括≥10 ℃積溫、≥0 ℃積溫兩個(gè)環(huán)境因子，其中≥10 ℃積溫的因子載荷量是最大的，是影響紅芪生長(zhǎng)的最顯著環(huán)境因子，因此，≥10 ℃積溫是甘肅地區(qū)栽培紅芪時(shí)首要考慮的因素；地理因子包括海拔、緯度兩個(gè)環(huán)境變量，緯度的因子載荷量是最高的，是影響紅芪生長(zhǎng)的第二因素；降水因子中只包含了年平均降水量，因此將年平均降水量作為影響紅芪生長(zhǎng)的第三因素；大氣因子包含了年平均氣溫、平均無(wú)霜期，其中平均無(wú)霜期的因子載荷量最大，是影響紅芪生長(zhǎng)的第四因素；光照因子包括年均日照數(shù)、年輻射量，兩環(huán)境因子中年均日照數(shù)的因子載荷量是最大的，因此年均日照數(shù)事影響紅芪生長(zhǎng)的第五因素。據(jù)此，影響紅芪生長(zhǎng)的主要環(huán)境因素為≥10℃積溫、緯度、平均降水量、平均無(wú)霜期、年均日照數(shù)，這些環(huán)境因子的確定有利于掌握影響紅芪分布的關(guān)鍵環(huán)境因子，為紅芪產(chǎn)地適宜性的研究奠定了基礎(chǔ)。同時(shí)因子分析結(jié)果還表明，紅芪分布趨于熱量條件高，緯度較低，具有一定降水量，光照充足的地區(qū)。

[1] 楊扶德，羅文蓉，林麗，等.甘肅紅芪質(zhì)量標(biāo)準(zhǔn)研究[J].中國(guó)藥房，2008，19(18)：1391-1393.

[2] 馬真金.甘肅人工栽培紅芪的質(zhì)量分析[J].甘肅中醫(yī)，2011，24(8)：15-16.

[3] 劉麗華，李加林.黃芪與紅芪的高效液相色譜鑒定[J].時(shí)珍國(guó)醫(yī)國(guó)藥，2010，21(5)：1277-1278.

[4] 胡芳弟，封士蘭，趙健雄.RP-HPLC法測(cè)定紅芪中兩種異黃酮類成分的含量[J].中藥材，2003，26(9)：634-636.

[5] 李冰，封士蘭，劉小花等.HPLC測(cè)定紅芪藥材中紅芪多糖的含量[J].中成藥，2008，30(5)：716-718.

[6] 董玉梅，王小虎，寇煒等.紅芪多糖聯(lián)合 X 線對(duì) HepG -2 細(xì)胞DNA 損傷的影響[J].實(shí)用腫瘤雜志，2012，27(4)：344-349.

[7] 李世剛，張永琦.紅芪多糖誘導(dǎo)人肝癌HEP-G2細(xì)胞凋亡的作用機(jī)制研究[J].中藥材，2009，32(8)：1249-1250.

[8] 鄭海生，金智生，劉 凱，等.紅芪多糖對(duì)2型糖尿病胰島素抵抗大鼠胰島素敏感性影響的研究[J].中華中醫(yī)藥學(xué)刊，2010，28(7)：1516-1518.

[9] 惠和平，封士蘭，胡芳弟，等.紅芪多糖的體外抗氧化活性研究[J].安徽農(nóng)業(yè)科學(xué)，2010，38(8)：4056-4057.

[10] 湯姿.基于因子分析和聚類分析的生態(tài)環(huán)境質(zhì)量評(píng)價(jià)[J].哈爾濱商業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版)，2006，22(6)：121-124.

[11] 盧穎，劉仁權(quán)，王文全，等.應(yīng)用因子分析探討影響甘草地理分布的環(huán)境因素[J].中草藥，2011，42(9)：1822-1827.

[12] 張鵬，潘水站，張杰，等.隴南紅芪無(wú)公害栽培技術(shù)[J].甘肅農(nóng)業(yè)科技，2005，(8)：68-69.

StudyonEnvironmentalFactorswithInfluencingEffectsonGeographicalDistributionofHedysariRadixinGansuProvinceBasedonFactorAnalysis

LIChengyi，QIANGZhengze，WANGYan，WANGMingwei，LIShuo*

(DepartmentofPharmacy，GansuCollegeofTraditionalChineseMedicine，Lanzhou730000，China)

Objective：Application of factor analysis to studying ecological dominant factors with influencing effects on geographical distribution of Hedysari Radix Methods.Methods：The statistics of the regions in Gansu with Hedysari Radix distribution and the geographical environment factors of them，including elevation，longitude and latitude，annual averrage temperature，July average temperature，extreme maximum temperature，extreme minimum temperature，frostless season，≥10 ℃ accumulative temperature days，≥0 ℃ accumulative temperature days，≥10 ℃ accumulative temperature，≥0 ℃ accumulative temperature，precipitation，annual sunshine hours，annual global solar radiation and annual aridity，are analyzed by the factor analysis.And then SPSS 21.0 software is used in the factor analysis.Results：Five common factors，which are heat factor，geographical factor，precipitation factor，air factor and illumination factor，are extracted.The heat factor is the key ecological factor influencing on geographical distribution of Hedysari Radix.The results also show that ≥10 ℃accumulative temperature，latitude，precipitation，frostless season and annual sunshine hours are important environmental factors impacting on the growth of Hedysari Radix.Conclusion：The heat factor is the environment dominant factor influencing on geographical distribution of licorice，the region with high factor total score is more appropriate for the growth and distribution of Hedysari Radix.

Hedysari Radix；geographical distribution；factor analysis；environment factors；environment dominant factor

2014-11-06)

國(guó)家自然科學(xué)基金委員會(huì)資助項(xiàng)目(地區(qū)科學(xué)基金項(xiàng)目)：甘肅地產(chǎn)藥材紅芪道地性與生態(tài)因子的相關(guān)性研究(81360621)

*

李碩，碩士；研究方向：中藥品種與質(zhì)量研究；Tel：(0931)8765385，E-mail：ice_ls@126.com；29068323@qq.com

10.13313/j.issn.1673-4890.2015.5.013