黃山松研究綜述(1960－2014)及其知識圖譜分析

蘇松錦，劉金福，蘭思仁，洪 偉，李文周

(1.福建農(nóng)林大學，福建 福州350002;2.三明林業(yè)學校，福建 三明365001;3.福建農(nóng)林大學海峽自然保護區(qū)研究中心，福建福州350002;4.戴云山國家級自然保護區(qū)管理局福建泉州362503)

黃山松(Pinus taiwanensis Hsia)為中國特有樹種，是我國東部亞熱帶高海拔山地綠化、造林和生態(tài)恢復(fù)的重要樹種.在中山氣候惡劣環(huán)境下，黃山松具有馬尾松無法企及的生態(tài)演替功能;在高山植被恢復(fù)、固碳、水源涵養(yǎng)、針葉藥效開發(fā)以及應(yīng)對全球氣候變化等方面顯示出巨大的生態(tài)、社會與經(jīng)濟價值.戴云山擁有6400 hm2原生性黃山松群落，是世界最大的黃山松種質(zhì)基因基地，也是亞熱帶中山地區(qū)植被演替研究的理想地.為讓國內(nèi)外學者對黃山松最新研究動態(tài)有所了解，筆者對其分類與資源分布、遺傳學、生態(tài)學、生理生化、森林培育和其它方面進行梳理.另一方面，在全球氣候變暖、生物多樣性監(jiān)測以及森林生態(tài)系統(tǒng)定位監(jiān)測研究日益深入背景下，對黃山松進行知識圖譜分析可快速挖掘或預(yù)測其研究前沿，并追蹤其發(fā)展動向，可為后續(xù)研究提供合作契機和參考依據(jù).

1 黃山松研究進展

1.1 分類與資源分布

黃山松首個模式標本出自黃山西海門[1].不同學者對黃山松與臺灣松分類爭議較大，鄭萬鈞等認為二者過渡性變異尚未達到質(zhì)的飛躍，提議將其合并為黃山松，拉丁學名為Pinus taiwanensis Hsia，得到眾多學者的支持[2];但從二者的外部形態(tài)、木材解剖超微結(jié)構(gòu)、樹脂道結(jié)構(gòu)、細胞和遺傳學特征等看，不少學者認為黃山松是獨立種系[3，4].在馬尾松和黃山松基因漸滲區(qū)存在杉松、武陵松等形態(tài)特征介于二者之間的變種.不同學者對黃山松與其它松樹的分類關(guān)系可用圖1表示.

圖1 黃山松與其它松樹分類關(guān)系Fig.1 Relationship between Pinus-taiwanensis and other pine tree species

黃山松在豫、皖、浙、閩、臺、贛、湘、鄂、晉、桂、黔、滇等省均有分布[2].從水平分布看，可將其分為臺灣、南部、西南部、東南沿海、中西部、東北和北部大別山7個氣候生態(tài)區(qū)[5]和閩南區(qū)、閩贛區(qū)、湘西南區(qū)、皖浙區(qū)、鄂湘贛和豫皖大別山區(qū)6個地理類型[6];從垂直分布看，黃山松在長江流域以南和以北海拔分布的上下限分別約為 800－1800 m 和 500－1700 m[5].

1.2 遺傳學

中國大陸天然黃山松不同居群間遺傳一致度變化范圍為0.5－0.8，遺傳距離與地理距離呈顯著“S”型關(guān)系[6].從核型看，黃山松有4－5對具二次縊痕染色體，而臺灣松僅第2對染色體具二次縊痕[7];從遺傳多樣性看，黃山松核苷酸多樣性與種內(nèi)遺傳分化顯著高于馬尾松;二者在C3H座位上存在共享單倍型，種間差異GI大于C3H;GI基因座進化速度相對較快[8].在確定黃山松等松屬植物DNA多態(tài)性方面，采用整合微流控芯片優(yōu)于凝膠電泳方法[9].從譜系生物地理學看，黃山松、琉球松和臺灣松種間、種內(nèi)均存在大量花粉流;mtDNA主要遺傳變異發(fā)生在3種松樹居群間，而cpDNA遺傳變異則主要集中在群體內(nèi)部;黃山松群體歷史擴張事件可能發(fā)生時間早于最近一次第四紀盛冰川期;3種松樹在第四紀冰期可能存在多個避難所[4].黃山松與馬尾松在海拔垂直分布上存在自然雜交和基因漸滲區(qū)，在基因漸滲過程中，馬尾松基因向雜交群體漸滲頻率高于黃山松[10].為揭示黃山松與馬尾松的生殖隔離遺傳機制，李淑嫻[11]系統(tǒng)地開展了AFLP實驗程序優(yōu)化及引物組合篩選;遺傳圖譜構(gòu)建及標記偏分離和偏分離遺傳位點定位及相關(guān)基因組區(qū)域確定研究.

1.3 生態(tài)學

1.3.1 空間格局 植物空間格局研究有助于解決營造林中植株配置和采伐問題，是物種共存、森林更新與演替等研究的基礎(chǔ)[12].植物空間格局與尺度、計量參數(shù)和研究方法等有關(guān).同一物種在不同地區(qū)空間分布差異較大，如武夷山黃山松種群在高海拔呈聚集分布，而在低海拔呈隨機分布[13];福建壽寧縣黃山松種群呈聚集分布，聚集強度大樹＞幼樹＞中徑樹，并隨取樣尺度增加而減少[14];黃山地區(qū)黃山松種群幼年期呈聚集分布，隨著年齡增長逐漸向隨機分布過渡[15];而福建屏南縣黃山松種群則呈隨機分布[14].對于空間格局潛在生態(tài)學過程的解釋，多數(shù)學者認為黃山松總體上呈聚集分布，并隨年齡或海拔增加趨于隨機分布，與其自身特性、競爭及所處的土壤、地形和氣候等環(huán)境條件有關(guān)[15].負二項參數(shù)、擴散型指標Iδ、分形關(guān)聯(lián)維數(shù)和計盒維數(shù)等傳統(tǒng)計量參數(shù)與研究方法未考慮樣點位置信息，且只能分析單一尺度空間格局;而O-ring統(tǒng)計彌補了傳統(tǒng)指標的不足，且克服了Ripley's K函數(shù)隨尺度變化出現(xiàn)累積效應(yīng)的問題[16]，是當前植被空間格局研究的主流方法.空間格局不僅因物種而異，而且同一物種在不同發(fā)育階段、不同生境、不同林層也有明顯區(qū)別[12].隨著生物多樣性監(jiān)測網(wǎng)絡(luò)不斷拓展與研究的日益深入，采用O-ring統(tǒng)計在武夷山、戴云山、黃山和天柱山等地開展1－50 hm2大樣地尺度黃山松種群空間格局是未來值得關(guān)注的方向之一.

1.3.2 生態(tài)位 生態(tài)位是生態(tài)學研究核心內(nèi)容之一.戴云山黃山松8個群叢中91.78%種對生態(tài)位重疊值為0－0.5;與黃山松更新生態(tài)位重疊值較大的樹種有窄基紅褐柃(Eurya rubiginosa)和短尾越桔(Vaccinium carlesii)[17].武夷山黃山松林多數(shù)種群在坡度、坡向和海拔上生態(tài)適應(yīng)能力相似性，生態(tài)位普遍重疊，資源共享趨勢明顯[18].古田山黃山松、木荷(Schima superba)、野漆樹(Rhus sylvestris)和甜櫧(Castanopsis eyrei)等生態(tài)位寬度和重疊度均較大[19].種群生態(tài)位寬的物種，其生態(tài)位重疊不一定高，如閩西北黃山松生態(tài)位寬度較大(0.5－ 0.6)，但其與杉木(Cunninghamia lanceolata)生態(tài)位重疊值僅 0.3424[20].

1.3.3 種群結(jié)構(gòu)與動態(tài) 不同林型和同一林型不同密度下斷面積最大增長速度徑級范圍差異顯著[21].邏輯斯諦或廣義邏輯斯諦模型在擬合黃山松種群數(shù)量動態(tài)規(guī)律時，效果良好[22].從結(jié)構(gòu)與數(shù)量特征看，不同地區(qū)黃山松種群差異明顯:黃山黃山松種群死亡高峰期為幼苗1級向2級的過渡期，其幼苗在海拔較高的山頂、山脊和陡坡等生境下自然更新較好[15].天柱山黃山松種群隨海拔升高，林分密度增大，徑級分布范圍減小;種群存活曲線由Deevey-Ⅱ型向Deevey-Ⅰ型過渡，群落由穩(wěn)定型向增長型過渡，與物種競爭引起的自疏和它疏作用有關(guān)[23].而福建龍棲山和屏南縣的黃山松種群存活曲線整體上在DeeveyⅡ－Ⅲ型之間[22].

1.3.4 生態(tài)系統(tǒng)定位監(jiān)測 森林生態(tài)系統(tǒng)定位監(jiān)測主要包括水文、土壤、氣象、生物等內(nèi)容.林盛松[24]觀測了閩東北鷲峰山天然黃山松林的降雨分配格局;董金相[25]估測了戴云山黃山松林植被層、0－60 cm土壤層和枯枝落葉層的碳儲量，并分析其地形分異規(guī)律.韓勇等[26]研究了武夷山黃山松凋落葉在不同植物群落中的分解動態(tài);吳澤民等[27]分析了黃山景區(qū)黃山松年輪S及Pb、Mn、Zn和Cu等元素的動態(tài)特征.林益明和林鵬等[28，29]則系統(tǒng)地報道了武夷山黃山松各器官灰分含量、干重熱值、去灰分熱值以及N、P、K、Na累積與循環(huán)規(guī)律.朱德煌[17]采用除趨勢典范對應(yīng)分析法分析了環(huán)境因子對黃山松群落格局的影響.鄭仕平[30]發(fā)現(xiàn)黃山松林密度為3200株·hm－2的生產(chǎn)力和養(yǎng)分利用率最大.黃山松作為中國東部亞熱帶中山地區(qū)典型代表性群系，在全球變暖背景下，開展黃山松集水區(qū)與嵌套流域、蒸散量、根際微生態(tài)、小氣候、碳通量、痕量氣體與氣溶膠、粗木質(zhì)殘體、穩(wěn)定同位素和氮沉降等方面定位監(jiān)測研究，有助于揭示其生態(tài)系統(tǒng)結(jié)構(gòu)和功能的時空變化規(guī)律和對環(huán)境的響應(yīng)與適應(yīng)機制[31].

1.3.5 群落特征與演替規(guī)律 不同地區(qū)黃山松群落特征差異大.天柱山黃山松群落物種組成共有52科105屬，具有北亞熱帶常綠闊葉林和暖溫帶落葉闊葉林特征[32];古田山黃山松群落生活型以高位芽植物為主，群落葉以小型葉、單葉、草質(zhì)和非全緣為主[33];湖南平江幕阜山黃山松群落種子植物約40科66屬70種[34].從物種多樣性看，天柱山黃山松群落物種Shannon-Wiener指數(shù)按灌木層、草本層和喬木層降序排列[32];古田山黃山松群落物種多樣性灌木層＞喬木層＞草本層[33].黃山松群落多樣性與海拔、中度旅游干擾等有關(guān)[36，37].黃山松耐逆向人為干擾能力弱，且正向人為干擾下恢復(fù)進程緩慢[37].從演替規(guī)律看，黃山生境優(yōu)越地區(qū)，黃山松純林朝針、闊混交林演變，并最終被落葉闊葉林群落取代;而在裸巖、陡坡、峰嶺或梁崗等生境中，黃山松可形成穩(wěn)定地形頂級群落[15，38].在平江幕阜山，黃 山松群落演替規(guī)律為裸地→黃山松林→黃山松+ 落 葉闊葉樹混交林→落 葉闊葉林→常綠闊葉林[34].從干擾角度看，按5.3%/10 a林冠干擾強度算，西天目山黃山松闊葉林在2051－2071年可能完全被闊葉純林取代[39].從更新角度看，吳澤民等[40]發(fā)現(xiàn)黃山松林隙形成與更新發(fā)生時間相近，40 m2林隙或110 m2擴展林隙的光照能保證黃山松幼苗的正常生長.

1.3.6 競爭 安徽天柱山黃山松種內(nèi)競爭強度隨對象木胸徑增大而降低，采用逐步擴大范圍法可有效地確定其鄰體競爭范圍[41].為精確預(yù)測林分生長，在構(gòu)建林分競爭生長模型時需綜合考慮林木胸徑、冠幅和林齡等.黃山松種群自疏效應(yīng)明顯，利用枝條著枝深度和林木變量建立的樹冠形狀模型及密度效應(yīng)模型，可測算其生產(chǎn)彈性值和邊際產(chǎn)量[13].競爭格局與種群自身、環(huán)境等有關(guān)，對不同階段黃山松幼樹進行適當透光伐，可提高其與闊葉樹的競爭能力[42].黃山松種內(nèi)競爭在壽寧縣小于屏南縣;在閩西北，天然黃山松個體空間關(guān)聯(lián)較強、競爭激烈和空間占據(jù)程度高的海拔范圍為1300－1310 m和1400 m[43].林分密度控制對林分經(jīng)營具有科學指導(dǎo)意義，考慮競爭因素，福建壽寧縣黃山松合理栽植密度應(yīng)控制在4500株·hm－2[44].

1.4 生理生化

黃山松適應(yīng)高山嚴酷環(huán)境能力強，其生理生態(tài)穩(wěn)定性主要通過改變表型來實現(xiàn).不同海拔黃山松針葉長、單根葉重和比葉面積差異顯著[45].黃山松中生型、邊生型、內(nèi)生型和鞘內(nèi)型針葉樹脂道與海拔關(guān)系密切，針葉樹脂道數(shù)目主要受≥10℃積溫、日照時數(shù)和降水量制約，年輪結(jié)構(gòu)和軸向樹脂道緯向變異明顯[6].從光合特性看，黃山松不同年齡段1年生針葉時蒸騰強度、水分利用率日變化均呈雙峰分布;黃山松幼齡受氣孔限制明顯，中幼齡抗旱能力最強[46].黃山松光合速率夏季＞秋季＞春季＞冬季，春、夏、秋3季有光合午休現(xiàn)象;其光合速率主要受細胞間隙CO2濃度和氣孔導(dǎo)度等影響[47].由于光補償點和光飽和點均較高，黃山松幼苗在林下競爭中處于劣勢，更新較困難.從生化角度看，低溫處理下黃山松幼苗浸出液相對電導(dǎo)率、超氧化物歧化酶活性及丙二醛質(zhì)量摩爾濃度均隨海拔升高而降低[48].針葉葉綠素和蛋白質(zhì)含量可反映黃山松的衰老狀況[49].

1.5 森林培育

1.5.1 生長過程分析 生長過程是植物定向培育研究的理論基礎(chǔ).Logistic模型可較好地模擬黃山松生長過程.黃山松天然林胸徑生長量比人工林低0.47 cm·a－1;樹高連年生長量最高峰(第20年)和材積快速生長點(第40年)比人工林分別遲5年和20年[50].黃山松林分形高隨林分樹高增大而增大，而隨胸高斷面積或疏密度增大而減少[51].從生物量看，密度為1530株·hm－2的10年生黃山松人工林生產(chǎn)力約1.675 t/(hm2·a)[52].疏伐可顯著提高球果生物量，當黃山松種子園郁閉度為0.5時，球果產(chǎn)量最高[53].采用經(jīng)驗公式W=a+bD+cD2預(yù)測黃山松生物量可達到較高精度，但速生階段宜用W=a+bD[54].以立地質(zhì)量為啞變量、間隔期為自變量的黃山松林分收獲模型，可有效解決天然異齡林生長動態(tài)預(yù)測中林齡不易確定問題，并為擇伐收益法在資產(chǎn)評估中的推廣應(yīng)用提供技術(shù)支撐[55].

1.5.2 種子園建設(shè) 黃山松優(yōu)樹選擇可采用數(shù)量與形質(zhì)指標，個體樹與對比樹相結(jié)合方法[56];高海拔黃山松采種最佳時間段為霜降前后，采種時間越早，種子品質(zhì)越差[57].黃山松較低的種子發(fā)芽率(21.6%)與綜合逆境活力(15.5%)一定程度上影響了其繁殖能力[15];在陰天對高約12 cm黃山松苗進行切根，可促進其徑粗生長與根系發(fā)達，并顯著提高其菌根感染率與造林成活率[58].黃山松無性系分株主要依靠嫁接來完成，嫁接時應(yīng)優(yōu)先考慮嫁接部位和接穗粗度.浙江天臺縣華頂林場黃山松種子園無性系物侯期和生長節(jié)律、土壤改良、樹體管理、定砧與嫁接、花粉管理及病蟲害防治等方面營建技術(shù)先進，效果良好，值得示范推廣[59].

1.5.3 撫育與造林成效 撫育與造林成效分析是黃山松林集約化經(jīng)營的重要環(huán)節(jié).黃山松幼林生長主要受坡位、種源和造林密度影響，而整地方式、施肥和幼林撫育對其影響較小[60].黃山松與光皮樺2:1混交，林分生長最好，并可有效解決杉木采伐跡地更新困難問題[61].編制黃山松數(shù)量化地位指數(shù)表可擴展地位質(zhì)量評定適用范圍[62].高海拔山地杉木、黃山松人工混交林對幼齡杉木個體樹高和中齡杉木個體胸徑有促進作用，其立木蓄積、胸高斷面積和生物量明顯高于杉木純林[63].

1.5.4 病蟲害防治 對黃山松有威脅的病蟲害主要有:松材線蟲(Bursaphelenchus xylophilus)、松毛蟲(Dendrulimus marmortus)、小卷蛾(Petrova cristata)、松果梢斑螟(Dioryctria pryeri)、松針褐斑病(Lecanosticta acicola)、細紋新須螨(Cenopalpus lineola)、小環(huán)綿盤菌(Cyclaneusmam inus)、黃山松黑松葉蜂(Nesodiprion.SP)、松斑螟(Salebria laruata)、松櫟銹病、赤枯病、紫色病(缺鎂)等[38，64].其中學者最為關(guān)注:黃山松是否為松材線蟲感病樹種[65－67]?黃山松感染松材線蟲可能與供試苗、接種源以及接種氣候條件等有關(guān)，自然狀態(tài)下受松材線蟲侵染的黃山松罹病木呈現(xiàn)出針葉部分變紅、木材有松香味、有松褐天牛危害痕跡等特點;黃山松感病木僅分布在其海拔下限區(qū)和人工栽植的林分中，而較高海拔天然黃山松林中未見有感染松材線蟲報道[65].篩選病蟲害抗性候補木是種子園建設(shè)的重要工作.陳鳳毛等[66]發(fā)現(xiàn)擬松材線蟲CFS4，SLD3，ZJ1和KOR1四個株系對黃山松苗有致病作用，其中接種ZJ1或SLD3，黃山松苗死亡率均為16.7%.徐六一等[67]于2005－2007年間篩選了57個家系、216株黃山松抗性候補木，隨后對篩選出的抗性候選個體開展MuPS標記識別分析.

1.6 其它

樹輪氣候?qū)W:樹木年輪分析對了解影響植物生長關(guān)鍵因子和重建歷史氣候資料意義重大.不同地區(qū)黃山松直徑生長及其影響因子差異較大，黃山黃山松徑向生長主要受當年＞10℃積溫的太陽輻射和積溫影響[68];大別山黃山松徑向生長主要受當年2－7月平均氣溫限制[69];九龍山自然保護區(qū)和安徽省南部牯牛降的高海拔黃山松樹輪徑向生長主要受夏季(上年6－7月和當年6月)水熱條件影響[70];而西大別山黃山松樹輪寬度生長主要受9、10月的水熱組合影響[71].

針葉的開發(fā)與利用:黃山松松針功能性成分研究目前僅見于馬文秀[72]鑒定出8，15-異海松二烯酸、8，15-海松二烯酸和濕地松酸;劉曉庚等[73]報道了黃山松的粗蛋白質(zhì)、多聚戊糖、纖維素、灰分、葉綠素和胡蘿卜素等理化成分含量.黃山松松針在保健食品和藥業(yè)生產(chǎn)上的開發(fā)與應(yīng)用有待深入.

災(zāi)害分析:黃山松受冰雪災(zāi)害，易發(fā)生樹干彎曲、斷干、掘根、斷冠和凍死等現(xiàn)象，樹高徑粗易斷大枝和斷冠;樹矮徑細、冠幅偏大易斷干或被掘根;幼樹易被凍死;尖削度小、樹冠較高且不均勻的黃山松抵御冰雪災(zāi)害能力較弱[74].

2 黃山松文獻知識圖譜分析

以“黃山松”為主題詞，Web of Science和中國知網(wǎng)CNKI為數(shù)據(jù)庫，對1960－2014年與林學相關(guān)學科領(lǐng)域文獻進行檢索，共得到461條記錄;去重、辨別(刪除書評及會議等)和整理后，獲得文獻220篇.采用CitespaceⅢ軟件[75]對文獻進行格式轉(zhuǎn)換，并對其關(guān)鍵詞、作者、機構(gòu)和期刊進行分析.

2.1 時間分布

黃山松受到關(guān)注和重視始于20世紀90年代，1995年以后，其發(fā)文數(shù)量呈波動上升.圖2中發(fā)文數(shù)量較多的年份為2000年(14篇)、2006年(12篇)和2012年(13篇).220篇文獻中僅有6篇被SCI收錄，表明黃山松研究國際社會影響力較弱.

圖2 1960－2014年黃山松林學相關(guān)文獻年篇數(shù)Fig.2 Papers of forestry literature on Pinus taiwanensis per year from 1960 to 2014

圖3 1960－2014年黃山松研究關(guān)鍵詞的知識圖譜Fig.3 Mapping knowledge domain of keywords in P.taiwanensis research from 1960 to 2014

2.2 關(guān)鍵詞分析

將CitespaceⅢ軟件中時區(qū)設(shè)為3年，以關(guān)鍵詞為網(wǎng)絡(luò)節(jié)點，基于各時區(qū)中高被引前10個關(guān)鍵詞，生成關(guān)鍵詞共被引聚類知識圖譜(圖3).關(guān)鍵詞知識圖譜網(wǎng)絡(luò)中共有118個節(jié)點，107個連接，網(wǎng)絡(luò)密度為0.0155.網(wǎng)絡(luò)中節(jié)點大小表示文獻數(shù)量，連線表示節(jié)點間合作關(guān)系，連線粗細表示合作的頻次，中心性表示節(jié)點在網(wǎng)絡(luò)中的重要性.從圖3可知，黃山松知識圖譜網(wǎng)絡(luò)中重要性大于0.08的關(guān)鍵詞為黃山松、馬尾松、樹脂道、天然更新和遺傳多樣性;頻次大于6的關(guān)鍵詞為黃山松、馬尾松、松屬、林分、生物量、樹脂道、天然更新和種群(表1).

表1 1960－2014年黃山松研究前24個高頻關(guān)鍵詞Table 1 Top 24 keywords of P.taiwanensis research appearing frequently from 1960 to 2014

2.3 發(fā)文作者分析

作者共被引聚類知識圖譜網(wǎng)絡(luò)中共有327個節(jié)點，251個連接，網(wǎng)絡(luò)密度為0.0047.節(jié)點越大，表明作者發(fā)文頻數(shù)越高.洪偉在圖譜中出現(xiàn)頻次最高，主要從事黃山松種群生態(tài)學研究;其次為范義榮、吳承禎、吳澤民和羅世家等(圖4)，分別從事林木遺傳與良種選育、種群生態(tài)學、森林生態(tài)學和植物學研究.發(fā)文較多作者之間合作關(guān)系密切，如吳承禎與洪偉，其它作者以獨自研究狀態(tài)為主.老師及其指導(dǎo)的學生是研究團體的核心部分.如福建農(nóng)林大學洪偉與吳承禎、封磊、鄭世群、劉金福和宋萍等;浙江農(nóng)林大學范義榮與其學生，安徽農(nóng)業(yè)大學吳澤民與其學生.發(fā)文數(shù)量在4篇以上作者共28人(表2).

圖4 1960－2014年黃山松研究作者的知識圖譜Fig.4 Mapping knowledge domain of authors in P.taiwanensis research(1960－2014)

2.4 發(fā)文期刊分析

參照《北京大學中文核心期刊目錄(2011年版)》，黃山松研究前19位期刊118篇文獻有78篇為農(nóng)業(yè)科學類核心期刊.220篇文獻中SCI期刊、核心期刊、博碩士論文分別有6、121和14篇，占總數(shù)的64.09%，表明發(fā)文質(zhì)量相對較高(表3).

表2 1960－2014年黃山松研究前28個作者Table 2 Top 28 authors in P.taiwanensis research from 1960 to 2014

黃山松發(fā)文機構(gòu)共被引聚類知識圖譜網(wǎng)絡(luò)中共有145個節(jié)點，48個連接，網(wǎng)絡(luò)密度為0.0046.發(fā)文頻次最高3家機構(gòu)分別為福建農(nóng)林大學、安徽農(nóng)業(yè)大學和浙江農(nóng)林大學(圖5).福建農(nóng)林大學與戴云山國家級自然保護區(qū);黃山學院與黃山風景區(qū)園林局;南京林業(yè)大學和安徽省林業(yè)勘察設(shè)計院、湖北民族學院等之間合作關(guān)系較緊密，其它機構(gòu)研究較獨立，如安徽農(nóng)業(yè)大學、浙江農(nóng)林大學、福建省壽寧縣林業(yè)局等.

表3 1960－2014年黃山松研究前19位期刊Table 3 Top 19 journals in Pinus taiwanensis research from 1960 to 2014

圖5 1960－2014年黃山松研究發(fā)文機構(gòu)的知識圖譜Fig.5 Mapping knowledge domain of institutes in P.taiwanensis research from 1960 to 2014

表4 1960－2014年黃山松研究前17個機構(gòu)Table 4 Top 17 institutes of P.taiwanensis research from 1960 to 2014

3 結(jié)論與展望

黃山松遺傳學和生理生化研究結(jié)果至今仍未能較好地解決黃山松分類問題，但多數(shù)學者傾向于將黃山松作為獨立種系，最常用學名為P.taiwanensis Hayata.目前生物多樣性、生態(tài)系統(tǒng)定位觀測、植被維持機制及其與環(huán)境因子耦合關(guān)系等研究，多基于1－50 hm2大樣地定位觀測數(shù)據(jù).戴云山國家級自然保護區(qū)6400 hm2黃山松群落是中國東南沿海針葉林生態(tài)系統(tǒng)的典型代表，其科研重點應(yīng)放在針葉林大樣地建設(shè)、空間格局和生態(tài)系統(tǒng)定位監(jiān)測上.黃山松人工林培育過程中應(yīng)重視良種選育以及病蟲害防治.黃山松徑向生長及其對氣候變化響應(yīng)結(jié)果區(qū)域差異大.因此，戴云山黃山松樹輪寬度氣候?qū)W分析是一個極有意義的研究點.黃山松在全國分布范圍較廣，但其景觀生態(tài)學研究幾乎空白，在景觀、區(qū)域或更大尺度上開展其碳儲量(生物量)、葉面積指數(shù)LAI、植被凈生產(chǎn)力NPP的反演研究，可豐富陸地生態(tài)系統(tǒng)和生物多樣性與全球變化的相互作用機制研究.從關(guān)鍵詞、作者和發(fā)文機構(gòu)知識圖譜看，黃山松林業(yè)領(lǐng)域研究以獨立研究狀態(tài)為主，合作研究較少，建議今后應(yīng)加強南北方黃山松林景觀或區(qū)域等大尺度合作研究.

[1]夏緯英.黃山松[J].中國植物學雜志，1936(1):17.

[2]中國科學院中國植物志編輯委員會.中國植物志(第七卷)[M].北京:科學出版社，1978.

[3]鄒惠渝，黎志，羅世家，等.黃山松的分類學研究[J].植物研究，2003，23(3):278－279.

[4]余雷.黃山松的譜系生物地理學研究[D].合肥:安徽農(nóng)業(yè)大學，2012.

[5]童再康，范義榮.黃山松分布區(qū)氣候生態(tài)區(qū)劃[J].中南林學院學報，1993，13(1):81－87.

[6]石雷.中國大陸地區(qū)黃山松地理變異的形態(tài)和RAPD分析[D].合肥:安徽農(nóng)業(yè)大學，2012.

[7]邢有華，方永鑫，歐陽家安.黃山松的核型研究[J].安徽農(nóng)學院學報，1985(2):72－76.

[8]張利銳，彭艷玲，任廣朋，等.馬尾松和黃山松兩個核基因位點的群體遺傳多樣性和種間分化[J].植物生態(tài)學報，2011，35(5):531－538

[9]唐娟娟，吳江林，劉金華，等.用整合微流控芯片系統(tǒng)分析黃山松RAPD片段多態(tài)性[J].遺傳，2003，25(4):419－424.

[10]羅世家，鄒惠渝，梁師文.黃山松與馬尾松基因漸滲的研究[J].林業(yè)科學，2001，37(6):118－122.

[11]李淑嫻.馬尾松和黃山松物種分化遺傳機制研究[D].南京:南京林業(yè)大學，2012.

[12]HOU JH，Mi XC，LIU CR，et al.Spatial patterns and associations in a Quercus-Betula forest in northern China[J].Journal of Vegetation Science，2004，15:407－414.

[13]蔡小英.武夷山黃山松種群結(jié)構(gòu)與動態(tài)研究[D].福州:福建農(nóng)林大學，2008.

[14]宋萍，洪偉，吳承禎，等.天然黃山松種群空間格局的分形特征——關(guān)聯(lián)維數(shù)[J].武漢植物學研究，2005，23(1):32－36.

[15]羅世家.黃山松的研究[D].南京:南京林業(yè)大學，1998.

[16]Wiegend T，Moloney K A.Rings，circles and null models for point pattern analysis in ecology[J].Oikos，2004，104:209－229.

[17]朱德煌.戴云山黃山松種群特征及其環(huán)境因子研究[D].福州:福建農(nóng)林大學，2013.

[18]封磊，洪偉，吳承禎，等.武夷山黃山松林主要種群生態(tài)位特征[J].山地學報，2008，26(3):308－316.

[19]胡正華，于明堅，彭傳正，等.古田山國家自然保護區(qū)黃山松林主要種群生態(tài)位研究[J].生態(tài)環(huán)境，2004，13(4):619－621，629.

[20]林思祖，黃世國，洪偉，等.杉闊混交林主要種群多維生態(tài)位特征[J].生態(tài)學報，2002，22(6):962－968.

[21]吳承楨，洪偉，林成來.黃山松種群數(shù)量動態(tài)研究[J].浙江林學院學報，1998，15(3):274－279.

[22]畢曉麗，洪偉，吳承禎.黃山松種群統(tǒng)計分析[J].林業(yè)科學，2002，38(1):61－67.

[23]段仁燕，黃敏毅，吳甘霖.黃山松種群結(jié)構(gòu)與動態(tài)研究[J].生物學雜志，2009，26(1):35－37，34.

[24]林盛松.閩東北鷲峰山天然黃山松林內(nèi)降雨分配的研究[J].福建林業(yè)勘察設(shè)計，2007(2):31－34.

[25]董金相.戴云山黃山松林碳儲量及其影響因子研究[D].福州:福建農(nóng)林大學，2012.

[26]韓勇，徐憲根，阮宏華，等.武夷山黃山松凋落葉在不同植物群落中的分解動態(tài)[J].南京林業(yè)大學學報:自然科學版，2010，34(3):141－145.

[27]吳澤民，高健，黃成林，等.黃山松年輪硫及重金屬元素含量動態(tài)特征[J].應(yīng)用生態(tài)學報，2005，16(5):820－824.

[28]林益明，林鵬.福建武夷山2個典型植物群落建群種的熱值研究[J].武夷科學，1999，15:118－123.

[29]林益明，李振基，楊志偉，等.福建武夷山黃山松群落的氮、磷累積和循環(huán)[J].熱帶亞熱帶植物學報，1997，5(2):26－32.

[30]鄭仕平.黃山松人工林養(yǎng)分循環(huán)與密度關(guān)系的研究[J].防護林科技，2014(8):13－16.

[31]丁訪軍.森林生態(tài)系統(tǒng)定位研究標準體系構(gòu)建[D].中國林業(yè)科學研究院，2011.

[32]周葆華，余世金.天柱山黃山松群落特征及其環(huán)境功能評價[J].地理研究，2008，27(2):257－265.

[33]胡正華，于明堅，丁炳揚，等.古田山國家級自然保護區(qū)黃山松群落特征及物種多樣性研究[J].生態(tài)環(huán)境，2003，12(4):436－439.

[34]李家湘，趙麗娟，黃展鵬.平江幕阜山黃山松群落特征及其演替規(guī)律的探討[J].湖南林業(yè)科技，2004，31(5):16－18.

[35]任國學，劉金福，徐道煒，等.戴云山國家級自然保護區(qū)黃山松群落類型與物種多樣性分析[J].植物資源與環(huán)境學報，2011，20(3):82－88.

[36]吳甘霖，黃敏毅，段仁燕，等.不同強度旅游干擾對黃山松群落物種多樣性的影響[J].生態(tài)學報，2006，26(12):3924－3930.

[37]姚劍飛，方樂金，吳俊，等.干擾對黃山高海拔木本植被生長的影響[J].東北林業(yè)大學學報，2012，40(3):28－31.

[38]季春峰，鄒惠渝，向其柏.黃山松研究進展(綜述)[J].安徽農(nóng)業(yè)大學學報，2004，31(1):111－114.

[39]夏冰，蘭濤，賀善安，等.西天目山黃山松闊葉林的冠層干擾與動態(tài)推測[J].植物資源與環(huán)境，1995，4(3):15－20.

[40]吳澤民，黃成林，韋朝領(lǐng).黃山松群落林隙光能效應(yīng)與黃山松的更新[J].應(yīng)用生態(tài)學報，2000，11(1):13－18.

[41]段仁燕，黃敏毅，吳甘霖，等.黃山松種群鄰體范圍與鄰體競爭強度的研究[J].廣西植物，2009，29(1):111－115.

[42]龔道懷.次生林內(nèi)黃山松幼樹的生長及其量化表達[J].林業(yè)勘察設(shè)計(福建)，2009(1):102－105.

[43]畢曉麗，洪偉，吳承禎，等.黃山松林不同樹種樹冠分形特征研究[J].福建林學院學報，2001，21(4):347－350.

[44]張進德.黃山松造林密度與生物量的關(guān)系[J].林業(yè)科技通訊，2000(1):7－10.

[45]馬元燦，王中生，余華，等.不同海拔梯度下黃山松與馬尾松針葉形態(tài)·光合生理特性的研究[J].安徽農(nóng)業(yè)科學，2012，40(29):14155－14158，14173.

[46]武金翠，吳澤民，洪淑媛，等.不同年齡黃山松蒸騰速率及水分利用效率的比較研究[J].安徽農(nóng)業(yè)大學學報，2007，34(2):239－243.

[47]洪淑媛，吳澤民，劉西軍.黃山松生理生態(tài)特性研究[J].安徽農(nóng)業(yè)大學學報，2006，33(4):462－467.

[48]李淑嫻，譚艷，陳穎，等.黃山松不同生理生化指標隨海拔高度變化趨勢[J].東北林業(yè)大學學報，2010，38(6):9－10，34.

[49]李東林，嚴景華，曹恒生，等.黃山松不同齡階針葉衰老指標的比較研究[J].林業(yè)科學研究，1998，11(2):218－221.

[50]盧琦，趙體順，羅天祥，等.黃山松天然林與人工林物種多樣性和林分生長規(guī)律的比較研究[J].林業(yè)科學研究，1996，9(3):273－277.

[51]劉欽.黃山松人工林形高與斷面積和疏密度的關(guān)系[J].林業(yè)科技，2003，17(增刊):19－20.

[52]胡道連，李志輝，謝旭東.黃山松人工林生物產(chǎn)量及生產(chǎn)力的研究[J].中南林學院學報，1998，18(1):60－64.

[53]陳見陽.黃山松種子園球果產(chǎn)量與郁閉度的研究[J].浙江林業(yè)科技，2014，34(1):58－60.

[54]鄭德華.黃山松人工林生物量測定方法的比較[J].福建林業(yè)科技，1999，26(增刊):36－39.

[55]華偉平，池上評，李曉偉，等.黃山松天然林收獲模型及在擇伐收益法中的應(yīng)用[J].福建林學院學報，2014，34(3):220－224.

[56]劉通高，張進德，范志平.黃山松天然林優(yōu)樹選擇標準和方法的研究[J].福建林業(yè)科技，1999，26(增刊):54－57.

[57]王國英，占曉云，范義榮.不同采收期對黃山松種子品質(zhì)的影響[J].浙江林業(yè)科技，1998，18(1):43－45.

[58]林盛松.黃山松大田苗切根培育技術(shù)研究[J].福建林業(yè)勘察設(shè)計，2006(2):122－125.

[59]范義榮，毛迎春，方陸明，等.黃山松種子園營建技術(shù)及初步效果[J].浙江林學院學報，1997，14(2):111－119.

[60]陳蘭仙.造林技術(shù)措施對黃山松幼林生長的效應(yīng)[J].福建林學院學報，2000，20(4):357－360.

[61]黃紹平，陶義俊，王興平，等.慶元林場黃山松與光皮樺人工林混交營造比例研究[J].現(xiàn)代農(nóng)業(yè)科技，2011(11):217，223.

[62]章禮拐.皖西黃山松人工林地位質(zhì)量數(shù)量化評定的研究[J].安徽林業(yè)科技，1995(2):20－22.

[63]徐小牛，李宏開.杉木、黃山松人工混交林的生長特性[J].林業(yè)科學研究，1996，9(3):278－283.

[64]方樂金，李豐伯，萬志兵，等.黃山松細紋新須螨發(fā)生與氣象因素關(guān)系的分析[J].西南林業(yè)大學學報，2014，34(2):104－107.

[65]蔣麗雅，江順利，汪振宇，等.黃山松自然狀態(tài)下感染松材線蟲病的初步調(diào)查[J].安徽農(nóng)業(yè)大學學報，2006，33(1):5－8.

[66]陳鳳毛，史延梅，王姝穎，等.不同株系擬松材線蟲對黃山松和黑松苗木的致病性[J].林業(yè)科學，2011，46(12):86－90.

[67]徐六一，高景斌，蔡衛(wèi)兵，等.黃山松松材線蟲病抗性候補木的選拔及評價研究[J].安徽農(nóng)業(yè)科學，2008，36(32):14093－14094，14115.

[68]吳澤民，黃成林，馬青山.黃山松年輪生長和氣候的關(guān)系[J].應(yīng)用生態(tài)學報，1999，10(2):147－150.

[69]鄭永宏，張永，邵雪梅，等.大別山地區(qū)黃山松和油松樹輪寬度的氣候意義[J].地理科學進展，2012，31(1):72－77.

[70]李玲玲，史江峰，侯鑫源，等.中國東南高海拔黃山松生長對氣候的響應(yīng)——以浙江省九龍山和安徽省枯牛降為例[J].應(yīng)用生態(tài)學報，2014，25(7):1849－1856.

[71]彭劍峰，李國棟，李玲玲.西大別山小林海黃山松樹輪寬度的氣候意義[J].應(yīng)用生態(tài)學報，2014，25(7):1857－1862.

[72]馬文秀.黃山松松脂樹脂酸少量成份研究[J].南京林業(yè)大學學報，1989，13(4):69－73.

[73]劉曉庚，陳梅梅.我國松針的開發(fā)利用研究及進展[J].林產(chǎn)化工通訊，2003，37(4):29－33.

[74]張志祥，劉鵬，邱志軍，等.浙江九龍山自然保護區(qū)黃山松種群冰雪災(zāi)害干擾及其受災(zāi)影響因子分析[J].植物生態(tài)學學報，2010，34(2):223－232.

[75]CHEN C.CiteSpace II:Detecting and visualizing emerging trends and transient patterns in scientific literature[J].Journal of the American Society for Information Science and Technology，2006，57(3):359－377.?