摘 要：" 為探討地形對(duì)林下灌草層植物生物量的影響，該研究采用嵌套設(shè)計(jì)法在東北虎豹國(guó)家公園調(diào)查了138個(gè)密林下樣地共1 685個(gè)植物樣方，通過嵌套方差分析與有序邏輯斯蒂回歸模型對(duì)林下灌草層植物生物量受地形的影響進(jìn)行了分析。結(jié)果表明：（1）不同坡位之間，谷底的灌草層植物生物量高于坡上，坡上高于坡下（Plt;0.01）；不同坡向之間，陰坡灌草層植物生物量低于陽坡及平地（Plt;0.01），后二者間無顯著差異；不同坡度之間，平坡灌草層植物生物量高于陡坡，陡坡高于緩坡（Plt;0.01）。（2）坡位與坡向的交互作用顯著，坡下平地、坡上平地、坡上陽坡與谷底的所有坡位灌草層植物生物量最高，坡下陰坡、坡下陽坡及坡上陰坡之間無顯著差異。（3）研究區(qū)現(xiàn)行狀態(tài)下，有序邏輯斯蒂回歸結(jié)果顯示，灌草層植物生物量在不同海拔、坡位及坡向坡度組合下不同。坡位、坡向及坡度對(duì)林下灌草層植物生物量有顯著影響，3個(gè)坡位等級(jí)間谷底最高而坡下最低，3個(gè)坡度等級(jí)間陡坡最高而緩坡最低，不同坡向比較，陰坡最低。（4）在不排除人為干擾、森林放牧的現(xiàn)實(shí)情況下，谷底、陡坡地帶灌草層植物生物量概率最高。該研究結(jié)果可為準(zhǔn)確估計(jì)東北虎豹國(guó)家公園林下灌草層植物對(duì)虎豹獵物種群的承載力提供重要參考，從而為瀕?；⒈谋Ｗo(hù)和管理提供科學(xué)依據(jù)。

關(guān)鍵詞： 林下灌草層植物， 有蹄類動(dòng)物食物資源， 地形， 東北虎豹國(guó)家公園， 有序邏輯斯蒂回歸

中圖分類號(hào)：" Q948

文獻(xiàn)標(biāo)識(shí)碼：" A

文章編號(hào)：" 1000-3142（2024）08-1512-12

Impacts of topography on spatial variability of temperate understory biomass in closed forest： a case study of Northeast China Tiger and Leopard National Park

WANG Le1，2， MOU Pu2，3， WANG Tianming2，3

（ 1. Institute of Ecological Protection and Restoration， Chinese Academy of Forestry， Beijing 100091， China; 2. State Forestry and GrasslandAdministration Key Laboratory for Conservation Ecology of Northeast China Tiger and Leopard National Park， Beijing100875， China; 3. College of Life Sciences， Beijing Normal University， Beijing 100875， China ）

Abstract：" To investigate the impact of topography on" understory biomass， this study employed a nested design and conducted surveys in the Northeast China Tiger and Leopard National Park， involving a total of 138 plots within closed forests， comprising 1 685 plant quadrats. The analysis of topographical effect on understory biomass was carried out using nested analysis of variance and an ordinal logistic regression model. The results were as follows： （1） Among different slope positions， the understory biomass" in valleys was higher than on upper slopes， and the footslopes where the understory biomass was the lowest （Plt;0.01）. The understory biomass" on shady slopes was lower than on sunny slopes and flats （Plt;0.01）， with no significant difference between the latter two. In terms of different slopes， the understory biomass on steep slopes was higher than on flats， and flats had higher biomass than gentle slopes （Plt;0.01）. （2） There was a significant interaction between slope position and slope aspect. The understory biomass" was the highest on flats at the footslopes， flats on upper slopes， upper slopes and sunny slopes， and in valleys among all slope positions. There was no significant difference in understory biomass among shady slopes at the footslopes， sunny slopes at the footslopes， and shady slopes on upper slopes. （3） Ordinal logistic regression analysis of the current state of the study area showed that the probability of high or low understory biomass varied across different combinations of altitude， slope position， and slope aspect. Slope position， slope aspect， and slope had significant effects on the understory biomass. Among the three slope positions， the understory biomass was the highest in valleys and the lowest on footslopes. Among the three slopes， steep slopes had the highest understory biomass， while gentle slopes had the lowest. Shady slopes had the lowest understory biomass compared to other slope aspects. （4） Considering the realistic situation of human disturbance and forest grazing， steep slope areas in valleys had the highest probability of high understory biomass. This study can provide important references for accurately estimating the carrying capacity of understory for preys populations of tiger and leopard in Northeast China Tiger and Leopard National Park， thereby offering a scientific basis for the conservation and management of endangered tigers and leopards.

Key words： understory， forage resource of ungulates， topography， Northeast China Tiger and Leopard National Park， ordinal logistic regression

在東北虎豹國(guó)家公園及鄰近地區(qū)，大中型有蹄類動(dòng)物如梅花鹿（Cervus nippon）、狍（Capreolus pygargus）和野豬（Sus scrofa）等是東北虎（Panthera tigris altaica）和東北豹（P. pardus orientalis）的主要獵物（Kerley et al.， 2015; Yang et al.， 2018; 王天明等， 2020），這些有蹄類動(dòng)物主要以林下灌草層植物為食（Massei et al.， 1996; Baskin amp; Danell， 2003; Gordon amp; Prins， 2008; Barancˇeková et al.， 2010; 王樂等， 2019a），因此林下灌草層植物生物量是影響生態(tài)系統(tǒng)食草動(dòng)物承載力的關(guān)鍵因素之一（Whittaker amp; Likens， 1973; Gilliam amp; Roberts， 2003）。早春時(shí)節(jié)，由于林冠樹葉稀疏，林下光照充裕，灌草層植物在較強(qiáng)的光照條件下迅速生長(zhǎng)，但隨著喬木冠層的逐漸郁閉，林下光照逐漸減少（王淼等， 2006），特別是林冠郁閉后，灌草層生物量的生長(zhǎng)開始減緩（Jolly et al.， 2004），此時(shí)地形因素對(duì)林下灌草層生產(chǎn)力的相對(duì)重要性明顯上升（Gracia et al.， 2007）。

復(fù)雜的地形產(chǎn)生了各種各樣的微環(huán)境，影響著局地微氣候和土壤，從而對(duì)局部生境中的光、溫、水、肥等植物生長(zhǎng)關(guān)鍵因素產(chǎn)生較大影響（龔艷賓， 2016; 申文輝等， 2017; 李芹等， 2019; 孫莉英等， 2020; 李娟和龔純偉， 2021）。地形因子如海拔、坡位、坡度、坡向，通過調(diào)節(jié)光照、溫度、水分及土壤中養(yǎng)分等生態(tài)因子來影響植物的生長(zhǎng)（李俊清， 2010；秦隨濤等， 2019），形成了灌草層植物資源的空間異質(zhì)分布格局（龔艷賓， 2016）。在東北虎豹國(guó)家公園的鄰近地區(qū)，地形對(duì)森林生物量的影響研究較少。在長(zhǎng)白山金溝嶺林場(chǎng)，不同海拔與坡度的森林地上生物量差異顯著，而不同坡向差異不大（姜慧泉， 2011）。海拔、坡度與坡向?qū)Υ笈d安嶺森林生物量空間異質(zhì)性影響顯著（黎如， 2010）。坡位、坡向是影響小興安嶺喬、灌層生物量的重要地形因子，生物量隨坡位上升呈先增加后減少的趨勢(shì)；坡向則隨陰坡、半陰坡、半陽坡、陽坡、平地有遞減趨勢(shì)（王曉莉， 2013）。黑龍江帽兒山的森林生物量與海拔和坡度呈顯著的相關(guān)性，與坡向則無顯著相關(guān)性（王維芳等， 2010）。上述研究結(jié)果說明海拔、坡位、坡向與坡度等地形因子對(duì)局地森林生物量的影響非常復(fù)雜。雖然局部生境的異質(zhì)性導(dǎo)致的森林生物量不同的研究已有許多（Gordon amp; Prins， 2008），但對(duì)林下灌草層植物生物量影響的研究，國(guó)內(nèi)外相關(guān)報(bào)道較少（de Castilho et al.， 2006; Nie et al.， 2019）。然而，探討地形對(duì)林下灌草層植物生物量的影響，是掌握灌草層生物量分異規(guī)律的基礎(chǔ)，也是研究生態(tài)系統(tǒng)食草動(dòng)物承載力的重要基礎(chǔ)。

本研究以東北虎豹國(guó)家公園為研究區(qū)域，針對(duì)園區(qū)內(nèi)大中型有蹄類動(dòng)物的主要食物來源——密林下灌草層植物，采用嵌套方差分析與有序多分類邏輯斯蒂回歸的方法，通過海拔、坡位、坡向坡度對(duì)研究區(qū)密林下灌草層植物生物量進(jìn)行分析與計(jì)算，擬探討以下問題：（1）地形因子對(duì)密林下灌草層植物生物量分異規(guī)律的影響；（2）預(yù)測(cè)高產(chǎn)灌草層植物生物量區(qū)域的地形條件。旨在掌握灌草層植物生物量在不同地形下的分異規(guī)律，并創(chuàng)新性地使用地形對(duì)高產(chǎn)灌草層植物生物量區(qū)域進(jìn)行預(yù)測(cè)。本研究結(jié)果可為國(guó)家公園中有蹄類動(dòng)物承載力的準(zhǔn)確估算提供一手?jǐn)?shù)據(jù)，進(jìn)而為公園內(nèi)虎豹種群的管理與保護(hù)提供依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)位于東北虎豹國(guó)家公園（129°10′—131°20′ E、42°30′—44°10′ N）的東部，是東北虎、東北豹在中國(guó)的主要分布區(qū)和優(yōu)先保護(hù)區(qū)（Hebblewhite et al.， 2014; Wang et al.， 2016; Feng et al.， 2017）。研究區(qū)屬于半濕潤(rùn)大陸性季風(fēng)氣候，1月均溫 -15 ℃，7月均溫20 ℃，年降水（30年記錄）365.0～842.9 mm，雨熱同季，冬季積雪較深。研究區(qū)地貌為山地丘陵，海拔2～1 486 m，主要河流有琿春河、圖們江、嘎呀河及綏芬河。山岳丘陵多形成各種地形，沿琿春河谷地，地勢(shì)平坦，海拔低，山勢(shì)也緩和，陡峭險(xiǎn)峻地勢(shì)少，而入深山區(qū)，海拔漸高，山勢(shì)亦漸相對(duì)高聳陡峭。該區(qū)地帶性土壤為暗棕壤，地帶性植被為闊葉紅松（Pinus koraiensis）林，由于人們的經(jīng)濟(jì)活動(dòng)和森林采伐，現(xiàn)已多被闊葉混交林、蒙古櫟（Quercus mongolica）林及山楊（Populus davidiana）林等不同演替階段的次生森林所取代（周以良和李景文， 1964）。林下灌木層植物以杜鵑屬（Rhododendron）、山梅花屬（Philadelphus）、溲疏屬（Deutzia）、五味子屬（Schisandra）、懸鉤子屬（Rubus）、榛屬（Corylus）和忍冬屬（Lonicera）為主，草本植物為薹草屬（Carex）、臭草屬（Melica）、委陵菜屬（Potentilla）、唐松草屬（Thalictrum）、落新婦屬（Astilbe）、堇菜屬（Viola）、茜草屬（Rubia）、蚊子草屬（Filipendula）、香茶菜屬（Isodon）、蟹甲草屬（Parasenecio）、野豌豆屬（Vicia）、蒿屬（Artemisia）和蕨屬（Pteridium）等植物（Wang et al.， 2023）。區(qū)內(nèi)分布有東北虎、東北豹、梅花鹿、狍和野豬等動(dòng)物（肖文宏等， 2014; Zhang et al.， 2014; Wang et al.， 2016）。有蹄類動(dòng)物主要取食的植物有青楷槭（Acer tegmentosum）、臭冷杉（Abies nephrolepis）、紫椴（Tilia amurensis）、白樺（Betula platyphylla）、五角槭（Acer pictum subsp. mono）、毛榛（Corylus mandshurica）、珍珠梅（Sorbaria sorbifolia）、粗莖鱗毛蕨（Dryopteris crassirhizoma）、薹草（Carex spp.）等（王祥生， 2019）。

1.2 數(shù)據(jù)采集

于2015—2016年植物生長(zhǎng)旺季（6—8月）隨機(jī)選擇了140個(gè)樣地進(jìn)行野外植被調(diào)查。在每個(gè)樣地中心，向北、東、西、南4個(gè)方向各拉50 m長(zhǎng)樣線。每條樣線上選取4個(gè)隨機(jī)點(diǎn)，每個(gè)樣點(diǎn)隨機(jī)選取1個(gè)1 m × 1 m的樣方，對(duì)樣方內(nèi)的草本植物（不包括苔蘚）的地上部分及高度2.2 m以下的喬木及灌木當(dāng)年生嫩枝葉的生物量進(jìn)行收獲并裝入紙袋，風(fēng)干后快遞寄回實(shí)驗(yàn)室。

在每個(gè)樣地內(nèi)，使用GPS（GARMIN， rino 530HCx， 美國(guó)）記錄樣地中心點(diǎn)的海拔，之后用羅盤測(cè)量并記錄每條樣線所在位置的坡向（0°～360°），按陰坡（坡向0°～135°和315°～360°）、陽坡（135°～315°）（張建亮等， 2016）和平地3類進(jìn)行記錄。每個(gè)樣地以500 m為界分為高海拔、低海拔2類。每條樣線上的4個(gè)樣方，以坡上、坡中（分析時(shí)歸入坡上）、坡下及谷底記錄坡位，并按平坡（坡度≤5°）、緩坡（5°lt;坡度≤25°）、陡坡（坡度gt;25°）3類測(cè)量并記錄樣方的坡度。由于懸崖等危險(xiǎn)地形造成樣線和樣方的缺失，因此共調(diào)查樣線496條、樣方1 948個(gè)，在樣方中用冠層分析儀（WinScanopy）測(cè)量林冠層郁閉度，生境開闊區(qū)域（郁閉度低于0.7）如林窗、林緣及河邊旁樣方（263個(gè)）增加了光照或水分的影響，不在本研究的考慮范圍內(nèi)，最終確定郁閉度在0.7～0.9之間的1 685個(gè)密林下植物樣方作為本研究的對(duì)象，它們歸屬于138個(gè)樣地（王樂等， 2019a）。

1.3 數(shù)據(jù)分析

1.3.1 嵌套方差分析 雖然林下1 685個(gè)樣方基本獨(dú)立，并且其環(huán)境信息基本準(zhǔn)確到樣方水平，但是考慮到樣方的環(huán)境因子多嵌套在138個(gè)樣地之上，為了防止可能出現(xiàn)的空間自相關(guān)導(dǎo)致的自由度過高，以海拔、坡位、坡向及坡度4個(gè)分類變量為自變量，將樣地作為嵌套因子（方差分析主要側(cè)重于地形因素的影響），分別以林下灌草層所有植物資源的總生物量為因變量，使用線性混合模型進(jìn)行嵌套多因素單變量（univariate）方差分析，并采用Tukey檢驗(yàn)對(duì)方差分析結(jié)果進(jìn)行后檢驗(yàn)。在數(shù)據(jù)分析之前，先對(duì)不同調(diào)查時(shí)間（年份、月份）之間林下灌草層所有植物總生物量進(jìn)行方差分析，結(jié)果顯示不存在顯著差異，可進(jìn)行下一步分析。上述所有統(tǒng)計(jì)計(jì)算過程采用R語言（Ver. 4.2.2）中的multcomp、lattice、lme4和lsmeans程序包進(jìn)行（R Core Team， 2021）。

1.3.2 有序邏輯斯蒂回歸計(jì)算過程

由于一些變量分類多于二類，采用有序邏輯斯蒂回歸對(duì)不同地形類組合的灌草層植物生物量進(jìn)行概率預(yù)測(cè)。有序多分類回歸模型依次將響應(yīng)變量按不同的取值水平分割成3個(gè)等級(jí)，對(duì)這3個(gè)等級(jí)建立的反應(yīng)變量為三分類的邏輯斯蒂回歸模型（張文彤， 2002）。計(jì)算公式如下（公式1～公式7）：

發(fā)生比odds = P/（1-P）（1）

兩邊取對(duì)數(shù) ln（odds） = ln ［ P/（1-P）］（2）

有序邏輯斯蒂回歸的連接函數(shù)logit （P1） = α1- （β1x1+…+βnxn） = ln ［ P1/（1-P1）］（3）

則logit （P1+P2） = α2- （β1x1+…+βnxn）=ln ［ （P1+P2）/（1-P1-P2）］=ln ［ （P1+P2）/P3 ］（4）

故P1 = elogit （P1）/ （1+ elogit （P1））= exp ［α1- （β1x1+…+βnxn）］/ ［ 1+ exp （α1- （β1x1+…+βnxn）） ］（5）

P2 = （P1+P2）-P1 = exp ［α2- （β1x1+…+βnxn）］/ ［1+ exp （α2- （β1x1+…+βnxn）） ］ -exp ［α1- （β1x1+…+βnxn）］/ ［ 1+ exp （α1- （β1x1+…+βnxn）） ］（6）

P3 = 1- （P1+P2） =1- exp ［α2- （β1x1+…+βnxn）］/ ［1+ exp （α2- （β1x1+…+βnxn）） ］（7）

式中： P表示事件發(fā)生的概率（P的取值范圍為0～1，P1、P2、P3分別表示三分類反應(yīng)變量發(fā)生的概率），x表示自變量（x1…xn分別表示1…n個(gè)自變量），α表示常數(shù)，β表示自變量的系數(shù)。

1.3.3 有序邏輯斯蒂回歸分析 為了對(duì)研究區(qū)林下灌草層植物生物量進(jìn)行預(yù)測(cè)，本研究選擇方差分析中使用的自變量，進(jìn)行變量組合。由于現(xiàn)實(shí)中有些組合不存在，有些組合可能理論上存在，但在研究區(qū)內(nèi)不存在，或樣本量極低，如平地緩坡、平地陡坡等，谷底的陡坡理論、現(xiàn)實(shí)上均可存在，但數(shù)據(jù)中僅有1個(gè)樣方，因此在計(jì)算時(shí)將坡向與坡度合并為陰坡緩坡、陽坡緩坡、陰坡陡坡、陽坡陡坡、平坡5類，將谷底陡坡（陰陽坡）與谷底平坡并入1類。以上分類變量，在進(jìn)行邏輯斯蒂回歸之前，轉(zhuǎn)換成數(shù)字啞變量（表1）。

將密林下灌草層植物生物量分為3個(gè)等級(jí)［低（生物量≤20 g·m-2）、中（20 g·m-2lt;生物量≤40 g·m-2）、高（生物量gt; 40 g·m-2）］。使用有序多分類邏輯斯蒂回歸，以海拔、坡位、坡向坡度這3個(gè)分類自變量對(duì)研究區(qū)林下灌草層植物生物量（分3級(jí)）進(jìn)行概率計(jì)算。計(jì)算使用IBM SPSS Statistics 20軟件（v20.0， SPSS Company， Chicago， USA）進(jìn)行。

2 結(jié)果與分析

2.1 不同地形密林下灌草層植物群落物種組成

根據(jù)調(diào)查，灌草層片除了44種常見針、闊葉樹種的幼樹幼苗外，還包括23種灌木、7種木質(zhì)藤本、3種禾草、5種莎草、83種雜類草及7種蕨類植物，共172種分屬于58科121屬（表2）。灌木以薔薇科（Rosaceae）和忍冬科（Caprifoliaceae）植物為主，草本植物以百合科（Liliaceae）、唇形科（Lamiaceae）、豆科（Fabaceae）、堇菜科（Violaceae）、菊科（Asteraceae）、毛茛科（Ranunculaceae）、傘形科（Apiaceae）、莎草科（Cyperaceae）植物為主。

山坡上部高于500 m的高海拔地區(qū)林型是以紅松和云冷杉為主的針闊混交林和以白樺、山楊、紫椴和蒙古櫟混交為主的落葉闊葉混交林，林下灌木主要有榛（Corylus heterophylla）、華北繡線菊（Spiraea fritschiana）、刺五加（Eleutherococcus senticosus）等，蓋度多在10%～20%之間，草本以朝鮮唐松草（Thalictrum ichangense var. coreanum）、廣布野豌豆（Vicia cracca）、興安獨(dú)活（Heracleum dissectum）等為主，蓋度多在30%[JP+3]～45%[JP+3]之間。山坡上部低于500 m的低海拔地區(qū)以純蒙古櫟林或以蒙古櫟、五角槭、花曲柳（Fraxinus chinensis subsp. rhynchophylla）、遼椴（Tilia mandshurica）等為主的混交林，灌木以大字杜鵑（Rhododendron schlippenbachii）、興安杜鵑（R. dauricum）、薄葉山梅花（Philadelphus tenuifolius）、東北溲疏（Deutzia parviflora var. amurensis）、胡枝子（Lespedeza bicolor）為主，蓋度多在20%～30%之間；林下草本層以山羅花（Melampyrum roseum）、歐鈴蘭（Convallaria majalis）、落新婦（Astilbe chinensis）、耳葉蟹甲草（Parasenecio auriculatus）、山尖子（P. hastatus）、山茄子（Brachybotrys paridiformis）、單穗升麻（Actaea simplex）以及大披針薹草（Carex lanceolata）、羊須草（C. callitrichos）等草本植物為主，草本植物蓋度多在20%～30%之間。

山坡下部及谷底在海拔高于500 m的區(qū)域其林型以針闊混交林為主，在低于500 m的區(qū)域以雜木林為主，主要樹種有黃花落葉松（Larix olgensis）、臭冷杉、遼東榿木（Alnus hirsuta）、水曲柳（Fraxinus mandshurica）、榆樹（Ulmus pumila）、青楷槭、紫花槭（Acer pseudosieboldianum）等，灌木以雞樹條（Viburnum opulus subsp. calvescens）、珍珠梅、遼東丁香（Syringa villosa subsp. wolfii）為主，蓋度多在30%～40%之間，林下草本主要有敗醬（Patrinia scabiosifolia）、水金鳳（Impatiens noli-tangere）、展枝沙參（Adenophora divaricata）、蹄葉橐吾（Ligularia fischeri）、驢蹄草（Caltha palustris）、毛茛（Ranunculus japonicus）、烏頭（Aconitum carmichaelii）、蚊子草（Filipendula palmata）、狹葉蕁麻（Urtica angustifolia）、荷青花（Hylomecon japonica）、拂子茅（Calamagrostis epigeios）、扁桿薹草（Carex planiculmis）、寬葉薹草（C. siderosticta）、蕨屬植物，草本植物蓋度多在30%～45%之間。

2.2 地形對(duì)密林下灌草層植物生物量的影響

2.2.1 密林下灌草層植物生物量本底情況 密林下1 685個(gè)樣方數(shù)據(jù)結(jié)果顯示，灌木層植物生物量為（12.26±0.30） g·m-2（x[TX-1mm]±sx[TX-0.5mm]），第一四分位數(shù) （Q1）為3.25 g·m-2，中位數(shù)（Q2）為9.11 g·m-2，第三四分位數(shù) （Q3）為17.51 g·m-2；草本層植物生物量為（17.84±0.39） g·m-2（x[TX-1mm]±sx[TX-0.5mm]），Q1為5.44 g·m-2， Q2為13.46 g·m-2，Q3為25.87 g·m-2；灌草層地上總生物量為（30.15±0.44） g·m-2（x[TX-1mm]±sx[TX-0.5mm]），Q1為16.50 g·m-2，Q2為26.45 g·m-2，Q3為40.40 g·m-2（圖1）。

2.2.2地形對(duì)密林下灌草層植物總生物量的影響 混合線性模型嵌套多因素方差分析結(jié)果（表3）顯示，坡位、坡向、坡度對(duì)密林下灌草層植物生物量均影響顯著（Plt;0.01），海拔無顯著影響（Pgt;0.05）。坡位與坡向交互作用顯著（表 3），其他交互作用均不顯著。Tukey多元比較結(jié)果顯示：整體而言，3個(gè)坡位等級(jí)間差異顯著，灌草層植物生物量由高到低依次為谷底、坡上及坡下（圖2：A）；陰坡灌草層植物生物量顯著低于陽坡及平地，后二者間無顯著差異（圖2：A）；3個(gè)坡度間差異顯著，灌草層植物生物量由高到低依次為陡坡、平坡及緩坡（圖2：A）； 坡位與坡向的交互作用（圖2：B）顯示，坡下平地、坡上平地、坡上陽坡與谷底所有坡位的灌草層植物生物量顯著最高，之間無顯著差異，坡下陰坡、坡下陽坡及坡上陰坡顯著最低，之間無顯著差異。

2.2.3 有序邏輯斯蒂回歸計(jì)算不同地形組合的林下灌草總生物量概率 有序邏輯斯蒂回歸結(jié)果顯示，模型擬合良好，通過平行線檢驗(yàn)（P=0.72）。海拔、坡位、坡向坡度3個(gè)分類變量所對(duì)應(yīng)的概率P值均小于顯著水平（Plt;0.05），進(jìn)入回歸方程（表4）。

由發(fā)生比歸納可知，低海拔密林下灌草層植物生物量是高海拔的0.76；坡下是谷底的0.47，坡上是谷底的0.69；以陽坡陡坡為準(zhǔn)， 陰坡緩坡的灌

草層植物生物量是陽坡陡坡的0.35，陽坡緩坡為陽坡陡坡的0.49，陰坡陡坡為陽坡陡坡的0.66，平坡為陽坡陡坡的0.63（表4）。

依據(jù)該有序邏輯斯蒂回歸的結(jié)果，得到連接函數(shù)Link1（公式8）和Link2（公式9）如下：

Link1=- 1.82 - （α海拔 +β坡位 +γ坡向坡度）（8）

Link2=- 0.03 - （α海拔 +β坡位 +γ坡向坡度）（9）

式中： α、β、γ分別為自變量海拔、坡位與坡向坡度的系數(shù)。

海拔、坡位及坡向坡度三類分類變量在公式8和公式9取值規(guī)則如下：（1）海拔是二分類變量，α取值-0.27，低海拔取值1，高海拔取值2；（2）坡位是三分類，按表5系數(shù)列按坡位（1=坡下，2=坡上，3=谷底）選擇相對(duì)應(yīng)β系數(shù)乘1進(jìn)入公式；（3）坡向坡度是五分類，按表5系數(shù)列按坡向坡度（1=陰坡緩坡；2=陽坡緩坡；3=陰坡陡坡；4=平坡；5=陽坡陡坡）選擇相對(duì)應(yīng)γ系數(shù)乘1進(jìn)入公式。每個(gè)因素水平中有1個(gè)分類是參數(shù)對(duì)照（其系數(shù)為0）（張文彤， 2002）。連接函數(shù)計(jì)算完畢，代入預(yù)測(cè)公式進(jìn)行下一步計(jì)算。

林下灌草層植物低、中、高生物量概率預(yù)測(cè)公式如下：

P低生物量=1/ ［1+exp （-Link1）］（10）

P中生物量=1/ ［1+exp （-Link2）］ -P低生物量（11）

P高生物量=1 -P低生物量 -P中生物量（12）

根據(jù)有序邏輯斯蒂回歸預(yù)測(cè)公式（公式10～公式12），計(jì)算各自變量組合下林下灌草層植物低、 中、 高生物量的概率（表6）。由表6可知， 低海拔的坡下區(qū)域低生物量概率普遍較大，大部分的組合中生物量概率較高，高生物量概率最高僅出現(xiàn)在高海拔坡上區(qū)域的陽坡陡坡。

3 討論

3.1 地形對(duì)密林下灌草層植物生物量的影響

研究區(qū)密林下灌草層植物總生物量約為0.30 t·hm-2，與Liu等（2019）在河南省雞公山的研究結(jié)果一致（林下總生物量變化范圍為0.08～0.69 t·hm-2， 平均為0.30 t·hm-2）。研究區(qū)密林下草本層植物生物量為0.18 t·hm-2，低于同緯度地帶40°～60° N森林平均地上部分草本層植物生物量的0.49 t·hm-2，以及北緯30°～40°地帶的0.42 t·hm-2（de Angelis et al.， 1981; Gilliam， 2003）。這可能是因?yàn)椴煌瑓^(qū)域群落物種組成不同導(dǎo)致，另外，林下植被生物量在大尺度上受到緯度、年均溫和年降水的顯著影響（Jin et al.， 2022）。

本研究區(qū)尺度較小，灌草層植物生物量空間分異顯著，海拔、坡位、坡向與坡度均顯示出較強(qiáng)的影響（張婷等， 2023）。高海拔處灌草層植物生物量高于低海拔處（約10%），通過野外走訪與調(diào)查分析認(rèn)為可能是由于低海拔區(qū)域人為干擾（如放牛、人參種植、林蛙養(yǎng)殖、山貨采集等）相對(duì)較高，但不同海拔數(shù)據(jù)存在較高的組內(nèi)變異（變異系數(shù)為60%），其統(tǒng)計(jì)結(jié)果并不顯著。細(xì)尺度上變化較大的坡位、坡度和坡向?qū)嗖輰又参锷锪烤嬖陲@著影響。坡位和坡度影響局域土壤的厚度，隨之土壤物化性質(zhì)和水肥條件也受到一系列影響（Power et al.， 1981; Woods amp; Schuman， 1988; Pachepsky et al.， 2001; 柴勇等， 2017）。一般而言，灌草生產(chǎn)力在谷底（多為塌積或沖積土壤）較高，而陡坡高于緩坡雖不符合預(yù)期但值得討論。首先，灌草層植物在坡度較大的林地比相同條件的緩坡、平地得到更多來自坡向方向的側(cè)光光照（陡坡效應(yīng)），會(huì)導(dǎo)致陡坡灌草生產(chǎn)力較高，同時(shí)陡坡土層較薄導(dǎo)致喬木密度較低、樹木個(gè)體矮小等客觀現(xiàn)象 （野外調(diào)查數(shù)據(jù)）， 會(huì)加劇陡坡效應(yīng)， Mou和Warrillow（2000）在研究冰雨對(duì)坡地森林影響時(shí)對(duì)此有所共識(shí)。其次， 陡坡條件下地勢(shì)積累的枯落物相對(duì)較少，這也是其灌草層生物量較高的一個(gè)可能原因（徐來仙等， 2022）。平地灌草層生物量高于坡地應(yīng)歸因于較好的土壤因素（Marques et al.， 2004），而陰、陽坡間的差異則可能既有光、溫差異的影響，也有水分條件的限制（如陽坡較干旱）。在溫帶濕潤(rùn)地區(qū)，陽坡光照、溫度條件均好于陰坡，水分條件則相反，水熱條件間的相互作用影響樹木生長(zhǎng)，也影響林下植物的生長(zhǎng)（李俊清， 2010）。陽坡灌草層植物生物量顯著高于陰坡，也可說明此處光照對(duì)于灌草層的重要性較高，尤其是在研究區(qū)里存在著大面積的中幼齡次生闊葉林，由于林分密度大、郁閉度高，平地、緩坡林下遮陰更甚。在長(zhǎng)白山地區(qū)，紅松針闊混交林生態(tài)系統(tǒng)喬、灌、草全年平均光合速率分別占到整個(gè)森林生態(tài)系統(tǒng)光合作用的89.9%、3.4%、6.7%（王淼等， 2006），可見灌草層的限制因素更多是光照，這可以間接地用來解釋本研究結(jié)果。

微地形對(duì)灌草生物量的影響普遍存在（Gilliam， 2003），坡位和坡向?qū)嗖萆锪匡@著的交互作用表明地形作用的復(fù)雜性。谷底平地與不同坡向灌草層植物生物量均無顯著差異，是由于在該地形下陰坡、陽坡、平坡基本處在密林冠之下，其林下光照條件并無明顯差異。此外，該區(qū)域土壤一般深厚，相對(duì)凹陷的局域地貌使這里水分條件良好，濕潤(rùn)環(huán)境指示物種——水金鳳是這里的常見物種（Hiratsuka amp; Inoue， 2013），并不局限于任一坡向。坡上林下灌草層植物生物量變化與坡向變化一致，光照與水分對(duì)林下灌草層植物生物量的相對(duì)重要性表現(xiàn)一致，但顯然前者作用應(yīng)該更大一些。坡下平地林下灌草層植物生物量高于坡下陰陽坡，也高于其他所有坡位坡向組合，對(duì)此本研究尚無法給出令人信服的解釋，這一問題亟待進(jìn)一步研究。

3.2 有序邏輯斯蒂回歸分析討論

鑒于研究區(qū)密林下灌草層生長(zhǎng)、分布的復(fù)雜性，利用地形因素作為綜合影響因子計(jì)算林下灌草層植物生物量的概率是一種比較可行的方法。通過對(duì)1 685個(gè)樣方調(diào)查數(shù)據(jù)的有序邏輯斯蒂回歸結(jié)果的分析，可以對(duì)研究區(qū)灌草層現(xiàn)狀歸納出下列趨勢(shì)和特點(diǎn)：（1）谷底各地形組合高灌草層植物生物量的概率高；（2）坡下各地形組合低灌草層植物生物量概率高；（3）陽坡陡坡高灌草層植物生物量概率高；（4）以上趨勢(shì)在高低海拔間差異較小。

部分結(jié)果與預(yù)期有些差異，如環(huán)境條件較差的陡坡地段有相對(duì)高的高生物量概率（趙磊磊等， 2012）。實(shí)際上，高海拔坡上的陰坡陡坡和陽坡陡坡的高灌草層植物生物量都有最高的概率，而一般環(huán)境較好的低海拔坡下，卻是低灌草層植物生物量概率高的地方。如果不考慮干擾因素，較好光照的陡坡效應(yīng)是一個(gè)合理解釋，并且陰坡陡坡和陽坡陡坡在高低海拔、不同坡位上較為一致的差異也對(duì)光照說形成了一定的支持。然而，整個(gè)陡坡的放牧強(qiáng)度相對(duì)較低，人為干擾相對(duì)較少也是事實(shí)，尤其在高海拔區(qū)域（Feng et al.， 2021）。因此，如何具體量化這些影響還有待進(jìn)一步的控制實(shí)驗(yàn)研究。

一般而言，坡下土壤深厚、水肥條件較好（Chamran et al.， 2002），植物的生產(chǎn)力和生物量均較高（Liu et al.， 2018），然而事實(shí)與模型預(yù)測(cè)結(jié)果均與此相悖。坡下好的水肥溫度條件對(duì)喬木生長(zhǎng)的促進(jìn)作用應(yīng)該更大，但實(shí)地調(diào)查顯示喬木平均胸徑相對(duì)較低（王樂等， 2019b）。此地形區(qū)域距離村鎮(zhèn)道路近，并且立地條件好，家畜（牛）活動(dòng)強(qiáng)度較高，人類活動(dòng)密度也遠(yuǎn)較高海拔、深山區(qū)的高（Feng et al.， 2021）。因此，在未排除干擾影響的調(diào)查數(shù)據(jù)基礎(chǔ)上計(jì)算的概率，對(duì)灌草生產(chǎn)力及質(zhì)量高低進(jìn)行判斷存在不足。未來將設(shè)置更加系統(tǒng)的圍欄調(diào)查實(shí)驗(yàn)，進(jìn)而對(duì)灌草層生產(chǎn)力進(jìn)行更為準(zhǔn)確地計(jì)算。

4 結(jié)論

地形因素（坡位、坡向及坡度）對(duì)密林下灌草層植物生物量有顯著影響，灌草層植物生物量在不同地形組合下有顯著差異。谷底的灌草層植物生物量高于坡上，坡上高于坡下（Plt;0.01）；陰坡灌草層植物生物量顯著低于陽坡及平地（Plt;0.01），后二者間無顯著差異；不同坡度，陡坡灌草層植物生物量高于平坡，平坡高于緩坡（Plt;0.01）。坡位與坡向的交互作用顯著，坡下平地、坡上平地、坡上陽坡與谷底的所有坡位灌草層植物生物量顯著最高，坡下陰坡、坡下陽坡及坡上陰坡之間無顯著差異均顯著最低。坡位、坡向及坡度對(duì)林下灌草層植物生物量有顯著影響，3個(gè)坡位等級(jí)間谷底最高而下坡位最低，3個(gè)坡度等級(jí)間陡坡最高而緩坡最低，不同坡向比較陰坡最低。

在研究區(qū)現(xiàn)行狀態(tài)，灌草層植物生物量高低在不同海拔、坡位及坡向坡度組合有不同的概率。在不排除人為干擾、森林放牧的情況下，谷底、陡坡地帶灌草層植物生物量高產(chǎn)概率最高。

致謝 野外調(diào)查工作得到了國(guó)家林業(yè)和草原局（國(guó)家公園管理局）、吉林省林業(yè)和草原局、琿春國(guó)家級(jí)自然保護(hù)區(qū)管理局、琿春市林業(yè)局、汪清國(guó)家級(jí)自然保護(hù)區(qū)管理局和汪清縣林業(yè)局、中國(guó)人民解放軍某邊防部隊(duì)大力支持。同時(shí)感謝東北虎豹國(guó)家公園保護(hù)生態(tài)學(xué)國(guó)家林草局重點(diǎn)實(shí)驗(yàn)室虎豹團(tuán)隊(duì)所有成員，以及張全國(guó)老師在使用R語言進(jìn)行嵌套方差分析上的指導(dǎo)和幫助。

參考文獻(xiàn)：

BARANACˇU3EKOV M， KROJEROV-PROKEOV J， USTR P， et al.， 2010. Annual changes in roe deer （Capreolus capreolus L.） diet in the Bohemian Forest， Czech Republic/Germany" ［J］. Eur J Wildl Res， 56： 327-333.

BASKIN L， DANELL K， 2003. Ecology of ungulates： a handbook of species in Eastern Europe and Northern and Central Asia" ［M］. Heidelberg： Springer-Verlag： 15-122.

CHAI Y， MENG GT， HE LP， et al.， 2017. Distribution pattens of tree species and its influencing factors of topography in a mid-montane humid ever-green broad-leaved forest in Gaoligong Mountains" ［J］. Guihaia， 37（12）： 1508-1520. ［柴勇， 孟廣濤， 和麗萍， 等， 2017. 高黎貢山中山濕性常綠闊葉林樹種的分布格局及其地形影響因子" ［J］. 廣西植物， 37（12）： 1508-1520.］

CHAMRAN F， GESSLER PE， CHADWICK OA， 2002. Spatially explicit treatment of soil-water dynamics along a semiarid catena" ［J］. Soil Sci Soc Am J， 66（5）： 1571-1583.

DE ANGELIS DL， GARDNER RH， SHUGSRT HH， 1981. Productivity of forest ecosystems studied during the IBP： the woodlands data set"" ［M］//REICHLE DE. Dynamics properties of forest ecosystems. Cambridge： Cambridge University Press： 567-672.

DE CASTILHO CV， MAGNUSSON WE， DE ARAU＇JO RNO， et al.， 2006. Variation in aboveground tree live biomass in a central Amazonian Forest： effects of soil and topography" ［J］. For Ecol Manage， 234（1/2/3）： 85-96.

FENG LM， SHEVTSOVA E， VITKALOVA A， et al.， 2017. Collaboration brings hope for the last Amur leopards" ［J］. Cat News， 65： 20.

FENG RN， L XY， XIAO WH， et al.， 2021. Effects of free-ranging livestock on sympatric herbivores at fine spatiotemporal scales" ［J］. Landscape Ecol， 36： 1441-1457.

GILLIAM FS， 2003. The herbaceous layer in forests of eastern North America" ［M］. New York： Oxford University Press： 24-26.

GONG YB， 2016. Understory species diversity and biomass of natural broad-leaved Korean pine forest" ［D］. Beijing： Beijing Forestry University." ［龔艷賓， 2016. 天然闊葉紅松林下植被物種多樣性及生物量研究" ［D］. 北京： 北京林業(yè)大學(xué).］

GORDON IJ， PRINS HHT， 2008. The ecology of browsing and grazing" ［M］. Heidelberg： Springer-Verlag： 126-177.

GRACIA M， MONTAN F， PIQU J， et al.， 2007. Overstory structure and topographic gradients determining diversity and abundance of understory shrub species in temperate forests in central Pyrenees （NE Spain）" ［J］. For Ecol Manage， 242（2/3）： 391-397.

HEBBLEWHITE M， MIQUELLE DG， ROBINSON H， et al.， 2014. Including biotic interactions with ungulate prey and humans improves habitat conservation modeling for endangered Amur tigers in the Russian Far East" ［J］. Biol Conserv， 178： 50-64.

HIRATSUKA A， INOUE O， 2013. Dispersability of dimorphic seeds in Impatiens noli-tangere and I. textori （Balsaminaceae）" ［J］. Ecol Rev， 21（3）： 157-161.

JIANG HQ， 2011. Spatial-temporal analysis and prediction of forest above-ground biomass in Jingouling forest farm" ［D］. Beijing： Chinese Academy of Forestry." ［姜慧泉， 2011. 金溝嶺林場(chǎng)森林地上生物量時(shí)空動(dòng)態(tài)分析與預(yù)測(cè)" ［D］. 北京： 中國(guó)林業(yè)科學(xué)研究院.］

JIN Y， LIU C， QIAN SS， et al.， 2022. Large-scale patterns of understory biomass and its allocation across Chinas forests ［J］. Sci Total Environ， 804： 150169.

JOLLY WM， NEMANI R， RUNNING SW， 2004. Enhancement of understory productivity by asynchronous phenology with overstory competitors in a temperate deciduous forest" ［J］. Tree Physiol， 24（9）： 1069-1071.

KERLEY LL， MUKHACHEVA AS， MATYUKHINA DS， et al.， 2015. A comparison of food habits and prey preference of Amur tiger （Panthera tigris altaica） at three sites in the Russian Far East" ［J］. Integr Zool， 10（4）： 354-364.

LI J， GONG CW， 2021. Effects of terrain factors on vegetation cover change on national park of Qilian Mountains" ［J］. Bull Soil Water Conserv， 41（3）： 228-237." ［李娟， 龔純偉， 2021. 祁連山國(guó)家公園植被覆蓋變化地形分異效應(yīng)" ［J］. 水土保持通報(bào)， 41（3）： 228-237.］

LI JQ， 2010. Forest ecology" ［M］. 2nd ed. Beijing： Higher Education Press： 16-96." ［李俊清， 2010. 森林生態(tài)學(xué)" ［M］. 2版. 北京： 高等教育出版社： 16-96.］

LI Q， RONG L， WANG M， 2019. Effects of topography on diversity and distribution pattern of plant species in karst mountains area" ［J］. Bull Soil Water Conserv， 39（6）： 27-34." ［李芹， 容麗， 王敏， 2019. 地形對(duì)喀斯特山地植物物種多樣性及分布格局的影響" ［J］. 水土保持通報(bào)， 39（6）： 27-34.］

LI R， 2010. Biomass estimates and its spatial pattern of a primary temperate forest in Xiaoxingan Mountains， Northeast China" ［D］. Harbin： Northeast Forestry University." ［黎如， 2010. 小興安嶺原始闊葉紅松林生物量及其空間分布格局" ［D］. 哈爾濱： 東北林業(yè)大學(xué).］

LIU HT， YAO L， LIN CW， et al.， 2018. 18-year grass hedge effect on soil water loss and soil productivity on sloping cropland" ［J］. Soil Tillage Res， 177： 12-18.

LIU Y， SHANG Q， WANG L，et al.， 2019. Effects of understory shrub biomass on variation of soil respiration in a temperate-subtropical transitional oak forest" ［J］. Forests， 10（2）： 88.

MARQUES JJ， SCHULZE DG， CURI N， et al.， 2004. Major element geochemistry and geomorphic relationships in Brazilian Cerrado soils" ［J］. Geoderma， 119（3/4）： 179-195.

MASSEI G， GENOV PV， STAINES BW， 1996. Diet， food availability and reproduction of wild boar in a Mediterranean coastal area" ［J］. Acta Theriol， 41（3）： 307-320.

MOU P， WARRILLOW MP， 2000. Ice storm damage to a mixed hardwood forest and its impacts on forest regeneration in the ridge and valley region of southwestern Virginia" ［J］. J Torrey Bot Soc， 127（1）： 66-82.

NIE X， GUO W， HUANG B， et al.， 2019. Effects of soil properties， topography and landform on the understory biomass of a pine forest in a subtropical hilly region" ［J］. Catena， 176： 104-111.

PACHEPSKY YA， TIMLIN DJ， RAWLS WJ， 2001. Soil water retention as related to topographic variables" ［J］. Soil Sci Soc Am J， 65（6）： 1787-1795.

POWER JF， SANDOVAL FM， RIES RE， et al.， 1981. Effects of topsoil and subsoil thickness on soil water content and crop production on a disturbed soil" ［J］. Soil Sci Soc Am J， 45（1）： 124-129.

QIN ST， LONG CL， WU BL， 2019. Niche characteristics of dominant tree populations in karst forest at different topography sites in Maolan National Natural Reserve of Guizhou Province" ［J］. Guihaia， 39（5）： 681-689." ［秦隨濤， 龍翠玲， 吳邦利， 2019. 茂蘭喀斯特森林不同地形部位優(yōu)勢(shì)喬木種群的生態(tài)位研究" ［J］. 廣西植物， 39（5）： 681-689.］

R Core Team， 2021. R： A language and environment for statistical computing. R Foundation for Statistical Computing， Vienna， Austria" ［EB/OL］. https：//www.R-project.org/.

SHEN WH， OU ZY， PANG SL， et al.， 2017. Relationship between the spatial distributions of dominant arbor species and environmental factors of Excentrodendron hsienmu community in Southwest Guangxi， China" ［J］. Guihaia， 37（6）： 694-701." ［申文輝， 歐芷陽， 龐世龍， 等， 2017. 桂西南蜆木群落優(yōu)勢(shì)樹種分布與環(huán)境因子的關(guān)系" ［J］. 廣西植物， 37（6）： 694-701.］

SUN LY， LI QY， PEI L， et al.， 2020. Effects of topographic factors on soil physical and chemical properties and plant species" ［J］. J Irrig Drain， 39（7）： 120-127." ［孫莉英， 栗清亞， 裴亮， 等， 2020. 地形因子對(duì)土壤理化性質(zhì)和植物種類的影響" ［J］. 灌溉排水學(xué)報(bào)， 39（7）： 120-127.］

WANG L， FENG JW， AMARSAIKHAN T， et al.， 2019a. Forest cattle grazing affects understory food resource of ungulates in the eastern part of the Northeast Tiger and Leopard National Park" ［J］. Acta Theriol Sin， 39（4）： 386-396." ［王樂， 馮佳偉， AMARSAIKHAN T， 等， 2019a. 森林放牧對(duì)東北虎豹國(guó)家公園東部有蹄類動(dòng)物灌草層食物資源的影響" ［J］. 獸類學(xué)報(bào)， 39（4）： 386-396.］

WANG L， FENG JW， MOU P， et al.， 2023. Relative abundance of Roe deer （Capreolus pygargus） related to overstory structure and understory food resources in Northeast China" ［J］. Glob Ecol Conserv， 46： e02542.

WANG L， YANG LM， SAI JM， et al.， 2019b. The quantity and quality of understory forages of the ungulates’ habitat in the eastern part of Northeast Tiger and Leopard National Park" ［J］. Acta Theriol Sin， 39（4）： 373-385." ［王樂， 楊麗萌， 賽九瑪， 等， 2019b. 東北虎豹國(guó)家公園東部有蹄類棲息地的灌草層食物資源的量與質(zhì) "［J］. 獸類學(xué)報(bào)， 39（4）： 373-385.］

WANG M， GUAN DX， WANG YS， et al.， 2006. Estimation of ecosystem productivity of red pine needle broad-leaved mixed forest in Changbai Mountain" ［J］. Sci Chin Ser D Earth Sci， 36（Suppl. Ⅰ）： 70-82." ［王淼， 關(guān)德新， 王躍思， 等， 2006. 長(zhǎng)白山紅松針闊葉混交林生態(tài)系統(tǒng)生產(chǎn)力的估算 ［J］. 中國(guó)科學(xué)D輯 地球科學(xué)， 36（增刊 Ⅰ）： 70-82.］

WANG TM， FENG LM， MOU P， et al.， 2016. Amur tigers and leopards returning to China： direct evidence and a landscape conservation plan" ［J］. Landscape Ecol， 31： 491-503.

WANG TM， FENG LM， YANG HT， et al.， 2020. An introduction to Long-term Tiger-Leopard Observation Network based on camera traps in Northeast China" ［J］. Biodivers Sci， 28（9）： 1059-1066." ［王天明， 馮利民， 楊海濤， 等， 2020. 東北虎豹生物多樣性紅外相機(jī)監(jiān)測(cè)平臺(tái)概述" ［J］. 生物多樣性， 28（9）： 1059-1066.］

WANG WF， SONG LN， SUI X， 2010. Estimation of forest biomass and its temporal and spatial distribution patterns in Maoershan Forest Farm， Heilongjiang Province" ［J］. J NE For Univ， 38（1）： 47-49." ［王維芳， 宋麗楠， 隋欣， 2010. 帽兒山林場(chǎng)森林生物量估測(cè)及時(shí)空動(dòng)態(tài)格局分析" ［J］. 東北林業(yè)大學(xué)學(xué)報(bào)， 38（1）： 47-49.］

WANG XL， 2013. Biomass allocation and spatial pattern of major forest ecosystems in Great Xingan Mountains， Heilongjiang Province" ［D］. Beijing： University of Chinese Academy of Sciences." ［王曉莉， 2013. 黑龍江省大興安嶺主要森林生態(tài)系統(tǒng)生物量分配及空間格局研究" ［D］. 北京： 中國(guó)科學(xué)院大學(xué).］

WANG XS， 2019. Comparison of nutritional strategies between Siberian roe deer and sika deer during winter in the north of Changbai Mountains， China" ［D］. Harbin： Northeast Forestry University." ［王祥生， 2019. 長(zhǎng)白山北部山地東北梅花鹿與同域分布西伯利亞狍冬季營(yíng)養(yǎng)對(duì)策比較研究" ［D］. 哈爾濱： 東北林業(yè)大學(xué).］

WHITTAKER RH， LIKENS GE， 1973. Primary production： the biosphere and man" ［J］. Hum Ecol， 1： 357-369.

WOODS LE， SCHUMAN GE， 1988. Cultivation and slope position effects on soil organic matter" ［J］. Soil Sci Soc Am J， 52（5）： 1371-1376.

XIAO WH， FENG LM， ZHAO XD， et al.， 2014. Distribution and abundance of Amur tiger， Amur leopard and their ungulate prey in Hunchun National Nature Reserve， Jilin" ［J］. Biodivers Sci， 22（6）： 717-724." ［肖文宏， 馮利民， 趙小丹， 等， 2014. 吉林琿春自然保護(hù)區(qū)東北虎和東北豹及其有蹄類獵物的多度與分布" ［J］. 生物多樣性， 22（6）： 717-724.］

XU LX， YAO L， AI XR， et al.， 2022. Effects of vegetation and litter on natural seedling germination of Metasequoia glyptostroboides" ［J］. J For Environ， 42（5）： 456-464." ［徐來仙， 姚蘭， 艾訓(xùn)儒， 等， 2022. 植被和凋落物對(duì)水杉天然幼苗萌發(fā)的影響" ［J］. 森林與環(huán)境學(xué)報(bào)， 42（5）： 456-464.］

YANG HT， DOU HL， BANIYA RK， et al.， 2018. Seasonal food habits and prey selection of Amur tigers and Amur leopards in Northeast China" ［J］. Sci Rep-UK， 8： 6930.

ZHANG JL， LIU FZ， CUI GF， 2016. Spatio-temporal variation of vegetation and analysis of its driving factors in Changbai Mountain National Nature Reserve" ［J］. Acta Ecol Sin， 36（12）： 3525-3536." ［張建亮， 劉方正， 崔國(guó)發(fā)， 2016. 長(zhǎng)白山國(guó)家級(jí)自然保護(hù)區(qū)植被時(shí)空變化及其驅(qū)動(dòng)因子" ［J］. 生態(tài)學(xué)報(bào)， 36（12）： 3525-3536.］

ZHANG L， CHENG YC， FENG LM， et al.， 2014. Abundance and occupancy of leopard and their prey in Wangqing Leopard Reserve， China" ［C］. 99th ESA Annual Meeting （August 10-15， 2014）. ［2018-08-14］.

ZHANG T， ZHANG JL， YANG T， et al.， 2023. Influence of topographic factors on plant diversity and spatial distribution of karst around FAST" ［J］. Guihaia， 42（3）： 473-483." ［張婷， 張建利， 楊濤， 等， 2023. 地形因子對(duì)FAST周邊喀斯特植物多樣性及空間分布的影響" ［J］. 廣西植物， 42（3）： 473-483.］

ZHANG WT， 2004. SPSS statistical analysis of advanced tutorial ［M］. Beijing： Higher Education Press： 189-193. ［張文彤， 2004. SPSS統(tǒng)計(jì)分析高級(jí)教程 ［M］. 北京： 高等教育出版社： 189-193.］

ZHAO LL， ZHU QK， NIE LS， et al.， 2012. Soil moisture variation patterns of steep slope in the loess region in northern Shaanxi Province" ［J］. Ecol Environ Sci， 21（2）： 253-259." ［趙磊磊， 朱清科， 聶立水， 等， 2012. 陜北黃土區(qū)陡坡土壤水分變異規(guī)律研究" ［J］. 生態(tài)環(huán)境學(xué)報(bào)， 21（2）： 253-259.］

ZHOU YL， LI JW， 1964. Characteristics and distribution of main vegetation types in eastern mountains of Northeast China" ［J］. Chin J Plant Ecol， 2（2）： 190-206." ［周以良， 李景文， 1964. 中國(guó)東北東部山地主要植被類型的特征及其分布規(guī)律" ［J］. 植物生態(tài)學(xué)報(bào)， 2（2）： 190-206.］

（責(zé)任編輯 鄧斯麗）

基金項(xiàng)目：" 中央級(jí)公益性科研院所基本科研業(yè)務(wù)費(fèi)專項(xiàng)（CAFYBB2023MA017）； 國(guó)家自然科學(xué)基金（31971539）； 科技部國(guó)家重點(diǎn)研發(fā)計(jì)劃（2019FY101702）； 中國(guó)林業(yè)科學(xué)研究院生態(tài)保護(hù)與修復(fù)研究所橫向項(xiàng)目（202204021-4102）。

第一作者： 王樂（1989—），博士，助理研究員，主要從事植被生態(tài)學(xué)與生態(tài)系統(tǒng)生態(tài)學(xué)研究，（E-mail）wangle@caf.ac.cn。

*通信作者：" 牟溥，博士，教授，主要從事陸地生態(tài)學(xué)與景觀生態(tài)學(xué)研究，（E-mail）ppmou@bnu.edu.cn。