2008特大冰凍災害后大明山常綠闊葉林林冠結構動態

周曉果, 溫遠光,2,3,*, 朱宏光,2,3, 王 磊

1 廣西大學林學院，廣西木林生態與保育重點實驗室培育基地，南寧 530004 2 亞熱帶農業生物資源保護與利用國家重點實驗室,南寧 530004 3 廣西友誼關森林生態系統定位觀測研究站,憑祥 532600

2008特大冰凍災害后大明山常綠闊葉林林冠結構動態

周曉果1, 溫遠光1,2,3,*, 朱宏光1,2,3, 王 磊1

1 廣西大學林學院，廣西木林生態與保育重點實驗室培育基地，南寧 530004 2 亞熱帶農業生物資源保護與利用國家重點實驗室,南寧 530004 3 廣西友誼關森林生態系統定位觀測研究站,憑祥 532600

林冠結構是研究森林生態系統眾多關鍵生態功能和過程的重要參數,常綠闊葉林是亞熱帶林區具有代表性的森林類型,對其林冠結構及動態特征的研究還很不深入。在廣西大明山中山區選擇了一個斜坡水平長200 m、寬160 m的典型坡面,在整個坡面建立了80個20 m×20 m的樣地,將樣地均勻分為5個坡段,每個坡段包含16個連續的樣地,在2009—2012年的生長季測定了林冠高度(CH)、林冠體積(CV)、林冠覆蓋度(CC)、林冠上/下冠蓋比(HLr)和林冠葉面積指數(LAI),分析了各林冠結構指標的坡位及年際動態,揭示了亞熱帶常綠闊葉林的林冠結構特征及短期動態規律。研究結果表明,大明山常綠闊葉林林冠結構的一般特征是：平均CH(12.09±0.05) m,平均CV(2642.51±278.33) m3(每 400 m2樣地),平均CC(59.90±3.29)%,平均HLr 2.48±0.23,平均LAI 2.00±0.06。大明山常綠闊葉林的林冠結構存在多層性,上層林冠覆蓋度平均為42.20%,中層為30.35%,下層為18.05%。大明山常綠闊葉林的林冠結構存在坡面和年際差異,坡面變異系數為CV(29.84%—55.89%)> HLr(32.90%—53.52%)> LAI(22.48%—43.89%)> CC(16.61%—25.74%)> CH(8.26%—12.77%)；年際變異系數為HLr(47.33%—57.00%)> CV(39.70%—49.06%)> LAI(21.58%—48.13%)> CC(20.35%—24.15%)> CH(9.19%—12.59%),表明CH有較強的穩定性。林冠LAI存在明顯的坡面尺度效應,即向下順坡每滑動100 m冠層LAI升高0.34。坡位對CH、HLr有顯著(P=0.022)和極顯著(P<0.001)影響；年份對HLr有顯著影響(P=0.013),對CV和CC有極顯著影響(P<0.001)；坡位×年份對CV和LAI的交互作用顯著(P=0.016,P=0.017)。回歸分析發現樹冠面積與林木胸徑呈極顯著的線性關系。此研究結果表明大明山常綠闊葉林冠層高度較低、林冠體積較小、林冠覆蓋度不高、上/下冠蓋比和葉面積指數偏小,這與研究區域的海拔較高(934—1223 m),土層淺薄(30—45 cm)以及經常受到冰凍災害(特別是2008年的特大冰凍災害)的影響有關,是山地常綠闊葉林樹冠結構與山地環境條件長期適應的結果。

常綠闊葉林；葉面積指數；林冠覆蓋度；林冠體積；上/下冠蓋比；動態

林冠是指森林中地表以上的所有葉片、枝條、小枝、各種附生生物及其枯死殘留物的總和[1- 4],林冠結構是林冠要素(如冠高、葉片、枝條等)在空間和時間上的組成、結構及動態[5- 6]。林冠高度(CH)、林冠體積(CV)、林冠覆蓋度(CC)、上/下冠蓋比(HLr)、林冠葉面積指數(LAI)等都是林冠結構的重要參數[5- 7],其三維空間的異質性及垂直分層使林冠蘊含了全球40% 的陸地生物多樣性[8],被視為“最后的生物前沿(the last biological frontier)”和生物圈中物種豐富卻最鮮為人知的生境之一[9]。林冠結構及其變化直接控制著森林生態系統與大氣的物質和能量交換,與森林小氣候、森林水文、森林養分循環密切相關[10-14],成為森林生態系統結構、功能及關鍵生態學過程的重要組分[6- 7],對森林更新、生物多樣性、生物生產力、碳固持及全球氣候變化有著重要的影響[1,15-16],林冠結構的多樣性和生態過程仍被視為“科學黑箱”[4,17],成為全球變化生態學研究的熱點[2,18]。

目前,有關林冠結構特征及動態的研究,主要集中在熱帶和溫帶地區的同齡林和人工林[19-24],研究的參數主要是LAI、林冠開度等[20- 24]少數指標,大都是對參數的季節變化[20-24]和坡向、坡位等立地特征差異的簡單比較上[24-26]。2008年特大冰凍災害后,有關森林林冠受損的研究明顯增多[27- 28],近年有學者對六盤山南坡華北落葉松人工林冠層LAI的坡面尺度效應進行了研究[19],但對災后常綠闊葉林林冠結構及動態的研究仍十分有限[29-30],尤其缺乏對原生性頂極森林及坡面較大面積林冠結構動態變化的研究[29]。

本研究以廣西大明山國家級自然保護區原生性常綠闊葉林為對象,對2008年特大冰凍災害后常綠闊葉林的林冠結構及動態進行了連續4年的監測,以期為深入認識南亞熱帶山地常綠闊葉林的林冠結構動態變化規律,為該地區林分、坡面、景觀和區域尺度上碳、水分和通量等方面的模擬提供基礎數據,也為豐富林冠學理論提供資料。

1 研究區概況

大明山國家級自然保護區(23°10′—23°38′ N,108°18′—108°45′ E),地處廣西中南部,是北熱帶與南亞熱帶的分界,北回歸線橫貫中部,屬南亞熱帶季風氣候區。年均溫為15.1℃,7月份平均氣溫21.9℃,1月份平均氣溫5.8℃,極端高溫為28.6℃,極端低溫-6.0℃；日均溫≥10℃的積溫為5047.7℃；年均降水量為2630.3 mm[31]。大明山地理位置獨特,地層古老,地貌發育,形成了復雜多樣的生態環境,孕育出豐富的生物多樣性,保存了較完好的原生性森林,共有維管植物234科918屬2374種[32]。茂密的天然森林資源使大明山成為周邊地區重要的物種基因庫、水庫、碳庫、氧庫,成為區域生態安全的重要屏障[33]。

2 研究方法

2.1 樣地構建

于2009年5月,在大明山常綠闊葉林中建立了3.2 hm2的長期監測樣地。樣地面積為200 m×160 m, 海拔934—1223 m之間,坡度范圍30—35°。采用全站儀將3.2 hm2樣地劃分成80個20 m×20 m的樣方,再把每個樣方劃分成4個10 m×10 m和16個5 m×5 m的小樣方,在每個20 m×20 m及10 m×10 m樣方的4個角用PVC管作永久標記,并標明樣方號及相應的坐標。由下至上將樣地分為5個連續的坡位,每個坡位包括16個20 m×20 m的樣方(表1)。調查區域森林植被的優勢種是云貴山苿莉(Huodendronbiaristatum,重要值: 33.1±37.1,下同)、羅浮槭(Acerfabri,20.6±11.2)、栓皮木姜子(Litseasuberosa,18.5±16.6)、刨花潤楠(Machiluspauhoi,17.2±12.3)、羅浮柿(Diospyrosmorrisiana,9.2±9.3)、黃杞(Engelhardtiaroxburghiana,9.0±8.8)和天目紫莖(Stewartiagemmata,8.7±10.5)。

表1 不同坡位常綠闊葉林樣地的基本特征

2.2 林冠調查

在樣地構建好的基礎上,于2009—2012年生長季,以5 m×5 m的小樣方為基本調查單元,對每個胸徑≥1.0 cm的木本植物掛鋁牌標記,記錄樹木的編號、物種名、坐標、胸徑、樹高、枝下高、冠幅等信息,每次調查需要25—30 d。同時采用CI- 110型數字植物冠層圖象分析儀(CID Inc.,Camas, Washington, USA),在每個樣方中心位置采集冠層影像。

2.3 數據處理與統計分析

大明山常綠闊葉林林冠結構圖解見圖1。將每個樣方的林木按樹高分上(樹高>8 m)、中(4—8 m)、下(0—4 m)3層。林冠高度(CH)以每個樣方最上層林木的平均高代表。林冠體積(CV)以每個樣方所有林木的樹冠體積之和表示,其中林木的樹冠體積采用圓柱體公式計算。林冠覆蓋度(CC)為每個樣方所有林木的垂直投影占樣方面積的百分比。上/下冠蓋比(HLr)為上層與下層林木林冠覆蓋度的比值。采用冠層分析儀自帶的冠層分析軟件對冠層影像進行分析,計算每個樣方的葉面積指數(LAI)。

采用單因素方差分析(one-way ANOVA)檢驗各調查年度、各坡位CH、CV、CC、HLr、LAI的差異顯著性,采用兩因素方差分析(two-way ANOVAs)檢驗坡位、調查時間及其互作對CH、CV、CC、HLr、LAI等林冠結構特征值影響的差異程度,并采用LSD法進行多重比較。

采用SPSS 20.0軟件對所有數據進行處理和統計分析,顯著性水平設置為P<0.05。數據繪圖由Sigmaplot 11.0軟件完成。

圖1 大明山常綠闊葉林的林冠結構圖解Fig.1 An illustration of the forest canopy structure in evergreen broadleaved forest at Damingshan

3 結果與分析

3.1 常綠闊葉林林冠高度、林冠體積及其變化

2009—2012年,大明山常綠闊葉林的林冠高度(CH)平均值變化在12.00—12.21 m之間(平均值為12.09 m),總體呈現微弱增加趨勢(表2)。大明山常綠闊葉林不同坡位的CH存在顯著差異(P=0.022),然而,4年內,年份對CH無顯著影響,坡位×年份的交互作用也不顯著(P=0.686、P=0.656)(表2,表3),其中2009年和2010年5個坡位的CH均無顯著差異,2011年下坡的CH顯著高于中坡,2012年顯著高于中下坡、中坡和上坡,其余坡位間差異不顯著(圖2)。2009—2012年,大明山常綠闊葉林不同坡位的林冠體積(CV)差異不顯著(P=0.964),每400 m2樣地的CV變化在2173.57—3283.90 m3之間,平均值為2624.51 m3,年份對CV有極顯著影響(P<0.001),坡位×年份的交互作用顯著(P=0.016)(表2，表3),2011年和2012年的CV顯著高于2009年和2010年(表2)。2009年和2010年5個坡位間CV無顯著差異,但2011年下坡和中下坡顯著高于中坡,2012年下坡顯著高于上坡(圖2)。CV的坡面變異性(29.84%—55.89%)和年際變化性(39.70%—49.06%)均明顯高于CH(8.26%—12.77%,9.19%—12.59%)。

3.2 常綠闊葉林林冠覆蓋度、上/下冠蓋比及其變化

2009—2012年,大明山常綠闊葉林林冠覆蓋度(CC)的平均值變化在53.55%—66.92%,平均值為59.90%,呈現逐漸恢復增長趨勢,坡位對CC影響不顯著(P=0.103),年份對CC有極顯著影響(P<0.001),坡位×年份的交互作用不顯著(P=0.380)(表2,表3)。2011年和2012年的CC顯著高于2009年和2010年,2012年比2011年出現輕微下降,但差異不顯著(表2)。大明山常綠闊葉林的CC在不同的坡位有所變化,2009年和2010年,中坡的CC顯著高于下坡和中下坡,2012年中坡顯著高于中上坡,2011年沒有差異(圖3)。

表2 2009—2012年大明山常綠闊葉林林冠結構特征均值及其變化

CH: 林冠高度canopy height；CV: 林冠體積canopy volume；CC: 林冠覆蓋度canopy cover；HLr: 上/下冠蓋比ratio of high to low cover；LAI: 葉面積指數leaf area index; 不同小寫字母表示該指標在不同年份間差異顯著(P<0.05);表中數據為平均值±標準誤(n=80)

表3 坡位和年份及其交互作用對常綠闊葉林林冠結構的影響

Table 3FandPvalues for the effects of slope position, year and their interaction on the forest canopy structures in evergreen broadleaved forest

因子Factorsdf林冠高度CH林冠體積CV林冠覆蓋度CC上/下冠蓋比HLr葉面積指數LAIFPFPFPFPFP坡位Position42．890．0220．150．9641．950．1037．41<0．0012．400．050年份Years30．490．68617．51<0．00119．40<0．0013．670．0132．250．083坡位×年份Position×years120．790．6562．120．0161．080．3801．230．2652．100．017

圖2 2009—2012年大明山常綠闊葉林不同坡位的林冠高度和林冠體積特征Fig.2 The forest canopy height and canopy volume in different positions in evergreen broadleaved forest at Damingshan during 2009—2012

2009—2012年,大明山常綠闊葉林林冠的上/下冠蓋比(HLr)平均值為2.10—3.14之間,呈現先降后升的趨勢,坡位和年份對HLr有極顯著和顯著影響(P<0.001,P=0.013),但兩者的交互作用不顯著(P=0.265)。2012年的HLr顯著高于前3年,而前3年間無顯著差異(表2)。2009—2012年,不同坡位的HLr存在明顯變化。2009年下坡的HLr顯著高于其它坡位,中下坡顯著高于中上坡和上坡；2010年下坡和中下坡的HLr顯著高于中上坡和上坡；2011年下坡的HLr顯著高于上坡；2012年中下坡的HLr顯著高于其它坡位,下坡顯著高于中坡、中上坡和上坡(圖3)。

圖3 2009—2012年大明山常綠闊葉林不同坡位的林冠覆蓋度和上/下冠蓋比特征Fig.3 The forest canopy cover and ratio of high to low cover in different positions in evergreen broadleaved forest at Damingshan during 2009—2012

3.3 常綠闊葉林林冠LAI及其變化

圖4 2009—2012年大明山常綠闊葉林不同坡位林冠LAI比較 Fig.4 LAI in different positions in evergreen broadleaved forest at Damingshan during 2009—2012

2009—2012年,大明山常綠闊葉林林冠葉面積指數(LAI)平均值變化在1.88—2.15之間(平均值2.00 m2/m2),總體呈現微弱增加趨勢,坡位和年份對LAI有一定影響,但均不顯著(P=0.050,P=0.083),但兩者的交互作用影響顯著(P=0.017)(表3)。其中2011年的林冠LAI顯著高于2009年,其它年份間差異不顯著(表2)。2009—2012年,大明山常綠闊葉林不同坡位的林冠LAI隨年份不同而有所變化。2009年下坡的林冠LAI顯著高于上坡,2011年為中坡顯著高于其它坡位,2012年則是中下坡顯著高于其它坡位,2010年各坡位差異不顯著(圖4)。

3.4 常綠闊葉林林冠垂直結構及其變化

根據三個層次的劃分統計結果,大明山常綠闊葉林林冠垂直結構呈現上層林冠(>8 m)(林冠覆蓋度42.20%)> 中層林冠(4—8 m)(30.35%)> 下層林冠(0—4 m)(18.05%),上層林冠覆蓋度極顯著高于中層林冠和下層林冠,中層林冠極顯著高于下層林冠(P<0.001)。

2009—2012年,大明山常綠闊葉林林冠垂直結構變化見圖5。從圖5可以看出,不同年份不同坡位的林冠垂直結構存在一定的變化。2009年,下層林冠覆蓋度以下坡最小,上坡最大,上坡顯著高于下坡和中下坡,中坡和中上坡顯著高于中下坡；中層林冠覆蓋度也是以下坡最小和上坡最大,上坡和中坡顯著高于下坡和中下坡,其它坡位差異不顯著；上層林冠的覆蓋度在不同坡位間差異均不顯著(圖5)。2010年,下層林冠覆蓋度為中坡、中上坡和上坡顯著高于下坡和中下坡；中層林冠覆蓋度為中坡和上坡顯著高于下坡和中下坡；上層林冠的覆蓋度在不同坡位間差異也不顯著(圖5)。2011年,下層林冠和中層林冠的林冠覆蓋度都是以上坡的最高,下坡最低,上層林冠則以下坡最高,中坡最低,但3個林冠層次在不同坡位林冠覆蓋度差異均不顯著(圖5)。2012年,下層林冠覆蓋度的坡位變化與2010年的相似；中層林冠覆蓋度除中坡顯著高于下坡外,其余坡位間差異不顯著；上層林冠覆蓋度為下坡、中下坡和中坡顯著高于上坡(圖5)。

3.5 樹冠面積與林木胸徑的關系

2009—2012年,大明山常綠闊葉林樹冠面積與林木胸徑有極顯著的相關性(P<0.001)(圖6),表明樹冠面積隨著林木胸徑的增加而增加,且隨著恢復時間的增加其增幅加大,以2011年的增幅最大(圖6)。

4 討論與結論

4.1 大明山常綠闊葉林林冠結構的基本特征

本研究表明,大明山常綠闊葉林林冠結構特征是：平均CH為(12.09±0.05)m,平均CV為(2642.51±278.33)m3(每 400 m2樣地),平均CC為(59.90±3.29)%,平均HLr為2.48±0.23,平均LAI 2.00±0.06。這一結果表明,大明山常綠闊葉林的冠層高度較低,林冠體積較小,林冠覆蓋度不高,上/下冠蓋比和葉面積指數偏小,這與研究區域的海拔較高(934—1223 m),土層淺薄(30—45 cm),以及經常受到冰凍災害(特別是2008年的特大冰凍災害)的影響有關。大明山常綠闊葉林的林冠結構存在多層性,上層林冠覆蓋度平均為42.20%,中層為30.35%,下層為18.05%,但與亞熱帶森林林冠的多層性和葉面積指數較大[29]的特點不完全一致。本研究中觀察到的LAI變化在1—6 之間,但平均值只有2.0左右,大明山常綠闊葉林林冠的LAI明顯小于哀牢山常綠闊葉林(LAI為5.60)[29]、銀瓶山常綠闊葉林(2.56)[34],而與桉樹人工林相近(1.5—2.5)[35]。此種林冠結構特征的形成不僅與2008年的特大冰凍災害有關,還與災后受損林木的逐漸死亡和氣候波動有關。據災后(2009年)的調查,在3.2 hm2固定樣地中,DBH≥1 cm的林木受災率為51.8%[36],樣地有85個林隙,其中68.40%的林隙形成于2008—2009年,擴展林隙和冠空隙的面積分別占71.70%和52.90%,每個擴展林隙的面積為200—600 m2[37]；同時,由于山地氣候冰凍干擾頻繁,形成狹窄的圓柱形樹冠[31],這是造成大明山常綠闊葉林CH、CV、LAI偏小的主要原因。誠然,大明山常綠闊葉林林冠LAI較小可能還與CI- 110的測定結果偏小有關[21]。

圖5 2009—2012年大明山常綠闊葉林林冠垂直結構變化Fig.5 The forest canopy vertical structures in evergreen broadleaved forest at Damingshan during 2009—2012

圖6 2009—2012年大明山常綠闊葉林林木胸徑與樹冠面積的關系Fig.6 The relationship between DBH and crown area in evergreen broadleaved forest at Damingshan during 2009—2012

4.2 災后林冠結構的坡位變化

有研究表明,華北落葉松人工林受坡面環境條件變化的影響,林冠LAI存在著坡面變化與尺度效應,在一個水平長為 398.2 m、寬度為 30 m 的坡面上,LAI變化在 2.66—3.49之間(平均3.11),LAI總體上隨坡頂向下的坡長增加呈波動性增大趨勢,在坡面中部(水平坡長188.45 m時)達到最大,之后稍微減小,森林冠層LAI存在著坡面空間尺度效應,即向下順坡每滑動100 m冠層LAI升高0.12[19]。關于原生性常綠闊葉林林冠LAI的坡面變化尚未見報道。本研究表明,在一個水平坡長為200 m、寬度為160 m的坡面上,常綠闊葉林林冠LAI也存在著一定的坡面差異,且有顯著的坡位×年份的交互作用(P=0.017),表明LAI的坡面變化規律會隨著年份的不同而不同。在本研究中,2009年為下坡的林冠LAI顯著高于上坡,2010年坡位間無顯著差異,2011年為中坡顯著高于其它坡位,2012年為中下坡顯著高于其它坡位。常綠闊葉林林冠LAI也存在坡面尺度效應,且在坡面中部達到最大,之后稍微減小,即向下順坡每滑動100 m冠層LAI升高0.34,與華北落葉松人工林林冠LAI的尺度效應相似,但強度更強。這與原生性常綠闊葉林的林冠結構和空間異質性更高有關。因此,對于原生性常綠闊葉林更應注意LAI的空間和年際變化,提高生態模型中LAI 這個參數的準確性在未來研究中值得引起注意。

常綠闊葉林CV、CC、HLr隨坡位的變化鮮有報道。本研究發現,這些林冠結構指標在不同坡面也存在較大的變異性(分別為39.70%—49.06%、20.35%—24.15%、47.33%—57.00%),特別是上層林冠的短期波動更明顯。作者認為冰凍干擾是影響林冠結構變化最直接的驅動力,不同坡面受災程度不同[33,37],導致不同坡面林冠結構的差異。

4.3 災后林冠結構的年際動態

研究表明,干擾會改變森林冠層原有的結構,增加冠層結構的空間異質性[22,30]。干擾過后森林在恢復過程中冠層結構的空間異質性還在不斷地變化,并對森林的格局與生態學過程有決定作用[30]。區余端和蘇志堯對粵北山地常綠闊葉林3a的監測表明,災后森林恢復初期樣地林冠開度的總體面積在逐年縮小,森林郁閉度越來越高[30]。本研究發現,2008年冰凍災害后,隨著森林的恢復,大明山常綠闊葉林的CH、CV、CC、LAI均呈現不同程度的增加,災后第3年的增幅最大,尤以CV和CC增幅最明顯,分別比2009年提高50.65%和24.97%,之后各指標又有所回落。出現此種狀況可能是災害發生后的前兩年(2009和2010年),受2009年和2010年春西南持續干旱影響,常綠闊葉林的林冠結構恢復較慢；而2012年恢復速率出現下降,如CV、CC、LAI分別比2011年減小10.69%、4.33%和5.58%,這與2012年春大明山常綠闊葉林又遭受較強的冰凍干擾,部分林木受損有關。由于山地常綠闊葉林區域自然干擾時有發生,常綠闊葉林的林冠結構總是處于不斷變化之中,因此,在利用LAI研究很多關鍵生態過程,如林冠光合作用、呼吸作用、冠層對降水的截留及水分蒸散、冠層中太陽輻射的傳遞、冠層下土壤水分的蒸發等時,除了考慮LAI的空間異質性外,還應考慮不同年份的變化。因為在基于過程的生態系統模型模擬和定量遙感分析中,描述LAI動態的能力直接決定著植被與大氣界面物質、能量交換的可靠性和準確性[38-39]。

4.4 樹冠面積與林木胸徑的關系

研究表明,樹冠面積與林木胸徑存在顯著的正相關關系[40]。張麗楠等的研究認為華北落葉松人工林樹冠面積與林木胸徑關系以二次方的擬合效果最佳[41]。Hemery等的研究認為,當林木胸徑為20—50 cm時,樹冠直徑(K)與胸徑(d)為線性相關[40]。在本研究中,樹冠面積與胸徑表現為極顯著的線性關系,與Hemery等的結果一致[40]。主要原因是大明山常綠闊葉林的林木胸徑較小,90%以上的林木胸徑小于50 cm；其次,正如本文中所指出的,大明山常綠闊葉林的樹冠面積較小,樹冠常表現為圓柱形。作者認為,這是山地常綠闊葉林樹冠結構與山地環境條件長期適應的結果。

[1] Pan Y D, Birdsey R A, Phillips O L, Jackson R B. The structure, distribution, and biomass of the world′s forests. Annual Review of Ecology, Evolution, and Systematics, 2013, 44: 593- 622.

[2] Leiterer R, Furrer R, Schaepman M E, Morsdorf F. Forest canopy-structure characterization: A data-driven approach. Forest Ecology and Management, 2015, 358: 48- 61.

[3] Lowman M D, Rinker H B. Forest Canopies. 2nd ed. London: Elsevier Academic Press, 2004.

[4] Lowman M D, Schowalter T D. Plant science in forest canopies-the first 30 years of advances and challenges (1980- 2010). New Phytologist, 2012, 194(1): 12- 27.

[5] Franklin J F, Van Pelt R. Spatial aspects of structural complexity in old-growth forests. Journal of Forestry, 2004, 102(3): 23- 28.

[6] Hansen A J, Phillips L B, Dubayah R, Goetz S, Hofton M. Regional-scale application of lidar: Variation in forest canopy structure across the southeastern US. Forest Ecology and Management, 2014, 329: 214- 226.

[7] McElhinny C, Gibbons P, Brack C, Bauhus J. Forest and woodland stand structural complexity: Its definition and measurement. Forest Ecology and Management, 2005, 218(1/3): 1- 24.

[8] Hopkin M. Biodiversity and climate form focus of forest canopy plan. Nature, 2005, 436(7050): 452- 452.

[9] Stork N E. Editorial: dynamics and processes in the canopy of an Australian tropical rainforest. Austral Ecology, 2007, 32(1): 2- 3.

[10] Brodersen C, Pohl S, Lindenlaub M, Leibundgut C, Von Wilpert K. Influence of vegetation structure on isotope content of throughfall and soil water. Hydrological Processes, 2000, 14(8): 1439- 1448.

[11] Hobbie S E. Effects of plant species on nutrient cycling. Trends in Ecology & Evolution, 1992, 7(10): 336- 339.

[12] Frolking S, Palace M W, Clark D B, Chambers J Q, Shugart H H, Hurtt G C. Forest disturbance and recovery: A general review in the context of spaceborne remote sensing of impacts on aboveground biomass and canopy structure. Journal of Geophysical Research, 2009, 114(G2): G00E02.

[13] Asner G P, Powell G V N, Mascaro J, Knapp D E, Clark J K, Jacobson J, Kennedy-Bowdoin T, Balaji A, Paez-Acosta G, Victoria E, Secada L, Valqui M, Hughes R F. High-resolution forest carbon stocks and emissions in the Amazon. Proceedings of the National Academy of Sciences of the United States of America, 2010, 107(38): 16738- 16742.

[14] Whitehurst A S, Swatantran A, Blair J B, Hofton M A, Dubayah R. Characterization of canopy layering in forested ecosystems using full waveform LiDAR. Remote Sensing, 2013, 5: 2014- 2036.

[15] Gao T, Hedblom M, Emilsson T, Nielsen A B. The role of forest stand structure as biodiversity indicator. Forest Ecology and Management, 2014, 330: 82- 93.

[16] Graf R, Mathys L, Bollmann K. Habitat assessment for forest dwelling species using LiDAR remote sensing: Capercaillie in the Alps. Forest Ecology and Management, 2009, 257(1): 160- 167.

[17] Parker G G, Smith A P, Hogan K P. Access to the upper forest canopy with a large tower crane. BioScience, 1992, 42(9): 664- 670.

[18] 吳毅, 劉文耀, 宋亮, 陳曦, 盧華正, 李蘇, 石賢萌. 基于林冠塔吊的附生植物生態學研究進展. 植物生態學報, 2016, 40(5): 508- 522.

[19] 王云霓, 鄧秀秀, 王彥輝, 曹恭祥, 于澎濤, 熊偉, 徐麗宏. 六盤山南坡華北落葉松人工林冠層LAI的坡面尺度效應. 生態學報, 2016, 36(12): 3564- 3571.

[20] 劉志理, 戚玉嬌, 金光澤. 小興安嶺谷地云冷杉林葉面積指數的季節動態及空間格局. 林業科學, 2013, 49(8): 58- 64.

[21] 蘇宏新, 白帆, 李廣起. 3類典型溫帶山地森林的葉面積指數的季節動態: 多種監測方法比較. 植物生態學報, 2012, 36(3): 231- 242.

[22] 陳夏, 桑衛國. 暖溫帶地區3種森林群落葉面積指數和林冠開闊度的季節動態. 植物生態學報, 2007, 31(3): 431- 436.

[23] 李根柱, 王賀新, 朱教君. 遼東次生林區主要闊葉林型葉面積指數季節動態. 生態學雜志, 2008, 27(12): 2049- 2055.

[24] 呂瑜良, 劉世榮, 孫鵬森, 張國斌, 張瑞蒲. 川西亞高山暗針葉林葉面積指數的季節動態與空間變異特征. 林業科學, 2007, 43(8): 1- 7.

[25] 童鴻強, 王玉杰, 王彥輝, 于澎濤, 熊偉, 徐麗宏, 周楊. 六盤山疊疊溝華北落葉松人工林葉面積指數的時空變化特征. 林業科學研究, 2011, 24(1): 13- 20.

[26] 張佳寧, 袁金國, 張莎. 2002—2011年河北省植被LAI時空變化特征. 南京林業大學學報: 自然科學版, 2015, 39(1): 88- 92.

[27] 曼興興, 米湘成, 馬克平. 雪災對古田山常綠闊葉林群落結構的影響. 生物多樣性, 2011, 19(2): 197- 205.

[28] 蘇志堯, 劉剛, 區余端, 戴朝暉, 李鎮魁. 車八嶺山地常綠闊葉林冰災后林木受損的生態學評估. 植物生態學報, 2010, 34(2): 213- 222.

[29] 鞏合德, 楊國平, 張一平, 劉玉洪, 鄭征, 甘建民. 哀牢山4類植物群落葉面積指數比較. 東北林業大學學報, 2007, 35(3): 34- 36.

[30] 區余端, 蘇志堯. 粵北山地常綠闊葉林自然干擾后冠層結構空間異質性動態. 植物科學學報, 2012, 30(3): 223- 229.

[31] 溫遠光, 元昌安, 李信賢, 和太平, 賴家業, 黃棉. 大明山中山植被恢復過程植物物種多樣性的變化. 植物生態學報, 1998, 22(1): 33- 40.

[32] 溫遠光, 和太平, 譚偉福. 廣西熱帶和亞熱帶山地的植物多樣性及群落特征. 北京: 氣象出版社, 2004.

[33] 溫遠光, 李婉舒, 朱宏光, 周曉果, 葉鐸, 王磊. 特大冰凍干擾對大明山常綠闊葉林樹冠及林冠層狀況的影響. 廣西科學, 2014, 21(5): 454- 462.

[34] 崔佳玉, 曾煥忱, 王永強, 張毅, 胡益珩, 蘇志堯. 銀瓶山自然保護區闊葉林冠層結構與輻射消減效應. 西北林學院學報, 2015, 30(4): 45- 49.

[35] 溫遠光. 連栽桉樹人工林植物多樣性與生態系統功能關系的長期實驗研究. 成都: 四川大學, 2006.

[36] 朱宏光, 李燕群, 溫遠光, 梁宏溫, 林建勇, 農友. 特大冰凍災害后大明山常綠闊葉林結構及物種多樣性動態. 生態學報, 2011, 13(19): 5571- 5577.

[37] 溫遠光, 林建勇, 朱宏光, 周曉果, 葉鐸, 王磊. 南亞熱帶山地常綠闊葉林林隙及其自然干擾特征研究. 廣西科學, 2014, 21(5): 447- 453.

[38] Asner G P, Scurlock J M O, Hicke J A. Global synthesis of leaf area index observations: Implications for ecological and remote sensing studies. Global Ecology and Biogeography, 2003, 12(3): 191- 205.

[39] Waring R H, Running S W. Forest Ecosystems: Analysis at Multiple Scales. 3rd ed. Burlington: Elsevier Academic Press, 2007.

[40] Hemery G E, Savill P S, Pryor S N. Applications of the crown diameter-stem diameter relationship for different species of broadleaved trees. Forest Ecology and Management, 2005, 215(1/3): 285- 294.

[41] 張麗楠, 王得祥, 郝亞中, 柴宗政, 張洋, 吳昊. 陜西寧東林業局華北落葉松人工林最適經營密度研究. 西北林學院學報, 2013, 28(1): 146- 150.

Short-term dynamics of canopy structure of evergreen broadleaved forest after a freezing disaster in 2008 in Damingshan, Southern China

ZHOU Xiaoguo1, WEN Yuanguang1,2,3,*, ZHU Hongguang1,2,3, WANG Lei1

1GuangxiKeyLaboratoryofForestEcologyandConservation,CollegeofForestry,GuangxiUniversity,Nanning530004,China2StateKeyLaboratoryforConservationandUtilizationofSubtropicalAgro-bioresources,Nanning530004,China3GangxiYouyiguangForestEcosystermResearchStation,Pingxiang532600,China

Canopy structure and dynamics are critical components in the functioning and key ecological processes of forest ecosystems. Evergreen broadleaved forests are the representative community type in the subtropics, but the characteristic of canopy structure of this forest is poorly understood. In the present study, 80 permanent plots (20 m × 20 m for each plot) on a typical slope (200 m × 160 m) in the mid-mountain region of Damingshan Mountain were built and equally divided into 5 groups on the slope (16 continuous plots on each slope segment). In order to reveal the feature and short-term dynamic regulation of the canopy structure in this evergreen broadleaved forest, we investigated and analyzed the slope effects and annual dynamics of canopy structure indexes, including canopy height (CH), canopy volume (CV), canopy cover (CC), ratio of high to low cover (HLr) and leaf area index (LAI) during the growing seasons from 2009 to 2012. Mean CH, CV, CC, HLr, and LAI of this evergreen broadleaved forest averaged (12.09±0.05) m, (2642.51±278.33) m3(in each plot), (59.90±3.29)%, 2.48±0.23, and 2.00±0.06, respectively. The canopy structure was multilayered with CC averaging 42.20% in the upper layer, 30.35% in the middle layer, and 18.05% in the lower layers. Meanwhile, the canopy structure showed differences between slopes and growth years. For different slopes, the coefficient of variation of the index ranked as follows: CV (29.84%—55.89%) > HLr (32.90%—53.52%) > LAI (22.48%—43.89%) > CC (16.61%—25.74%) > CH (8.26%—12.77%). For different growth years, the coefficient of variation of the index ranked as HLr (47.33%—57.00%) > CV (39.70%—49.06%) > LAI (21.58%—48.13%) > CC (20.35%—24.15%) > CH (9.19%—12.59%), showing that CH had relatively strong stability. We found significant effects of slope scale on LAI with an elevation of 0.34 when sliding down 100 m in the downslope direction. The results of two-way ANOVAs showed that slope position had significant and extremely significant effects on CH (P=0.022) and HLr (P<0.001), respectively, while year of growth had significant effects on HLr (P=0.013) and extremely significant effects on CV and CC (bothP<0.001). The interaction effects of slope × year on CV and LAI were significant (P=0.016 andP=0.017, respectively). Correlation analysis showed that there was extremely significant positive correlation between DBH and canopy area. Our results indicated that the canopy structure of evergreen broadleaved forests on Damingshan Mountain have the characteristics of relatively lower CH, CC, HLr, and LAI, and smaller CV. This could be attributed to the relatively higher elevation (934—1223 m), shallow soil (30—45 cm), and the frequency of freezing disturbances, especially the severe ice storm damage in 2008 on Damingshan Mountain. This could also be the result of the long-term adaptation of canopy structure in montane evergreen broadleaved forest to environmental conditions.

evergreen broadleaved forest; leaf area index; canopy cover; canopy volume; ratio of high to low cover; dynamics

國家自然科學基金項目(31460121,39460022, 30860059)；廣西研究生教育創新計劃項目(YCBZ2014008)

2016- 06- 22;

2016- 11- 04

10.5846/stxb201606221221

*通訊作者Corresponding author.E-mail: wenyg@263.net

周曉果, 溫遠光, 朱宏光, 王磊.2008特大冰凍災害后大明山常綠闊葉林林冠結構動態.生態學報,2017,37(4):1137- 1146.

Zhou X G, Wen Y G, Zhu H G, Wang L.Short-term dynamics of canopy structure of evergreen broadleaved forest after a freezing disaster in 2008 in Damingshan, Southern China.Acta Ecologica Sinica,2017,37(4):1137- 1146.