• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    1981—2010年長(zhǎng)江中下游地區(qū)單季稻生殖生長(zhǎng)期對(duì)氣候變化和技術(shù)進(jìn)步的響應(yīng)

    2023-03-23 07:07:52劉二華周廣勝武炳義宋艷玲何奇瑾呂曉敏周夢(mèng)子
    作物學(xué)報(bào) 2023年5期
    關(guān)鍵詞:長(zhǎng)江中下游地區(qū)單季稻乳熟期

    劉二華 周廣勝,,* 武炳義 宋艷玲 何奇瑾 呂曉敏 周夢(mèng)子

    1981—2010年長(zhǎng)江中下游地區(qū)單季稻生殖生長(zhǎng)期對(duì)氣候變化和技術(shù)進(jìn)步的響應(yīng)

    劉二華1,2,3周廣勝1,2,3,4,*武炳義2宋艷玲1,3何奇瑾4呂曉敏1,3周夢(mèng)子1,3

    1中國(guó)氣象科學(xué)研究院, 北京 100081;2復(fù)旦大學(xué)大氣科學(xué)研究院, 上海 200439;3中國(guó)氣象科學(xué)研究院與鄭州大學(xué)生態(tài)氣象聯(lián)合實(shí)驗(yàn)室, 河南鄭州 450001;4中國(guó)農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院, 北京 100193

    作物生殖生長(zhǎng)期長(zhǎng)度與作物產(chǎn)量和品質(zhì)密切相關(guān)。為深入探究作物生殖生長(zhǎng)期長(zhǎng)度(reproductive growth period lengths, RGLs)對(duì)氣候變化和技術(shù)進(jìn)步的響應(yīng), 基于1981—2010年長(zhǎng)江中下游地區(qū)單季稻生殖生長(zhǎng)期和氣象數(shù)據(jù), 量化不同RGLs (孕穗期—抽穗期(booting to heading, BDHD)、抽穗期—乳熟期(heading to milking, HDMS)、乳熟期—成熟期(milking to maturity, MSMD)和孕穗期—成熟期(booting to maturity, BDMD))對(duì)平均溫度(mean temperature, TEM)、累積降水量(cumulative precipitation, PRE)和累積日照時(shí)數(shù)(cumulative sunshine duration, SSD)的敏感性, 并分離氣候變化和技術(shù)進(jìn)步對(duì)不同RGLs的影響。結(jié)果表明, 1981—2010年長(zhǎng)江中下游地區(qū)單季稻BDMD呈延長(zhǎng)趨勢(shì)(0.24 d a–1), 其中, HDMS延長(zhǎng)趨勢(shì)最明顯(0.16 d a–1)。氣候因子中高溫和寡照不利于單季稻不同RGLs延長(zhǎng), 其中, TEM對(duì)BDHD、HDMS和MSMD變化趨勢(shì)的平均相對(duì)貢獻(xiàn)分別為–50.0%、–50.7%和–21.9%, SSD對(duì)BDHD、HDMS和MSMD變化趨勢(shì)的平均相對(duì)貢獻(xiàn)分別為–47.2%、–48.7%和–67.6%。技術(shù)進(jìn)步彌補(bǔ)了氣候變化對(duì)不同RGLs變化趨勢(shì)的不利影響。研究表明, 技術(shù)進(jìn)步可能是當(dāng)前單季稻穩(wěn)產(chǎn)高產(chǎn)和趨利避害的主要手段, 未來(lái)可以采用較長(zhǎng)生殖生長(zhǎng)期和耐熱性品種來(lái)適應(yīng)持續(xù)的氣候變化。

    單季稻; 氣候變化; 生殖生長(zhǎng)期; 技術(shù)進(jìn)步

    作物生育期是重要的農(nóng)業(yè)信息, 對(duì)合理安排作物布局、調(diào)控作物生長(zhǎng)發(fā)育和預(yù)測(cè)產(chǎn)量具有重要意義[1-2]。同時(shí), 作物生育期在陸面動(dòng)態(tài)過(guò)程、生物物理過(guò)程中扮演重要角色[3]。作物關(guān)鍵生育期的相對(duì)長(zhǎng)短是決定作物產(chǎn)量的關(guān)鍵因素[4]。氣候變暖普遍加快作物生長(zhǎng)發(fā)育, 縮短作物生長(zhǎng)期, 縮短作物生物量的積累時(shí)間, 威脅作物的生產(chǎn)力和糧食安全[5]。因此, 隨氣候變暖的加劇[6], 作物生育期對(duì)氣候變化的響應(yīng)已然成為研究的熱點(diǎn)[7]。

    作物生育期變化反映了作物生長(zhǎng)過(guò)程對(duì)自然因素和人為因素的綜合響應(yīng)[8-9], 自然因素主要為氣候因素(包括溫度、水分和光照等), 人為因素表現(xiàn)為技術(shù)進(jìn)步(包括品種、種植方式、栽培管理和生產(chǎn)投入等)。作物生育期的變化表現(xiàn)為提前或推遲, 生長(zhǎng)期縮短或延長(zhǎng)[10]。研究表明中國(guó)單季稻生育期呈延長(zhǎng)趨勢(shì)[11-13], 尤其以生殖生長(zhǎng)期(reproductive growth periods lengths, RGLs)顯著延長(zhǎng)為主[14]。同時(shí), 中國(guó)玉米、小麥、大豆以及棉花等作物RGLs普遍延長(zhǎng)[5,15-19]。已有大量研究通過(guò)統(tǒng)計(jì)模型和作物模擬模型等多種方法廣泛評(píng)價(jià)了作物生育期對(duì)氣候變化和技術(shù)進(jìn)步的響應(yīng)[20-22]。例如, Zhang等[14]的研究表明水稻RGLs由溫度和品種特性共同決定, 溫度升高加快水稻抽穗期, 品種改變延緩水稻抽穗期。而Hu等[23]認(rèn)為中國(guó)水稻生育期延長(zhǎng)主要由播種期和品種改變調(diào)節(jié), Bai等[24]進(jìn)一步研究表明播種期顯著影響水稻營(yíng)養(yǎng)生長(zhǎng)期, 而對(duì)RGLs影響不明顯, 品種改變使得水稻RGLs每10年延長(zhǎng)0.2~2.8 d。然而, 氣候因素(如降水量和日照時(shí)數(shù))對(duì)作物生育期的影響也至關(guān)重要, 但以往研究通常忽略這些氣候因素[19,25-26]。因此, 厘清氣候變化和技術(shù)進(jìn)步對(duì)作物生育期的影響程度對(duì)于準(zhǔn)確實(shí)施適應(yīng)氣候變化的管理措施具有重要參考價(jià)值[5,25,27-28]。由于技術(shù)的不斷進(jìn)步, 作物生育期對(duì)氣候變化的響應(yīng)難以評(píng)估[8]。當(dāng)前已有一些研究開(kāi)展了關(guān)于量化氣候變化和作物管理對(duì)作物生育期的影響研究[23,26], 但是, 鮮少有關(guān)于氣候變化和技術(shù)進(jìn)步對(duì)作物不同RGLs變化趨勢(shì)的綜合定量評(píng)價(jià)研究。

    水稻(L.)是世界上一半以上人口的主食[29], 長(zhǎng)江中下游地區(qū)是中國(guó)最古老最重要的水稻主產(chǎn)地之一[30], 水稻種植面積和總產(chǎn)量均占全國(guó)45%以上[31]。長(zhǎng)江中下游地區(qū)水稻種植制度為單季稻和早、晚雙季稻[32]。1980—2008年呈現(xiàn)出單季稻種植面積增加和雙季稻種植面積減少的態(tài)勢(shì)[33-35], 這表明單季稻在長(zhǎng)江中下游地區(qū)水稻生產(chǎn)中占有越來(lái)越重要的地位。本研究假設(shè)單季稻不同RGLs是氣候變化和技術(shù)進(jìn)步綜合作用的結(jié)果?;谝陨霞僭O(shè), 本研究的目標(biāo)是(1) 闡明長(zhǎng)江中下游地區(qū)單季稻不同RGLs時(shí)空變化特征; (2) 單季稻不同RGLs對(duì)平均溫度(mean temperature, TEM)、累積降水量(cumulative precipitation, PRE)和累積日照時(shí)數(shù)(cumulative sunshine duration, SSD)的敏感性; (3) 分離和量化氣候變化和技術(shù)進(jìn)步對(duì)單季稻RGLs變化趨勢(shì)影響的貢獻(xiàn)。

    1 材料與方法

    1.1 研究區(qū)域與資料來(lái)源

    研究區(qū)域位于中國(guó)長(zhǎng)江中下游地區(qū)(26°~35°N, 103°~122°E), 該地區(qū)屬于雨熱同期的亞熱帶季風(fēng)氣候區(qū)。本研究涉及的地區(qū)主要包括江蘇、河南南部、湖北、重慶和四川東部(圖1)。本研究所用資料主要包括氣象資料和單季稻生育期資料, 資料均來(lái)源于中國(guó)氣象局。氣象資料主要包括1981—2010年10個(gè)氣象站點(diǎn)觀測(cè)的逐日平均溫度、降水量和日照時(shí)數(shù)。根據(jù)長(zhǎng)期觀測(cè)資料缺失情況, 選取的10個(gè)農(nóng)業(yè)氣象站位于與上述氣象站相同的地點(diǎn)。單季稻生育期主要包括分蘗期、孕穗期、抽穗期、乳熟期和成熟期, 本研究將生殖生長(zhǎng)階段進(jìn)一步劃分為孕穗期—抽穗期(booting to heading, BDHD)、抽穗期—乳熟期(heading to milking, HDMS)、乳熟期—成熟期(milking to maturity, MSMD)和孕穗期—成熟期(booting to maturity, BDMD), 同時(shí)探討了分蘗期—孕穗期(tillering to booting, TDBD)的變化特征。

    圖1 10個(gè)氣象站點(diǎn)地理位置

    JXH、JZJ、HGS、HXY、HFX、HZX、SDZ、SGP、CFD和CYY分別表示興化站、鎮(zhèn)江站、固始站、信陽(yáng)站、房縣站鐘祥站、大竹站、高坪站、豐都站和酉陽(yáng)站。矢量地圖數(shù)據(jù)來(lái)源: 全國(guó)地理信息資源目錄服務(wù)系統(tǒng)(https://www.webmap.cn/)。

    JXH, JZJ, HGS, HXY, HFX, HZX, SDZ, SGP, CFD, and CYY represent Xinghua station, Zhenjiang station, Gushi station, Xinyang station, Fangxian station, Zhongxiang station, Dazhu station, Gaoping station, Fengdu station, and Youyang station, respectively. The above vector map is from the National Catalogue Service for Geographic Information (https://www.webmap.cn/).

    1.2 分析方法

    利用線性回歸方法分析單季稻關(guān)鍵生殖生長(zhǎng)期和不同RGLs (Booting、Heading、Milking、Maturity、TDBD、BDHD、HDMS、MSMD和BDMD)以及相應(yīng)時(shí)期氣候因子(TEM、PER和SSD)的變化趨勢(shì), 采用雙尾t檢驗(yàn)分析確定統(tǒng)計(jì)顯著性。

    Y=T×X+B(1)

    式(1)中,Y是第年各變量的值; T為線性回歸斜率;X是年份(= 1, 2, 3, …, 30); B是截距。

    面板回歸模型法[36]已被證明適合氣候變化影響的研究[37]。本研究利用面板回歸模型法得到單季稻RGLs對(duì)氣候因子的響應(yīng)及其敏感性。為降低氣候因子間的多重共線性, 本研究根據(jù)相關(guān)性分析并結(jié)合前人研究[25,38], 選取TEM、PRE和SSD三個(gè)氣候因子作為預(yù)測(cè)變量。

    首先, 構(gòu)建了5個(gè)單季稻RGLs (TDBD、BDHD、HDMS、MSMD和BDMD)多元面板回歸模型。模型形式如下:

    RGLt,s=β0,s+β1×TEMt,s+β2×PREt,s+β3×SSDt,s+(2)

    式(2)中, RGLt,s表示第s站點(diǎn)和第t年的生殖生長(zhǎng)期長(zhǎng)度(d), TEMt,s、PREt,s和SSDt,s分別表示第s站點(diǎn)和第t年的生殖生長(zhǎng)階段的TEM (℃)、PRE (mm)和SSD (h)。β0表示站點(diǎn)的固定效應(yīng)。β1、β2和β3分別表示RGLs對(duì)TEM、PRE和SSD的敏感性(d ℃–1、d mm–1和d h–1)。ε是一個(gè)隨機(jī)誤差項(xiàng), 表示氣候變化以外的因素對(duì)作物生育期的影響。

    其次, 將各個(gè)站點(diǎn)單季稻分蘗期、孕穗期、抽穗期、乳熟期和成熟期的日期確定為1981—2010年實(shí)際觀測(cè)日期平均值, 統(tǒng)計(jì)1981—2010年各站點(diǎn)氣候因子在不同RGLs的變化趨勢(shì), 據(jù)此分離氣候變化趨勢(shì)對(duì)單季稻RGLs變化趨勢(shì)的影響[39]。因此, 僅受氣候變化影響的單季稻RGLs變化趨勢(shì)的估算方法如下:

    TRGL_cli=β1×TTEM+β2×TPRE+β3×TSSD(3)

    式(3)中, TRGL_cli表示氣候變化引起的生殖生長(zhǎng)期變化趨勢(shì)。TTEM、TPRE和TSSD分別表示特定生殖生長(zhǎng)階段的TEM、PRE和SSD的變化趨勢(shì)(℃ a–1、mm a–1和h a–1)。TTEM、TPRE和TSSD均由公式(1)計(jì)算得到。

    第三, 從氣候變化和技術(shù)進(jìn)步的綜合影響中消除氣候變化的影響, 得到技術(shù)進(jìn)步對(duì)單季稻RGLs變化趨勢(shì)的影響:

    TRGL_tec=TRGL_cli_tec–TRGL_cli(4)

    式(4)中, TRGL_tec表示技術(shù)進(jìn)步對(duì)單季稻RGLs變化趨勢(shì)的影響(d a–1), TRGL_cli_tec表示氣候變化和技術(shù)進(jìn)步對(duì)單季稻RGLs變化趨勢(shì)的綜合影響(d a–1), 即實(shí)際觀測(cè)得到的不同RGLs變化趨勢(shì)。TRGL_cli_tec可由公式(1)計(jì)算得到。

    最后, 分離出氣候變化和技術(shù)進(jìn)步導(dǎo)致單季稻RGLs變化趨勢(shì)的貢獻(xiàn), 計(jì)算式如下:

    式中, CRGL_cli和CRGL_tec分別表示氣候變化和技術(shù)進(jìn)步對(duì)RGLs變化趨勢(shì)的貢獻(xiàn)(%), RCRGL_TEM表示TEM對(duì)單季稻RGLs變化趨勢(shì)的平均相對(duì)貢獻(xiàn)(%), 同樣,利用公式(7)計(jì)算PRE和SSD對(duì)單季稻RGLs變化趨勢(shì)的平均相對(duì)貢獻(xiàn)(RCRGL_PRE和RCRGL_SSD)。參考Li等[10]研究, 計(jì)算氣候變化和技術(shù)進(jìn)步對(duì)單季稻RGLs變化趨勢(shì)的平均貢獻(xiàn), 以及各個(gè)氣候因子對(duì)單季稻RGLs變化趨勢(shì)的平均相對(duì)貢獻(xiàn)。

    2 結(jié)果與分析

    2.1 單季稻生殖生長(zhǎng)期時(shí)空變化

    1981—2010年長(zhǎng)江中下游地區(qū)單季稻分蘗期、孕穗期、抽穗期、乳熟期和成熟期以及不同RGLs的時(shí)空變化趨勢(shì)如圖2所示。研究發(fā)現(xiàn), 分蘗期以0.37 d a–1的趨勢(shì)提前, 其中10個(gè)站點(diǎn)中9個(gè)呈提前趨勢(shì), 7個(gè)顯著提前。不同于其他站點(diǎn), 興化站分蘗期呈顯著推遲趨勢(shì)(0.39 d a–1)。孕穗期和抽穗期時(shí)空變異較大, 平均變化趨勢(shì)分別為0.08 d a–1和0.11 d a–1。乳熟期和成熟期主要呈不同程度的推遲趨勢(shì), 分別為0.27 d a–1和0.32 d a–1。以上結(jié)果表明生殖生長(zhǎng)期前生育期呈提前趨勢(shì), 隨生育進(jìn)程推進(jìn), 各個(gè)生育期逐漸推遲。

    圖2-f~j分別表示長(zhǎng)江中下游地區(qū)單季稻TDBD、BDHD、HDMS、MSMD和BDMD時(shí)空變化趨勢(shì)。分蘗期顯著提前導(dǎo)致TDBD以0.45 d a–1的速率延長(zhǎng), 其中10個(gè)站點(diǎn)中9個(gè)呈延長(zhǎng)趨勢(shì), 7個(gè)顯著延長(zhǎng), 而鐘祥站TDBD以0.34 d a–1速率縮短。單季稻BDHD變化不明顯(0.03 d a–1)。單季稻HDMS以0.16 d a–1速率延長(zhǎng), MSMD變化不明顯(0.05 d a–1), BDMD以0.24 d a–1速率延長(zhǎng)。綜合而言, 長(zhǎng)江中下游地區(qū)單季稻RGLs呈延長(zhǎng)趨勢(shì), 其中以HDMS延長(zhǎng)為主。

    2.2 單季稻生殖生長(zhǎng)期氣候因子變化趨勢(shì)

    1981—2010年長(zhǎng)江中下游地區(qū)單季稻不同RGLs的TEM、PRE和SSD特征及變化趨勢(shì)如下圖3。研究發(fā)現(xiàn), 長(zhǎng)江中下游地區(qū)單季稻分蘗期到成熟期的TEM先升高后降低(圖3-a), BDHD的TEM達(dá)到最高(27.2℃), MSMD的TEM最低(24.5℃)。生殖生長(zhǎng)期前(TDBD)的平均PRE (120.5 mm)和平均SSD (245.5 h)最高, BDHD的平均PRE (26.6 mm)和平均SSD (60.1 h)最低(圖3-b, c)。就單季稻整個(gè)生殖生長(zhǎng)階段(BDMD)而言, TEM增加趨勢(shì)不明顯, 而MSMD的TEM有略微下降趨勢(shì)(–0.017℃ a–1), 其他RGLs的TEM均以增加為主。單季稻MSMD的TEM降低可能是乳熟期和成熟期推遲導(dǎo)致TEM季節(jié)性降低引起的。單季稻BDHD的TEM增加趨勢(shì)最明顯(0.027℃ a–1), 其次為HDMS的TEM (0.013℃ a–1) (圖3-d)。單季稻BDMD的平均PRE呈增加趨勢(shì)(1.27 mm a–1), 其中, HDMS的平均PRE增加趨勢(shì)最明顯(0.71 mm a–1), BDHD和MSMD的平均PRE變化趨勢(shì)不明顯(圖3-e)。單季稻BDMD的平均SSD以減小為主(-0.17 h a–1), 其中, BDHD和MSMD的平均SSD呈減小趨勢(shì)(–0.16 h a–1和–0.45 h a–1), 而HDMS的平均SSD呈增加趨勢(shì)(0.44 h a–1) (圖3-f)。由此可見(jiàn), 單季稻BDHD的氣候變暖速率最強(qiáng), 平均PER變化不明顯,平均SSD呈減小趨勢(shì), HDMS的氣候變暖速率次強(qiáng), 平均PRE和平均SSD均主要呈增加趨勢(shì), 其他生殖生長(zhǎng)階段氣候因子變化趨勢(shì)各異。

    圖2 1981–2010年長(zhǎng)江中下游地區(qū)單季稻生殖生長(zhǎng)期和生殖生長(zhǎng)階段長(zhǎng)度變化趨勢(shì)

    a、b、c、d、e、f、g、h、i和j分別表示分蘗期、孕穗期、抽穗期、乳熟期、成熟期、分蘗期–孕穗期、孕穗期–抽穗期、抽穗期–乳熟期、乳熟期–成熟期和孕穗期–成熟期。各指標(biāo)變化趨勢(shì)由1981–2010年觀測(cè)的各生殖生長(zhǎng)期時(shí)間序列依據(jù)公式(1)計(jì)算得到。*表示在0.05概率水平差異顯著, no表示統(tǒng)計(jì)不顯著。矢量地圖數(shù)據(jù)來(lái)源: 全國(guó)地理信息資源目錄服務(wù)系統(tǒng)(https://www.webmap.cn/)。

    a, b, c, d, e, f, g, h, i, and j represent tillering stage, booting stage, heading stage, milking stage, maturity stage, tillering-booting, booting-heading, heading-milking, milking-maturity, and booting-maturity, respectively. Trends of indexes were calculated according to formula (1) for each reproductive growth period observed from 1981 to 2010. *:< 0.05; no: no significant difference. The above vector map is from the National Catalogue Service for Geographic Information (https://www.webmap.cn/).

    2.3 單季稻生殖生長(zhǎng)期對(duì)氣候因子的敏感性分析

    單季稻不同RGLs對(duì)氣候因子(TEM、PRE和SSD)的敏感性如表1所示。面板模型回歸系數(shù)β1、β2和β3分別表示不同RGLs對(duì)TEM、PRE和SSD的敏感性, 決定系數(shù)(2)表示模型回歸系數(shù)對(duì)因變量的解釋率。結(jié)果表明, TEM與RGLs均呈負(fù)相關(guān)關(guān)系, TEM每升高1℃, BDHD、HDMS和MSMD的平均長(zhǎng)度分別縮短0.77、1.00和1.17 d, 且RGLs對(duì)TEM的敏感性隨生育進(jìn)程增加。PRE和SSD與不同RGLs均呈正相關(guān)關(guān)系, PRE每增加100 mm, BDHD、HDMS和MSMD的平均長(zhǎng)度分別延長(zhǎng)1.24、0.86和0.71 d, 且RGLs對(duì)PRE的敏感性隨生育進(jìn)程減小; SSD每增加100 h, BDHD、HDMS和MSMD的平均長(zhǎng)度分別延長(zhǎng)7.46、8.06和6.88 d??傮w而言, TEM增加不利于延長(zhǎng)RGLs, PRE和SSD增加有利于延長(zhǎng)RGLs。標(biāo)準(zhǔn)化處理不同RGLs對(duì)氣候因子的敏感性表明BDHD、HDMS和MSMD均對(duì)SSD的敏感性最強(qiáng), 對(duì)TEM的敏感性次之, 對(duì)PRE的敏感性最小(表2)。

    圖3 不同生殖生長(zhǎng)期氣候因子變化特征

    a、b和c分別表示不同生殖生長(zhǎng)期的平均溫度、累積降水量和累計(jì)日照時(shí)數(shù); d、e和f分別表示平均溫度、累積降水量和累積日照時(shí)數(shù)的變化趨勢(shì)。TEM、PRE和SSD分別表示平均溫度、累積降水量和累積日照時(shí)數(shù)。TDBD、BDHD、HDMS、MSMD和BDMD分別表示分蘗期–孕穗期、孕穗期–抽穗期、抽穗期–乳熟期、乳熟期–成熟期和孕穗期–成熟期。各指標(biāo)變化趨勢(shì)由1981–2010年觀測(cè)的氣象數(shù)據(jù)時(shí)間序列依據(jù)公式(1)計(jì)算得到。

    a, b, and c represent mean temperature, cumulative precipitation, and cumulative sunshine duration of different reproductive growth periods lengths, respectively. d, e, and f represent trends of mean temperature, cumulative precipitation, and cumulative sunshine duration, respectively. TEM, PRE, and SSD represent mean temperature, cumulative precipitation, and cumulative sunshine duration, respectively. TDBD, BDHD, HDMS, MSMD, and BDMD represent tillering to booting, booting to heading, heading to milking, milking to maturity, and booting to maturity, respectively. The change trends of indexes were calculated by the time sequence of meteorological data observed from 1981 to 2010 based on formula (1).

    表1 面板模型中不同生殖生長(zhǎng)期對(duì)平均溫度、累計(jì)降水量和累計(jì)日照時(shí)數(shù)的相關(guān)系數(shù)

    *:< 0.05;**:< 0.01.

    表2 標(biāo)準(zhǔn)化面板模型中不同生殖生長(zhǎng)期對(duì)平均溫度、累計(jì)降水量和累計(jì)日照時(shí)數(shù)的相關(guān)系數(shù)

    *表示在0.05概率水平差異顯著,**表示在0.01概率水平差異顯著。*:< 0.05;**:< 0.01.

    2.4 氣候變化和技術(shù)進(jìn)步對(duì)生殖生長(zhǎng)期變化趨勢(shì)的影響以及貢獻(xiàn)

    氣候變化和技術(shù)進(jìn)步對(duì)單季稻不同RGLs變化趨勢(shì)的單獨(dú)影響如圖4。研究發(fā)現(xiàn), 氣候變化不利于RGLs延長(zhǎng), 對(duì)不同RGLs均有縮短效應(yīng)。技術(shù)進(jìn)步有利于延長(zhǎng)不同RGLs, 對(duì)HDMS延長(zhǎng)效應(yīng)最大, 延長(zhǎng)趨勢(shì)為0.22 d a–1, 對(duì)全生殖生長(zhǎng)期的延長(zhǎng)趨勢(shì)為0.33 d a–1。氣候變化和技術(shù)進(jìn)步對(duì)RGLs變化趨勢(shì)的綜合作用與技術(shù)進(jìn)步對(duì)RGLs變化趨勢(shì)的單獨(dú)作用相似。

    分離后的氣候變化和技術(shù)進(jìn)步對(duì)單季稻RGLs變化趨勢(shì)的平均貢獻(xiàn)如圖5。研究發(fā)現(xiàn), 氣候變化和技術(shù)進(jìn)步對(duì)單季稻BDMD變化趨勢(shì)的平均貢獻(xiàn)分別為–20.0%和80.0% (圖5-a)。然而, 氣候變化和技術(shù)進(jìn)步的平均貢獻(xiàn)在單季稻BDHD、HDMS和MSMD存在差異。對(duì)于BDHD, 氣候變化的平均貢獻(xiàn)(–65.1%)大于技術(shù)進(jìn)步的平均貢獻(xiàn)(34.9%), 對(duì)于HDMS, 氣候變化的平均貢獻(xiàn)(–24.7%)小于技術(shù)進(jìn)步的平均貢獻(xiàn)(75.3%), 而對(duì)于MSMD, 氣候變化的平均貢獻(xiàn)(–48.1%)略微小于技術(shù)進(jìn)步的平均貢獻(xiàn)(51.9%)。由此可見(jiàn), 技術(shù)進(jìn)步對(duì)RGLs變化趨勢(shì)的延長(zhǎng)起主導(dǎo)作用, 技術(shù)進(jìn)步對(duì)RGLs變化趨勢(shì)的貢獻(xiàn)彌補(bǔ)了氣候變化的不利影響。

    TEM、PRE和SSD對(duì)單季稻RGLs變化趨勢(shì)的平均相對(duì)貢獻(xiàn)如圖5-b。PRE對(duì)不同RGLs變化趨勢(shì)的平均相對(duì)貢獻(xiàn)絕對(duì)值最小。TEM和SSD是影響RGLs變化趨勢(shì)的主要?dú)夂蛞蜃? 兩者對(duì)不同RGLs變化趨勢(shì)的平均相對(duì)貢獻(xiàn)差異大。對(duì)于TDBD, TEM的平均相對(duì)貢獻(xiàn)最大(–69.9%), SSD的平均相對(duì)貢獻(xiàn)次之(–27.7%), 表明該階段單季稻生長(zhǎng)發(fā)育主要受TEM的影響。對(duì)于BDHD, TEM和SSD的平均相對(duì)貢獻(xiàn)相當(dāng)(–50.0%和–47.2%), HDMS變化趨勢(shì)受氣候因子的影響與BDHD類似, TEM和SSD的平均相對(duì)貢獻(xiàn)相當(dāng)(–50.7%和–48.8%), 表明這2個(gè)階段單季稻生長(zhǎng)發(fā)育受制于TEM和SSD。對(duì)于MSMD, SSD的平均相對(duì)貢獻(xiàn)最大(–67.6%), TEM的平均相對(duì)貢獻(xiàn)次之(–21.9%), 表明該階段單季稻生長(zhǎng)發(fā)育主要受SSD的影響。因此, TEM和SSD均為影響長(zhǎng)江流域地區(qū)單季稻RGLs變化趨勢(shì)的主要?dú)夂蛞蜃印?/p>

    圖4 氣候變化和技術(shù)進(jìn)步對(duì)單季稻不同生殖生長(zhǎng)期變化趨勢(shì)的影響

    T_cli、T_tec、T_cli_tec分別表示氣候因素、技術(shù)進(jìn)步、氣候因素和技術(shù)進(jìn)步綜合作用對(duì)不同生殖生長(zhǎng)期長(zhǎng)度的影響趨勢(shì)。生殖生長(zhǎng)期長(zhǎng)度縮寫同圖3。

    T_cli, T_tec, and T_cli_tecindicate the effect trends of climatic factors, technological progress, and the combined of climatic factors and technological progress on different reproductive growth periods lengths, respectively. Abbreviations of different reproductive growth periods lengths are the same as those given in Fig. 3.

    圖5 氣候變化和技術(shù)進(jìn)步對(duì)單季稻生殖生長(zhǎng)期長(zhǎng)度的平均貢獻(xiàn)

    a表示氣候變化和技術(shù)進(jìn)步對(duì)不同生殖生長(zhǎng)期的平均貢獻(xiàn); b表示氣候因子對(duì)不同生殖生長(zhǎng)期的平均相對(duì)貢獻(xiàn)。CRGL_cli和CRGL_tec分別表示氣候因素和技術(shù)進(jìn)步對(duì)不同生殖生長(zhǎng)期長(zhǎng)度的平均貢獻(xiàn); RCRGL_TEM、RCRGL_PRE和RCRGL_SSD分別表示平均溫度、累積降水量和累積日照時(shí)數(shù)對(duì)不同生殖生長(zhǎng)期長(zhǎng)度的平均相對(duì)貢獻(xiàn)。生殖生長(zhǎng)期長(zhǎng)度縮寫同圖3。

    a: the mean contributions of climate change and technological progress to different reproductive growth periods lengths; b: the mean relative contributions of climatic factors to different reproductive growth periods lengths. CRGL_cliand CRGL_tecindicate mean contribution of climatic factors and technological progress to different reproductive growth periods lengths, respectively; RCRGL_TEM, RCRGL_PRE, and RCRGL_SSDindicate mean relative contribution of mean temperature, cumulative precipitation, and cumulative sunshine duration to different reproductive growth periods lengths, respectively. Abbreviations of different reproductive growth periods lengths are the same as those given in Fig. 3.

    3 討論

    3.1 單季稻生殖生長(zhǎng)期對(duì)氣候變化和技術(shù)進(jìn)步的響應(yīng)

    中國(guó)不同地區(qū)水稻生育期和生長(zhǎng)階段長(zhǎng)度的變化趨勢(shì)各不相同[38]。研究發(fā)現(xiàn), 長(zhǎng)江中下游單季稻乳熟期和成熟期顯著推遲, 生殖生長(zhǎng)期尤其是抽穗期到乳熟期顯著延長(zhǎng)(圖1), 該研究結(jié)果與Tao等[32]的結(jié)果一致。抽穗期到乳熟期延長(zhǎng)為灌漿期提供了充足的時(shí)間, 從而有利于生物量累積和產(chǎn)量提高[40]。本研究結(jié)果表明, 過(guò)去30年長(zhǎng)江中下游地區(qū)農(nóng)民采用了長(zhǎng)生殖生長(zhǎng)期單季稻品種, 這也表明延長(zhǎng)單季稻生殖生長(zhǎng)期是適應(yīng)氣候變化的有效措施[38]。

    氣候變化對(duì)作物生育期的影響效應(yīng)因研究對(duì)象和生育期時(shí)段不同而存在差異[12]。本研究表明單季稻RGLs隨平均溫度的升高而縮短, 隨累積降水量和累積日照時(shí)數(shù)的增加而延長(zhǎng), 該結(jié)果與以往一些研究結(jié)果一致[11,25], 而與Liu等[38]的研究不同。其他作物生育期對(duì)溫度、累積降水量和累積日照時(shí)數(shù)的響應(yīng)與本研究結(jié)論一致[7,18,41]。本研究揭示了孕穗期至抽穗期、抽穗期至乳熟期和乳熟期至成熟期對(duì)累積日照時(shí)數(shù)的敏感性最強(qiáng), 對(duì)平均溫度的敏感性次之, 對(duì)累積降水量的敏感性最低。平均溫度和累積日照時(shí)數(shù)對(duì)不同生殖生長(zhǎng)期變化趨勢(shì)的平均相對(duì)貢獻(xiàn)占比均較大, 且兩者的長(zhǎng)期變化趨勢(shì)均不利于生殖生長(zhǎng)期延長(zhǎng), 這進(jìn)一步表明高溫和寡照是長(zhǎng)江中下游地區(qū)單季稻生殖生長(zhǎng)期的主要限制因子。本研究發(fā)現(xiàn)氣候變化對(duì)單季稻生殖生長(zhǎng)期的貢獻(xiàn)小于技術(shù)進(jìn)步, 并且技術(shù)進(jìn)步在延長(zhǎng)生殖生長(zhǎng)期變化趨勢(shì)上占主導(dǎo)作用(圖5), 該結(jié)果與Li等[10]、Wang等[12]、Liu和Dai[27]利用統(tǒng)計(jì)方法和作物模型法的研究結(jié)果一致。因此, 采用長(zhǎng)生殖生長(zhǎng)期品種, 可緩解氣候變化帶來(lái)的生育期縮短的負(fù)面影響[42]。

    3.2 單季稻生殖生長(zhǎng)期適應(yīng)氣候變化的措施

    實(shí)施作物管理有利于協(xié)調(diào)水稻生育進(jìn)程, 充分利用水稻生長(zhǎng)季的有效光溫資源, 獲得最大的光合生產(chǎn)量, 提高結(jié)實(shí)率和千粒重, 實(shí)現(xiàn)高產(chǎn)、優(yōu)質(zhì)和高效栽培[23]。1981—2010年中國(guó)單季稻主產(chǎn)區(qū)高溫持續(xù)時(shí)間和高溫強(qiáng)度顯著上升[43]。湖北中東部、安徽中部、江蘇西部、浙江北部和湖南東北部均易發(fā)生熱脅迫, 在未來(lái)的2021—2050年和2071—2100年2個(gè)時(shí)段, 整個(gè)區(qū)域的熱脅迫顯著增加[44]。長(zhǎng)江中下游地區(qū)熱脅迫主要發(fā)生在7月至8月, 該時(shí)間段正值單季稻生殖生長(zhǎng)階段, 尤其是單季稻開(kāi)花期和灌漿期與最高氣溫時(shí)期重疊[45-46], 張倩等[47]研究也表明長(zhǎng)江中下游地區(qū)熱脅迫主要影響單季稻孕穗期、開(kāi)花期和灌漿期。由本研究各生殖生長(zhǎng)階段氣候因子變化特征可知, 長(zhǎng)江中下游地區(qū)單季稻產(chǎn)量形成關(guān)鍵期處于熱脅迫容易發(fā)生的時(shí)段。且研究表明單季稻乳熟期前后為溫度敏感期[48], 本研究結(jié)果也證實(shí)了這一點(diǎn)(表2)。而這些關(guān)鍵階段的熱脅迫是導(dǎo)致水稻減產(chǎn)的主要原因[29]。因此單季稻生殖生長(zhǎng)期極易遭受熱脅迫的影響, 而避免高溫是預(yù)防熱脅迫的主要策略[47], 培育耐熱品種是保障糧食安全的最有效措施之一[49]。優(yōu)化播種期是緩解氣候變暖不利影響的另一有效措施, Bai等[24]研究表明調(diào)整播種期主要改變了單季稻的開(kāi)花期和成熟期, 而對(duì)生殖生長(zhǎng)期長(zhǎng)度沒(méi)有影響, 因此調(diào)整播種期獲得不同的光溫配置是提高水稻資源利用率和適應(yīng)氣候變化的能力的重要途徑[50]。長(zhǎng)江中下游地區(qū)適當(dāng)早播可提高單季稻品種的光溫資源利用[51], 合理的播種期科學(xué)配套適宜單季稻品種類型對(duì)充分利用當(dāng)?shù)毓鉁刭Y源、發(fā)揮品種潛力、確保高產(chǎn)穩(wěn)產(chǎn)具有十分重要的意義。因此, 單季稻提前播種以避開(kāi)花期高溫、轉(zhuǎn)向現(xiàn)有耐高溫品種以及培育長(zhǎng)生殖生長(zhǎng)期優(yōu)良品種等3項(xiàng)適應(yīng)性措施可在不同程度上有效抵消氣候變化的負(fù)面影響[52]。另外, 合理灌溉、施肥和噴灑外源化學(xué)物質(zhì)可以緩解外界環(huán)境的不利影響。例如, 灌溉是降低水稻花期冠層溫度的一種有效的實(shí)時(shí)栽培措施[46], 合理灌溉可以緩解因氣候變暖導(dǎo)致的開(kāi)花期提前[53];灌漿期熱脅迫影響單季稻品質(zhì), 主要通過(guò)降低葉片抗氧化酶活性, 從而降低光合作用, 導(dǎo)致光合同化物降低等[54], 灌漿期前噴灑外源亞精胺可顯著提高超氧化物歧化酶和過(guò)氧化物酶的活性, 降低丙二醛的積累, 提高熱脅迫下水稻葉片的可溶性糖含量, 維持葉片的滲透平衡, 提高光合和蒸騰速率[55]; 施氮可以增強(qiáng)單季稻生殖生長(zhǎng)期對(duì)熱脅迫的耐受性, 進(jìn)而增加葉片壽命和生長(zhǎng)期長(zhǎng)度, 然而, 需肥量應(yīng)根據(jù)當(dāng)?shù)厮靖魃A段實(shí)際日照狀況計(jì)算得到, 日照不足時(shí)應(yīng)少施肥。例如, 長(zhǎng)江中下游地區(qū)分蘗期到孕穗期和抽穗期到乳熟期日照時(shí)數(shù)呈增加趨勢(shì)(圖3), 這2個(gè)關(guān)鍵生長(zhǎng)階段可以適當(dāng)施肥。

    3.3 問(wèn)題和挑戰(zhàn)

    研究氣候變化和技術(shù)進(jìn)步對(duì)作物生育期長(zhǎng)度的影響十分關(guān)鍵。氣候變化和技術(shù)進(jìn)步對(duì)作物生育期的定量影響最大的問(wèn)題仍然是有效資料的匱乏, 這在理解生育期響應(yīng)氣候變化機(jī)制方面仍存在很大的不確定性。現(xiàn)有長(zhǎng)期作物生育期觀測(cè)資料是不同站點(diǎn)不同技術(shù)措施下得到的, 每個(gè)站點(diǎn)品種類型、施肥量、灌溉量以及農(nóng)藥使用等相關(guān)資料很難準(zhǔn)確獲得, 因此技術(shù)進(jìn)步的貢獻(xiàn)僅是一個(gè)綜合結(jié)果, 不同管理措施難以定量化, 這可能限制了作物適應(yīng)氣候變化的精細(xì)化措施的實(shí)施。以往研究通常用品種替換代替技術(shù)進(jìn)步的影響, 一定程度上可能高估了品種改變對(duì)作物生育期的影響。而基于田間試驗(yàn)的有關(guān)品種變化對(duì)作物生育期影響的資料是有限的[5]。因此, 今后應(yīng)利用田間試驗(yàn)、原位觀測(cè)和模型模擬等多源數(shù)據(jù)開(kāi)展研究[9,18,25]。

    4 結(jié)論

    本研究詳細(xì)地分析了1981—2010年長(zhǎng)江中下游地區(qū)單季稻不同生殖生長(zhǎng)期長(zhǎng)度時(shí)空變化特征, 以及單季稻不同生殖生長(zhǎng)期長(zhǎng)度對(duì)氣候變化和技術(shù)進(jìn)步的響應(yīng)。近30年來(lái)長(zhǎng)江中下游地區(qū)單季稻乳熟期和成熟期推遲, 全生殖生長(zhǎng)期延長(zhǎng), 其中以抽穗期到乳熟期延長(zhǎng)為主。平均溫度升高縮短了單季稻生殖生長(zhǎng)期長(zhǎng)度, 累積降水量和累積日照時(shí)數(shù)增加延長(zhǎng)了單季稻生殖生長(zhǎng)期長(zhǎng)度, 且生殖生長(zhǎng)期長(zhǎng)度對(duì)累積日照時(shí)數(shù)的敏感性最強(qiáng)、對(duì)平均溫度的敏感性次之, 對(duì)降水量的敏感性最小。長(zhǎng)江中下游地區(qū)高溫和寡照不利于單季稻生殖生長(zhǎng)期的延長(zhǎng), 而技術(shù)進(jìn)步彌補(bǔ)了氣候變化對(duì)生殖生長(zhǎng)期的縮短效應(yīng)。因此, 氣候變暖背景下, 長(zhǎng)江中下游地區(qū)采取強(qiáng)耐熱型和長(zhǎng)生殖生長(zhǎng)期單季稻品種仍然是適應(yīng)氣候變化的關(guān)鍵措施。

    [1] Liu Y J, Bachofen C, Wittwer R, Silva D G, Sun Q, Klaus V H, Buchmann N. Using PhenoCams to track crop phenology and explain the effects of different cropping systems on yield., 2022, 195: 103306.

    [2] van Vliet A J H, Bron W A, Mulder S, van der Slikke W, Odé B. Observed climate-induced changes in plant phenology in the Netherlands., 2013, 14: 997–1008.

    [3] Liu F S, Chen Y, Shi W J, Zhang S, Tao F L, Ge Q S. Influences of agricultural phenology dynamic on land surface biophysical process and climate feedback., 2017, 27: 1085–1099.

    [4] Mohapatra S, Tripathy S K, Mohanty A K, Tripathy S. Effect of heat stress on yield and economics of rice (L.) cultivars under different sowing dates., 2021, 23: 38–45.

    [5] Wang Z B, Chen J, Tong W J, Xu C C, Fu C. Impacts of climate change and varietal replacement on winter wheat phenology in the North China Plain., 2018, 12: 251–263.

    [6] Piao S L, Ciais P, Huang Y, Shen Z H, Peng S S, Li J S, Zhou L P, Liu H Y, Ma Y C, Ding Y H, Friedlingstein P, Liu C Z, Tan K, Yu Y Q, Zhang T Y, Fang J Y. The impacts of climate change on water resources and agriculture in China., 2010, 467: 43–51.

    [7] Liu Y J, Qin Y, Ge Q S, Dai J H, Chen Q M. Reponses and sensitivities of maize phenology to climate change from 1981 to 2009 in Henan province, China., 2017, 27: 1072–1084.

    [8] He L, Asseng S, Zhao G, Wu D R, Yang X Y, Zhuang W, Jin N, Yu Q. Impacts of recent climate warming, cultivar changes, and crop management on winter wheat phenology across the Loess Plateau of China., 2015, 200: 135–143.

    [9] Wang Z, Chen J, Li Y, Li C, Zhang L, Chen F. Effects of climate change and cultivar on summer maize phenology., 2016, 10: 509–526.

    [10] Li N, Li Y, Biswas A, Wang J H, Dong H Z, Chen J H, Liu C C, Fan X Y. Impact of climate change and crop management on cotton phenology based on statistical analysis in the main-cotton- planting areas of China., 2021, 298: 126750.

    [11] Bai H Z, Xiao D P. Spatiotemporal changes of rice phenology in China during 1981?2010., 2020, 140: 1483–1494.

    [12] Wang X H, Ciais P, Li L, Ruget F, Vuichard N, Viovy N, Zhou F, Chang J F, Wu X C, Zhao H F. Management outweighs climate change on affecting length of rice growing period for early rice and single rice in China during 1991–2012., 2017, 233: 1–11.

    [13] Ye T, Zong S, Kleidon A, Yuan W P, Wang Y, Shi P J. Impacts of climate warming, cultivar shifts, and phenological dates on rice growth period length in China after correction for seasonal shift effects., 2019, 155: 127–143.

    [14] Zhang S, Tao F L, Zhang Z. Rice reproductive growth duration increased despite of negative impacts of climate warming across China during 1981?2009., 2014, 54: 70–83.

    [15] Li Z G, Yang P, Tang H J, Wu W B, Yin H, Liu Z H, Zhang L. Response of maize phenology to climate warming in Northeast China between 1990 and 2012., 2013, 14: 39–48.

    [16] Mo F, Sun M, Liu X Y, Wang J Y, Zhang X C, Ma B L, Xiong Y C. Phenological responses of spring wheat and maize to changes in crop management and rising temperatures from 1992 to 2013 across the Loess Plateau., 2016, 196: 337–347.

    [17] Wang H L, Gan Y T, Wang R Y, Niu J Y, Zhao H, Yang Q G, Li G C. Phenological trends in winter wheat and spring cotton in response to climate changes in northwest China., 2008, 148: 1242–1251.

    [18] Wang Z B, Chen J, Xing F F, Han Y C, Chen F, Zhang L F, Li Y B, Li C D. Response of cotton phenology to climate change on the North China Plain from 1981 to 2012., 2017, 7: 6628.

    [19] Xiao D P, Zhao Y X, Bai H Z, Hu Y K, Cao J S. Impacts of climate warming and crop management on maize phenology in northern China., 2019, 11: 892–903.

    [20] Fatima Z, Ur-Rehman A, Abbas G, Iqbal P, Zakir I, Khan M A, Kamal G M, Ahmed M, Ahmad S. Quantification of climate warming and crop management impacts on phenology of pulses- based cropping systems., 2020, 15: 107–123.

    [21] Li K N, Yang X G, Tian H Q, Pan S F, Liu Z J, Lu S. Effects of changing climate and cultivar on the phenology and yield of winter wheat in the North China Plain., 2016, 60: 21–32.

    [22] Liu Y J, Chen Q M, Ge Q S, Ge Q S, Dai J H, Qin Y, Dai L, Zou X T, Chen J. Modelling the impacts of climate change and crop management on phenological trends of spring and winter wheat in China., 2018, 248: 518–526.

    [23] Hu X Y, Huang Y, Sun W J, Yu L F. Shifts in cultivar and planting date have regulated rice growth duration under climate warming in China since the early 1980s., 2017, 247: 34–41.

    [24] Bai H Z, Xiao D P, Zhang H, Tao F L, Hu Y H. Impact of warming climate, sowing date, and cultivar shift on rice phenology across China during 1981–2010., 2019, 63: 1077–1089.

    [25] Chen J, Liu Y J, Zhou W M, Zhang J, Pan T. Effects of climate change and crop management on changes in rice phenology in China from 1981 to 2010., 2021, 101: 6311–6319.

    [26] Zhang T Y, Huang Y, Yang X G. Climate warming over the past three decades has shortened rice growth duration in China and cultivar shifts have further accelerated the process for late rice., 2013, 19: 563–570.

    [27] Liu Y J, Dai L. Modelling the impacts of climate change and crop management measures on soybean phenology in China., 2020, 262: 121271.

    [28] Wu D R, Wang P J, Jiang C Y, Yang J Y, Huo Z G, Yu Q. Measured phenology response of unchanged crop varieties to long- term historical climate change., 2018, 13: 47–58.

    [29] Shi P, Tang L, Wang L, Sun T, Liu L, Cao W, Zhu Y. Post-heading heat stress in rice of south China during 1981–2010., 2015, 10: e0130642.

    [30] Rehmani M I A, Ding C Q, Li G H, Ata-Ul-Karim S T, Hadifa A, Bashir M A, Hashem M, Alamri S, Al-Zubair F, Ding Y F. Vulnerability of rice production to temperature extremes during rice reproductive stage in Yangtze River valley, China., 2021, 33: 101599.

    [31] 劉書通. 長(zhǎng)江中下游地區(qū)水稻生產(chǎn)能力分析. 中國(guó)農(nóng)業(yè)科學(xué)院碩士學(xué)位論文, 北京, 2014.

    Liu S T. Analysis of Rice Production Capacity in the Middle and Lower Reaches of the Yangtze River. MS Thesis of Chinese Academy of Agricultural Sciences, Beijing, China, 2014 (in Chinese with English abstract).

    [32] Tao F L, Zhang Z, Shi W J, Liu Y J, Xiao D P, Zhang S, Zhu Z, Wang M, Liu F S. Single rice growth period was prolonged by cultivars shifts, but yield was damaged by climate change during 1981–2009 in China, and late rice was just opposite., 2013, 19: 3200–3209.

    [33] Jiang Y L, Yin X G, Wang X H, Zhang L, Lu Z, Lei Y D, Chu Q Q, Chen F. Impacts of global warming on the cropping systems of China under technical improvements from 1961 to 2016., 2021, 113: 187–199.

    [34] 陳風(fēng)波, 馬志雄, 陳培勇. 農(nóng)戶水稻種植模式選擇現(xiàn)狀及其影響因素分析: 對(duì)長(zhǎng)江中下游地區(qū)四個(gè)省的調(diào)查. 農(nóng)業(yè)經(jīng)濟(jì)與管理, 2011, (3): 62–73.

    Chen F B, Ma Z X, Chen P Y. Selecting status of farmers’ rice cropping patterns and its impact factor analysis: four provinces survey along the Yangtze River., 2011, (3): 62–73 (in Chinese with English abstract).

    [35] 王旭, 張?zhí)焐? 劉懂偉, 王曉亮, 李輝. 長(zhǎng)江中下游水稻種植影響因素及變化分析. 現(xiàn)代農(nóng)業(yè)科技, 2016, (6): 54–55.

    Wang X, Zhang T S, Liu D W, Wang X L, Li H. Analysis of factors influencing rice cultivation and changes in the middle and lower reaches of the Yangtze River., 2016, (6): 54–55 (in Chinese with English abstract).

    [36] Lobell D B, Burke M B. On the use of statistical models to predict crop yield responses to climate change., 2010, 150: 1443–1452.

    [37] Zhang S, Tao F L, Zhang Z. Changes in extreme temperatures and their impacts on rice yields in southern China from 1981 to 2009., 2016, 189: 43–50.

    [38] Liu Y J, Zhou W M, Ge Q S. Spatiotemporal changes of rice phenology in China under climate change from 1981 to 2010., 2019, 157: 261–277.

    [39] Zhang Y, Zhao Y X, Sun Q. Increasing maize yields in Northeast China are more closely associated with changes in crop timing than with climate warming., 2021, 16: 054052.

    [40] Zhu P, Jin Z N, Zhuang Q L, Ciais P, Bernacchi C, Wang X H, Makowski D, Lobell D. The important but weakening maize yield benefit of grain filling prolongation in the US Midwest., 2018, 24: 4718–4730.

    [41] Guo E L, Zhang J Q, Wang Y F, Alu S, Wang R, Li D J, Ha S. Assessing non-linear variation of temperature and precipitation for different growth periods of maize and their impacts on phenology in the midwest of Jilin province, China., 2017, 132: 685–699.

    [42] Tao F L, Zhang S L, Zhang Z, R?tter R P. Maize growing duration was prolonged across China in the past three decades under the combined effects of temperature, agronomic management, and cultivar shift., 2014, 20: 3686–3699.

    [43] Sun X S, Long Z W, Song G P, Chen C Q. High-temperature episodes with spatial-temporal variation impacted middle-season rice yield in China., 2018, 110: 961?969.

    [44] Zhang L, Yang B Y, Li S, Hou Y Y, Huang D P. Potential rice exposure to heat stress along the Yangtze River in China under RCP8.5 scenario., 2018, 248: 185–196.

    [45] Meng L, Wang C Y, Zhang J Q. Heat injury risk assessment for single-cropping rice in the middle and lower reaches of the Yangtze River under climate change., 2016, 30: 426–443.

    [46] Wang Y L, Wang L, Zhou J X, Hu S B, Chen H Z, Xiang J, Zhang Y K, Zeng Y J, Shi Q H, Zhu D F, Zhang Y P. Research progress on heat stress of rice at flowering stage., 2019, 26: 1–10.

    [47] 張倩, 趙艷霞, 王春乙. 長(zhǎng)江中下游地區(qū)高溫?zé)岷?duì)水稻的影響. 災(zāi)害學(xué), 2011, 26(4): 58–62.

    Zhang Q, Zhao Y X, Wang C Y. Study on the impact of high temperature damage to rice in the lower and middle reaches of the Yangtze River., 2011, 26(4): 58–62. (in Chinese with English abstract).

    [48] Sanchez B, Rasmussen A, Porter J R. Temperatures and the growth and development of maize and rice: a review., 2014, 20: 408–417.

    [49] Wu C, Cui K H, Li Q, Li L Y, Wang W C, Hu Q Q, Ding Y F, Li G H, Fahad S, Huang J L, Nie L X, Peng S B. Estimating the yield stability of heat-tolerant rice genotypes under various heat conditions across reproductive stages: a 5-year case study., 2021, 11: 13604.

    [50] 鄧南燕. 中國(guó)水稻產(chǎn)量差評(píng)估及長(zhǎng)江中下游地區(qū)增產(chǎn)途徑探究.華中農(nóng)業(yè)大學(xué)博士學(xué)位論文, 湖北武漢, 2018.

    Deng N Y. Estimation of Rice Yield Gaps in China and Exploration of Approaches to Improve Yield in Middle and Lower Reaches of Yangtze River Valley. PhD Dissertation of Huazhong Agricultural University, Wuhan, Hubei, China, 2018 (in Chinese with English abstract).

    [51] 王文婷. 沿江地區(qū)溫度要素對(duì)優(yōu)質(zhì)粳稻產(chǎn)量與品質(zhì)的影響研究. 揚(yáng)州大學(xué)博士學(xué)位論文, 江蘇揚(yáng)州, 2021.

    Wang W T. Study on the Effect of Temperature-Light Factors on the Formation of Quality and Yield ofRice along the Yangtze River. PhD Dissertation of Yangzhou University, Yangzhou, Jiangsu, China, 2021 (in Chinese with English abstract).

    [52] Xu C C, Wu W X, Ge Q S, Zhou Y, Lin Y M, Li Y M. Simulating climate change impacts and potential adaptations on rice yields in the Sichuan Basin, China., 2015, 22: 565–594.

    [53] Muleke A, Harrison M, de Voil P, Hunt L, Liu K, Yanotti M, Eisner R. Earlier crop flowering caused by global warming alleviated by irrigation., 2022, 17: 044032.

    [54] Liu Q H, Wu X, Ma J Q, Li T, Zhou X B, Guo T. Effects of high air temperature on rice grain quality and yeld under field condition., 2013, 105: 446.

    [55] Tang S, Zhang H X, Li L L, Liu X, Chen L, Chen W Z, Ding Y F. Exogenous spermidine enhances the photosynthetic and antioxidant capacity of rice under heat stress during early grain-filling period., 2018, 45: 911–921.

    Response of reproductive growth period length to climate warming and technological progress in the middle and lower reaches of the Yangtze River during 1981–2010 in single-cropping rice

    LIU Er-Hua1,2,3, ZHOU Guang-Sheng1,2,3,4,*, WU Bing-Yi2, SONG Yan-Ling1,3, HE Qi-Jin4, LYU Xiao-Min1,3, and ZHOU Meng-Zi1,3

    1State Key Laboratory of Severe Weather, Chinese Academy of Meteorological Sciences, Beijing 100081, China;2Institute of Atmospheric Sciences, Fudan University, Shanghai 200439, China;3Joint Eco-Meteorological Laboratory of Chinese Academy of Meteorological Sciences and Zhengzhou University, Zhengzhou 450001, Henan, China;4College of Resources and Environmental Sciences, China Agricultural University, Beijing 100193, China

    Crop growth period length is closely linked to climate change and technological progress. Even though the extensive researches conducting on crop growth period length variation and its response to climate change, particularly temperature change, the response of reproductive growth periods lengths (RGLs) to climate change and technological progress remains unclear. Based on the reproductive growth periods and meteorological data of single-cropping rice in the middle and lower reaches of the Yangtze River (MLYR) during 1981?2010, the trends of the RGLs (including booting to heading (BDHD), heading to milking (HDMS), milking to maturity (MSMD), and booting to maturity (BDMD)) and climatic variations were analyzed. In addition, to explore the confounding effects of climate change and technological progress on the RGLs, the sensitivities of the RGLs to mean temperature (TEM), cumulative precipitation (PRE), and cumulative sunshine duration (SSD) were measured. The results showed that the BDMD had an extension trend (0.24 d a–1), among which the extension trend in the HDMS (0.16 d a–1) was the most obvious in RGLs, while the extension trends of BDHD (0.03 d a–1) and MSMD (0.05 d a–1) were not significant. High temperature and low sunshine duration were unfavorable to the extension of the RGLs. The mean relative contributions of TEM to the BDHD, HDMS, and MSMD were –50.0%, –50.7%, and –21.9%, which were –47.2%, –48.7%, and –67.6% in terms of SSD, respectively. Technological progress compensated for the adverse impacts of climate change on the trends of different RGLs. These results suggested that cultivar selection and agronomic management were the effective adaptation strategies benefiting for the stable and high yield of single-cropping rice. Single-cropping rice cultivar with longer RGLs and heat-tolerant may be suitable to cope with the continuous climate change in the future.

    single-cropping rice; climate change; reproductive growth period lengths; technological progress

    10.3724/SP.J.1006.2023.22026

    本研究由國(guó)家重點(diǎn)研究計(jì)劃項(xiàng)目(2018YFA0606103), 國(guó)家自然科學(xué)基金重點(diǎn)項(xiàng)目(42130514, 4213000565)和中國(guó)氣象科學(xué)研究院基本科研業(yè)務(wù)費(fèi)專項(xiàng)(2020Z004)資助。

    This study is supported by the National Key Research and Development Program of China (2018YFA0606103), the National Natural Science Foundation of China (42130514, 4213000565), and the Basic Research Fund of Chinese Academy of Meteorological Sciences (2020Z004).

    周廣勝, E-mail: zhougs@cma.gov.cn

    E-mail: leh4179@163.com

    2022-04-30;

    2022-10-10;

    2022-10-31.

    URL: https://kns.cnki.net/kcms/detail/11.1809.S.20221028.1445.002.html

    This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

    猜你喜歡
    長(zhǎng)江中下游地區(qū)單季稻乳熟期
    秸稈還田與灌溉方式對(duì)春玉米根冠抗氧化能力的影響
    飼用小黑麥在濟(jì)源地區(qū)的刈割期及刈割次數(shù)研究
    2017年長(zhǎng)江中下游地區(qū)一次持續(xù)性異常降水過(guò)程分析
    攀西地區(qū)冬閑田種植晚熟型燕麥的最佳刈割期研究
    長(zhǎng)江中下游地區(qū)梅雨期異常年降水及大氣熱源分布特征
    青藏高原春季土壤濕度對(duì)長(zhǎng)江中下游地區(qū)初夏短期氣候影響的數(shù)值模擬
    2015年貴池區(qū)單季稻肥料配方驗(yàn)證試驗(yàn)
    論青藏高原溫度對(duì)長(zhǎng)江中下游地區(qū)降水的影響
    增效控失肥在粳雜甬優(yōu)8號(hào)上的施肥效果研究
    貴池區(qū)水稻機(jī)插秧高效水肥運(yùn)籌試驗(yàn)研究
    1024香蕉在线观看| 少妇精品久久久久久久| 99精国产麻豆久久婷婷| 亚洲欧美清纯卡通| 亚洲成国产人片在线观看| 香蕉丝袜av| 十八禁高潮呻吟视频| 女人爽到高潮嗷嗷叫在线视频| 日韩一卡2卡3卡4卡2021年| 亚洲精品中文字幕一二三四区 | 欧美久久黑人一区二区| 99久久精品国产亚洲精品| 日本一区二区免费在线视频| 国产在线视频一区二区| www.999成人在线观看| 亚洲伊人色综图| 亚洲成人免费av在线播放| 在线观看www视频免费| a 毛片基地| 亚洲色图综合在线观看| 亚洲少妇的诱惑av| 久久久国产一区二区| 最近中文字幕2019免费版| 免费在线观看视频国产中文字幕亚洲 | 日本五十路高清| 另类精品久久| 精品少妇一区二区三区视频日本电影| 色婷婷久久久亚洲欧美| 亚洲综合色网址| 国产精品亚洲av一区麻豆| 高清黄色对白视频在线免费看| 涩涩av久久男人的天堂| 18禁观看日本| 女人被躁到高潮嗷嗷叫费观| 黄色 视频免费看| 9191精品国产免费久久| 国产成人啪精品午夜网站| 亚洲精品国产一区二区精华液| 各种免费的搞黄视频| 少妇的丰满在线观看| 人妻久久中文字幕网| 一级片'在线观看视频| 欧美国产精品一级二级三级| 久久久国产一区二区| 精品一区在线观看国产| 国产日韩一区二区三区精品不卡| 侵犯人妻中文字幕一二三四区| 人人妻,人人澡人人爽秒播| 欧美老熟妇乱子伦牲交| 亚洲精品日韩在线中文字幕| 成人国产av品久久久| 十八禁网站网址无遮挡| 亚洲一码二码三码区别大吗| 99久久国产精品久久久| 成人三级做爰电影| 久久精品国产综合久久久| 精品福利观看| 久久精品aⅴ一区二区三区四区| 亚洲av电影在线进入| 国产精品成人在线| 成人18禁高潮啪啪吃奶动态图| 夜夜骑夜夜射夜夜干| 国产无遮挡羞羞视频在线观看| videos熟女内射| 国产一区二区三区综合在线观看| 国产黄色免费在线视频| 午夜激情久久久久久久| 国产老妇伦熟女老妇高清| 丝瓜视频免费看黄片| h视频一区二区三区| 久久中文字幕一级| 青草久久国产| av一本久久久久| 少妇猛男粗大的猛烈进出视频| 日韩一区二区三区影片| 两个人免费观看高清视频| 国产有黄有色有爽视频| 亚洲欧美成人综合另类久久久| 久久影院123| 五月天丁香电影| 成年av动漫网址| 久久狼人影院| 亚洲精品一区蜜桃| 日韩欧美免费精品| 成人亚洲精品一区在线观看| 黑丝袜美女国产一区| 国产免费现黄频在线看| 亚洲国产精品一区二区三区在线| 中文欧美无线码| 亚洲一卡2卡3卡4卡5卡精品中文| 99精品久久久久人妻精品| 亚洲国产av影院在线观看| 首页视频小说图片口味搜索| 一级片'在线观看视频| 午夜成年电影在线免费观看| 国产精品秋霞免费鲁丝片| 少妇人妻久久综合中文| 色婷婷av一区二区三区视频| 啪啪无遮挡十八禁网站| 亚洲国产欧美一区二区综合| 亚洲av电影在线进入| 一二三四社区在线视频社区8| 亚洲 国产 在线| 久久久久国产精品人妻一区二区| 亚洲精品av麻豆狂野| 国产成+人综合+亚洲专区| 久久精品亚洲熟妇少妇任你| 免费一级毛片在线播放高清视频 | av天堂久久9| 岛国在线观看网站| 18禁观看日本| 亚洲久久久国产精品| 啪啪无遮挡十八禁网站| 成人国产av品久久久| 多毛熟女@视频| 久久精品国产亚洲av高清一级| av网站免费在线观看视频| 青草久久国产| 久久久久国产一级毛片高清牌| 亚洲一卡2卡3卡4卡5卡精品中文| 免费观看a级毛片全部| 一本—道久久a久久精品蜜桃钙片| 久久99热这里只频精品6学生| 亚洲欧美日韩另类电影网站| 一二三四社区在线视频社区8| 亚洲人成电影观看| 脱女人内裤的视频| 97精品久久久久久久久久精品| 男女无遮挡免费网站观看| 国产野战对白在线观看| 狠狠狠狠99中文字幕| 性少妇av在线| 日日摸夜夜添夜夜添小说| 悠悠久久av| 日韩 欧美 亚洲 中文字幕| 在线精品无人区一区二区三| 国产男女超爽视频在线观看| 亚洲人成77777在线视频| 搡老乐熟女国产| 亚洲伊人久久精品综合| 狠狠狠狠99中文字幕| 电影成人av| 亚洲精品久久久久久婷婷小说| 后天国语完整版免费观看| 老司机午夜十八禁免费视频| 99热国产这里只有精品6| 欧美日韩亚洲国产一区二区在线观看 | 午夜免费鲁丝| 免费在线观看视频国产中文字幕亚洲 | 亚洲精品国产一区二区精华液| 制服人妻中文乱码| 欧美乱码精品一区二区三区| 午夜激情av网站| 少妇粗大呻吟视频| 在线观看一区二区三区激情| 黄色 视频免费看| 婷婷色av中文字幕| 18禁观看日本| 午夜免费观看性视频| 精品人妻一区二区三区麻豆| 成人国产一区最新在线观看| 黑人巨大精品欧美一区二区mp4| 日韩制服丝袜自拍偷拍| 亚洲精品一区蜜桃| 新久久久久国产一级毛片| 一个人免费看片子| 亚洲成av片中文字幕在线观看| 香蕉丝袜av| 大香蕉久久成人网| 色精品久久人妻99蜜桃| 国产精品偷伦视频观看了| 老鸭窝网址在线观看| 中文欧美无线码| 成人三级做爰电影| 精品卡一卡二卡四卡免费| 精品少妇一区二区三区视频日本电影| 国产深夜福利视频在线观看| 国产免费av片在线观看野外av| avwww免费| 久久精品人人爽人人爽视色| av网站免费在线观看视频| av片东京热男人的天堂| 夜夜夜夜夜久久久久| 国产高清国产精品国产三级| 国产av国产精品国产| 欧美精品啪啪一区二区三区 | 国产欧美日韩一区二区三区在线| 天天影视国产精品| 国产在线视频一区二区| 大型av网站在线播放| 一区二区av电影网| av天堂久久9| 色婷婷久久久亚洲欧美| 18禁观看日本| 国产精品久久久人人做人人爽| 欧美激情高清一区二区三区| 久久香蕉激情| 国产精品 国内视频| 久久久精品国产亚洲av高清涩受| 精品人妻在线不人妻| 美女扒开内裤让男人捅视频| 自拍欧美九色日韩亚洲蝌蚪91| 人妻久久中文字幕网| 99九九在线精品视频| 亚洲三区欧美一区| 久久精品国产亚洲av香蕉五月 | 91字幕亚洲| 欧美黑人精品巨大| 国产av一区二区精品久久| 51午夜福利影视在线观看| 久久九九热精品免费| 国产精品欧美亚洲77777| 99香蕉大伊视频| 国产欧美日韩一区二区精品| bbb黄色大片| 国产精品 欧美亚洲| 中文字幕高清在线视频| 丝袜在线中文字幕| a级片在线免费高清观看视频| 纵有疾风起免费观看全集完整版| 亚洲熟女精品中文字幕| 欧美中文综合在线视频| 中文字幕精品免费在线观看视频| 又紧又爽又黄一区二区| avwww免费| 中文字幕人妻丝袜制服| 一区二区av电影网| 国产精品一区二区免费欧美 | 国产91精品成人一区二区三区 | 少妇被粗大的猛进出69影院| 欧美日韩亚洲综合一区二区三区_| 日本av免费视频播放| 高清在线国产一区| 老司机亚洲免费影院| 亚洲国产毛片av蜜桃av| 精品熟女少妇八av免费久了| 国产精品香港三级国产av潘金莲| 狠狠婷婷综合久久久久久88av| 久久性视频一级片| 日日摸夜夜添夜夜添小说| 新久久久久国产一级毛片| 免费观看人在逋| 亚洲综合色网址| 老司机亚洲免费影院| 国产免费av片在线观看野外av| 免费看十八禁软件| 精品一区二区三区av网在线观看 | 国产欧美日韩综合在线一区二区| 亚洲精品av麻豆狂野| 国产在线一区二区三区精| 十八禁网站网址无遮挡| 大型av网站在线播放| 美女主播在线视频| 看免费av毛片| 少妇猛男粗大的猛烈进出视频| 一区二区三区激情视频| 欧美日本中文国产一区发布| 国产国语露脸激情在线看| 99久久综合免费| 18禁观看日本| 欧美xxⅹ黑人| 黑人猛操日本美女一级片| 日韩中文字幕欧美一区二区| 亚洲,欧美精品.| 中文字幕人妻熟女乱码| 亚洲av成人一区二区三| 日韩视频在线欧美| 麻豆av在线久日| 色视频在线一区二区三区| 久久久精品区二区三区| 99国产精品免费福利视频| 中文字幕人妻熟女乱码| 欧美日韩av久久| 国产精品自产拍在线观看55亚洲 | 日本猛色少妇xxxxx猛交久久| 乱人伦中国视频| 成人黄色视频免费在线看| 国产精品一区二区精品视频观看| 大型av网站在线播放| 99九九在线精品视频| 精品少妇黑人巨大在线播放| 天堂8中文在线网| 80岁老熟妇乱子伦牲交| 国产成人啪精品午夜网站| 91成人精品电影| 精品熟女少妇八av免费久了| 亚洲国产欧美日韩在线播放| 如日韩欧美国产精品一区二区三区| 丝袜脚勾引网站| 中文字幕另类日韩欧美亚洲嫩草| 99久久99久久久精品蜜桃| 一本一本久久a久久精品综合妖精| 熟女少妇亚洲综合色aaa.| 男女午夜视频在线观看| 日本vs欧美在线观看视频| 一二三四在线观看免费中文在| 黄频高清免费视频| 国产亚洲精品一区二区www | 美女高潮喷水抽搐中文字幕| 午夜成年电影在线免费观看| 在线观看一区二区三区激情| 可以免费在线观看a视频的电影网站| 亚洲全国av大片| 国产成人影院久久av| 丝袜脚勾引网站| 51午夜福利影视在线观看| 9191精品国产免费久久| 免费高清在线观看日韩| 久久精品国产综合久久久| 国产精品九九99| 国产人伦9x9x在线观看| 考比视频在线观看| 久久久久精品国产欧美久久久 | 一区二区av电影网| 国产深夜福利视频在线观看| 日本wwww免费看| 在线亚洲精品国产二区图片欧美| 男人操女人黄网站| 悠悠久久av| 亚洲精品久久久久久婷婷小说| 97人妻天天添夜夜摸| 91av网站免费观看| av不卡在线播放| 精品久久久久久久毛片微露脸 | 欧美黄色片欧美黄色片| 国产精品一区二区在线不卡| 国产一卡二卡三卡精品| 两个人看的免费小视频| 国产精品国产av在线观看| 国产激情久久老熟女| 巨乳人妻的诱惑在线观看| 性少妇av在线| 老熟女久久久| 菩萨蛮人人尽说江南好唐韦庄| 男女之事视频高清在线观看| 精品一区在线观看国产| 国产真人三级小视频在线观看| 18禁裸乳无遮挡动漫免费视频| 精品亚洲乱码少妇综合久久| 国产97色在线日韩免费| 国产亚洲av片在线观看秒播厂| 老司机午夜十八禁免费视频| 精品亚洲成国产av| 亚洲伊人久久精品综合| 在线观看人妻少妇| 一本—道久久a久久精品蜜桃钙片| 十八禁高潮呻吟视频| 午夜老司机福利片| 新久久久久国产一级毛片| 菩萨蛮人人尽说江南好唐韦庄| 亚洲专区中文字幕在线| 女警被强在线播放| av线在线观看网站| 桃花免费在线播放| 亚洲第一av免费看| 国产成人啪精品午夜网站| 成年动漫av网址| av网站在线播放免费| 天天影视国产精品| av网站在线播放免费| 亚洲国产看品久久| 亚洲欧洲日产国产| 蜜桃在线观看..| av天堂久久9| 1024香蕉在线观看| 久久久久久久精品精品| 美女高潮到喷水免费观看| 啦啦啦中文免费视频观看日本| 精品人妻熟女毛片av久久网站| 国产高清国产精品国产三级| 男女床上黄色一级片免费看| 日日爽夜夜爽网站| 大香蕉久久网| 午夜91福利影院| 90打野战视频偷拍视频| 夜夜骑夜夜射夜夜干| 免费观看av网站的网址| 两性午夜刺激爽爽歪歪视频在线观看 | 日本wwww免费看| 日韩制服丝袜自拍偷拍| 首页视频小说图片口味搜索| 精品国产一区二区三区四区第35| 夫妻午夜视频| 国产欧美日韩一区二区三 | 美女中出高潮动态图| 狂野欧美激情性xxxx| 免费少妇av软件| 欧美另类亚洲清纯唯美| 蜜桃在线观看..| 丰满迷人的少妇在线观看| 亚洲精品一二三| 正在播放国产对白刺激| 精品人妻熟女毛片av久久网站| 一级毛片女人18水好多| 久久女婷五月综合色啪小说| 国产精品久久久久成人av| 久久精品国产亚洲av高清一级| 亚洲一区中文字幕在线| 亚洲精品久久久久久婷婷小说| 建设人人有责人人尽责人人享有的| 91国产中文字幕| 老司机靠b影院| 99热国产这里只有精品6| 99久久人妻综合| 十分钟在线观看高清视频www| 国产在线观看jvid| www.熟女人妻精品国产| 一本一本久久a久久精品综合妖精| 免费在线观看日本一区| 成人免费观看视频高清| 亚洲精品国产av成人精品| 久久精品国产亚洲av高清一级| 欧美日韩av久久| 18在线观看网站| 不卡一级毛片| 日本一区二区免费在线视频| 丝袜喷水一区| 国产精品 欧美亚洲| 亚洲成人国产一区在线观看| 日日爽夜夜爽网站| www.av在线官网国产| 欧美另类亚洲清纯唯美| 国产在线一区二区三区精| av超薄肉色丝袜交足视频| 免费在线观看日本一区| 久久精品熟女亚洲av麻豆精品| 午夜老司机福利片| 男人爽女人下面视频在线观看| 大香蕉久久网| 多毛熟女@视频| 精品一区二区三卡| 亚洲精品乱久久久久久| 12—13女人毛片做爰片一| 叶爱在线成人免费视频播放| 制服人妻中文乱码| 狂野欧美激情性xxxx| 国产精品久久久久久精品电影小说| 国产亚洲欧美在线一区二区| 人人澡人人妻人| 日韩欧美一区二区三区在线观看 | 黄色毛片三级朝国网站| 亚洲成人国产一区在线观看| 婷婷丁香在线五月| 另类亚洲欧美激情| 黑人欧美特级aaaaaa片| 两个人免费观看高清视频| 两人在一起打扑克的视频| av在线老鸭窝| 欧美黄色淫秽网站| 日本五十路高清| 操美女的视频在线观看| 一区二区日韩欧美中文字幕| 国产成人av激情在线播放| www.av在线官网国产| 国产亚洲精品一区二区www | 欧美激情高清一区二区三区| 亚洲专区国产一区二区| 别揉我奶头~嗯~啊~动态视频 | 多毛熟女@视频| 欧美黑人欧美精品刺激| 国产成人a∨麻豆精品| 欧美激情高清一区二区三区| 侵犯人妻中文字幕一二三四区| 国产色视频综合| 欧美+亚洲+日韩+国产| 久久久久视频综合| 国产又色又爽无遮挡免| 91精品国产国语对白视频| 99精国产麻豆久久婷婷| 狂野欧美激情性bbbbbb| 操出白浆在线播放| 午夜影院在线不卡| 丰满迷人的少妇在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 777米奇影视久久| 少妇被粗大的猛进出69影院| 亚洲色图综合在线观看| 亚洲欧美精品自产自拍| 激情视频va一区二区三区| 大香蕉久久成人网| 亚洲精品一区蜜桃| 免费av中文字幕在线| www.999成人在线观看| 丰满少妇做爰视频| 中国美女看黄片| 欧美 日韩 精品 国产| 首页视频小说图片口味搜索| 1024香蕉在线观看| 热re99久久国产66热| 69av精品久久久久久 | 亚洲国产精品999| 看免费av毛片| 香蕉丝袜av| 水蜜桃什么品种好| 操出白浆在线播放| 国产精品香港三级国产av潘金莲| 亚洲情色 制服丝袜| 黑丝袜美女国产一区| av免费在线观看网站| 久久久国产成人免费| av在线老鸭窝| 9色porny在线观看| 成在线人永久免费视频| 无限看片的www在线观看| 亚洲中文av在线| av电影中文网址| av超薄肉色丝袜交足视频| 亚洲精品中文字幕在线视频| 在线av久久热| 亚洲精品久久成人aⅴ小说| 久久影院123| 久久国产精品大桥未久av| 国产精品1区2区在线观看. | 19禁男女啪啪无遮挡网站| 99国产精品一区二区三区| 欧美日韩成人在线一区二区| 亚洲少妇的诱惑av| 成人影院久久| 久久国产精品大桥未久av| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲精品成人av观看孕妇| 国产成人一区二区三区免费视频网站| 久热这里只有精品99| 久久国产亚洲av麻豆专区| 欧美xxⅹ黑人| 黄色 视频免费看| 亚洲国产欧美在线一区| 一本一本久久a久久精品综合妖精| av电影中文网址| 欧美日韩成人在线一区二区| 精品亚洲成国产av| 青青草视频在线视频观看| 国产有黄有色有爽视频| 性色av一级| 99热网站在线观看| 欧美激情 高清一区二区三区| 99九九在线精品视频| 亚洲精品久久久久久婷婷小说| 十八禁高潮呻吟视频| av天堂久久9| 俄罗斯特黄特色一大片| 少妇 在线观看| 日韩一卡2卡3卡4卡2021年| 一区二区日韩欧美中文字幕| 丝袜美足系列| 男女边摸边吃奶| 日韩制服丝袜自拍偷拍| 欧美精品av麻豆av| 成在线人永久免费视频| 亚洲中文字幕日韩| 狂野欧美激情性xxxx| 十八禁网站网址无遮挡| 美女国产高潮福利片在线看| 欧美人与性动交α欧美精品济南到| 日韩视频在线欧美| 国产成人免费无遮挡视频| 69av精品久久久久久 | 国产精品一区二区在线观看99| a 毛片基地| 国产麻豆69| 精品人妻在线不人妻| 亚洲av成人一区二区三| 桃红色精品国产亚洲av| 黄色 视频免费看| 亚洲精品国产一区二区精华液| av有码第一页| 亚洲精品第二区| 国产97色在线日韩免费| 99国产精品免费福利视频| 女性被躁到高潮视频| 操出白浆在线播放| 久久毛片免费看一区二区三区| 咕卡用的链子| 两性夫妻黄色片| 男人添女人高潮全过程视频| 91精品伊人久久大香线蕉| 麻豆av在线久日| 丁香六月天网| 亚洲欧美日韩高清在线视频 | av网站在线播放免费| 国产亚洲欧美在线一区二区| 亚洲五月色婷婷综合| 免费不卡黄色视频| 91精品三级在线观看| 老熟妇乱子伦视频在线观看 | 动漫黄色视频在线观看| av网站在线播放免费| 成人亚洲精品一区在线观看| 亚洲中文字幕日韩| 中文精品一卡2卡3卡4更新| 日韩熟女老妇一区二区性免费视频| 国产视频一区二区在线看| 一本—道久久a久久精品蜜桃钙片| 国产在线视频一区二区| 一二三四在线观看免费中文在| 国产区一区二久久| 久久久水蜜桃国产精品网| 久久久久国产精品人妻一区二区| 日本一区二区免费在线视频| kizo精华| 色婷婷久久久亚洲欧美| 国产精品秋霞免费鲁丝片| 国产伦理片在线播放av一区| 免费在线观看视频国产中文字幕亚洲 | 欧美日韩av久久| 亚洲av男天堂| 欧美日本中文国产一区发布| 一级黄色大片毛片| 亚洲第一欧美日韩一区二区三区 | 纵有疾风起免费观看全集完整版| 乱人伦中国视频| 亚洲欧美日韩高清在线视频 | av一本久久久久| 天堂中文最新版在线下载| 欧美+亚洲+日韩+国产|