直線型線截抽樣在森林資源綜合監(jiān)測中的應用

羅仙仙

（泉州師范學院 數(shù)學與計算機科學學院，福建 泉州 362000）

線截抽樣是稀疏總體的調(diào)查方法，是一種有放回不等概抽樣方法[1-2]，是在蒲豐隨機投針問題（Buf fon’s needle problem）基礎上推導而來，具有完善的理論基礎，1941年由Canfield首次提出，用來估計植被分布密度[3]， 先后用于采伐剩余物的調(diào)查[4]， 道路長度估計[5]， 可燃物抽樣調(diào)查[6]， 法國國家森林清查中樹籬長度的估計[7]、美國堪薩斯州林帶面積估計[8]和干旱地區(qū)生物量估計[9]。近年來，線截抽樣技術在冠空隙[10]、森林生物量[11]等方面得到應用。劉素青等[11-12]利用線截抽樣進行林分株數(shù)和林分蓄積量調(diào)查。目前，線截抽樣被廣泛應用于林地粗木質(zhì)物殘體、倒木、野生動物、漁業(yè)和其他流動總體抽樣調(diào)查[2]。線截抽樣截線的形狀類型及其抽樣協(xié)議見文獻[1-2]。粗木質(zhì)物殘體（coarse woody debris，CWD）是指完好的和處于不同腐解時期的直徑（通常指粗頭部分）D≥10 cm，長度L≥1 m的倒木、枯立木、大凋落枝，以及直徑D≥10 cm，長度L＜1 m的根樁和直徑D＞1 cm的地下粗根殘體[13]。目前，對粗木質(zhì)物殘體的研究大多為生態(tài)功能研究[14-16]，尚未見有關調(diào)查方法及統(tǒng)計估計方法的研究。在森林資源監(jiān)測逐漸向與林業(yè)可持續(xù)發(fā)展相適應的森林生態(tài)系統(tǒng)監(jiān)測發(fā)展過程中，許多林業(yè)發(fā)達國家均增加了粗木質(zhì)物殘體的調(diào)查與監(jiān)測[17]。如美國森林資源清查與監(jiān)測體系采用Y型線截抽樣進行粗木質(zhì)殘體的調(diào)查[18]。為適應森林生態(tài)系統(tǒng)監(jiān)測需求，利用直線型線截抽樣原理對試驗區(qū)的倒木進行調(diào)查，驗證線截抽樣調(diào)查方法的可行性。同時，探討不同截線長度對調(diào)查精度的影響，從而確定最佳的截線長度，以期為森林資源綜合監(jiān)測地面樣地設計提供抽樣技術指導。

1 倒木調(diào)查的直線型線截抽樣原理

線截抽樣的截線布設方式有2種：隨機布設和系統(tǒng)布設，對總體參數(shù)的估計方法可分為條件估計方法和無條件估計方法。條件估計是指截線在隨機布設下的估計；無條件是指線截抽樣中每條截線的方向在區(qū)間[0，π]內(nèi)被隨機一致布設。因為無法從相對精度的角度確定條件估計與無條件估計2種方法的優(yōu)劣，從實踐角度上說，盡量采取能夠減少外業(yè)工作量的估計方法[2]。本研究采用是無條件估計方法對倒木的蓄積、長度、條數(shù)進行估計。假定τy為總體總值，則總體密度為，總體均值為，其中A為總體區(qū)域面積，N為區(qū)域內(nèi)總體單元數(shù)。L為直線型截線的長度。如圖1所示：μk為倒木，Dk為大頭直徑，dk為小頭直徑，lk為倒木長度。則個體μk構成最小封閉凸面的圍長ck≈2lk+dk+Dk，因此，個體μk的包含概率為 πk=（ck·L）/（πA）=（2lk+dk+Dk）L/（πA）。 第s條截線屬性為 yk的 τy無條件估計為：

圖1 直線型截線與直徑為dk，Dk，長為lk倒木相交Figure1 A straight line transect intersects a log with diameters dk， Dk， and length lk

式（1）中Ls為第s條截線。

當?shù)鼓局睆较鄬﹂L度較小時即dk，Dk遠小于lk時，式（1）近似公式為：

1.1 倒木蓄積總值無條件估計

若有m條截線，則所有倒木蓄積總值τy的HT（Horvitz-Thompson，赫維茨-湯普森）估計值為：

1.2 倒木蓄積密度無條件估計

m條截線所截倒木蓄積密度無條件估計為：

1.3 倒木長度無條件估計

當估計區(qū)域Φ中的倒木總長度時，此時，yk=lk，第s條截線所截倒木長度估計為：

式（7）中：ns為第s條截線所截倒木的條數(shù)。

m條截線所截倒木總長度無條件估計為：

1.4 倒木總條數(shù)無條件估計

當估計區(qū)域Φ中的倒木總條數(shù)時，此時，yk＝1，區(qū)域Φ中第s條截線的倒木條數(shù)無條件估計為：

m條截線所截倒木總條數(shù)無條件估計為：

1.5 方差估計

1.6 估計精度

τy無條件估計的估計精度為：

λy無條件估計的估計精度為：

1.7 區(qū)間估計

τy無條件估計的估計區(qū)間為：

λy無條件估計的估計區(qū)間為：

1.8 樣線數(shù)的確定

假定個體分布為正態(tài)、獨立、同分布。無條件估計的樣線數(shù)為：

式（17）中：E為相對誤差限。

2 研究區(qū)概況

1987年，北京林業(yè)大學和吉林省汪清林業(yè)局合作在金溝嶺林場建立檢查法試驗項目。項目的宗旨是預定和提交林分級尺度上森林生態(tài)系統(tǒng)；增強林地的承載力和生產(chǎn)力；提高森林經(jīng)營水平。試驗區(qū)屬長白山系， 位于 43°22′N， 130°10′E， 總面積為340.9 hm2；各個大區(qū)分5個小區(qū)，共15個小區(qū)。本研究試驗選擇在檢查法試驗I大區(qū)4小區(qū)，有關I大區(qū)樣地設計基本情況見表1。

表1 I大區(qū)各小區(qū)樣地面積和樣地數(shù)量表Table1 Area and quantity of plots in the first section

3 研究方法

2008年8月，在檢查法的試驗I大區(qū)4小區(qū)進行線截抽樣的試驗研究。I大區(qū)4小區(qū)共有22個樣地，各個樣地大小均為0.04 hm2，樣地布設方式見圖2。其中東西方向為坡上坡下方向，為減少坡度的影響，在各個樣地西樁處布設1條南北方向測繩，以西樁為中心，控制測繩在南北方向各50 m長，南方向起點為0 m，測繩總長100 m，各截線長度統(tǒng)計示意圖所圖3所示，共有22條截線。為探討直線型線截抽樣調(diào)查方法的可行性，只對大頭直徑Dk≥10 cm，長度L≥1 m的倒木進行調(diào)查。調(diào)查對象為與截線完全相交或部分相交的倒木，調(diào)查的指標為倒木的大頭直徑、小頭直徑、長度、分解等級。分解等級分為新近死亡、開始分解、大量分解、幾乎完全分解4個等級。同時記錄倒木截點處位置，便于不同截線長度線截抽樣的倒木對象的統(tǒng)計分析。

圖2 I大區(qū)4小區(qū)樣點分布圖Figure2 Plots distribution in the fourth small section of the first section

圖3 各截線長度統(tǒng)計示意圖Figure3 Statistical hint of different transect length

4 結果與分析

利用前面式（4）（6）（13），不同截線長度倒木蓄積量無條件估計的抽樣效果見表2。根據(jù)公式（10）（16）不同截線長度倒木條數(shù)無條件估計的抽樣效果見表3。不同截線長度所截取倒木條數(shù)及估計值見表4。截線長度是影響精度的主要因子[19]，通過不同截線長度對倒木估計的影響，從而確定最優(yōu)的截線長度。

4.1 倒木蓄積總量的估計精度與截線的長度關系

如圖4所示：倒木蓄積總量的估計精度隨截線的長度增加而增加。截線長度是影響精度的主要因子[19]。當截線長度處于區(qū)間[10，50]m時，隨截線長度的增加，抽樣精度明顯提高；當截線長度達到50 m，增加截線長度對抽樣精度影響不大，即當截線長度為50 m時，抽樣精度為90.05%，截線長度取50 m直線型線截抽樣既滿足抽樣精度要求，同時也減少外業(yè)工作量。

4.2 倒木總條數(shù)估計值與截線的長度關系

如圖5所示：倒木總條數(shù)估計值隨截線的長度增加先增大后減少，當截線長度區(qū)間為[50，100]m時，倒木總條數(shù)估計值在（5000，6000）相對穩(wěn)定變化。

4.3 倒木蓄積密度與截線的長度關系

如圖6所示：倒木蓄積密度隨截線的長度增加先增大后減少，當截線長度區(qū)間為[50，100]m時，倒木蓄積密度在（25，30）相對穩(wěn)定變化。

4.4 樣本倒木條數(shù)與截線的長度關系

如圖7所示：所選樣本倒木條數(shù)隨截線的長度增加而近線性增加，說明倒木在研究區(qū)域中的分布情況相對均勻。

綜上各指標與截線長度的關系，在一定區(qū)域內(nèi)對倒木進行調(diào)查時，截線長度取50 m既可以滿足抽樣精度要求，又可以減少外業(yè)工作量，因此是直線型線截抽樣的最優(yōu)截線長度。此時，I大區(qū)4小區(qū)倒木蓄積總量為536.85 m3，蓄積密度為27.53 m3·hm-2，倒木總條數(shù)為5887條，抽樣精度達90.05%（可靠性為95%）。

表2 不同截線長度倒木蓄積量無條件估計的抽樣效果Table2 Sampling effect of different transect length on the volume of log

表3 不同截線長度倒木條數(shù)無條件估計的抽樣效果Table3 Sampling effect of different transect length on the number of log

5 結論與討論

5.1 結論

本研究在檢查法的基礎上，根據(jù)現(xiàn)有森林生態(tài)系統(tǒng)監(jiān)測需求，應用直線型線截抽樣原理對實驗區(qū)進行了倒木的調(diào)查。探討了不同截線長度對精度的影響，確定了倒木調(diào)查的最佳截線長度為50 m，此時，倒木蓄積總量為536.85 m3，蓄積密度為27.53 m3·hm-2，倒木總條數(shù)為5887條，抽樣精度達90.05%（可靠性為95%）。

5.2 討論

森林類型和干擾歷史決定了粗木質(zhì)物殘體數(shù)量。本研究只對大頭直徑Dk≥10 cm，長度L≥1 m的倒木進行系統(tǒng)布設的直線型線截抽樣調(diào)查與研究，對其他粗木質(zhì)物殘體采用其他截線類型線截抽樣可進一步研究。然而在粗木質(zhì)物殘體隨機分布的立地中，采用復雜類型線截抽樣不如采用單條直線型線截抽樣[18]。

表4 不同截線長度所截取倒木條數(shù)及估計值Table4 Number of log and its estimation of different transect length

圖4 倒木蓄積總量的估計精度與截線長度的關系Figure4 Relationship between total volume of log and transect length

圖5 倒木總條數(shù)估計值與截線長度的關系Figure5 Relationship between total number of dead tree and transect length

調(diào)查中沒有考慮坡度的影響；截線南北方向不一定十分準確，這在一定程度上影響抽樣的概率。

圖6 倒木蓄積密度與截線長度的關系Figure6 Relationship between volume density of dead tree and transect length

圖7 樣本倒木條數(shù)與截線的長度關系Fiure 7 Relationship between numbers of dead tree of dead tree and transect length

本研究調(diào)查是在檢查法的基礎上進行的，截線的布設方式是系統(tǒng)布設的，截線的條數(shù)與檢查法的樣地數(shù)一致。截線數(shù)量對倒木估計精度的影響可進一步研究。

此次調(diào)查研究中，所采用的截線長度是等長的?？梢愿鶕?jù)粗木質(zhì)物殘體的分布狀況，采用不等長的截線長度進行調(diào)查，同時亦可采用多階線截抽樣進行其他監(jiān)測指標的調(diào)查。

本研究所采用方法只針對研究區(qū)域進行，對于其他區(qū)域或其他森林類型的倒木的調(diào)查直線型截線長度是否50 m就能達到90.05%精度有待試驗與研究。

[1]GREGOIRE T G， VALENTINE H T.Line intersect sampling： EII-shaped transects and multiple intersections [J].Environ Ecol Stat， 2003， 10 （2）： 263-279.

[2]GREGOIRE T G， VALENTINE H T.Sampling Strategies for Natural Resources and the Environment[M].Boca Raton：Chapman&Hall/CRC， 2008： 279-326.

[3]CANFIELD R H.Application of the line interception method in sampling range vegetation [J].J For， 1941， 39 （4）：34-40.

[4]WARREN W G， OLSEN P F.A line intersect technique for assessing logging waste [J].For Sci， 1964， 10 （3）： 267-276.

[5]MATéRN B.A method of estimating the total length of roads by means of a line survey [J].Stud For Suec， 1964， 18：68-70.

[6]BROWN T K.A plannar intersect method for sampling fuel volume and surface area [J].For Sci， 1971， 17 （1）： 96-102.

[7]CHEVROU R.Inventaire des haies [J].Rev For Franc， 1973， 25： 47-53.

[8]HANSEN M H.Line intersect sampling of wooded strips [J].For Sci， 1985， 31 （2）： 282-288.

[9]de VRIES P G.Sampling Theory for Forest Inventory [M].Berlin： Springer-Verlag， 1986： 242-279.

[10]BATTLES J J， DUSHOFF J G， FAHEY T J.Line intersect sampling of forest canopy gaps [J].For Sci， 1996， 42 （2）：131-138.

[11]NELSON R， SHORT A， VALENTI M.Measuring biomass and carbon in delaware using an airborne profiling LIDAR[J].Scandinavian J For Res， 2004， 19 （6）： 500-511.

[11]劉素青，洪家勝，彭世揆.林分蓄積量調(diào)查的線截抽樣原理[J].東北林業(yè)大學學報，1999，27（2）：60-64.LIU Suqing， HONG Jiasheng， PENG Shikui.The principle of line intersect sampling for surveying stand volume[J].J Northeast For Univ， 1999， 27 （2）： 60-64.

[12]劉素青，洪家勝，彭世揆.線截抽樣與線狀物體調(diào)查原理研究[J].江西農(nóng)業(yè)大學學報，1998，20（3）：376-380.LIU Suqing， HONG Jiasheng， PENG Shikui.Line intersect sampling and surveying on linearshape object （road net，forest belt etc.） [J].Acta Agric Univ Jiangxi， 1998， 20 （3）： 376-380.

[13]閆恩榮，王希華，黃建軍.森林粗死木質(zhì)殘體的概念及其分類[J].生態(tài)學報，2005，25（1）：158-167.YAN Enrong， WANG Xihua， HUANG Jianjun.Concept and classification of coarse woody debris in forest ecosystems[J].Acta Ecol Sin， 2005， 25 （1）： 158-167.

[14]侯平，潘存德.森林生態(tài)系統(tǒng)中的粗死木質(zhì)殘體及其功能[J].應用生態(tài)學報，2001，12（2）：309-314.HOU Ping， PAN Cunde.Coarse woody debris and its function in forest ecosystem [J].Chin J Appl Ecol， 2001， 12 （2）：309-314.

[15]劉惠英，張思玉，吉霞.粗死木質(zhì)殘體的水土保育功能[J].世界林業(yè)研究，2004，17（3）：25-28.LIU Huiying， ZHANG Siyu， JI Xia.The soil and water conservation function of coarse woody debris [J].World For Res，2004， 17 （3）： 25-28.

[16]趙玉濤，余新曉，程根偉，等.粗死木質(zhì)殘體的水文生態(tài)功能[J].山地學報，2002，20（1）：12-18.ZHAO Yutao， YU Xinxiao， CHENG Genwei， et al.A slighting tache in field of forest hydrology research： hydrological effects of coarse woody debris （CWD） [J].J Mount Res， 2002， 20 （1）： 12-18.

[17]羅仙仙，亢新剛.森林資源綜合監(jiān)測研究綜述[J].浙江林學院學報，2008，25（6）：803-809.LUO Xianxian， KANG Xingang.Progress in research on the comprehensive monitoring of forest resources [J].J Zhejiang For Coll， 2008， 25 （6）： 803-809.

[18]du CROS R T， LOPEZ S.Preliminary study on the assessment of deadwood volume by the French national forest inventory [J].Ann For Sci， 2009， 66 （3）： 1-10.

[19]WOLDENDORP G， KEENAN R J， BARRY S， et al.Analysis of sampling methods for coarse woody debris [J].For Ecol Manage， 2000， 198 （1/3）： 133-148.