中國東北地區(qū)近50年凈生態(tài)系統(tǒng)生產(chǎn)力的時空動態(tài)

李 潔，張遠東，顧峰雪，*，黃 玫，郭 瑞，郝衛(wèi)平，夏 旭

(1.中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)環(huán)境與可持續(xù)發(fā)展研究所，農(nóng)業(yè)部旱作節(jié)水農(nóng)業(yè)重點實驗室，北京 100081;2.中國林業(yè)科學(xué)研究院森林生態(tài)環(huán)境與保護研究所，國家林業(yè)局森林生態(tài)環(huán)境重點實驗室，北京 100091;3.中國科學(xué)院地理科學(xué)與資源研究所生態(tài)系統(tǒng)觀測與模擬重點實驗室，北京 100101)

陸地生態(tài)系統(tǒng)在全球碳循環(huán)中具有重要作用［1］，在一定程度上具有減緩大氣CO2濃度增加和氣候變化的作用［2］。由于京都議定書已將陸地生態(tài)系統(tǒng)碳吸收納入到溫室氣體減排的框架之中［3］，因此對陸地碳匯時空變化的研究對京都議定書的執(zhí)行具有重要意義［4-5］。氣候變化是影響陸地生態(tài)系統(tǒng)碳循環(huán)的因素之一［6-8］。在氣候變暖條件下，中高緯度地區(qū)生態(tài)系統(tǒng)NPP通常表現(xiàn)為增加，土壤呼吸作用增強［9］，這兩個過程都是由溫度升高所導(dǎo)致，對陸地碳匯卻具有相反的作用。因此，綜合分析生態(tài)系統(tǒng)碳循環(huán)過程對氣候變化的響應(yīng)，對于量化和理解陸地碳匯的區(qū)域分布和年際變化具有重要意義［5］。

作為凈初級生產(chǎn)力NPP和異養(yǎng)呼吸HR的差值，凈生態(tài)系統(tǒng)生產(chǎn)力(NEP)直接定性定量地描述陸地生態(tài)系統(tǒng)碳源/匯的性質(zhì)和能力。在不考慮各種自然和人為擾動的情況下，NEP可近似看作是陸地生態(tài)系統(tǒng)與大氣系統(tǒng)之間的凈碳交換量［10］，最近十幾年來，越來越受到科學(xué)家的關(guān)注［11-13］。以往研究表明，從20世紀(jì)80年代到90年代，陸地生態(tài)系統(tǒng)的碳蓄積量明顯增加了［2，14-15］，并且北半球中高緯度陸地生態(tài)系統(tǒng)是其中一個巨大碳匯［2，12，16］，但對陸地碳匯時空格局的變化及其與氣候變化關(guān)系的方面仍存在很大不確定性［17-18］。Cao等研究了氣候變化對中國陸地生態(tài)系統(tǒng)碳蓄積的影響，在氣候變化和CO2濃度增加的共同作用下，中國陸地生態(tài)系統(tǒng)總體上具有碳匯作用，并且碳源匯分布存在區(qū)域差異，中國大部分地區(qū) NEP接近于零［6］。Piao基于ORCHIDEE模型研究了氣候變化對青藏地區(qū)草地生態(tài)系統(tǒng)NPP、NEP的時空變化的影響，結(jié)果表明降水對NEP增加的貢獻為36%，是引起NEP年際變化的主要驅(qū)動因素，溫度對該區(qū)凈碳蓄積增加沒有明顯作用［13］。Lu等基于TEM模型在美國中西部地區(qū)的研究表明，美國中西部地區(qū)陸地生態(tài)系統(tǒng)NEP具有0.62TgC/a的增加趨勢，年平均NEP與降水呈正相關(guān)，與溫度的相關(guān)性不顯著［8］。以上研究結(jié)果表明，氣候變化對NEP的空間格局變化的影響存在區(qū)域差異，而且研究時段主要在1980年以后的20a，而陸地碳匯的變化是長期氣候變化的結(jié)果［2，5］，因此有必要在更長時間尺度上對NEP的時空動態(tài)及其對氣候變化的響應(yīng)進行研究。

東北地區(qū)處于我國最高緯度地區(qū)，是中國最大的天然林區(qū)和重要的商品糧生產(chǎn)基地，也是世界三大黑土帶之一［19-21］。東北地區(qū)是全球氣候變化的敏感區(qū)域之一［22］。鑒于東北地區(qū)特殊的地理位置以及在全球變化中的重要作用，國內(nèi)許多學(xué)者開展了對東北地區(qū)碳收支的研究［20，23-25］，但研究大多是對NPP和土壤呼吸這兩個碳收支過程的單獨研究，而且主要集中于森林和農(nóng)田生態(tài)系統(tǒng)［20，24-25］，不能反映整個東北地區(qū)的凈碳收支狀況。目前僅有的對東北地區(qū)NEP的研究，研究時間段短，并且沒有對該地區(qū)NEP的時空動態(tài)變化以及與氣候變化的關(guān)系進行深入探討。王軍邦等模擬了東北地區(qū)NPP和NEP的空間格局［26］。趙俊芳等人模擬了當(dāng)前氣候狀況和不同氣候情景對東北森林生態(tài)系統(tǒng)NPP和NEP的影響，認為NPP和NEP對溫度升高的反應(yīng)比對降水的反應(yīng)更為敏感［27］?；诖?，本研究應(yīng)用高分辨率的氣候數(shù)據(jù)庫和生態(tài)系統(tǒng)過程模型CEVSA，對1961—2010年東北地區(qū)凈生態(tài)系統(tǒng)生產(chǎn)力(NEP)的時空格局及變化趨勢進行分析，并討論氣候變化與NEP的關(guān)系，闡明影響東北區(qū)域NEP變化的關(guān)鍵氣候因子。

1 資料與方法

1.1 模型簡介

CEVSA建立于1998年，該模型是一個基于生理生態(tài)過程模擬植物-土壤-大氣系統(tǒng)能量交換和水碳氮耦合循環(huán)的生物地球化學(xué)循環(huán)模型［28-30］。CEVSA模型包含3個子模型:植被子模型用來計算植被的NPP、分配和凋落物產(chǎn)量;生物物理子模型計算蒸散、土壤水分動態(tài)和氣孔導(dǎo)度;生物地球化學(xué)子模型計算凋落物和土壤有機碳的分解和傳輸［28］。CEVSA模型已應(yīng)用于區(qū)域和全球尺度，模擬分析了陸地生態(tài)系統(tǒng)碳循環(huán)對氣候變化的響應(yīng)［5，28-36］。以下是對模型主要過程模擬方法的簡要描述，關(guān)于模型更詳細的結(jié)構(gòu)可以參考相關(guān)文獻［28，37-38］。

1.1.1 光合和氣孔導(dǎo)度

光合作用速率決定于葉肉組織光合酶對CO2的利用效率和CO2向葉肉組織的擴散速率，氣孔導(dǎo)度的變化采用改進的Ball-Berry模型模擬。

式中，Ab代表由酶系統(tǒng)活性決定的碳同化速率，gs為氣孔傳導(dǎo)度，Wc由Rubisco所決定的羧化速率，Wj是由電子傳遞速度決定的羧化速率，Wp決定于光合反應(yīng)過程對磷酸丙糖(Triose phosphate)的利用效率。Po和Pc分別是葉肉組織中 O2和 CO2的分壓。τ是Rubisco的特異因子，Rd為白晝非光合呼吸速率。go是在光補償點下光合速率為零時的氣孔傳導(dǎo)度，g1是靈敏度參數(shù)，A為實際光合速率，T為絕對溫度，kg(ws)為氣孔導(dǎo)度對土壤含水量ws的反應(yīng)函數(shù)。

1.1.2 氮吸收和自養(yǎng)呼吸

植物的氮含量決定了植物的最大光合速率，氮吸收對溫度的依賴性(NT)可表達為:

式中，Amax為最大光飽和的光合速率，N為某一層葉片的氮含量，NT則為植物的總氮吸收量，u1、u2、u3和kT(T)均為溫度的響應(yīng)函數(shù)，Tk為絕對溫度。葉片的氮含量(N)還決定了葉片的維持呼吸(Rml)［39］:

式中，r1(T)和r2(T)分別是溫度的響應(yīng)函數(shù)。

1.1.3 土壤碳氮動態(tài)

CEVSA模型將土壤有機質(zhì)分為8個碳庫，每個碳庫都具有特定的分解速率，并與碳庫本身特性、溫度、濕度、氮的有效性和土壤質(zhì)地有關(guān)。氮礦化量(Na)和土壤異養(yǎng)呼吸(HR)可分別用下式表達:

式中，i表示不同的碳庫。OM是每個碳庫的大小，K是潛在衰變率，CN是碳氮比，Kag(i)是實際各個土壤有機碳庫的衰變率，ε為同化效率。

1.2 輸入數(shù)據(jù)來源、參數(shù)獲取和模型運行

模型中所使用的氣象數(shù)據(jù)來自于國家氣象信息中心，包括1961—2010年全國756個氣象臺站的每旬平均氣溫，降水量，云量和相對濕度。使用ANUSPLIN4.1插值軟件的樣條函數(shù)插值法對氣象數(shù)據(jù)進行內(nèi)插，得到全國1961—2010年的0.1°氣象柵格數(shù)據(jù)。大氣 CO2濃度資料來源于美國夏威夷Mauna Loa觀測所(http://co2now.org/Current-CO2/CO2-Now/noaa-mauna-loa-co2-data.html)。土壤類型和質(zhì)地資料取自1∶100萬土壤類型圖和第二次土壤普查數(shù)據(jù)，對其進行數(shù)字化并重采樣到0.1°。植被數(shù)據(jù)來自于(European Commission，Joint Research Centre，2003. http://bioval.jrc.ec.europa.eu/products/glc2000/glc2000.php)，共有22種土地覆被類型，在ArcGIS中進行插值重采樣，匹配于氣象數(shù)據(jù)的分辨率。

CEVSA模型包含43個參數(shù)和11個初始的狀態(tài)變量。在CEVSA模型應(yīng)用于區(qū)域和全球的模擬時，初始的狀態(tài)變量，如植被碳、土壤碳、土壤碳氮比、初始的土壤含水量等來源于靜態(tài)模型的輸出。與土壤質(zhì)地有關(guān)的參數(shù)包含田間持水量、飽和含水量、萎蔫系數(shù)、土壤分解的最佳持水量等參數(shù)，這些土壤水分參數(shù)由土壤顆粒組成確定，依據(jù)Raich等［40］、McGuie等［41］和 Cao 和 Woodward［28］等文獻中提供的數(shù)據(jù)進行估算。CEVSA最初包含13種土地覆蓋類型，與植被類型相關(guān)的參數(shù)主要是根、莖、葉的周轉(zhuǎn)時間、根、莖、葉的分配比例和比葉面積等，主要來源于已發(fā)表文獻的估計。在本研究中使用的22種土地覆蓋類型，如與最初模型中已有土地覆蓋類型直接對應(yīng)的，就使用原有參數(shù)，不直接對應(yīng)的植被類型，本研究通過查閱文獻和模型調(diào)試的方法獲取相對應(yīng)植被類型的相關(guān)參數(shù)。詳細的參數(shù)列表和取值方法可以參考已發(fā)表的文章，如 Cao 和 Woodward［28］、Zhang 等［38］、Gu 等［42］。

本研究中，模型運行的時間分辨率為10d，空間分辨率為0.1°。首先應(yīng)用1961—2010年的平均氣候數(shù)據(jù)運行模型至生態(tài)系統(tǒng)平衡態(tài)，即各個狀態(tài)變量如植被、土壤碳貯量以及土壤含水量等年際變化量小于0.1%，且NPP、凋落物產(chǎn)生量(LT)與HR相等，然后用1961—2010年每旬資料進行動態(tài)模擬，并反復(fù)運行模型以消除假定的生態(tài)系統(tǒng)狀態(tài)變量的初始值(即平衡態(tài)假設(shè))對模擬結(jié)果的影響。從最終模擬結(jié)果中提取中國東北地區(qū)的數(shù)據(jù)進行分析。

1.3 數(shù)據(jù)分析

(1)時間序列變化分析:對年平均溫度、降水、NEP序列做一元線性回歸分析，線性方程的斜率代表溫度、降水、NEP的變化趨勢。(2)空間序列變化分析:逐個格點計算NEP一元線性回歸的回歸系數(shù)，得到NEP的空間變化趨勢;(3)由于本研究中的NEP的計算結(jié)果，是通過溫度與降水等因子在模型中運算得到的，NEP與溫度和降水之間不具有獨立性，但出于分析的需要，本研究不考慮這一關(guān)系，采用相關(guān)分析的方法，研究溫度、降水對NEP的影響程度。NEP與溫度降水的相關(guān)分析:用PEARSON相關(guān)系數(shù)法，計算NEP與年平均溫度、降水量的相關(guān)系數(shù)，并進行顯著性檢驗，當(dāng)相關(guān)系數(shù)通過顯著性水平(P＜0.05)時，認為NEP與溫度、降水有顯著相關(guān)關(guān)系。

1.4 研究區(qū)域概況

本研究中，東北地區(qū)(北緯 38°48'—53°33'N，東經(jīng) 115°31'—135°05'E)主要包括黑龍江、吉林、遼寧以及內(nèi)蒙古地區(qū)東北部，北部與俄羅斯接壤，東南部與朝鮮半島相接，南部濱臨中國渤海和黃海［43］。本研究中統(tǒng)計的總面積為11.73×105km2，1961—2010東北地區(qū)年平均溫度變化在2.49—6.02℃之間，從南到北隨著緯度升高，年均溫逐漸降低，年降水量變化在427—680mm，從東南向西北，降水量逐漸減少。冬季寒冷干燥，夏季炎熱多雨，特殊的地理位置和氣候特點，形成了東北地區(qū)獨特的植被分布格局。從南向北隨著熱量的變化，出現(xiàn)暖溫帶、溫帶和寒溫帶，植被也相應(yīng)形成了暖溫帶落葉闊葉林帶，溫帶針闊混交林帶和寒溫帶針葉林帶;從東南向西北隨著濕度的分異，出現(xiàn)濕潤地區(qū)、半濕潤地區(qū)和半干旱地區(qū)，植被則有森林、森林草原、典型草原植被［44］。

2 結(jié)果分析

2.1 模型驗證

CEVSA模型所應(yīng)用的生物和生態(tài)學(xué)原理、計算方程和參數(shù)均取自于大量的實驗室和野外試驗觀測，曾應(yīng)用葉片和植株生理試驗、樣點初級生產(chǎn)力觀測和遙感反演數(shù)據(jù)進行了驗證，模型估算的NPP、葉面積指數(shù)及植被和土壤中的碳儲量與實地調(diào)查和測定所得的結(jié)果有很好的一致性［28，30］。但 CEVSA模型自建立后一直缺乏在生態(tài)系統(tǒng)尺度上的驗證。20世紀(jì)90年代后大型環(huán)境控制試驗和渦度相關(guān)技術(shù)的應(yīng)用提供了連續(xù)、長期的生態(tài)系統(tǒng)尺度碳通量觀測數(shù)據(jù)［45-47］。它們?yōu)檠芯可鷳B(tài)系統(tǒng)尺度的碳交換［48-49］以及模型的檢驗測試［50］提供了大量有用的數(shù)據(jù)。為檢驗該CEVSA模型對生態(tài)系統(tǒng)碳通量季節(jié)動態(tài)的模擬效果，評價模型對生態(tài)系統(tǒng)碳循環(huán)過程機理定量描述的準(zhǔn)確性，顧峰雪等［51］在2007年曾應(yīng)用多個站點的渦度相關(guān)觀測數(shù)據(jù)對模型的模擬結(jié)果進行了驗證和評價，包括位于東北地區(qū)的長白山站，結(jié)果表明CEVSA模型能夠較好地模擬長白山針闊混交林碳交換的季節(jié)和年際動態(tài)(圖1)。

圖1 長白山站NEP模擬值與通量觀測值的比較NEP觀測值與模擬值的季節(jié)動態(tài);NEP觀測值與模擬值的1∶1線Fig.1 Comparison between modeled and measured monthly net ecosystem productivity(NEP)in Changbaishan station from 2003 to 2005 simulated and observed monthly NEP;correlation of simulated and observed monthly NEP

另外，為了進一步驗證模型在區(qū)域模擬時對碳收支的模擬效果，本研究應(yīng)用東北地區(qū)81個實測站點地面NPP數(shù)據(jù)［52-54］對模型模擬的NPP進行了驗證，這81個實測站點包含了落葉闊葉林、落葉松林、針闊混交林等東北地區(qū)的多種森林類型。原數(shù)據(jù)中以生物量表示的NPP，本文中取碳轉(zhuǎn)化率0.45，得到以碳含量表示的實測NPP。81個實測站點NPP數(shù)據(jù)的平均值為491 gC·m-2·a-1，根據(jù)樣點坐標(biāo)提取的1961—2010年模擬NPP年均值為479 gC·m-2·a-1，平均相對誤差為2.2%。NPP年均模擬值與對應(yīng)樣點NPP實測值間的相關(guān)系數(shù)為r=0.584(P＜0.001)，NPP實測值與模擬值基本吻合(圖2)。同時，本研究將CEVSA模型模擬的NPP與已有研究進行了比較，結(jié)果表明，本研究模擬的東北地區(qū)平均的NPP為439 gC·m-2·a-1，在以往其他研究結(jié)果的范圍內(nèi)(400.85—445 gC·m-2·a-1)［23，55-56］。另外，本研究收集了東北地區(qū)17個實測土壤異養(yǎng)呼吸數(shù)據(jù)［57-66］，這17個實測數(shù)據(jù)包含了農(nóng)田、森林、草地等多種植被類型，模擬值與對應(yīng)經(jīng)緯度土壤異養(yǎng)呼吸實測值間存在顯著相關(guān)性，相關(guān)系數(shù)為r=0.880(P＜0.01)(圖2)。

圖2 NPP，HR模擬值與地面實測值的比較Fig.2 Comparison of modeled NPP and HR with the observed values in the Northeast

2.2 NEP的年際及年代際變化

近50年中，東北地區(qū)陸地生態(tài)系統(tǒng)是大氣CO2的匯，但碳吸收的年際間變異很大。1961—2010年，東北地區(qū)年NEP總量最低值為-0.094PgC/a，最高值為0.117PgC/a，年平均為0.026PgC/a。年NEP單位平均值變化范圍為-79.0—99.8 gC·m-2·a-1，平均為22.5 gC·m-2·a-1，變異系數(shù)為206%。研究時段內(nèi)年NEP總量具有下降趨勢，但線性變化趨勢不顯著(P=0.240)。在年代際總量上，20世紀(jì)80年代的碳匯強度最大(52.47TgC/a)，占研究時段內(nèi)總碳吸收量的40%。60年代和70年代的碳吸收處于較高水平(27.73TgC/a，26.11TgC/a)，分別占研究時段總量的21%和20%。90年代和2000—2010年后碳吸收量急劇減少(13.02TgC/a和10.51TgC/a)，分別占研究時段總量的10%和8%。研究時段內(nèi)，僅有30%的年份表現(xiàn)為碳源，且碳源強度較弱，除20世紀(jì)90年代后期以及21世紀(jì)初期向大氣釋放了少量碳外，中國東北地區(qū)陸地生態(tài)系統(tǒng)在其他時期主要表現(xiàn)為碳匯，尤其是80年代(圖3)。

圖3 研究時段內(nèi)東北地區(qū)凈生態(tài)系統(tǒng)生產(chǎn)力NEP的年際變化及年代際變化Fig.3 Interannual variations and decadal variations in NEP during 1961—2010 in the Northeast

研究時段內(nèi)，東北地區(qū)的年平均溫度變化在2.49—6.02℃之間，溫度升高明顯(P＜0.05)，年降水量427.26—680.72mm，趨勢傾向率為-5.4mm/10a，變化趨勢不顯著，降水的年際變化較大，變異系數(shù)為9.62(圖4)。NEP的年際變化與溫度呈顯著負相關(guān)(r=-0.343，P＜0.05)，與降水相關(guān)性極顯著(r=0.859，P＜0.01)，同時東北地區(qū) NEP和年降水量的變化規(guī)律基本一致，即同期上升或達到最高值。因此，溫度和降水共同作用導(dǎo)致東北地區(qū)NEP的年際變化，而年降水量的變化對 NEP年際變化起主要作用。

圖4 1961—2010年東北地區(qū)年平均溫度、降水量的年際變化及年代際變化Fig.4 Interannual variations and decadal variations in mean annual temperature and mean annual precipitation during 1961—2010 in the Northeast

2.3 NEP的空間分布及變化趨勢

東北地區(qū)NEP呈現(xiàn)出東部高，西部和中部低，北部高，南部低的空間分布特點。NEP正值代表碳匯，NEP負值代表碳源。東北平原的部分地區(qū)，內(nèi)蒙古高原部分地區(qū)以及三江平原小部分地區(qū)為碳源區(qū)，碳源區(qū)僅占全區(qū)總面積的15.7%。碳匯值較高的地區(qū)分布在大、小興安嶺和長白山的森林分布區(qū)，約為45—60gC·m-2·a-1，部分地區(qū)達到60—100 gC·m-2·a-1。碳匯較低的地區(qū)主要分布在遼寧，三江平原，內(nèi)蒙古與遼寧、吉林交界處的農(nóng)田分布區(qū)，碳匯強度約為0—30gC·m-2·a-1(圖5 a)。

圖5 1961—2010年東北地區(qū)NEP的空間分布(a)及其變化趨勢(b)Fig.5 Spatial distribution(a)and trend(b)of annual NEP during 1961—2010 in the Northeast

從NEP的變化趨勢來看，東北地區(qū)NEP的變化趨勢與其空間分布有良好的一致性。NEP增加地區(qū)的面積占研究區(qū)域總面積的31.8%，而顯著增加的面積僅占總面積29.3%，主要集中在三江平原地區(qū)，內(nèi)蒙古高原等碳源區(qū)。NEP減少地區(qū)的面積占研究區(qū)域總面積的66.8%，其中顯著減少的面積占研究區(qū)域總面積的55%，主要在長白山西麓，大、小興安嶺等碳匯區(qū)(圖5 b)?？傮w來看，碳源區(qū)向大氣釋放的碳量在減少，碳匯區(qū)從大氣吸收的碳也在減少。

東北地區(qū)NEP與溫度呈顯著負相關(guān)(P＜0.05)的面積占31.6%，主要分布在大興安嶺東麓，小興安嶺以及長白山地區(qū)(圖6)。NEP與降水呈顯著正相關(guān)(P＜0.05)的面積占研究區(qū)域總面積的91.5%，只有大、小興安嶺和長白山的部分地區(qū)與降水的相關(guān)性不顯著(圖6)。因此，在空間上，東北地區(qū)NEP空間分布主要由降水決定，只有一小部分地區(qū)的NEP受溫度控制。

3 討論

1961—2010年，東北地區(qū)凈生態(tài)系統(tǒng)生產(chǎn)力(NEP)單位平均值為22.5gC·m-2·a-1，低于在西南高山地區(qū)NEP值的29.7gC·m-2·a-1［36］。東北地區(qū)固定的總碳量為0.026PgC/a，占全國NEP總量的15%—37%［6，67-68］。在年代際變化量上，Cao 在中國區(qū)域內(nèi)的研究結(jié)果表明，與20世紀(jì)80年代相比，東北地區(qū)NEP總量在90年代的減少，抵消了南部地區(qū)NEP的增加，本研究中，80年代到90年代，東北地區(qū)NEP總量減少了39.45TgC/a，與Cao的研究結(jié)果一致［6］，而且在2000年以后，東北地區(qū)NEP總量仍在減少，但變化量很小，為2.51TgC/a。

NEP的年際變化與氣候所導(dǎo)致的凈初級生產(chǎn)力(NPP)和異養(yǎng)呼吸(HR)的變化密切相關(guān)。研究時段內(nèi)，東北地區(qū)增溫明顯，年平均NPP與溫度不相關(guān)，HR與溫度有極顯著的正相關(guān)關(guān)系(r=0.505，P＜0.01)，溫度每升高1℃，東北地區(qū)由土壤異養(yǎng)呼吸向大氣釋放的碳增加14.85Tg(圖7)。NPP與降水有極顯著的正相關(guān)關(guān)系(r=0.942，P＜0.01)，降水量每增加100mm，東北地區(qū)的植被凈初級生產(chǎn)力增加89.4Tg，而HR與降水的相關(guān)關(guān)系不顯著(圖7)。這表明東北地區(qū)降水不足，是植被生長的限制因子，這與國志興等基于EOS/MOS衛(wèi)星遙感資料的研究，分析東北地區(qū)年均NPP與溫度和降水量的相關(guān)性結(jié)論一致［23］。

圖6 1961—2010年東北地區(qū)NEP與溫度、降水的相關(guān)關(guān)系Fig.6 Correlation coefficient between annual NEP and mean annual temperature and mean annual precipitation during 1961—2010 in the Northeast

圖7 東北地區(qū)NPP、HR與溫度、降水的相關(guān)關(guān)系Fig.7 Sensitivity to temperature and precipitation of NPP and HR in the Northeast

NPP和HR共同決定了NEP的年際變化，年平均NEP與NPP和HR的相關(guān)系數(shù)分別為0.867(P＜0.05)和-0.338(P＜0.05)。這表明東北地區(qū)的陸地碳匯主要是由于降水量的變化影響植被生產(chǎn)力引起的，與之相比，溫度對土壤呼吸作用對陸地碳匯的貢獻相對較小，這與在東亞地區(qū)的研究存在差異［12］，在東亞地區(qū)，溫度和太陽輻射的年際變化是導(dǎo)致區(qū)域碳收支年際變化的主要原因，可能原因是東亞大部分地區(qū)處于濕潤的季風(fēng)氣候區(qū)，降水量的變化和水分脅迫沒有顯著影響，NPP對降水的敏感性很低，而東北地區(qū)處于歐亞大陸的東緣，自東向西，氣候上從濕潤區(qū)、半濕潤區(qū)過渡到半干旱區(qū)，降水是植被生長的限制因子。

東北地區(qū)氣候變暖導(dǎo)致HR增加，而NPP隨著降水量的年際變化呈波動趨勢(圖8)，兩者的綜合作用導(dǎo)致東北地區(qū)NEP的年際變化沒有明顯的線性變化趨勢，而是與降水量的年際波動趨勢相一致。20世紀(jì)60和70年代，東北地區(qū)碳吸收量分別占研究時段總量的21%和20%，80年代的碳吸收量顯著增加，這與80年代東北地區(qū)降水量增加導(dǎo)致的NPP增加量大于溫度升高造成土壤呼吸的增加量有關(guān)，最終導(dǎo)致NEP呈增長趨勢(圖8)。進入90年代，尤其是在90年代后期，東北地區(qū)由于降水量下降NPP減少趨勢明顯，同期的土壤呼吸量卻顯著增加，碳吸收強度有所下降。21世紀(jì)的前10a，東北地區(qū)的溫度繼續(xù)升高，降水顯著減少，碳吸收總量較90年代減少了2.51TgC/a，這一時期中的 2000年和2007年，是研究時段內(nèi)凈生態(tài)系統(tǒng)生產(chǎn)力(NEP)最低的兩年，降水量分別為465mm和411mm，低于降水量的多年平均值，因而NPP水平很低，而且由于溫度的升高伴隨著土壤呼吸的增加，導(dǎo)致這兩年東北地區(qū)凈碳釋放量達-0.094PgC和-0.093PgC，約是1961—2010年年平均碳吸收的3.6倍，但這期間東北地區(qū)陸地生態(tài)系統(tǒng)仍起著碳匯作用。因此，近20年，升溫伴隨降水減少是東北區(qū)域碳匯能力減弱的重要原因。

圖8 1961—2010年東北地區(qū)NPP和HR的年際變化Fig.8 Interannual variations in NPP and HR in the Northeast

NEP的空間分布和變化趨勢與溫度和降水也有密切關(guān)系。東北平原、內(nèi)蒙古高原以及三江平原的部分區(qū)域，是東北地區(qū)主要的碳源區(qū)，在過去50a，由于降水與凈初級生產(chǎn)力呈正相關(guān)關(guān)系，降水的減少造成NPP減少，同時溫度的降低使得土壤呼吸作用減弱，但降水引起的NPP減少量小于溫度降低造成的土壤呼吸減少量，使得過去50a這些地區(qū)的NEP呈增加趨勢，即碳源區(qū)向大氣釋放的碳量在減少。大、小興安嶺和長白山等主要碳匯區(qū)，大部分地區(qū)NEP與溫度呈極顯著負相關(guān)，與降水正相關(guān)，低溫或多雨是這些地區(qū)形成碳匯的主要原因。研究時段內(nèi)，碳匯區(qū)溫度升高，降水減少，導(dǎo)致土壤呼吸增高，NPP減少，NEP呈減少趨勢，即碳匯區(qū)從大氣吸收的碳量在減少。

4 結(jié)論

(1)過去50a東北區(qū)域NEP沒有明顯的線性變化趨勢，NEP和年降水量的變化規(guī)律基本一致，即同期上升或達到最高值。東北地區(qū)的陸地碳匯主要是由于降水量的變化影響植被生產(chǎn)力引起的，與之相比，溫度對土壤呼吸作用對陸地碳匯的貢獻相對較小。

(2)NEP的空間分布呈現(xiàn)出東部高，西部和中部低，北部高，南部低的空間趨勢，并隨降水量和溫度分布呈現(xiàn)區(qū)域異質(zhì)性。降水是決定NEP空間分布的最主要因子。

(3)升溫伴隨降水增加導(dǎo)致1961—1990年NEP呈增加趨勢，而其后升溫伴隨降水減少則是近20年東北區(qū)域碳匯能力減弱的重要原因。

［1］Canadell J G，Mooney H A，Baldocchi D D，Berry J A，Ehleringer J R，F(xiàn)ield C B，Gower S T，Hollinger D Y，Hunt J E，Jackson R B，Running S W，Shaver G R，Steffen W，Trumbore S E，Valentini R，Bond B Y.Carbon metabolism of the terrestrial biosphere: A multitechnique approach for improved understanding.Ecosystems，2000，3(2):115-130.

［2］Schimel D S，House I J，Hibbard K A，Bousquet P，Ciais P，Peylin P，Braswell B H，Apps M J，Baker D，Bondeau A，Canadell J，Churkina G，Cramer W，Denning A S，F(xiàn)ield C B，F(xiàn)riedlingstein P，Goodale C，Heimann M，Houghton R A，Melillo J M，MooreⅢ B，Murdiyarso D，Noble I，Pacala S W，Prentice I C，Raupach M R，Rayner P J，Scholes R J，Steffen W L，Wirth C.Recent patterns and mechanisms of carbon exchange by terrestrial ecosystems.Nature，2001，414(6860):169-172.

［3］IGBP Terrestrial Carbon Working Group.The terrestrial carbon cycle:implications for the Kyoto Protocol.Science，1998，280(5368):1393-1394.

［4］Take K R，Scott N A，Parshotam A，Brown L，Wilde R H，Giltrap D J，Trustrum N A，Gomez B，Ross D J.A multi-scale analysis of a terrestrial carbon budget:Is New Zealand a source or sink of carbon?.Agriculture，Ecosystems ＆ Environment，2000，82(1/3):229-246.

［5］Cao M K，Prince S D，Tao B，Small J，Li K R.Regional pattern and interannual variations in global terrestrial carbon uptake in response to changes in climate and atmospheric CO2.Tellus B，2005，57(3):210-217.

［6］Cao M K，Prince S D，Li K R，Tao B，Small J，Shao X M.Response of terrestrial carbon uptake to climate interannual variability in China.Global Change Biology，2003，9(4):536-546.

［7］Mohamed M A A，Babiker I S，Chen Z M，Ikeda K，Ohta K，Kato K.The role of climate variability in the inter-annual variation of terrestrial net primary production(NPP).Science of the Total Environment，2004，332(1):123-137.

［8］Lu X L，Zhuang Q L.Evaluating climate impacts on carbon balance of the terrestrial ecosystems in the Midwest of the United States with a process-based ecosystem model.Mitigation and Adaptation strategies Global Change，2010，15(5):467-487.

［9］Xu X F，Tian H Q，Wan S Q.Climate warming impacts on carbon cycling in terrestrial ecosystems.Journal of Plant Ecology(Chinese Version)，2007，31(2):175-188.

［10］Woodwell G M，Whitaker R H，Reiners W A，Likens G E，Delwich C C，Botkin D B.Biota and the world carbon budget.Science;(United States)，1978，199(4325):141-146.

［11］Tian H Q，Melillo J M，Kicklighter D W，McGuire A D，HelfrichⅢ J V K，MooreⅢ B，Vorosmarty C J.Effect of interannual climate variability on carbon storage in Amazonian ecosystems.Nature，1998，396(6712):664-667.

［12］Ito A.The regional carbon budget of East Asia simulated with a terrestrial ecosystem model and validated using AsiaFlux data.Agricultural and Forest Meteorology，2008，148(5):738-747.

［13］Piao S L，Tan K，Nan H J，Ciais P，F(xiàn)ang J Y，Wang T.Impacts of climate and CO2changes on the vegetation growth and carbon balance of Qinghai-Tibetan grasslands over the past five decades.Global and Planetary Change，2012，98-99:73-80.

［14］Gurney K R，Law R M，Denning A S，Rayner P J，Baker D，Bousquet P，Bruhwiler L，Chen Y H，Ciais P，F(xiàn)an S M，F(xiàn)ung I Y，Gloor M，Heimann M，Higuchi K，John J，Maki T，Maksyutov S，Masarie K，Peylin P，Prather M，Pak B C，Randerson J，Sarmiento J，Taguchi S，Takahashi T，Yuen C W.Towards robust regional estimates of CO2sources and sinks using atmospheric transport models.Nature， 2002， 415(6872):626-630.

［15］Houghton R A.Why are the estimates of the terrestrial carbon balance so different?.Global Change Biology，2003，9(4):500-509.

［16］Pacala S W，Hurtt G C，Baker D，Peylin P，Houghton R A，Birdsey R A，Heath L，Sundquist E T，Stallard R F，Ciais P，Moorcroft P，Caspersen J P，Shevliakova E，Moore B，Kohlmaier G，Holland E，Gloor M，Harmon M E，F(xiàn)an S M，Sarmiento J L，Goodale C L，Schimel D，F(xiàn)ield C B.Consistent land-and atmosphere-based U.S.carbon sink estimates.Science，2001，292(5525):2316-2320.

［17］Tao B，Cao M K，Li K R，Gu F X，Ji J J，Huang M，Zhang L M.Spatial patterns of terrestrial net ecosystem productivity in China during1981—2000.Science in China (SeriesD:Earth Sciences)，2007，50(5):745-753.

［18］Liu W G.Research of carbon and NPP changes of the land ecosystem in XinJing［D］.Wulumuqi:Xinjiang University，2007.

［19］Yu L，Zhang B.The degradation situations of black soil in China and its prevention and counter measures.Journal of Arid Land Resources and Environment，2004，18(1):99-103.

［20］Zhu F，LiuZ M，Wang Z M，Song K S.Temporal-spatial characteristics and factors influencing crop NPP across NortheasternChina.ResourcesScience， 2010， 32(11):2079-2084.

［21］Yu Y，F(xiàn)an W Y，Li M Z.Forest carbon rates at different scales in Northeast China forest area.Chinese Journal of Applied Ecology，2012，23(2):341-346.

［22］Wang S Q，Zhou C H，Liu J Y，Li K R，Yang X M.Stimulation analyses of terrestrial carbon cycle balance model in Northeast China.Acta Geographica Sinica，2001，56(4):390-400.

［23］Guo Z X，Wang Z M，Zhang B，Liu D W，Yang G，Song K S，Li F.Analysis oftemporal-spatialcharacteristics and factors influencing vegetation NPP in Northeast China from 2000 to 2006.Resources Science，2008，30(8):1226-1235.

［24］Liu X，Guo Q X，Liu J W.Simulation of spatial pattern of soil respiration in forests in the eastern part of Northeast China.Journal of Northeast Forestry University.2010，28(12):58-64.

［25］Liu X，Guo Q X，Liu J W.Sensitivity analysis for main factors influencing NPP of forests simulated by IBIS in the eastern area of Northeast China.ActaEcologicaSinica， 2011， 31(7):1772-1782.

［26］Wang J B.Chinese terrestrial net ecosystem productive model applied remote sensing data［D］. Hangzhou:Zhejiang University，2004.

［27］Zhao J F，Yan X D，Jia G S.Simulating the responses of forest net primary productivity and carbon budget to climate change in Northeast China.Acta Ecologica Sinica，2008，28(1):92-102.

［28］Cao M K，Woodward F I.Net primary and ecosystem production and carbon stocks of terrestrial ecosystems and their responses to climate change.Global Change Biology，1998，4(2):185-198.

［29］Cao M K，Woodward F I.Dynamic responses of terrestrial ecosystem carbon cycling to global climate change.Nature，1998，393(6682):249-252.

［30］Woodward F I，Smith T M，Emanuel W R.A global land primary productivityandphytogeographymodel.Globalbiogeochemical cycles，1995，9(4):471-490.

［31］Cao M K，Prince S D，Shugart H H.Increasing terrestrial carbon uptake from the 1980s to the 1990s with changes in climate and atmospheric CO2.GlobalBiogeochemicalcycles， 2002， 16(4):1069.

［32］Cao M K，Zhang Q F，Shugart H H.Dynamic responses of African ecosystem carbon cycling to climate change.Climate Research，2001，17(2):183-193.

［33］Cao M K，Tao B，Li K R，Shao X M，Prience S D.Interannual variation in terrestrial ecosystem carbon fluxes in China from 1981—1998.Acta Botanica Sinica(Chinese Edition)，2003，45(5):552-560.

［34］Tao B.Modeling research on terrestrial net primary productivity＆net ecosystem productivity in China［D］.Beijing:Graduate University of Chinese Academy of Sciences，2003.

［35］Li K R，Wang S Q，Cao M K.Carbon storage of vegetation and soil in china.Science in China(Series D)，2003，33(1):72-80.

［36］Pang R，Gu F X，Zhang Y D，Hou Z H，Liu S R.Temporalspatial variations of net ecosystem productivity in alpine area of southwestern China.Acta Ecologica Sinica，2012，32(24):7844-7856.

［37］Gu F X，Cao M K，Wen X F，Liu Y F，Tao B.A comparison between simulated and measured CO2and water flux in a subtropical coniferous forest.Science in China(Series D:Earth Sciences)，2006，49(SupplⅡ):241-251.

［38］Zhang L，Yu G R，Gu F X，He H L，Zhang L M，Han S J.Uncertainty analysis of modeled carbon fluxes for a broad-leaved Korean pine mixed forest using a process-based ecosystem model.Journal of Forest Research，2012，17(3):268-282.

［39］Harley P C，Thomas R B，Reynolds J F，Strain B R.Modelling photosynthesis of cotton grown in elevated CO2.Plant，Cell ＆Environment，1992，15(3):271-282.

［40］Raich J W，Rastetter E B，Melillo J M，Kicklighter D W，Steudler P A，Peterson B J，Grace A L，MooreⅢ B，Vorosmarty C J.Potentialnetprimaryproductivityin south America:Application of a global model.Ecological Applications，1991，1(4):399-429.

［41］McGuire A D，Melillo J M，Joyce L A，Kicklighter D W，Grace A L，MooreⅢ B，Vorosmarty C J.Interactions between carbon and nitrogen dynamics in estimating net primary productivity for potential vegetation in North America.Global Biogeochemical Cycles，1992，6(2):101-124.

［42］Gu F X，Zhang Y D，Tao B，Wang Q F，Yu G R，Zhang L M，Li K R.Modeling the effects of nitrogen deposition on carbon budget in two temperate forests.Ecological Complexity，2010，7(2):139-148.

［43］Mao D H，Wang Z M，Luo L，Han J X.Spatio-temporal pattern of net primary productivity and its driven factors in Northeast China in 1982—2010.Scientia Geographica Sinica，2012，32(9):1106-1111.

［44］Xu W D，He X Y，Chen W，Liu C F，Zhao G L，Zhou Y.Ecological division of vegetation in Northeast China.Chinese Journal of Ecology，2008，27(11):1853-1860.

［45］Lloyd J，Grace J，Miranda A C，Meir P，Wong S C，Miranda H S，Wright I R，Gash J H C，McIntyre J.A simple calibrated model of Amazon rainforest productivity based on leaf biochemical properties.PlantCellandEnvironment， 1995， 18(10):1129-1145.

［46］Baldocchi D D.Assessing the eddy covariance technique for evaluating carbon dioxide exchange rates of ecosystems:past，present and future.Global Change Biology，2003，9(4):479-492.

［47］Yu G R，Zhang L M，Sun X M，F(xiàn)u Y L，Li Z G.Advances in carbon flux observation and research in Asia.Science in China(Series D:Earth Sciences)，2005，48(SupplⅠ):1-16.

［48］Baldocchi D，F(xiàn)alge E，Gu L H，Olson R，Hollinger D，Running S，Anthoni P.FLUXNET:A new tool to study the temporal and spatial variability of ecosystem-scale carbon dioxide，water vapor，and energy flux densities.Bulletin of the American Meteorological Society，2001，82(11):2415-2434.

［49］Falge E，Baldocchi D，Tenhunen J，Aubinet M，Bakwin P，Berbigier P，Bernhofer C.Seasonality of ecosystem respiration and gross primary production as derived from FLUXNET measurements.Agricultural and Forest Meteorology，2002，113(1):53-74.

［50］Baldocchi D D，Wilson K B.Modeling CO2and water vapor exchange of a temperate broadleaved forest across hourly to decadal time scales.Ecological Modelling，2001，142(1):155-184.

［51］Gu F X，Cao M K，Yu G R，Tao B，Wen X F，Liu Y F，Zhang L M.Modeling carbon exchange in different forest ecosystems by CEVSA model:Comparison with eddy covariance measurements.Advances in Earth Science，2007，22(3):223-234.

［52］Luo T X.Patterns of net primary productivity for Chinese major forest types and their mathematical models［D］.Beijing:Chinese Academy of Sciences，1996.

［53］Zhao M.Zhou G S.A new methodology for estimating forest NPP based on forest inventory data ——A case study of Chinese pine forest.Journal of Forestry Research，2004，15(2):93-100.

［54］Wang Y H，Zhou G S，Jiang Y L，Yang Z Y.Estimating biomass and NPP ofLarix forests using forest inventory data(FID).Acta Phytoecologica Sinica.2001，25(4):420-425.

［55］Zhang F，Zhou G S.Spatial-temporal variation in net primary productivity along Northeast China Transect(NECT)from 1982 to 1999.Journal of Plant Ecology(Chinese Version).2008，32(4):798-809.

［56］Zhao G S，Wang J B，F(xiàn)an W Y，Ying T Y.Vegetation net primary productivity in NortheastChina in 2000—2008:Simulation and seasonal change.Chinese Journal of Applied Ecology，2011，22(3):621-630.

［57］Zhu J J，Yan Q L，F(xiàn)an A N.Soil carbon flux for secondary forest ecosystems in the northeast mountains，China//Proceeding of the Fifth Youth EcologistsConference.Guangdong:Guangdong Province Association for Science and Technology，2008:112-120.

［58］Jiang C S，Hao Q J，Song C C，Hu B Q.Effects of marsh reclamation on soil respiration in the Sanjiang Plain.Acta Ecologica Sinica.2010，30(17):4539-4548.

［59］Liu Y.Soil carbon efflux in the typical forest ecosystems of Changbai Mountains［D］.Beijing:Graduate University of Chinese Academy of Sciences，2007.

［60］Liu Q J，Zhang G C，Xu Q Q，Wang Y D，Wang H M.Soil respiration of two dominant tundra communities under controlled temperatures in ChangbaiMountain，NortheastChina.Acta Ecologica Sinica，2011，31:137-44.

［61］Han G X，Zhou G S，Xu Z Z.Seasonal dynamics of soil respiration and carbon budget of maize(Zea mays L.)farmland ecosystem.Chinese Journal of Eco-Agriculture，2009，17(5):874-879.

［62］Wang X，Zhou G S，Jiang Y L，Li F.Comparison of soil respiration in broad-leaved Korean pine forest and reclaimed cropland in Changbai Mountains，China.Journal of Plant Ecology.2006，30(6):887-893.

［63］Yang J Y，Wang C K.Soil carbon storage and flux of temperate forest ecosystems in northeastern China.Acta Ecologica Sinica，2005，25(11):2875-2881.

［64］Lu B，Wang S H，Mao Z J，Sun T，Jia G M，Jin S B，Sun P F，Cheng C X.Soil respiration characteristics of four primary Korean pine communities in growing season at Xiaoxing'an Mountain，China.Acta Ecologica Sinica，2010，30(15):4065-4074.

［65］Qi Y C，Dong Y S，Manfred D，Geng Y B，Liu L X，Liu X R.Comparison of CO2effluxes and their driving factors between two temperate steppesin InnerMongolia，China.Advancesin Atmospheric Sciences，2006，23(5):726-736.

［66］Qi Y C，Dong Y S，Liu L X，Liu X R，Peng Q，Xiao S S，He Y T.Spatial-temporal variation in soil respiration and its controlling factors in three steppes of Stipa L.in Inner Mongolia，China.Science China(Series D:Earth Sciences)，2010，53(5):683-693.

［67］Mu Q Z，Zhao M S，Running S W，Liu M L，Tian H Q.Contribution of increasing CO2and climate change to the carbon cycle in China's ecosystems.Journal of Geophysical Research:Biogeosciences(2005—2012)，2008，113(G01018):1-15.

［68］Fang J Y，Guo Z D，Piao S L，Chen A P.Terrestrial vegetation carbon sinks in China from 1981 to 2000.Science in China(Series D:Earth Sciences)，2007，37(6):804-812.

參考文獻:

［9］徐小鋒，田漢勤，萬師強.氣候變化對生態(tài)系統(tǒng)碳循環(huán)的影響.植物生態(tài)學(xué)報，2007，31(2):175-188.

［17］陶波，曹明奎，李克讓，顧峰雪，季勁鈞，黃玫，張雷明.1981—2000年中國陸地凈生態(tài)系統(tǒng)生產(chǎn)力空間格局及其變化.中國科學(xué)(D輯)，2006，36(12):1131-1139.

［18］劉衛(wèi)國.新疆陸地生態(tài)系統(tǒng)凈初級生產(chǎn)力和碳時空變化研究［D］.烏魯木齊:新疆大學(xué)，2007.

［19］于磊，張柏.中國黑土退化現(xiàn)狀與防治對策.干旱區(qū)資源與環(huán)境，2004，18(1):99-103.

［20］朱鋒，劉志明，王宗明，宋開山.東北地區(qū)農(nóng)田凈初級生產(chǎn)力時空特征及其影響因素分析.資源科學(xué)，2010，32(11):2079-2084.

［21］于穎，范文義，李明澤.東北林區(qū)不同尺度森林的含碳率.應(yīng)用生態(tài)學(xué)報，2012，23(2):341-346.

［22］王紹強，周成虎，劉紀(jì)遠，李克讓，楊曉梅.東北地區(qū)陸地碳循環(huán)平衡模擬分析.地理學(xué)報，2001，56(4):390-400.

［23］國志興，王宗明，張柏，劉殿偉，楊桄，宋開山，李方.2000年—2006年東北地區(qū)植被NPP的時空特征及影響因素分析.資源科學(xué)，2008，30(8):1226-1235.

［24］劉曦，國慶喜，劉經(jīng)偉.東北東部森林土壤呼吸空間分布格局的模擬.東北林業(yè)大學(xué)學(xué)報.2010，28(12):58-64.

［25］劉曦，國慶喜，劉經(jīng)偉.IBIS模擬東北東部森林NPP主要影響因子的敏感性.生態(tài)學(xué)報，2011，31(7):1772-1782.

［26］王軍邦.中國陸地凈生態(tài)系統(tǒng)生產(chǎn)力遙感模型研究［D］.杭州:浙江大學(xué)，2004.

［27］趙俊芳，延曉冬，賈根鎖.東北森林凈第一性生產(chǎn)力與碳收支對氣候變化的響應(yīng).生態(tài)學(xué)報，2008，28(1):92-102.

［34］陶波.中國陸地生態(tài)系統(tǒng)凈初級生產(chǎn)力和凈生態(tài)系統(tǒng)生產(chǎn)力模擬研究［D］.北京:中國科學(xué)院研究生院，2003.

［35］李克讓，王紹強，曹明奎.中國植被和土壤碳儲量.中國科學(xué)(D輯)，2003，33(1):72-80.

［36］龐瑞，顧峰雪，張遠東，侯振宏，劉世榮.西南高山地區(qū)凈生態(tài)系統(tǒng)生產(chǎn)力時空動態(tài).生態(tài)學(xué)報，2012，32(24):7844-7856.

［37］顧峰雪，曹明奎，溫學(xué)發(fā)，劉允芬，陶波.亞熱帶針葉林水碳通量的模擬及其與觀測的對比研究.中國科學(xué)(D輯)，2006，36(增刊I):224-233.

［43］毛德華，王宗明，羅玲，韓佶興.1982—2010年中國東北地區(qū)植被NPP時空格局及驅(qū)動因子分析.地理科學(xué)，2012，32(9):1106-1111.

［44］徐文鐸，何興元，陳瑋，劉常富，趙桂玲，周園.中國東北植被生態(tài)區(qū)劃.生態(tài)學(xué)雜志，2008，27(11):1853-1860.

［51］顧峰雪，曹明奎，于貴瑞，陶波，溫學(xué)發(fā)，劉允芬，張雷明.典型森林生態(tài)系統(tǒng)碳交換的機理模擬及其與觀測的比較研究.地理科學(xué)進展，2007，22(3):223-234.

［52］羅天祥.中國主要森林類型生物生產(chǎn)力格局及其數(shù)學(xué)模型［D］.北京:中國科學(xué)院地理科學(xué)與資源研究所，1996.

［54］王玉輝，周廣勝，蔣延玲，楊正宇.基于森林資源清查資料的落葉松林生物量和凈生長量估算模式.植物生態(tài)學(xué)報，2001，25(4):420-425.

［55］張峰，周廣勝.中國東北樣帶植被凈初級生產(chǎn)力時空動態(tài)遙感模擬.植物生態(tài)學(xué)報，2008，32(4):798-809.

［56］趙國帥，王軍邦，范文義，應(yīng)天玉.2000—2008年中國東北地區(qū)植被凈初級生產(chǎn)力的模擬及季節(jié)變化.應(yīng)用生態(tài)學(xué)報，2011，22(3):621-630.

［57］朱教君，閆巧玲，范阿南.東北山區(qū)次生林生態(tài)系統(tǒng)土壤碳通量研究//第五屆中國青年生態(tài)學(xué)工作者學(xué)術(shù)研討會論文集.廣東:廣東省科學(xué)技術(shù)協(xié)會，2008:112-120.

［58］江長勝，郝慶菊，宋長春，胡必琴.墾殖對沼澤濕地土壤呼吸速率的影響.生態(tài)學(xué)報，2010，30(17):4539-4548.

［59］劉穎.長白山典型森林生態(tài)系統(tǒng)土壤碳通量研究［D］.北京:中國科學(xué)院研究生院，2007.

［61］韓廣軒，周廣勝，許振柱.玉米農(nóng)田生態(tài)系統(tǒng)土壤呼吸作用季節(jié)動態(tài)與碳收支初步估算.中國生態(tài)農(nóng)業(yè)學(xué)報，2009，17(5):874-879.

［62］王旭，周廣勝，蔣延玲，李峰.長白山紅松針闊混交林與開墾農(nóng)田土壤呼吸作用比較.植物生態(tài)學(xué)報，2006，30(6):887-893.

［63］楊金艷，王傳寬.東北東部森林生態(tài)系統(tǒng)土壤碳貯量和碳通量.生態(tài)學(xué)報，2005，25(11):2875-2881.

［64］陸斌，王淑華，毛子軍，孫濤，賈桂梅，靳世波，孫鵬飛，程春香.小興安嶺4種原始紅松林群落類型生長季土壤呼吸特征.生態(tài)學(xué)報，2010，30(15):4065-4074.

［66］齊玉春，董云社，劉立新，劉杏認，彭琴，肖勝生，何亞婷.內(nèi)蒙古錫林河流域主要針茅屬草地土壤呼吸變化及其主導(dǎo)因子.中國科學(xué)(D輯)，2010，40(3):341-351.

［68］方精云，郭兆迪，樸世龍，陳安平.1981—2000年中國陸地植被碳匯的估算.中國科學(xué)(D輯)，2007，37(6):804-812.