紅花玉蘭潛在引種氣候適生區研究

施曉燈， 馬履一， 段 劼， 桑子陽， 朱仲龍， 姚 娜， 周明明， 賈忠奎

（1．北京林業大學 省部共建森林培育與保護教育部重點實驗室,北京100083；2.北京林業大學 國家能源非糧生物質原料研發中心，北京100083；3.五峰博翎紅花玉蘭科技發展有限公司，湖北宜昌443400）

木蘭科Magnoliaceae植物在全世界有近11屬246種，主要分布于亞洲東南部溫帶及熱帶、美洲中部等地［1－2］，木蘭屬Magnolia是木蘭科中最大的屬，中國是世界上木蘭屬植物最多的國家，被譽為 “木蘭王國”。從中國西南部省區到東南部省區均有木蘭屬樹種分布［3－4］。紅花玉蘭Magnolia wufengensis及其變種多瓣紅花玉蘭Magnolia wufengensisvar.multitepala是在湖北五峰發現的木蘭科木蘭屬玉蘭亞屬Yulania Subgenus新種［5－6］。其樹形高大，冠型優美，花被片數目、顏色、形狀以及花型等存在豐富的變異。花被片數9～46，花色從深紅到粉紅具有7個色系［7］，花被片形狀有卵形、倒卵形、狹倒卵形和長披針形等，花型有荷花型、月季型、菊花型和牡丹型等［7］。紅花玉蘭可作為城市綠化和山地景觀造林樹種，具有極高的觀賞價值和巨大的開發利用潛力［8］。模糊數學是研究和揭示模糊現象的定量處理方法，廣泛應用于模糊聚類分析、模糊模式識別、模糊綜合評判、模糊線性規劃等方面，在決策、土木、農業、林業、氣象、風險評估等領域［9－15］已有相關研究成果。模糊相似優先比法是模糊數學中的重要內容，它是以成對的樣品與一個固定樣品同時作比較，以確定哪一個樣品與固定樣品最相似，在所有樹木引種區劃的數學方法中是較為合理的且應用較多的方法［16－17］。根據定義優先比的方法不同分為海明距離模糊相似優先比法和相對歐氏距離模糊相似優先比法［16］等。陸鼎煌等［18］曾采用海明距離（絕對值距離）模糊相似優先比法分析了油橄欖Olea europaea在中國的引種區劃，區劃結果與油橄欖實際調查引種現狀幾乎吻合。但袁嘉祖等［16］同時在2個生態環境相似性分析中指出海明距離模糊相似優先比法的不足之處在于分析的數據太多，計算繁瑣。而相對歐氏距離模糊相似優先比法計算起來相對簡便，結果與前者基本一致，如果再考慮各個因素的權重，劃分出的相似性結果會更可靠［19］。有祥亮等［20］利用相對歐氏距離模糊相似優先比法對中國引種栓皮槭Acer campestre進行了生態適生區劃分，并分別提出了相應的栽培措施。同時還有日本甜柿Diospyros kaki［21］， 美國黃松Pinus ponderosa［22］， 梭梭Haloxylon ammodendron［23］等均采用此方法取得了較好的適生區劃分結果。為保證紅花玉蘭能夠成功引種到原生地以外區域，本研究通過應用模糊數學原理［9］，采用相對歐氏距離模糊相似優先比法分析和劃分紅花玉蘭引種氣候適生區，為相應植物引種提供科學依據。

1 材料與方法

1.1 氣候因子選擇

查閱影響植物生長發育的主要氣候因子［24－27］，并依據紅花玉蘭生長特性對氣候的需求［28］，選取年平均氣溫T1（℃），最熱月平均氣溫T2（℃）， 最冷月平均氣溫T3（℃），極端最高氣溫T4（℃），極端最低氣溫T5（℃）， 年平均降水量T6（mm）， 年平均風速T7（m·s－1）， 年平均相對濕度T8（%）， 年日照時數T9（h）等 9個氣候因子作為模糊相似優先比分析的因素。

1.2 資料收集與整理

按照全面、詳細、有代表性的原則，選取全國32個省（市、自治區）196個氣象觀測站點進行模糊相似優先比排序分析。紅花玉蘭是多年生樹種，生長發育受到長期氣候影響，因此，選擇的各個氣象觀測站點的資料為1981－2010年30 a的平均值（表1）。固定樣地的氣象數據采用紅花玉蘭原生地湖北省五峰土家族自治縣1981-2010年相關氣象數據均值。氣象數據大部分來自中國氣象局氣象數據中心（http：//data.cma.cn/site/index.html），個別站點數據來源于各省（市、自治區）氣象網站。

表1 全國196個臺站的1981-2010年主要氣候要素值Table 1 Values of the main climatic elements of 196 stations of China （the average value from 1981 to 2010）

1.3 數據處理與分析方法

試驗采用相對歐氏距離模糊相似優先比法，用Excel 2007和SPSS 22.0軟件對各地氣候數據進行分析，再依據分析的結果，劃分紅花玉蘭引種氣候適生區。

x={x1，x2，…， xn，x0}，n=1， 2，…，196;T={T1，T2，…，Tm}， m=1，2， …，9;xi={xi1，xi2，…，xim}， i=1， 2， …， n;xj={xj1， xj2， …， xjm}， j=1， 2， …， n。 其中： x1， x2， …， xn表示全國 196個氣象觀測站樣本；T1，T2，…，Tm為9個氣候因子；xi，xj為其中任意2個樣本；x0表示湖北五峰氣象觀測站，即固定樣本。

先對中國氣象局氣象數據中心下載的數據進行標準化處理，去除量綱對數據的影響，采用極差正規化方法。公式如下：

式（1）中：ximax，ximin分別是第m個因子中的極大值和極小值，x′im∈［0， 1］。

根據層次分析法［29］對因子之間相對重要性標度（表2）進行打分,然后用求平均及四舍五入的方法構造判斷矩陣（表3），求出各因素的權重值Ws。用歸一后的各氣候因子值x′im數據與相對應的權重 Ws相乘， 得到 x″im。 即： x″im=x′im×Ws。

用相對歐氏距離計算被選樣品 xi與固定樣品x0之間的差距（Di0）。

表2 相對重要性標度Table 2 Means of relative importance scale

表3 紅花玉蘭氣候因子判斷矩陣Table 3 Estimation matrix of climatic elements of Magnolia wufengensis

計算相似優先比（rij）：

式（2）中：rji=1－rij，當i=j時，rij=1。另外，令rij=rji=1，其中，i，j為備選樣品，k為固定樣品。

取λ水平截集，計算出全國196個站點與固定樣本的相似順序。即從表3矩陣中的元素，求出每行的最小值組成一列，將此列中的最大值所在的行所對應的樣本記為1，即該樣本與固定樣本最相似，然后刪除該行及相對應的列，重復上述過程，再依次求取相似樣品，序號為2，3，4，…，196等。

2 結果與分析

2.1 確定氣候因子權重

依據表3可得出各個氣候因子權重值（表4），并對表3進行一致性檢驗得出其最大特征根λmax=9.729 4，一致性指標IC=0.091 18，平均隨機一致性指標IR=1.46，隨機一致性比率RC=0.062 45＜0.1，表明該矩陣具有令人滿意的一致性，得出的權重值不需要調整。

表4 各個氣候因子權重值Table 4 Weights of each climatic factors

2.2 建立模糊相似優先比矩陣

通過式（3）求得模糊相似優先比矩陣R，R=［rij］n×n（表 5）。 其中： 元素r（xi，xj）表示xi，xj與固定樣地x0比較時，xi比xj或者xj比xi的優越程度，若rij∈［0.5，1.0］，則表示xi比xj優先； 若rij∈［0， 0.5］， 則表示xj比xi優先。

2.3 劃分適宜區

根據表5紅花玉蘭模糊相似優先比矩陣采用λ水平截集的計算結果，將各臺站按序號值從小到大進行排列（表6），表中序號值的大小基本反映了全國196個氣象觀測站綜合氣候生態條件與紅花玉蘭原產地五峰相似程度的大小：序號越小，表示相似程度越高；序號越大，表示相似程度越低。根據引種的氣候相似性原理，將中國引種的紅花玉蘭的氣候適生區大致劃分如表6所示，詳細氣候因子劃分范圍見表7。

2.3.1 適宜區 該區序號值最靠前，且相對歐氏距離值均小于等于0.02，與紅花玉蘭原產地五峰氣候條件最相似（表8），主要包括：①華東、華中地區：從地域上看，該區包含地點分別屬于黃淮平原、長江中下游平原、江南丘陵、浙閩丘陵西北部。該區主要以亞熱帶溫濕季風氣候為主，地處濕潤區，與原產地湖北五峰有非常相似的水熱條件，四季分明，降水豐富，冬季氣溫較溫和，夏季氣溫偏高，雨熱同期，降水量多集中在夏季，容易形成高溫高濕的環境，造成圃地內澇，發生病蟲害，其中紅花玉蘭主要病害是葉片表面蓋滿煤污狀黑霉的煤污病和從根莖部開始腐爛的根腐病，主要蟲害是刺吸性害蟲，如紅蜘蛛Tetranychus cinnbarinus，斑潛蠅Liriomyza，粉蚧Pseudococcus等。因此在引種栽植紅花玉蘭時應選擇適宜的立地條件，若在低洼地段應做高床，夏季及時做好圃地清理及排水工作。發生煤污病需噴質量分數為0.3%～0.5%波爾多液進行防治，根腐病可采用500倍多菌靈噴灑，連續打藥2～3次，基本上可控制病情發展，刺吸類害蟲用氧化樂果500～1 000倍液噴霧。②華南、西南地區：從地域上看，該區包含地點分別屬于四川盆地東南部、云貴高原東北部、浙閩丘陵東南部和兩廣丘陵北部。主要以亞熱帶季風氣候為主，地處濕潤區，氣候溫和，濕度較大，例如四川萬源，云南昆明、大理等地四季如春，但是需要注意一些地區降雨季節分配不均，夏秋季引種容易遭受較長期的干旱，紅花玉蘭較難成活，因此這些地區引種栽植紅花玉蘭時，遇到較大旱情需要及時增加灌溉等抗旱措施。

表5 紅花玉蘭模糊相似優先比矩陣Table 5 Fuzzy similitude priority ratio matrix of Magnolia wufengensis

表6 各站點引種紅花玉蘭適應性排序Table 6 Order of stations of China according to the suitability of ordinal numbers

表6 （續）Table 6 Continued

表7 各類適生區氣候因子對比Table 7 Comparison of climate factor types for each suitability county （city）

表8 紅花玉蘭引種適宜區Table 8 Suitable area for the introduction of Magnolia wufengensis

2.3.2 次適宜區 該區序號值排在適宜區之后，且歐氏距離為0.02～0.04，主要是適宜區向四周擴散的區域，綜合氣候條件與紅花玉蘭原產地五峰有一定差異（表9）。主要包括：①華南地區、西南區。從地域上看，該區包含地點分別屬于四川盆地西北部、云貴高原西南部，兩廣丘陵以南，雷州半島，海南島和臺灣島。該區以亞熱帶季風氣候為主，部分高海拔地區為山地氣候，地跨濕潤區、半濕潤區。與適宜區相比，西南山區海拔較高，氣候變化幅度較大，引種紅花玉蘭需要栽植在海拔較低的區域，同時在個別干熱河谷，區域內光熱資源豐富，氣候炎熱少雨，水土流失嚴重，生態十分脆弱，寒、旱、風、蟲、草、火等自然災害特別突出，需要注意這些地區不適宜引種紅花玉蘭。華南地區夏季相對適宜區更加高溫高濕，太陽輻射更強，雨熱同期時間更長，臺風、洪澇災害頻發。這些地區引種紅花玉蘭需要注意地形地勢的選擇，遠離洪澇頻發區域，發生病蟲害應及時進行治理，以防蔓延，出現內澇及時疏通清理，減少引種紅花玉蘭的死亡。②華東、華北、東北地區。從地域上看，該區包含地點分別屬于晉冀魯山地丘陵地帶，華北平原以及遼東半島，以暖溫帶大陸性季風氣候為主，地處半濕潤區。該區光照充足，冬季寒冷干燥，風速較大，夏季高溫炎熱，年平均降水量偏少，最冷月均溫較低，苗齡低于3～4 a的紅花玉蘭無法自然越冬，因此，在引種此類苗齡紅花玉蘭時冬季需要采取適當的防寒措施，保護紅花玉蘭安全越冬。目前，園林綠化方面應用比較成熟的防寒措施有搭建防風障和使用保溫棉進行樹木包裹并覆土。夏季高日照會引起紅花玉蘭葉片日灼等傷害，需要搭建遮陽網，同時注意高溫干旱會引起紅花玉蘭苗木生長不適等狀況，需及時進行人工灌溉。③西北地區。從地域上看，該區包含地點分別屬于黃土高原中部，青藏高原東南部地區。該區以溫帶季風性氣候為主，地處半濕潤區，水、熱條件基本能滿足紅花玉蘭生長需求，但是最冷月均溫較低，3～4年生紅花玉蘭無法自然越冬，需要對紅花玉蘭進行越冬保護，如搭建防風障，使用保溫棉進行苗木包裹并覆土等，還要做好灌溉設施配套和管理，充分利用好水資源，常年對紅花玉蘭進行人工灌溉。同時在引進紅花玉蘭不同品種時，盡量挑選抗旱、抗寒性強的品種。

表9 紅花玉蘭引種次適宜區Table 9 Sub-suitable area for the introduction of Magnolia wufengensis

2.3.3 不適宜區 該區序號排在最后，且相對歐氏距離大于0.04，綜合氣候條件與紅花玉蘭原產地五峰具有較大差異（表10）。主要包括：①東北地區。該區從地域上看屬于東北平原，主要以溫帶季風性氣候為主，地跨濕潤區、半濕潤區。該區光照充足，但是最冷月均溫遠低于紅花玉蘭原產地五峰，冬春低溫干旱時間過長，不適宜紅花玉蘭生長。②西北地區。該區從地域上看屬于黃土高原西北部，青藏高原西北部，新疆，內蒙古等地，主要以溫帶大陸性氣候為主，地跨半干旱區和旱區。該區光照條件豐富，冬季寒冷干燥，夏季高溫，但是水、熱條件較差，最冷月均溫遠低于紅花玉蘭原產地，風、旱、雪、凍災害頻發，土地荒漠化嚴重，大部分是干旱地區，因此，在該區不適合引種紅花玉蘭。

表10 紅花玉蘭引種不適宜區Table 10 Un-suitable area for the introduction of Magnolia wufengensis

3 結論

紅花玉蘭在中國引種的適宜區包括黃淮平原、長江中下游平原、江南丘陵、浙閩丘陵、四川盆地東南部、云貴高原東北部和兩廣丘陵北部。該區受亞熱帶季風影響，夏季溫度偏高，雨熱同期，降水量多集中在夏季，容易形成圃地內澇和病蟲害，因此該區引種紅花玉蘭在夏季要及時做好圃地清理、病蟲害防治及排水工作。目前，生產上已經在上述引種適宜區中的安徽滁州，江蘇常州、南京，上海，河南南陽，湖南婁底，廣西柳州，貴州貴陽和云南昆明等地進行了小范圍引種并且生長表現良好。建議這些地區可以擴大引種。

紅花玉蘭在中國引種的次適宜區包括四川盆地西北部、云貴高原西南部、兩廣丘陵以南、雷州半島、海南島和臺灣島、晉冀魯山地丘陵地帶、華北平原、遼東半島、黃土高原中部及青藏高原東南部地區。該區南部臺風、洪澇等自然災害較多，引種紅花玉蘭需要注意地形地勢的選擇。北部地區持續低溫時間較長且極端氣溫較低，栽植苗齡為3～4年生及以下的紅花玉蘭易受凍害，應進行越冬防寒，西北部還需常年進行灌溉。目前，生產上已經在上述引種次適宜區的遼寧丹東、甘肅天水、河北石家莊、天津薊縣、北京海淀、山東臨沂等地也進行了小范圍引種紅花玉蘭。在夏季這些地區容易干旱導致紅花玉蘭生長不良，所以應增加灌溉，在冬季發現這些地區紅花玉蘭若不采取保護措施會發生生理干旱導致枯稍，因此，冬季應對紅花玉蘭進行防寒。

紅花玉蘭在中國引種的不適宜區包括東北平原、黃土高原西北部、青藏高原西北部、新疆和內蒙古等地。目前，陜西西安進行了紅花玉蘭少量引種，但是此地屬于紅花玉蘭引種不適宜區，若需種植必須采取保護地栽培。

4 討論

湖北、湖南、河南、浙江、福建、江西、安徽、云南、貴州、四川等省為紅花玉蘭的引種適宜區。這些地區屬于亞熱帶季風氣候，水熱條件較好，寒害凍害風險低，影響紅花玉蘭生長的極端氣候因子較少，非常適合紅花玉蘭的生長。該區天然分布木蘭屬樹種占到全國種的93.9%［30］，而且這些地區的木蘭屬樹種與變異最為豐富［31］，其中云南（10屬69種）、四川（5屬24種）、廣西（7屬36種）、廣東（5屬27種）等地是木蘭屬樹種的現代分布中心［32－33］，說明該區劃結果可靠，可為紅花玉蘭的大范圍引種提供科學依據。

紅花玉蘭的半致死溫度為－15℃，但是在適宜區中的江蘇沐陽、徐州、泗洪，安徽亳州，河南駐馬店等地極端氣溫低于－15℃，如果單從極端最低氣溫角度考慮，這些地區是不適宜引種紅花玉蘭的。因此像這樣的站點在引種紅花玉蘭時必須挑選適當的品種，充分考慮引種地小氣候條件，針對限制因子采取相對應的栽培措施。

在植物潛在引種氣候適生區區劃研究中，指標體系的選擇構建、指標權重確定對區劃結果的可靠性有直接影響。在構建區劃指標體系中，本研究綜合考慮了紅花玉蘭生長過程中所有可能限制紅花玉蘭分布的氣候因子，并在指標權重方面，考慮到各個因子對紅花玉蘭生長發育所起的作用（即權重）是不同的，因此，采用專家打分法和層次分析法相結合的方法確定權重，避免了以往單純依賴專家主觀打分帶來的局限性，更加注重內在邏輯性，使權重確定更加具有科學合理性。本研究僅分析了植物引種氣候相似性原理中的氣候因素，但是引種紅花玉蘭能否成功是其遺傳基因和環境共同作用的結果，此過程不僅受氣候因子的影響，當地的地形、海拔、土壤條件以及其他非生物因素等都可能影響其成活和生長［21］，例如胡建忠［24］在植物引種栽培試驗研究方法中所提到的風土馴化理論就考慮了不只氣候——“風”這一個因素，還有 “土”，即指土壤也影響著植物引種的成功。

5 參考文獻

［1］ 劉玉壺.木蘭科分類系統的初步研究［J］.植物分類學報，1984，22（2）：89－109.LAW Yuhwu.A preliminary study on the taxonomy of the family Magnoliaceae ［J］.Acta Phytotaxon Sin,1984,22（2）：89－109.

［2］ 張慶寶，申亞梅，范義榮.木蘭屬（Magnolia）觀賞植物育種現狀及育種策略［J］.江蘇林業科技，2008，35（6）：46－48,52.ZHANG Qingbao， SHEN Yamei， FAN Yirong.Current breeding situation and breeding strategy of ornamental plant breeding inMagnolia［J］.J Jiangsu For Sci Technol， 2008， 35（6）： 46 － 48， 52.

［3］ 王亞玲，李勇，張壽洲，等.用matK序列分析探討木蘭屬植物的系統發育關系 ［J］.植物分類學報，2006，44（2）： 135 － 147.WANG Yaling， LI Yong， ZHANG Shouzhou,et al.The utility ofmatKgene in the phylogenetic analysis of the genusMagnolia［J］.Acta Phytotaxon Sin， 2006， 44（2）： 135 － 147.

［4］ 林祁，段林東，袁瓊，等.單性木蘭屬（木蘭科）植物的分類學訂正［J］.武漢植物研究，2005，23（3）：236－238.LIN Qi， DUAN Lindong， YUAN Qiong，et al.Taxonomic notes on the genusKmeria（Pierre） Dandy （Magnoliaceae） ［J］.J Wuhan Bot Res， 2005， 23（3）： 236 － 238.

［5］ 馬履一，王羅榮，賀隨超，等.中國木蘭科木蘭屬一新種［J］.植物研究，2006，26（1）：4-7.MA Lüyi， WANG Luorong， HE Suichao，et al.A new species ofMagnolia（Magnoliaceae） from Hubei， China ［J］.Bull Bot Res， 2006， 26（1）： 4 － 7.

［6］ 馬履一，王羅榮，賀隨超，等.中國木蘭科木蘭屬一新變種［J］.植物研究，2006，26（5）：516-519.MA Lüyi， WANG Luorong， HE Suichao，et al.A new variety ofMagnolia（Magnoliaceae） from Hubei， China ［J］.Bull Bot Res， 2006， 26（5）： 516 － 519.

［7］ 梁大偉.紅花玉蘭優樹選擇與類型劃分［D］.北京：北京林業大學，2010.LINAG Dawei.Select on Superior Tree Magnolia wufengensis and Cultivar Classification［D］.Beijing：Beijing Forestry University,2010.

［8］ 桑子陽，馬履一，陳發菊，等.五峰紅花玉蘭種質資源保護現狀與開發利用對策［J］.湖北農業科學，2011，50（8）： 1564 － 1567.SANG Ziyang， MA Lüyi， CHEN Faju，et al.Protection status and utilization countermeasure of germplasm resources of theMagnolia wufengensisin Wufeng County ［J］.Hubei Agric Sci， 2011， 50（8）： 1564 － 1567.

［9］ 賀仲雄.模糊數學及其應用［M］.天津：天津科學技術出版社，1983.

［10］ ZHANG Huimin.Group decision making based on incomplete multiplicative and fuzzy preference relations ［J］.Appl Soft Comput， 2016， 48： 735 － 744.

［11］ 李整，文雨松，羅藝，等.基于典型斷面雨量-水位函數的中小橋水害預測［J］.中南大學學報（自然科學版），2010， 41（4）： 1597 － 1603.LI Zheng， WEN Yusong， LUO Yi，et al.Flood forecast method for middle ＆ small bridges based on rainfall-water level functions of representative riverbed sections ［J］.J Cent South Univ Sci Technol， 2010， 41（4）： 1597 － 1603.

［12］ 李麗純，陳家金，林晶，等.基于模糊相似優先比的閩臺兩地果樹氣候相似性分析［J］.中國生態農業學報，2013， 21（9）： 1149 － 1156.LI Lichun， CHEN Jiajin， LIN Jing，et al.Using fuzzy similarity priority ratio to analyse climate similarity between Fujian and Taiwan for fruit tree cultivation ［J］.Chin J Eco-Agric,2013， 21（9）： 1149 － 1156.

［13］ 韓蓓蓓.陜西省5種主要經濟林樹種氣候區劃研究［D］.楊凌：西北農林科技大學，2010.HAN Beibei.Study on Suitability Division of Five Main Economic Forest in Shaanxi Province［D］.Yangling：Northwest A＆F University,2010.

［14］ 張京偉，張德順，有祥亮，等.基于模糊相似優先比法劃分與上海氣候相似的美國區域［J］.江西農業大學學報， 2009， 31（6）： 1178 － 1182.ZHANG Jingwei， ZHANG Deshun， YOU Xiangliang，et al.Zoning of USA for Shanghai introduction with the method of fuzzy analogy peferred ratio ［J］.Acta Agric Univ Jiangxi， 2009， 31（6）： 1178 － 1182.

［15］ RADIONOVS A， UZHGA-REBROV O.Application of fuzzy logic for risk assessment ［J］.Inf Technol Manage Sci，2015， 17（1）： 50 － 54.

［16］ 袁嘉祖，馮晉層.模糊數學及其在林業中的應用［M］.北京：中國林業出版社，1988.

［17］ 袁嘉祖，鄧朝經.模糊相似優先比法的探討［J］.四川林業科技，1987,8（2）：93－96.YUAN Jiazu， DENG Chaojing.Discussion on method of fuzzy similarity priority ratio ［J］.J Sichuan For Sci Technol,1987,8（2）： 93 － 96.

［18］ 陸鼎煌，袁嘉祖.應用模糊相似優選比區劃中國油橄欖引種分布［J］.北京林學院學報，1982，4（4）：1－13.LU Dinghuang，YUAN Jiazu.Division of introduction regions of olive tree in China by means of indistinct simila-r optimum seeking ratio ［J］.J Beijing For Coll， 1982， 4（4）： 1 － 13.

［19］ 袁嘉祖.2個生態環境相似性計算方法的探討［J］.生態學雜志，1987，6（6）：55－56,45.YUAN Jiazu.Discussion on calculation method of similarity of two ecological environment ［J］.J Ecol,1987,6（6）：55－56,45.

［20］ 有祥亮，余新曉，陳吉虎，等.我國引種栓皮槭的氣候生態適生區劃分［J］.中國水土保持科學，2006，4（5）： 83 － 91.YOU Xiangliang， YU Xinxiao， CHEN Jihu，et al.Division on suitable climatic and ecological region for introducing ofAcer campestreL.to China ［J］.Sci Soil Water Conserv， 2006， 4（5）： 83 － 91.

［21］ 陳紅星，王勁風，龔榜初，等.我國引種日本甜柿氣候生態適生區初步分析［J］.林業科學研究，2000，13（3）： 323 － 327.CHEN Hongxing， WANG Jinfeng， GONG Bangchu，et al.A preliminary analysis on the suitable climatic and ecological region for introducing Japanese sweet persimmon ［J］.For Res， 2000， 13（3）： 323 － 327.

［22］ 潘靜，唐德瑞，車少輝，等.我國引種太平洋黃松的氣候適生區區劃［J］.西北林學院學報，2008，23（3）：80－84.PAN Jing， TANG Derui， CHE Shaohui，et al.Division on suitable climatic and ecological region for introducing ofPinus ponderosavar.ponderosato China ［J］.J Northwest For Univ， 2008， 23（3）： 80 － 84.

［23］ 王德英，丁國棟，趙媛媛，等.中國珍稀荒漠植物梭梭在半干旱沙地引種適應性［J］.北京林業大學學報，2015， 37（4）： 74 － 81.WANG Deying， DING Guodong， ZHAO Yuanyuan，et al.On suitability for introduction of Chinese rare desert plantHaloxylon ammodendronin semi-arid sandy regions of northwestern China ［J］.J Beijing For Univ， 2015， 37（4）： 74 － 81.

［24］ 胡建忠.植物引種栽培試驗研究方法［M］.鄭州：黃河水利出版社，2002.

［25］ 田士林，李莉.香樟在我國中部引種適應性研究［J］.安徽農業科學，2006，34（11）：2403－2445.TIAN Shilin， LI Li.Research on the adaptability of fragrant camphor tree in the central section in China ［J］.J Anhui Agric Sci， 2006， 34（11）： 2403 － 2445.

［26］ 李耀明，王玉杰，儲小院，等.降雨因子對縉云山地區典型森林植被類型地表徑流的影響［J］.水土保持研究， 2009， 16（4）： 244 － 249.LI Yaoming， WANG Yujie， CHU Xiaoyuan，et al.Effects of the rainfall factors on surface runoff of typical forest vegetations in three gorges region ［J］.Res Soil Water Conserv， 2009， 16（4）： 244 － 249.

［27］ 普宗朝，張山清，賓建華,等.氣候變化對烏魯木齊市冬小麥產量的影響［J］.西北農林科技大學學報（自然科學版）， 2013， 41（3）： 115 － 123.PU Zongchao， ZHANG Shanqing， BIN Jianhua，et al.Impact of climate change on winter wheat yield in Urumqi，Xinjiang ［J］.J Northwest A&F Univ Nat Sci Ed， 2013， 41（3）： 115 － 123.

［28］ 桑子陽，馬履一，陳發菊.干旱脅迫對紅花玉蘭幼苗生長和生理特性的影響［J］.西北植物學報，2011，31（1）： 109 － 115.SANG Ziyang， MA Lüyi， CHEN Faju.Growth and physiological characteristics ofMagnolia wufengensisseedlings under drought stress ［J］.Acta Bot Boreal-Occident Sin， 2011， 31（1）： 109 － 115.

［29］ 劉孟霞，張延龍，牛立新，等.運用層次-關聯分析法綜合評價加拿大引種草本花卉［J］.西北農業學報，2009， 18（4）： 261 － 266.LIU Mengxia， ZHANG Yanlong， NIU Lixin，et al.The comprehensive evaluation of the ornamental flower plants introduced from Canada by the application of hierarchy-relation analysis ［J］.Acta Agric Boreal-Occident Sin，2009， 18（4）： 261 － 266.

［30］ 閆雙喜，李永華，位鳳宇.中國木蘭科植物的地理分布［J］.武漢植物學研究，2008，26（4）：379－384.YAN Shuangxi， LI Yonghua， Wei Fengyu.Geographical distribution of Magnoliaceae plants in China ［J］.J Wuhan Bot Res， 2008， 26（4）： 379 － 384.

［31］ 何彥峰.我國木蘭屬植物研究進展［J］.北方園藝， 2010（3）： 186－190.HE Yanfeng.The research progress ofMagnoliaplant in China ［J］.Northern Hortic， 2010（3）： 186 － 190.

［32］ 劉玉壺，夏念和，楊惠秋.木蘭科（Magnoliaceae）的起源、進化和地理分布［J］.熱帶亞熱帶植物學報，1995,3（4）： 1 － 12.LIU Yuhu， XIA Nianhe， YANG Huiqiu.The origin， evolution and phytogeography of Magnoliaceae ［J］.J Trop Subtrop Bot， 1995， 3（4）： 1 － 12.

［33］ 吳征鎰.論中國植物區系的分區問題［J］.云南植物研究，1979，1（1）：1－20.WU Zhengyi.On the zoning of the flora of China ［J］.Acta Bot Yunan,1979,1（1） ：1 － 20.