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

    華北平原冬小麥主要發(fā)育階段日數(shù)對(duì)溫度變化的敏感性分析*

    2016-11-16 09:02:43靜鄔定榮王培娟陳京華閆峰趙煜飛王佳強(qiáng)
    中國(guó)農(nóng)業(yè)氣象 2016年4期
    關(guān)鍵詞:發(fā)育階段發(fā)育期日數(shù)

    高 靜鄔定榮王培娟陳京華閆 峰趙煜飛王佳強(qiáng)

    ?

    華北平原冬小麥主要發(fā)育階段日數(shù)對(duì)溫度變化的敏感性分析*

    高 靜1,鄔定榮2,3**,王培娟2,3,陳京華1,閆 峰4,趙煜飛1,王佳強(qiáng)1

    (1.國(guó)家氣象信息中心,北京 100081;2.中國(guó)氣象科學(xué)研究院,北京 100081;3.中國(guó)氣象科學(xué)研究院災(zāi)害天氣國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100081;4.河北省易縣氣象局,保定 074200)

    氣候變化背景下,作物發(fā)育階段日數(shù)發(fā)生了明顯變化,其對(duì)溫度變化的敏感性也成為氣候變化對(duì)農(nóng)業(yè)影響研究的重要內(nèi)容。針對(duì)現(xiàn)有溫度敏感性研究的不足,本文提出相對(duì)敏感性的概念,并以此分析華北平原冬小麥敏感性的時(shí)空分布特征。利用65個(gè)農(nóng)業(yè)氣象觀測(cè)站1980-2012年冬小麥發(fā)育期觀測(cè)資料及同期逐日平均氣溫?cái)?shù)據(jù),將發(fā)育期分為播種-越冬、越冬-返青、返青-抽穗和抽穗-成熟4個(gè)主要發(fā)育階段,采用線性回歸方法計(jì)算各階段發(fā)育期日數(shù)對(duì)溫度變化的相對(duì)敏感性,并用GIS空間插值方法研究相應(yīng)指標(biāo)的時(shí)空變化特征。結(jié)果表明:(1)各發(fā)育階段平均溫度均呈上升趨勢(shì),不同發(fā)育階段日數(shù)具有明顯的區(qū)域分布規(guī)律;(2)播種-越冬期對(duì)溫度的敏感性為-0.113~0.029℃-1,區(qū)域平均-0.040℃-1,該期敏感性呈中間高兩端低的區(qū)域分布特征,且敏感性具有一定的普遍性;(3)越冬-返青期敏感性為-0.081~0.091℃-1,平均0.013℃-1。該期敏感性區(qū)域波動(dòng)大,但無(wú)明顯區(qū)域特征,敏感性的穩(wěn)定程度也較弱;(4)返青-抽穗期敏感性在-0.112~-0.035℃-1,平均-0.074℃-1。敏感性無(wú)明顯區(qū)域特征,但敏感性的穩(wěn)定性高;(5)抽穗-成熟期敏感性在-0.114~0.014℃-1,平均-0.042℃-1,呈顯著的南高北低的空間分布特征,在區(qū)域上具有很強(qiáng)的穩(wěn)定性??傮w上,不同發(fā)育階段對(duì)溫度的敏感性差異較大,且敏感性具有明顯的區(qū)域分布特征。

    氣候變化;發(fā)育階段日數(shù);溫度的相對(duì)敏感性;區(qū)域分布規(guī)律;時(shí)間分布特征

    氣候變化對(duì)農(nóng)業(yè)的影響一直是農(nóng)業(yè)領(lǐng)域國(guó)內(nèi)外研究關(guān)注的焦點(diǎn)[1]。作物物候?qū)夂蜃兓捻憫?yīng)是農(nóng)業(yè)應(yīng)對(duì)氣候變化的核心研究?jī)?nèi)容之一。近幾十年來(lái),為準(zhǔn)確評(píng)估氣候變化對(duì)農(nóng)業(yè)生產(chǎn)的影響,國(guó)內(nèi)外學(xué)者對(duì)歷史氣象條件下作物發(fā)育期對(duì)氣候變化的響應(yīng)進(jìn)行了廣泛研究。研究表明,氣候變暖使國(guó)外多種作物的發(fā)育期呈提前趨勢(shì)[2]。在中國(guó)華北平原,夏玉米生育期天數(shù)增加速率為0.34d·a-1[3]。在東北,除抽穗期及成熟期外, 水稻各生育期均隨著年代的增加而提前[4]。河西走廊春小麥各發(fā)育期提前明顯,其中抽穗期提前0.36d·a-1,開(kāi)花期提前0.41d·a-1,成熟期提前0.36d·a-1,而玉米在春夏季發(fā)育期提前,秋季乳熟期、成熟期的變化趨勢(shì)不明顯[5]。在關(guān)注作物發(fā)育期隨年份變化的同時(shí)[6-13],國(guó)內(nèi)外學(xué)者重點(diǎn)關(guān)注了發(fā)育期對(duì)溫度的敏感性[13-15]。德國(guó)的研究表明,溫度每升高1℃,燕麥成熟期提前約4d[14]。在國(guó)內(nèi),過(guò)去30a的氣候變暖對(duì)華北平原的作物物候產(chǎn)生了重要影響,發(fā)育期內(nèi)平均溫度每升高1℃,冬小麥返青-開(kāi)花期和開(kāi)花-成熟期分別縮短3.8d和1.2d[15-16],夏玉米全生育期和生殖生長(zhǎng)期天數(shù)分別平均縮短2.71d 和1.07d[3]。

    上述研究極大提升了對(duì)農(nóng)作物發(fā)育期與氣候變化及品種改進(jìn)關(guān)系的理解,為氣候變化背景下農(nóng)業(yè)適應(yīng)的對(duì)策研究提供了堅(jiān)實(shí)的科學(xué)基礎(chǔ)。然而,現(xiàn)有研究的結(jié)果多以發(fā)育期日數(shù)隨溫度的絕對(duì)變化速率來(lái)表達(dá),即溫度每上升1℃作物發(fā)育階段日數(shù)的變化天數(shù)[12,15-16]。但絕對(duì)變化速率受該發(fā)育階段長(zhǎng)度的影響,長(zhǎng)度越長(zhǎng),變化速率越大,因此,這種表達(dá)方式不適合在不同作物、不同站點(diǎn)和不同發(fā)育期之間進(jìn)行發(fā)育速率對(duì)溫度敏感性的比較。

    基于此,本文提出相對(duì)敏感性的概念,即將發(fā)育階段的日數(shù)除以該階段多年平均日數(shù),在此基礎(chǔ)上計(jì)算其與同期平均溫度的關(guān)系,并以華北平原(北京、天津、河北、山東、河南)多個(gè)冬小麥農(nóng)業(yè)氣象試驗(yàn)站點(diǎn)冬小麥發(fā)育期和逐日氣象觀測(cè)資料為基礎(chǔ),在站點(diǎn)和區(qū)域尺度上研究冬小麥不同發(fā)育階段日數(shù)對(duì)溫度的相對(duì)敏感性及敏感性的區(qū)域差異,以期為氣候變化背景下作物產(chǎn)量響應(yīng)的區(qū)域模擬研究提供科學(xué)依據(jù)。

    1 資料與方法

    1.1 資料

    選用1980-2012年華北平原冬小麥連續(xù)種植超過(guò)10a的65個(gè)農(nóng)業(yè)氣象觀測(cè)站的資料進(jìn)行分析(圖1)。冬小麥生育期觀測(cè)資料與氣象資料均來(lái)自國(guó)家氣象信息中心。發(fā)育期資料包括播種、開(kāi)始越冬、返青、抽穗和成熟等5個(gè)發(fā)育期。氣象資料為與農(nóng)業(yè)氣象觀測(cè)站對(duì)應(yīng)的65個(gè)氣象臺(tái)站1980-2012年的逐日平均氣溫?cái)?shù)據(jù)。

    1.2 研究方法

    與冬小麥5個(gè)發(fā)育期對(duì)應(yīng)的是4個(gè)發(fā)育階段,分別為播種-越冬(S-O,Sowing to Overwintering)、越冬-返青(O-T,Overwintering to Turning green)、返青-抽穗(T-H,Turning green to Heading)和抽穗-成熟(H-M,Heading to Maturity)。按此劃分方式,分別計(jì)算各站歷年各發(fā)育階段的日數(shù)和同期平均溫度,采用線性回歸方法計(jì)算日數(shù)和平均溫度隨年份變化的趨勢(shì),并進(jìn)行顯著性檢驗(yàn)[17]。

    以往研究中,發(fā)育階段對(duì)溫度的敏感性通過(guò)計(jì)算發(fā)育階段日數(shù)與平均溫度的回歸關(guān)系得到[3, 14-16],即

    式中,D為該站完成某發(fā)育階段所需的日數(shù)(d),T為同期的平均溫度(℃),a和b為回歸系數(shù)。b值代表發(fā)育期內(nèi)平均溫度每升高1oC時(shí)發(fā)育階段變化的天數(shù)。當(dāng)發(fā)育階段經(jīng)歷的天數(shù)較多時(shí)(如返青-成熟),作物暴露在溫度升高環(huán)境中的時(shí)間更長(zhǎng),此時(shí)b值會(huì)比天數(shù)更短的發(fā)育階段(如返青-拔節(jié))更大。因此,式(1)中的b值不能用于不同發(fā)育階段日數(shù)溫度敏感性的比較。由于不同站點(diǎn)同一發(fā)育階段的日數(shù)也不同,因此也不能進(jìn)行不同站點(diǎn)同一發(fā)育階段敏感性的比較。鑒于此,本文將D值除以日數(shù)的平均值,提出發(fā)育階段日數(shù)受溫度影響相對(duì)速率的概念,以便在不同發(fā)育階段之間進(jìn)行比較,公式為

    式中,Daver是該站完成某發(fā)育階段所需日數(shù)的多年平均值(d)。式(2)中的b值即為發(fā)育階段日數(shù)受溫度影響的相對(duì)變化量,可用于不同發(fā)育階段日數(shù)對(duì)溫度敏感性和不同區(qū)域?qū)ν话l(fā)育階段敏感性的比較。

    采用點(diǎn)面結(jié)合的方法分析發(fā)育階段日數(shù)對(duì)溫度的相對(duì)敏感性。首先計(jì)算各站4個(gè)主要發(fā)育階段日數(shù)及同期平均溫度,求其與年份的線性回歸系數(shù)和顯著性水平。在此基礎(chǔ)上,統(tǒng)計(jì)成區(qū)域平均值,進(jìn)行發(fā)育階段日數(shù)和溫度變化的背景分析。隨后根據(jù)式(2)計(jì)算各站4個(gè)主要發(fā)育階段日數(shù)對(duì)溫度的敏感性,并利用ArcGIS空間分析工具中的反距離權(quán)重法進(jìn)行敏感性區(qū)域分布規(guī)律的展示和分析。

    2 結(jié)果與分析

    2.1 發(fā)育階段平均溫度變化趨勢(shì)

    區(qū)域差值結(jié)果表明,不同發(fā)育階段的溫度差異較大,同一發(fā)育階段不同站點(diǎn)的溫度差異也很大。其中,越冬-返青期的平均溫度呈顯著的南高北低的趨勢(shì),其余3個(gè)發(fā)育階段則相反,呈南低北高的趨勢(shì)。

    表1是對(duì)不同發(fā)育階段平均溫度特征的統(tǒng)計(jì)結(jié)果。由表可見(jiàn),氣候變暖背景下,全區(qū)發(fā)育階段的溫度均呈上升趨勢(shì),但不同發(fā)育階段上升速率不同。其中播種-越冬的升溫速率最低,僅為0.03℃·10a-1,抽穗-成熟期的升溫速率為0.17℃·10a-1,越冬-返青期升溫速率較大,為0.35℃·10a-1,返青-抽穗期升溫速率最大,達(dá)0.61℃·10a-1。但上述升溫速率通過(guò)顯著性檢驗(yàn)的站點(diǎn)并不多,播種-越冬期僅有5個(gè)站點(diǎn)的升溫速率達(dá)顯著性水平,增溫最顯著的返青-抽穗期也僅有22個(gè)站點(diǎn)達(dá)顯著性水平(P<0.05)。

    表1 四個(gè)發(fā)育階段平均溫度及其隨年份變化趨勢(shì)的區(qū)域平均值

    2.2 發(fā)育階段日數(shù)的變化特征

    表2是對(duì)不同發(fā)育階段日數(shù)特征的統(tǒng)計(jì),由表可見(jiàn),不同發(fā)育階段日數(shù)相差較大。越冬-返青期最長(zhǎng),達(dá)70.0d,播種-越冬期次之,抽穗-成熟期最短,僅41.1d。不同發(fā)育階段日數(shù)多具有顯著的區(qū)域分布特征(圖略)。播種-越冬期日數(shù)呈由南向北遞減的趨勢(shì),越冬-返青期日數(shù)趨勢(shì)則相反,即由南向北遞增,返青-抽穗期日數(shù)呈南北兩端低中間高的趨勢(shì),抽穗-成熟期日數(shù)呈由南向北遞增減的趨勢(shì)。

    由表2還可見(jiàn),研究時(shí)段內(nèi),營(yíng)養(yǎng)生長(zhǎng)期日數(shù)均呈現(xiàn)隨年份縮短的趨勢(shì)。其中返青-抽穗期縮短較多,達(dá)2.98d·10a-1,但這種回歸關(guān)系僅6個(gè)站點(diǎn)通過(guò)P<0.05水平的顯著性檢驗(yàn),僅占研究區(qū)域站點(diǎn)總數(shù)的9.23%。越冬-返青期縮短速率次之,而播種-越冬期日數(shù)變化不大。生殖生長(zhǎng)期日數(shù)隨年份有延長(zhǎng)的趨勢(shì),延長(zhǎng)速率達(dá)1.39d·10a-1,且有近半數(shù)的站點(diǎn)通過(guò)P<0.05水平的顯著性檢驗(yàn),其中14個(gè)站的顯著性水平達(dá)到P<0.001,表明生殖生長(zhǎng)期日數(shù)隨年份的增加而延長(zhǎng)的趨勢(shì)在全研究區(qū)域普遍存在。

    表2 四個(gè)發(fā)育階段日數(shù)及其隨年份變化趨勢(shì)的區(qū)域平均值

    2.3 發(fā)育階段日數(shù)對(duì)溫度變化的相對(duì)敏感性

    2.3.1 不同發(fā)育階段敏感性的統(tǒng)計(jì)

    由圖2可見(jiàn),不同發(fā)育階段的日數(shù)對(duì)溫度變化的相對(duì)敏感性各不相同,且同一發(fā)育期不同站點(diǎn)的敏感性也具有較大波動(dòng)。除越冬期外,其它發(fā)育階段日數(shù)均隨溫度的上升而縮短。其中,播種-越冬期相對(duì)敏感性的區(qū)域平均值為-0.040℃-1,即溫度每升高1℃,日數(shù)縮短4%。返青-抽穗期是冬小麥營(yíng)養(yǎng)生長(zhǎng)的主要階段,各站對(duì)溫度均較為敏感,區(qū)域平均值為-0.074℃-1,即溫度每升高1℃,日數(shù)縮短7.4%。抽穗-成熟期為冬小麥的生殖生長(zhǎng)階段,該階段對(duì)溫度較不敏感,溫度每升高1℃,日數(shù)縮短4.2%。

    總體上,返青-抽穗期對(duì)溫度最敏感,其次是抽穗-成熟期。播種-越冬期對(duì)溫度較不敏感。以往基于絕對(duì)速率的敏感性研究多表明營(yíng)養(yǎng)生長(zhǎng)期對(duì)溫度的敏感性強(qiáng)于生殖生長(zhǎng)期。如據(jù)報(bào)道,平均溫度每升高1℃,華北平原冬小麥返青-開(kāi)花期和開(kāi)花-成熟期分別縮短3.8d和1.2d[15]。本文基于相對(duì)速率計(jì)算的二者的敏感性分別為-0.074℃-1和-0.042℃-1,雖然營(yíng)養(yǎng)生長(zhǎng)期的敏感性仍強(qiáng)于生殖生長(zhǎng)期,但二者差距大為縮小。

    注:分別代表最大值和最小值;分別代表第99和第1分位數(shù);框圖中的上、中、下三橫分別為上四分位數(shù)、中位數(shù)和下四分位數(shù);為平均值

    Note:is the maximum and minimum value of relative sensitivity;is the 99th and 1st percentile; The upper, middle and lower horizontal line inare the upper quartile, median and lower quartile, respectively;is the average value

    2.3.2 不同發(fā)育階段敏感性的區(qū)域分布

    圖3表明,播種-越冬期日數(shù)對(duì)溫度的相對(duì)敏感性在-0.113~0.029℃-1區(qū)間變化,平均-0.040℃-1。區(qū)域分布上,該生育階段敏感性呈顯著的中部區(qū)域高兩端低的趨勢(shì)(圖3a)。平原中南部的日數(shù)對(duì)溫度較敏感,普遍在-0.08℃-1以上。隨著緯度的增加或降低,敏感性均呈逐漸降低的趨勢(shì)。在平原南北兩端,如北京密云和河南固始等地,隨著溫度的升高,播種-越冬期生育階段日數(shù)基本不變或略有增加。華北平原約50%站點(diǎn)的敏感性通過(guò)顯著性檢驗(yàn),表明這種敏感性具有一定的普遍性。

    越冬-返青期日數(shù)對(duì)溫度的相對(duì)敏感性在-0.08~0.09℃-1區(qū)間變化,平均為0.013℃-1。該生育階段敏感性在河南波動(dòng)較大(圖3b),如河南西南部站點(diǎn)隨溫度的升高越冬-返青期日數(shù)呈明顯減小趨勢(shì),但在河南東北部則相反,即隨著溫度的升高日數(shù)呈延長(zhǎng)趨勢(shì)。平原中偏南部多數(shù)站點(diǎn)敏感性為正值,表明隨著溫度升高,日數(shù)呈延長(zhǎng)趨勢(shì),全區(qū)域的平均值也為正值,這可能是氣候變化在不同季節(jié)的非均一性所致。僅有河南鄭州、許昌和方城及山東濟(jì)寧的敏感性通過(guò)顯著性檢驗(yàn),表明這種敏感性的穩(wěn)定程度很弱。

    返青-抽穗期日數(shù)對(duì)溫度的敏感性在-0.112~-0.035℃-1區(qū)間變化,平均值為-0.074℃-1。該生育階段各區(qū)域敏感性在數(shù)值上均高于其余3個(gè)發(fā)育階段(圖3c),表明該發(fā)育階段日數(shù)最易受溫度的影響。敏感性在區(qū)域上沒(méi)有顯著的分布趨勢(shì)。除潢川外,其它所有站點(diǎn)的敏感性均通過(guò)顯著性檢驗(yàn),表明該階段的溫度敏感性非常穩(wěn)定。返青-抽穗期是冬小麥主要的營(yíng)養(yǎng)生長(zhǎng)階段,此階段敏感性較強(qiáng)的區(qū)域,如河南方城和南陽(yáng)一帶,以及山東莒縣、泰安和文登一帶等,在氣候變暖背景下,返青-抽穗期日數(shù)縮短較多,不利于干物質(zhì)的積累,對(duì)后期產(chǎn)量的形成也造成了不利影響。而河北黃驊、河間和欒城等敏感性較弱的區(qū)域,在氣候變暖背景下,返青-抽穗期日數(shù)縮短較少,因而比敏感性強(qiáng)的區(qū)域更有利于干物質(zhì)的積累。

    抽穗-成熟期日數(shù)對(duì)溫度的敏感性在-0.114~0.014℃-1區(qū)間變化,平均值為-0.042℃-1,低于返青-抽穗期(-0.074),表明在生殖生長(zhǎng)階段,溫度對(duì)發(fā)育期的作用弱于營(yíng)養(yǎng)生長(zhǎng)階段(圖3d)。其中,約有80%站點(diǎn)的敏感性通過(guò)顯著性檢驗(yàn),表明該生育階段敏感性在區(qū)域上具有較強(qiáng)的穩(wěn)定性。敏感性在區(qū)域分布上呈顯著的南高北低的趨勢(shì),表明隨著溫度的上升,南部站點(diǎn)抽穗-成熟期日數(shù)縮短比北部多,尤其是河南西南部,達(dá)到-0.06℃-1以上,而北部站點(diǎn)約為-0.02℃-1。抽穗-成熟期是產(chǎn)量形成的關(guān)鍵時(shí)期,平原南部快速降低的抽穗-成熟期日數(shù)將可能不利于冬小麥的產(chǎn)量形成,而平原北部受影響的程度較弱。

    3 結(jié)論與討論

    (1) 播種-越冬期、返青-抽穗期、抽穗-成熟期日數(shù)均隨溫度的上升而縮短,相對(duì)敏感性的區(qū)域平均值依次為-0.040℃-1、-0.074℃-1和-0.042℃-1??梢?jiàn),華北平原冬小麥返青-抽穗期對(duì)溫度最敏感,其次是抽穗-成熟,播種-越冬較為不敏感。雖然營(yíng)養(yǎng)生長(zhǎng)期的敏感性仍高于生殖生長(zhǎng)期,但二者之間的差距較前人的研究已大為縮小[15-16]。

    (2) 不同發(fā)育階段日數(shù)對(duì)溫度變化具有明顯不同的敏感性,不同區(qū)域同一發(fā)育期的敏感性也各不相同。區(qū)域分布上,播種-越冬期的敏感性呈中間高兩端低的趨勢(shì),且敏感性具有一定的普遍性;越冬-返青期的敏感性區(qū)域波動(dòng)大,但無(wú)明顯區(qū)域特征,敏感性的穩(wěn)定程度也較弱;返青-抽穗期的敏感性無(wú)明顯區(qū)域特征,敏感性的穩(wěn)定性很高;抽穗-成熟期的敏感性呈顯著的南高北低的趨勢(shì),在區(qū)域上具有較強(qiáng)的穩(wěn)定性。多數(shù)發(fā)育階段內(nèi)的敏感性區(qū)域差異較大,呈負(fù)值和正值的站點(diǎn)均有,但是敏感性為正值的站點(diǎn)無(wú)一通過(guò)顯著性檢驗(yàn)。導(dǎo)致這種現(xiàn)象的原因還有待進(jìn)一步研究。

    近30a來(lái),華北平原冬小麥各發(fā)育階段內(nèi)平均溫度均呈明顯上升趨勢(shì)。不同站點(diǎn)同一發(fā)育階段天數(shù)各不相同,同一站點(diǎn)不同發(fā)育期階段日數(shù)也存在顯著差異。因此,傳統(tǒng)的表示發(fā)育階段對(duì)日數(shù)敏感性的方法,不適于不同發(fā)育階段之間的比較,也不適于不同站點(diǎn)之間的比較。本文提出相對(duì)敏感性的概念及計(jì)算方法,可供不同發(fā)育期和不同站點(diǎn)之間進(jìn)行比較。

    與基于絕對(duì)速率的敏感性相比[15-16],本文提出的相對(duì)敏感性可較好地反映站點(diǎn)和區(qū)域尺度上冬小麥不同發(fā)育階段日數(shù)對(duì)溫度的敏感性,為分析敏感性的區(qū)域差異及隨發(fā)育期的變化特征提供了有利工具,尤其適于歷史氣候變化及品種改進(jìn)背景下溫度對(duì)發(fā)育期影響的研究。除溫度之外,土壤水分狀況和耕作制度等對(duì)發(fā)育期也有重要影響,本研究?jī)H考慮了溫度這一主要因子,尚不能完整反映其它環(huán)境因子對(duì)發(fā)育期的影響,還需開(kāi)展多因子對(duì)發(fā)育速率的綜合影響研究,以使研究結(jié)果更具應(yīng)用價(jià)值。

    參考文獻(xiàn)References

    [1]IPCC.Climate change 2007:the physical science basis[M]. Cambridge:Cambridge University Press,2007.

    [2]Garcia-Mozo H,Mestre A,Galan C.Phenological trends in southern Spain:a response to climate change[J].Agricultural and Forest Meteorology,2010,150:575-580.

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

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

    [4]王萍,李廷全,閆平,等.黑龍江省近35年氣候變化對(duì)粳稻發(fā)育期及產(chǎn)量的影響[J].中國(guó)農(nóng)業(yè)氣象,2008,29(3):268 -271.

    Wang P,Li T Q,Yan P,et al.Impacts of climate change over last 35 years on rice growing period and yield in Heilongjiang province[J].Chinese Journal of Agrometeorology,2008,29(3): 268 -271.(in Chinese)

    [5]蔣菊芳,王鶴齡,魏育國(guó),等.河西走廊東部不同類(lèi)型植物物候?qū)夂蜃兓捻憫?yīng)[J].中國(guó)農(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)

    [6]Tao F L,Zhang S,Zhang Z,et al.Maize growing duration was prolonged across China in the past three decades under the combined effects of temperature,agronomic management,and cultivar shift[J].Global Change Biology,2014,20:3686-3699.

    [7]Wang J,Wang E L,Feng L P,et al.Phenological trends of winter wheat in response to varietal and temperature changes in the North China Plain[J].Field Crops Research,2013,144: 135-144.

    [8]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[J].Global Change Biology,2013,19:563-570.

    [9]張宇,王石立,王馥棠.氣候變化對(duì)我國(guó)小麥發(fā)育及產(chǎn)量可能影響的模擬研究[J].應(yīng)用氣象學(xué)報(bào),2000,11:264-270.

    Zhang Y,Wang S L,Wang F T.Research on the possible effects of climate change on growth and yield of wheat in China[J].Quarterly Journal of Applied Meteorology,2000, 11:264-270.(in Chinese)

    [10]楊建瑩,梅旭榮,劉勤,等.氣候變化背景下華北地區(qū)冬小麥生育期的變化特征[J].植物生態(tài)學(xué)報(bào),2011,35(6):623-631.

    Yang J Y,Mei X R,Liu Q,et al.Variations of winter wheat growth stages under climate changes in northern China[J]. Chinese Journal of Plant Ecology,2011,35(6):623- 631.(in Chinese)

    [11]邱美娟,劉建棟,鄔定榮,等.華北冬小麥開(kāi)花及成熟期變化特征分析[J].干旱地區(qū)農(nóng)業(yè)研究,2014,32(3):250-255.

    Qiu M J,Liu J D,Wu D R,et al.Characteristics of variations in growth stages of winter wheat in North China Plain[J]. Agricultural Research in the Arid Areas,2014,32(3):250-255. (in Chinese)

    [12]陳書(shū)濤,王讓會(huì),許遐禎,等.氣溫及降水變化對(duì)江蘇省典型農(nóng)業(yè)區(qū)冬小麥、水稻生育期的影響[J].中國(guó)農(nóng)業(yè)氣象,2011,2(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,2(2): 235-239.(in Chinese)

    [13]余衛(wèi)東,趙國(guó)強(qiáng),陳懷亮.氣候變化對(duì)河南省主要農(nóng)作物生育期的影響[J].中國(guó)農(nóng)業(yè)氣象,2007,28(1): 9-12.

    Yu W D,Zhao G W,Chen H L.Impacts of climate change on growing stages of main crops in Henan province[J].Chinese Journal of Agrometeorology,2007,28(1):9-12.(in Chinese)

    [14]Estrella N,Sparks T H,Menzel A.Trends and temperature response in the phenology of crops in Germany[J].Global Change Biology,2007,13:1737-1747.

    [15]肖登攀,陶福祿.過(guò)去30年氣候變化對(duì)華北平原冬小麥物候的影響研究[J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào),2012,20(11): 1539-1545.

    Xiao D P,Tao F L.Impact of climate change in 1981-2009 on winter wheat phenology in the North China Plain[J].Chinese Journal of Eco-Agriculture,2012,20(11):1539-1545.(in Chinese)

    [16]肖登攀,陶福祿,沈彥俊,等.華北平原冬小麥對(duì)過(guò)去30 年氣候變化響應(yīng)的敏感性研究[J].中國(guó)生態(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)

    [17]魏鳳英.現(xiàn)代氣候統(tǒng)計(jì)診斷與預(yù)測(cè)技術(shù)[M].北京:氣象出版社,2007:37-39.

    Wei F Y.Modern climatic statistical diagnosis and prediction technology[M].Beijing:China Meteorological Press,2007:37- 39.(in Chinese)

    Relative Sensitivity of Main Growth Durations to Temperature for Winter Wheat in North China Plain

    GAO Jing1, Wu Ding-rong2,3, WANG Pei-juan2,3,CHEN Jing-hua1, YAN Feng4, ZHAO Yu-fei1, WANG Jia-qiang1

    (1.National Meteorological Information Center,Beijing 100081,China; 2.Chinese Academy of Meteorological Sciences, Beijing 100081; 3.State Key Laboratory of Severe Weather, Chinese Academy of Meteorological Sciences, Beijing 100081; 4.Meteorological Bureau of Yixian, Hebei, Baoding 074200)

    Under the condition of climate change, crop growth duration changed substantially. As a result, sensitivity of crop growth duration to temperature becomes an important part of the research of climate change impact on agriculture. Aiming to overcome shortcoming of the existing research about sensitivity, this paper put forward the concept of relative sensitivity and applied in North China Plain (NCP) to analyze the spatial-temporal characteristics of winter wheat. Using winter wheat phenology observation and daily-average temperature data from 65 agricultural meteorological observation stations in the NCP during 1980-2012, the changes of growth duration and average temperature during 4 major growth periods were calculated, including sowing to the start of overwintering (S-O), the start of overwintering to turning green (O-T), turning green to hearing (T-H) and heading to maturity (H-M). Based on the calculation, relative sensitivity of growth duration to temperature was calculated. Linear regression and GIS spatial interpolation were used to explain the results. Results showed that: (1) average temperature increased during all four growth period. Growth duration varied substantially spatially and temporally; (2) For the period S-O, sensitivity varied from -0.113 to 0.029℃-1, with an average of -0.040℃-1. Sensitivity is fairly stable during this period. Generally, middle plain is more sensitive than the north and south plain; (3) For the period O-T, sensitivity varied from -0.081 to 0.091℃-1, with an average of 0.013℃-1. Sensitivity is not stable for most stations and varied much among different stations; (4) For the period T-H, sensitivity varied from -0.112 to -0.035℃-1, with an average of -0.074℃-1. Value of sensitivity has no obvious regional distribution, but the stability of sensitivity is extremely high; (5) For the period H-M, sensitivity varied from -0.114 to 0.014℃-1, with an average of -0.042℃-1. The south plain is more sensitivity than the north. Stability of sensitivity is very high during this period. Temperature sensitivity is varied substantially among different developmental stages, and also varied significantly in different region. This study promotes understanding of regional crop response to climate change and varieties shift, and provides scientific basis for simulating phenology and hence yield response to future climate change.

    Climate change; Growth duration; Relative sensitivity to temperature; Spatial pattern; Temporal pattern

    10.3969/j.issn.1000-6362.2016.04.007

    2016-01-04

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

    國(guó)家氣象信息中心青年科技基金(NMICQJ201603);國(guó)家自然科學(xué)基金(41371410)

    高靜(1985-),女,碩士,工程師,從事氣象數(shù)據(jù)分析服務(wù)工作。E-mail:duoduolaiba@163.com

    猜你喜歡
    發(fā)育階段發(fā)育期日數(shù)
    漢江上游漢中區(qū)域不同等級(jí)降水日數(shù)的氣候變化特征分析
    綠色科技(2022年16期)2022-09-15 03:04:46
    小麥生殖發(fā)育階段對(duì)低溫的敏感性鑒定
    天津市濱海新區(qū)塘沽地域雷暴日數(shù)變化規(guī)律及特征分析
    天津科技(2020年2期)2020-03-03 05:09:48
    設(shè)施葡萄果實(shí)發(fā)育期管理技術(shù)要點(diǎn)
    對(duì)森工林區(qū)在商品林基地培養(yǎng)速生楊樹(shù)探討
    海南省雷暴日數(shù)年代際變化特征
    基于累積熱量單位的甜瓜幼苗發(fā)育期模擬
    ESSENTIAL NORMS OF PRODUCTS OF WEIGHTED COMPOSITION OPERATORS AND DIFFERENTIATION OPERATORS BETWEEN BANACH SPACES OF ANALYTIC FUNCTIONS?
    大花黃牡丹葉片發(fā)育過(guò)程中氣孔密度和氣孔指數(shù)的動(dòng)態(tài)變化
    1980—2011年阿勒泰市光熱資源變化對(duì)春小麥發(fā)育期的影響
    免费大片黄手机在线观看| 黑人高潮一二区| 国产精品嫩草影院av在线观看| av一本久久久久| 午夜福利在线在线| 亚洲第一区二区三区不卡| 成人综合一区亚洲| 波多野结衣巨乳人妻| 看黄色毛片网站| 不卡视频在线观看欧美| 欧美xxⅹ黑人| 两个人的视频大全免费| 国产精品1区2区在线观看.| 欧美区成人在线视频| 日日啪夜夜撸| 色综合亚洲欧美另类图片| 菩萨蛮人人尽说江南好唐韦庄| 色吧在线观看| 一边亲一边摸免费视频| 青青草视频在线视频观看| 18禁在线无遮挡免费观看视频| 男人和女人高潮做爰伦理| 91精品伊人久久大香线蕉| 噜噜噜噜噜久久久久久91| 日韩精品有码人妻一区| 日本-黄色视频高清免费观看| 成人午夜精彩视频在线观看| 夜夜爽夜夜爽视频| 丝瓜视频免费看黄片| 男人狂女人下面高潮的视频| 午夜福利在线观看吧| 高清日韩中文字幕在线| 日韩精品青青久久久久久| 国产 一区精品| 中文欧美无线码| 精华霜和精华液先用哪个| 国产永久视频网站| 在线天堂最新版资源| 日韩国内少妇激情av| 韩国高清视频一区二区三区| 老司机影院成人| 超碰av人人做人人爽久久| av播播在线观看一区| 全区人妻精品视频| 黄色日韩在线| 亚洲人成网站在线播| 一级黄片播放器| 少妇裸体淫交视频免费看高清| 亚洲av.av天堂| 纵有疾风起免费观看全集完整版 | 国产精品国产三级国产专区5o| 一个人观看的视频www高清免费观看| 国产高清国产精品国产三级 | 精品不卡国产一区二区三区| 国产黄频视频在线观看| 国产不卡一卡二| 精品一区二区免费观看| 777米奇影视久久| 日日摸夜夜添夜夜爱| 美女高潮的动态| 精品久久久久久久久久久久久| 欧美97在线视频| 精品久久久噜噜| 又黄又爽又刺激的免费视频.| 国产成人aa在线观看| 国产精品蜜桃在线观看| 成人亚洲精品av一区二区| 日韩一区二区三区影片| 三级毛片av免费| 亚洲三级黄色毛片| www.av在线官网国产| 久久99精品国语久久久| 69人妻影院| 国产69精品久久久久777片| 国产色婷婷99| 天堂影院成人在线观看| 非洲黑人性xxxx精品又粗又长| 亚洲精品国产av成人精品| 久久久精品欧美日韩精品| 能在线免费观看的黄片| 久久久欧美国产精品| 久久热精品热| 亚洲av成人av| 成人国产麻豆网| 真实男女啪啪啪动态图| 国产精品三级大全| 亚洲国产精品sss在线观看| 国产一级毛片七仙女欲春2| 亚洲国产av新网站| 亚洲自拍偷在线| 午夜日本视频在线| 国产免费一级a男人的天堂| 国产精品福利在线免费观看| av国产免费在线观看| 亚洲av福利一区| 夫妻午夜视频| 免费看美女性在线毛片视频| 91久久精品国产一区二区成人| 精品一区在线观看国产| 在线免费观看的www视频| 七月丁香在线播放| 久久久久免费精品人妻一区二区| 国产午夜精品一二区理论片| 黄色一级大片看看| 国内精品美女久久久久久| 日韩一本色道免费dvd| 中文字幕av成人在线电影| 亚洲成人精品中文字幕电影| 久久99蜜桃精品久久| 夫妻性生交免费视频一级片| 久久99精品国语久久久| 亚洲国产av新网站| 亚洲伊人久久精品综合| 日本猛色少妇xxxxx猛交久久| 亚洲av成人精品一二三区| 国产精品久久久久久精品电影| 久久久久久久久久久丰满| 日韩欧美精品免费久久| 水蜜桃什么品种好| 熟妇人妻不卡中文字幕| 精品99又大又爽又粗少妇毛片| 国产成人午夜福利电影在线观看| 国产精品一区二区性色av| 2022亚洲国产成人精品| 免费播放大片免费观看视频在线观看| 久久99精品国语久久久| 婷婷六月久久综合丁香| 国产成人精品婷婷| 听说在线观看完整版免费高清| 毛片一级片免费看久久久久| 99热这里只有是精品在线观看| 欧美 日韩 精品 国产| 老女人水多毛片| 国产欧美日韩精品一区二区| 极品少妇高潮喷水抽搐| 人人妻人人看人人澡| 在线a可以看的网站| 国产精品一区二区三区四区免费观看| 免费大片黄手机在线观看| 在线天堂最新版资源| 日韩av不卡免费在线播放| 日产精品乱码卡一卡2卡三| 草草在线视频免费看| 噜噜噜噜噜久久久久久91| 淫秽高清视频在线观看| 国产高清国产精品国产三级 | 国产综合懂色| 天堂av国产一区二区熟女人妻| 久久久欧美国产精品| 老司机影院毛片| 深夜a级毛片| 久久精品久久久久久噜噜老黄| 国产成年人精品一区二区| 欧美性感艳星| 男女视频在线观看网站免费| 亚洲精品国产av成人精品| 亚洲婷婷狠狠爱综合网| 欧美97在线视频| 日韩精品青青久久久久久| 亚洲av电影在线观看一区二区三区 | 高清av免费在线| 建设人人有责人人尽责人人享有的 | av国产久精品久网站免费入址| 国产日韩欧美在线精品| 精华霜和精华液先用哪个| 日本一二三区视频观看| 国产亚洲精品av在线| 看非洲黑人一级黄片| 丰满少妇做爰视频| 日韩在线高清观看一区二区三区| 亚洲综合色惰| 少妇人妻精品综合一区二区| 久久久久久久久久人人人人人人| 狂野欧美激情性xxxx在线观看| 国产精品蜜桃在线观看| 亚洲国产高清在线一区二区三| 亚洲国产欧美人成| 午夜福利成人在线免费观看| 身体一侧抽搐| 欧美一区二区亚洲| 一级爰片在线观看| 亚洲精品成人久久久久久| 成人无遮挡网站| 成人欧美大片| 日韩制服骚丝袜av| 51国产日韩欧美| 天天一区二区日本电影三级| 亚洲18禁久久av| 青春草国产在线视频| 久久久久久久久久久免费av| av在线天堂中文字幕| 别揉我奶头 嗯啊视频| ponron亚洲| 欧美激情在线99| 欧美97在线视频| 亚洲精品日韩av片在线观看| 国产美女午夜福利| 波野结衣二区三区在线| 成人特级av手机在线观看| 久久久久久久大尺度免费视频| 色网站视频免费| 国产黄a三级三级三级人| 欧美3d第一页| 我的老师免费观看完整版| 我要看日韩黄色一级片| 日日啪夜夜爽| 亚洲精品日韩在线中文字幕| av国产久精品久网站免费入址| 国产精品人妻久久久久久| 日韩制服骚丝袜av| 欧美 日韩 精品 国产| 国产精品一区二区性色av| 91在线精品国自产拍蜜月| 亚洲av成人精品一区久久| 亚洲成人精品中文字幕电影| 久久久久九九精品影院| 亚洲性久久影院| av黄色大香蕉| 十八禁网站网址无遮挡 | 成年av动漫网址| 久久久久久九九精品二区国产| 69av精品久久久久久| 国产色婷婷99| 久久久久久久久大av| 一级黄片播放器| 岛国毛片在线播放| 欧美另类一区| 国产亚洲av片在线观看秒播厂 | 18禁裸乳无遮挡免费网站照片| 99久久精品一区二区三区| 少妇被粗大猛烈的视频| 日韩亚洲欧美综合| 亚洲av二区三区四区| 国产精品日韩av在线免费观看| 免费黄网站久久成人精品| 日韩人妻高清精品专区| 国产精品熟女久久久久浪| 欧美+日韩+精品| 爱豆传媒免费全集在线观看| 久久草成人影院| 中文天堂在线官网| 2021少妇久久久久久久久久久| 国内精品一区二区在线观看| 免费大片18禁| 国语对白做爰xxxⅹ性视频网站| 免费黄网站久久成人精品| 亚洲av中文av极速乱| 国产在视频线精品| 免费大片黄手机在线观看| 国产亚洲91精品色在线| 狠狠精品人妻久久久久久综合| 91午夜精品亚洲一区二区三区| 六月丁香七月| 免费看光身美女| 一边亲一边摸免费视频| 国产日韩欧美在线精品| 国产精品一及| 亚洲精品乱码久久久v下载方式| 日韩精品青青久久久久久| 亚洲精品国产av成人精品| 亚洲性久久影院| 亚洲18禁久久av| 中文字幕亚洲精品专区| 神马国产精品三级电影在线观看| 乱系列少妇在线播放| 岛国毛片在线播放| 日本爱情动作片www.在线观看| a级一级毛片免费在线观看| 丰满人妻一区二区三区视频av| 成人二区视频| 亚洲精品乱码久久久久久按摩| 久久韩国三级中文字幕| 看黄色毛片网站| 亚洲一级一片aⅴ在线观看| 国产乱来视频区| 欧美日韩一区二区视频在线观看视频在线 | 亚洲不卡免费看| 99久国产av精品| 老女人水多毛片| 中文字幕av在线有码专区| 国产探花极品一区二区| 永久网站在线| 国产精品无大码| 国产 一区精品| 国内揄拍国产精品人妻在线| 日韩av免费高清视频| 熟女人妻精品中文字幕| 秋霞伦理黄片| 婷婷色麻豆天堂久久| 夫妻午夜视频| 人妻系列 视频| 午夜免费观看性视频| 国内揄拍国产精品人妻在线| videossex国产| 熟女人妻精品中文字幕| 亚洲欧美一区二区三区国产| 天堂俺去俺来也www色官网 | 亚洲最大成人中文| 亚洲欧美成人精品一区二区| 国产精品一区二区三区四区免费观看| 中文天堂在线官网| 精品久久久久久久久av| 嫩草影院入口| 嘟嘟电影网在线观看| 精品少妇黑人巨大在线播放| 亚洲经典国产精华液单| 亚洲精品乱码久久久久久按摩| 国产精品国产三级国产av玫瑰| 精品久久久久久久末码| 成年免费大片在线观看| 秋霞伦理黄片| 91午夜精品亚洲一区二区三区| 亚洲av电影在线观看一区二区三区 | av卡一久久| 啦啦啦韩国在线观看视频| 国产女主播在线喷水免费视频网站 | 久久久久性生活片| 久久久久免费精品人妻一区二区| 久久久久精品久久久久真实原创| 国语对白做爰xxxⅹ性视频网站| 欧美潮喷喷水| 美女脱内裤让男人舔精品视频| 精品一区在线观看国产| 亚洲国产成人一精品久久久| 精品国产一区二区三区久久久樱花 | 一级片'在线观看视频| 亚洲成人一二三区av| www.av在线官网国产| 久久久久久久久大av| 99久久精品国产国产毛片| 97超碰精品成人国产| 边亲边吃奶的免费视频| 99热网站在线观看| 日本一二三区视频观看| 搞女人的毛片| 最新中文字幕久久久久| 男插女下体视频免费在线播放| 国产一区有黄有色的免费视频 | 最新中文字幕久久久久| 高清欧美精品videossex| 免费少妇av软件| 高清午夜精品一区二区三区| 九色成人免费人妻av| 成年人午夜在线观看视频 | 99久久精品一区二区三区| 成人亚洲精品一区在线观看 | 免费观看性生交大片5| videos熟女内射| 欧美不卡视频在线免费观看| 成人毛片a级毛片在线播放| 日日撸夜夜添| 噜噜噜噜噜久久久久久91| 国产精品久久久久久精品电影小说 | 精品国产露脸久久av麻豆 | 久久精品国产自在天天线| 国产成人91sexporn| 少妇裸体淫交视频免费看高清| 欧美日韩亚洲高清精品| 亚洲av在线观看美女高潮| 黄色欧美视频在线观看| 亚洲精品久久久久久婷婷小说| 日韩精品有码人妻一区| 热99在线观看视频| 综合色av麻豆| 麻豆精品久久久久久蜜桃| av在线亚洲专区| 日本wwww免费看| 搡老妇女老女人老熟妇| 男人和女人高潮做爰伦理| 成人欧美大片| 久久热精品热| 国产 一区精品| 欧美一区二区亚洲| 成人性生交大片免费视频hd| 欧美一级a爱片免费观看看| 女人久久www免费人成看片| ponron亚洲| 观看美女的网站| 少妇熟女欧美另类| 久久鲁丝午夜福利片| 天天躁夜夜躁狠狠久久av| 日日撸夜夜添| 丰满少妇做爰视频| 联通29元200g的流量卡| av在线观看视频网站免费| 色综合亚洲欧美另类图片| 亚洲国产精品成人综合色| 久热久热在线精品观看| 免费无遮挡裸体视频| 日韩人妻高清精品专区| 亚洲欧洲国产日韩| 六月丁香七月| 在线观看人妻少妇| 久久精品夜色国产| 综合色丁香网| 成人美女网站在线观看视频| 国产精品综合久久久久久久免费| 性插视频无遮挡在线免费观看| 联通29元200g的流量卡| 久久精品夜夜夜夜夜久久蜜豆| 国产人妻一区二区三区在| 80岁老熟妇乱子伦牲交| 91久久精品电影网| 少妇丰满av| 日本av手机在线免费观看| 三级国产精品片| 中文资源天堂在线| 边亲边吃奶的免费视频| 一级毛片 在线播放| 国产老妇伦熟女老妇高清| 日韩一区二区视频免费看| 你懂的网址亚洲精品在线观看| 久久久精品欧美日韩精品| 直男gayav资源| 男女啪啪激烈高潮av片| 欧美一区二区亚洲| 国产成人a区在线观看| 亚洲激情五月婷婷啪啪| 18禁在线无遮挡免费观看视频| 大又大粗又爽又黄少妇毛片口| 亚洲欧美精品专区久久| 国产在视频线精品| 小蜜桃在线观看免费完整版高清| 亚洲成人中文字幕在线播放| 91精品伊人久久大香线蕉| 亚洲欧美成人精品一区二区| 国产 亚洲一区二区三区 | 深爱激情五月婷婷| 中文字幕制服av| 麻豆久久精品国产亚洲av| 91在线精品国自产拍蜜月| 中文在线观看免费www的网站| 国产一区二区三区av在线| 午夜福利视频1000在线观看| 国产色爽女视频免费观看| 中文天堂在线官网| 国精品久久久久久国模美| 亚洲精品一区蜜桃| 亚洲最大成人中文| 免费av毛片视频| 久热久热在线精品观看| 精品熟女少妇av免费看| 偷拍熟女少妇极品色| 亚洲精品国产av蜜桃| 伦理电影大哥的女人| 一级a做视频免费观看| 男人和女人高潮做爰伦理| 日韩欧美一区视频在线观看 | 夫妻性生交免费视频一级片| 免费av毛片视频| 2021少妇久久久久久久久久久| 成人综合一区亚洲| 久久久久久九九精品二区国产| 99热这里只有是精品在线观看| 晚上一个人看的免费电影| 国产黄色免费在线视频| 97精品久久久久久久久久精品| 成人毛片60女人毛片免费| 五月天丁香电影| 亚洲av日韩在线播放| 91精品伊人久久大香线蕉| 小蜜桃在线观看免费完整版高清| 国产亚洲91精品色在线| 又黄又爽又刺激的免费视频.| 国产不卡一卡二| 大片免费播放器 马上看| 久久精品熟女亚洲av麻豆精品 | 国产精品人妻久久久影院| 大香蕉97超碰在线| 老师上课跳d突然被开到最大视频| 一个人看的www免费观看视频| 欧美潮喷喷水| 伊人久久精品亚洲午夜| 久久久久久久亚洲中文字幕| 中国美白少妇内射xxxbb| 夫妻性生交免费视频一级片| 色吧在线观看| 亚洲最大成人av| 人妻系列 视频| 国内精品美女久久久久久| 听说在线观看完整版免费高清| 精品国内亚洲2022精品成人| 爱豆传媒免费全集在线观看| 18禁动态无遮挡网站| 国内少妇人妻偷人精品xxx网站| 欧美另类一区| 国产精品美女特级片免费视频播放器| videos熟女内射| 汤姆久久久久久久影院中文字幕 | 国产视频首页在线观看| 国产精品女同一区二区软件| 如何舔出高潮| 大话2 男鬼变身卡| 插阴视频在线观看视频| 国产午夜精品一二区理论片| 国内揄拍国产精品人妻在线| 午夜爱爱视频在线播放| 夫妻性生交免费视频一级片| 欧美成人a在线观看| 国产免费福利视频在线观看| 又爽又黄无遮挡网站| www.色视频.com| 亚洲av中文字字幕乱码综合| 久久精品久久精品一区二区三区| 永久网站在线| 免费看日本二区| 高清日韩中文字幕在线| 午夜免费男女啪啪视频观看| 毛片一级片免费看久久久久| 天天躁夜夜躁狠狠久久av| 午夜福利成人在线免费观看| 久久人人爽人人爽人人片va| 亚洲最大成人中文| 亚洲精品日本国产第一区| 国产三级在线视频| 最近最新中文字幕大全电影3| 激情 狠狠 欧美| 夫妻性生交免费视频一级片| 一级黄片播放器| 九九久久精品国产亚洲av麻豆| 日韩欧美一区视频在线观看 | 亚洲欧美精品自产自拍| 亚洲一区高清亚洲精品| 日韩三级伦理在线观看| 80岁老熟妇乱子伦牲交| 美女内射精品一级片tv| 中文字幕av成人在线电影| 熟妇人妻不卡中文字幕| 一级毛片电影观看| av天堂中文字幕网| 精品欧美国产一区二区三| 哪个播放器可以免费观看大片| 观看免费一级毛片| 国产精品精品国产色婷婷| 在线天堂最新版资源| 天堂俺去俺来也www色官网 | 99热这里只有精品一区| 夫妻午夜视频| 夫妻性生交免费视频一级片| 三级国产精品片| 99久久中文字幕三级久久日本| 小蜜桃在线观看免费完整版高清| 久久午夜福利片| 国产精品国产三级专区第一集| 午夜久久久久精精品| 国产精品精品国产色婷婷| 夜夜看夜夜爽夜夜摸| 亚洲成人一二三区av| 亚洲av福利一区| 免费大片黄手机在线观看| 国产人妻一区二区三区在| 成人国产麻豆网| 波多野结衣巨乳人妻| videossex国产| 亚洲真实伦在线观看| 免费黄频网站在线观看国产| 天美传媒精品一区二区| 国产高潮美女av| 中文字幕亚洲精品专区| 欧美另类一区| 精品国产露脸久久av麻豆 | 天天躁夜夜躁狠狠久久av| 别揉我奶头 嗯啊视频| 成人特级av手机在线观看| 大陆偷拍与自拍| 国产69精品久久久久777片| 亚洲国产欧美在线一区| 国产精品一区二区三区四区免费观看| 最近2019中文字幕mv第一页| 色综合亚洲欧美另类图片| 色5月婷婷丁香| 国产精品嫩草影院av在线观看| 日韩伦理黄色片| 久久久久久久久中文| 亚洲综合精品二区| 欧美日韩在线观看h| 一级爰片在线观看| av在线播放精品| 你懂的网址亚洲精品在线观看| 精品人妻熟女av久视频| 国产午夜福利久久久久久| 波野结衣二区三区在线| 日韩亚洲欧美综合| 日日啪夜夜撸| 亚洲精品成人av观看孕妇| 亚洲熟女精品中文字幕| 高清在线视频一区二区三区| 久久人人爽人人爽人人片va| av在线亚洲专区| 嫩草影院入口| 日本黄大片高清| 99热全是精品| 久久久久精品久久久久真实原创| videos熟女内射| 国产一区二区亚洲精品在线观看| 日韩成人伦理影院| 精品久久久久久久末码| 精品熟女少妇av免费看| 午夜激情福利司机影院| 中文欧美无线码| 欧美成人午夜免费资源| 国产成人免费观看mmmm| 亚洲三级黄色毛片| 亚洲av电影在线观看一区二区三区 | 成年人午夜在线观看视频 | 69人妻影院| 国内精品美女久久久久久| 国产一级毛片在线| 免费电影在线观看免费观看| 日韩欧美一区视频在线观看 | 亚洲激情五月婷婷啪啪| av在线老鸭窝| 久久韩国三级中文字幕| 街头女战士在线观看网站| 国产高清三级在线|