內(nèi)蒙大興安嶺林區(qū)年降水量變化特征及周期分析

邊玉明, 代海燕, 王 冰, 張 健, 巴 特, 張秋良

(1.內(nèi)蒙古農(nóng)業(yè)大學(xué) 林學(xué)院, 呼和浩特 010010 ; 2.內(nèi)蒙古生態(tài)與農(nóng)業(yè)氣象中心,呼和浩特 010051; 3.興安盟林業(yè)科學(xué)研究所, 內(nèi)蒙古 興安盟 137400)

內(nèi)蒙大興安嶺林區(qū)年降水量變化特征及周期分析

邊玉明1, 代海燕2, 王 冰1, 張 健3, 巴 特1, 張秋良1

(1.內(nèi)蒙古農(nóng)業(yè)大學(xué) 林學(xué)院, 呼和浩特 010010 ; 2.內(nèi)蒙古生態(tài)與農(nóng)業(yè)氣象中心,呼和浩特 010051; 3.興安盟林業(yè)科學(xué)研究所, 內(nèi)蒙古 興安盟 137400)

利用內(nèi)蒙古大興安嶺林區(qū)11個氣象站1971—2010年降水觀測資料，采用一元線性回歸、累積距平法、曼—肯德爾檢驗(Mann-Kendall)法和Morlet小波分析方法對內(nèi)蒙古大興安嶺林區(qū)近40年年降水量變化進(jìn)行了分析，以揭示該地區(qū)在氣候變暖大背景下年降水量的變化特征。結(jié)果表明：1971—2010年內(nèi)蒙古大興安嶺林區(qū)年降水量整體在波動中呈微弱下降的趨勢，傾向率為-7.4 mm/10 a，降水量年際變化幅度較大，降水氣候傾向率具有很強(qiáng)的區(qū)域特點，整體經(jīng)歷了豐—枯—豐—枯的變化過程；結(jié)合累積距平分析與曼—肯德爾檢驗(Mann-Kendall)法兩種方法判斷出降水量發(fā)生突變的大致時間為1975—1979年與1999年；Morlet小波分析結(jié)果表明，近40年內(nèi)蒙古大興安嶺林區(qū)年降水量包含多個不同時間尺度的周期，其中27 a左右的周期振蕩最為強(qiáng)烈，對降水作用最為顯著，且在2014—2020年可能出現(xiàn)降水偏多期。

內(nèi)蒙古大興安嶺; 林區(qū); 降水; Morlet小波分析

根據(jù)政府間氣候變化專門委員會(Intergovernmental Panel on Climate Change，IPCC)第一工作組第五次評估報告，1951—2012年全球平均地表溫度的升溫速率是1880年以來升溫速率的兩倍[1]；中國陸地表面平均溫度上升趨勢與全球基本一致[2]。氣候變暖帶來的氣溫、降水等氣象因子的變化直接或間接影響了植物生長和生態(tài)平衡，進(jìn)而引起一系列生態(tài)環(huán)境問題如河流斷流、濕地萎縮、森林草場退化等。相對溫度而言降水的區(qū)域性差異更明顯[3]，降水量的多寡與干旱的發(fā)生也有直接而密切的聯(lián)系，可以說是干旱能否發(fā)生的關(guān)鍵性因素之一。

降水是描述一個地區(qū)氣候變化的重要指標(biāo)之一。在氣候變暖大背景下有很多專家學(xué)者做了降水事件的研究。張代青等[4]利用Morlet小波變換，對中國大陸近54 a降水事件進(jìn)行多時間尺度周期性分析，其結(jié)果表明：中國大陸年降水量變化呈現(xiàn)多尺度性、主周期性、同步性和總趨勢減少的規(guī)律；李邦東等[5]對東北地區(qū)降水事件進(jìn)行了分析，得出東北地區(qū)年降水量總體有減少趨勢，但線性趨勢不顯著；董旭光等[6]分析了山東省近50 a來的降水日數(shù)和降水強(qiáng)度的氣候變化特征，認(rèn)為近50 a山東省年降水日數(shù)和強(qiáng)度存在著明顯的年代際振蕩，降水日數(shù)減少趨勢達(dá)到極顯著水平；曲學(xué)斌等[7]研究結(jié)果表明：呼倫貝爾全市年降水量存著在10 a，26 a間尺度上的周期性變化，整體未發(fā)現(xiàn)減少趨勢，與整個東北地區(qū)表現(xiàn)不一致?？偠灾?，過去鮮有對典型的區(qū)域氣候進(jìn)行研究，從而忽略了區(qū)域氣候條件和區(qū)域差異的影響。內(nèi)蒙古大興安嶺林區(qū)具有典型的區(qū)域氣候條件，是亞歐大陸北方森林帶的重要組成部分，大興安嶺山脈貫穿整個林區(qū)，其森林生態(tài)系統(tǒng)在涵養(yǎng)水源、保持水土、保護(hù)生物多樣性等方面發(fā)揮著不可替代的重要作用，被譽(yù)為“祖國北方的重要生態(tài)屏障”。同時，也是我國“氣候脆弱區(qū)”和氣候變化響應(yīng)最為敏感區(qū)之一[8]，其一旦破壞，不可復(fù)制。本文利用內(nèi)蒙古大興安嶺林區(qū)11個氣象站點1971—2010年降水的連續(xù)觀測數(shù)據(jù)對年降水量進(jìn)行趨勢分析、突變檢驗和降水時間序列變化特征多時間尺度周期性分析，力圖揭示在氣候變暖大背景下內(nèi)蒙古大興安嶺林區(qū)降水變化規(guī)律，這對林區(qū)的生態(tài)建設(shè)和水資源保護(hù)及配置具有極其重要的現(xiàn)實意義和應(yīng)用價值；也為該地區(qū)水資源脆弱性評價提供科學(xué)的參考價值。

1 資料與方法

1.1 研究區(qū)概況

內(nèi)蒙古大興安嶺林區(qū)位于我國內(nèi)蒙古自治區(qū)東北部，地理坐標(biāo)119°36′20″—125°20′50″E,46°08′40″—53°20′00″N，地處我國高緯地帶，東部與黑龍江省接壤，南至洮兒河，西部與呼倫貝爾大草原相連，西北部與俄羅斯毗鄰，地跨呼倫貝爾市、興安盟等9個旗市，是我國最大的集中連片的國有林區(qū)。林區(qū)屬寒溫帶大陸性季風(fēng)氣候，冬季漫長而寒冷，夏季短暫而炎熱；年降水量300～500 mm，集中于7—8月份；年平均氣溫為-3.2℃，極端氣溫達(dá)-49.8℃。

1.2 資料選取及處理

降水資料為內(nèi)蒙古自治區(qū)11個(呼倫貝爾市9個，興安盟2個)站點逐月降水觀測數(shù)據(jù)，數(shù)據(jù)資料由內(nèi)蒙古自治區(qū)氣象局提供。分析時間定在1971—2010年，這一時間內(nèi)降水?dāng)?shù)據(jù)完整，沒有缺失值。借用VisualFoxpro 6.0軟件建立林區(qū)降水因子的時間序列資料，并進(jìn)行時間一致性檢驗。

1.3 研究方法

本文主要對內(nèi)蒙古大興安嶺林區(qū)降水事件年際變化進(jìn)行趨勢分析、突變檢驗和多時間尺度的復(fù)雜結(jié)構(gòu)分析研究。方法介紹如下：氣候要素的趨勢變化一般采用一元線性回歸分析，即：Y(t)=kx+b，式中：Y(t)為氣候要素，t為時間，k為線性回歸系數(shù)；k為正值，表示氣候要素呈線性增加的趨勢，反之呈線性減少的趨勢。這里采用10k表示降水要素每10 a的變化率[9]。

累積距平法是一種比較常用的通過曲線來直觀判斷變化趨勢的一種分析方法。根據(jù)N個時刻的累積距平值繪制累積距平曲線對時間序列的變化趨勢進(jìn)行分析，曲線由上升變?yōu)橄陆祷蛴上陆底優(yōu)樯仙礊檗D(zhuǎn)折點，轉(zhuǎn)折點對應(yīng)的時間為突變大致時間[9]。

曼—肯德爾檢驗(Mann-Kendall)法是一種非參數(shù)統(tǒng)計檢驗方法，其樣本可以服從任何分布，也不受其他異常值的干擾，比較適合順序類型的檢驗，同時可以明確指出突變的大致開始時間，目前已廣泛應(yīng)用于降水、氣溫等方面研究[10-13]。

小波分析最初由法國Morlet首次提出的，本文在Matlab 7.0軟件平臺的支持下，采用Morlet小波計算小波系數(shù)，利用Surfer 12.0等繪圖軟件繪制小波變換系數(shù)實部等值線、小波變換系數(shù)模及小波方差圖。目前，小波分析已經(jīng)在氣象領(lǐng)域廣泛應(yīng)用，已有不少學(xué)者利用小波分析方法研究了不同地區(qū)降水量多時間尺度的變化特征[14-16]。

2 結(jié)果與分析

2.1 年降水量變化趨勢分析

從時間尺度來看，在1971—2010年內(nèi)蒙古大興安嶺林區(qū)年降水量整體在波動中呈微弱下降的趨勢，變化傾向率為-7.4 mm/10 a，沒有通過顯著性水平檢驗(圖1)，這與賈建英等[17]研究結(jié)果一致。其多年均值為452 mm，最大值出現(xiàn)在1998年為740 mm，最小值出現(xiàn)在2007年為300 mm，兩者相差達(dá)440 mm，說明年降水量年際變化幅度較大。從空間尺度來看，降水傾向率具有很強(qiáng)的區(qū)域性(圖2)。其中趨勢下降的站點全部位于西部和中部林區(qū)，而增加的站點全部位于東部，東部林區(qū)毗鄰東北平原，是否受東北平原季風(fēng)和太平洋季風(fēng)影響其原因有待進(jìn)一步分析。

圖11971-2010年內(nèi)蒙古大興安嶺林區(qū)

年降水量線性變化趨勢

圖21971-2010年內(nèi)蒙古大興安嶺林區(qū)降水氣候傾向率分布

2.2 年降水量突變檢驗與累積距平

采用曼—肯德爾檢驗(Mann-Kendall)突變檢驗法對內(nèi)蒙古大興安嶺林區(qū)1971—2010年年降水量進(jìn)行突變檢驗分析(圖3)。由圖3可以看出UFk曲線整體上經(jīng)歷了上升—下降—上升—下降的趨勢，可判斷林區(qū)降水經(jīng)歷了豐—枯—豐—枯的過程。在1977—2003年一直處于上升趨勢，持續(xù)時間比較長，并且在1983—1986年與1989—1995年上升趨勢通過了0.05的顯著性水平檢驗，但均未通過0.001的顯著性水平檢驗。UFk與UBk曲線在置信區(qū)間以內(nèi)大致有4個交點，其中在2009—2010年出現(xiàn)兩個交點，距離比較近，可近似看做為一個交點，由此推斷在1971—2010年年降水量可能發(fā)生了3次突變，大致時間分別為1975年、1999年、2009年。累積距平曲線分析表明年降水量大體經(jīng)歷了的下降—上升—下降三個階段(圖3)，1971—1979年降水量呈緩慢下降趨勢，負(fù)距平占89%，表明降水量低于歷史平均水平，為少雨期，在1979年達(dá)到歷史最低水平，1980—1999年呈現(xiàn)明顯的上升趨勢，正距平占70%，為多雨期，且持續(xù)時間較長，且1999年為明顯的一次多雨波動，距平值為300 mm，達(dá)到歷史最高水平，2000—2010年降水量呈下降趨勢，負(fù)距平占64%，其中2000—2007年下降趨勢明顯，2008—2010年下降趨勢非常緩慢，由此判斷降水量發(fā)生突變的年份大致為1979年與1999年。綜合以上兩種方法的結(jié)果判斷出內(nèi)蒙古大興安嶺林區(qū)年降水量發(fā)生突變的大致時間為1975—1979年與1999年。

2.3 年降水量周期性變化分析

實部等值線圖中虛線表示小波系數(shù)實部小于0，即降水量偏低(枯)；實線表示小波系數(shù)實部大于0，即降水量偏高(豐)。0值所對應(yīng)的時刻即枯—豐轉(zhuǎn)折點。圖4A可以看出近40 a內(nèi)蒙古大興安嶺林區(qū)年降水量包含了多個不同尺度的周期性變化，形成各種尺度正負(fù)相間的振蕩中心，存在明顯的年代和年際變化。年際變化中具有5個明顯的4 a，6 a，10 a，17 a，27 a尺度周期的嵌套結(jié)構(gòu)，其中27 a尺度主要經(jīng)歷了5個枯—豐交替，1977年、1985年、1998年、2001年為枯—豐交替的轉(zhuǎn)折點，并且具有全域性；17 a尺度主要經(jīng)歷了7個明顯枯—豐交替，枯—豐轉(zhuǎn)折點為1976年、1981年、1988年、1995年、2001年、2006年，且具有全域性；10 a，6 a，4 a分別經(jīng)歷了12個、9個、14個枯—豐交替，但是4 a與6 a尺度上的周期變化不具有全域性，在Matlab軟件中處理發(fā)現(xiàn)4 a尺度振蕩主要集中在1986—2010年，6 a尺度振蕩主要集中在1977—1998年。這表明林區(qū)降水變化受大、中、小多重時間尺度周期的共同作用。

年降水量小波方差圖(圖4B)顯示小波方差共出現(xiàn)5個峰值，對應(yīng)時間尺度從小到大分別為4 a，6 a，17 a，10 a，27 a，表明這5個周期在近40 a中對降水量起主要作用。其中最大峰值對應(yīng)著27 a的時間尺度，可視為第一主周期，說明對降水影響最為顯著，判讀其等值線閉合情況可知2013年左右為枯—豐交替點，在2013年以后的大約8年間(2014—2020年)可能會出現(xiàn)降水偏多期；10 a時間尺度對應(yīng)第二峰值，可視為第二主周期，判讀其等值線閉合情況可知在2010年以后的大約8 a間(2011—2018年)也可能會出現(xiàn)降水偏多期，17 a，6 a，4 a的時間尺度分別為第三、四、五主周期；小波系數(shù)的模值是反映不同時間尺度變化周期所對應(yīng)的能量密度在時間域的分布情況，系數(shù)的模值越大表明其所對應(yīng)時間尺度的周期性規(guī)律越強(qiáng)。通過判讀年降水量Morlet小波系數(shù)模時頻分布圖(圖4C)可得最顯著性周期出現(xiàn)在25～32 a尺度上，其中27 a左右時間尺度上的周期振蕩最強(qiáng)烈，而且主要集中在1973年左右和2010年以后，7～11 a的時間尺度周期振蕩次之，其他時間尺度周期性變化較弱。

圖31971-2010年內(nèi)蒙古大興安嶺林區(qū)年降水量突變分析和累積距平曲線

圖4內(nèi)蒙古大興安嶺林區(qū)年降水量Morlet小波變換系數(shù)實部時頻分布、小波方差圖和小波變換系數(shù)模時頻分布

3 結(jié)論與討論

3.1 結(jié) 論

(1) 在1971—2010年內(nèi)蒙古大興安嶺林區(qū)年降水量整體在波動中呈下降趨勢，與整個中國的趨勢是一致的，變化傾向率為-7.4 mm/10 a，降水量年際變化幅度較大，但沒有顯著性變化趨勢；降水氣候傾向率具有很強(qiáng)的區(qū)域特點。

(2) 近40 a內(nèi)蒙古大興安嶺林區(qū)年降水經(jīng)歷了豐—枯—豐—枯的變化過程；年降水量發(fā)生突變的大致時間在1975—1979年與1999年。

(3) 內(nèi)蒙古大興安嶺林區(qū)40 a年降水量包含了多個不同時間尺度的周期性變化，其年際變化中存在著5個明顯周期變化規(guī)律：4 a，6 a，10 a，17 a，27 a，在全球氣候變暖背景下27 a左右的周期振蕩最為強(qiáng)烈，為第一主周期，對降水作用最為顯著，而且主要集中在1973年左右和2010年以后，在2014—2020年可能出現(xiàn)降水偏多期。

3.2 討 論

降水量的多少是影響一個地區(qū)水資源的主要因素之一，而水資源狀況影響著地區(qū)植被分布、植被種類、植被多樣性和生態(tài)環(huán)境的好壞。研究[18-19]表明大興安嶺林區(qū)是整個內(nèi)蒙古水分條件最好的地區(qū)，近年來隨著氣候變暖，降水量減少，林區(qū)有進(jìn)一步暖干化趨勢，林業(yè)干旱發(fā)生更加頻繁，林區(qū)內(nèi)濕地、草原等生態(tài)系統(tǒng)存在退化趨勢，大興安嶺林區(qū)生態(tài)穩(wěn)定性的維持風(fēng)險加劇，內(nèi)蒙古大興安嶺林區(qū)作為我國重要的生態(tài)屏障，一旦遭到破壞，不僅影響的是整個東北地區(qū)，甚至對中國北方生態(tài)環(huán)境也造成不利影響；此外，一些研究結(jié)果[20]表明，氣象和氣候條件嚴(yán)重制約著森林火災(zāi)的發(fā)生發(fā)展，尤其是降水量的多寡對林區(qū)火災(zāi)的發(fā)生、蔓延有著重要影響，于文穎等[21]研究表明降水量與過火次數(shù)呈顯著的正相關(guān)。從本文研究來看雖然在未來幾年的時間尺度上林區(qū)降水會增多，但林區(qū)降水量還會呈現(xiàn)減小的趨勢，未來還要逐漸提高火險等級預(yù)警，加大森林火災(zāi)防控力度；加強(qiáng)對林業(yè)干旱的預(yù)報預(yù)測，及時采取積極有效的措施應(yīng)對林業(yè)干旱帶來的不利影響。

在氣候變暖的大背景下，內(nèi)蒙古大興安嶺林區(qū)年降水量呈下降趨勢，與中國和整個東北降水量趨勢基本一致，有些研究結(jié)果表明整個東北地區(qū)年降水量下降趨勢達(dá)到了顯著性水平[17]，而林區(qū)年降水量下降卻沒達(dá)到顯著性水平，這可能受研究區(qū)特殊的區(qū)域環(huán)境影響或因為時間序列和所選站點不同而造成的差異，這有待我們進(jìn)一步探討研究。另外，關(guān)于林區(qū)各季降水特點，在這里沒有進(jìn)行分析，季節(jié)降水的分配也應(yīng)該引起我們的注意。還有生長季雨日，大中小雨日數(shù)、雨強(qiáng)和頻率的變化等對林業(yè)干旱的影響還存在很多的不確定，仍需要進(jìn)一步研究。同樣等量級的降水，其頻率和強(qiáng)度不同，對土壤水分的有效性也不同，如小于10 mm降水的日數(shù)較多，雨量比例高，對林地補(bǔ)充土壤水的作用效果較差，同樣大雨暴雨次數(shù)多，地表徑流多，只有中雨日數(shù)和量級比例高才會對林地的作用較大，年降雨量整體變化趨勢不明顯，但有效降水仍是該區(qū)域關(guān)注的重點。內(nèi)蒙古大興安嶺林區(qū)生態(tài)修復(fù)和保護(hù)，離不開主要的限制因子水分條件的研究，同時相關(guān)的蒸發(fā)、生長期雨量變化、雨熱系數(shù)等也是以后研究的重點問題。

[1] IPCC. Summary for Policymakers of Climate Change 2013:The Physical Science Basis. Contribution of Working Group1to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change[R]. Cambridge:Cambridge University Press, 2013.

[2] 丁一匯,任國玉,石廣玉.氣候變化國家評估報告(I):中國氣候變化的歷史和未來趨[J].氣候變化研究進(jìn)展,2006,2(1):3-8.

[3] Liang L, Li L, Liu Q. Precipitation variability in Northeast China from 1961 to 2008[J]. Journal of Hydrology, 2011,404(1):67-76.

[4] 張代青,梅亞東,楊娜,等.中國大陸近54年降水量變化規(guī)律的小波分析[J].武漢大學(xué)學(xué)報:工學(xué)版,2010,43(3):278-282.

[5] 李邦東,趙中軍,舒黎忠,等.1961—2010年東北地區(qū)降水事件時空均勻性研究[J].氣象與環(huán)境學(xué)報,2014,30(3):52-58.

[6] 董旭光,顧偉宗,孟祥新,等.山東省近50年來降水事件變化特征[J].地理學(xué)報,2014,69(5):661-671.

[7] 曲學(xué)斌,吳昊.呼倫貝爾市53年氣候變化特征分析[J].水土保持研究,2014,21(1):178-182.

[8] 姜曉艷,劉樹華,馬明敏,等.東北地區(qū)近百年降水時間序列變化規(guī)律的小波分析[J].地理研究,2009,28(2):354-362.

[9] 魏鳳英.現(xiàn)代氣候統(tǒng)計診斷與預(yù)測技術(shù)[M].2版.北京:氣象出版社,2007:175-181.

[10] 陳少勇,王勁松,石圓圓,等.中國東部季風(fēng)區(qū)1961—2006年年平均氣溫變化特征[J].資源科學(xué),2009,31(3):462-471.

[11] 劉兆飛,王翊晨,姚治君,等.太湖流域降水、氣溫與徑流變化趨勢及周期分析[J].自然資源學(xué)報,2011,26(9):1575-1584.

[12] 伍紅雨,杜堯東.1961—2008年華南地區(qū)極端氣溫的非對稱變化[J].生態(tài)學(xué)雜志,2010,29(11):2241-2248.

[13] 劉宇峰,孫虎,原志華.1960—2007年汾河流域氣溫年際和季節(jié)性變化特征分析[J].資源科學(xué),2011,33(3):489-495.

[14] 李淼,夏軍,陳社明,等.北京地區(qū)近300年降水變化的小波分析[J].自然資源學(xué)報,2011,26(6):1001-1011.

[15] 張利平,夏軍,胡志芳.華北地區(qū)降水多時間尺度演變特征[J].氣候變化研究進(jìn)展,2008,4(3):140-144.

[16] 胡永寧,王林和,張國盛,等.毛烏素沙地1969—2009年主要氣候因子時間序列小波分析[J].中國沙漠,2013,33(2):390-395.

[17] 賈建英,郭建平.東北地區(qū)近46年氣候變化特征分析[J].干旱區(qū)資源與環(huán)境,2011,25(10):109-115.

[18] 張存厚,王明玖,李興華,等.近30年來內(nèi)蒙古地區(qū)氣候干濕狀況時空分布特征[J].干旱區(qū)資源與環(huán)境,2011,28(5):70-75.

[19] 包云,李曉兵,黃玲梅,等.1961—2007年內(nèi)蒙古降水時空分布[J].干旱區(qū)地理,2011,34(1):52-61.

[20] 徐明超,馬文婷.干旱氣候因子與森林火災(zāi)[J].冰川凍土.2012,34(3):603-608.

[21] 于文穎,周廣勝,趙先麗,等.大興安嶺林區(qū)火災(zāi)特征及影響因子[J].氣象與環(huán)境學(xué)報,2009,25(4):1-5.

SpatiotemporalDistributionandPeriodicityofAnnualPrecipitationinDaxing′anlingForestRegionofInnerMongolia

BIAN Yuming1, DAI Haiyan2, WANG Bing1, ZHANG Jian3, BA Te1, ZHANG Qiuliang1

(1.ForestryCollegeofInnerMongoliaAgriculturalUniversity,Hohhot010010,China;2.InnerMongoliaEcologyandAgriculturalMeteorologyCentre,Hohhot010051,China;3.XinganMengForestryScienceResearchInstitute,HingganLeague,InnerMongolia137400,China)

Based on the observational data of precipitation in 11 meteorological stations from 1951 to 2010 in forest region of Inner Mongolia Daxing′anling, we used a linear regression, accumulative anomaly method, Mann Kendall test (Mann-Kendall) method and Morlet wavelet analysis method to analyze the annual precipitation changes in forest region of Inner Mongolia Daxing′anling in recent 40 years and to reveal the variation characteristics of annual precipitation under the climate warming background. The results showed as follows： precipitation in forest region of Inner Mongolia Daxing′anling showed a slight downward with fluctuation, the trend tendency rate was -7.4 mm/decade, inter annual variation precipitation was significant, precipitation of climate trend rate had very obvious regional features, which experienced the wet—dry—wet—dry process; the time of precipitation mutation occurred roughly from 1975 to 1979 and in 1999 based on the cumulative distance average analysis and Mann Kendall (Mann-Kendall) test method. Morlet wavelet analysis results showed that annual precipitation in recent 40 years in Daxing′anling forest region of Inner Mongolia contained a number of periods with different time scales, the 27-year periodic oscillation was the strongest, and had the most significant effect on precipitation, and high precipitation period will occur in the years 2014 to 2020.

Daxing′anling of Inner Mongolia; forest region; precipitation; Morlet wavelet analysis

2016-09-04

：2016-10-30

中國科學(xué)院野外站聯(lián)盟項目“東北地區(qū)生態(tài)變化評估”(KFJ-SW-YW026)

邊玉明(1992—),男,河北保定人,碩士研究生,研究方向:森林資源經(jīng)營與管理。E-mail:bianyuming0801@126.com

張秋良(1960—),男,內(nèi)蒙古托縣人,博士生導(dǎo)師,博士,研究方向:森林資源經(jīng)營管理與森林生態(tài)。E-mail:18686028468@163.com

P426.6

：A

：1005-3409(2017)03-0146-05