孫文娟，張晨希，劉曉萌，劉超，陳健，胡正華*

（1. 中國科學院植物研究所，植被與環(huán)境變化國家重點實驗室，北京 100093；2. 中國科學院大學，北京 100049；3. 南京信息工程大學應用氣象學院，江蘇 南京 210044）

大氣CO2濃度升高是全球重要的環(huán)境問題。特別是近10年來，大氣CO2濃度增速加快，2015年大氣CO2濃度為400 μmol/mol，約為工業(yè)革命前的1.44倍[1]。在溫室氣體排放中等穩(wěn)定化情景下，估計2100年大氣CO2濃度將達到538～670 μmol/mol[2]。研究顯示：CO2濃度升高可促進光合作用，從而提高作物生產(chǎn)力[3-4]，其“施肥效應”能在一定程度上補償由于氣候變化導致的糧食減產(chǎn)[5]。

對全球文獻數(shù)據(jù)的統(tǒng)計分析顯示，在CO2濃度升高300 μmol/mol條件下，水稻和小麥葉片光合速率分別提高47%±5%和67%±9%，生物量分別增加 36%±2% 和 34%±2%（http://www.co2science.org/index.php）。Kimball等[6]對全球 FACE（Free Air CO2Enrichment）試驗結(jié)果進行了總結(jié)，發(fā)現(xiàn)在不同的氮肥水平和水分條件下，小麥地上生物量和籽粒產(chǎn)量分別增加3%～17%和5%～25%；氮肥供應充足條件下，水稻地上生物量和籽粒產(chǎn)量分別增加8%～17%和3%～14%。我國的相關研究始于上世紀90年代。王春乙等[7]利用自行設計的OTC-1型開頂式氣室，對冬小麥、棉花、玉米和大豆開展了連續(xù)3年的CO2濃度倍增（350 μmol/mol vs. 700 μmol/mol）試驗，結(jié)果表明CO2濃度倍增使得作物發(fā)育進程加快，株高增加，經(jīng)濟產(chǎn)量和生物產(chǎn)量增長明顯，且C3作物的增長幅度大于C4作物。Cai等[8]在FACE條件下對稻麥作物研究發(fā)現(xiàn)，CO2濃度增加到500 μmol/mol，可顯著增加作物群體結(jié)實率和產(chǎn)量，但穗數(shù)和籽粒單重未增加。對雜交和常規(guī)秈稻的研究結(jié)果也顯示，F(xiàn)ACE處理均使得其產(chǎn)量較對照有顯著增加[9-11]。王斌等[12]采用開頂箱設計，發(fā)現(xiàn)CO2濃度升高60 μmol/mol均可促進早稻和晚稻的營養(yǎng)生長和產(chǎn)量形成。采用室內(nèi)氣體熏蒸平臺的研究顯示[13-14]，CO2濃度升高200 μmol/mol可使得水稻拔節(jié)、抽穗和灌漿期凈光合速率增加，控制實驗CO2濃度增加從 60 μmol/mol到 200 μmol/mol不等[8,15]，但相同試驗條件下不同CO2濃度水平對作物生長和產(chǎn)量影響的報道還不多見。

水稻在我國糧食生產(chǎn)中占據(jù)重要的位置，明確CO2濃度升高對我國水稻產(chǎn)量的影響，是預測和評估未來全球變化背景下糧食產(chǎn)量的基礎。水稻籽粒產(chǎn)量的大部分來自開花后葉片光合產(chǎn)物，籽粒的灌漿特性對產(chǎn)量和品質(zhì)均具有重要的影響，而劍葉光合產(chǎn)物是籽粒灌漿物質(zhì)的重要來源[16-17]。CO2濃度升高，可提高水稻葉片的光合速率，加速葉片的衰老，從而影響籽粒灌漿和產(chǎn)量形成[18]。不同CO2濃度升高水平下，劍葉葉綠素含量對灌漿速率的影響可能存在差異。因此，本文試圖通過對兩個CO2濃度升高水平下水稻穗重變化和劍葉葉綠素含量的分析，探討CO2濃度升高水平對水稻灌漿速率的影響，為未來CO2濃度升高背景下估算作物產(chǎn)量提供一定的數(shù)據(jù)資料，并為農(nóng)業(yè)有效應對氣候變化提供實驗依據(jù)。

1 材料與方法

1.1 實驗設計

田間試驗在南京信息工程大學農(nóng)業(yè)氣象與生態(tài)試驗站（32.16°N，118.86°E）進行，該試驗站位于南京市浦口區(qū)盤城鎮(zhèn)落橋村，屬于亞熱帶季風氣候區(qū)，年平均降水量為1 100 mm，多年平均溫度為15.6 oC，平均日照時數(shù)超過1 900 h，無霜期為237 d。供試土壤為潴育型水稻土亞類，灰馬肝土種，耕層土壤質(zhì)地為壤質(zhì)黏土，黏粒含量為26.1%，0～20 cm土壤容重為1.57 g/cm3，有機碳、全氮含量分別為11.95 g/kg和1.19 g/kg，pH（H2O）值為6.3。

試驗以背景大氣CO2濃度為對照（CK），設置低CO2濃度升高，即在背景大氣CO2濃度基礎上增加 40 μmol/mol（CK40），高 CO2濃度升高，即在背景大氣CO2濃度基礎上增加200 μmol/mol（CK200）兩個水平處理，每個處理4個重復。采用開頂式氣室（OTC）模擬大氣CO2濃度升高，共12個OTC按隨機區(qū)組分布。OTC外形呈正八邊形棱柱狀，高3 m，對邊直徑3.75 m，底面積約為10 m2。側(cè)面為鋁合金框架，高透光性的普通玻璃。為聚攏氣體并使氣體散失速度放緩，OTC頂部開口向內(nèi)斜45°。OTC內(nèi)CO2濃度采用自動控制系統(tǒng)進行快速反饋調(diào)節(jié)，通過接收氣室內(nèi)CO2傳感器信號，根據(jù)實驗設定濃度范圍閾值，每天24 h連續(xù)供氣，實時向氣室補充CO2氣體，使得OTC內(nèi)CO2氣體濃度達到目標濃度。CO2氣源為鋼瓶氣（高壓液態(tài)CO2，純度為99%），CO2傳感器采用芬蘭Vaisala公司的GMM222 傳感器（Vaisala Inc.，Helsinki，F(xiàn)inland），量程 0～2 000 μmol/mol，精度經(jīng)標定 1 000 μmol/mol范圍以內(nèi)可達±20 μmol/mol，響應時間為30 s。水稻生長季內(nèi)CO2濃度增加時段為移栽至成熟期。

實驗選取的水稻品種為南粳9108，屬常規(guī)粳稻，于2016年5月20日播種（大田育秧），6月20日移栽，栽插密度為30穴/m2。生長季內(nèi)分三個時期施肥，總施肥量176 kgN/hm2，各時期氮肥施用的百分比為基肥∶蘗肥∶穗肥=40%:30%:30%。基肥采用復合肥（N∶P2O5∶K2O=15%:15%:15%），分蘗肥和穗肥均為尿素，水分和其他田間管理措施統(tǒng)一按當?shù)爻Ｒ?guī)栽培要求執(zhí)行。

1.2 取樣和測定方法

在水稻抽穗開花期，于各OTC中選擇同日始穗的單莖至少100個，掛牌并標注日期[19]。將中部穎花盛花期定為全穗的開花期，開花后第10 d開始取樣，前中期（9/5—9/25）每5 d取樣一次；后期（9/26—10/21）每隔8～10 d取樣一次。每個處理隨機從掛牌的穗中取10穗，取樣時間在上午10: 00～12: 00。共取樣8次。取樣后的稻穗于105℃殺青30 min后70 ℃烘干至恒重。全文統(tǒng)一以10穗重為計量單位。

為避免受實驗區(qū)域所限造成的取樣頻次和取樣量的不足，采用Logistic方程對各OTC內(nèi)的觀測數(shù)據(jù)進行擬合[20]，以期更客觀的反應灌漿期內(nèi)稻穗重量的變化：

式中：Wmax代表理論終極生長量（g/10穗），W0為初始穗重（g/10穗），k為生長速率參數(shù)；DAYs為取樣時穗開花后的天數(shù)（d）。灌漿速率（GR）用單位時間穗重的增加速率表示，計算方法為：

在每次取穗的同時，選取同穗的劍葉葉片，用手持式葉綠素測量儀（SPAD-502Plus，Konica Minolta Ltd.，Japan）測定其葉綠素含量（全文葉綠素含量用SPAD讀數(shù)來表達）。每葉片從葉尖至葉柄均勻記錄10個數(shù)據(jù)，取平均值作為該葉的葉綠素含量。

1.3 統(tǒng)計分析

數(shù)據(jù)采用Microsoft Excel 2010進行處理，統(tǒng)計分析用SYSTAT 10.0統(tǒng)計軟件對不同處理間灌漿速率的差異進行ANOVA分析，采用LSD法進行差異顯著性檢驗。

2 結(jié)果與分析

2.1 不同處理下水稻穗重的變化

對照和2個CO2濃度處理下，穗重均隨時間延長呈增加趨勢，且開花后前3周增加最快，并于開花后6周趨于平緩（圖1）。灌漿初期每10株穗重為 15.7±2.4 g，至成熟期達到 44.2±1.5 g，增加約 3倍。灌漿前期（開花后2周左右），CK200處理的穗重明顯高于CK和CK40處理（P＜0.05）；開花后20 d至水稻成熟，各處理與對照均無顯著差異（P＞0.05）。與開花后10 d的穗重相比，成熟期CK的最大穗重增加最多，約為開花后10 d的3.3倍，其次為CK40處理（3.0倍），而CK200處理增加最少（2.4倍）。

對不同處理下水稻開花后天數(shù)與穗重分別用Logistic方程（方程1）進行擬合，所得方程參數(shù)見表1。各處理的Wmax和k值在不同處理間均沒有顯著差異，CO2濃度升高對稻穗理論最大生長量和生長速率影響不大。CK200的開花后初始穗重W0顯著高于CK40，但與CK無顯著差異，說明與低濃度CO2升高處理相比，高濃度CO2處理增加了初始稻穗的重量。

圖1 水稻開花后穗重變化Fig. 1 Changes in panicle weight after anthesis. Bars are mean±SE, the same below

表1 不同處理下的水稻穗重Logistic方程參數(shù)估計Table 1 Estimates of parameters for Logistic equation of panicle weight in rice under different treatments

將上述參數(shù)帶入方程（1），計算得到逐日穗重，從而模擬水稻開花至成熟期穗重的變化動態(tài)。將各處理4個OTC內(nèi)的穗重模擬值取平均，并與觀測平均值進行對比。可以看出，采用Logistic方程能很好的模擬CK、CK40、CK200處理下水稻灌漿期穗重的變化，模擬值與觀測值具有良好的一致性（圖2）。各處理模擬穗重最大值均出現(xiàn)在開花后6～7周左右，隨后趨于平緩，與觀測結(jié)果相吻合。

2.2 不同處理下水稻灌漿速率的變化

由圖2計算得到的逐日穗重代入方程（2），求得逐日灌漿速率。結(jié)果顯示，CK、CK40和CK200處理下灌漿速率的最小值分別為0.14±0.02、0.09±0.01和 0.08±0.01 g/(10穗 .d)；最大值分別為1.20±0.04、1.28±0.03 和 1.22±0.04 g/(10 穗 .d)，分別是最小值的8.6、14.2和15.3倍。進一步分析表明，CK200處理的最大灌漿速率在抽穗后12～14 d，比CK和CK40處理提前了2～3 d。參照沈直等[21]對灌漿活躍期的定義（W達到最終穗重95%所經(jīng)歷的時間），各處理灌漿活躍期為開花至花后35～38d，不同處理間活躍灌漿期的平均灌漿速率無顯著差異（P＞0.05）。將各處理下由方程（2）計算得到的灌漿速率與實測值進行比較（圖3），可以看出，模擬值和觀測值之間具有較好的一致性（R2=0.75，P＜0.001）。

圖2 水稻開花后穗重觀測值與模擬值的比較Fig. 2 Comparison between observed and estimated panicle weights in rice after anthesis

圖3 灌漿速率觀測值與模擬值的比較Fig. 3 Comparison between observed and estimated grain filling rates in rice

2.3 劍葉葉綠素含量對灌漿速率的影響

對各處理不同測定階段植株劍葉葉綠素含量的分析表明，葉綠素含量自水稻開花后約2周左右達到最大值（SPAD值接近50），其后逐漸降低，至成熟期下降到39.3，最大值為最小值的1.3倍（圖4）。開花后6周內(nèi)各處理間劍葉葉綠素含量均無顯著差異（P＞0.05），而6周后CK葉綠素含量下降較CK40和CK200更快，6周后葉綠素含量平均值表現(xiàn)為CK200＞CK40＞CK。表明在CO2濃度升高200 μmol/mol下，水稻臨近成熟期其劍葉葉綠素含量較高，其光合作用的能力仍高于CK和CK40處理。將各處理的葉綠素含量與稻穗灌漿速率相比較發(fā)現(xiàn)，開花后葉綠素含量的變化趨勢與稻穗灌漿速率的變化具有明顯的一致性，均表現(xiàn)為開花后2周左右達到最大值，隨后降低，至成熟期最低（圖4）。

圖4 開花后灌漿速率與劍葉葉綠素含量動態(tài)Fig. 4 Dynamics of grain filling rate of flag leaves and chlorophyll content (SPAD reading) after anthesis

進一步分析表明，不同處理下的灌漿速率隨劍葉葉綠素含量的增加而增加，二者呈現(xiàn)極顯著線性正相關關系（P＜0.001），說明劍葉葉綠素含量是決定灌漿速率的主要因素（圖5）。CK40和CK200處理下，單位葉綠素含量對灌漿速率的貢獻較為接近，分別為0.136和0.131，且均高于對照（k=0.093)，表明CO2濃度升高可提高劍葉葉綠素對水稻籽粒灌漿的貢獻。

圖5 灌漿速率與劍葉葉綠素含量的關系Fig. 5 Relationship between grain filling rate and chlorophyll content (SPAD reading) of flag leaves

3 討論

本研究采用實驗觀測和經(jīng)驗模型相結(jié)合的方法，模擬灌漿期內(nèi)水稻穗重的變化，并由此計算灌漿速率。由于開頂箱實驗的小區(qū)面積較小（10 m2），受實驗區(qū)域的限制，用于研究灌漿速率的稻穗每次取10穗，共取8次。這與其他同類研究的取樣量和取樣頻次[22-23]相比略少。同樣由于取樣量所限，個別處理的穗重出現(xiàn)隨時間增加反而降低的現(xiàn)象，即灌漿速率表現(xiàn)為負值（如圖2的CK40處理最后兩次穗重），這顯然有悖于實際情況。因此，采用作物生長方程擬合的方法，能更客觀的反映灌漿期內(nèi)穗重的變化。趙黎明等[22]和李俊周等[24]也采用生長方程的方法對籽粒灌漿過程的籽粒增重觀測數(shù)據(jù)進行擬合，分析最終生長量，并計算灌漿速率動態(tài)。本研究采用Logistic生長方程對每個處理的觀測值進行擬合，并確定方程參數(shù)（表1），可得到逐日穗重的變化量。結(jié)果顯示方程的模擬效果較好，與觀測結(jié)果的變化趨勢吻合（圖2）。

李軍營等[25]利用我國第一個FACE平臺，對武香粳14的研究結(jié)果顯示，CO2濃度升高200 μmol/mol可顯著促進其灌漿速率，并使灌漿速率提前3 d到最大值。本研究也發(fā)現(xiàn)，CO2濃度增加200 μmol/mol處理下，水稻灌漿前期穗重高于CK40處理和CK（圖1），且其最大灌漿速率比CK和CK40處理提前了2～3 d。但值得注意的是，與CK200不同，低CO2濃度升高處理（CK40）對籽粒灌漿速率和灌漿進程的影響不顯著。全球大部分針對大氣CO2濃度升高對作物生理生態(tài)和產(chǎn)量影響的控制實驗研究都是建立在一個CO2濃度增加較高的水平（如，比背景大氣增加200 μmol/mol）下，但自然狀態(tài)下大氣CO2濃度升高是一個緩慢的過程，升高200 μmol/mol需要若干年[2]。因此，低CO2濃度升高下的研究結(jié)果，對理解未來大氣CO2濃度升高的生態(tài)學效應具有更重要的借鑒意義。

Klironomos等[26]以燕麥草（B. inermis）為對象，研究了CO2濃度增加200 μmol/mol和每隔15周增加 10 μmol/mol逐漸增加至 200 μmol/mol的影響，發(fā)現(xiàn)CO2驟然增加和逐漸增加200 μmol/mol，對雀麥草生長的影響不同。高濃度CO2驟然增加的影響要高于低濃度逐漸增加的影響。本研究結(jié)果顯示，低CO2濃度增加對水稻籽粒灌漿進程的影響與背景大氣無顯著性差異，但低于高濃度CO2增加的影響，與Klironomos等[26]的結(jié)果具有一致性。由于本研究的試驗只進行了一年，低CO2濃度增加是否存在多年累積效應，以及籽粒灌漿對長期高CO2濃度增加是否存在適應性還需進一步開展研究，以便為更客觀評估未來CO2濃度升高的生態(tài)效應提供參考。

水稻劍葉光合作用對產(chǎn)量的貢獻最大[16]，而開花后至成熟期劍葉葉綠素含量的變化反映了水稻光合能力的變化。籽粒灌漿啟動后劍葉開始逐漸衰老，葉綠素含量降低[16,27]，本研究顯示了類似的趨勢（圖4）。Goufo等[28]采用開頂箱對水稻的研究發(fā)現(xiàn)，CO2濃度升高175 μmol/mol，水稻葉片的葉綠素含量與對照相比略有下降。而胡健等[29]發(fā)現(xiàn)，不同施氮水平下，高濃度CO2處理下的劍葉葉綠素含量在抽穗后5～10 d均有所增加，但在抽穗后15～25 d下降。王惠貞等[30]研究了CO2濃度600 μmol/mol下水稻劍葉葉綠素含量，發(fā)現(xiàn)抽穗期和乳熟期劍葉總?cè)~綠素含量較背景大氣下有顯著增加。

本研究結(jié)果顯示，高低兩個CO2濃度升高水平下，劍葉葉綠素含量與對照相比均無顯著差異（P＞0.05），但CO2濃度升高顯著提高了單位葉綠素含量對灌漿速率的貢獻。這可能與CO2濃度升高促進了葉片衰老[31]，從而加速其灌漿進程[18]有關。歐陽杰等[32]對5個秈稻品種研究發(fā)現(xiàn)，水稻灌漿中后期，劍葉葉綠素含量與單穗粒重關系密切，且抽穗后15～20 d，劍葉中葉綠素降解速率與單穗粒重顯著負相關，這與本研究葉綠素含量對灌漿速率影響的結(jié)果有一致性。

4 結(jié)論

高CO2濃度升高下，水稻灌漿期前期穗重高于低CO2濃度升高和背景大氣，其后至成熟期，CO2濃度升高對其無顯著影響；高CO2濃度升高下，水稻最大灌漿速率出現(xiàn)時間較低CO2濃度升高和背景大氣下有所提前，而活躍灌漿期的平均灌漿速率無顯著差異。植株劍葉葉綠素含量影響灌漿速率，葉綠素含量越高，灌漿速率越大，CO2濃度升高下，葉綠素含量對灌漿速率的貢獻高于背景大氣。

致謝：感謝江蘇省農(nóng)業(yè)科學院王才林研究員提供水稻種子。

[1] World Meteorological Organization. The state of greenhouse gases in the atmosphere based on global observations through 2013[M]//WMO Greenhouse Gas Bulletin, Switzerland, 2014.

[2] IPCC. Summary for Policymakers[M]//Stocker T F, Qin D,Plattner G K. Climate Change 2013: The Physical Science Basis.Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge University Press, Cambridge, United Kingdom and New York,NY, USA, 2013.

[3] Norby R J, Zak D R. Ecological lessons from free-air CO2enrichment (FACE) experiments[J]. Annual Review of Ecology,Evolution, and Systematics, 2011, 42: 181-203.

[4] Rosenzweig C, Elliott J, Deryng D, et al. Assessing agricultural risks of climate change in the 21st century in a global gridded crop model intercomparison[J]. Proceedings of the National Academy of Sciences USA, 2014, 111(9): 3268-3273.

[5] Long S P. Virtual Special Issue on food security-greater than anticipated impacts of near-term global atmospheric change on rice and wheat[J]. Global Change Biology, 2012, 18: 1489-1490.

[6] Kimball B A, Morris C F, Pinter Jr P J, et al. Elevated CO2,drought and soil nitrogen effects on wheat grain quality[J]. New Phytologist, 2001, 150: 295-303.

[7] 王春乙, 潘亞茹, 白月明, 等. CO2濃度倍增對中國主要作物影響的試驗研究[J]. 氣象學報, 1997, 55(1): 86-94.Wang C Y, Pan Y R, Bai Y M, et al. The experiment study of effects doubled CO2concentration on several main crops in China[J]. Acta Meteorologica Sinica, 1997, 55(1): 86-94.

[8] Cai C, Yin X, He S, et al. Responses of wheat and rice to factorial combinations of ambient and elevated CO2and temperature in FACE experiments[J]. Global Change Biology, 2016, 22: 856-874.

[9] 黃建曄, 楊洪建, 董桂春, 等. 開放式空氣CO2濃度增高對水稻產(chǎn)量形成的影響[J]. 應用生態(tài)學報, 2002, 13(10): 1210-1214.Huang J Y, Yang H J, Dong G C, et al. Effects of free-air CO2enrichment (FACE) on yield formation in rice (Oryza sativa)[J].Chinese Journal of Applied Ecology, 2002, 13(10): 1210-1214.

[10] Yang L X, Huang J Y, Yang H J, et al. The impact of free-air CO2enrichment (FACE) and N supply on yield formation of rice crops with large panicle[J]. Field Crops Research, 2006, 98: 141-150.

[11] 劉紅江, 楊連新, 黃建曄, 等. FACE對常規(guī)秈稻揚稻6號產(chǎn)量形成的影響[J]. 農(nóng)業(yè)環(huán)境科學學報, 2009, 28(2): 299-304.Liu H J, Yang L X, Huang J Y, et al. Effect of free air CO2enrichment（FACE）on yield formation of Indica rice cultivar Yangdao 6[J]. Journal of Agro-Environment Science, 2009, 28(2):299-304.

[12] 王斌, 萬運帆, 郭晨, 等. 模擬大氣溫度和CO2濃度升高對雙季稻氮素利用的影響[J]. 作物學報, 2015, 41(8): 1295-1303.Wang B, Wan Y F, Guo C, et al. Effects of elevated air temperature and carbon dioxide concentration on nitrogen use of double rice(Oryza sativa L.) in open-top chambers[J]. Acta Agronomica Sinica, 2015, 41(8): 1295-1303.

[13] 邵在勝, 趙軼鵬, 宋琪玲, 等. 大氣CO2和O3濃度升高對水稻‘汕優(yōu)63’葉片光合作用的影響[J]. 中國生態(tài)農(nóng)業(yè)學報, 2014,22(4): 422-429.Shao Z S, Zhao Y P, Song Q L, et al. Impact of elevated atmospheric carbon dioxide and ozone concentrations on leaf photosynthesis of ‘Shanyou 63’ hybrid rice[J]. Chinese Journal of Eco-Agriculture, 2014, 22(4): 422-429.

[14] 趙軼鵬, 邵在勝, 王云霞, 等. 大氣CO2和O3濃度升高對汕優(yōu)63生長動態(tài)、物質(zhì)生產(chǎn)和氮素吸收的影響[J]. 生態(tài)學報，2015, 35(24): 8128-8138.Zhao Y P, Shao Z S, Wang Y X, et al. Impacts of elevated atmospheric carbon dioxide and ozone concentration on growth dynamic, dry matter production, and nitrogen uptake of hybrid rice Shanyou 63[J]. Acta Ecologic Sinica, 2015, 35(24): 8128-8138.

[15] Chen C, Jiang Q, Ziska L H, et al. Seed vigor of contrasting rice cultivars in response to elevated carbon dioxide[J]. Field Crops Research, 2015, 178: 63-68.

[16] Haque M M, Khanam A H M, Biswas D K, et al. The effect of elevated CO2concentration on leaf chlorophyll and nitrogen contents in rice during post-flowering phases[J]. Biologia Plantarum, 2006, 50: 69-73.

[17] 金松恒, 蔣德安, 王品美, 等. 水稻孕穗期不同葉位葉片的氣體交換與葉綠素熒光特性[J]. 中國水稻科學, 2004, 18(5): 443-448.Jin S H, Jiang D A, Wang P M, et al. Characteristics of gas exchange and chlorophyll fluorescence in different position leaves at booting stage in rice plants[J]. Chinese Journal of Rice Science,2004, 18(5): 443-448.

[18] Zhu C, Zhu J, Zeng Q, et al. Elevated CO2accelerates flag leaf senescence in wheat due to ear photosynthesis which causes greater ear nitrogen sink capacity and ear carbon sink limitation[J].Functional Plant Biology, 2009, 36(4): 291-299.

[19] 朱慶森, 曹顯祖, 駱亦琪. 水稻籽粒灌漿的生長分析[J]. 作物學報, 1988, 14(3): 182-189.Zhu Q S, Cao X Z, Luo Y Q. Growth analysis on the process of grain filling in rice[J]. Acta Agronomica Sinica, 1988, 14(3):182-189.

[20] 姜樂, 趙洪亮, 謝立勇, 等. CO2濃度增高對北方水稻籽粒灌漿特性的影響[J]. 華北農(nóng)學報, 2014, 29(3): 136-141.Jiang L, Zhao H L, Xie L Y, et al. Effect of CO2enrichment on rice grain filling under FACE system in north China[J]. Acta Agriculturae Boreali-Sinica, 2014, 29(3): 136-141.

[21] 沈直, 唐設, 張海祥, 等. 灌漿期開放式增溫對水稻強勢粒和弱勢粒淀粉代謝關鍵酶相關基因表達水平的影響[J]. 南京農(nóng)業(yè)大學學報, 2016, 39(6): 898-906.Shen Z, Tang S, Zhang H Y, et al. Effect of T-FACE high temperature on genes expression level of key enzymes involved in starch metabolism in superior spikelets and inferior spikelets of rice during grain filling period[J]. Journal of Nanjing Agricultural University, 2016, 39(6): 898-906.

[22] 趙黎明, 李明, 鄭殿峰, 等. 灌溉方式對寒地水稻產(chǎn)量及籽粒灌漿的影響[J]. 中國農(nóng)業(yè)科學, 2015, 48(22): 4493-4506.Zhao L M, Li M, Zheng D F, et al. Effects of irrigation regimes on yield and grain filling of rice (Oryza sativa L.) in cold region[J].Scientia Agricultura Sinica, 2015, 48(22): 4493-4506.

[23] 姜元華, 張洪程, 趙可, 等. 長江下游地區(qū)不同類型水稻品種產(chǎn)量及其構(gòu)成因素特征的研究[J]. 中國水稻科學, 2014, 28(6):621-631.Jiang Y H, Zhang H C, Zhao K, et al. Difference in yield and its components characteristics of different type rice cultivars in the lower reaches of the Yangtze River[J]. Chinese Journal of Rice Science, 2014, 28(6): 621-631.

[24] 李俊周, 李磊, 孫傳范, 等. 水氮互作對水稻籽粒充實及產(chǎn)量的影響[J]. 中國農(nóng)業(yè)大學學報, 2011, 16(3): 42-47.Li J Z, Li L, Sun C F, et al. Effects of water-nitrogen interaction on rice grain plumpness and yield[J]. Journal of China Agricultural University, 2011, 16(3): 42-47.

[25] 李軍營, 徐長亮, 謝輝, 等. CO2濃度升高加快水稻灌漿前期籽粒的生長發(fā)育進程[J]. 作物學報, 2006, 32(6): 905-910.Li J Y, Xu C L, Xie H, et al. Acceleration of grain growth and development process by FACE during early grain filling stage of rice (Oryza sativa L.)[J]. Acta Agronomica Sinica, 2006, 32(6):905-910.

[26] Klironomos J N, Allen M F, Rillig M C, et al. Abrupt rise in atmospheric CO2overestimates community response in a model plant–soil system[J]. Nature, 2005, 433: 621-624.

[27] Panda D, Sarkar R K. Natural leaf senescence: Probed by chlorophyll fluorescence, CO2photosynthetic rate and antioxidant enzyme activities during grain filling in different rice cultivars[J].Physiology and Molecular Biology of Plants, 2013, 19(1): 43-51.

[28] Goufo P, Pereira J, Moutinho-Pereira J, et al. Rice (Oryza sativa L.) phenolic compounds under elevated carbon dioxide (CO2)concentration[J]. Environmental and Experimental Botany, 2014,99: 28-37.

[29] 胡健, 楊連新, 周娟, 等. 開放式空氣CO2濃度增高(FACE)對水稻灌漿動態(tài)的影響[J]. 中國農(nóng)業(yè)科學, 2007, 40(11): 2443-2451.Hu J, Yang L X, Zhou J, et al. Effect of free-air CO2enrichment(FACE) on grain filling dynamics of rice[J]. Scientia Agricultura Sinica, 2007, 40(11): 2443-2451.

[30] 王惠貞, 趙洪亮, 馮永祥, 等. 北方水稻生育后期劍葉可溶性物質(zhì)含量及植株生產(chǎn)力對CO2濃度增高的響應[J]. 作物學報,2014, 40(2): 320-328.Wang H Z, Zhao H L, Feng Y X, et al. Response of soluble substances content in flag leaves during late growth stage and plant productivity of rice to elevated CO2in north China[J]. Acta Agronomica Sinica, 2014, 40(2): 320-328.

[31] Zhu C, Ziska L, Zhu J, et al. The temporal and species dynamics of photosynthetic acclimation in flag leaves of rice (Oryza sativa)and wheat (Triticum aestivum) under elevated carbon dioxide[J].Physiologia Plantarum, 2012, 145: 395-405.

[32] 歐陽杰, 王楚桃, 何光華, 等. 水稻灌漿中后期功能葉中葉綠素含量及其變化趨勢與谷物產(chǎn)量關系研究[J]. 西南農(nóng)業(yè)學報,2012, 25(4): 1201-1204.Ouyang J, Wang C T, He G H, et al. Study on relationship between different functional leaf chlorophyll content and its trends in mid and late period of rice filling and grain yield[J]. Southwest China Journal of Agricultural Sciences, 2012, 25(4): 1201-1204.