李友邁 趙光影

摘要 天然林保護(hù)工程的實(shí)施改變了森林景觀日益破碎化的格局，因此通過(guò)研究天保工程實(shí)施背景下森林景觀格局的變化，可以說(shuō)明天保工程背景下經(jīng)營(yíng)模式的轉(zhuǎn)變是否對(duì)森林景觀變化趨勢(shì)起到積極的作用；同時(shí)為今后天然林保護(hù)工作措施的改進(jìn)提供科學(xué)依據(jù)。通過(guò)研究利用景觀格局指數(shù)計(jì)算方法，對(duì)1987—2017年露水河林業(yè)局實(shí)行天然林保護(hù)工程前后其森林景觀格局的變化進(jìn)行分析，結(jié)果表明：研究區(qū)域內(nèi)林地為優(yōu)勢(shì)景觀；非林地面積占比較小且面積逐年遞減，林地內(nèi)部相互轉(zhuǎn)化，其中大約26.2%～38.6%的針闊葉林轉(zhuǎn)變?yōu)榛旖涣帧Ｑ芯繀^(qū)域內(nèi)NP指數(shù)下降到7 695， SHDI指數(shù)始終保存在1.1～1.4。研究區(qū)域內(nèi)商品林區(qū)NP指數(shù)由9 905下降到9 225，公益林區(qū)NP指數(shù)由從13 474下降到11 001。商品林區(qū)，公益林區(qū)SHDI指數(shù)在0.9～1.2浮動(dòng)，SHEI指數(shù)較高。景觀類型總體上分布均勻，多樣性水平較低。

關(guān)鍵詞 景觀格局空間變化；景觀破碎化；景觀多樣性

中圖分類號(hào) S718.5? 文獻(xiàn)標(biāo)識(shí)碼 A

文章編號(hào) 0517-6611（2023）14-0109-05

doi：10.3969/j.issn.0517-6611.2023.14.027

作者簡(jiǎn)介 李友邁（1998—），男，遼寧沈陽(yáng)人，碩士研究生，研究方向：生態(tài)學(xué)。*通信作者，教授，從事地表生態(tài)與生態(tài)調(diào)控研究。

景觀是各種具有特定結(jié)構(gòu)、功能和動(dòng)態(tài)特征的生態(tài)系統(tǒng)鑲嵌組成的異質(zhì)性區(qū)域［ 1］。景觀的異質(zhì)性導(dǎo)致各景觀要素的構(gòu)成具有多樣性的特征。景觀的多樣性揭示了景觀的復(fù)雜性，對(duì)景觀功能和過(guò)程有重要的影響［ 2-3］。大小形狀各異的景觀要素的空間分布為景觀格局，景觀格局隨時(shí)間的變化稱為景觀動(dòng)態(tài)變化［ 4］。景觀動(dòng)態(tài)變化是地理學(xué)和景觀生態(tài)學(xué)領(lǐng)域中長(zhǎng)期關(guān)注的熱點(diǎn)問(wèn)題［ 5］，其中森林景觀格局變化及其驅(qū)動(dòng)力研究歷來(lái)備受關(guān)注［ 6］。森林景觀格局的變化會(huì)對(duì)森林生態(tài)系統(tǒng)內(nèi)的物質(zhì)循環(huán)和能量流動(dòng)產(chǎn)生影響，制約多種生態(tài)過(guò)程［ 7］。森林景觀格局變化的原因在于外界的干擾作用［ 8］，外界的干擾作用總體包括自然因素和人為因素；其中人為因素是導(dǎo)致景觀出現(xiàn)變化的主要原因［ 9］。景觀格局變化的主要表現(xiàn)為破碎化現(xiàn)象，景觀破碎化會(huì)引起生態(tài)系統(tǒng)的一系列變化，進(jìn)而影響生物多樣性和生態(tài)功能的發(fā)揮，因此景觀破碎化是研究景觀格局變化的重要內(nèi)容之一［ 10-11］。

長(zhǎng)白山地區(qū)是我國(guó)森林的重要區(qū)域，由于受到多年的人為開發(fā)和破壞，森林質(zhì)量明顯下降，森林景觀破碎化趨勢(shì)日益嚴(yán)重，對(duì)長(zhǎng)白山地區(qū)森林生態(tài)系統(tǒng)的穩(wěn)定性造成嚴(yán)重威脅［ 12］。從1998年開始，長(zhǎng)白山重點(diǎn)國(guó)有林區(qū)成為天然林保護(hù)工程的試點(diǎn)區(qū)域［ 13］。因此，筆者在天然林保護(hù)工程實(shí)施的背景下，從公益林和商品林的角度，通過(guò)RS和GIS技術(shù)反映1987—2017年長(zhǎng)白山地區(qū)露水河林業(yè)局景觀格局動(dòng)態(tài)變化過(guò)程，探討森林經(jīng)營(yíng)方式的變化對(duì)森林景觀格局演化的作用，為今后長(zhǎng)白山地區(qū)森林資源的保護(hù)及林業(yè)經(jīng)營(yíng)措施的實(shí)施和改進(jìn)提供科學(xué)依據(jù)。

1 研究區(qū)概況

露水河林業(yè)局地處長(zhǎng)白山脈西北部，全域經(jīng)營(yíng)面積120 934 hm2，包含8個(gè)林場(chǎng)。該林業(yè)局主要從事木材生產(chǎn)、林產(chǎn)品加工等活動(dòng)。20世紀(jì)80年代以前，露水河林業(yè)局主要采取皆伐的經(jīng)營(yíng)方式，曾短暫實(shí)行過(guò)采育兼顧、徑級(jí)伐的經(jīng)營(yíng)方式，之后又重新采用皆伐方式，導(dǎo)致林業(yè)局森林環(huán)境出現(xiàn)退化的危險(xiǎn)。自1998年開始，露水河林業(yè)局作為天然林保護(hù)工程的試點(diǎn)單位，逐漸改變了林業(yè)經(jīng)營(yíng)策略。因此，以露水河林業(yè)局作為研究對(duì)象說(shuō)明林業(yè)經(jīng)營(yíng)政策的改變引起的森林景觀格局變化現(xiàn)象具有代表性，同時(shí)可為該地森林經(jīng)營(yíng)模式的實(shí)施和改進(jìn)提供科學(xué)依據(jù)。

2 研究方法

2.1 影像分類

選取1987年7月、1997、2008、2017年8月4期影像，選用Landsat8OLI傳感器7個(gè)波段；Landsat5TM（專題制圖儀）數(shù)據(jù)（表1）。應(yīng)用多光譜影像數(shù)據(jù)的波段組合顯示，確定解譯標(biāo)志，完成監(jiān)督分類。由于所下載的Landsat-5、Landsat-8數(shù)據(jù)屬于L1TP級(jí)別，其是以地面控制點(diǎn)和DEM數(shù)據(jù)完成幾何精校正，因此該級(jí)別的數(shù)據(jù)無(wú)需再進(jìn)行幾何校正。數(shù)據(jù)預(yù)處理流程包括輻射定標(biāo)、大氣校正、圖像配準(zhǔn)、鑲嵌、裁剪等步驟。

參考露水河二類調(diào)查數(shù)據(jù)，將研究區(qū)域劃分為闊葉林、針葉林、混交林、灌木林、非林地5個(gè)景觀類型。采用ENVI監(jiān)督分類功能（supersied classification）最大似然法（maximum likelihood）進(jìn)行分類。依照景觀類型在遙感影像上反映的紋理、顏色等特點(diǎn)選擇訓(xùn)練樣本，根據(jù)選好的訓(xùn)練樣本在Classification模塊下進(jìn)行監(jiān)督分類。運(yùn)用混淆矩陣對(duì)分類結(jié)果進(jìn)行檢驗(yàn)，經(jīng)過(guò)檢驗(yàn)可知，4期數(shù)據(jù)的總體精度均在98%以上，4期數(shù)據(jù)的Kappa系數(shù)均在0.9以上，符合該研究的精度要求（表2）。采取Majority/Minority功能和聚類處理（clump）功能對(duì)分類結(jié)果進(jìn)行處理。

2.2 景觀格局指數(shù)選取

景觀格局指數(shù)是高度濃縮的景觀格局信息，是反映景觀結(jié)構(gòu)組成和空間配置某些特征方面的簡(jiǎn)單量化指標(biāo)。通過(guò)景觀格局指數(shù)的計(jì)算，可以分析不同景觀類型的分布結(jié)構(gòu)和組成特征。截至目前，用于景觀分析的景觀指數(shù)有很多，而且各景觀指數(shù)直接具有密切的相關(guān)性。因此，可以根據(jù)不同的研究需要和側(cè)重點(diǎn)選擇相應(yīng)的景觀指數(shù)。根據(jù)該研究區(qū)的特點(diǎn)結(jié)合前人研究經(jīng)驗(yàn)，選擇以下景觀指數(shù)進(jìn)行研究［ 14-15］：

3 結(jié)果與分析

3.1 景觀類型面積現(xiàn)狀

該研究區(qū)域景觀基質(zhì)為混交林類型。露水河林業(yè)局1987—2017年景觀類型面積表現(xiàn)為混交林>針葉林>闊葉林>非林地>灌木林（表3，圖1）。該研究區(qū)域中混交林面積最大，是該區(qū)域中各景觀類型的重要組成部分。該區(qū)域的景觀結(jié)構(gòu)屬于典型的基質(zhì)-廊道-斑塊結(jié)構(gòu)，基質(zhì)為混交林，其他景觀斑塊鑲嵌于基質(zhì)中，道路作為廊道連接各斑塊。

3.2 景觀類型變化特征

1987—2017年該研究區(qū)域內(nèi)非林地景觀類型向林地景觀類型轉(zhuǎn)變，林地內(nèi)部景觀類型相互轉(zhuǎn)換，針闊葉林面積呈縮小的趨勢(shì)。1987—2008年，針闊葉林向混交林轉(zhuǎn)變，針闊葉林面積開始縮小，非林地面積開始萎縮并向林地方向轉(zhuǎn)變。1997—2008年，非林地類型向林地類型轉(zhuǎn)移面積有所減小，林地內(nèi)部變化程度相比10年前也有所變化，混交林向針闊葉林轉(zhuǎn)移面積有所減小，針闊葉林向混交林轉(zhuǎn)移面積減小，使得針闊葉林面積小幅回升。2008—2017年非林地類型轉(zhuǎn)移林地面積為近30年最大值；林地內(nèi)部混交林向針闊葉林轉(zhuǎn)移的面積為30年來(lái)最低，針闊葉林面積呈繼續(xù)萎縮的變化趨勢(shì)。灌木林面積在近30年來(lái)變化不明顯，其他景觀類型向灌木林轉(zhuǎn)移的趨勢(shì)亦不明顯（表4，圖1）。

3.3 景觀格局變化特征

3.3.1 景觀指數(shù)變化。以林業(yè)局作為單位來(lái)看（表5，圖2），1987—1997年NP、PD呈增長(zhǎng)趨勢(shì)，說(shuō)明在前10年林業(yè)局景觀破碎化程度加重；1997—2017年NP、PD值均呈下降趨勢(shì)，說(shuō)明在后20年該地區(qū)景觀破碎化情況得到改善。1987—2017年ED呈現(xiàn)逐漸上升后下降趨勢(shì)，說(shuō)明在這30年間該地區(qū)斑塊形狀變化程度呈先加重后緩解的趨勢(shì)，PLADJ指數(shù)和COHESION指數(shù)變化不大且始終維持高水平，說(shuō)明該地區(qū)斑塊連通性較好。SHDI指數(shù)較高，說(shuō)明該地區(qū)景觀多樣性水平較低，但斑塊整體分布較為均勻。

3.3.2 不同經(jīng)營(yíng)區(qū)劃內(nèi)景觀格局變化特征。以林業(yè)局內(nèi)部經(jīng)營(yíng)區(qū)作為單位進(jìn)行分析（表6，圖2），商品林區(qū)、公益林區(qū)NP、PD指數(shù)均呈先上升后下降的變化趨勢(shì)。說(shuō)明林業(yè)局各經(jīng)營(yíng)區(qū)內(nèi)部景觀破碎化程度趨于緩解。ED指數(shù)與NP、PD指數(shù)變化趨勢(shì)相同，說(shuō)明各經(jīng)營(yíng)區(qū)斑塊形狀逐漸趨于穩(wěn)定。各經(jīng)營(yíng)區(qū)PLADJ指數(shù)在2008—2017年呈上升趨勢(shì)，COHESION指數(shù)變化趨勢(shì)不明顯且指數(shù)較高，說(shuō)明斑塊連通性得到改善。商品林區(qū)、公益林區(qū)SHDI指數(shù)較低，說(shuō)明這2個(gè)經(jīng)營(yíng)區(qū)景觀多樣性水平較低；SHEI指數(shù)較高，說(shuō)明各經(jīng)營(yíng)區(qū)內(nèi)景觀類型分布較均勻。

4 討論

4.1 景觀類型面積轉(zhuǎn)移分析

經(jīng)營(yíng)模式的變化影響森林景觀類型面積的分布。該研究結(jié)果表明，露水河林業(yè)局內(nèi)占主要成分的景觀類型為混交林、針葉林、闊葉林，這與長(zhǎng)白山自然保護(hù)區(qū)內(nèi)整體景觀面積分布情況相同［ 16］。混交林、針葉林、闊葉林在近30年共同組成為該區(qū)域的復(fù)合景觀，它們對(duì)該研究區(qū)域的生態(tài)環(huán)境具有調(diào)控作用，并對(duì)當(dāng)?shù)厣鷳B(tài)環(huán)境演變施加影響［ 17］，主導(dǎo)著當(dāng)?shù)鼐坝^整體上的結(jié)構(gòu)和動(dòng)態(tài)變化［ 18］。林地面積基本保持不變，林地內(nèi)部之間轉(zhuǎn)換，主要為針葉林，闊葉林向混交林轉(zhuǎn)移的變化趨勢(shì)。這一變化趨勢(shì)與天然林保護(hù)工程的經(jīng)營(yíng)策略方向相符合，該經(jīng)營(yíng)策略要求撫育工作要做到人工林天然化，在進(jìn)行人工林培育時(shí)可以進(jìn)行帶狀或團(tuán)狀的培育，向針闊混交的方向進(jìn)行培育［ 19］。說(shuō)明目前的景觀類型分布格局基本符合天然林保護(hù)工程的要求。2008年之前針葉林、闊葉林向混交林轉(zhuǎn)移的速度相比2008年之后較為緩慢，這與森林經(jīng)營(yíng)模式改革逐漸深化過(guò)程密切相關(guān)。非林地面積逐漸減小，其主要原因?yàn)殡S著天然林保護(hù)工程的實(shí)施，林業(yè)局對(duì)經(jīng)營(yíng)模式作出了調(diào)整，改變育林思路，加大了育林力度等措施［ 20］。說(shuō)明該林業(yè)局實(shí)施天保工程的效果較為顯著。通過(guò)以上分析可以明確該區(qū)域森林景觀類型分布情況的變化與經(jīng)營(yíng)模式的轉(zhuǎn)變密切相關(guān)。

4.2 林業(yè)局景觀格局變化分析

從整體上看，經(jīng)營(yíng)方式的轉(zhuǎn)變是導(dǎo)致林業(yè)局景觀格局變化的主要因素。1987—1997年該地區(qū)景觀破碎化，斑塊程度愈加嚴(yán)重，主要原因是在90年代期間該林業(yè)局仍采取傳統(tǒng)的皆伐方式，形成采伐比重過(guò)大，忽略林業(yè)撫育的局面［ 21-22］。同一時(shí)期該地區(qū)開展了道路修建等經(jīng)濟(jì)生產(chǎn)活動(dòng)［ 23-24］，導(dǎo)致該區(qū)域的森林景觀結(jié)構(gòu)整體呈現(xiàn)破碎化的趨勢(shì)。1997—2017年景觀破碎化現(xiàn)象逐漸緩解。主要原因是隨著天然林保護(hù)工程的展開，林業(yè)局改變傳統(tǒng)的經(jīng)營(yíng)方式，采取分類經(jīng)營(yíng)模式。依照分類經(jīng)營(yíng)的理念，林業(yè)局管轄區(qū)域內(nèi)自然保護(hù)區(qū)、水源地、生態(tài)環(huán)境較為脆弱的林地以及立地條件較差，交通管理不便的林地劃分為公益林區(qū)；林業(yè)局管轄區(qū)域內(nèi)立地條件較好，交通較為便利的林地劃分為商品林［ 19］。經(jīng)過(guò)以上經(jīng)營(yíng)策略的實(shí)施，林業(yè)局內(nèi)景觀破碎化程度逐漸得到遏制，但出現(xiàn)了景觀多樣性水平較低且變化不明顯的現(xiàn)象。根據(jù)管東生等［ 25-26］研究，多樣性水平變化不明顯且水平較低的原因可能是經(jīng)營(yíng)區(qū)內(nèi)景觀類型面積比例沒(méi)有發(fā)生太大變化，所占面積較大的斑塊未發(fā)生較大變化，林場(chǎng)內(nèi)部景觀類型面積所占比例差異懸殊導(dǎo)致的。這與該研究關(guān)于景觀類型面積的研究結(jié)果相一致。結(jié)合以上分析可以有理由認(rèn)為，該林業(yè)局轉(zhuǎn)變經(jīng)營(yíng)方式對(duì)維持該地區(qū)森林生態(tài)環(huán)境起著積極作用。

4.3 經(jīng)營(yíng)區(qū)景觀格局變化分析

經(jīng)過(guò)對(duì)該林業(yè)局實(shí)施天保工程前后森林景觀格局變化進(jìn)行分析，對(duì)商品林區(qū)、公益林區(qū)內(nèi)部景觀格局變化情況進(jìn)行了分析。商品林區(qū)、公益林區(qū)在1987—1997年破碎化程度較高，斑塊形狀較復(fù)雜，經(jīng)過(guò)前人研究表明，該研究區(qū)1985—1999年森林景觀呈破碎化趨勢(shì)［ 27］，這與該研究結(jié)果相符，說(shuō)明在前10年經(jīng)營(yíng)區(qū)受外部干擾較為嚴(yán)重，由于1987—1997年該區(qū)域尚未完全普及分類經(jīng)營(yíng)制度，部分區(qū)域仍實(shí)行傳統(tǒng)的經(jīng)營(yíng)方式，即未節(jié)制地砍伐活動(dòng)，導(dǎo)致商品林、公益林區(qū)域破碎化嚴(yán)重；1997—2008年商品林區(qū)破碎化趨勢(shì)受到遏制，這種現(xiàn)象與天保工程的實(shí)施密切相關(guān)，隨著天保工程的實(shí)施及禁伐區(qū)的設(shè)立，同時(shí)放棄了皆伐的生產(chǎn)方式，導(dǎo)致商品林區(qū)采伐量減少，進(jìn)而使得過(guò)去破碎化嚴(yán)重的現(xiàn)象得到緩解。2008—2017年商品林區(qū)景觀破碎化有再度發(fā)展的趨勢(shì)，其主要原因是在天保工程實(shí)施的背景下，商品林區(qū)的采伐量大幅下降，與此同時(shí)商品林區(qū)的蓄積量也大幅下降，因此從總體上看商品林區(qū)的相對(duì)采伐強(qiáng)度增大，進(jìn)而導(dǎo)致商品林區(qū)的破碎化趨勢(shì)有增強(qiáng)趨勢(shì)［ 28］。1997—2017年公益林區(qū)破碎化趨勢(shì)逐漸得到遏制。這種現(xiàn)象與天保工程的實(shí)施密切相關(guān)。針對(duì)公益林區(qū)的撫育工作，該林業(yè)局遵守嚴(yán)格的撫育次數(shù)和年限；針對(duì)該區(qū)域中重點(diǎn)公益林部分則嚴(yán)格采取禁伐措施；針對(duì)該區(qū)域中生態(tài)公益林部分則采取了適當(dāng)采伐的方式且采伐目的是以改善森林系統(tǒng)為主導(dǎo)。經(jīng)過(guò)以上的舉措，公益林區(qū)的破碎化趨勢(shì)得到很大改善。1987—2017年商品林區(qū)、公益林區(qū)斑連通性始終保持較高水平，說(shuō)明商品林區(qū)和公益林區(qū)的森林資源環(huán)境并未受到大規(guī)模的破壞和退化，這與之前得出的破碎化程度先嚴(yán)重后逐漸趨緩的結(jié)論相矛盾。這可能是由于經(jīng)營(yíng)區(qū)內(nèi)的景觀破碎化現(xiàn)象只是發(fā)生在小范圍區(qū)域內(nèi)，并未威脅到林地景觀斑塊整體的聚集程度，說(shuō)明經(jīng)營(yíng)區(qū)內(nèi)森林環(huán)境目前沒(méi)有發(fā)生大規(guī)模的退化。商品林區(qū)和公益林區(qū)景觀多樣性水平較低，景觀分布均勻且隨時(shí)間變化不明顯。目前露水河林業(yè)局中商品林和公益林分別占森林面積的64.3%和35.7%，在這種分布情況下商品林區(qū)和公益林區(qū)的SHDI指數(shù)仍維持在較低水平，說(shuō)明該區(qū)域中優(yōu)勢(shì)景觀類型所占面積沒(méi)有發(fā)生太大變化，由此可知各經(jīng)營(yíng)區(qū)內(nèi)景觀格局水平?jīng)]有發(fā)生大規(guī)模的退化，同時(shí)各景觀之間的分布較為均勻，說(shuō)明維持目前的經(jīng)營(yíng)區(qū)劃對(duì)于維持森林生態(tài)環(huán)境穩(wěn)定具有積極作用［ 28］。

5 結(jié)論

該研究結(jié)果表明，露水河林業(yè)局中混交林為優(yōu)勢(shì)景觀， 其他類型斑塊鑲嵌其中，整體呈現(xiàn)典型的基質(zhì)—廊道—斑塊景觀類型格局。非林地類型向林地類型轉(zhuǎn)變，林地面積基本保持不變。實(shí)施天保工程后，林業(yè)局整體景觀破碎化程度趨緩，但破碎化現(xiàn)象依然嚴(yán)峻，斑塊連通性較好，景觀多樣性仍保持較低水平。商品林區(qū)景觀破碎化現(xiàn)象明顯，但已經(jīng)出現(xiàn)好轉(zhuǎn)的跡象，公益林區(qū)景觀破碎化程度減弱 。在實(shí)施天保工程后露水河林業(yè)局森林生態(tài)環(huán)境整體上呈現(xiàn)逐漸改善的趨勢(shì)。未來(lái)10年內(nèi)露水河林業(yè)局森林生態(tài)環(huán)境呈現(xiàn)繼續(xù)向好的發(fā)展趨勢(shì)，因此目前的經(jīng)營(yíng)方式對(duì)于維持露水河林業(yè)局森林生態(tài)環(huán)境具有積極作用。

參考文獻(xiàn)

［1］楊國(guó)靖，肖篤寧.中祁連山淺山區(qū)山地森林景觀空間格局分析［J］.應(yīng)用生態(tài)學(xué)報(bào)，2004，15（2）：269-272.

［2］ 肖德榮，田昆，袁華，等.滇西北高原典型退化濕地納帕海植物群落景觀多樣性［J］.生態(tài)學(xué)雜志，2007，26（8）：1171-1176.

［3］ 公霞，常禹，布仁倉(cāng)，等.呼中林業(yè)局森林采伐方式對(duì)森林景觀格局的長(zhǎng)期影響［J］.生態(tài)學(xué)雜志，2006，25（7）：805-812.

［4］ 賈寶全，慈龍駿，任一萍.綠洲景觀動(dòng)態(tài)變化分析［J］.生態(tài)學(xué)報(bào)，2001，21（11）：1947-1951.

［5］ 劉吉平，趙丹丹，田學(xué)智，等.1954—2010年三江平原土地利用景觀格局動(dòng)態(tài)變化及驅(qū)動(dòng)力［J］.生態(tài)學(xué)報(bào)，2014，34（12）：3234-3244.

［6］ 李月輝，胡遠(yuǎn)滿，常禹，等.大興安嶺呼中林業(yè)局森林景觀格局變化及其驅(qū)動(dòng)力［J］.生態(tài)學(xué)報(bào)，2006，26（10）：3347-3357.

［7］ 廖芳均，趙東升.南嶺國(guó)家級(jí)自然保護(hù)區(qū)森林景觀格局變化與動(dòng)態(tài)模擬［J］.地理科學(xué)，2014，34（9）：1099-1107.

［8］ 肖篤寧，趙羿，孫中偉，等.沈陽(yáng)西郊景觀格局變化的研究［J］.應(yīng)用生態(tài)學(xué)報(bào)，1990，1（1）：75-84.

［9］ 劉紅玉，呂憲國(guó)，張世奎，等.三江平原流域濕地景觀破碎化過(guò)程研究［J］.應(yīng)用生態(tài)學(xué)報(bào)，2005，16（2）：289-295.

［10］ 馬克明，傅伯杰，周華峰.景觀多樣性測(cè)度：格局多樣性的親和度分析［J］.生態(tài)學(xué)報(bào)，1998，18（1）：76-81.

［11］ 李月輝，常禹，胡遠(yuǎn)滿，等.人類活動(dòng)對(duì)森林景觀影響研究進(jìn)展［J］.林業(yè)科學(xué)，2006，42（9）：119-126.

［12］ 李國(guó)偉，趙偉，魏亞偉，等.天然林資源保護(hù)工程對(duì)長(zhǎng)白山林區(qū)森林生態(tài)系統(tǒng)服務(wù)功能的影響［J］.生態(tài)學(xué)報(bào)，2015，35（4）：984-992.

［13］ 張佩昌.試論天然林保護(hù)工程［J］.林業(yè)科學(xué)，1999，35（2）：124-131.

［14］ 高照忠，魏海霞，黃鐵蘭.粵港澳大灣區(qū)土地覆蓋及景觀格局時(shí)空變化分析［J］.測(cè)繪通報(bào)，2021（5）：25-29.

［15］ 白少云，劉斌，余珮珩，等.融合多尺度影像數(shù)據(jù)的杞麓湖流域景觀格局分析［J］.測(cè)繪通報(bào)，2020（2）：90-95.

［16］ 王美玲，張繼超，王舶鑒，等.長(zhǎng)白山區(qū)森林景觀格局動(dòng)態(tài)［J］.生態(tài)學(xué)雜志，2017，36（11）：3138-3147.

［17］ 趙福強(qiáng)，代力民，于大炮，等.長(zhǎng)白山露水河林業(yè)局森林景觀格局動(dòng)態(tài)［J］.應(yīng)用生態(tài)學(xué)報(bào)，2010，21（5）：1180-1184.

［18］ 岳剛，楊華，亢新剛，等.基于GIS的長(zhǎng)白山森林景觀格局分析［J］.中南林業(yè)科技大學(xué)學(xué)報(bào)，2013，33（7）：35-39.

［19］ 羅菊春.論我國(guó)天然林保護(hù)的技術(shù)措施及其基本原理［J］.北京林業(yè)大學(xué)學(xué)報(bào)，2000，22（4）：128-129.

［20］ 王志新.實(shí)施森林可持續(xù)經(jīng)營(yíng)是林業(yè)發(fā)展的必然選擇：關(guān)于吉林省森林經(jīng)營(yíng)工作的調(diào)查與思考［J］.吉林林業(yè)科技，2013，42（1）：30-36.

［21］ 李秀娟. 吉林省國(guó)有林區(qū)經(jīng)濟(jì)社會(huì)環(huán)境系統(tǒng)協(xié)調(diào)發(fā)展評(píng)價(jià)研究［D］.北京：北京林業(yè)大學(xué)，2008.

［22］ 周瑜，岳彩榮.香格里拉縣森林景觀格局分析［J］.林業(yè)建設(shè)，2009（6）：25-29.

［23］ 段春霞，胡遠(yuǎn)滿，李月輝，等.大興安嶺北部林區(qū)景觀格局變化及其影響分析［J］.生態(tài)學(xué)雜志，2004，23（2）：133-135.

［24］? 余世孝，郭濼.山東泰山地區(qū)景觀結(jié)構(gòu)變化及其分形分析［J］.生態(tài)學(xué)報(bào)，2005，25（1）：129-134.

［25］ 管東生，鐘曉燕，鄭淑穎.廣州地區(qū)森林景觀多樣性分析［J］.生態(tài)學(xué)雜志，2001，20（4）：9-12.

［26］ 于德永，王艷艷，郝占慶，等.吉林省露水河地區(qū)森林景觀格局變化［J］.資源科學(xué)，2005，27（4）：147-153.

［27］ 于大炮，周旺明，包也，等.天保工程實(shí)施以來(lái)東北闊葉紅松林的可持續(xù)經(jīng)營(yíng)［J］.生態(tài)學(xué)報(bào)，2015，35（1）：10-17.

［28］ 侯瑩. 吉林省森林分類經(jīng)營(yíng)的研究［D］.長(zhǎng)春：東北師范大學(xué)，2008.安徽農(nóng)業(yè)科學(xué)，J.Anhui Agric.Sci. 2023，51（14）：168-177