王 磊，崔 明，劉玉國，周夢玲，武建宏，周桃龍

（1.中國林業(yè)科學研究院 荒漠化研究所，北京 100091；2.北京林業(yè)大學 水土保持學院，北京 100083；3.河南省淅川縣林業(yè)局，河南 淅川 474450）

林分生物量是森林生態(tài)系統(tǒng)對有機質(zhì)所積累的質(zhì)量，是森林生態(tài)系統(tǒng)最基本的特征數(shù)據(jù)，對研究森林生態(tài)系統(tǒng)物質(zhì)循環(huán)、能量流動具有重要的作用[1?3]。林分生物量的動態(tài)變化不僅能夠反映物質(zhì)與能量在生態(tài)系統(tǒng)中變化的復(fù)雜關(guān)系，同時在一定程度上能夠反映生態(tài)系統(tǒng)的健康水平[4?5]。生物量的積累表現(xiàn)為林分生長。立地條件與林分生長密切相關(guān)，好的立地條件為林分生長提供充足的養(yǎng)分，加快林分生長[6]。林分結(jié)構(gòu)和樹種多樣性也會影響林分有機質(zhì)的積累速度。巖溶地區(qū)高溫多雨且分布不均，基巖裸露率高、土壤瘠薄、季節(jié)性巖溶干旱明顯[7?8]。巖溶區(qū)的適生植被不僅需具有石生性、喜鈣性、耐干旱和耐貧瘠的特點，而且需具有發(fā)達的根系。植被的恢復(fù)和重建是巖溶退化區(qū)的重點和難點[9]，生物量作為生態(tài)系統(tǒng)生產(chǎn)力的基礎(chǔ)和功能的表現(xiàn)形式，是退化區(qū)植被恢復(fù)成功與否的評價指標[10]。目前對巖溶生態(tài)系統(tǒng)植被生物量的研究集中在年齡序列、演替階段等單一因素對生物量的影響[11?13]。然而，生態(tài)系統(tǒng)中的影響是一個復(fù)雜的過程，由很多相互的作用組成[14?15]，因此綜合探討林分生物量與多種影響因素之間的關(guān)系具有必要性和現(xiàn)實性。傳統(tǒng)的路徑分析及主成分分析不能同時處理多個因變量與自變量的聯(lián)系，不能提供完整的信息[16]。結(jié)構(gòu)方程模型是一種多變量的統(tǒng)計方法，能夠同時研究自然系統(tǒng)中多個因子間的相互關(guān)系[6,17]。相較于傳統(tǒng)的相關(guān)性分析，結(jié)構(gòu)方程模型是基于研究者的先驗知識預(yù)先設(shè)定系統(tǒng)內(nèi)因子間的依賴關(guān)系，能夠判別各因子間的關(guān)系強度，同時還能作出整體擬合和判斷，反映自變量與因變量之間的直接和間接關(guān)系，幫助研究者更全面掌握自然系統(tǒng)中的真實機制[6,18]。近年來結(jié)構(gòu)方程模型在生態(tài)學和林學上得到廣泛應(yīng)用[19?20]。河南省淅川縣是南北過渡帶和巖溶石漠化集連片區(qū)的最北端，巖溶面積達348.82 km2，又是南水北調(diào)中線工程核心水源地、主要淹沒區(qū)及渠首所在地，地理位置特殊[21]。調(diào)查發(fā)現(xiàn)：巖溶區(qū)林分質(zhì)量明顯低于非巖溶區(qū)，生態(tài)資源環(huán)境問題嚴重，影響著南水北調(diào)水中線安全。提高該區(qū)域生態(tài)環(huán)境，保障南水北調(diào)中線安全勢在必行，林分質(zhì)量提升已經(jīng)迫在眉睫。本研究以淅川縣巖溶區(qū)天然次生林為研究對象，利用結(jié)構(gòu)方程模型，探討巖溶區(qū)建群種生物量、林分總生物量與基巖裸露度、土壤厚度以及林分密度和物種多樣性之間的耦合關(guān)系，明確各關(guān)系之間影響強度，以期為河南省淅川縣巖溶區(qū)植被的經(jīng)營和管理提供理論依據(jù)，以期改善淅川縣巖溶區(qū)生態(tài)環(huán)境。

1 研究地區(qū)概況與研究方法

1.1 研究地區(qū)概況

淅川縣位于河南省西南部，豫、鄂、陜?nèi)〗Y(jié)合處，32°55′～33°23′N, 110°58′～111°53′E，地處伏牛山余脈，擁有亞洲最大的淡水湖，屬于南水北調(diào)中線工程水源地所在地。淅川縣氣候適宜，屬于北亞熱帶向暖溫帶過渡的季風性氣候區(qū)，四季分明。年平均降水量為804.3 mm，6?9月降水占全年降水量的58.0%以上，最高達79.8%。年平均氣溫為15.8 ℃，無霜期約228 d。巖溶地貌在淅川分布廣泛，主要集中在丹江口庫區(qū)周邊和荊關(guān)、西簧、金河、大石橋等12個鄉(xiāng)鎮(zhèn)，巖溶區(qū)現(xiàn)有植被以次生林為主。其中主要的喬木樹種為栓皮櫟Quercusvariabilis、黃連木Pistaciachinensis、化香Platycaryastrobilacea、側(cè)柏Platycladusorientalis、女貞Ligustrumlucidum等；主要的灌木樹種為小果薔薇Rosacymosa、菝葜Smilaxchina、黃荊Vitexnegundo、莢蒾Viburnumdilatatum、連翹Hypericumascyron等；主要的草本為碎米莎草Cyperusiria、綬草Spiranthessinensis、天門冬Asparaguscochinchinensis、藎草Arthraxon hispidus等。

1.2 調(diào)查方法與數(shù)據(jù)處理

2018年7?9月分別在淅川縣金河鎮(zhèn)、西簧鄉(xiāng)、荊關(guān)鎮(zhèn)巖溶區(qū)選擇6個樣點18個20 m×30 m樣地，在每個樣地中劃分6個10 m×10 m小樣方，共108個樣方進行數(shù)據(jù)采集。在每塊樣地進行每木檢尺，確定喬木層樹種組成，測定樹高、胸徑、枝下高和冠幅。同時測量每塊樣地的坡度、坡向、海拔、基巖裸露度、土壤厚度，測定土壤含水率和土壤容重等。

采用中國森林生態(tài)系統(tǒng)碳儲量—生物量方程[22]計算林分生物量和建群種生物量。

采用Shannon-Wiener指數(shù)計算各樣地樹種多樣性[23]：

式(1)中：H為樹種多樣性；B為林分胸徑面積總和；Bi為樹種i的胸徑斷面積總和；N為樹種數(shù)。結(jié)果見表1。

1.3 模型構(gòu)建原理

結(jié)構(gòu)方程模型是一種建立、估計和檢驗因果關(guān)系的模型，結(jié)合了因素分析和路徑分析2種統(tǒng)計方法，并且同時檢驗?zāi)Ｐ椭械娘@性變量、潛在變量、干擾或誤差變量間的關(guān)系。觀測變量是直接測得的變量、潛在變量是不可直接測得的變量，但是可由觀測變量表示。如基巖裸露度、土壤厚度、海拔等為可以直接測量的指標，這些指標能夠反映不能夠直接測量的指標立地條件(潛在變量)，本研究中采用能夠直接測量的指標構(gòu)建結(jié)構(gòu)方程模型。

結(jié)構(gòu)方程模型由測量模型和結(jié)構(gòu)模型兩部分組成。測量模型是潛在變量與顯性變量之間關(guān)系的模型。

式(2)～(3)中：X為外生顯性變量向量；Y為內(nèi)生顯性變量向量；Λx和Λy為指標變量(X,Y)的因素負荷量；δ、ε為外生顯性變量和內(nèi)生顯性變量的測量誤差；ξ為外生潛在變量；η為內(nèi)生潛在變量[24]。

立地條件決定林分的生長狀況，直接影響林分的結(jié)構(gòu)和樹木的多樣性以及林分的生物量，立地條件包含地形因子、土壤因子以及氣候因子[18]，森林的功能與林分結(jié)構(gòu)密切相關(guān)，林分結(jié)構(gòu)表征了林木的分布與生長狀況。林分多樣性影響樹種間的差異，并且影響林分的生長[25]。由于巖溶區(qū)基巖裸露度高，土壤瘠薄且不連續(xù)，對巖溶區(qū)植被有著最直接的影響。故本研究中將基巖裸露度(x1)、土壤厚度(x2)、樹種多樣性(x3)、林分密度(x4)等觀測變量作為探討與林分生物量關(guān)系的顯變量，對應(yīng)誤差為δ1、δ2、δ3、δ4。建群種生物量(η1)和林分總生物量(η2)作為與內(nèi)生潛在變量(生物量η)有關(guān)的觀測變量，ε1、ε2為對應(yīng)的誤差，構(gòu)建各顯變量與內(nèi)生潛變量(生物量η)結(jié)構(gòu)方程模型。運用SPSS 18.0軟件對數(shù)據(jù)進行處理，運用Amos 23.0軟件構(gòu)建結(jié)構(gòu)方程模型進行林分生長與影響因子間關(guān)系耦合。

表 1 樣地概況Table 1 Sample information

2 模型構(gòu)建與修正

本研究建立基巖裸露度(x1)、土壤厚度(x2)、樹種多樣性(x3)、林分密度(x4)等觀測變量與建群種生物量和林分生物量耦合關(guān)系的結(jié)構(gòu)方程模型。根據(jù)擬合結(jié)果，逐步的調(diào)整路徑關(guān)系，刪除關(guān)系不顯著(P＞0.05)的路徑，得到初始結(jié)構(gòu)方程模型1(圖1A)。初始模型χ2=13.236，df=7，χ2/df=1.891＜3，顯著性概率值P=0.067＞0.05，接受虛無假設(shè)。但模型適配度指數(shù)IGF=0.829＜0.90，模型規(guī)準適配指數(shù)(INF)、相對適配指數(shù)(IRF)、非規(guī)準適配指數(shù)(ITL)均較小(表2)，模型擬合結(jié)果需進一步進行模型修正。修正后得到最優(yōu)模型2(圖1B)。

最優(yōu)模型χ2=5.051，df=6，χ2/df=0.842＜3，顯著性概率值P=0.537＞0.05，接受虛無假設(shè)，觀測數(shù)據(jù)所導(dǎo)出的方差協(xié)方差矩陣與假設(shè)模型導(dǎo)出的方差協(xié)方差矩陣適配的假設(shè)獲得支持，假設(shè)模型與觀測數(shù)據(jù)適配。模型中各項檢驗指標均達到標準，如表2。同時漸進殘差均方和平方根RMESA=0＜0.08、以及赤池信息量AIC(預(yù)設(shè)模型)=35.051＜42.000(飽和模型)＜123.965(獨立模型)，假設(shè)模型與數(shù)據(jù)適配良好。

3 結(jié)果與分析

圖 1 初始與修正模型路徑圖Figure 1 Initial and correction model path map

表 2 林分生長結(jié)構(gòu)方程模型的擬合參數(shù)Table 2 Fitting parameters of structure equation model for the growth of forest

由表3所示：基巖裸露度對林分密度、樹種多樣性有正影響，對建群種生物量有負影響，其路徑系數(shù)分別為0.198、0.519、?0.722，即每當基巖裸露度變化1時，其林分密度、樹種多樣性和建群種生物量分別變化0.198、0.519、?0.722?；鶐r裸露度對樹種多樣性的的正影響較大，對建群種生物量的負影響最大。建群種物量受到土壤厚度、林分密度和物種多樣性的正影響，受到基巖裸露度的負影響，其路徑系數(shù)分別為0.258、0.096、0.084、?0.722，即當土壤厚度、林分密度、物種多樣性和基巖裸露度分別變化1時，建群種生物量分別變化0.258、0.096、0.084、?0.722。土壤厚度對建群種生物量有最大的正影響。林分總生物量受到樹種多樣性和建群種生物量的正影響，其路徑系數(shù)分別為0.126、1.014。林分總生物量的變化幾乎完全受林分建群種生物量的影響。

表 3 結(jié)構(gòu)方程模型標準化影響系數(shù)Table 3 Standard impact coefficient of SEM

河南省淅川縣巖溶區(qū)建群種生物量和林分總生物量主要受到基巖裸露度、土壤厚度、林分密度和物種多樣性的直接或間接影響。基巖裸露度對建群種生物量的總影響系數(shù)為?0.659，直接影響系數(shù)為?0.722，間接影響系數(shù)為0.063；對林分總生物量的總影響系數(shù)為?0.604，間接影響系數(shù)為?0.604，直接影響系數(shù)為0。因為基巖裸露度對林分總生物量直接影響系數(shù)可忽略，在結(jié)構(gòu)方程模型中將該路徑刪除。當基巖裸露度變化1時，建群種生物量和總生物量分別變化?0.659、?0.604。即隨著基巖裸露度增加，建群種生物量和林分總生物量下降。土壤厚度對建群種生物量的總影響系數(shù)或直接影響系數(shù)為0.258；對林分總生物量的總影響系數(shù)或間接影響系數(shù)0.262。隨著土壤厚度的增加，林分建群種生物量和總生物量增加。樹種多樣性對林分建群種生物量的總影響系數(shù)或直接影響系數(shù)為0.084，對林分總生物量影響系數(shù)為0.211，直接影響系數(shù)為0.126，間接影響系數(shù)為0.085。隨著樹種多樣性的增加建群種生物量和總生物量緩慢增加。林分密度對林分建群種生物量的總影響系數(shù)或直接影響系數(shù)為0.096，對林分總生物量的總影響系數(shù)或間接影響系數(shù)為0.098。建群種生物量對林分總生物量的影響系數(shù)為1.014，林分總生物量幾乎完全受建群種生物量的決定。巖溶區(qū)土壤瘠薄、土壤水肥條件較差、基巖裸露度較大，立地條件差，對植被的適應(yīng)性要求較高，建群種在立地競爭中占據(jù)絕對優(yōu)勢。巖溶區(qū)建群種生物量占林分總生物量的77.50%以上(表1)，提升建群種的生長條件，能夠直接提升林分的總生物量。

基巖裸露度對樹種多樣性的總影響系數(shù)或直接影響系數(shù)為0.519，對林分密度的總影響系數(shù)或直接影響系數(shù)為0.198，即隨著基巖裸露度的增加，物種的多樣性和林分密度增加。巖溶區(qū)基巖的裸露不僅僅造成地被的不連續(xù)，而且造成明顯的小生境，高度異質(zhì)的巖溶區(qū)物種多樣性隨著基巖裸露度增加而增加。

4 討論與結(jié)論

4.1 討論

河南省淅川縣巖溶區(qū)天然次生林基巖裸露度、土壤厚度、樹種多樣性、林分密度與建群種生物量和總生物量存在著復(fù)雜的耦合關(guān)系?；鶐r裸露度與林分密度為正相關(guān)關(guān)系，與樹種多性為極顯著正相關(guān)關(guān)系(P＜0.01)，與建群種生物量為極顯著負相關(guān)關(guān)系(P＜0.01)；土壤厚度與建群種生物量顯著正相關(guān)關(guān)系(P＜0.05)，樹種多樣性、林分密度與建群種生物量為正相關(guān)關(guān)系；建群種生物量與林分總生物量為顯著正相關(guān)關(guān)系(P＜0.01)。根據(jù)模型擬合結(jié)果表明淅川縣巖溶區(qū)主要影響林分生物量的觀測變量基巖裸露度和土壤厚度。基巖裸露度增加，造成地被破碎，土壤不連續(xù)，植物生存空間變少、生長的立地條件變差，樹種多樣性增加，林分密度增加，但樹木的胸徑樹高減少，林分生物量反而下降。土壤厚度是土壤肥力存在和植物生長的重要物質(zhì)基礎(chǔ)，是植被重建最基本的條件[26?27]。土壤厚度增加能夠促進植物的生長，增加林分生物量。

基巖裸露度和土壤厚度等作為立地質(zhì)量的評價指標，能反映立地條件的優(yōu)劣。一般隨立地條件的改善，林分的物種多樣性越來越高，結(jié)構(gòu)功能越來越好[28]。同時生態(tài)學理論認為群落的穩(wěn)定性與植物的多樣性關(guān)系密切，群落的物種多樣性是衡量生態(tài)系統(tǒng)穩(wěn)定性的指標[28?29]。淅川縣巖溶區(qū)立地條件提升時樹種多樣性會降低，可能是立地質(zhì)量改善促進了建群種個體的生長，在競爭中占據(jù)優(yōu)勢地位，導(dǎo)致其他物種減少，林分復(fù)雜程度降低，群落的穩(wěn)定性降低。巖溶區(qū)植被受到強烈的特殊地質(zhì)環(huán)境脅迫，植物群落簡單，多樣性低，并且?guī)r溶區(qū)生長植物既要有石生性、喜鈣性以及耐干旱的特點，又必須具有發(fā)達、強壯的根系[30]，所以巖溶區(qū)的適生種、優(yōu)勢種具有明顯的特殊性[31?32]。立地條件提升能顯著促進建群種的生長，使得林分復(fù)雜程度越來越低，多樣性降低。

生態(tài)系統(tǒng)的變化是一個復(fù)雜的過程，并不是簡單的一一對應(yīng)關(guān)系[15]。研究系統(tǒng)復(fù)雜的生態(tài)過程用傳統(tǒng)的方法有一定局限性，傳統(tǒng)的生長模型大多為直接模型，無法反映隱含的間接關(guān)系[33?34]。本研究用結(jié)構(gòu)方程模型來探討淅川縣巖溶區(qū)基巖裸露度、土壤厚度、林分密度、樹種多樣性和建群種生物量、林分總生物量之間的復(fù)雜關(guān)系。根據(jù)模型計算結(jié)果可知：因子間既有直接的影響關(guān)系，也有間接的關(guān)系。基巖裸露度對建群種生物量有直接作用，同時通過影響樹種多樣性和建群種生物量對林分生物量產(chǎn)生間接作用。土壤厚度、物種多樣性、林分密度對建群種生物量有直接作用，并且通過影響建群種生物量對林分總生物量產(chǎn)生間接作用。保護巖溶區(qū)生態(tài)環(huán)境，防止水土流失，保持土壤厚度以及穩(wěn)定或降低基巖裸露度能夠提高巖溶區(qū)次生林的生物量，增加干物質(zhì)積累。

本研究對象是淅川縣巖溶區(qū)天然次生林，由于巖溶區(qū)特殊的地質(zhì)水文結(jié)構(gòu)，林分的樹種多樣性低，穩(wěn)定性較差。同時所選取的樣地數(shù)量相對較少，模型的最終適配結(jié)果雖然存在一定的不足，但是也能反映出巖溶區(qū)林分生長與影響因子間的關(guān)系。合理改善立地條件，降低基巖裸露度和增加土壤厚度，及時撫育提升樹種多樣性，增加林分密度，有利于增加巖溶區(qū)生物量的積累，促進林分生長，逐步提高巖溶區(qū)林分的生態(tài)功能。

4.2 結(jié)論

基巖裸露度與林分生物量關(guān)系復(fù)雜?；鶐r裸露度能直接降低林分生物量，但又能通過增加樹種多樣性和林分密度降低間接提升林分生物量。土壤厚度、物種多樣性和林分密度直接影響林分的生物量。河南省淅川縣巖溶區(qū)林分生長與影響因子間的關(guān)系耦合結(jié)果表明：巖溶區(qū)建群種生物量和林分總生物量受到基巖裸露度的極顯著負影響(P＜0.01)，其影響系數(shù)分別為?0.659、?0.604；受到土壤厚度的顯著正影響(P＜0.05)，影響系數(shù)分別為0.258、0.262；受到樹種多樣性和林分密度正影響，影響系數(shù)分別為0.084、0.211，0.096、0.098。在巖溶區(qū)防治水土流失，穩(wěn)定土壤厚度及基巖裸露度，能改善林分結(jié)構(gòu)等立地條件，增加建群種和林分總生物量，提升巖溶區(qū)的生態(tài)環(huán)境質(zhì)量。