• <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)
    亚洲国产精品成人久久小说| 久久综合国产亚洲精品| 国产精品成人在线| 91老司机精品| 男女国产视频网站| 欧美另类亚洲清纯唯美| 亚洲avbb在线观看| 国产成人av激情在线播放| 黄色毛片三级朝国网站| 亚洲欧美清纯卡通| 欧美精品一区二区大全| 久久久国产精品麻豆| 国产在线观看jvid| 久久久欧美国产精品| 久久精品国产亚洲av香蕉五月 | 免费不卡黄色视频| 一级毛片精品| 飞空精品影院首页| 亚洲成av片中文字幕在线观看| 中国美女看黄片| 国产亚洲欧美精品永久| av不卡在线播放| 国产野战对白在线观看| 大陆偷拍与自拍| 国产成人av教育| 国产高清国产精品国产三级| av免费在线观看网站| 国产日韩欧美亚洲二区| 黑丝袜美女国产一区| 亚洲国产av影院在线观看| 亚洲精品av麻豆狂野| 国产精品久久久久久精品古装| 午夜福利一区二区在线看| 欧美在线一区亚洲| 高清av免费在线| 欧美性长视频在线观看| 亚洲欧美色中文字幕在线| 叶爱在线成人免费视频播放| 亚洲av男天堂| 国产亚洲精品一区二区www | 亚洲成人免费电影在线观看| 国产精品香港三级国产av潘金莲| 在线天堂中文资源库| 免费观看人在逋| 日韩熟女老妇一区二区性免费视频| 999久久久国产精品视频| 亚洲国产毛片av蜜桃av| 久久精品aⅴ一区二区三区四区| 99精品欧美一区二区三区四区| 成人亚洲精品一区在线观看| 在线观看免费视频网站a站| 人妻人人澡人人爽人人| 亚洲欧美清纯卡通| 性色av一级| 日本91视频免费播放| 国产97色在线日韩免费| 一级毛片电影观看| 我的亚洲天堂| 国产成人免费无遮挡视频| 国产1区2区3区精品| 黑人操中国人逼视频| 男人爽女人下面视频在线观看| 亚洲五月色婷婷综合| 亚洲少妇的诱惑av| 欧美老熟妇乱子伦牲交| 老汉色av国产亚洲站长工具| 另类亚洲欧美激情| 丝袜在线中文字幕| 免费黄频网站在线观看国产| 免费在线观看完整版高清| 精品免费久久久久久久清纯 | 久久久久久免费高清国产稀缺| 亚洲成人免费电影在线观看| 欧美少妇被猛烈插入视频| www.精华液| 国产一区二区 视频在线| 制服人妻中文乱码| 黄色a级毛片大全视频| 在线观看舔阴道视频| a级毛片黄视频| 亚洲专区中文字幕在线| 99国产精品免费福利视频| 男人操女人黄网站| 黄色毛片三级朝国网站| 97在线人人人人妻| 免费观看人在逋| 色综合欧美亚洲国产小说| 久久久久网色| 精品一区二区三卡| 人人妻人人澡人人爽人人夜夜| 男女之事视频高清在线观看| 欧美精品一区二区免费开放| 亚洲伊人久久精品综合| 美女扒开内裤让男人捅视频| 久久人人爽av亚洲精品天堂| 中文字幕高清在线视频| 美女国产高潮福利片在线看| 激情视频va一区二区三区| av在线app专区| 久热这里只有精品99| 老司机影院毛片| 亚洲国产av影院在线观看| 亚洲欧美色中文字幕在线| 亚洲欧美一区二区三区黑人| 男男h啪啪无遮挡| 亚洲熟女毛片儿| 亚洲少妇的诱惑av| 午夜免费观看性视频| 国产精品久久久人人做人人爽| 亚洲一码二码三码区别大吗| 91麻豆av在线| 人妻人人澡人人爽人人| 精品国产超薄肉色丝袜足j| 亚洲第一青青草原| 国产黄色免费在线视频| 久久久国产一区二区| 亚洲 欧美一区二区三区| 手机成人av网站| 亚洲国产欧美日韩在线播放| 亚洲精品一二三| 免费看十八禁软件| 黑人猛操日本美女一级片| 日韩 亚洲 欧美在线| 久久久久国产精品人妻一区二区| 国产日韩欧美在线精品| 人人妻人人澡人人爽人人夜夜| 久久这里只有精品19| 99久久99久久久精品蜜桃| 丝瓜视频免费看黄片| 成年人免费黄色播放视频| 国产成人影院久久av| 欧美大码av| 亚洲午夜精品一区,二区,三区| 建设人人有责人人尽责人人享有的| 日本91视频免费播放| 久久人妻熟女aⅴ| e午夜精品久久久久久久| 一本久久精品| 午夜精品国产一区二区电影| 视频区欧美日本亚洲| 国产精品成人在线| 中文字幕人妻熟女乱码| 国产男女内射视频| 亚洲精品久久午夜乱码| 国产精品熟女久久久久浪| 韩国精品一区二区三区| 久久狼人影院| 日韩欧美免费精品| 亚洲精品久久久久久婷婷小说| 成人国产av品久久久| 99国产综合亚洲精品| a在线观看视频网站| 亚洲人成电影观看| 91国产中文字幕| 操出白浆在线播放| 我的亚洲天堂| 亚洲精品成人av观看孕妇| 在线观看人妻少妇| 另类亚洲欧美激情| 亚洲一码二码三码区别大吗| 夜夜夜夜夜久久久久| 亚洲欧洲精品一区二区精品久久久| 91九色精品人成在线观看| 王馨瑶露胸无遮挡在线观看| 亚洲国产欧美日韩在线播放| 巨乳人妻的诱惑在线观看| www.999成人在线观看| 欧美在线一区亚洲| 日本wwww免费看| 法律面前人人平等表现在哪些方面 | 王馨瑶露胸无遮挡在线观看| 亚洲精品第二区| 国产日韩欧美在线精品| 99精国产麻豆久久婷婷| 成人黄色视频免费在线看| 成年女人毛片免费观看观看9 | 欧美黑人精品巨大| 亚洲中文av在线| 一区二区日韩欧美中文字幕| 国产精品九九99| 麻豆av在线久日| 亚洲国产精品一区二区三区在线| 亚洲一码二码三码区别大吗| 日韩精品免费视频一区二区三区| 人妻 亚洲 视频| 成年动漫av网址| 亚洲精品国产精品久久久不卡| 欧美xxⅹ黑人| 久久精品久久久久久噜噜老黄| 久久ye,这里只有精品| 亚洲国产中文字幕在线视频| av不卡在线播放| 亚洲av成人一区二区三| 精品少妇一区二区三区视频日本电影| 欧美 日韩 精品 国产| 亚洲,欧美精品.| 人人妻人人爽人人添夜夜欢视频| 亚洲七黄色美女视频| 成人三级做爰电影| 国产在视频线精品| 涩涩av久久男人的天堂| 一个人免费在线观看的高清视频 | 亚洲欧美一区二区三区久久| 亚洲国产欧美日韩在线播放| 久久久欧美国产精品| 男女床上黄色一级片免费看| 亚洲精华国产精华精| 欧美黄色淫秽网站| 亚洲精品国产色婷婷电影| 亚洲欧美一区二区三区久久| 亚洲成人手机| 一二三四在线观看免费中文在| tube8黄色片| 亚洲av日韩在线播放| 女人精品久久久久毛片| 高清欧美精品videossex| 免费日韩欧美在线观看| 侵犯人妻中文字幕一二三四区| 在线观看人妻少妇| 美女脱内裤让男人舔精品视频| 美女脱内裤让男人舔精品视频| 成人国语在线视频| 热99国产精品久久久久久7| 少妇被粗大的猛进出69影院| 老鸭窝网址在线观看| 男女高潮啪啪啪动态图| 欧美乱码精品一区二区三区| 新久久久久国产一级毛片| 免费在线观看完整版高清| 视频在线观看一区二区三区| 欧美少妇被猛烈插入视频| 午夜激情av网站| 久久精品国产亚洲av香蕉五月 | 女人高潮潮喷娇喘18禁视频| 宅男免费午夜| 亚洲精华国产精华精| 夜夜夜夜夜久久久久| 美国免费a级毛片| 亚洲伊人色综图| 国产精品自产拍在线观看55亚洲 | 亚洲 国产 在线| 极品少妇高潮喷水抽搐| 考比视频在线观看| 男女免费视频国产| 午夜成年电影在线免费观看| 久久久久久免费高清国产稀缺| tube8黄色片| 日本a在线网址| 人妻久久中文字幕网| 国产精品熟女久久久久浪| 精品人妻在线不人妻| 这个男人来自地球电影免费观看| 人人妻人人添人人爽欧美一区卜| 蜜桃国产av成人99| 久久亚洲国产成人精品v| 波多野结衣一区麻豆| 50天的宝宝边吃奶边哭怎么回事| 午夜视频精品福利| 美女高潮到喷水免费观看| 三级毛片av免费| 热99国产精品久久久久久7| 1024香蕉在线观看| 一进一出抽搐动态| 热99国产精品久久久久久7| 国产欧美日韩一区二区三 | 中国国产av一级| 建设人人有责人人尽责人人享有的| 国产欧美日韩一区二区三 | 免费在线观看视频国产中文字幕亚洲 | 肉色欧美久久久久久久蜜桃| 成人三级做爰电影| 日韩一区二区三区影片| 亚洲精品美女久久av网站| 亚洲性夜色夜夜综合| 黑丝袜美女国产一区| 大码成人一级视频| 国产精品一二三区在线看| 国产精品九九99| 一本色道久久久久久精品综合| 亚洲欧美清纯卡通| 国产免费av片在线观看野外av| 日本一区二区免费在线视频| 色老头精品视频在线观看| 亚洲精品粉嫩美女一区| av网站在线播放免费| 免费在线观看影片大全网站| 欧美日韩亚洲综合一区二区三区_| 免费久久久久久久精品成人欧美视频| 欧美成人午夜精品| 一个人免费在线观看的高清视频 | 人妻人人澡人人爽人人| 日本精品一区二区三区蜜桃| 亚洲免费av在线视频| 激情视频va一区二区三区| 黄频高清免费视频| 黄片播放在线免费| 丁香六月天网| 亚洲色图 男人天堂 中文字幕| 国产亚洲精品一区二区www | 免费久久久久久久精品成人欧美视频| 婷婷成人精品国产| 叶爱在线成人免费视频播放| 久久人人爽av亚洲精品天堂| 亚洲avbb在线观看| 国产亚洲一区二区精品| 中文字幕另类日韩欧美亚洲嫩草| 久久久久精品国产欧美久久久 | 99精国产麻豆久久婷婷| 制服人妻中文乱码| 人人妻,人人澡人人爽秒播| 日韩人妻精品一区2区三区| 丁香六月欧美| 亚洲精品一区蜜桃| 中文字幕人妻丝袜制服| 精品一区二区三区av网在线观看 | 亚洲国产精品成人久久小说| 亚洲精品粉嫩美女一区| √禁漫天堂资源中文www| 国产欧美日韩一区二区三区在线| 一级,二级,三级黄色视频| 在线观看免费高清a一片| 久久精品国产亚洲av香蕉五月 | 亚洲精品美女久久av网站| 婷婷丁香在线五月| 黄色 视频免费看| 69av精品久久久久久 | 又紧又爽又黄一区二区| 99精国产麻豆久久婷婷| 高潮久久久久久久久久久不卡| 久久人人97超碰香蕉20202| 交换朋友夫妻互换小说| 亚洲欧洲精品一区二区精品久久久| 亚洲全国av大片| 欧美亚洲 丝袜 人妻 在线| 中文字幕av电影在线播放| 午夜免费观看性视频| 国产国语露脸激情在线看| 欧美日韩av久久| 91精品三级在线观看| 国产精品秋霞免费鲁丝片| 亚洲精品一卡2卡三卡4卡5卡 | 热99re8久久精品国产| 女性被躁到高潮视频| kizo精华| 欧美黄色淫秽网站| 一区二区三区激情视频| 亚洲人成77777在线视频| 国产精品欧美亚洲77777| 老司机靠b影院| 天天操日日干夜夜撸| 亚洲色图综合在线观看| 精品一区二区三卡| 欧美日韩视频精品一区| 中文字幕精品免费在线观看视频| 人人妻人人爽人人添夜夜欢视频| 国产97色在线日韩免费| 国产亚洲一区二区精品| 国产福利在线免费观看视频| 精品少妇久久久久久888优播| 亚洲 国产 在线| 伊人亚洲综合成人网| 日本av手机在线免费观看| 亚洲国产毛片av蜜桃av| 后天国语完整版免费观看| 欧美日韩精品网址| 日韩视频一区二区在线观看| 国产高清国产精品国产三级| 一个人免费看片子| 国产av又大| 国产成人免费无遮挡视频| 欧美激情 高清一区二区三区| 成年女人毛片免费观看观看9 | 欧美黑人欧美精品刺激| 夜夜夜夜夜久久久久| 国产成人欧美| 亚洲国产欧美一区二区综合| 王馨瑶露胸无遮挡在线观看| 午夜福利,免费看| 欧美中文综合在线视频| 成年av动漫网址| 十八禁人妻一区二区| 亚洲精品国产av成人精品| 久久久欧美国产精品| 最新在线观看一区二区三区| 三上悠亚av全集在线观看| 母亲3免费完整高清在线观看| 99九九在线精品视频| 亚洲男人天堂网一区| 亚洲精品乱久久久久久| 欧美成人午夜精品| 免费观看av网站的网址| 人人妻人人添人人爽欧美一区卜| 国产一区有黄有色的免费视频| 91国产中文字幕| 免费在线观看视频国产中文字幕亚洲 | 精品一区二区三卡| 国产精品成人在线| 亚洲一区中文字幕在线| 91精品三级在线观看| 大片电影免费在线观看免费| 欧美日韩福利视频一区二区| 久久人人爽人人片av| av网站免费在线观看视频| 色视频在线一区二区三区| 香蕉国产在线看| 国产成人欧美| 亚洲国产av新网站| 午夜福利,免费看| 高潮久久久久久久久久久不卡| 日本黄色日本黄色录像| 精品国产国语对白av| 精品国产乱码久久久久久男人| 国产成+人综合+亚洲专区| 久久精品国产亚洲av高清一级| 亚洲一区中文字幕在线| 亚洲人成电影观看| 在线亚洲精品国产二区图片欧美| e午夜精品久久久久久久| 国产在线观看jvid| av一本久久久久| 国产欧美日韩一区二区三 | 午夜影院在线不卡| 夫妻午夜视频| www.999成人在线观看| 黄色毛片三级朝国网站| 国内毛片毛片毛片毛片毛片| 国产精品99久久99久久久不卡| 国产成人影院久久av| 丝袜在线中文字幕| 99久久精品国产亚洲精品| 亚洲第一av免费看| 欧美大码av| 老汉色av国产亚洲站长工具| 99九九在线精品视频| 国产免费福利视频在线观看| 黑人巨大精品欧美一区二区mp4| 91国产中文字幕| 欧美日韩亚洲高清精品| 夜夜夜夜夜久久久久| 美女中出高潮动态图| 夫妻午夜视频| 久久久久精品人妻al黑| 一本久久精品| 中国美女看黄片| 蜜桃国产av成人99| 亚洲国产av新网站| 男女无遮挡免费网站观看| 成人18禁高潮啪啪吃奶动态图| 精品福利观看| 久久人妻福利社区极品人妻图片| 一级毛片女人18水好多| h视频一区二区三区| 国产精品秋霞免费鲁丝片| 亚洲国产欧美日韩在线播放| 嫁个100分男人电影在线观看| 啦啦啦啦在线视频资源| 久久这里只有精品19| 搡老岳熟女国产| 欧美另类一区| 在线天堂中文资源库| 久久 成人 亚洲| 精品人妻熟女毛片av久久网站| 18禁裸乳无遮挡动漫免费视频| 黄色a级毛片大全视频| 久热这里只有精品99| 高清av免费在线| 99国产精品99久久久久| 国产在线视频一区二区| 一区二区三区四区激情视频| 美女高潮喷水抽搐中文字幕| 一本—道久久a久久精品蜜桃钙片| 久久青草综合色| h视频一区二区三区| 精品一区二区三区av网在线观看 | 国产一区有黄有色的免费视频| svipshipincom国产片| 操美女的视频在线观看| 成年动漫av网址| av天堂久久9| 成人免费观看视频高清| 美女国产高潮福利片在线看| 国产精品偷伦视频观看了| 少妇粗大呻吟视频| 大码成人一级视频| 欧美黑人欧美精品刺激| 一本久久精品| 这个男人来自地球电影免费观看| www.自偷自拍.com| 国产亚洲欧美精品永久| 十八禁网站免费在线| 女人爽到高潮嗷嗷叫在线视频| 最近最新免费中文字幕在线| 成人av一区二区三区在线看 | 久久亚洲精品不卡| 啦啦啦中文免费视频观看日本| 精品久久久久久久毛片微露脸 | 免费观看人在逋| 亚洲欧美日韩另类电影网站| 亚洲欧美成人综合另类久久久| 免费少妇av软件| 国产一区二区 视频在线| 国产免费现黄频在线看| 曰老女人黄片| 一区在线观看完整版| 在线观看人妻少妇| 久久热在线av| 捣出白浆h1v1| 中文字幕高清在线视频| 天天操日日干夜夜撸| 亚洲欧美激情在线| 操美女的视频在线观看| 日韩精品免费视频一区二区三区| 人人妻人人添人人爽欧美一区卜| 肉色欧美久久久久久久蜜桃| 一二三四在线观看免费中文在| 欧美人与性动交α欧美精品济南到| 久久精品国产a三级三级三级| 欧美国产精品一级二级三级| 黄色视频不卡| 一本大道久久a久久精品| 超色免费av| 自线自在国产av| 爱豆传媒免费全集在线观看| 黄片小视频在线播放| 女人爽到高潮嗷嗷叫在线视频| 亚洲激情五月婷婷啪啪| 99国产极品粉嫩在线观看| 美女扒开内裤让男人捅视频| 亚洲中文av在线| 欧美激情久久久久久爽电影 | 99久久综合免费| 亚洲精华国产精华精| 999久久久国产精品视频| 美女扒开内裤让男人捅视频| 五月开心婷婷网| 9热在线视频观看99| av福利片在线| 国产淫语在线视频| 欧美久久黑人一区二区| 国产精品av久久久久免费| 精品久久久精品久久久| 天天添夜夜摸| 中文字幕精品免费在线观看视频| 老汉色av国产亚洲站长工具| 麻豆av在线久日| 老司机在亚洲福利影院| av有码第一页| 国产在视频线精品| 老鸭窝网址在线观看| 老汉色∧v一级毛片| 菩萨蛮人人尽说江南好唐韦庄| 国产在线观看jvid| 国产成人欧美在线观看 | av在线播放精品| 日韩视频一区二区在线观看| 欧美日韩av久久| 免费观看人在逋| 在线 av 中文字幕| 99精国产麻豆久久婷婷| a级毛片黄视频| 免费人妻精品一区二区三区视频| 在线亚洲精品国产二区图片欧美| 成人国产av品久久久| 丝袜在线中文字幕| 久久香蕉激情| 欧美+亚洲+日韩+国产| 精品欧美一区二区三区在线| 精品一品国产午夜福利视频| 国产精品一区二区免费欧美 | 亚洲av片天天在线观看| 中文字幕人妻丝袜制服| 国产精品一区二区精品视频观看| 男女边摸边吃奶| 精品国产超薄肉色丝袜足j| 嫁个100分男人电影在线观看| 亚洲精品国产区一区二| 久久久精品国产亚洲av高清涩受| 天天躁狠狠躁夜夜躁狠狠躁| 午夜成年电影在线免费观看| 蜜桃在线观看..| 亚洲 欧美一区二区三区| 午夜福利免费观看在线| 亚洲专区字幕在线| a级毛片在线看网站| 男人舔女人的私密视频| 91大片在线观看| 男女之事视频高清在线观看| 午夜精品国产一区二区电影| 国产人伦9x9x在线观看| 欧美xxⅹ黑人| 亚洲专区中文字幕在线| av天堂久久9| 亚洲精品中文字幕一二三四区 | 亚洲av电影在线进入| 丰满饥渴人妻一区二区三| 成年女人毛片免费观看观看9 | 十八禁高潮呻吟视频| 成人18禁高潮啪啪吃奶动态图| 久久精品国产亚洲av高清一级| 亚洲 国产 在线| 黑人欧美特级aaaaaa片| 大型av网站在线播放| 日韩欧美国产一区二区入口| 精品第一国产精品| 交换朋友夫妻互换小说| 国产日韩欧美亚洲二区| 久久久久久久精品精品| 亚洲精华国产精华精| 国产精品麻豆人妻色哟哟久久| 久久 成人 亚洲| 国产黄频视频在线观看| 精品视频人人做人人爽| 色播在线永久视频| 伊人亚洲综合成人网| 婷婷成人精品国产|