舒海燕, 江洪,,3*, 陳曉峰, 孫文文, 馬錦麗, 張敏霞

?

安吉毛竹林非生長季CO2通量與其光合響應關系探究

舒海燕1, 江洪1,2,3*, 陳曉峰2, 孫文文2, 馬錦麗2, 張敏霞2

1. 西南大學三峽庫區(qū)生態(tài)環(huán)境教育部重點實驗室, 重慶400715 2. 浙江農林大學林業(yè)與生物技術學院, 浙江臨安311300 3. 南京大學國際地球系統(tǒng)科學研究所, 江蘇南京210046

利用渦度相關法和生理生態(tài)法(光合)對安吉毛竹林生態(tài)系統(tǒng)非生長季(2014.10—2015.4, 無2月)的CO2通量進行測定, 結果表明, 兩種方法具有良好的響應關系, 基于渦度的生態(tài)系統(tǒng)凈交換量(和基于光合的生態(tài)系統(tǒng)凈初級生產力(日變化均呈單峰型曲線, 受光合有效輻射()的直接影響, 在午間12:00前后有最大值, 且兩者均表現(xiàn)為秋季最大(10月最高), 春季次之, 冬季最低。就各月而言, 基于光合測定的日最大值和月總量均低于渦度相關技術白天的測定值, 但均占到后者的40%和65%以上, 在毛竹光合能力比較強的10、11和4月, 所占比例更大。表明在白天, 安吉毛竹林生態(tài)系統(tǒng)的主要取決于毛竹的光合作用。此外, 觀測期間, 毛竹林生態(tài)系統(tǒng)在的凈初級生產力總量為1365.59 gCO2·m–2, 吸收的碳約372.44 gC·m–2, 整體表現(xiàn)為較強的碳匯。

毛竹林生態(tài)系統(tǒng); 非生長季; 渦度相關技術; 光合作用; 凈初級生產力

1 前言

隨著社會的進步、經(jīng)濟水平的提升以及工業(yè)產業(yè)迅速發(fā)展, 越來越多的環(huán)境問題日益突出, 尤其是全球氣候變暖吸引著全社會人類的目光, 不止對我國的經(jīng)濟、社會及生態(tài)環(huán)境造成了重大的影響, 也是世界社會經(jīng)濟可持續(xù)發(fā)展所面臨的嚴峻挑戰(zhàn), 同時嚴重威脅著人類的生活及安全, 正如所指: “溫室效應”是工業(yè)革命以來影響范圍最廣的生態(tài)環(huán)境問題之一[1]。大氣中CO2等溫室氣體濃度不斷升高是引起氣候變暖的重要原因, 對全球變暖的貢獻率達60%[2]。森林生態(tài)系統(tǒng)作為陸地生態(tài)系統(tǒng)的主體, 覆蓋面積約1/3左右[3], 因其高大的植被特性和復雜的結構, 使之具有更高的碳密度貯量和最高的生產力以及生物量積累, 使得森林生態(tài)系統(tǒng)成為吸收和儲存CO2的有效碳庫, 具有很強的碳匯能力[4]。我國自20世紀70年代以來, 人工林生態(tài)系統(tǒng)面積逐漸增加且居世界首位, 已成為我國森林生態(tài)系統(tǒng)的重要組成部分[5–8], 所以量化人工林生態(tài)系統(tǒng)的固碳對研究陸地生態(tài)系統(tǒng)碳循環(huán)以及全球碳的源/匯具有相當重要的意義。毛竹林作為典型的人工林生態(tài)系統(tǒng), 具有較高的固碳能力[9]。毛竹還具有生長快、產量高、用途廣、效益大和持續(xù)更新等特點, 在竹類資源中占有最大比例, 現(xiàn)成為我國分布最多、面積最大、利用領域最廣的經(jīng)濟竹種[10–11]。此外, 由于毛竹林根鞭發(fā)達, 更新頻繁, 對區(qū)域防風固沙、固碳排氧和調節(jié)氣候等方面發(fā)揮著巨大的生態(tài)保護作用, 具有較高的生態(tài)效益[12–13]。但是, 在非生長季, 主要由毛竹的光合作用轉向植被、凋落物和土壤微生物等的呼吸作用, 使得森林碳吸收降低, 甚至表現(xiàn)為碳釋放, 很大程度上影響了毛竹林生態(tài)系統(tǒng)全年的碳源/匯強度。因此, 研究毛竹林非生長季的CO2通量顯得尤為重要。

關于生態(tài)系統(tǒng)CO2通量研究的方法較多, 主要的是以微氣象學為基礎的渦度協(xié)方差法和生理生態(tài)學方法。其中最為直接的可連續(xù)測定的方法是渦度協(xié)方差法, 1986年被首次用于森林生態(tài)系統(tǒng)CO2通量的研究中[14], 發(fā)展至今, 盡管還存在著一些不足之處, 該方法仍然作為現(xiàn)今碳通量研究的一個標準方法[15–17]; 而生理生態(tài)學方法主要是使用測定室對森林生態(tài)系統(tǒng)不同組分分別進行測定并通過數(shù)學方法進行總量估計的方法, 該法能夠細化生態(tài)系統(tǒng)的碳平衡于不同組分之間[18]。關于渦度協(xié)方差法, 運用較廣, 如草原[19]、荒漠[20]、稻田[21]、針葉林[22–25]、針闊混交林[26–29]、闊葉混交林[30–32]、熱帶季雨林[33]、雷竹林[34]、毛竹林[9][13]等等, 但是將生理生態(tài)法與微氣象法同時結合運用的研究較少。本研究將以安吉毛竹林為研究對象, 同時結合兩種研究方法研究其非生長季(2014.10—2015.04)的CO2通量, 探究毛竹林非生長季CO2通量與其光合的響應關系, 明確毛竹光合在生態(tài)系統(tǒng)凈交換量(net ecosystem exchange,), 乃至凈初級生產力(net ecosystem productivity,)中的貢獻, 以便揭示毛竹林生態(tài)系統(tǒng)的生理生態(tài)學基礎, 從而解釋安吉毛竹林生態(tài)系統(tǒng)的碳源/匯能力。

2 研究地區(qū)與研究方法

2.1 研究地區(qū)概況

研究區(qū)位于浙江省湖州市安吉縣山川鄉(xiāng)(30°28′34.5″N, 119°40′25.7″E), 為典型的亞熱帶季風氣候, 全年氣候溫和、雨水充沛、光照充足、四季分明。歷史年均氣溫14.6 ℃, 年日照1980—2050 h, 年無霜期203—226 d, 年降雨日140—170 d, 降水量達1100—1900 mm。降水量由北至南遞增, 6月和8月雨量最多, 12月雨量最少, 具明顯的雙峰型降雨的特點, 各年份梅雨量的實際差異大, 是造成旱澇的原因[13]。該地區(qū)森林覆蓋率達到71%, 擁有山林13.2萬hm2, 其中竹林面積6.33萬hm2, 毛竹林面積為4.99萬hm2, 占森林總面積37.8%, 竹林總面積的78.8%。觀測站點海拔高度380 m, 站點下墊面坡度為2.5°—14.0°, 以觀測塔為中心1000 m范圍內主要森林類型為毛竹林, 土壤類型主要為黃壤。試驗區(qū)內的毛竹林為粗放型人工純林, 竹林密度為4500株·hm–2, 胸徑12—18 cm, 高度16—20 m, 冠幅3.0 m × 2.5 m, 蓋度90%, 郁閉度0.9, 毛竹節(jié)間短, 壁厚, 一般3—4月為發(fā)筍期, 5月毛竹開始抽葉, 6月份大量展開, 7—9月份為主要生長季, 生長能力較強。

2.2 研究方法

2.2.1 CO2通量的測定

本研究采用了開路式渦度相關系統(tǒng)對安吉毛竹林生態(tài)系統(tǒng)的CO2通量進行連續(xù)測定。觀測林地微氣象觀測塔高40 m, 開路渦度相關系統(tǒng)(OPEC)探頭的安裝高度為距地面38 m, 由三維超聲風速儀(CAST3, Campbell Inc., USA)和開路CO2/H2O分析儀(Li-7500, Li-Cor Inc., USA)組成, 前者用于測定風速和溫度脈動, 后者可測定CO2/H2O密度脈動, 并在線計算并存儲30 min CO2通量(), 其原始采樣頻率10 Hz, 通過數(shù)據(jù)采集器(CR1000, Campbell Inc., USA)存儲數(shù)據(jù)。此外, 常規(guī)氣象觀測系統(tǒng)包括7層風速(010C, metone, USA)、7層大氣溫度和濕度(HMP45C, Vaisala, Helsinki, Finland), 安裝高度分別為1 m、7 m、11 m、17 m 、23 m、30 m和38 m; 光合有效輻射和凈輻射傳感器安裝高度38 m, 用于采集光合有效輻射、上行/下行的長波/短波輻射、凈輻射數(shù)據(jù)。常規(guī)氣象觀測系統(tǒng)數(shù)據(jù)采樣頻率為5 Hz, 通過數(shù)據(jù)采集器(CR1000, Campbell Inc., USA)每30 min自動記錄平均風速、濕度、溫度、氣壓、凈輻射等常規(guī)氣象信息。本研究中空氣溫度、空氣濕度與飽和水汽壓差選用23 m高冠層氣象數(shù)據(jù)。儀器設備詳情可參考孫成[9]、劉玉莉[12]和楊爽[13]相關研究報道。

關于植被與大氣間CO2通量的計算, 用植被與大氣間凈生態(tài)系統(tǒng)CO2交換量(, mg·m–2·s–1)表示, 定義為:

(2)

式中:c(mg·m–2·s–1)為CO2湍流通量, 通過利用三維超聲風速儀測定三維風速, 以及紅外線CO2/H2O氣體分析儀測定CO2密度脈動, 再采用EdiRe ()軟件進行協(xié)方差運算, 并進行二次坐標軸旋轉、WPL()修正等, 最終將利用數(shù)據(jù)采集器記錄并儲存的10Hz的原始數(shù)據(jù)轉化為30 min數(shù)據(jù)而得到。具體數(shù)據(jù)處理方法詳見文獻[13]。s(mg·m–2·s–1)為觀測高度以下的碳儲存通量;為觀測高度(m);c為高度h處前后兩次相鄰時間測定的CO2濃度差值(mg·m–3);t為前后兩次測定的時間間隔(s)。當CO2從生態(tài)系統(tǒng)排放到大氣時, 定義符號為正, 當CO2從大氣進入到生態(tài)系統(tǒng)時, 定義為負。

2.2.2 葉片光響應曲線的測定

利用便攜式光合測定儀LI-6400XT(Li-cor, USA)對毛竹林冠層葉片的光響應曲線進行測定。于2014年10月—2015年4月間每月選定一個晴天, 就毛竹林冠層中上部的葉片進行光合測定, 測量時, 人在鐵塔上, 選用LED 2 cm×3 cm葉室。每次至少測定5個重復, 最后求取平均值。得到的光響應曲線采用IBM SPSS Statistics 21.0軟件通過非直角雙曲線公式進行擬合:

式中:為光合有效輻射(μmol·m–2·s–1),n為葉片在某一值下的凈光合速率(μmol·m–2·s–1),max為光飽和時最大凈光合速率(μmol·m–2·s–1),d為植物的暗呼吸速率(μmol·m–2·s–1)。為表觀量子效率(μmol·μmol–1),為曲線曲率系數(shù), 在0—1之間。用公式(3)擬合出這些參數(shù)之后, 再根據(jù)通量塔得到的值, 計算出各月白天(≥2 μmol·m–2·s–1)n的平均日變化。

根據(jù)公式(1)可以計算出林分的值, 在生理生態(tài)學方法中, 與之對應的是凈生態(tài)系統(tǒng)生產力(, mg·m–2·s–1)[18]。符號的定義與植被/大氣間凈生態(tài)系統(tǒng)CO2交換量()正好相反, 當CO2從大氣進入到生態(tài)系統(tǒng)時, 定義為正; 當CO2從生態(tài)系統(tǒng)排放到大氣時, 定義符號為負。但它們絕對值相等, 即。在此,可以粗略地計算為毛竹葉片的凈光合速率與當月葉面積指數(shù)(leaf index area,)的乘積。即:

式中0.044為單位換算系數(shù), 即1 μmol·m–2·s–1=0.044 mg·m–2·s–1。這樣, 白天n的平均日變化就可以轉化為白天的平均日變化, 從而與通量的日變化進行比較。然而, 因為溫度和天氣等原因, 2015年2月數(shù)據(jù)空缺。

2.3 數(shù)據(jù)處理

首先, 植被和大氣間凈生態(tài)系統(tǒng)CO2交換()在白天(≥2 μmol·m–2·s–1)僅僅來源于生態(tài)系統(tǒng)植被的光合作用[35]。同植物葉片的光響應曲線一樣,對也存在一定的響應關系, 且符合直角雙曲線方程, 其函數(shù)通?？梢杂肕ichaelis-Menten模型表達[36–37]:

公式(5)中,為白天凈生態(tài)系統(tǒng)交換量(mg·m–2·s–1), 為表觀初始初始光能利用率(表觀量子效率, mg·μmol–1),為光合有效輻射(μmol·m–2·s–1),max為光飽和時最大凈光合速率(mg·m–2·s–1),d為光照強度趨近于0時植物的呼吸速率(暗呼吸, mg·m–2·s–1)。

其次, 利用IBM SPSS Statistics 21.0軟件擬合毛竹葉片的光響應曲線, 進而擬合對的響應曲線, 并得到d、max和等參數(shù)。

最后, 采用Excel 2003對數(shù)據(jù)進行處理, 同時運用統(tǒng)計分析軟件IBM SPSS Statistics 21.0對數(shù)據(jù)進行顯著性分析。作圖采用Origin 8.6軟件。數(shù)據(jù)以平均數(shù)±標準差(Mean ± SD)形式呈現(xiàn)。

3 結果與分析

3.1 毛竹林生態(tài)系統(tǒng)氣候因子變化特征

試驗區(qū)屬于亞熱帶季風氣候, 全年氣候溫和、雨水充沛、光照充足、四季分明。就2014年10月—2015年4月毛竹林生態(tài)系統(tǒng)的光合有效輻射(photosynthetic active radiation,)、空氣溫度(air temperature,a, 23m)、空氣相對濕度(relative humidity,, 23m)和飽和水汽壓差(vapor pressure deficit,, 23m)的動態(tài)變化特征進行分析, 如圖1所示。

由圖1(a)可知, 2014年10月—2015年4月的光合有效輻射()先降低后升高, 但因為2014年11月陰雨天氣較多, 光照被削弱,值偏低; 此外, 冬季值相對較低, 1月最低, 平均為181.74± 268.01 μmol·m–2·s–1, 10月和4月相對較高, 平均分別達到311.10±434.98 μmol·m–2·s–1和349.45± 440.77 μmol·m–2·s–1, 入春之后,值繼續(xù)增大?？諝鉁囟?)變化趨勢呈“U”型曲線(圖1(b)), 10月最高, 為17.53±1.18 ℃, 最小值出現(xiàn)在12月, 僅4.23± 1.78 ℃, 入春后,a迅速升高。相對濕度()和飽和水汽壓差()變化趨勢剛好相反(圖1(c, d))。受值的影響,和在12月存在轉折, 12月最低, 平均僅47.66%±5.99%, 表明12月份降雨最低, 符合冬季寒冷少雨特征; 而該月的稍微升高, 達到4.59±1.15 kPa。另外, 受和a共同影響,在2月最低, 為2.81±1.003 kPa, 2月以后逐漸上升, 表明2月之后, 毛竹林的水分利用能力(或者蒸騰作用)逐漸增強, 主要因為2015年為毛竹生理大年, 出筍量較大, 水分消耗相應較高。

3.2 毛竹林生態(tài)系統(tǒng)CO2通量的變化特征

統(tǒng)計得到2014年10月—2015年4月(不包括2015年2月, 下同)同時刻的數(shù)據(jù), 并進行平均處理, 即得到該時間段毛竹林生態(tài)系統(tǒng)CO2凈交換量()的月平均日變化, 如圖2。同時, 圖3以2014年10月和2015年1月為例, 比較了隨光照強度()的改變而發(fā)生的變化。

由圖2可知, 2014年10月—2015年4月各月的變化規(guī)律性較強, 均呈“U”型曲線變化, 各月CO2交換速率最大值出現(xiàn)在午間12:00前后, 夜間為正值(釋放CO2), 白天為負(吸收CO2), 冬季數(shù)值較低。各月均值也為負, 表明毛竹林生態(tài)系統(tǒng)在非生長季各月均表現(xiàn)為碳吸收, 10月吸收速率最大, 為0.12±0.29 mg·m–2·s–1, 1月吸收速率最低, 僅0.04±0.17 mg·m–2·s–1, 但各月間沒有顯著差異(>0.05); 表明毛竹林生態(tài)系統(tǒng)非生長季總體表現(xiàn)為一定的碳匯作用, 10月的碳匯作用較強, 也證明了毛竹林具有較強的固碳能力[9], 甚至在冬季還能夠表現(xiàn)為碳吸收。另外, 因為11月份的值偏低(圖1(a)), 限制了植被的光合作用, 使得其數(shù)值略低于12月。

圖3結合日出和日落顯示了與的響應變化, 可知, 每天的值呈現(xiàn)倒“V”變化, 并且與有著良好的負相關關系。綜合圖2和圖3, 可以看出早晨6:00(冬季6:30)之前和傍晚18:00(冬季17:30)以后,值為0, 此時, 毛竹林不能依靠光能進行光合作用吸收CO2, 而呼吸作用仍在進行, 導致繼續(xù)向外界釋放CO2, 使得毛竹林生態(tài)系統(tǒng)為正值; 然而, 受低溫、黑暗等的影響, 各月該時間段值較低, 維持在0—0.15 mg·m–2·s–1之間變化, 尤其是12—1月更加明顯, 其中1月最低, 平均僅0.08 mg·m–2·s–1, 而春秋兩季相對較高。說明該生態(tài)系統(tǒng)在此時間段(夜間)表現(xiàn)為較弱的碳源作用, 且有明顯的季節(jié)變化, 即春秋季高于冬季。早晨6:00(或6:30)以后, 伴隨著日出,值逐漸增大, 毛竹林開始通過光合作用從外界吸收CO2, 與此同時, 毛竹林生態(tài)系統(tǒng)慢慢從正值向負值轉變, 表明該生態(tài)系統(tǒng)由碳源逐漸變?yōu)樘紖R; 午間12:00前后,達到峰值, 同時, CO2的吸收速率也達到一天的最大值, 10月最高, 達到0.73 mg·m–2·s–1, 而1月僅0.38 mg·m–2·s–1, 表明此時毛竹林的光合作用為當天最強, 吸收的CO2最多; 再之后, 隨著時間的推移, 太陽西下,值逐漸減小, 毛竹光合作用逐漸減弱,也相應降低, 并慢慢在18:00(或17:30)后變?yōu)樨撝? 毛竹林又逐漸變成碳源。分析原因, 生態(tài)系統(tǒng)的變化主要是由植物的光合作用和生態(tài)系統(tǒng)呼吸引起的, 而白天最主要的是植物的光合作用, 它是植物從大氣中吸收固定CO2的主要途徑, 而值是植物光合作用的必須條件, 制約著光合作用的進行,值為0時, 植物沒有光合作用, 此時為正值; 而增大時, 植物的光合作用增強,逐漸變?yōu)樨撝? 因此說毛竹林生態(tài)系統(tǒng)取決于毛竹林的光合作用, 也間接取決于值, 從另一角度表明的變化基本取決于林木光合作用的變化[35]。

3.3 毛竹葉片光響應曲線及的變化特征

毛竹葉片的光合作用強度影響著該生態(tài)系統(tǒng)的變化, 而光合作用強度又主要受控于環(huán)境光照強度, 光響應曲線(light response curve)可以直接反映其凈光合速率與光照強度的關系, 從而使得光合作用和光照強度之間的關系變得易于理解和研究。圖4分析比較了2014年10月—2015年4月(不包括2015年2月)毛竹葉片的光響應曲線。

可知, 10月光合作用最高, 平均凈光合速率為1.41±0.82 μmol·m–2·s–1, 光補償點較低(<10 μmol·m–2·s–1),光飽和點較高(>500 μmol·m–2·s–1), 呼吸較低 (0.17 μmol·m–2·s–1), 而1月光合作用最低, 平均僅0.59±0.52 μmol·m–2·s–1。毛竹葉片12—4月暗呼吸偏高, 約0.55 μmol·m–2·s–1; 此外, 各月最大凈光合速率(max)由大到小依次為: 10月>11月>4月>3月>12月>1月, 其中最大值為2.15 μmol·m–2·s–1(10月), 約是最小值(1月, 約0.96 μmol·m–2·s–1)的2.2倍。表明冬季和早春時節(jié), 毛竹的光合作用較弱, 尤其冬季最低, 主要是因為此時間段環(huán)境溫度低, 影響著一些酶的生物活性, 從而限制著植物的生理活動; 此外, 入冬后, 光照時間縮短, 光照強度減弱, 也削弱了植物的光合作用。總之, 對比圖2可知毛竹葉片的光合作用強度與毛竹林生態(tài)系統(tǒng)有著一致的變化, 凈光合速率越大,數(shù)值也越高(負向越大), 反之, 也越低。

同時, 運用公式(3)和(4)計算出該非生長季各月白天(≥2 μmol·m–2·s–1)的凈初級生產力(), 設定夜間(=0)的為0, 如圖5?？梢姼髟掳滋斓木鶠檎? 并呈現(xiàn)倒“U”型變化, 直接受控于值的改變。因為受到飽和光照的影響, 盡管中午時分光照強烈, 但維持在最大值而沒有劇烈起伏變化, 10月最高, 1月最低, 各月依次為0.43、0.36、0.23、0.19、0.25和0.31 mg·m–2·s–1, 均低于渦度相關系統(tǒng)觀測結果。主要是因為渦度相關系統(tǒng)觀測范圍較廣, 觀測結果較為全面, 而基于光合的測定組分較為狹窄, 不包括除了毛竹葉片以外的樹干、枝條等進行光合作用吸收CO2的部分。另外, 受光照、溫度、濕度等環(huán)境因子的影響, 各月總體表現(xiàn)為秋季最高, 春季居中, 冬季最低, 與的季節(jié)變化一致, 表明兩者響應關系較好。

3.4 毛竹林生態(tài)系統(tǒng)和對的響應

由2.2可知,與有著顯著的負相關關系, 而這種關系可以通過Michaelis-Menten方程來描述和體現(xiàn)。就2014年10月—2015年4月(無2月)白天(≥2 μmol·m–2·s–1)的對的響應曲線進行擬合, 結果如圖6所示。同時利用SPSS軟件對白天和進行擬合, 所得擬合參數(shù)結果見表1。

據(jù)圖6可知, 12月和1月份曲線比較接近,最大數(shù)值分別是0.44和0.38 mg·m–2·s–1, 相對應光合有效輻射分別是944.90和734.77 μmol·m–2·s–1; 秋季和春季的數(shù)值相對較高, 10、11和4月近似線性變化, 10月有最大值, 達到0.73 mg·m–2·s–1, 對應的光照強度為1244.45 μmol·m–2·s–1, 4月達到0.60 mg·m–2·s–1, 對應值為1203.83 μmol·m–2·s–1?？梢姸镜墓庹諒姸热? 日照時間短, 環(huán)境溫度低, 毛竹林吸收CO2的能力也相應減弱, 隨著春季的到來, 光照強度又逐漸增加, 毛竹林吸收CO2的能力也逐漸增強。

采用Michaelis-Menten方程, 非線性回歸擬合毛竹林生態(tài)系統(tǒng)非生長季不同月份白天對值的響應曲線, 擬合系數(shù)2均在0.93以上, 且均達到極顯著水平(圖6)。2014年10月—2015年4月生態(tài)系統(tǒng)暗呼吸速率d較低, 平均值變化為0.004—0.015 mg·m–2·s–1, 11月份偏高(表1), 主要是該月陰雨天較多, 云層較厚, 白天云層吸收熱量, 在夜間釋放出來, 使得生態(tài)系統(tǒng)內部溫度稍高, 利于植物的呼吸[38]; 12月之后偏低, 因為環(huán)境溫度較低, 直接影響著植物的生理活動, 導致其呼吸速率也相應較低, 而隨著溫度的升高, 毛竹的呼吸速率逐漸增大; 最大凈光合作用速率先降低后升高, 而初始光能利用率則先升高后下降, 4月最小, 為0.0005mg·μmol–1Photons, 表明光照強度較弱時, 毛竹的光合作用也較弱, 但是毛竹的初始光能利用率反而較高, 能捕獲更多的光能進行光合作用, 從而維持自身的生長。相關研究表明[18], 陰天植物的光能利用率更高, 主要因為陰天有較多的散射光及植物對散射光較高的光合利用效率; 相應地, 外界光照強度較大時, 毛竹的光合作用較強, 但是初始光能利用率卻相對降低, 因為此時有多余的光能供利用。

表1 白天NEE/NEP—PAR擬合參數(shù)

用非直角雙曲線擬合基于生理生態(tài)學方法測得的安吉毛竹林生態(tài)系統(tǒng)白天凈初級生產力與的響應曲線而得到d、max和(表1), 可知, 各月d先升高后降低, 變化范圍為0.002—0.012 mg·m–2·s–1, 低于Michaelis-Menten方程擬合的結果; 與d變化相反,max和先降低再升高, 明顯低于max, 但高于。主要是因為利用封閉葉室在可控條件下進行測定, 排除了其他因子的干擾, 就只有毛竹單一物種的光合, 對象是毛竹葉片, 且經(jīng)過光誘導, 在給予其不同的光照強度時, 其光能利用效率較高; 此外, 冬季寒冷干燥, 溫度和水分同時控制的植物的呼吸, 使得毛竹葉片的呼吸相對較低。而的測定是在自然條件下測定, 存在著諸多不可控因素, 就算是白天, 除了毛竹自身的CO2釋放外, 還有林下其他植物的光合和呼吸作用以及土壤呼吸等, 因此d和max較高?？傊?的測定的成分復雜, 要想確定各個組分, 須結合更詳細的生理生態(tài)方法, 具體有待進一步研究。

4 討論

利用渦度相關法和生理生態(tài)法對安吉毛竹林非生長季(2014年10月—2015年4月, 無2月)的CO2通量進行測定, 總體上, 渦度相關法與生理生態(tài)法(光合)響應關系較好。

4.1 日最大值比較

設定用渦度相關技術觀測得的凈初級生產力為1, 基于光合測定的為2。通過2.2和2.3知道兩種方法測定毛竹林非生長季各月白天的日變化均呈現(xiàn)倒“U”型變化, 中午有最大值, 但由于有飽和光強,2達到飽和光強后, 會穩(wěn)定在最大值附近波動, 直到光強低于飽和值。表2比較了兩種方法測定的各月日變化的最大值?？芍? 各月日2最大值均低于1, 即生理生態(tài)法低估了生態(tài)系統(tǒng)日變化的最大值。分析原因,1是利用通量塔進行測定, 能夠在一個較大的時間和空間尺度上估計系統(tǒng)的總生產力, 除了受到系統(tǒng)邊緣農田作物的CO2吸收的影響外, 在系統(tǒng)內部測定的組分較多, 也比較全面, 克服了通過測定室測定后累加所導致的不確定性[18]; 而2是基于毛竹葉片的光合測定的, 測定組分較為狹窄, 不包括除了毛竹葉片以外的樹干、枝條等進行光合作用吸收CO2的部分; 而且雖然此生態(tài)系統(tǒng)中近似毛竹純林, 但還是避免不了存在著其他的一些物種, 如林下草本、苔蘚類植物、蕨類等, 能夠吸收CO2進行自身的生理活動, 而這些都是2所不包括的; 此外, 儀器配置的光源為人工冷光源, 雖盡可能模仿自然光, 但是散射光對葉片光合的影響很難模擬[18]; 加上毛竹的葉面積指數(shù)利用魚眼鏡頭拍攝分析, 受光照影響大, 也很大程度地影響著2, 最終導致低估該生態(tài)系統(tǒng)2日變化的最大值。另一方面,2雖較低, 但占到1的40%以上, 在毛竹光合作用較旺盛的10月、11月和4月, 所占比例更大。表明毛竹林通過光合作用固定的CO2對生態(tài)系統(tǒng)有著較大的貢獻, 即在白天, 安吉毛竹林生態(tài)系統(tǒng)的還是主要取決于毛竹的光合作用。

表2 基于渦度相關技術和光合的日NEP最大值比較

注:1是渦度相關技術觀測得的凈初級生產力,2是基于光合測定的凈初級生產力。

4.2 月總量比較

設定用渦度相關技術觀測得的各月總凈初級生產力為1T, 白天的凈初級生產力為1D, 基于光合測定的凈初級生產力為2。表3比較了基于兩種方法測定的毛竹林非生長季(2014年10月—2015年4月, 無2月)各月的總量。

可知, 毛竹林生態(tài)系統(tǒng)非生長季的凈初級生產力總量(1T為1365.59 g·m–2, 吸收的碳約372.44 gC·m–2,整體表現(xiàn)較強的碳匯, 高于相似環(huán)境中的次生櫟林生態(tài)系統(tǒng)(為碳源)[39], 也高于環(huán)境迥異的長白山紅松針闊葉混交林(–127 gC·m–2, 為碳源)[28]。其中10月最高, 達到330.35g·m–2, 占總量的24.19%, 而1月最低, 僅94.38g·m–2, 占總量的6.91%, 3、4月又高于1月, 平均為269.88 g·m–2, 從而表現(xiàn)出明顯的季節(jié)變化, 即秋末最大, 早春次之, 冬季最低。

就白天而言,1D總量為2237.40 g·m–2,2為1900.33 g·m–2, 并占到1D的84.93%, 比例較高。各月2均低于1D, 主要是因為2測定組分單一, 不包括除毛竹葉片以外的其他能夠吸收消耗CO2的部分, 王文杰等[18]的研究也表明, 需要考慮到林下植物的光合能力, 才能夠是渦度協(xié)方差法計算的結果與生理生態(tài)學方法測定結果吻合。而且測量時不能模擬散射光, 以及葉面積指數(shù)偏低也會導致2被低估; 但是, 各月2所占比例偏高, 尤其是毛竹光合能力比較強的10、11和4月, 占到96%以上, 甚至4月高達98.06%, 分析原因, 首先, 利用渦度相關技術觀測CO2通量, 雖然比較全面, 受環(huán)境干擾小, 時間分辨率也較高, 但是地勢不平坦、平流干擾都會影響其觀測精度[30]。其次, 通量數(shù)據(jù)的處理, 包括一些極端天氣數(shù)據(jù)的收集與整理, 會影響最終得到的1D。第三, 葉片光響應曲線的測定均是在毛竹生長良好的晴天進行的, 且采用的是封閉葉室, 受外界環(huán)境影響較小, 沒有陰雨天的干擾。最后, 擬合光響應曲線的時候, 常常預測值略高于實測值, 甚至實測值會有“光抑制”, 從而使得2被高估。

表3 基于渦度相關技術和光合的NEP總量的比較

注:1T是渦度相關技術觀測得的凈初級生產力總量,1D是基于渦度相關技術測定的白天的凈初級生產力,2是基于光合測定的凈初級生產力。

另外, 從表3還可以看出, 非生長季各月的1T均低于其1D和2, 主要是因為在非生長季, 植被、凋落物及土壤微生物的呼吸作用較為旺盛, 甚至森林主要的CO2交換由光合吸收轉變?yōu)楹粑尫? 冬季最為明顯, 從而對毛竹林生態(tài)系統(tǒng)凈CO2的吸收具有抑制作用。11月份受陰雨天氣影響, 土壤濕度大, 其土壤呼吸更強, 這在Michael等[40]以及王春林等[27]的研究中也得到證實。總之, 毛竹林生態(tài)系統(tǒng)的CO2通量變化受系統(tǒng)內多個組分和過程的影響, 而不僅僅是毛竹的光合作用, 要明確毛竹林生態(tài)系統(tǒng)凈交換量的生理生態(tài)學基礎, 還需要對系統(tǒng)內多個組分進行分別觀測。

5 結論

利用渦度相關法和生理生態(tài)法對安吉毛竹林非生長季(2014年10月—2015年4月, 無2月)的CO2通量進行測定, 可知, 渦度相關法與生理生態(tài)法(光合)響應關系較好, 均具有明顯的季節(jié)變化: 秋季最高, 春季次之, 冬季最低。兩種方法測定的毛竹林生態(tài)系統(tǒng)凈初級生產力()日變化均呈倒“U”型曲線變化, 受值的直接影響, 在午間12:00前后有最大值, 10月最高。但由于光合作用有飽和光強, 故基于光合的2達到飽和光強后, 會穩(wěn)定在最大值附近波動, 直到光強低于飽和值。

就各月而言, 生理法測定的2日變化最大值均低于渦度相關技術測定的1, 但均占到1的40%以上;2月總量也均低于渦度相關技術白天測定的1D, 但2所占比例均在65%以上, 尤其是毛竹光合能力比較強的10、11和4月, 所占比例更大。表明毛竹林通過光合作用固定的CO2對生態(tài)系統(tǒng)有著較大的貢獻, 即在白天, 安吉毛竹林生態(tài)系統(tǒng)的還是主要取決于毛竹的光合作用。另外, 觀測期間, 毛竹林仍有一定的碳吸收能力, 該生態(tài)系統(tǒng)的凈初級生產力總量(1T)為1365.59 gCO2·m–2, 吸收的碳約372.44 gC·m–2, 整體表現(xiàn)為較強的碳匯(碳吸收)。

[1] 方精云. 全球生態(tài)學: 氣候變化的生態(tài)響應[M]. 北京: 高等教育出版社, 2000: 56.

[2] GIFFORD D M. The global carbon cycle: a viewpoint on the missing sink[J]. Australian Journal of Plant Physiology, 21(1): 1–15.

[3] 王效科, 馮宗煒. 森林生態(tài)系統(tǒng)生物碳儲存量的研究歷史[M]//王如松. 現(xiàn)代生態(tài)學的熱點問題研究. 北京: 中國科學技術出版社, 1996: 335–347.

[4] 崔崧. 基于微氣象學方法的落葉松人工林CO2通量研究[D]. 哈爾濱: 東北林業(yè)大學, 2007.

[5] BONAN G B. Forests and climate change: forcings, feedbacks, and the climate benefits of forests[J]. Science, 2008, 320(5882): 1444–1449.

[6] O'GRANDY A P, TISSUE D T, BEADLE C L. Canopy processes in a changing climate[J]. Tree Physiology, 2011, 31(9): 887–892.

[7] FANG Jingyun, CHEN Anping, PENG Changhui, et al. Changes in forest biomass carbon storage in China between 1949 and 1998[J]. Science, 2001, 292(5525): 2320–2322.

[8] PIAO Shilong, FANG Jingyun, CIAIS P, et al. The carbon balance of terrestrial ecosystems in China[J]. Nature, 2009, 458 (7241): 1009–1013.

[9] 孫成, 江洪, 周國模, 等. 我國亞熱帶毛竹林CO2通量的變異特征[J]. 應用生態(tài)學報, 2013, 24(10): 2717–2724.

[10] 李睿, 維爾格, 鐘章成. 施肥對毛竹竹筍生長的影響[J]. 植物生態(tài)學報, 1997, 21(1): 19–26.

[11] 馬少杰, 李正才, 王剛, 等. 集約和粗放經(jīng)營下毛竹林土壤活性有機碳的變化[J]. 植物生態(tài)學報, 2011, 35(5): 551–557.

[12] 劉玉莉, 江洪, 陳健, 等. 安吉毛竹林碳水通量及水分利用效率的日動態(tài)研究[J]. 生態(tài)科學, 2015, 34(4): 43–51.

[13] 楊爽. 浙江安吉毛竹林生態(tài)系統(tǒng)CO2通量觀測研究[D]. 臨安: 浙江農林大學, 2012.

[14] VERMA S B, BALDOCCHI D D, ANDERSON D E, et al. Eddy fluxes of CO2, water vapor, and sensible heat over a deciduous forest[J]. Boundary Layer Meteorology, 1986, 36: 71–91.

[15] BALDOCCARDO D, VALENTINI R, RUNNING S, et al. Strategies for measuring and modeling carbon dioxide and water vapor fluxes over terrestrial ecosystem [J]. Global Change Biology, 1996, 2(3): 159–168.

[16] LEE X H. On micrometeorological observation of surface- air exchange over tall vegetation [J]. Agricultural and Forest Meteorology, 1998, 91: 39–49.

[17] FALGE E, BALDOCCHI D, OLSON R, et al. Gap filling strategies for defensible annual sums of net ecosystem exchange[J]. Agricultural and Forest Meteorology, 2001, 107(1): 43–69.

[18] 王文杰, 祖元剛, 王輝民, 等. 基于渦度協(xié)方差法和生理生態(tài)法對落葉松CO2通量的初步研究[J]. 植物生態(tài)學報, 2007, 31(1): 118–128.

[19] 張憲洲, 石培禮, 劉允芬, 等. 青藏高原高寒草原生態(tài)系統(tǒng)土壤CO2排放及其碳平衡[J]. 中國科學D輯地球科學, 2004, 34(增刊II): 193–199.

[20] 王玉輝, 井長青, 白潔, 等. 亞洲東部干旱區(qū)3個典型生態(tài)系統(tǒng)生長季水碳通量特征[J]. 植物生態(tài)學報, 2014, 38(8): 795–808.

[21] 朱詠莉, 童成立, 吳金水, 等. 亞熱帶稻田生態(tài)系統(tǒng)CO2通量的季節(jié)變化特征[J]. 環(huán)境科學, 2007, 28(2): 283–288.

[22] 周麗艷, 賈丙瑞, 周廣勝, 等. 中國北方針葉林生長季碳交換及其調控機制[J]. 應用生態(tài)學報, 2010, 21(10): 2449–2456.

[23] 劉允芬, 宋霞, 孫曉敏, 等. 千煙洲人工針葉林CO2通量季節(jié)變化及其環(huán)境因子的影響[J]. 中國科學D輯地球科學, 2004, 34(增刊II): 109–117.

[24] 邱嶺, 祖元剛, 王文杰, 等. 帽兒山地區(qū)落葉松人工林CO2通量特征及對林分碳收支的影響[J]. 應用生態(tài)學學報, 2011, 22(1): 1–8.

[25] 周麗艷, 曹丙瑞, 曾偉, 等. 原始興安落葉松林生長季凈生態(tài)系統(tǒng)CO2交換及其光響應特征[J]. 生態(tài)學報, 2010, 30(24): 6919–6926.

[26] 王春林, 周國逸, 王旭, 等. 鼎湖山針闊葉混交林冠層下方CO2通量及其環(huán)境響應[J]. 生態(tài)學報, 2007, 27(3): 846–854.

[27] 王春林, 于貴瑞, 周國逸, 等. 鼎湖山常綠針闊葉混交林CO2通量估算[J]. 中國科學D輯地球科學, 2006, 36(增刊I): 119–129.

[28] 吳家兵, 關德新, 施婷婷, 等. 非生長季長白山紅松針闊葉混交林CO2通量特征[J]. 林業(yè)科學, 2006, 4(9): 1–6.

[29] 沈艷, 繆啟龍, 劉允芬. 亞熱帶紅壤丘陵人工混交林區(qū)CO2源匯及變化[J]. 生態(tài)學報, 2005, 25(6): 1371–1375.

[30] 于成龍, 劉丹. 小興安嶺天然闊葉混交林生長季CO2通量特征分析[J]. 中國農業(yè)氣象, 2011, 32(4): 525–529.

[31] 焦振, 王傳寬, 王興昌. 溫帶落葉闊葉林冠層CO2濃度的時空變異[J]. 植物生態(tài)學報, 2011, 35(5): 512–522.

[32] 趙曉松, 關德新, 吳家兵, 等. 長白山闊葉紅松林CO2通量與溫度的關系[J]. 生態(tài)學報, 2006, 26(4): 1088–1095.

[33] 趙雙菊, 張一平, 于貴瑞, 等. 西雙版納熱帶季雨林晴天CO2交換的日變化和季節(jié)變化特征[J]. 植物生態(tài)學報, 2006, 30(2): 295–301.

[34] 陳云飛, 江洪, 周國模, 等. 人工高效經(jīng)營雷竹林CO2通量估算及季節(jié)變化特征[J]. 生態(tài)學報, 2013, 33(11): 3434–3444.

[35] 劉允芬, 于貴瑞, 溫學發(fā), 等. 千煙洲中亞熱帶人工林生態(tài)系統(tǒng)CO2通量的季節(jié)變異特征[J]. 中國科學D輯地球科學, 2006, 36(增刊I): 91–102.

[36] HOLLINGER D Y, KELLIHER F M, BYERS J N, et al. Carbon dioxide exchange between an undisturbed old- growth temperate forest and the atmosphere[J]. Ecology, 1994, 75: 134–150.

[37] LOESCHER H W, OBERBAUER S F, GHOLZ H L, et al. Environmental controls on net ecosystem-level carbon exchange and productivity in a Central American tropical wet forest[J]. Global Change Biology, 2003, 9: 396–412.

[38] 藺恩杰, 江洪, 陳云飛. 太湖源雷竹林水汽通量變化及其對凈輻射的響應[J]. 浙江農林大學學報, 2013, 30(3): 313–318

[39] 劉乙, 胡海波, 劉準橋. 北亞熱帶次生櫟林生態(tài)系統(tǒng)非生長季CO2通量特征[J]. 東北林業(yè)大學學報, 2013, 41(7): 22–27.

[40] MICHAEL L G, BRUCE C D, FAN S M, et al. Physiological responses of a black spruce forest to weather[J]. Journal of geophysical research, 1997, 102: 28987–28996.

Study on the response relationship between CO2flux and photosynthesis offorest during non-growing season in Anji, Zhejiang

SHU Haiyan1, JIANG Hong1,2,3*, CHEN Xiaofeng2, SUN Wenwen2, MA Jinli2, ZHANG Minxia2

1. Key Laboratory of Eco-environments of Three Gorges Reservoir Region, Ministry of Education, Southwest University, Chongqing 400715, China 2. College of Forestry & Bio-technology, Zhejiang A & F University, Lin’an, Zhejiang 311300, China 3. International Institute for Earth System Science, Nanjing University, Nanjing 210046, China

CO2fluxes offorest ecosystem during non-growing season (from October in 2014 to April in 2015, but no February) in Anji, Zhejiang Province were measured by eddy covariance micrometeorological and eco-physiological method (photosynthesis of leaves). Both the results showed a good relationship. Diurnal variations of the net ecosystem exchange () based on eddy covariance micrometeorological method and the net ecosystem productivity () based on the eco-physiological method were single peak curves, and the maximum values were at around 12:00 noon, which were influenced bydirectly. Moreover,andwere the highest in autumn (October was the highest) and the lowest in winter. For months, the daily maximums and monthly totals ofbased on photosynthesis method were lower than those measured by the eddy covariance technique during the day, but each accounted to over 40% and 65% out of the latter respectively. Especially in October, November and April, the proportions were greater because of the high photosynthesis of. What is more, during the non-growing season, the totalofforest ecosystem was 1365.59 gCO2·m–2, and the absorption of carbon was about 372.44 gC·m–2. That was, the ecosystem was a moderate carbon sink outside the growing season.

forest ecosystem; non-growing season; eddy covariance technique; photosynthesis; net ecosystem productivity

10.14108/j.cnki.1008-8873.2017.02.003

Q945.79 , S718.5

A

1008-8873(2017)02-015-11

2015-09-02;

2015-12-03

國家自然科學基金(61190114和41171324); 上海市科學委員會項目(12231205101); 科技部國家科技基礎條件平臺項目“地球系統(tǒng)科學數(shù)據(jù)共享網(wǎng)”(2005DKA32300)

舒海燕(1990—), 女, 四川廣元人, 碩士生, 主要從事植物生理生態(tài)學研究, E-mail: shu001297@163.com

江洪, 男, 博士, 教授, 主要從事全球變化等方向研究, E-mail: jianghong_china@hotmail.com

舒海燕,江洪,陳曉峰, 等. 安吉毛竹林非生長季CO2通量與其光合響應關系探究[J]. 生態(tài)科學, 2017, 36(2): 15-25.

SHU Haiyan, JIANG Hong, CHEN Xiaofeng, et al. Study on the response relationship between CO2flux and photosynthesis offorest during non-growing season in Anji, Zhejiang[J]. Ecological Science, 2017, 36(2): 15-25.