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

    華南主產(chǎn)區(qū)雙季水稻物候變化及其與氣候條件的關(guān)系*

    2021-11-30 01:46:04李偉光侯美亭張京紅車秀芬陳小敏
    中國農(nóng)業(yè)氣象 2021年12期
    關(guān)鍵詞:物候晚稻早稻

    李偉光,侯美亭,張京紅,車秀芬,陳小敏

    華南主產(chǎn)區(qū)雙季水稻物候變化及其與氣候條件的關(guān)系*

    李偉光1,2,侯美亭3,張京紅1,2,車秀芬1,2,陳小敏1,2

    (1.海南省氣候中心,???570203;2.海南省南海氣象防災(zāi)減災(zāi)重點(diǎn)實(shí)驗(yàn)室,???570203;3.中國氣象局干部學(xué)院,北京 100081)

    為探索華南水稻生長發(fā)育對氣候變化的敏感性,估算氣候變化、品種更替對水稻物候期變化的貢獻(xiàn),選取1981-2013年華南雙季水稻種植區(qū)物候和氣象觀測數(shù)據(jù),采用Theil-Sen估計線性傾向率、Mann-Kendall趨勢顯著性檢驗(yàn)等方法,分析雙季水稻物候期的時空變化趨勢。結(jié)果表明:與其它地區(qū)不同,華南早稻播種期提前、晚稻成熟期推遲趨勢較小;早稻、晚稻營養(yǎng)生長期存在縮短趨勢,而生殖生長期存在延長趨勢。早稻、晚稻的營養(yǎng)生長期、生殖生長期均隨著降水量和日照時數(shù)的增加而延長,隨著氣溫上升而縮短;平均氣溫每上升1℃,生長期長度縮短0.5~4.7d。氣候變化與品種對水稻物候期影響的貢獻(xiàn)率分析表明,品種更替對早稻和晚稻生長期長度的影響大于氣候變化。研究期內(nèi)品種更替中,早稻傾向于選擇營養(yǎng)生長期和生殖生長期均長的品種;晚稻傾向于選擇營養(yǎng)生長期短、生殖生長期長的品種。

    水稻;物候;發(fā)育期;氣候變化;貢獻(xiàn)率

    植物物候是一年中植物生長階段的時間記錄,它反應(yīng)了植物對周圍環(huán)境條件變化的響應(yīng),是全球氣候變化的有力證據(jù)[1]。天然植被物候的變化已經(jīng)在全球各地觀測到,主要表現(xiàn)為春季提前、秋季推遲[2]。在農(nóng)業(yè)生產(chǎn)中,人們會通過改良品種、改進(jìn)栽培措施來延長生育期,以充分利用熱量資源,提升產(chǎn)量。因此,觀測到的物候變化是氣候變化和品種更新共同作用的結(jié)果[3]。國內(nèi)外已有眾多學(xué)者對糧食作物開展物候變化的研究,并取得大量成果。不同地區(qū)的冬小麥對氣候變化的反應(yīng)基本一致,播種期推遲[4],小麥營養(yǎng)生長期和生長季都在縮短,而生殖生長期長度呈現(xiàn)延長趨勢[5]。中國春玉米物候期提前而夏玉米和春夏播玉米物候期推遲,營養(yǎng)生長期呈縮短、生殖生長期呈延長趨勢[6]。

    水稻作為中國產(chǎn)量最大的糧食作物之一,針對其物候期變化的相關(guān)研究也已開展。Tao等發(fā)現(xiàn)中國的單季稻移栽期提前、成熟期推遲,營養(yǎng)生長期、生殖生長期及全生育期均延長;而長江中下游的雙季水稻物候期提前,營養(yǎng)生長期、全生育期縮短,生殖生長期延長[7]。Bai等發(fā)現(xiàn)早稻生育期提前、晚稻生育期延遲,單季稻和早稻的營養(yǎng)生長、生殖生長和總生育期均延長,而晚稻表現(xiàn)為縮短[8]。東北地區(qū)一季稻在2001?2012年間生育期表現(xiàn)出縮短而1992?2001年的生育期表現(xiàn)出延長現(xiàn)象[3]。Hu等研究發(fā)現(xiàn)氣候變暖使單季稻、早稻和晚稻生育期均呈現(xiàn)縮短趨勢[9]。氣象等環(huán)境因子也影響著水稻生育期。Liu等認(rèn)為氣溫升高、日照時數(shù)增加使生育期縮短,而降水增加對延長生育期有幫助[10]。水稻分布廣泛,有秈稻、粳稻兩個亞種,地區(qū)間一季水稻、雙季水稻、早稻晚稻等種植制度不同,加之氣象條件的空間分布差異、季節(jié)增溫的不對稱性[11]共同導(dǎo)致了水稻物候變化的時空不一致。

    華南地區(qū)水熱資源充沛,限制及影響水稻物候變化的環(huán)境條件與其它地區(qū)差異較大。因此,本文對華南雙季水稻主產(chǎn)區(qū)早稻、晚稻的物候變化及其與氣象條件的關(guān)系展開研究,以期發(fā)現(xiàn)雙季水稻適應(yīng)氣候變化的規(guī)律,提出品種適應(yīng)當(dāng)前氣候變化的針對性措施,進(jìn)而指導(dǎo)該地區(qū)農(nóng)業(yè)研究和生產(chǎn)。

    1 數(shù)據(jù)與方法

    1.1 農(nóng)業(yè)氣象站和數(shù)據(jù)

    水稻物候觀測數(shù)據(jù)來源于中國農(nóng)作物生長發(fā)育數(shù)據(jù)集(http://data.cma.cn/data/cdcdetail/dataCode/AGME_ AB2_CHN_TEN.html),數(shù)據(jù)時間范圍為1981-2013年。數(shù)據(jù)的觀測記錄依據(jù)《農(nóng)業(yè)氣象觀測規(guī)范》[12],水稻物候期包括播種、出苗、三葉、移栽、返青、分蘗、拔節(jié)、孕穗、抽穗、乳熟、成熟共11個。營養(yǎng)生長期(Vegetative growth period,VGP)為播種-抽穗階段,生殖生長期(Reproductive growth period,RGP)為抽穗-成熟階段,生長季(Whole growth period,WGP)為播種-成熟階段。

    氣候觀測數(shù)據(jù)來自中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)提供的中國地面氣候資料日值數(shù)據(jù)集(http://data. cma.cn/data/cdcdetail/dataCode/SURF_CLI_CHN_MUL_DAY.html),選擇日平均氣溫、降水和日照時數(shù)進(jìn)行計算。

    數(shù)據(jù)預(yù)處理過程分三步。第一步,選取有連續(xù)物候觀測的農(nóng)業(yè)氣象觀測站點(diǎn),對數(shù)據(jù)進(jìn)行質(zhì)量控制,僅選取數(shù)據(jù)有效率90%以上站點(diǎn),有23個站點(diǎn)入選,分布于廣東、廣西及海南省的海南島,詳細(xì)分布見圖1。第二步,將物候期日期轉(zhuǎn)換為年日序(Day of year,DOY)。第三步,按站點(diǎn)和年份篩選各物候期的DOY值,分別計算生長季、營養(yǎng)生長期、生殖生長期內(nèi)的光、溫、水和積溫等氣候要素數(shù)值。

    圖1 研究區(qū)域位置及站點(diǎn)分布

    注:基于國家測繪地理信息局標(biāo)準(zhǔn)地圖服務(wù)網(wǎng)站下載的審圖號為GS( 2019) 1698號的標(biāo)準(zhǔn)地圖制作,底圖無修改。

    Note:Based on the standard map No. GS (2019) 1698, downloaded from National Bureau of surveying, mapping and geographic information, and the base map was not modified.

    1.2 方法

    1.2.1 積溫計算

    積溫(Growing degree days, GDD)是指某個時間段內(nèi)日平均氣溫大于基準(zhǔn)溫度部分的累積值,具有重要的生物學(xué)意義,可以用來估計作物在生長期內(nèi)的生長發(fā)育情況[13]。計算式為

    式中,Dstart為某發(fā)育階段的起始日期,Dend是該發(fā)育階段結(jié)束的日期,TEMd是d日內(nèi)的日平均氣溫,TEMbase為水稻發(fā)育的基準(zhǔn)溫度,即發(fā)育溫度下限,設(shè)置為10℃[14]。

    1.2.2傾向率及顯著性檢驗(yàn)

    早、晚稻物候期和生長期長度的線性傾向率采用Theil-Sen方法計算。Theil-Sen估計較普通最小二乘(OLS)對數(shù)據(jù)中異常值的容錯性更高[15]。計算式為

    式中,n為分析時段的年數(shù);Pi,k和Pj,k分別表示觀測站點(diǎn)在第i年和第j年的k物候期(用年日序表示,Day of year, DOY)或生長期長度(d);Sk表示k物候期或生長期長度的傾向率(d·a-1),median為對斜率序列取中位數(shù)。Sk大于0,表示物候期或生長期長度推遲或延長,反之表示提前或縮短。

    根據(jù)歷年生長期初始日和終止日觀測值,統(tǒng)計作物生長發(fā)育過程中的平均氣溫、降水量、日照時數(shù)和有效積溫。以上氣候要素的傾向率也用式(2)計算。

    為檢驗(yàn)變化趨勢的顯著性,進(jìn)一步結(jié)合Mann-Kendall法[16]進(jìn)行分析。Mann-Kendall法是一種非參數(shù)統(tǒng)計檢驗(yàn)方法,其優(yōu)點(diǎn)是不需要測量值服從正態(tài)分布,也不要求趨勢是線性的,并且不受缺失值和異常值的影響,廣泛應(yīng)用于長時間序列數(shù)據(jù)的趨勢顯著檢驗(yàn)。計算方法為對于時間序列Xi,i=1,2,3,…,n定義標(biāo)準(zhǔn)化檢驗(yàn)統(tǒng)計量Z。

    1.2.3 水稻生長期長度對各氣候因子變化的敏感度分析

    水稻生長期長度對生長期內(nèi)平均氣溫、降水和日照時數(shù)變化的敏感度計算式為[17]

    1.2.4 氣候變化及品種對發(fā)育期長度變化的貢獻(xiàn)率計算

    觀測到的作物發(fā)育期變化是由作物品種更替、氣候變化以及其它農(nóng)事活動共同作用的結(jié)果。作物發(fā)育期長度的趨勢可以通過品種與氣候因素引發(fā)的發(fā)育期趨勢表示[9,18-19],即

    氣候變化導(dǎo)致的水稻生長期長度變化趨勢[18,20]為

    式中,Tphe-cli表示由氣候變化導(dǎo)致的作物發(fā)育期長度變化趨勢。Ttem、Tpre、Tssd分別為氣溫、降水、日照時數(shù)的變化趨勢。a、b和c為研究區(qū)域內(nèi)的所有站點(diǎn)水稻生長期長度對平均氣溫、降水和日照時數(shù)敏感度的平均值。

    作物品種引發(fā)的水稻生長期長度變化趨勢為

    式中,Tphe-GDD表示由作物品種變化導(dǎo)致的作物發(fā)育期長度變化趨勢,TGDD表示水稻積溫隨時間變化的趨勢,TEMavg表示發(fā)育期內(nèi)的平均氣溫,TEMbase表示水稻發(fā)育的基準(zhǔn)溫度(10℃)。

    氣候變化和品種更替引發(fā)的發(fā)育期長度趨勢變化的貢獻(xiàn)率計算式[18]分別為

    2 結(jié)果與分析

    2.1 雙季水稻播種期和成熟期變化趨勢

    利用1981?2013年華南雙季水稻播種及成熟期數(shù)據(jù),經(jīng)過傾向率估計及顯著性檢驗(yàn),雙季水稻播種期和成熟期的變化趨勢見圖2。由圖可見,1981?2013年華南雙季水稻種植區(qū)早稻播種期提前、晚稻成熟期推遲的趨勢較小,而早稻成熟期提前、晚稻播種期推遲的趨勢明顯。早稻播種期平均每10年提前0.3d,成熟期提前1.2d。23個站點(diǎn)中早稻播種期呈現(xiàn)提前趨勢的站點(diǎn)有9個,其中7個通過0.05水平的顯著性檢驗(yàn);推遲的站點(diǎn)有10個,其中5個通過0.05水平的顯著性檢驗(yàn)(圖2a)。早稻成熟期呈現(xiàn)提前趨勢的站點(diǎn)有14個,其中6個站點(diǎn)通過0.05水平的顯著性檢驗(yàn);推遲的站點(diǎn)有5個,沒有站點(diǎn)通過顯著性檢驗(yàn)(圖2b)。晚稻播種、成熟期變化幅度更大,播種期平均每10年推遲4.8d,收獲期推遲0.6d。晚稻播種期表現(xiàn)為推遲的站點(diǎn)有21個,其中15個站點(diǎn)通過0.05水平的顯著性檢驗(yàn);僅有1個站點(diǎn)表現(xiàn)為播種期提前(圖2c)。晚稻收獲期有8個站點(diǎn)表現(xiàn)出提前,其中3個站點(diǎn)通過0.05水平的顯著性檢驗(yàn);14個站點(diǎn)推遲,其中7個站點(diǎn)通過0.05水平的顯著性檢驗(yàn)(圖2d)。

    圖2 1981?2013年雙季水稻播種期和成熟期線性傾向率的空間分布

    注:帶〇的站點(diǎn)表示通過0.05水平的顯著性檢驗(yàn)。下同。

    Note: The site with 〇 passed 0.05 level significance test. The same as below.

    早稻成熟、晚稻播種期變化趨勢空間分布上較為一致(圖2b、圖2c),從南到北普遍出現(xiàn)提前、推遲的趨勢,相反表現(xiàn)的站點(diǎn)呈零星分布。早稻播種、晚稻成熟變化趨勢的空間一致性較差(圖2a、圖2d),表現(xiàn)出推遲、提前的站點(diǎn)零星分布,沒有規(guī)律性的聚集。

    2.2 雙季水稻生育期長度的變化趨勢

    利用雙季水稻營養(yǎng)生長期、生殖生長期及生長季的長度數(shù)據(jù),經(jīng)過線性傾向率估計及顯著性檢驗(yàn),得到1981?2013年華南雙季水稻生育期長度的變化趨勢(圖3)。由圖3可見,從早稻、晚稻各生長期長度變化趨勢的空間分布來看,早晚稻生長季普遍存在縮短的趨勢;營養(yǎng)生長期、生殖生長期分開來看,前者呈現(xiàn)出縮短而后者呈現(xiàn)出延長的趨勢。早稻全生育期平均以0.7d·10a-1的速率縮短,23個站點(diǎn)中早稻全生育期縮短的有12個站點(diǎn),其中5個通過0.05水平的顯著性檢驗(yàn);全生育期延長的有8個站點(diǎn),僅2個通過0.05水平的顯著性檢驗(yàn)。早稻營養(yǎng)生長期縮短趨勢更大,約以1.1d·10a-1的速率縮短。23個站點(diǎn)中早稻營養(yǎng)生長期縮短的有15個站點(diǎn),其中5個通過0.05水平的顯著性檢驗(yàn);延長的站點(diǎn)僅4個。早稻生殖生長期平均以0.6d·10a-1的速率延長,半數(shù)的站點(diǎn)呈現(xiàn)出延長趨勢,其中4個通過0.05水平的顯著性檢驗(yàn);9個站點(diǎn)未表現(xiàn)出上升或下降趨勢,3個站點(diǎn)縮短。

    圖3 雙季水稻生長期長度線性傾向率的空間分布

    Note: VGP is vegetative growth period, RGP is reproductive growth period, WGP is whole growth period. The same as below.

    晚稻全生育期、營養(yǎng)生長期長度變化的趨勢更為一致,普遍表現(xiàn)出縮短的趨勢,分別以4.5d·10a-1、4.4d·10a-1的速率縮短。生殖生長期表現(xiàn)出小幅度的延長趨勢,平均以0.3d·10a-1的速率延長。從變化趨勢的空間一致性來看,趨勢也較為統(tǒng)一。全生育期、營養(yǎng)生長期縮短的站點(diǎn)分別有16、19個,分別有12個、11個站點(diǎn)通過0.05水平的顯著性檢驗(yàn),但表現(xiàn)為上升的站點(diǎn)均未通過顯著性檢驗(yàn)。超過2/3的站點(diǎn)表現(xiàn)出生殖生長期延長,但延長幅度較小。

    2.3 雙季水稻生長季內(nèi)主要?dú)夂蛞蜃拥淖兓厔?/h3>

    利用1981?2013年雙季水稻生長季內(nèi)的平均氣溫、降水量、日照時數(shù)及積溫數(shù)據(jù),經(jīng)過線性傾向率及顯著性檢驗(yàn),得到華南雙季水稻生長季內(nèi)主要?dú)夂蛞蜃拥淖兓厔?。由圖4a可見,華南地區(qū)早稻生長季內(nèi)降水量變化趨勢不統(tǒng)一,表現(xiàn)出上升和下降趨勢的站點(diǎn)分別為12個和11個。晚稻生長季內(nèi)大部分站點(diǎn)表現(xiàn)出降水減少的趨勢(圖4b),23個站點(diǎn)中19個站點(diǎn)下降,總體平均以31.6mm·10a-1的速率下降。降水量減少的原因可能是因?yàn)橥淼旧L季縮短顯著。早稻生長季內(nèi)多數(shù)站點(diǎn)(圖4c)日照時數(shù)表現(xiàn)出增加的趨勢,其中3個站點(diǎn)通過0.05水平的顯著性檢驗(yàn);總體以16.1h·10a-1的速率增加。而晚稻生長季內(nèi)有14個站點(diǎn)表現(xiàn)出日照時數(shù)下降(圖4d),其中4個站點(diǎn)通過0.05水平的顯著性檢驗(yàn),平均以13.7h·10a-1的速率減少。

    早稻、晚稻生育期內(nèi)的平均氣溫上升趨勢均較明顯(圖4e、圖4f),上升速率分別為0.2℃·10a-1和0.1℃·10a-1。早稻生長季內(nèi)有15個站點(diǎn)平均氣溫上升,上升速率超過0.3℃的站點(diǎn)有11個,9個站點(diǎn)上升趨勢通過顯著性檢驗(yàn);呈現(xiàn)下降的站點(diǎn)有8個,呈零星分布。晚稻生長季內(nèi)平均氣溫上升速率略小,有14個站點(diǎn)平均氣溫上升,其中6個站點(diǎn)上升趨勢通過顯著性檢驗(yàn)。早稻生長季的有效積溫普遍表現(xiàn)出增加趨勢,所有站點(diǎn)平均以42℃·10a-1的速率增加,12個站點(diǎn)的增加趨勢通過顯著性檢驗(yàn),僅有3個站點(diǎn)表現(xiàn)出減少趨勢。與早稻不同,晚稻生長季內(nèi)多半站點(diǎn)有效積溫表現(xiàn)出輕微下降趨勢(22℃·10a-1),23個站點(diǎn)中有15個站點(diǎn)表現(xiàn)為下降,7個站點(diǎn)下降趨勢通過顯著性檢驗(yàn)。下降站點(diǎn)在兩廣地區(qū)分布相對集中,而有效積溫增加站點(diǎn)零星分布。有效積溫減少的原因可能也是晚稻生長季縮短引起的。

    2.4 雙季水稻生長期長度對氣候因子變化的敏感性

    利用雙季水稻營養(yǎng)生長期、生殖生長期長度與相應(yīng)時段內(nèi)的降水量、日照時數(shù)、平均氣溫數(shù)據(jù),建立線性回歸模型,得到1981?2013年華南雙季水稻種植區(qū)生長期長度對相應(yīng)時段內(nèi)平均氣溫、降水量和日照時數(shù)變化的平均敏感度,結(jié)果見表1。由表可見,構(gòu)建的水稻生長期長度對平均氣溫、降水量和日照時數(shù)變化的敏感度方程平均決定系數(shù)R2較大,達(dá)到0.380~0.762(P<0.05),說明構(gòu)建的方程具有較高的可靠性。各生長期長度對平均氣溫的敏感度為負(fù)值,說明早稻、晚稻的營養(yǎng)生長和生殖生長期長度均隨氣溫升高而縮短。具體數(shù)值顯示,平均氣溫每升高1℃,早稻營養(yǎng)生長期縮短4.7d,生殖生長期縮短2.1d;晚稻營養(yǎng)生長期縮短4.3d,生殖生長期縮短0.5d。與晚稻相比,早稻對平均氣溫更加敏感。早稻和晚稻營養(yǎng)生長期對平均氣溫的敏感度顯著高于生殖生長期。

    各生長期長度對日照時數(shù)變化的敏感度均為正值,說明各生長期長度均隨日照時數(shù)增加而延長。早稻、晚稻不同生長期長度對日照的敏感性相近,日照時數(shù)每增加100h,生育階段長度延長4.7~7.8d。

    各生長期長度對降水量變化的敏感度數(shù)值均為正值,說明早稻、晚稻的營養(yǎng)生長和生殖生長期長度均隨降水增多而延長。營養(yǎng)生長期內(nèi)降水量每增加100mm生長期延長1.1~1.7d,生殖生長期內(nèi)降水量每增加100mm生長期延長0.9~2.1d。

    2.5 品種更替及氣候變化對生長期長度變化的貢獻(xiàn)率

    積溫是植物內(nèi)在生命活動周期的度量,積溫的變化趨勢能夠在一定程度上反應(yīng)出品種的變遷。統(tǒng)計1981?2013年水稻各生育期內(nèi)積溫的平均傾向率表明(表2),華南雙季水稻主產(chǎn)區(qū)的早稻營養(yǎng)生長期和生殖生長期以及晚稻生殖生長期的積溫均

    圖4 雙季水稻生長季內(nèi)氣象要素線性傾向率的空間分布

    表1 雙季水稻生長期長度對各氣候因子變化的平均敏感度及顯著性

    表2 雙季水稻積溫變化趨勢折算生長期長度趨勢(CT)

    表現(xiàn)出上升趨勢。早稻營養(yǎng)生長期、生殖生長期積溫分別以25.2、21.8℃·d·10a-1的速率增加。而晚稻營養(yǎng)生長期積溫表現(xiàn)為明顯減少的趨勢,速率為43.5d·10a-1;生殖生育期呈現(xiàn)增加趨勢,速率為26.6℃·d·10a-1,增加幅度與早稻生殖生長期相近。根據(jù)積溫趨勢與平均氣溫折算生長期趨勢,早稻營養(yǎng)生長期、生殖生長期、晚稻生殖生長期均表現(xiàn)延長趨勢,延長速率分別為1.9、1.2、2.0℃·d·10a-1,晚稻營養(yǎng)生長期表現(xiàn)為以2.5℃·d·10a-1速率縮短。

    對1981?2013年23個站點(diǎn)氣候變化及品種更替對水稻生長期長度變化的平均貢獻(xiàn)率分析發(fā)現(xiàn)(圖5),氣候變化對華南雙季水稻生長期長度變化的貢獻(xiàn)較??;品種對生長期長度變化的貢獻(xiàn)較大且一般為正效應(yīng),即品種更替延長水稻生長期長度;僅在晚稻營養(yǎng)生長期兩者表現(xiàn)出負(fù)效應(yīng)。在早稻和晚稻的生殖生長期,品種因素均表現(xiàn)出超過50%的延長貢獻(xiàn)率,而氣候變化貢獻(xiàn)率均不足5%。營養(yǎng)生長期,品種更替對早稻與晚稻生長期長度的作用相反,對早稻表現(xiàn)為延長,對晚稻表現(xiàn)為縮短。氣候變化對生長期影響最大的時期為晚稻營養(yǎng)生長期,貢獻(xiàn)率達(dá)到45%,表現(xiàn)為縮短生長期長度。

    圖5 氣候變化及品種因素對水稻生長期長度變化的貢獻(xiàn)率

    3 結(jié)論與討論

    3.1 結(jié)論

    (1)1981?2013年華南雙季水稻種植區(qū)早稻播種期、成熟期分別以平均0.3、1.2 d·10a-1的速率提前,晚稻播種期、成熟期分別以4.8、0.6 d·10a-1的速率推遲。早稻成熟期、晚稻播種期變化趨勢的空間分布相對一致,區(qū)域內(nèi)大部分站點(diǎn)整體提前。早稻和晚稻的營養(yǎng)生長期均存在縮短,生殖生長期均存在延長的趨勢,總生育期也存在縮短的趨勢。

    (2)早稻生長季內(nèi)降水量變化趨勢不明顯,而晚稻生長季內(nèi)降水量普遍表現(xiàn)出減少的趨勢。早稻生長季內(nèi)多數(shù)地區(qū)日照時數(shù)增加而晚稻生長季內(nèi)大部分區(qū)域日照時數(shù)下降。早稻、晚稻生長季內(nèi)平均氣溫普遍呈上升趨勢。早稻生長季有效積溫表現(xiàn)出增加趨勢,而晚稻生長季有效積溫下降。

    (3)早稻、晚稻的營養(yǎng)生長期和生殖生長期長度隨降水量和日照時數(shù)增加而延長,隨平均氣溫的上升而縮短,平均氣溫每上升1℃,生長期縮短0.5~4.7d。通過積溫反應(yīng)的品種變化引起的生長期長度變化來看,早稻營養(yǎng)生長期、生殖生長期及晚稻生殖生長期存在較大的延長趨勢(1.2~2.0d·10a-1),晚稻營養(yǎng)生長期長度趨勢為以2.5d·10a-1速率縮短。貢獻(xiàn)率分析表明品種對生長期長度變化趨勢的影響大于氣候變化。品種更替過程中,早晚稻均傾向于選擇生殖生長期長的品種。但在晚稻營養(yǎng)生長期,品種與氣候變化均表現(xiàn)為負(fù)效應(yīng),且貢獻(xiàn)了90%以上的作用。

    3.2 討論

    1981?2013年華南雙季水稻種植區(qū),早稻播種期提前、晚稻成熟期推遲的趨勢較小,這與該研究區(qū)域以北地區(qū)水稻物候變化趨勢存在明顯不同[7]。其原因可能是華南地區(qū)水熱資源較北方水稻種植區(qū)豐富,水稻界限溫度上的日數(shù)已經(jīng)足夠雙季水稻發(fā)育。水稻種植上利用氣候變化帶來的增溫需求不強(qiáng)烈,所以早稻播種期提前、晚稻成熟期推遲趨勢較小。

    華南雙季水稻生育期長度變化趨勢與其它地區(qū)[7-8,10]基本一致,表現(xiàn)為營養(yǎng)生長期和生長季普遍縮短,而生殖生長期延長。從作物產(chǎn)量形成的角度看,生殖生長期越長,接收的光合輻射越多,積累的干物質(zhì)越多,產(chǎn)量越高[21-22]。所以延長作物生殖生長期是育種家追求的重要目標(biāo),這一點(diǎn)在華南早稻和晚稻中均有確實(shí)證據(jù)。從驅(qū)動水稻物候變化因素的貢獻(xiàn)率看,氣候變化對華南雙季水稻的影響遠(yuǎn)小于品種更替因素,這可能是因?yàn)槿A南水熱充足,氣候條件的波動仍然在水稻適宜范圍內(nèi)。晚稻營養(yǎng)生長期,品種更替及氣候變化一致表現(xiàn)為縮短貢獻(xiàn),產(chǎn)生這個現(xiàn)象的原因可能是晚稻營養(yǎng)生長期是一年中氣溫最高的時段,環(huán)境氣溫可能已經(jīng)接近水稻生長的最適溫度[23]。在氣溫升高的背景下,農(nóng)戶主動選擇營養(yǎng)生長期短的品種,作物被動縮短生育期來共同應(yīng)對高溫脅迫。

    農(nóng)作物的物候變化是氣候變化和農(nóng)事管理共同作用的結(jié)果,區(qū)分二者貢獻(xiàn)對客觀評價氣候變化帶來的影響具有重要意義。通常采用的方法有統(tǒng)計[24]和作物模型[25],但兩種方法還沒有壓倒性的優(yōu)勢。作物模型的方法可以模擬作物、環(huán)境條件及管理措施之間的相互作用,它需要更多的信息進(jìn)行模型校準(zhǔn)和驗(yàn)證,而它的跨區(qū)域推廣性較差[26]。統(tǒng)計方法只是建立了作物表現(xiàn)與氣象因素之間的數(shù)理關(guān)系,而這種關(guān)系并不一定是因果關(guān)系。本研究早稻、晚稻不同發(fā)育階段對氣象因子的敏感系數(shù)大小相近,這也在一定程度上肯定了統(tǒng)計方法結(jié)果的可靠性。未來可以結(jié)合田間試驗(yàn),進(jìn)一步探索其中的生理基礎(chǔ),分析水稻對氣候變化的響應(yīng)機(jī)制,增強(qiáng)水稻適應(yīng)氣候變化的能力。

    [1] Ge Q,Wang H,Rutishauser T,et al.Phenological response to climate change in China:a meta-analysis[J].Glob Chang Biology,2015,21(1):265-274.

    [2] Chen X,Wang L,Niu Z,et al.The effects of projected climate change and extreme climate on maize and rice in the Yangtze River Basin,China[J].Agricultural and Forest Meteorology,2020:282-283.

    [3] Pachauri R K,Allen M R,Barros V R,et al.Climate change 2014:Synthesis report.Contribution of Working Groups I,II and III to the Fifth assessment report of the intergovernmental panel on climate change[M]. Geneva,Switzerland:IPCC, 2014.

    [4] 劉玉潔,陳巧敏,葛全勝,等.氣候變化背景下1981-2010中國小麥物候變化時空分異[J].中國科學(xué):地球科學(xué),2018, 48(7):888-898.

    Liu Y J,Chen Q M,Ge Q S,et al.Spatiotemporal differentiation of changes in wheat phenology in China under climate change from 1981 to 2010[J].Science China:Earth Sciences,2018,48(7):888-898.(in Chinese)

    [5] 張媛鈴,郭曉磊,王娜,等.基于決策系統(tǒng)模擬不同降水年型旱作冬小麥的最佳播期[J].中國農(nóng)業(yè)氣象,2021,42(6): 475-485.

    Zhang Y L,Guo X L,Wang N,et al.Analysis on the optimal sowing date of dry-land winter wheat under different precipitation pattern based on wheat decision system[J]. Chinese Journal of Agrometeorology,2021,42(6):475-485. (in Chinese)

    [6] 秦雅,劉玉潔,葛全勝.氣候變化背景下1981-2010年中國玉米物候變化時空分異[J].地理學(xué)報,2018,73(5):906-916.

    Qin Y, Liu Y J, Ge Q S. Spatiotemporal variations in maize phenology of China under climate change from 1981 to 2010[J].Acta Geographica Sinica,2018,73(5):906-916.(in Chinese)

    [7] Tao F L,Zhao Z,Shi W J, et al.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[J].Global change biology,2013,19(10): 3200-3209.

    [8] Bai H,Xiao D.Spatiotemporal changes of rice phenology in China during 1981-2010[J].Theoretical and Applied Climatology,2020,140(3):1483-1494.

    [9] Hu X,Huang Y,Sun W,et al.Shifts in cultivar and planting date have regulated rice growth duration under climate warming in China since the early 1980s[J].Agricultural and Forest Meteorology,2017,247:34-41.

    [10] Liu Y,Zhou W,Ge Q.Spatiotemporal changes of rice phenology in China under climate change from 1981 to 2010[J].Climatic Change,2019,157(2):261-277.

    [11] Cai Q,Liu Y,Wang Y,et al.Recent warming evidence inferred from a tree-ring-based winter-half year minimum temperature reconstruction in northwestern Yichang,South Central China,and its relation to the large-scale circulation anomalies[J].International Journal of Biometeorology,2016, 60(12):1885-1896.

    [12] 國家氣象局.農(nóng)業(yè)氣象觀測規(guī)范(水稻)[M].北京:氣象出版社,1993.

    National Meteorological Administration.Standard for agrometeorological observation(rice)[M].Beijing:China Meteorology Press,1993.(in Chinese)

    [13] Xiao L,Liu L,Asseng S,et al.Estimating spring frost and its impact on yield across winter wheat in China[J]. Agricultural and Forest Meteorology,2018,260-261:154- 164.

    [14] 黃維,吳炫柯,劉永裕,等.氣候變化對廣西雙季稻種植布局的影響[J].中國農(nóng)業(yè)氣象,2020,41(9):539-551.

    Huang W,Wu X K,Liu Y Y,et al.Impact of climate change on layout of double cropping rice in Guangxi[J].Chinese Journal of Agrometeorology,2020,41(9):539-551.(in Chinese)

    [15] Fernandes R,Leblanc S G.Parametric(modified least squares) and non-parametric(Theil–Sen) linear regressions for predicting biophysical parameters in the presence of measurement errors[J].Remote Sensing of Environment, 2005,95(3):149-153.

    [16] 黃家俊,張強(qiáng),張生,等.基于信息熵的新疆降水時空變異特征研究[J].生態(tài)學(xué)報,2017,37(13):4444-4455.

    Huang J J,Zhang Q,Zhang S,et al.Information entropy- based analysis of spatial and temporal variation in precipitation in Xinjiang[J].Acta Ecologica Sinica,2017, 37(13):4444-4455.(in Chinese)

    [17] Liu Y J, Qin Y,Ge Q S.Spatiotemporal differentiation of changes in maize phenology in China from 1981 to 2010[J].Journal of Geographical Sciences,2019,29(3):351- 362.

    [18] He L,Jin N,Yu Q.Impacts of climate change and crop management practices on soybean phenology changes in China[J].Science of the Total Environment,2020,707: 135638.

    [19] Verón S R,De A D,Lobell D B.Impacts of precipitation and temperature on crop yields in the Pampas[J].Climatic Change,2015,130(2):235-245.

    [20] Liu Y,Chen Q,Ge Q,et al.Modelling the impacts of climate change and crop management on phenological trends of spring and winter wheat in China[J].Agricultural and Forest Meteorology,2018,248:518-526.

    [21] Gautam P,Lal B,Nayak A K,et al.Inter-relationship between intercepted radiation and rice yield influenced by transplanting time, method, and variety[J].International Journal of Biometeorology,2019,63(3):337-349.

    [22] 吳俊,鄧啟云,袁定陽,等.超級雜交稻研究進(jìn)展[J].科學(xué)通報,2016,61(35):3787-3796.

    Wu J,Deng Q Y,Yuan D Y,et al.Progress of super hybrid rice research in China[J].Chinese Science Bulletin,2016, 61(35):3787-3796.(in Chinese)

    [23] 宋有金,吳超,李子煜,等.水稻產(chǎn)量對生殖生長階段不同時期高溫的響應(yīng)差異[J].中國水稻科學(xué),2021,35(2):177-186.

    Song Y J,Wu C,Li Z Y,et al.Differential responses of grain yields to high temperature in different stages of reproductive growth in rice[J].Chinese Journal of Rice Science,2021,35(2):177-186.(in Chinese)

    [24] Zhao H,Fu Y H,Wang X,et al.Timing of rice maturity in China is affected more by transplanting date than by climate change[J].Agricultural and Forest Meteorology, 2016,216:215-220.

    [25] He L,Asseng S,Zhao G,et al.Impacts of recent climate warming,cultivar changes,and crop management on winter wheat phenology across the Loess Plateau of China[J]. Agricultural and Forest Meteorology,2015,200:135-143.

    [26] Zhang Q,Zhang W,Li T,et al.Accuracy and uncertainty analysis of staple food crop modelling by the process-based Agro-C model[J].International Journal of Biometeorology, 2021,65:587-599.

    Spatiotemporal Changes of Double Cropping Rice Phenology in South China and Relationship with Climatic Conditions

    LI Wei-guang1,2, HOU Mei-ting3, ZHANG Jing-hong1,2, CHE Xiu-fen1,2, CHEN Xiao-min1,2

    (1.Climate Center of Hainan Province, Haikou 570203, China; 2. Key Laboratory of Meteorological Disaster Prevention and Reduction in South China Sea, Haikou 570203; 3. Cadre College of China Meteorological Administration, Beijing 100081)

    In order to explore the sensitivity of rice growth and development to climate change, and to estimate the contribution of climate change and varieties change to the change of rice phenology in South China, the phenological and meteorological observation data of double cropping rice planting in South China from 1981 to 2013 were selected for analysis. Theil-Sen estimation linear tendency rate and Mann-Kendall trend significance test were used to analyze the change trend and spatial distribution characteristics of double cropping rice phenology. Different from other regions, the sowing date of early rice in South China is ahead of schedule and the mature period of late rice is delayed slightly. The vegetative growth period of early and late rice was shortened, while the reproductive growth period was prolonged. The length of vegetative and reproductive growth periods of early rice and late rice increased with the increase of precipitation and sunshine hours, and decreased with the increase of temperature; the length of growth period was shortened by 0.5?4.7 days when the average temperature increased by 1℃. The analysis of the contribution rate of climate change and varieties to the impact of rice phenology showed that the influence of varieties change on the long-term length of early and late rice was greater than that of climate change. In the past 30 years, early rice tended to choose the varieties with long vegetative and reproductive growth period, late rice tends to choose varieties with short vegetative growth period and long reproductive growth period.

    Rice;Phenology;Growth phases;Climate change;Contribution

    10.3969/j.issn.1000-6362.2021.12.004

    李偉光,侯美亭,張京紅,等.華南主產(chǎn)區(qū)雙季水稻物候變化及其與氣候條件的關(guān)系[J].中國農(nóng)業(yè)氣象,2021,42(12):1020-1030

    2021?04?27

    中國氣象局創(chuàng)新發(fā)展專項(CXFZ2021J070);氣象軟科學(xué)項目(2021ZZXM16)

    李偉光,高級工程師,研究方向?yàn)闅夂蜃兓皻庀鬄?zāi)害對農(nóng)業(yè)的影響,E-mail:163great@163.com

    猜你喜歡
    物候晚稻早稻
    早稻迎豐收
    海南橡膠林生態(tài)系統(tǒng)凈碳交換物候特征
    2021年全國早稻總產(chǎn)量2802萬噸增長2.7%
    2021年全國早稻總產(chǎn)量2802萬噸 比2020年增產(chǎn)72.3萬噸
    宋代早稻若干問題探討
    ‘灰棗’及其芽變品系的物候和生育特性研究
    早熟蓮藕-晚稻栽培技術(shù)
    5種忍冬科植物物候期觀察和比較
    約旦野生二棱大麥在川西高原的物候期和農(nóng)藝性狀分析
    2014年晚稻新品種生產(chǎn)試驗(yàn)
    亚洲国产欧美一区二区综合| 欧美 日韩 精品 国产| 欧美人与性动交α欧美软件| 美国免费a级毛片| 国精品久久久久久国模美| 91老司机精品| 青青草视频在线视频观看| 国产成人av激情在线播放| 女人精品久久久久毛片| 欧美国产精品一级二级三级| 午夜免费鲁丝| 丁香六月欧美| 老司机影院成人| 麻豆精品久久久久久蜜桃| 自线自在国产av| 欧美精品高潮呻吟av久久| 新久久久久国产一级毛片| 国产xxxxx性猛交| 亚洲欧美日韩另类电影网站| 欧美乱码精品一区二区三区| 9热在线视频观看99| 另类精品久久| www.自偷自拍.com| 免费黄频网站在线观看国产| 三上悠亚av全集在线观看| 中文字幕制服av| 久久久久久久大尺度免费视频| 纯流量卡能插随身wifi吗| 免费黄网站久久成人精品| 国产一级毛片在线| 亚洲精品国产av蜜桃| 精品亚洲乱码少妇综合久久| 丝瓜视频免费看黄片| 无遮挡黄片免费观看| 在线观看免费视频网站a站| 亚洲人成网站在线观看播放| 午夜福利在线免费观看网站| 久久久精品国产亚洲av高清涩受| 久久 成人 亚洲| 一区二区三区激情视频| 男女床上黄色一级片免费看| 亚洲人成电影观看| 久久影院123| 曰老女人黄片| 欧美黑人精品巨大| 欧美xxⅹ黑人| 波多野结衣一区麻豆| 亚洲美女视频黄频| 成年女人毛片免费观看观看9 | 亚洲av电影在线观看一区二区三区| 日本猛色少妇xxxxx猛交久久| 中国国产av一级| 久久久久久久精品精品| 99九九在线精品视频| 色婷婷久久久亚洲欧美| 国产精品一二三区在线看| 在线亚洲精品国产二区图片欧美| 国产亚洲精品第一综合不卡| 国产免费又黄又爽又色| 捣出白浆h1v1| 色网站视频免费| 久久免费观看电影| 丝袜人妻中文字幕| 国产精品秋霞免费鲁丝片| 日韩不卡一区二区三区视频在线| 老鸭窝网址在线观看| 精品一区二区三区四区五区乱码 | 哪个播放器可以免费观看大片| 亚洲美女黄色视频免费看| 亚洲av欧美aⅴ国产| 亚洲精品视频女| 精品国产乱码久久久久久小说| 午夜久久久在线观看| 啦啦啦在线免费观看视频4| 国产老妇伦熟女老妇高清| 国产精品久久久久久精品古装| 啦啦啦啦在线视频资源| 久久久久久久大尺度免费视频| 午夜影院在线不卡| 亚洲,欧美精品.| 香蕉丝袜av| av网站免费在线观看视频| 欧美精品一区二区免费开放| 精品国产一区二区久久| 亚洲av福利一区| 我的亚洲天堂| 欧美日韩国产mv在线观看视频| 亚洲精品成人av观看孕妇| 另类精品久久| 亚洲视频免费观看视频| 成人亚洲欧美一区二区av| 麻豆av在线久日| 久久精品aⅴ一区二区三区四区| 精品一区二区免费观看| 国产视频首页在线观看| 亚洲自偷自拍图片 自拍| 18禁动态无遮挡网站| 欧美老熟妇乱子伦牲交| 美女大奶头黄色视频| 成人国产麻豆网| 国产亚洲午夜精品一区二区久久| av在线播放精品| 侵犯人妻中文字幕一二三四区| 免费黄频网站在线观看国产| 亚洲精品美女久久av网站| 亚洲人成电影观看| 啦啦啦在线观看免费高清www| 黄色怎么调成土黄色| 亚洲欧美日韩另类电影网站| av线在线观看网站| 99re6热这里在线精品视频| 亚洲七黄色美女视频| 中文欧美无线码| 99精品久久久久人妻精品| 精品一区二区免费观看| 亚洲欧洲精品一区二区精品久久久 | 精品人妻熟女毛片av久久网站| 成人亚洲精品一区在线观看| 久久精品亚洲av国产电影网| 一二三四中文在线观看免费高清| 日本vs欧美在线观看视频| 少妇猛男粗大的猛烈进出视频| av线在线观看网站| 午夜福利网站1000一区二区三区| 久久久精品94久久精品| 中文字幕亚洲精品专区| 美女脱内裤让男人舔精品视频| 国产成人91sexporn| netflix在线观看网站| 97在线人人人人妻| 一级毛片黄色毛片免费观看视频| 久久久久精品国产欧美久久久 | 国产精品亚洲av一区麻豆 | 免费人妻精品一区二区三区视频| 国产av国产精品国产| 成年av动漫网址| 波野结衣二区三区在线| 各种免费的搞黄视频| 男男h啪啪无遮挡| 丰满饥渴人妻一区二区三| av网站在线播放免费| 在线观看三级黄色| 久久久久国产一级毛片高清牌| 欧美日韩一级在线毛片| 99久久精品国产亚洲精品| 亚洲国产中文字幕在线视频| 国产免费一区二区三区四区乱码| 亚洲三区欧美一区| 国产精品麻豆人妻色哟哟久久| 九草在线视频观看| 久久免费观看电影| 91精品伊人久久大香线蕉| 老司机靠b影院| 19禁男女啪啪无遮挡网站| 亚洲欧洲国产日韩| 最近手机中文字幕大全| 日韩视频在线欧美| 如日韩欧美国产精品一区二区三区| av在线播放精品| 女人高潮潮喷娇喘18禁视频| 母亲3免费完整高清在线观看| 国产人伦9x9x在线观看| 欧美日韩亚洲高清精品| 亚洲精品aⅴ在线观看| 搡老乐熟女国产| 999精品在线视频| 亚洲精品美女久久av网站| 菩萨蛮人人尽说江南好唐韦庄| 女人爽到高潮嗷嗷叫在线视频| 日本欧美视频一区| 最新的欧美精品一区二区| 精品一区二区三区av网在线观看 | 亚洲,欧美,日韩| 欧美日韩福利视频一区二区| 日韩欧美精品免费久久| 秋霞伦理黄片| av福利片在线| 日韩欧美精品免费久久| 亚洲,欧美精品.| 视频区图区小说| 美国免费a级毛片| 久久精品aⅴ一区二区三区四区| 久久性视频一级片| 美国免费a级毛片| 精品人妻在线不人妻| 一区在线观看完整版| 人妻人人澡人人爽人人| 人人妻人人爽人人添夜夜欢视频| 国产97色在线日韩免费| 丰满饥渴人妻一区二区三| 成人亚洲精品一区在线观看| 国产人伦9x9x在线观看| 如何舔出高潮| 国产成人精品久久久久久| 国产视频首页在线观看| 国产野战对白在线观看| 免费人妻精品一区二区三区视频| 老司机深夜福利视频在线观看 | 国产在线一区二区三区精| 满18在线观看网站| 国产精品女同一区二区软件| 精品一区二区三卡| 99热国产这里只有精品6| av不卡在线播放| 国产伦理片在线播放av一区| 只有这里有精品99| 成人三级做爰电影| 一区二区三区精品91| e午夜精品久久久久久久| 自拍欧美九色日韩亚洲蝌蚪91| netflix在线观看网站| 伊人久久国产一区二区| 久久天堂一区二区三区四区| 黄频高清免费视频| av福利片在线| 两个人看的免费小视频| 狂野欧美激情性bbbbbb| 男人添女人高潮全过程视频| 免费观看av网站的网址| 日韩精品有码人妻一区| 精品一区二区三区av网在线观看 | 男男h啪啪无遮挡| 韩国精品一区二区三区| 亚洲人成77777在线视频| 成人国产麻豆网| 美女福利国产在线| 亚洲av在线观看美女高潮| 夫妻性生交免费视频一级片| 国产欧美日韩综合在线一区二区| 人妻人人澡人人爽人人| 精品少妇黑人巨大在线播放| 亚洲欧洲日产国产| 日本91视频免费播放| 色精品久久人妻99蜜桃| 老司机深夜福利视频在线观看 | 国产xxxxx性猛交| 国产精品无大码| 亚洲四区av| 欧美日本中文国产一区发布| 国产片内射在线| 国产爽快片一区二区三区| 成人毛片60女人毛片免费| 免费观看性生交大片5| 男女床上黄色一级片免费看| 一边亲一边摸免费视频| 国产乱人偷精品视频| 亚洲第一区二区三区不卡| 一本色道久久久久久精品综合| 国产精品久久久久成人av| av线在线观看网站| 大香蕉久久网| 欧美日韩成人在线一区二区| 亚洲国产欧美在线一区| 丰满少妇做爰视频| 亚洲av成人不卡在线观看播放网 | 看免费av毛片| 亚洲国产成人一精品久久久| 啦啦啦 在线观看视频| 日本av免费视频播放| 色视频在线一区二区三区| 美女视频免费永久观看网站| 亚洲国产日韩一区二区| 在线观看免费午夜福利视频| 久久 成人 亚洲| 乱人伦中国视频| 国产精品av久久久久免费| 国产精品亚洲av一区麻豆 | 你懂的网址亚洲精品在线观看| 国产精品亚洲av一区麻豆 | 纵有疾风起免费观看全集完整版| 国产亚洲最大av| videos熟女内射| 一区二区三区精品91| 超碰成人久久| 国产av精品麻豆| 久久影院123| 丁香六月天网| 国产熟女午夜一区二区三区| 亚洲精品第二区| 欧美日韩亚洲综合一区二区三区_| 久久久久久久久免费视频了| 亚洲精品久久成人aⅴ小说| 99香蕉大伊视频| 久久久久国产精品人妻一区二区| 香蕉丝袜av| 国产99久久九九免费精品| 秋霞伦理黄片| 色综合欧美亚洲国产小说| 女人被躁到高潮嗷嗷叫费观| 中文欧美无线码| 一边摸一边抽搐一进一出视频| 色94色欧美一区二区| 青春草亚洲视频在线观看| 国产免费又黄又爽又色| 成人国产av品久久久| 成年美女黄网站色视频大全免费| 9191精品国产免费久久| 捣出白浆h1v1| 色视频在线一区二区三区| 午夜免费观看性视频| 亚洲国产av新网站| 国产爽快片一区二区三区| 免费高清在线观看视频在线观看| 国语对白做爰xxxⅹ性视频网站| 国产精品二区激情视频| 免费av中文字幕在线| 免费黄频网站在线观看国产| 免费不卡黄色视频| 熟妇人妻不卡中文字幕| av免费观看日本| 777久久人妻少妇嫩草av网站| 黑丝袜美女国产一区| 久久青草综合色| 少妇被粗大的猛进出69影院| 日韩中文字幕视频在线看片| 无遮挡黄片免费观看| 一区在线观看完整版| 亚洲欧美日韩另类电影网站| 色精品久久人妻99蜜桃| 一级毛片我不卡| 成人午夜精彩视频在线观看| 婷婷色麻豆天堂久久| 国产精品女同一区二区软件| 高清在线视频一区二区三区| av网站免费在线观看视频| 啦啦啦在线免费观看视频4| 搡老岳熟女国产| 日本av免费视频播放| 大片免费播放器 马上看| 天天影视国产精品| 欧美精品人与动牲交sv欧美| 亚洲欧洲日产国产| 99久国产av精品国产电影| 青春草视频在线免费观看| 又大又爽又粗| 日本爱情动作片www.在线观看| 日日啪夜夜爽| 亚洲欧美精品自产自拍| 又大又黄又爽视频免费| 亚洲精品国产av成人精品| 久久久欧美国产精品| 人人妻人人澡人人看| 国产精品久久久久成人av| 亚洲视频免费观看视频| 在线观看免费午夜福利视频| 国产精品熟女久久久久浪| 人人妻人人爽人人添夜夜欢视频| e午夜精品久久久久久久| 亚洲欧美日韩另类电影网站| 精品福利永久在线观看| 亚洲精品日本国产第一区| 欧美日韩一区二区视频在线观看视频在线| 午夜福利影视在线免费观看| 免费高清在线观看日韩| 天天躁日日躁夜夜躁夜夜| 久久久久久久久久久免费av| 热re99久久国产66热| 欧美日本中文国产一区发布| 高清不卡的av网站| 色94色欧美一区二区| 一二三四中文在线观看免费高清| www.av在线官网国产| 久久久久久久精品精品| 久久久久精品人妻al黑| 人人妻人人添人人爽欧美一区卜| 亚洲精品美女久久av网站| 我的亚洲天堂| 国产欧美亚洲国产| 精品一区二区三区av网在线观看 | 成人国语在线视频| 午夜福利视频在线观看免费| 日韩中文字幕视频在线看片| 男人添女人高潮全过程视频| 欧美中文综合在线视频| 欧美变态另类bdsm刘玥| 男女高潮啪啪啪动态图| 免费看不卡的av| 夫妻午夜视频| 国产精品久久久人人做人人爽| √禁漫天堂资源中文www| 精品一品国产午夜福利视频| 国产男女内射视频| 国产又色又爽无遮挡免| 久久综合国产亚洲精品| 久久久久精品人妻al黑| 卡戴珊不雅视频在线播放| 国产成人精品福利久久| 在线观看免费高清a一片| 精品酒店卫生间| 久久精品久久精品一区二区三区| 国产成人精品福利久久| 9191精品国产免费久久| 中文字幕人妻丝袜制服| 国产精品久久久久成人av| 黄频高清免费视频| 9热在线视频观看99| av网站免费在线观看视频| 97人妻天天添夜夜摸| 国产精品久久久久久人妻精品电影 | 在线观看免费视频网站a站| 国产av码专区亚洲av| 一区二区三区乱码不卡18| 亚洲成人av在线免费| 丰满少妇做爰视频| 男女免费视频国产| 亚洲欧美中文字幕日韩二区| 这个男人来自地球电影免费观看 | 成人手机av| 婷婷色av中文字幕| 日本一区二区免费在线视频| 国产深夜福利视频在线观看| 操美女的视频在线观看| 欧美精品高潮呻吟av久久| 亚洲成国产人片在线观看| 国产一区有黄有色的免费视频| 在线观看三级黄色| 国产精品久久久人人做人人爽| 国产人伦9x9x在线观看| 亚洲精品美女久久久久99蜜臀 | 女人被躁到高潮嗷嗷叫费观| 国产精品秋霞免费鲁丝片| 高清在线视频一区二区三区| 一本—道久久a久久精品蜜桃钙片| 一级爰片在线观看| 国产精品三级大全| 亚洲男人天堂网一区| 久久久久精品国产欧美久久久 | 一边摸一边做爽爽视频免费| 男人操女人黄网站| 久久精品国产a三级三级三级| 操美女的视频在线观看| 欧美97在线视频| 夫妻午夜视频| 91精品三级在线观看| 亚洲,一卡二卡三卡| 精品少妇黑人巨大在线播放| 嫩草影院入口| 色网站视频免费| 日本欧美国产在线视频| 麻豆av在线久日| 久久精品久久久久久久性| 婷婷色综合大香蕉| 亚洲精品一区蜜桃| av网站在线播放免费| 80岁老熟妇乱子伦牲交| 久久久久久久精品精品| 欧美日韩综合久久久久久| 婷婷色综合大香蕉| 亚洲av综合色区一区| 国产精品成人在线| 亚洲av成人不卡在线观看播放网 | 波多野结衣一区麻豆| 国产日韩欧美视频二区| 别揉我奶头~嗯~啊~动态视频 | 午夜激情av网站| 日韩大码丰满熟妇| 纯流量卡能插随身wifi吗| 日本vs欧美在线观看视频| 国产一区有黄有色的免费视频| 国产成人啪精品午夜网站| 在线观看免费视频网站a站| 久久午夜综合久久蜜桃| 伦理电影大哥的女人| av电影中文网址| 国产av码专区亚洲av| 国产精品久久久久久精品电影小说| 又大又爽又粗| 天天躁夜夜躁狠狠躁躁| 成人国语在线视频| 国产午夜精品一二区理论片| 久久精品久久精品一区二区三区| 丝袜美腿诱惑在线| 水蜜桃什么品种好| 哪个播放器可以免费观看大片| 久久久久久人人人人人| 国产成人欧美在线观看 | 中文天堂在线官网| 婷婷色av中文字幕| 一级毛片 在线播放| 亚洲国产看品久久| 91老司机精品| 99精国产麻豆久久婷婷| 亚洲欧美精品综合一区二区三区| 一区二区三区精品91| 男人舔女人的私密视频| 中文精品一卡2卡3卡4更新| kizo精华| 中文字幕av电影在线播放| 色网站视频免费| 丝袜脚勾引网站| 91国产中文字幕| 成年av动漫网址| 亚洲av福利一区| 亚洲,欧美精品.| 亚洲成人手机| 亚洲av欧美aⅴ国产| 亚洲精品国产区一区二| 日韩大片免费观看网站| 亚洲视频免费观看视频| 少妇猛男粗大的猛烈进出视频| 精品免费久久久久久久清纯 | 久久精品aⅴ一区二区三区四区| 青青草视频在线视频观看| 亚洲人成77777在线视频| a级毛片黄视频| 亚洲美女黄色视频免费看| 亚洲美女搞黄在线观看| 午夜激情久久久久久久| 国产探花极品一区二区| 2021少妇久久久久久久久久久| 汤姆久久久久久久影院中文字幕| 老熟女久久久| 国产精品三级大全| 大香蕉久久网| 国产精品久久久久久精品电影小说| 国产一区二区三区av在线| 日韩制服丝袜自拍偷拍| 一本—道久久a久久精品蜜桃钙片| 久久久欧美国产精品| 国产成人a∨麻豆精品| 欧美激情高清一区二区三区 | 人人妻人人澡人人看| 欧美97在线视频| 午夜免费鲁丝| 亚洲欧美一区二区三区黑人| 国产无遮挡羞羞视频在线观看| 亚洲av男天堂| a级片在线免费高清观看视频| 久久97久久精品| 中文字幕另类日韩欧美亚洲嫩草| 亚洲欧美日韩另类电影网站| 久久这里只有精品19| 精品国产一区二区三区四区第35| 亚洲精品久久午夜乱码| 国产成人精品久久久久久| 久久久久久免费高清国产稀缺| 一级毛片 在线播放| 1024香蕉在线观看| 一级毛片我不卡| 亚洲精品久久久久久婷婷小说| 亚洲国产精品国产精品| 国产人伦9x9x在线观看| 久久国产精品大桥未久av| 777久久人妻少妇嫩草av网站| 免费在线观看黄色视频的| 日韩电影二区| 久久 成人 亚洲| av网站免费在线观看视频| 纯流量卡能插随身wifi吗| 一边亲一边摸免费视频| 伦理电影大哥的女人| 国产精品99久久99久久久不卡 | 国产成人系列免费观看| 黄色怎么调成土黄色| 狠狠婷婷综合久久久久久88av| 久久精品久久久久久久性| 各种免费的搞黄视频| 在线精品无人区一区二区三| 老司机亚洲免费影院| 国产熟女午夜一区二区三区| 美女中出高潮动态图| 久久97久久精品| 国产精品久久久久久人妻精品电影 | 久久精品国产亚洲av高清一级| 婷婷成人精品国产| 国产亚洲欧美精品永久| 成年人免费黄色播放视频| 亚洲一区中文字幕在线| 日韩精品有码人妻一区| 国产午夜精品一二区理论片| 欧美乱码精品一区二区三区| 91国产中文字幕| 国产精品久久久av美女十八| 久久99热这里只频精品6学生| 国产成人精品无人区| 精品国产超薄肉色丝袜足j| 国产欧美日韩综合在线一区二区| 狠狠精品人妻久久久久久综合| e午夜精品久久久久久久| 国产日韩欧美在线精品| 亚洲欧洲精品一区二区精品久久久 | xxx大片免费视频| 午夜影院在线不卡| 久久午夜综合久久蜜桃| 色94色欧美一区二区| 老司机亚洲免费影院| 国产极品天堂在线| 嫩草影视91久久| 一边亲一边摸免费视频| 亚洲久久久国产精品| netflix在线观看网站| 久久久久久久久久久久大奶| 久久精品亚洲熟妇少妇任你| 婷婷色麻豆天堂久久| 免费黄网站久久成人精品| 国产男人的电影天堂91| 久久亚洲国产成人精品v| 亚洲精品av麻豆狂野| 天天添夜夜摸| 亚洲精品第二区| 男人操女人黄网站| 亚洲欧美色中文字幕在线| 国产精品欧美亚洲77777| 欧美黑人欧美精品刺激| 亚洲精品av麻豆狂野| 成年美女黄网站色视频大全免费| 国产 精品1| 亚洲一卡2卡3卡4卡5卡精品中文| 国产亚洲精品第一综合不卡| svipshipincom国产片| 一本一本久久a久久精品综合妖精| 日日摸夜夜添夜夜爱| 黄频高清免费视频| 欧美成人午夜精品| 一级片免费观看大全|