劉建泉+李進(jìn)軍+郝虎

摘要 在樣地調(diào)查的基礎(chǔ)上，利用青海云杉胸徑—材積的一元模型、青海云杉林蓄積—生物量模型、相關(guān)分析和主成分分析法，對(duì)甘肅祁連山青海云杉林的蓄積量、生物量、碳儲(chǔ)量以及影響青海云杉林碳儲(chǔ)量的主要因子進(jìn)行了研究。結(jié)果表明，青海云杉林單位面積的蓄積量、生物量和碳儲(chǔ)量平均值分別為223.445 4 m3/hm2、136.075 5 t/hm2和96.186 4 t/hm2；總生物量和碳儲(chǔ)量分別達(dá)到2.59×107 t和1.52×107 t，其中中齡林和近熟林的生物量和碳儲(chǔ)量分別占總量的42.581 8%和41.134 8%，而幼齡林僅占1.080 9%；影響森林植被碳儲(chǔ)量的主要因子為環(huán)境水熱狀況，溫度是影響祁連山森林植被碳儲(chǔ)量最關(guān)鍵的因子。

關(guān)鍵詞 青海云杉林；蓄積量；生物量；碳儲(chǔ)量；影響因素；祁連山

中圖分類號(hào) S718.5 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 1007-5739（2017）12-0140-04

Abstract Based on sample survey，using a meta model of DBH-volume of Picea crassifolia，stand volume-biomass model of Picea crassifolia forest，correlation analysis and principal component analysis，the stand volume，biomass，carbon storage and the main factors affecting the carbon storage of Picea crassifolia forest in Gansu Qilian Mountains were studied.The results showed that the stand volume，biomass and carbon storage of per unit area averages was 223.445 4 m3/hm2，136.075 5 t/hm2 and 96.186 4 t/hm2.The total biomass and total carbon storage reached 2.59×107 t and 1.52×107 t，biomass and carbon storage of mid-age forest and near-mature forest accounted for 42.581 8% and 41.134 8% of the total，respectively，and young forest accounted for only 1.080 9%.The main factors that affectd the forest vegetation carbon storage was hydrothermal condition of the environment，the temperature was the most critical factor affecting forest vegetation carbon storage in Qilian mountains.

Key words Picea crassifolia forest；stand volume；biomass；carbon storage；influence factor；Qilian Mountains

碳源（carbon source）是指向大氣中釋放碳的過程、活動(dòng)或機(jī)制。碳匯（carbon sink）是指通過植樹造林、森林管理、植被恢復(fù)等措施，利用植物光合作用吸收大氣中的CO2，并將其固定在植被或土壤中，從而減少溫室氣體在大氣中濃度的過程、活動(dòng)或機(jī)制[1]。當(dāng)生態(tài)系統(tǒng)固定碳量大于排放的碳量，該系統(tǒng)就稱為大氣CO2的匯，反之，則為碳源[2]。森林生態(tài)系統(tǒng)是陸地生態(tài)系統(tǒng)的主體，它貯存了陸地生態(tài)系統(tǒng)76%～98%的有機(jī)碳；森林是陸地上最大的碳匯，森林碳匯（forest carbon sinks）是指森林植物吸收大氣中CO2并將其固定在植被或土壤中，從而減少該氣體在大氣中的濃度的過程、活動(dòng)或機(jī)制[3]。森林植被所固定的碳量約占陸地植被總固碳量的82.5%，是森林固碳能力的重要標(biāo)志和植被碳匯的主要載體[3-4]。中國森林碳匯占總碳匯的52.85%[1]，其在應(yīng)對(duì)氣候變化和增加大氣碳吸收中具有重要作用。

青海云杉（Picea crassifolia）林主要分布于祁連山和賀蘭山兩大山系，其中祁連山分布面積占總面積的94.6%[3]，是祁連山森林生態(tài)系統(tǒng)的主要組成成分，對(duì)發(fā)揮森林生態(tài)系統(tǒng)服務(wù)功能、維護(hù)山地—綠洲—荒漠生態(tài)系統(tǒng)的平衡、促進(jìn)西部經(jīng)濟(jì)社會(huì)的穩(wěn)定發(fā)展都具有重要意義。目前，對(duì)青海云杉林的研究主要集中在森林生態(tài)、森林水文、森林土壤等領(lǐng)域[5-10]，在生物量方面也有研究[11-14]，對(duì)其碳儲(chǔ)量和碳匯功能的研究相對(duì)較少[3-4，15-16]。本文對(duì)祁連山青海云杉林不同林組的生物量和碳儲(chǔ)量進(jìn)行估算，從減緩氣候變化方面論證青海云杉林的碳匯功能；對(duì)影響青海云杉林碳儲(chǔ)量的因素進(jìn)行分析，探討青海云杉林對(duì)氣候變化的響應(yīng)，為青海云杉林的保護(hù)、恢復(fù)、發(fā)展和森林碳匯功能的經(jīng)營管理提供參考。

1 研究區(qū)域與研究方法

1.1 研究區(qū)概況

研究區(qū)域位于甘肅省祁連山北坡寺大隆林區(qū)，99°31′～100°53′ E，38°14′～38°44′ N，海拔2 600～4 200 m，屬大陸性高寒半濕潤(rùn)山地氣候，年平均氣溫0.6 ℃，1 月平均氣溫 -13.1 ℃，7 月平均氣溫12.1 ℃；年平均降水量437.2 mm，主要集中在5—9月；年蒸發(fā)量1 066.2 mm，年相對(duì)濕度60%，年日照時(shí)數(shù)1 892.6 h，林木生長(zhǎng)期90～120 d；土壤和植被隨山地地形和氣候的差異形成明顯的垂直帶譜，分布有寒漠土、高山和亞高山草甸土、灌叢草甸土、山地栗鈣土、森林灰褐土以及沼澤土等土壤，發(fā)育有高山流石灘植被、高山和亞高山灌叢、溫性和寒溫性針葉林、山地草原等地帶性植被和沼澤等非地帶性植被，主要優(yōu)勢(shì)種有青海云杉、祁連圓柏（Sabina przewalskii）、金露梅（Potentilla fruticosa）、鬼箭錦雞兒（Caragana jubata）、吉拉柳（Salix gilashanica）、高山繡線菊（Spiraea alpina）、珠芽蓼（Polygonum viviparum）、針茅（Stipa spp.）、苔草（Carex spp.）、嵩草（Kobresia spp.）、委陵菜（Poten-tilla spp.）以及苔蘚等，生物多樣性十分豐富。寺大隆林區(qū)的青海云杉林主要分布在海拔2 500～3 300 m的陰坡和半陰坡，多為中齡林和近熟林，面積13 968.3 hm2。

1.2 研究方法

1.2.1 樣地調(diào)查。利用衛(wèi)星影像和林相圖，在研究區(qū)域沿海拔梯度每間隔100 m設(shè)置1塊調(diào)查樣地，共設(shè)6塊樣地，記錄樣地所在位置的地理坐標(biāo)；根據(jù)地理坐標(biāo)，選擇受干擾小、接近自然狀態(tài)的地塊設(shè)置樣地，樣地面積為20 m×20 m；進(jìn)行環(huán)境因子和樣地調(diào)查，并對(duì)樹高≥1.5 m的林木測(cè)量胸徑和株數(shù)，記錄樹高<1.5 m的幼樹株數(shù)，調(diào)查樣地均為中齡林。

1.2.2 蓄積量計(jì)算。林木材積與胸徑有密切的關(guān)系，國內(nèi)外林學(xué)家相繼提出了幾十個(gè)不同類型的一元材積模型[17]，根據(jù)文獻(xiàn)[6]資料，利用Microsoft Excel 2007提供的數(shù)據(jù)分析工具，對(duì)青海云杉胸徑與材積的一元模型進(jìn)行擬合，篩選出最佳單株材積回歸模型：

1.2.5 影響青海云杉林碳儲(chǔ)量的主要因素分析。祁連山區(qū)氣候垂直分帶非常明顯，海拔每升高100 m，降水量增加18.6 mm，平均氣溫降低0.58 ℃，由此可以推算出樣地的平均溫度和平均降水量，將此組數(shù)據(jù)補(bǔ)充到表1中，并將表1中的坡向、坡位進(jìn)行數(shù)字化（西北坡=1，北坡=2；上坡位=1，中坡位=2，下坡位=3），構(gòu)成17×6的數(shù)據(jù)矩陣，用相關(guān)分析（correlation analysis）和主成分分析（principal components an-alysis，PCA）法分析影響青海云杉林蓄積量的主要因素。全部計(jì)算用SPSS 19軟件在計(jì)算機(jī)上完成。

2 結(jié)果與分析

2.1 樣地調(diào)查結(jié)果

根據(jù)調(diào)查結(jié)果，樣地的基本情況見表1。

2.2 青海云杉林蓄積量與生物量的垂直分布

祁連山區(qū)的溫度和降水量隨著海拔的變化而發(fā)生明顯的規(guī)律性變化，深刻影響青海云杉林的分布及其蓄積量與生物量。根據(jù)公式（2）（5）（10）的計(jì)算結(jié)果（表2）表明，隨著海拔升高、降水增加、氣溫降低，青海云杉林內(nèi)樹高小于1.5 m的幼樹呈減少趨勢(shì)，林分的天然更新能力降低；林分的蓄積量和生物量隨著海拔的升高和氣溫的降低呈減少趨勢(shì)，海拔2 770 m的蓄積量、生物量和碳密度都達(dá)到最大，分別為309.197 5 m3/hm2、177.931 1 t/hm2和122.940 6 t/hm2，海拔 2 770～3 250 m的青海云杉林單位面積的蓄積量、生物量和碳儲(chǔ)量（碳密度）平均值分別為223.445 4 m3/hm2、136.075 5 t/hm2和96.186 4 t/hm2。

用海拔作為自變量，與青海云杉林蓄積量、生物量和碳儲(chǔ)量進(jìn)行回歸分析，擬合得到最佳的回歸方程（表3）；回歸方程符合多項(xiàng)式方程，相關(guān)系數(shù)的平方（復(fù)測(cè)定系數(shù)R2）為0.798 9～0.802 8；回歸方程揭示了青海云杉林蓄積量、生物量和碳儲(chǔ)量隨著海拔升高而逐漸降低的趨勢(shì)。

2.3 祁連山林區(qū)青海云杉林碳儲(chǔ)量估算

祁連山林區(qū)是我國重要的水源涵養(yǎng)林區(qū)和生態(tài)功能區(qū)，對(duì)應(yīng)對(duì)氣候變化有著舉足輕重的作用；青海云杉林是祁連山森林生態(tài)系統(tǒng)的主體，主要分布于甘肅和青海二省，甘肅省是青海云杉林的分布中心，分布面積13.42萬hm2。利用式（5）和式（10）對(duì)甘肅祁連山林區(qū)青海云杉林的生物量和碳儲(chǔ)量進(jìn)行估算，結(jié)果（表4）表明，青海云杉林的生物量和碳儲(chǔ)量分別達(dá)到2.59×107 t和1.52×107 t。其中，中齡林和近熟林的生物量和碳儲(chǔ)量分別占總量的42.581 8%和41.134 8%，而幼齡林只占1.080 9%。

2.4 青海云杉林碳儲(chǔ)量與樣地因子的相關(guān)性

相關(guān)分析結(jié)果（表5）表明，在選取的17個(gè)因子中，僅有郁閉度（x7）、草本高度（x11）與林分碳儲(chǔ)量有不顯著的負(fù)相關(guān)，幼樹株數(shù)（x13）與碳儲(chǔ)量有不顯著的正相關(guān)，其余14個(gè)因子與青海云杉林碳儲(chǔ)量均成顯著或極顯著的正相關(guān)或負(fù)相關(guān)，其中極顯著相關(guān)因子12個(gè)、顯著相關(guān)因子2個(gè)。因此，這14個(gè)因子是研究林分碳儲(chǔ)量的重要因子。

2.5 影響青海云杉林碳儲(chǔ)量的主要因素

主成分分析結(jié)果（表6）表明，①特征值≥1的主成分有3個(gè)，因而可以提取3個(gè)主成分；②前3個(gè)主成分的累計(jì)貢獻(xiàn)率已經(jīng)達(dá)到96.012 7%，可以基本反映影響林分碳儲(chǔ)量的全部指標(biāo)信息；③第一主成分上載荷較高的有平均氣溫、平均降水量、海拔、草本蓋度、坡度、坡位，基本反映了林地的水熱狀況；第二主成分上載荷較高的有草本高度、苔蘚蓋度、幼樹株數(shù)、大樹株數(shù)、平均胸徑，基本反映了林分結(jié)構(gòu)特點(diǎn)；第三主成分則是對(duì)第一、第二主成分的補(bǔ)充。綜上分析，林分環(huán)境的水熱狀況是影響青海云杉林碳儲(chǔ)量的主要因子，而且海拔、平均氣溫、平均降水量的載荷量分別達(dá)到0.9802、-0.985 0和0.980 2，是影響青海云杉林碳儲(chǔ)量的關(guān)鍵因子，特別是平均氣溫的影響最大。

3 結(jié)論與討論

甘肅祁連山林區(qū)青海云杉林的總生物量為2.59×107 t，平均生物量為193.2055 t/hm2，接近王金葉等[16]203.08 t/hm2、常學(xué)向等[12]242.98 t/hm2和彭守璋等[3]134.3～250.1 t/hm2的研究結(jié)果，大于張 雷等[11]128.61 t/hm2的研究結(jié)果；總碳儲(chǔ)量為1.52×107 t，平均碳儲(chǔ)量（碳密度）為96.186 4 t/hm2，平均碳儲(chǔ)量（碳密度）高于寧夏賀蘭山保護(hù)區(qū)青海云杉林64.51 t/hm2 [15]和青海大通縣林區(qū)青海云杉林65.68 t/hm2 [4]，接近祁連山青海云杉林70.4～131.1 t/hm2 [3]的研究結(jié)果。甘肅祁連山林區(qū)青海云杉林中齡林和近熟林的生物量和碳儲(chǔ)量分別占總量的42.581 8%和41.134 8%，至少能夠維持3個(gè)林組80～120年生物量和碳儲(chǔ)量的持續(xù)增長(zhǎng)，80～120年后也可以保持一定時(shí)期的穩(wěn)定，但幼齡林僅占1.080 9%，后備資源匱乏。因此，從可持續(xù)發(fā)展和持續(xù)應(yīng)對(duì)氣候變化角度出發(fā)，必須持續(xù)擴(kuò)大造林和封山育林面積，加強(qiáng)幼齡林撫育管理，有效地培育森林資源，以保證碳匯功能的持續(xù)發(fā)揮。

林分碳儲(chǔ)量是林木通過光合作用將CO2轉(zhuǎn)換為有機(jī)物儲(chǔ)存在林分中的產(chǎn)物，是林分生物量非常重要的表現(xiàn)形式。對(duì)影響甘肅祁連山青海云杉林碳儲(chǔ)量關(guān)鍵因子的分析證明，在祁連山林區(qū)，坡位和坡向的改變會(huì)影響林地的水熱條件；特別是隨著海拔升高，林地的水熱條件逐步發(fā)生變化：水分條件會(huì)顯著改善，而平均氣溫明顯下降，造成林木光合速率降低，吸收和固定CO2的能力減弱。因此，溫度成為影響祁連山森林植被碳儲(chǔ)量最關(guān)鍵的因子。

林分生物量的測(cè)定是估算林分碳儲(chǔ)量的基礎(chǔ)。標(biāo)準(zhǔn)木分段切割法是常用的生物量測(cè)定方法，具有測(cè)定結(jié)果相對(duì)準(zhǔn)確、人力物力消耗較大的特點(diǎn)，而且天然林資源保護(hù)工程和林地管理相關(guān)法規(guī)的實(shí)施對(duì)林木采伐有了嚴(yán)格的規(guī)定。因此，研究者把目光逐漸轉(zhuǎn)移到利用青海云杉林生物量—蓄積量模型研究林分的生物量[3，4，17-18]，具有良好的效果。

4 參考文獻(xiàn)

[1] 魯豐先，張艷，秦耀辰，等.中國省級(jí)區(qū)域碳源匯空間格局研究[J].地理科學(xué)進(jìn)展，2013，32（12）：1751-1759.

LU F X，ZHANG Y，QIN Y C，et al.Spatial patterns of provincial carbon source and sink in China[J].Progress in Geography，2013，32（12）：1751-1759.

[2] 方精云，郭兆迪，樸世龍，等.1981—2000年中國陸地植被碳匯的估算[J].中國科學(xué)D輯：地球科學(xué)，2007，37（6）：804-812.

FANG J Y，GUO Z D，PIAO S L，et al.Estimation of carbon sink in China land from 1981 to 2000[J].Science in China Series D，2007，37（6）：804-812.

[3] 彭守璋，趙傳燕，鄭祥霖，等.祁連山青海云杉林生物量和碳儲(chǔ)量空間分布特征[J].應(yīng)用生態(tài)學(xué)報(bào)，2011，22（7）：1689-1694.

PENG S Z，ZHAO C Y，ZHENG X L，et al.Spatial distribution characte-ristics of the biomass and carbon storage of Qinghai spruce（Picea crassi-folia） forests in Qilian Mountains[J].Chinese Journal of Ap-plied Ecology，2011，22（7）：1689-1694.

[4] 郭倩倩，賀康寧，劉碩，等.青海省大通縣青海云杉林碳儲(chǔ)量初步估算[J].西北林學(xué)院學(xué)報(bào)，2011，26（6）：51-55.

GUO Q Q，HE K N，LIU S，et al.Preliminary estimate of vegetation carbon storage of Picea crassifolia in Datong Country of Qinghai Province[J].Journal of Northwest Forestry Un-iversity，2011，26（6）：51-55.

[5] 趙維俊，劉賢德，徐麗恒，等.祁連山青海云杉林動(dòng)態(tài)監(jiān)測(cè)樣地土壤pH和養(yǎng)分的空間異質(zhì)性[J].干旱區(qū)地理，2015，38（6）：1179-1189.

ZHAO W J，LIU X D，XU L H，et al.Spatial heterogeneity of soil pH and nutrients of a dyn-amical monitoring plot of Picea crassifolia forest in Qilian Mountains[J].Arid Land Geography，2015，38（6）：1179-1189.

[6] 劉興聰.青海云杉[M].蘭州：蘭州大學(xué)出版社，1992：389-391.

[7] 劉建泉.影響東大河林區(qū)青海云杉林物種多樣性的因素[J].草業(yè)科學(xué)，2015，32（1）：28-35.

LIU J Q.Influence factors of species diversity of Picea crassifolia forest in Dongdahe forest area[J].Pratacultural Science，2015，32（1）：28-35.

[8] 彭煥華，趙傳燕，許仲林，等.祁連山青海云杉林冠層持水能力[J].應(yīng)用生態(tài)學(xué)報(bào)，2011，22（9）：2233-2239.

PENG H H，ZHAO C Y，XU Z L，et al.Water storage capacity of Qinghai spruce（Picea crassifolia） forest canopy in Qilian Mountains[J].Chinese Journal of Applied Ecology，2011，22（9）：2233-2239.

[9] 齊鵬，劉賢德，趙維俊，等.祁連山中段青海云杉林土壤養(yǎng)分特征[J].山地學(xué)報(bào)，2015，33（5）：538-545.

QI P，LIU X D，ZHAO W J，et al.Soil nutrient cha-racteristics of Picea crassifolia forest in the middle segment of Qilian Mo-untains[J].Journal of Mountain Science，2015，33（5）：538-545.

[10] 趙維俊，劉賢德，張學(xué)龍，等.祁連山青海云杉（Picea crassifolia）林土壤有機(jī)碳與化學(xué)性質(zhì)的相互關(guān)系[J].冰川凍土，2014，36（6）：1565-1571.

ZHAO W J，LIU X D，ZHANG X L，et al.Relationship between soil organic carbon content and chemical properties of Picea crassifolia for-est in the Qilian Mountains[J].Journal of Glaciology and Geocryology，2014，36（6）：1565-1571.

[11] 張雷，于澎濤，王彥輝，等.祁連山北坡青海云杉中齡林生物量隨海拔的變化[J].林業(yè)科學(xué)，2015，51（8）：1-7.

ZHANG L，YU P T，WANG Y H，et al.Biomass change of middle aged forest of Qinghai spruce along an altitudinal gradient on the north slope of Qilian Mountains[J].Scientia Silvae Sinicae，2015，51（8）：1-7.

[12] 常學(xué)向，車克鈞，宋彩福.祁連山林區(qū)青海云杉林群落生物量的初步研究[J].西北林學(xué)院學(xué)報(bào)，1995，11（1）：19-23.

CHANG X X，CHE K J，SONG C F.Preliminary studies on biomass of Picea crassifolia forest community in Qilian Mountains[J].Journal of Northwest Forestry University，1995，11（1）：19-23.

[13] 敬文茂，劉賢德，趙維俊，等.祁連山典型林分生物量與凈生產(chǎn)力研究[J].甘肅農(nóng)業(yè)大學(xué)學(xué)報(bào)，2011，46（6）：81-85.

JING W M，LIU X D，ZHAO W J，et al.Study on biomass and net productivity of typical forest stand in the Qilian Mountains[J].Journal of Gansu Agricultural University，2011，46（6）：81-85.

[14] 王金葉，車克鈞，傅輝恩，等.祁連山水源涵養(yǎng)林生物量的研究[J].福建林學(xué)院學(xué)報(bào)，1998，18（4）：319-323.

WANG J Y，CHE K J，F(xiàn)U H E，et al.Study on Biomass of Water Conservation Forest on North Slope of Qilian Mou-ntains[J].Journal of Fujian College Forestry，1998，18（4）：319-323.

[15] 季波，王繼飛，何建龍，等.寧夏賀蘭山自然保護(hù)區(qū)青海云杉林的有機(jī)碳儲(chǔ)量[J].草業(yè)科學(xué)，2014，31（8）：1445-1449.

JI B，WANG J F，HE J L，et al.Carbon storage of Qinghai spruce forest in Ningxia Helan Nature Reserve[J].Pratacultural Science，2014，31（8）：1445-1449.

[16] 王金葉，車克鈞，蔣志榮.祁連山青海云杉林碳平衡研究[J].西北林學(xué)院學(xué)報(bào)，2000，15（1）：9-14.

WANG J Y，CHE K J，JIANG Z R.A study on carbon balance of Picea crassifolia in Qilian Mountains[J].Journal of Northwest Forestry Unive-rsity，2000，15（1）：9-14.

[17] 張連翔，林陽.林木胸徑與材積的關(guān)系：Logistic衍生模型[J].東北林業(yè)大學(xué)學(xué)報(bào)，2001，29（2）：99-101.

ZHANG L X，LIN Y.Relationship between diameter at breast height and volume of forest-tree： Logistic derivation model[J].Journal of Northeast Forestry University，2001，29（2）：99-101.

[18] 王斌，劉某承，張彪.基于森林資源清查資料的森林植被凈生產(chǎn)量及其動(dòng)態(tài)變化研究[J].林業(yè)資源管理，2009（1）：35-43.

WANG B，LIU M C，ZHANG B.Dynamics of Net Production of Chinese Forest Vegetation Based on Forest Inventory Data[J].Forest Resources Management，2009（1）：35-43.