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

    雙季早晚稻不同發(fā)育階段日數(shù)對溫度變化的敏感性比較*

    2018-03-19 07:29:23王治海金志鳳鄔定榮毛智軍陳中赟
    中國農(nóng)業(yè)氣象 2018年3期
    關(guān)鍵詞:發(fā)育階段孕穗日數(shù)

    王治海,金志鳳,鄔定榮,毛智軍,陳中赟

    ?

    雙季早晚稻不同發(fā)育階段日數(shù)對溫度變化的敏感性比較*

    王治海1,金志鳳1,鄔定榮2**,毛智軍3,陳中赟4

    (1.浙江省氣候中心,杭州 310017;2.中國氣象科學(xué)研究院,北京 100081;3.龍游縣氣象局,衢州 324400;4.湖州市氣象局,湖州 313000)

    基于1981-2010年長江中下游地區(qū)38個(gè)農(nóng)業(yè)氣象觀測站雙季早、晚稻發(fā)育期及同期逐日氣象數(shù)據(jù),結(jié)合水稻特性將發(fā)育期劃分為出苗-移栽、移栽-返青、返青-孕穗、孕穗-成熟4個(gè)主要發(fā)育階段,比較雙季早、晚稻不同發(fā)育階段日數(shù)對溫度變化的相對敏感性,以探究不同作物、不同發(fā)育期對氣候變化響應(yīng)的差異。結(jié)果表明:水稻各發(fā)育階段日數(shù)隨溫度升高基本呈減少趨勢,但相對敏感性差異較大。早稻返青-孕穗期日數(shù)對溫度變化最敏感,相對敏感性平均為?0.094個(gè)點(diǎn)·℃?1(P<0.05),晚稻孕穗-成熟期日數(shù)對溫度變化最不敏感,相對敏感性平均為?0.045個(gè)點(diǎn)·℃?1(P<0.05)。相同發(fā)育階段內(nèi),返青前晚稻對溫度變化更敏感,返青后早稻更敏感。除移栽-返青期外,兩種作物其它營養(yǎng)生長階段敏感性均強(qiáng)于生殖生長。相對敏感性與階段日數(shù)變異系數(shù)呈負(fù)相關(guān)關(guān)系,日數(shù)波動越大,階段日數(shù)隨溫度升高而減少越明顯,尤以早稻返青-孕穗期相關(guān)系數(shù)達(dá)?0.761(P<0.001)。

    氣候變暖;發(fā)育階段日數(shù);相對敏感性;水稻;長江中下游地區(qū)

    IPCC第五次評估報(bào)告指出,1880-2012年全球平均地表溫度上升了0.85℃,1951-2012年增溫速率接近1880年以來的兩倍[1]。氣候變暖使作物生長季內(nèi)的光溫水等氣候資源發(fā)生變化,進(jìn)而影響作物的生長發(fā)育和產(chǎn)量形成[2-3]。作物生長模型是研究氣候變化對作物產(chǎn)量影響的主要工具,但不同模型在模擬產(chǎn)量時(shí)存在較大的不確定性,其生育期模塊采用的溫度響應(yīng)曲線的差異是不確定性的主要來源[4]。因此,科學(xué)分析作物發(fā)育期對氣候變化的響應(yīng)機(jī)制是提高作物生長模型預(yù)測產(chǎn)量精度的關(guān)鍵,對制定農(nóng)業(yè)發(fā)展策略、保障糧食安全具有重要意義。

    目前,圍繞氣候變暖對大宗作物物候的影響研究取得了不少成果。針對水稻作物,李正國等[5-6]分別在東北、華北稻區(qū)的研究表明,溫度升高引起水稻播種和移栽日期提前、生育期延長的趨勢;但也有學(xué)者在緯度較低的華中[7]、華南[8]以及江淮地區(qū)[9]研究發(fā)現(xiàn),水稻生育期因增溫而縮短,但縮短速率存在區(qū)域差異。圍繞麥類作物,肖登攀等[10]發(fā)現(xiàn)華北冬小麥營養(yǎng)生長階段以7.4d·10a?1的速率縮短,而生殖生長階段以1.7d·10a?1的速率延長,即氣候變化對不同發(fā)育階段日數(shù)的影響不同;蔣菊芳等[11]對河西走廊春小麥的物候研究也有類似結(jié)論,三葉-抽穗、拔節(jié)-開花期分別以0.3d·10a?1、0.26d·10a?1的速率縮短,而孕穗-乳熟期以0.23d·10a?1的速率延長。對于玉米作物,孟林等[12]指出,平均氣溫每升高1℃,華北夏玉米全生育期和生殖生長期分別縮短2.7d、1.1d。以上研究均通過發(fā)育階段絕對變化日數(shù),反映不同糧食作物物候期對氣候變化的響應(yīng)情況。但受地域、品種、種植制度差異的影響,不同地區(qū)、不同作物的不同發(fā)育期對氣候變暖的響應(yīng)不一致。溫度每升高1℃,經(jīng)歷天數(shù)較多的發(fā)育階段其絕對變化日數(shù)可能比天數(shù)短的發(fā)育期更多。此外,各地區(qū)為適應(yīng)氣候變暖對作物播期的調(diào)整也可引起生育期天數(shù)的變化。因此,基于絕對速率的敏感性不適于直接比較不同地區(qū)、不同作物、不同生育期的溫度敏感性。對此,高靜等[13]提出相對敏感性的概念,分析了華北冬小麥不同發(fā)育階段日數(shù)對溫度變化的敏感性。理論上,該指標(biāo)可用于比較不同作物、不同發(fā)育階段之間的溫度敏感性。

    水稻是中國最重要的糧食作物之一,長江中下游地區(qū)是中國最大的水稻產(chǎn)區(qū)[14]。水稻因種植制度不同,可分為雙季早稻和雙季晚稻,兩者生長發(fā)育對溫度變化的敏感性可能存在一定差異。此外,水稻移栽容易導(dǎo)致秧苗根系損傷,推遲整個(gè)發(fā)育進(jìn)程,這種影響也可能會在溫度敏感性上體現(xiàn)。因此,本研究擬以長江中下游移栽雙季稻為研究對象,利用相對敏感性分析方法,比較雙季早、晚稻不同發(fā)育階段日數(shù)對溫度變化的相對敏感性及其區(qū)域性差異,以期為氣候變化背景下作物生育期管理、合理布局和制定適應(yīng)性政策等提供科學(xué)依據(jù)。

    1 資料與方法

    1.1 研究區(qū)概況

    長江中下游地區(qū)包括江蘇、安徽、浙江、江西、湖北、湖南、上海7?。ㄊ校貏菸鞲邧|低,以平原和低山丘陵為主。地處亞熱帶溫暖濕潤季風(fēng)氣候區(qū),雨熱同期且降水充沛,年平均氣溫14~18℃,年降水量1000~1500mm,無霜期210~270d,10℃以上活動積溫達(dá)4500~5000℃·d[15]。根據(jù)水稻種植區(qū)劃[16],整個(gè)區(qū)域在1981-2010年內(nèi)以雙季稻為主[17]。

    1.2 資料

    選取1981-2010年長江中下游地區(qū)連續(xù)10a以上有雙季早、晚稻移栽種植的38個(gè)農(nóng)業(yè)氣象觀測站(圖1),各站早稻和晚稻均有25a觀測記錄。發(fā)育期資料包括雙季稻出苗、移栽、返青、孕穗、成熟日期;氣象資料為與農(nóng)業(yè)氣象觀測站對應(yīng)的38個(gè)氣象站1981-2010年的逐日平均氣溫。發(fā)育期和氣象數(shù)據(jù)均來源于國家氣象信息中心。

    1.3 方法

    多數(shù)作物的全生育期一般可分為營養(yǎng)生長(出苗-孕穗)和生殖生長(孕穗-成熟)兩個(gè)階段,但與直播作物不同的是,移栽水稻的營養(yǎng)生長階段增加了移栽和返青兩個(gè)生育期。移栽可以延長水稻生育期,并通過控制田間密度以保證穗粒數(shù),利于提高產(chǎn)量。但移栽也可引起水稻根系損傷,損傷程度因秧苗特性和天氣條件的不同而異,輕則水分失衡重則死亡[18]。移栽后下一個(gè)生育階段為返青,移栽-返青階段水稻生長對低溫環(huán)境敏感,溫度過低容易延遲秧苗返青甚至不能返青[19]。由此可見,移栽和返青是移栽水稻對環(huán)境溫度要求較為特別的生育期,需要從營養(yǎng)生長階段單獨(dú)劃分出來進(jìn)行分析。本研究將移栽水稻整個(gè)生育期劃分為4個(gè)主要階段,依次為出苗-移栽(E-T)、移栽-返青(T-R)、返青-孕穗(R-B)、孕穗-成熟(B-M)。前3個(gè)階段是營養(yǎng)生長階段,后1個(gè)是生殖生長階段。

    圖1 長江中下游地區(qū)38個(gè)雙季稻農(nóng)業(yè)氣象觀測站點(diǎn)分布

    以往研究中,作物發(fā)育期的溫度敏感性多以階段日數(shù)受溫度影響的絕對變化量來表達(dá)[20]。但絕對變化速率受該發(fā)育階段長度影響,無法用于不同作物、不同發(fā)育期之間的溫度敏感性比較。為此,高靜等[13]提出了相對敏感性的概念,即發(fā)育階段相對日數(shù)與同期日平均氣溫的線性回歸關(guān)系,計(jì)算式為

    式中,D為某站水稻完成某發(fā)育階段日數(shù)(d),Daver為該站水稻完成該發(fā)育階段日數(shù)的多年平均值(d),T為同期日平均氣溫(℃),a、b為待定參數(shù),其中b值為該發(fā)育階段日數(shù)對溫度變化的相對敏感性(個(gè)點(diǎn)·℃-1)。由于采用的是發(fā)育日數(shù)的相對值,理論上該指標(biāo)可用于比較不同作物以及同一作物不同發(fā)育階段之間的溫度敏感性差異。

    利用C#計(jì)算機(jī)編程語言,統(tǒng)計(jì)各站點(diǎn)4個(gè)主要發(fā)育階段日數(shù)及同期日平均氣溫,基于最小二乘法擬合式(1)中的參數(shù)b。借用SPSS軟件進(jìn)行顯著性檢驗(yàn)分析,基于ArcGIS軟件繪制各站點(diǎn)的平均氣溫、階段日數(shù)以及相對敏感性的空間分布圖。

    2 結(jié)果與分析

    2.1 雙季早稻、晚稻主要發(fā)育階段平均氣溫的變化

    1981-2010年統(tǒng)計(jì)結(jié)果顯示,長江中下游地區(qū)雙季早稻整個(gè)生育期(4-7月)歷年平均氣溫在22.2~24.7℃。由圖2a1-a4可見,出苗-移栽期(4-5月)歷年平均氣溫在16.1~20.0℃,平均18.2℃,相對較低;移栽-返青期(5月)歷年平均氣溫在18.5~22.7℃,平均21.2℃;返青-孕穗期(5-6月)歷年平均氣溫在22.1~24.9℃,平均23.5℃;孕穗-成熟期(6-7月)歷年平均氣溫在26.3~28.7℃,平均27.6℃。分析各站點(diǎn)關(guān)鍵生育期平均氣溫線性變化趨勢可知,34.2%的站點(diǎn)雙季早稻孕穗-成熟期平均氣溫呈顯著升高趨勢,這些站點(diǎn)平均增溫速率為1.11℃·10a?1(P<0.05),其它生育期內(nèi)也有不同比例站點(diǎn)的平均氣溫呈現(xiàn)顯著線性升高趨勢,其中出苗-移栽期和返青-孕穗期顯著增溫的站點(diǎn)比例分別為21.1%(P<0.05)、15.8%(P<0.05),而移栽-返青期僅有5.3%的站點(diǎn)增溫趨勢顯著(P<0.05)。

    雙季晚稻整個(gè)生育期(6-10月)歷年平均氣溫在23.4~27.1℃。由圖2b1-b4可見,出苗-移栽期(6-7月)歷年平均氣溫在27.0~29.7℃,平均28.4℃;移栽-返青期(7月)歷年平均氣溫在28.0~31.0℃,平均29.2℃;返青-孕穗期(7-9月)歷年平均氣溫在26.3~29.4℃,平均28.0℃;孕穗-成熟期(9-10月)歷年平均氣溫在18.7~23.9℃,平均22.1℃,相對較低。分析各站點(diǎn)雙季晚稻關(guān)鍵生育期平均氣溫線性變化趨勢可知,63.2%的站點(diǎn)孕穗-成熟期增溫趨勢顯著,氣溫平均升高速率為1.05℃·10a?1(P<0.05),該階段變暖特征最明顯;出苗-移栽期次之,36.8%的站點(diǎn)平均氣溫呈顯著升高趨勢,平均增溫速率為0.89℃·10a?1(P<0.05);返青-孕穗期和移栽-返青期增溫趨勢通過0.05水平顯著性檢驗(yàn)的站點(diǎn)僅占13.2%和2.6%。

    2.2 雙季早稻、晚稻主要發(fā)育階段持續(xù)日數(shù)的變化

    1981-2010年統(tǒng)計(jì)結(jié)果表明,長江中下游地區(qū)雙季早稻全生育期持續(xù)日數(shù)在97.6~117.7d,平均104.6d。由圖3a1-a4可見,返青-孕穗期持續(xù)日數(shù)最長,平均40.2d;孕穗-成熟期次之,平均34.0d;出苗-移栽期持續(xù)日數(shù)平均為25.8d;移栽-返青期最短,平均4.6d。分析區(qū)域內(nèi)各關(guān)鍵發(fā)育階段持續(xù)日數(shù)隨時(shí)間的變化趨勢可知,28.9%的站點(diǎn)出苗-移栽期持續(xù)日數(shù)呈顯著減少趨勢,平均線性傾向率為-4.57d·10a?1(P<0.05);而返青-孕穗期有26.3%站點(diǎn)的持續(xù)日數(shù)呈顯著增加趨勢,平均線性傾向率為5.84d·10a?1(P<0.05);孕穗-成熟期和移栽-返青期持續(xù)日數(shù)顯著變化的站點(diǎn)比例為34.2%、29.0%,其線性傾向率變化范圍分別為-10.00~4.32d·10a?1(P<0.05)、-2.64~1.61d·10a?1(P<0.05),其中為正值的站點(diǎn)分別占61.5%、36.4%,即這兩個(gè)階段的持續(xù)日數(shù)有增加也有減少趨勢,但無明顯的區(qū)域分布特征。

    圖2 1981-2010年雙季早稻、晚稻4個(gè)關(guān)鍵發(fā)育階段各站點(diǎn)平均氣溫及其線性變化趨勢

    Note: E-T is the stage from emergence to transplanting, T-R is the stage from transplanting to reviving, R-B is the stage from reviving to booting and B-M is the stage from booting to mature. The same as below

    雙季晚稻全生育期日數(shù)為106.6~135.5d,平均119.8d。由圖3b1-b4可見,孕穗-成熟期持續(xù)日數(shù)最長,平均48.5d;返青-孕穗期次之,平均38.7d;出苗-移栽期持續(xù)日數(shù)平均為28.7d;移栽-返青期最短,平均3.9d。分析各站點(diǎn)關(guān)鍵發(fā)育階段持續(xù)日數(shù)的線性變化趨勢可知,63.2%的站點(diǎn)出苗-移栽期持續(xù)日數(shù)呈顯著減少趨勢,平均線性傾向率為?5.83d·10a?1(P<0.05);孕穗-成熟期內(nèi)有36.8%的站點(diǎn)持續(xù)日數(shù)減少趨勢顯著,平均線性傾向率為?4.10d·10a?1(P<0.05);而在返青-孕穗期,34.2%的站點(diǎn)持續(xù)日數(shù)呈顯著增加趨勢,平均線性傾向率為3.84d·10a?1(P<0.05);移栽-返青期持續(xù)日數(shù)的線性變化趨勢通過顯著性檢驗(yàn)站點(diǎn)比例達(dá)42.1%,線性傾向率在?1.71~1.10d·10a?1(P<0.05),為負(fù)值的站點(diǎn)約占50%,即該階段日數(shù)有縮短也有延長趨勢,但無明顯的區(qū)域分布特征,這與各地移栽時(shí)間差異較大有關(guān)。

    2.3 早晚稻發(fā)育階段持續(xù)日數(shù)對溫度變化的相對敏感性分析

    2.3.1 相對敏感性比較

    由圖4可見,多數(shù)通過顯著性檢驗(yàn)站點(diǎn)的各生育階段日數(shù)對溫度變化的相對敏感性均為負(fù)值,主要是溫度升高加快了作物發(fā)育速率,導(dǎo)致發(fā)育階段日數(shù)隨之減少。但同一作物不同發(fā)育階段相對敏感性大小各不相同,同一發(fā)育期內(nèi)不同作物的相對敏感性也存在一定差異。早稻返青-孕穗期相對敏感性最強(qiáng),區(qū)域波動較大,39.5%的站點(diǎn)通過顯著性檢驗(yàn),這些站點(diǎn)的相對敏感性在?0.204~?0.042個(gè)點(diǎn)·℃?1(P<0.05),平均?0.094個(gè)點(diǎn)·℃?1(P<0.05),即溫度每升高1℃,階段日數(shù)平均減少9.4%;晚稻出苗-移栽期有42.1%的站點(diǎn)通過顯著性檢驗(yàn),相對敏感性在?0.162~?0.040個(gè)點(diǎn)·℃?1(P<0.05),平均?0.093個(gè)點(diǎn)·℃?1(P<0.05);早稻出苗-移栽期內(nèi),52.6%的站點(diǎn)通過顯著性檢驗(yàn),相對敏感性在?0.123~0.048個(gè)點(diǎn)·℃?1(P<0.05),平均?0.061個(gè)點(diǎn)·℃?1(P<0.05),其中僅浙江嘉興相對敏感性為正值(P<0.05);早稻移栽-返青與孕穗-成熟期相對敏感性平均值相近,在?0.053個(gè)點(diǎn)·℃?1左右(P<0.05),但通過顯著性檢驗(yàn)站點(diǎn)比例不同,分別為31.6%、60.5%,相對而言,早稻移栽-返青期相對敏感性在區(qū)域內(nèi)較不穩(wěn)定;晚稻孕穗-成熟期68.4%的站點(diǎn)通過顯著性檢驗(yàn),相對敏感性在?0.110~?0.021個(gè)點(diǎn)·℃?1(P<0.05),平均?0.045個(gè)點(diǎn)·℃?1(P<0.05),相對敏感性較弱,區(qū)域內(nèi)波動較小;晚稻返青-孕穗期相對敏感性通過顯著性檢驗(yàn)站點(diǎn)的比例為18.4%,相對敏感性在?0.072~0.141個(gè)點(diǎn)·℃?1(P<0.05),平均0.021個(gè)點(diǎn)·℃?1(P<0.05),其中僅江西的宜豐和樟樹、浙江的椒江和金華相對敏感性為正值(P<0.05);而晚稻移栽-返青期內(nèi)僅浙江紹興站相對敏感性通過了0.05水平的顯著性檢驗(yàn)。

    圖3 1981-2010年雙季早稻、晚稻4個(gè)關(guān)鍵發(fā)育階段各站點(diǎn)持續(xù)日數(shù)及其線性變化趨勢

    從全區(qū)平均情況看,早稻返青-孕穗期階段日數(shù)對溫度變化最敏感,晚稻孕穗-成熟期最不敏感。晚稻移栽-返青期和返青-孕穗期的相對敏感性通過顯著性檢驗(yàn)的站點(diǎn)比例較小,敏感性在區(qū)域內(nèi)均不穩(wěn)定。針對同一作物,移栽-返青期較其它階段相對敏感性最不穩(wěn)定,期間水稻生長發(fā)育對低溫敏感[21],加上該階段日數(shù)較短且受農(nóng)戶決策影響較大,因此,溫度升高對該階段日數(shù)的影響相對較小。

    2.3.2 相對敏感性差異分析

    為探究不同作物、不同發(fā)育階段對溫度變化相對敏感性差異的原因,統(tǒng)計(jì)雙季早稻、晚稻相對敏感性通過顯著性檢驗(yàn)站點(diǎn)的不同發(fā)育階段日數(shù)的變異系數(shù),分析相對敏感性與變異系數(shù)的相關(guān)關(guān)系。

    圖4 雙季早晚、晚稻4個(gè)關(guān)鍵發(fā)育階段日數(shù)對溫度變化相對敏感性的空間分布

    注: ?表示相對敏感性通過0.05水平的顯著性檢驗(yàn)

    Note: ? is P<0.05

    如圖5所示,除早稻移栽-返青期外,各階段相對敏感性均與持續(xù)日數(shù)變異系數(shù)呈顯著負(fù)相關(guān)關(guān)系,即持續(xù)日數(shù)波動越大,相對敏感性為負(fù)的絕對值越大,階段日數(shù)隨溫度升高而減少越明顯。這種關(guān)系以早稻返青-孕穗期最突出,相關(guān)系數(shù)絕對值最大,回歸關(guān)系通過了0.001水平的顯著性檢驗(yàn),變異系數(shù)較大,平均為0.14(P<0.05);早稻和晚稻的出苗-移栽期次之,相關(guān)系數(shù)在?0.60左右(P<0.05);早稻和晚稻的孕穗-成熟期相關(guān)系數(shù)約?0.50(P<0.05),日數(shù)變異系數(shù)在0.10左右(P<0.05),波動較小。晚稻移栽-返青及返青-孕穗期的相對敏感性通過顯著性檢驗(yàn)的站點(diǎn)較少,區(qū)域穩(wěn)定性很弱,相關(guān)系數(shù)也均未通過顯著性檢驗(yàn)。

    由上可見,各發(fā)育階段對溫度變化的相對敏感性與階段日數(shù)變異系數(shù)的相關(guān)性從強(qiáng)到弱排序?yàn)榉登?孕穗期>出苗-移栽期>孕穗-成熟期,這與其相對敏感性的對比結(jié)果基本一致。因此,各發(fā)育階段日數(shù)的變化差異可以在一定程度反映發(fā)育階段間的相對敏感性差異。這主要因?yàn)闇囟葪l件的區(qū)域差異可引起不同站點(diǎn)內(nèi)同一作物同一發(fā)育階段日數(shù)的不同,而發(fā)育階段日數(shù)的差異反過來影響該階段的平均溫度。早稻和晚稻的移栽-返青期相對敏感性與階段日數(shù)變異系數(shù)相關(guān)性均未通過顯著性檢驗(yàn),但該階段日數(shù)變異系數(shù)大于其它階段,表明該階段相對敏感性可能還受別的因素影響。移栽是水稻種植的關(guān)鍵環(huán)節(jié),適宜的移栽期有利于提高有效穗數(shù)和產(chǎn)量。水稻適宜移栽期除了受溫度的直接影響外,還受降水和日照影響。此外,土壤特性、品種更替、災(zāi)害防御等因素也可造成各地移栽期的不同[22]。因此,該階段敏感性是環(huán)境、品種、農(nóng)戶決策等多種因素綜合作用的結(jié)果。

    圖5 雙季早稻、晚稻不同發(fā)育階段相對敏感性與階段日數(shù)變異系數(shù)的相關(guān)性分析

    注: 圖中所有站點(diǎn)相對敏感性均通過0.05水平的顯著性檢驗(yàn)。晚稻移栽-返青及返青-孕穗相對敏感性通過顯著性檢驗(yàn)站點(diǎn)較少,且相關(guān)性不顯著(圖略)

    Note: The relative sensitivity of the stations was significant (P<0.05). There existed few staions whose relative sensitivity was significant (P<0.05) during T-R and R-B of late rice and its relationship was not significant (Figure was omitted)

    3 結(jié)論與討論

    (1)歷史資料統(tǒng)計(jì)分析表明,長江中下游地區(qū)雙季早、晚稻4個(gè)關(guān)鍵發(fā)育階段平均氣溫隨時(shí)間均呈升高趨勢,但增溫特征不均勻。早稻和晚稻的孕穗-成熟期增溫最明顯,氣溫平均升高速率均超過1.0℃·10a?1,水稻移栽前后增暖不明顯。氣候變化背景下,早稻和晚稻的出苗-移栽持續(xù)日數(shù)呈線性減少趨勢,這與同類研究成果基本一致[23]。而早稻和晚稻的返青-孕穗期持續(xù)日數(shù)均呈增加趨勢,這可能是由不同品種生長的最適宜溫度的差異所致。同一溫度條件下,有的品種最適宜生長,生長速度快,階段日數(shù)減少,但同樣溫度條件有可能就超出了其它品種的最適溫度范圍,導(dǎo)致發(fā)育速率隨溫度升高而降低,完成該發(fā)育階段所需日數(shù)反而增加[24]。

    (2)不同作物、不同發(fā)育階段的持續(xù)日數(shù)對溫度變化的相對敏感性存在一定差異,其中早稻返青-孕穗期對溫度變化最敏感,相對敏感性平均為?0.094個(gè)點(diǎn)·℃?1(P<0.05),晚稻孕穗-成熟期最不敏感,相對敏感性平均為?0.045個(gè)點(diǎn)·℃?1(P<0.05)。在同一發(fā)育階段,返青前晚稻更敏感,返青后早稻敏感性更強(qiáng);對于同一作物,除移栽-返青期外,其它營養(yǎng)生長階段日數(shù)對溫度變化的敏感性均強(qiáng)于生殖生長,這與前人研究結(jié)果一致[25]。由此可見,基于相對敏感性可以較準(zhǔn)確地反映雙季早稻、晚稻不同階段日數(shù)對溫度變化的敏感性差異,印證了該方法在比較不同作物、不同發(fā)育階段對溫度變化的敏感性中的可行性。

    (3)雙季早稻、晚稻主要發(fā)育階段相對敏感性與持續(xù)日數(shù)變異系數(shù)呈負(fù)相關(guān)關(guān)系,持續(xù)日數(shù)波動越大,相對敏感性為負(fù)的絕對值越大,階段日數(shù)隨溫度升高而減少越明顯,這種關(guān)系以早稻返青-孕穗期最突出。影響階段日數(shù)變化的因素除了溫度條件,可能還包括品種更替和農(nóng)戶決策等。從災(zāi)害防御實(shí)踐上看,雙季早稻有秧田期低溫和孕穗抽穗期高溫,雙季晚稻有孕穗抽穗期低溫和灌漿成熟期高溫。通過調(diào)整播期可以避免低溫影響,針對高溫危害防御,水稻種植則可移向高海拔山區(qū)或選擇耐高溫型晚熟品種[26],由此造成同一階段日數(shù)的區(qū)域差異。與此同時(shí),階段日數(shù)的不同也反過來影響發(fā)育期溫度條件的區(qū)域差異。因此,環(huán)境、品種、農(nóng)戶決策等多種影響因素的共同相互作用可能是相對敏感性無明顯空間分布特征的原因。

    晚稻返青-孕穗期的相對敏感性統(tǒng)計(jì)結(jié)果表明,該生育階段日數(shù)隨溫度升高而增加的現(xiàn)象在研究區(qū)域不普遍,但與其它發(fā)育階段相比仍較突出,這可能與品種更替適應(yīng)氣候變化有關(guān)。隨著氣候變暖,各地種植的水稻品種逐漸由早熟型向晚熟型方向轉(zhuǎn)變。水稻品種更替可以降低作物營養(yǎng)生長階段的基本早熟性,抵消溫度升高對水稻生育期的負(fù)面影響,穩(wěn)定或延長了水稻的生育期天數(shù)[27],引起該階段日數(shù)對溫度變化敏感性不明顯。對于晚稻移栽-返青期,僅1個(gè)站點(diǎn)相對敏感性通過了顯著性檢驗(yàn),區(qū)域內(nèi)該階段日數(shù)對溫度變化不敏感。在水稻模型ORYZA2000的發(fā)育期模塊中,移栽后水稻生長將停滯一段時(shí)間,期間發(fā)育速率定為0℃·d,該階段秧苗發(fā)育進(jìn)程對溫度變化不敏感[28],本研究結(jié)果與之符合。此外,移栽期受農(nóng)戶決策影響較大,使得移栽-返青期溫度敏感性不顯著且不穩(wěn)定。

    準(zhǔn)確刻畫作物生長發(fā)育對溫度變化的響應(yīng)是提高作物生長模型對產(chǎn)量模擬精準(zhǔn)度的關(guān)鍵。Wang等[29]指出,大多數(shù)作物生長模型對產(chǎn)量模擬的不確定可用發(fā)育速率與溫度變化的關(guān)系來解釋,通過改進(jìn)作物生長發(fā)育對溫度變化的響應(yīng)機(jī)制,產(chǎn)量模擬誤差可明顯降低。光溫水條件、種植制度、品種分布等在時(shí)空上發(fā)生了不同程度的變化,影響著作物生長發(fā)育對氣候變化的響應(yīng)[30]。本研究僅考慮了平均氣溫這一主要環(huán)境因子,尚不能充分反映品種、地形、土壤、氣候等因素對作物發(fā)育期的綜合影響。因此,作物發(fā)育期對氣候變化的敏感性研究還需要結(jié)合多種因子的影響進(jìn)一步完善,為提高作物模型模擬精度和農(nóng)業(yè)生產(chǎn)適應(yīng)能力提供科學(xué)參考。

    [1]IPCC.Climate change 2013:the physical science basis[M]. Cambridge,U K:Cambridge University Press,2013:1-300.

    [2]胡琦,潘學(xué)標(biāo),邵長秀,等.1961-2010年中國農(nóng)業(yè)熱量資源分布和變化特征[J].中國農(nóng)業(yè)氣象,2014,35(2):119-127.

    Hu Q,Pan X B,Shao C X,et al.Distribution and variation of China agricultural heat resources in 1961-2010[J].Chinese Journal of Agrometeorology,2014,35(2):119-127.(in Chinese)

    [3]譚凱炎,鄔定榮,趙花榮.氣候變暖背景下華北平原冬小麥生育期溫度條件變化趨勢分析[J].中國農(nóng)業(yè)氣象,2017,38(6): 333-341.

    Tan K Y,Wu D R,Zhao H R.Trend analysis of temperature conditions over different growth periods of winter wheat under climate warming in North China Plain[J].Chinese Journal of Agrometeorology,2017,38(6):333-341.(in Chinese)

    [4]Zhang T Y,Li T,Yang X G,et al.Model biases in rice phenology under warmer climates[J].Scientific Reports,2016,(6):1-9.

    [5]李正國,楊鵬,唐華俊,等.氣候變化背景下東北三省主要作物典型物候期變化趨勢分析[J].中國農(nóng)業(yè)科學(xué),2011,44(20): 4180-4189.

    Li Z G,Yang P,Tang H J,et al.Trend analysis of typical phenophases of major crops under climate change in the three provinces of Northeast China[J].Scientia Agricultura Sinica, 2011,44(20):4180-4189.(in Chinese)

    [6]薛昌穎,劉榮花,吳騫.氣候變暖對信陽地區(qū)水稻生育期的影響[J].中國農(nóng)業(yè)氣象,2010,31(3):353-357.

    Xue C Y,Liu R H,Wu Q.Effect of climate warming on rice growing stages in Xinyang[J].Chinese Journal of Agromete- orology,2010,31(3):353-357.(in Chinese)

    [7]葛道闊,金之慶.氣候及其變率變化對長江中下游稻區(qū)水稻生產(chǎn)的影響[J].中國水稻科學(xué),2009,23(1):57-64.

    Ge D K,Jin Z Q.Impacts of climate change and its variability on rice production in the middle and lower valley of the Yangtze River,China[J].Chinese Journal of Rice Science,2009, 23(1):57-64.(in Chinese)

    [8]丁麗佳,謝松元.氣候變暖對潮州水稻主要生育期的影響和對策[J].中國農(nóng)業(yè)氣象,2009,30(S1):97-102.

    Ding L J,Xie S Y.Impacts and countermeasures of climate warming on the main growth periods of rice in Chaozhou[J]. Chinese Journal of Agrometeorology,2009, 30(Supp1):97-102. (in Chinese)

    [9]董文軍,鄧艾興,張彬,等.開放式晝夜不同增溫對單季稻影響的試驗(yàn)研究[J].生態(tài)學(xué)報(bào),2011,31(8):2169-2177.

    Dong W J,Deng A X,Zhang B,et al.An experimental study on the effects of different diurnal warming regimes on single cropping rice with free air temperature increased (FATI) facility[J].Acta Ecologica Sinica,2011,31(8):2169-2177.(in Chinese)

    [10]肖登攀,陶福祿,沈彥俊,等.華北平原冬小麥對過去30 年氣候變化響應(yīng)的敏感性研究[J].中國生態(tài)農(nóng)業(yè)學(xué)報(bào),2014, 22(4): 430-438.

    Xiao D P,Tao F L,Shen Y J,et al.Sensitivity of response of winter wheat to climate change in the North China Plain in the last three decades[J].Chinese Journal of Eco-Agriculture, 2014,22(4):430-438.(in Chinese)

    [11]蔣菊芳,王鶴齡,魏育國,等.河西走廊東部不同類型植物物候?qū)夂蜃兓捻憫?yīng)[J].中國農(nóng)業(yè)氣象,2011,32(4):543- 549.

    Jiang J F,Wang H L,Wei Y G,et al.Response of different type plants’ phenology to climate change in east of Hexi Corridor[J].Chinese Journal of Agrometeorology,2011,32(4): 543-549.(in Chinese)

    [12]孟林,劉新建,鄔定榮,等.華北平原夏玉米主要生育期對氣候變化的響應(yīng)[J].中國農(nóng)業(yè)氣象,2015,36(4):375-382.

    Meng L,Liu X J,Wu D R,et al.Responses of summer maize main phenology to climate change in the North China Plain[J]. Chinese Journal of Agrometeorology,2015,36(4): 375-382.(in Chinese)

    [13]高靜,鄔定榮,王培娟,等.華北平原冬小麥主要發(fā)育階段日數(shù)對溫度變化的敏感性分析[J].中國農(nóng)業(yè)氣象,2016,37(4): 431-436.

    Gao J,Wu D R,Wang P J,et al.Relative sensitivity of main growth durations to temperature for winter wheat in the North China Plain[J].Chinese Journal of Agrometeorology, 2016,37(4):431-436.(in Chinese)

    [14]李勇,楊曉光,葉清,等.全球氣候變暖對中國種植制度可能影響Ⅸ:長江中下游地區(qū)單雙季稻高低溫災(zāi)害風(fēng)險(xiǎn)及其產(chǎn)量影響[J].中國農(nóng)業(yè)科學(xué),2013,46(19):3997-4006.

    Li Y,Yang X G,Ye Q,et al.The possible effects of global warming on cropping systems in China Ⅸ:the risk of high and low temperature disasters for single and double rice and its impacts on rice yield in the Middle-Lower Yangtze Plain[J].Scientia Agricultura Sinica,2013,46(19):3997-4006. (in Chinese)

    [15]王春乙,姚蓬娟,張繼權(quán),等.長江中下游地區(qū)雙季早稻冷害、熱害綜合風(fēng)險(xiǎn)評價(jià)[J].中國農(nóng)業(yè)科學(xué),2016,49(13):2469- 2483.

    Wang C Y,Yao P J,Zhang J Q,et al.Risk assessment of cold and hot damages for double-cropping early rice in Lower- Middle Reaches of the Yangtze River Basin[J]. Scientia Agricultura Sinica,2016,49(13):2469-2483.(in Chinese)

    [16]程式華,李建.現(xiàn)代中國水稻[M].北京:金盾出版社,2007: 116-240.

    Cheng S H,Li J.Modern Chinese rice[M].Beijing:Jindun Press,2007:116-240.(in Chinese)

    [17]黃亮,謝小萍,楊沈斌,等.SRES AlB情景下氣候變化對長江中下游稻區(qū)單雙季稻生產(chǎn)的影響[J].江蘇農(nóng)業(yè)科學(xué), 2013,41(9):36-41.

    Huang L,Xie X P,Yang S B,et al.Impacts of climate change on single and double cropping rice production in Lower-Middle reaches of the Yangtze River Basin under SRES AlB scenario[J].Jiangsu Agricultural Sciences,2013, 41(9):36-41.(in Chinese)

    [18]Singh H,Singh K N.Transplating shock in temperate rice and its influence on rooting characteristics and grain yield[J]. Indian J. Agric. Res.,2014,48(5):389-393.

    [19]Zhang Y P,Xiang J,Chen H Z,et al.Physiological characterization of green turning of rice seedlings at different temperatures[J]. Agricultural Science and Technology, 2015, 16(7):1390-1394.

    [20]陳書濤,王讓會,許遐禎,等.氣溫及降水變化對江蘇省典型農(nóng)業(yè)區(qū)冬小麥、水稻生育期的影響[J].中國農(nóng)業(yè)氣象,2011, 32(2):235-239.

    Chen S T,Wang R H,Xu X Z,et al.Impacts of variations in air temperature and precipitation on the growth stages of winter wheat and rice in typical agricultural zones of Jiangsu province[J].Chinese Journal of Agrometeorology,2011,32(2): 235-239.(in Chinese)

    [21]馬偉,霍治國,姬靜華,等.殼聚糖和水楊酸混配液對移栽返青期水稻幼苗抗寒性影響及機(jī)理[J].中國稻米,2015,21(4): 150-154.

    Ma W,Huo Z G,Ji J H,et al.Cold resistance effect and mechanism of mixed carboxymethl chitosan and salicylic acid on rice seedlings during transplanting and rehabilitation period[J].China Rice,2015,21(4):150-154.(in Chinese)

    [22]高蘋,徐敏,徐經(jīng)緯,等.基于大尺度海溫因子的江蘇省水稻適宜移栽期預(yù)報(bào)模型研究[J].氣象,2015,41(12):1538-1544.

    Gao P,Xu M,Xu J W,et al.Study of the prediction model of suitable transplanting date of rice in Jiangsu based on large- scale factors of sea surface temperature[J]. Meteorological Monthly,2015,41(12):1538-1544.(in Chinese)

    [23]楊沈斌,申雙和,趙小艷,等.氣候變化對長江中下游稻區(qū)水稻產(chǎn)量的影響[J].作物學(xué)報(bào),2010,36(9):1519-1528.

    Yang S B,Shen S H,Zhao X Y,et al.Impacts of climate changes on rice production in the middle and lower reaches of the Yangtze River[J].Acta Agronomica Sinica,2010,36(9): 1519-1528.(in Chinese)

    [24]Lu P L,Yu Q,Wang E L,et al.Effects of climatic variation and warming on rice development across South China[J].Climate Research,2008,36:79-88.

    [25]郭建平.氣候變化對中國農(nóng)業(yè)生產(chǎn)的影響研究進(jìn)展[J].應(yīng)用氣象學(xué)報(bào),2015,26(1):1-11.

    Guo J P.Advances in impacts of climate change on agricultural production in China[J].Journal of Applied Meteorological Sciences,2015,26(1):1-11.(in Chinese)

    [26]江敏,金之慶,石春林,等.福建省基于自適應(yīng)調(diào)整的水稻生產(chǎn)對未來氣候變化的響應(yīng)[J].作物學(xué)報(bào),2012,38(12): 2246-2257.

    Jiang M,Jin Z Q,Shi C L,et al.Response of rice production based on self-adaption to climate change in Fujian province[J]. Acta Agronomica Sinica,2012,38(12):2246- 2257.(in Chinese)

    [27]Liu L L,Wang E L,Zhu Y,et al.Contrasting effects of warming and autonomous breeding on single-rice productivity in China[J].Agriculture, Ecosystems and Environment,2012,149: 20-29.

    [28]Bouman B A M,Kropff M J,Tuong T P,et al.ORYZA2000: modeling lowland rice[M].Wageningen,the Netherlands: IRRI,Los Ba?os,Philippines and Wageningen University and Research Center,2001:34-35.

    [29]Wang E L,Martre P,Zhao Z G,et al.The uncertainty of crop yield projections is reduced by improved temperature response functions[J].Nature Plants,2017,(3):1-11.

    [30]周廣勝.氣候變化對中國農(nóng)業(yè)生產(chǎn)影響研究展望[J].氣象與環(huán)境科學(xué),2015,38(1):80-94.

    Zhou G S.Research prospect on impact of climate change on agricultural production in China[J].Meteorological and Environmental Sciences,2015,38(1):80-94.(in Chinese)

    Comparison of Temperature Sensitivity During Different Development Stages for Double-season Early and Late Rice

    WANG Zhi-hai1, JIN Zhi-feng1, WU Ding-rong2, MAO Zhi-jun3, CHEN Zhong-yun4

    (1. Zhejiang Climate Center, Hangzhou 310017, China; 2. Chinese Academy of Meteorological Sciences, Beijing 100081;3. Meteorological Bureau of Longyou, Quzhou 324400; 4. Meteorological Bureau of Huzhou, Huzhou 313000)

    In order to explore responses of different crops’ different development stages to climate change, the relative sensitivities of growth duration to temperature of early and late rice in 38 agrometeorological observation stations during 1981-2010 in the Middle-Lower Yangtze Plain were calculated. Both early and late rice’s growth stages were divided into four stages, i.e., emergence to transplanting(E-T), transplanting to reviving(T-R), reviving to booting(R-B) and booting to mature(B-M). Results showed that growth duration decreased with temperature increasing, but its relative sensitivity varied much between all the four stages. The R-B period of early rice was the most sensible to temperature with average relative sensitivity value of ?0.094 point·℃?1(P<0.05). The B-M period was less sensible than others, with average value of ?0.045 point·℃?1(P<0.05). For the same stage, duration of late rice was more sensible to temperature than early rice before reviving, but was less sensible after reviving. For the two crops, the vegetative growth was more sensible than reproductive growth except T-R period. Further analysis showed that the relative sensitivity was negatively correlated with variation coefficient of duration, indicating in such stations where growth duration varied substantially, the duration would decreased more sharply under warmer condition, especially in the R-B period of early rice, whose correlation coefficient was ?0.761 (P<0.001).

    Warming; Growth duration; Relative sensitivity to temperature; Rice; Middle-Lower Yangtze Plain

    10.3969/j.issn.1000-6362.2018.03.006

    王治海,金志鳳,鄔定榮,等.雙季早晚稻不同發(fā)育階段日數(shù)對溫度變化的敏感性比較[J].中國農(nóng)業(yè)氣象,2018,39(3):185-194

    收稿日期:2017-07-07

    通訊作者。E-mail:wudr@camscma.cn

    浙江省重點(diǎn)科技專項(xiàng)(2015C02048);浙江省氣象局重點(diǎn)專項(xiàng)(2016ZD09);華東區(qū)域氣象科技協(xié)同創(chuàng)新基金(QYHZ201609)

    王治海(1988-),碩士,工程師,主要從事農(nóng)業(yè)氣象服務(wù)相關(guān)工作。E-mail:shan_huhai@163.com

    猜你喜歡
    發(fā)育階段孕穗日數(shù)
    漢江上游漢中區(qū)域不同等級降水日數(shù)的氣候變化特征分析
    綠色科技(2022年16期)2022-09-15 03:04:46
    小麥生殖發(fā)育階段對低溫的敏感性鑒定
    磷酸二氫鉀葉面肥在水稻上的應(yīng)用效果
    江陰市“揚(yáng)麥25”高產(chǎn)高效栽培技術(shù)研究
    簡述寒地水稻科學(xué)施肥技術(shù)
    天津市濱海新區(qū)塘沽地域雷暴日數(shù)變化規(guī)律及特征分析
    天津科技(2020年2期)2020-03-03 05:09:48
    對森工林區(qū)在商品林基地培養(yǎng)速生楊樹探討
    海南省雷暴日數(shù)年代際變化特征
    ESSENTIAL NORMS OF PRODUCTS OF WEIGHTED COMPOSITION OPERATORS AND DIFFERENTIATION OPERATORS BETWEEN BANACH SPACES OF ANALYTIC FUNCTIONS?
    大花黃牡丹葉片發(fā)育過程中氣孔密度和氣孔指數(shù)的動態(tài)變化
    亚洲欧美激情在线| 国产精品98久久久久久宅男小说| 99国产极品粉嫩在线观看| 热re99久久国产66热| 正在播放国产对白刺激| 中文字幕另类日韩欧美亚洲嫩草| 动漫黄色视频在线观看| av天堂久久9| 国产精品自产拍在线观看55亚洲| 丝袜在线中文字幕| 老司机深夜福利视频在线观看| 看片在线看免费视频| 亚洲在线自拍视频| 欧美中文日本在线观看视频| 别揉我奶头~嗯~啊~动态视频| 男人操女人黄网站| 一区在线观看完整版| 免费人成视频x8x8入口观看| 久久人妻熟女aⅴ| 亚洲精品美女久久av网站| 婷婷六月久久综合丁香| 国产一区在线观看成人免费| 又紧又爽又黄一区二区| 国产伦一二天堂av在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 久9热在线精品视频| 五月开心婷婷网| 午夜a级毛片| 国产三级在线视频| 日本三级黄在线观看| 91国产中文字幕| 日韩大码丰满熟妇| 91精品三级在线观看| 亚洲成人精品中文字幕电影 | 99国产精品一区二区三区| 美国免费a级毛片| 午夜日韩欧美国产| 久久青草综合色| 91大片在线观看| 高清在线国产一区| 国产片内射在线| 亚洲国产中文字幕在线视频| 一进一出抽搐动态| 亚洲免费av在线视频| 成年女人毛片免费观看观看9| 久久九九热精品免费| av网站免费在线观看视频| 久久人妻熟女aⅴ| 久久亚洲真实| 国产午夜精品久久久久久| 十八禁网站免费在线| 国产精品国产av在线观看| 免费少妇av软件| 韩国精品一区二区三区| 亚洲专区中文字幕在线| 久久精品91蜜桃| 国产精品1区2区在线观看.| 精品国产国语对白av| 人妻丰满熟妇av一区二区三区| 最近最新中文字幕大全电影3 | 国产成+人综合+亚洲专区| 免费看a级黄色片| 视频在线观看一区二区三区| 中文字幕人妻熟女乱码| 视频区欧美日本亚洲| 少妇被粗大的猛进出69影院| 天天影视国产精品| 在线视频色国产色| 黑人巨大精品欧美一区二区蜜桃| 精品久久久久久电影网| 黄色女人牲交| 国产亚洲精品第一综合不卡| 99久久人妻综合| 午夜视频精品福利| 一区二区三区激情视频| 一级a爱视频在线免费观看| 在线观看舔阴道视频| 精品国内亚洲2022精品成人| 人妻久久中文字幕网| 人妻久久中文字幕网| 国产麻豆69| 欧美日韩福利视频一区二区| 母亲3免费完整高清在线观看| 国产精品免费一区二区三区在线| 99re在线观看精品视频| 欧美日韩黄片免| 久久精品亚洲熟妇少妇任你| 搡老熟女国产l中国老女人| 久久国产亚洲av麻豆专区| 人妻丰满熟妇av一区二区三区| 欧美久久黑人一区二区| 国产成人啪精品午夜网站| 女人被狂操c到高潮| 91大片在线观看| 国产成人一区二区三区免费视频网站| 黄色a级毛片大全视频| 熟女少妇亚洲综合色aaa.| 国产成人一区二区三区免费视频网站| 真人一进一出gif抽搐免费| 黄片大片在线免费观看| 母亲3免费完整高清在线观看| 中文字幕另类日韩欧美亚洲嫩草| 欧美激情极品国产一区二区三区| 99精品欧美一区二区三区四区| 亚洲精品一二三| videosex国产| 狠狠狠狠99中文字幕| 一二三四社区在线视频社区8| 久久精品aⅴ一区二区三区四区| 99国产极品粉嫩在线观看| 亚洲少妇的诱惑av| 日韩免费av在线播放| 99久久久亚洲精品蜜臀av| 一区二区日韩欧美中文字幕| 国产1区2区3区精品| 淫妇啪啪啪对白视频| 91大片在线观看| 亚洲国产欧美网| 脱女人内裤的视频| 黄色毛片三级朝国网站| 淫妇啪啪啪对白视频| 久久国产亚洲av麻豆专区| 久久久久久久久中文| 亚洲精品av麻豆狂野| 成人精品一区二区免费| 9191精品国产免费久久| 国产乱人伦免费视频| 日韩免费av在线播放| 久久国产精品男人的天堂亚洲| 欧美在线黄色| 黄色视频不卡| 久久人妻熟女aⅴ| 超色免费av| 久久久久久亚洲精品国产蜜桃av| 日韩大尺度精品在线看网址 | 免费久久久久久久精品成人欧美视频| 午夜久久久在线观看| 久久国产精品人妻蜜桃| 国产精品电影一区二区三区| 黄色毛片三级朝国网站| 高清欧美精品videossex| 亚洲三区欧美一区| 男人舔女人的私密视频| av在线播放免费不卡| 人人妻人人爽人人添夜夜欢视频| 美女高潮到喷水免费观看| 亚洲熟妇熟女久久| 国产成人精品久久二区二区91| 日韩一卡2卡3卡4卡2021年| 国产激情久久老熟女| 精品国产乱子伦一区二区三区| av超薄肉色丝袜交足视频| 美女高潮到喷水免费观看| 国产成人av教育| 黄片播放在线免费| 国产av精品麻豆| 国产亚洲精品久久久久5区| 国产亚洲欧美98| 一级毛片高清免费大全| 精品福利观看| 人人妻,人人澡人人爽秒播| 长腿黑丝高跟| 午夜亚洲福利在线播放| 91在线观看av| 在线观看免费高清a一片| 国产精品国产高清国产av| 日韩中文字幕欧美一区二区| 亚洲专区中文字幕在线| 欧美在线黄色| 亚洲精品粉嫩美女一区| 国产欧美日韩精品亚洲av| 精品免费久久久久久久清纯| 欧美乱妇无乱码| 1024视频免费在线观看| 日本a在线网址| 欧美成狂野欧美在线观看| 成人黄色视频免费在线看| 免费女性裸体啪啪无遮挡网站| 男女下面插进去视频免费观看| 一二三四社区在线视频社区8| 欧美激情高清一区二区三区| 深夜精品福利| 美女午夜性视频免费| 精品国产乱码久久久久久男人| 大香蕉久久成人网| 国产成人欧美在线观看| xxxhd国产人妻xxx| 久久久久国产精品人妻aⅴ院| 男女下面插进去视频免费观看| 后天国语完整版免费观看| aaaaa片日本免费| 亚洲一区二区三区不卡视频| 999久久久精品免费观看国产| 在线观看一区二区三区激情| 一二三四社区在线视频社区8| 在线看a的网站| 99久久精品国产亚洲精品| 丁香欧美五月| 五月开心婷婷网| 这个男人来自地球电影免费观看| 女人精品久久久久毛片| 高潮久久久久久久久久久不卡| 欧美日韩av久久| 搡老岳熟女国产| 极品教师在线免费播放| 无人区码免费观看不卡| 可以在线观看毛片的网站| 久久精品91蜜桃| 久久精品aⅴ一区二区三区四区| 极品教师在线免费播放| 久久香蕉激情| 最近最新中文字幕大全电影3 | 老熟妇乱子伦视频在线观看| 大型黄色视频在线免费观看| 亚洲精品一二三| 丝袜人妻中文字幕| aaaaa片日本免费| av国产精品久久久久影院| 变态另类成人亚洲欧美熟女 | 天天躁夜夜躁狠狠躁躁| 久久亚洲精品不卡| 亚洲av第一区精品v没综合| 亚洲精品国产色婷婷电影| 国产日韩一区二区三区精品不卡| 久久久久国产精品人妻aⅴ院| 国产精品影院久久| 国产午夜精品久久久久久| 日韩大码丰满熟妇| 国内久久婷婷六月综合欲色啪| 国产伦一二天堂av在线观看| 久久久久久免费高清国产稀缺| 亚洲第一av免费看| av在线天堂中文字幕 | 国产精品一区二区在线不卡| 老司机深夜福利视频在线观看| 波多野结衣高清无吗| 亚洲中文av在线| 国产极品粉嫩免费观看在线| 欧美大码av| 国产精品国产高清国产av| 91成人精品电影| 成人免费观看视频高清| 后天国语完整版免费观看| 黄色丝袜av网址大全| 欧美日韩精品网址| 在线观看一区二区三区| 国产精品自产拍在线观看55亚洲| 日本一区二区免费在线视频| 美女福利国产在线| 12—13女人毛片做爰片一| 久久中文看片网| 亚洲aⅴ乱码一区二区在线播放 | 久久久国产一区二区| 淫妇啪啪啪对白视频| 亚洲专区国产一区二区| 村上凉子中文字幕在线| 欧美激情 高清一区二区三区| 啪啪无遮挡十八禁网站| 女人精品久久久久毛片| 午夜91福利影院| 欧美+亚洲+日韩+国产| 日韩精品免费视频一区二区三区| 淫秽高清视频在线观看| x7x7x7水蜜桃| 久久精品国产综合久久久| 99国产极品粉嫩在线观看| 国产一区二区激情短视频| 操出白浆在线播放| 中文字幕高清在线视频| 一边摸一边做爽爽视频免费| 在线观看免费视频网站a站| 交换朋友夫妻互换小说| 精品久久久久久成人av| 9热在线视频观看99| 丁香欧美五月| 超色免费av| 99精品在免费线老司机午夜| 亚洲人成网站在线播放欧美日韩| 男女床上黄色一级片免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 久久久久久久久中文| 在线观看www视频免费| 欧美日韩亚洲综合一区二区三区_| 成熟少妇高潮喷水视频| 亚洲精品成人av观看孕妇| 午夜免费观看网址| 后天国语完整版免费观看| 精品国产国语对白av| 午夜影院日韩av| 精品福利永久在线观看| 丰满人妻熟妇乱又伦精品不卡| 变态另类成人亚洲欧美熟女 | 99国产极品粉嫩在线观看| 久久香蕉国产精品| 女警被强在线播放| 91在线观看av| 亚洲精品在线观看二区| av在线天堂中文字幕 | 麻豆一二三区av精品| 精品国产超薄肉色丝袜足j| 男女做爰动态图高潮gif福利片 | 一区福利在线观看| 中文字幕另类日韩欧美亚洲嫩草| 久久久精品国产亚洲av高清涩受| 丝袜美足系列| 久久香蕉精品热| 久久天堂一区二区三区四区| 又大又爽又粗| 国产精品秋霞免费鲁丝片| 国产日韩一区二区三区精品不卡| 手机成人av网站| 黄色女人牲交| 嫁个100分男人电影在线观看| 久久热在线av| 亚洲成人久久性| 十八禁网站免费在线| 欧美人与性动交α欧美精品济南到| 成人18禁高潮啪啪吃奶动态图| 丝袜在线中文字幕| 1024香蕉在线观看| 制服诱惑二区| 脱女人内裤的视频| 午夜视频精品福利| 91成年电影在线观看| 国产区一区二久久| 日韩国内少妇激情av| 精品国产一区二区三区四区第35| 中文字幕人妻丝袜制服| 亚洲一区高清亚洲精品| 曰老女人黄片| 国产免费av片在线观看野外av| 亚洲精华国产精华精| 正在播放国产对白刺激| 久热爱精品视频在线9| 亚洲自拍偷在线| 欧美成人性av电影在线观看| 成人亚洲精品av一区二区 | 精品免费久久久久久久清纯| 十分钟在线观看高清视频www| 日韩免费高清中文字幕av| 丁香欧美五月| 中亚洲国语对白在线视频| 久久国产精品人妻蜜桃| 亚洲全国av大片| 精品久久蜜臀av无| 一本综合久久免费| 在线永久观看黄色视频| 一区二区三区精品91| 热re99久久国产66热| 一区福利在线观看| 一级作爱视频免费观看| 免费日韩欧美在线观看| 久久精品影院6| 成人三级做爰电影| 麻豆成人av在线观看| 一个人免费在线观看的高清视频| 成人三级做爰电影| 丝袜在线中文字幕| 亚洲av第一区精品v没综合| 最新美女视频免费是黄的| 免费看a级黄色片| 午夜福利在线观看吧| 亚洲精品久久午夜乱码| 涩涩av久久男人的天堂| 亚洲五月婷婷丁香| 午夜日韩欧美国产| 久热这里只有精品99| 99香蕉大伊视频| 超碰成人久久| 在线免费观看的www视频| 国产精品美女特级片免费视频播放器 | 岛国在线观看网站| 亚洲av美国av| 两人在一起打扑克的视频| 老熟妇乱子伦视频在线观看| 国产精品国产av在线观看| 极品教师在线免费播放| 色婷婷av一区二区三区视频| 亚洲欧美精品综合一区二区三区| 成年人免费黄色播放视频| 99热国产这里只有精品6| 在线播放国产精品三级| 国产无遮挡羞羞视频在线观看| 亚洲熟妇中文字幕五十中出 | 变态另类成人亚洲欧美熟女 | 叶爱在线成人免费视频播放| 男女高潮啪啪啪动态图| 久久香蕉精品热| 免费高清在线观看日韩| 久久影院123| 美女午夜性视频免费| x7x7x7水蜜桃| 丰满饥渴人妻一区二区三| 90打野战视频偷拍视频| 多毛熟女@视频| 国产深夜福利视频在线观看| 国产黄色免费在线视频| 亚洲欧美一区二区三区久久| 少妇被粗大的猛进出69影院| 视频区图区小说| 精品少妇一区二区三区视频日本电影| 琪琪午夜伦伦电影理论片6080| 午夜成年电影在线免费观看| 嫩草影视91久久| 午夜精品在线福利| a级毛片黄视频| av在线天堂中文字幕 | 黄色成人免费大全| 国产免费现黄频在线看| 成熟少妇高潮喷水视频| 韩国精品一区二区三区| 电影成人av| 国产精品乱码一区二三区的特点 | 亚洲熟妇中文字幕五十中出 | 日韩精品免费视频一区二区三区| 一区在线观看完整版| 老司机在亚洲福利影院| 日韩视频一区二区在线观看| 搡老乐熟女国产| 日韩 欧美 亚洲 中文字幕| 亚洲人成电影免费在线| 成年版毛片免费区| 久久精品人人爽人人爽视色| 国产高清国产精品国产三级| 最好的美女福利视频网| 桃红色精品国产亚洲av| 亚洲专区字幕在线| 久久久水蜜桃国产精品网| www国产在线视频色| 韩国精品一区二区三区| 两性夫妻黄色片| 午夜激情av网站| 免费女性裸体啪啪无遮挡网站| 高潮久久久久久久久久久不卡| 国产xxxxx性猛交| 国产亚洲欧美98| 天天躁夜夜躁狠狠躁躁| 男人的好看免费观看在线视频 | 妹子高潮喷水视频| 性色av乱码一区二区三区2| 久久青草综合色| 丰满的人妻完整版| 国产单亲对白刺激| 自线自在国产av| 国产av又大| 中文字幕精品免费在线观看视频| 亚洲专区字幕在线| 欧美国产精品va在线观看不卡| 长腿黑丝高跟| www.熟女人妻精品国产| 亚洲自拍偷在线| 精品熟女少妇八av免费久了| 亚洲三区欧美一区| 男女之事视频高清在线观看| 丰满迷人的少妇在线观看| 搡老岳熟女国产| 国产成人一区二区三区免费视频网站| 久久香蕉精品热| 一个人免费在线观看的高清视频| 精品国产一区二区三区四区第35| 国产成人一区二区三区免费视频网站| 麻豆成人av在线观看| 精品国产乱码久久久久久男人| 国产成人免费无遮挡视频| 在线永久观看黄色视频| 91麻豆av在线| 午夜免费成人在线视频| 香蕉国产在线看| 不卡av一区二区三区| 99精国产麻豆久久婷婷| 美女福利国产在线| 99在线人妻在线中文字幕| 国产精品一区二区在线不卡| av在线播放免费不卡| 美女 人体艺术 gogo| 欧美亚洲日本最大视频资源| 最近最新免费中文字幕在线| 一本大道久久a久久精品| videosex国产| 国产成人av激情在线播放| 一区福利在线观看| 女人精品久久久久毛片| 亚洲成人久久性| 99riav亚洲国产免费| 欧美日韩av久久| 满18在线观看网站| 国产伦一二天堂av在线观看| 国产成人精品在线电影| 啦啦啦在线免费观看视频4| www.999成人在线观看| 丰满人妻熟妇乱又伦精品不卡| 精品一区二区三卡| 色综合欧美亚洲国产小说| 日韩精品青青久久久久久| 99精品久久久久人妻精品| 动漫黄色视频在线观看| 久久精品国产亚洲av高清一级| 成人精品一区二区免费| 深夜精品福利| 久久人妻av系列| 国产精品免费一区二区三区在线| 多毛熟女@视频| 国产又色又爽无遮挡免费看| 在线av久久热| 国产野战对白在线观看| bbb黄色大片| 亚洲一区二区三区色噜噜 | 夫妻午夜视频| 无遮挡黄片免费观看| 日本wwww免费看| 亚洲成人久久性| 欧美色视频一区免费| 精品国产亚洲在线| 丝袜人妻中文字幕| 国产视频一区二区在线看| 中文欧美无线码| 国产99白浆流出| 久久九九热精品免费| 女同久久另类99精品国产91| 18美女黄网站色大片免费观看| 757午夜福利合集在线观看| 午夜福利在线免费观看网站| 很黄的视频免费| 女人高潮潮喷娇喘18禁视频| 免费在线观看亚洲国产| 成人国产一区最新在线观看| 在线观看日韩欧美| 久久人人97超碰香蕉20202| 久99久视频精品免费| 久久久久久人人人人人| 69精品国产乱码久久久| 国产欧美日韩一区二区三区在线| 别揉我奶头~嗯~啊~动态视频| 国产高清视频在线播放一区| 后天国语完整版免费观看| 不卡一级毛片| 91国产中文字幕| 丰满饥渴人妻一区二区三| 丝袜人妻中文字幕| 亚洲欧美激情在线| 黄片大片在线免费观看| 激情视频va一区二区三区| 丝袜人妻中文字幕| 一边摸一边做爽爽视频免费| 久久精品亚洲av国产电影网| 九色亚洲精品在线播放| 一边摸一边抽搐一进一小说| 成人特级黄色片久久久久久久| 两性夫妻黄色片| 757午夜福利合集在线观看| 午夜福利在线观看吧| 亚洲va日本ⅴa欧美va伊人久久| 欧美日韩亚洲综合一区二区三区_| 搡老熟女国产l中国老女人| 成在线人永久免费视频| 一本大道久久a久久精品| 亚洲av熟女| 女性被躁到高潮视频| 亚洲一码二码三码区别大吗| 日韩国内少妇激情av| 男女下面插进去视频免费观看| 欧美日韩精品网址| 一个人观看的视频www高清免费观看 | 久久影院123| 午夜a级毛片| 精品一区二区三区四区五区乱码| 在线观看免费视频日本深夜| 免费女性裸体啪啪无遮挡网站| 亚洲 欧美一区二区三区| avwww免费| 亚洲第一欧美日韩一区二区三区| 亚洲精品中文字幕在线视频| av超薄肉色丝袜交足视频| 亚洲中文日韩欧美视频| 88av欧美| 伦理电影免费视频| 久久久国产欧美日韩av| 亚洲人成电影免费在线| 黄网站色视频无遮挡免费观看| 最近最新中文字幕大全免费视频| 中文字幕av电影在线播放| 欧美午夜高清在线| 欧美黄色淫秽网站| 好男人电影高清在线观看| 夜夜爽天天搞| 欧美乱码精品一区二区三区| 十八禁网站免费在线| 免费av中文字幕在线| 午夜免费成人在线视频| а√天堂www在线а√下载| 极品教师在线免费播放| 最近最新中文字幕大全电影3 | 成年人免费黄色播放视频| 人人澡人人妻人| 日本精品一区二区三区蜜桃| 午夜成年电影在线免费观看| 日韩 欧美 亚洲 中文字幕| 在线av久久热| 久久久国产欧美日韩av| 色婷婷av一区二区三区视频| 日韩免费av在线播放| 久9热在线精品视频| 在线观看www视频免费| 一边摸一边抽搐一进一小说| 精品人妻在线不人妻| 国内久久婷婷六月综合欲色啪| 久久精品aⅴ一区二区三区四区| 国产成人欧美在线观看| 国产野战对白在线观看| 亚洲中文日韩欧美视频| 亚洲色图av天堂| 中文字幕另类日韩欧美亚洲嫩草| 91在线观看av| 亚洲精品中文字幕一二三四区| 色老头精品视频在线观看|