• <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| 午夜两性在线视频| 80岁老熟妇乱子伦牲交| cao死你这个sao货| 精品无人区乱码1区二区| www.精华液| 午夜福利在线免费观看网站| 深夜精品福利| 日本黄色日本黄色录像| 免费高清在线观看日韩| 天天添夜夜摸| 亚洲中文av在线| 久久香蕉国产精品| 亚洲九九香蕉| 电影成人av| 亚洲熟女精品中文字幕| 最新在线观看一区二区三区| 夫妻午夜视频| av国产精品久久久久影院| 91九色精品人成在线观看| 又黄又爽又免费观看的视频| 建设人人有责人人尽责人人享有的| 乱人伦中国视频| 一区福利在线观看| 男女床上黄色一级片免费看| 99国产综合亚洲精品| 亚洲欧美激情综合另类| 在线十欧美十亚洲十日本专区| 国产亚洲精品第一综合不卡| 亚洲av日韩在线播放| 首页视频小说图片口味搜索| 激情在线观看视频在线高清 | av超薄肉色丝袜交足视频| 亚洲综合色网址| 国产免费男女视频| 天天影视国产精品| 老鸭窝网址在线观看| 欧美中文综合在线视频| 中文亚洲av片在线观看爽 | 窝窝影院91人妻| 欧美亚洲日本最大视频资源| 美女高潮到喷水免费观看| 69精品国产乱码久久久| 日日摸夜夜添夜夜添小说| 黄色 视频免费看| 一本一本久久a久久精品综合妖精| 最新的欧美精品一区二区| 校园春色视频在线观看| 午夜激情av网站| 亚洲国产看品久久| 最新的欧美精品一区二区| 欧美大码av| 国产精品亚洲一级av第二区| 亚洲精品美女久久av网站| 丝袜人妻中文字幕| 亚洲欧美激情综合另类| 一级毛片精品| 国产亚洲精品第一综合不卡| 国产精品 国内视频| 美女扒开内裤让男人捅视频| 丰满饥渴人妻一区二区三| 亚洲成人国产一区在线观看| 日韩有码中文字幕| 黄片播放在线免费| 下体分泌物呈黄色| 黑人操中国人逼视频| 一区在线观看完整版| 欧美丝袜亚洲另类 | 大型av网站在线播放| 亚洲黑人精品在线| 亚洲一卡2卡3卡4卡5卡精品中文| 大片电影免费在线观看免费| 一级毛片高清免费大全| 九色亚洲精品在线播放| 欧美+亚洲+日韩+国产| 国产免费av片在线观看野外av| 天天添夜夜摸| 麻豆国产av国片精品| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲在线自拍视频| 欧美日韩瑟瑟在线播放| 韩国精品一区二区三区| 国产一区有黄有色的免费视频| 精品国内亚洲2022精品成人 | 又黄又粗又硬又大视频| 老司机在亚洲福利影院| 好男人电影高清在线观看| 80岁老熟妇乱子伦牲交| 亚洲性夜色夜夜综合| 国产精华一区二区三区| 色综合欧美亚洲国产小说| 欧美最黄视频在线播放免费 | 日日爽夜夜爽网站| 97人妻天天添夜夜摸| 乱人伦中国视频| 天堂动漫精品| 中亚洲国语对白在线视频| 国产成人免费观看mmmm| 精品国内亚洲2022精品成人 | 亚洲精品美女久久久久99蜜臀| cao死你这个sao货| 十分钟在线观看高清视频www| 成人18禁高潮啪啪吃奶动态图| 成年版毛片免费区| 婷婷精品国产亚洲av在线 | 欧美成人免费av一区二区三区 | 久久久久久久国产电影| 一区二区日韩欧美中文字幕| 亚洲人成77777在线视频| 99riav亚洲国产免费| 人妻久久中文字幕网| 色综合婷婷激情| 一区二区三区激情视频| 国产xxxxx性猛交| 黑丝袜美女国产一区| 久久九九热精品免费| 国产精品欧美亚洲77777| 97人妻天天添夜夜摸| 欧美黄色片欧美黄色片| bbb黄色大片| 国产成人一区二区三区免费视频网站| 国产精品乱码一区二三区的特点 | 国产成人系列免费观看| 国产精品成人在线| 香蕉久久夜色| 欧美黑人欧美精品刺激| 女人被狂操c到高潮| 欧美激情久久久久久爽电影 | 精品一区二区三卡| 国产成人影院久久av| 黄色怎么调成土黄色| 黑人巨大精品欧美一区二区蜜桃| 久久精品91无色码中文字幕| 精品高清国产在线一区| 日韩欧美一区视频在线观看| 欧美乱妇无乱码| 亚洲成人手机| 在线永久观看黄色视频| 操出白浆在线播放| 亚洲精品久久午夜乱码| 国产成人欧美在线观看 | 精品卡一卡二卡四卡免费| 看黄色毛片网站| 少妇 在线观看| 亚洲国产看品久久| 俄罗斯特黄特色一大片| 99久久国产精品久久久| 国产精品一区二区免费欧美| 国产精品99久久99久久久不卡| 村上凉子中文字幕在线| 人人妻人人澡人人爽人人夜夜| 亚洲熟女精品中文字幕| av天堂久久9| 女性被躁到高潮视频| 色在线成人网| 黄色视频,在线免费观看| 一级片'在线观看视频| 亚洲中文字幕日韩| 少妇裸体淫交视频免费看高清 | 黑人操中国人逼视频| 久久精品熟女亚洲av麻豆精品| 国产精华一区二区三区| 亚洲一区中文字幕在线| 女性生殖器流出的白浆| 免费av中文字幕在线| 国产精品综合久久久久久久免费 | 啦啦啦在线免费观看视频4| 国产男靠女视频免费网站| 一级片免费观看大全| 我的亚洲天堂| 国产又色又爽无遮挡免费看| av福利片在线| 一级黄色大片毛片| 久久国产精品大桥未久av| svipshipincom国产片| 每晚都被弄得嗷嗷叫到高潮| 中文字幕人妻丝袜制服| 亚洲精品国产区一区二| 国产成人精品久久二区二区91| 91老司机精品| 久久精品熟女亚洲av麻豆精品| 后天国语完整版免费观看| 国产成人av激情在线播放| 欧美日韩亚洲高清精品| av网站免费在线观看视频| 黄色怎么调成土黄色| 亚洲国产精品sss在线观看 | 亚洲欧美激情综合另类| 建设人人有责人人尽责人人享有的| 国产一区二区三区在线臀色熟女 | 亚洲成人免费电影在线观看| 亚洲精品美女久久av网站| 国产精品成人在线| 国产乱人伦免费视频| 日韩欧美三级三区| 亚洲综合色网址| 亚洲专区字幕在线| 一区福利在线观看| 精品无人区乱码1区二区| 国产精品电影一区二区三区 | 欧美日韩亚洲国产一区二区在线观看 | 午夜视频精品福利| 亚洲国产欧美一区二区综合| 成人永久免费在线观看视频| 在线天堂中文资源库| 大码成人一级视频| 久久久精品区二区三区| 久久精品国产亚洲av高清一级| 午夜免费鲁丝| 欧美黑人精品巨大| e午夜精品久久久久久久| 国产亚洲av高清不卡| 久久国产精品大桥未久av| 成人影院久久| 乱人伦中国视频| 18在线观看网站| 黑人操中国人逼视频| 一进一出好大好爽视频| 99久久综合精品五月天人人| 日本五十路高清| 视频在线观看一区二区三区| 热99国产精品久久久久久7| 精品人妻熟女毛片av久久网站| 久久久久视频综合| 国产无遮挡羞羞视频在线观看| 91av网站免费观看| 美女福利国产在线| 国产精品二区激情视频| 欧美中文综合在线视频| 亚洲av成人不卡在线观看播放网| 18禁裸乳无遮挡免费网站照片 | 视频区欧美日本亚洲| 免费久久久久久久精品成人欧美视频| 国产亚洲精品第一综合不卡| 18禁黄网站禁片午夜丰满| 97人妻天天添夜夜摸| 亚洲av电影在线进入| 亚洲欧洲精品一区二区精品久久久| 亚洲伊人色综图| 在线观看免费日韩欧美大片| 在线观看66精品国产| 国产男女超爽视频在线观看| 精品少妇久久久久久888优播| 国产淫语在线视频| 十分钟在线观看高清视频www| 国产一区二区三区视频了| 亚洲一码二码三码区别大吗| 亚洲欧美一区二区三区久久| 国产三级黄色录像| 我的亚洲天堂| 亚洲第一青青草原| 色尼玛亚洲综合影院| 婷婷精品国产亚洲av在线 | 91成人精品电影| 两性夫妻黄色片| 成年人午夜在线观看视频| 国产精品九九99| av一本久久久久| videos熟女内射| 欧美色视频一区免费| 亚洲色图av天堂| 国产av精品麻豆| 高潮久久久久久久久久久不卡| 亚洲精品久久成人aⅴ小说| 日本撒尿小便嘘嘘汇集6| 亚洲精品在线美女| 高清毛片免费观看视频网站 | 婷婷成人精品国产| 在线观看免费日韩欧美大片| 亚洲五月色婷婷综合| 日韩欧美免费精品| 亚洲午夜理论影院| 国产三级黄色录像| 国精品久久久久久国模美| 国产免费av片在线观看野外av| 亚洲精品乱久久久久久| 涩涩av久久男人的天堂| 日韩一卡2卡3卡4卡2021年| 国产99白浆流出| 久久久久久久久久久久大奶| www日本在线高清视频| 国产精品一区二区在线观看99| 黄色丝袜av网址大全| 少妇粗大呻吟视频| 欧美激情高清一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 成人av一区二区三区在线看| 黄色视频,在线免费观看| 国产精品亚洲av一区麻豆| 亚洲午夜精品一区,二区,三区| 国产不卡av网站在线观看| 久久性视频一级片| 国产成人精品久久二区二区免费| 国产成人系列免费观看| 成人黄色视频免费在线看| 欧美日韩国产mv在线观看视频| 日韩 欧美 亚洲 中文字幕| 高清毛片免费观看视频网站 | √禁漫天堂资源中文www| 超碰97精品在线观看| av有码第一页| 精品久久蜜臀av无| 一二三四在线观看免费中文在| 亚洲男人天堂网一区| 首页视频小说图片口味搜索| 精品亚洲成a人片在线观看| 国产精品美女特级片免费视频播放器 | 最近最新中文字幕大全电影3 | 国产成人欧美在线观看 | 国产精品欧美亚洲77777| 天天影视国产精品| 色综合欧美亚洲国产小说| 亚洲美女黄片视频| 成人国语在线视频| 色在线成人网| av片东京热男人的天堂| 日本五十路高清| 无人区码免费观看不卡| 91字幕亚洲| 波多野结衣一区麻豆| 久久香蕉精品热| 最近最新免费中文字幕在线| 水蜜桃什么品种好| 精品免费久久久久久久清纯 | 女人高潮潮喷娇喘18禁视频| x7x7x7水蜜桃| av中文乱码字幕在线| 99精品久久久久人妻精品| 精品国产亚洲在线| 中文欧美无线码| 亚洲精品一卡2卡三卡4卡5卡| 99国产精品一区二区蜜桃av | 亚洲人成电影观看| 国产亚洲一区二区精品| 青草久久国产| 亚洲国产欧美日韩在线播放| 日韩有码中文字幕| 亚洲av成人不卡在线观看播放网| 日本vs欧美在线观看视频| 天天躁夜夜躁狠狠躁躁| 免费在线观看亚洲国产| av中文乱码字幕在线| 精品久久久久久电影网| videos熟女内射| aaaaa片日本免费| 久久久久久久国产电影| 亚洲人成伊人成综合网2020| 亚洲av成人不卡在线观看播放网| av片东京热男人的天堂| 午夜影院日韩av| 最新在线观看一区二区三区| 日本欧美视频一区| 久久久国产一区二区| 日本欧美视频一区| 法律面前人人平等表现在哪些方面| 精品视频人人做人人爽| av一本久久久久| 极品少妇高潮喷水抽搐| 国产欧美亚洲国产| 亚洲精品中文字幕一二三四区| 久久人妻av系列| 搡老岳熟女国产| 高清毛片免费观看视频网站 | 美女 人体艺术 gogo| 亚洲一区中文字幕在线| 99在线人妻在线中文字幕 | 精品亚洲成国产av| 这个男人来自地球电影免费观看| 中文字幕av电影在线播放| 久久精品熟女亚洲av麻豆精品| 欧美日韩视频精品一区| 无限看片的www在线观看| 国产一卡二卡三卡精品| 在线观看日韩欧美| av不卡在线播放| 老汉色av国产亚洲站长工具| 国产成人免费观看mmmm| 一区二区三区激情视频| 中文字幕人妻丝袜制服| 精品少妇久久久久久888优播| 欧美在线黄色| 久9热在线精品视频| 亚洲国产精品合色在线| av有码第一页| 亚洲精品国产色婷婷电影| av天堂在线播放| 最新的欧美精品一区二区| 好看av亚洲va欧美ⅴa在| 中文字幕人妻丝袜一区二区| 亚洲综合色网址| 嫩草影视91久久| 在线观看午夜福利视频| 亚洲成人免费电影在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 乱人伦中国视频| 亚洲一区高清亚洲精品| 亚洲精品久久成人aⅴ小说| 色尼玛亚洲综合影院| 1024香蕉在线观看| 99精品在免费线老司机午夜| 极品教师在线免费播放| 岛国在线观看网站| 又大又爽又粗| 每晚都被弄得嗷嗷叫到高潮| 色婷婷av一区二区三区视频| 中文字幕av电影在线播放| 国产亚洲欧美在线一区二区| 成人永久免费在线观看视频| 夜夜爽天天搞| 激情视频va一区二区三区| 波多野结衣一区麻豆| 久久人人97超碰香蕉20202| 美女 人体艺术 gogo| 村上凉子中文字幕在线| √禁漫天堂资源中文www| 亚洲精品一二三| 我的亚洲天堂| 国产精品自产拍在线观看55亚洲 | 最近最新免费中文字幕在线| 国产一区二区三区视频了| 日本一区二区免费在线视频| 欧美av亚洲av综合av国产av| 国产精品久久电影中文字幕 | 大片电影免费在线观看免费| 亚洲成人国产一区在线观看| 国产成人免费无遮挡视频| 露出奶头的视频| 一本综合久久免费| a级毛片在线看网站| 日本黄色视频三级网站网址 | 99re在线观看精品视频| 成人18禁高潮啪啪吃奶动态图| 捣出白浆h1v1| 韩国av一区二区三区四区| 亚洲七黄色美女视频| 一二三四在线观看免费中文在| 国产精品成人在线| av线在线观看网站| 高清av免费在线| 在线永久观看黄色视频| 亚洲精品av麻豆狂野| 五月开心婷婷网| 国产在视频线精品| 久久久久久久久久久久大奶| 美女视频免费永久观看网站| 性少妇av在线| 精品视频人人做人人爽| 夜夜躁狠狠躁天天躁| 视频区欧美日本亚洲| 99国产精品一区二区三区| 中亚洲国语对白在线视频| 超碰97精品在线观看| 精品无人区乱码1区二区| 日本欧美视频一区| 欧美久久黑人一区二区| 亚洲三区欧美一区| 欧美日韩成人在线一区二区| 女人被躁到高潮嗷嗷叫费观| 夜夜夜夜夜久久久久| 国产不卡一卡二| 国产亚洲欧美精品永久| 国产精品免费一区二区三区在线 | 国产精品久久久人人做人人爽| 成人影院久久| 国产一区二区三区综合在线观看| 香蕉久久夜色| 一边摸一边抽搐一进一小说 | 久久久久久亚洲精品国产蜜桃av| 午夜福利乱码中文字幕| 亚洲精品久久午夜乱码| 国产有黄有色有爽视频| 国产亚洲av高清不卡| 麻豆国产av国片精品| 国产男靠女视频免费网站| 国产免费男女视频| 欧美在线一区亚洲| 亚洲精品国产区一区二| 欧美丝袜亚洲另类 | 天天影视国产精品| 电影成人av| 天堂俺去俺来也www色官网| 国产成人一区二区三区免费视频网站| 国产淫语在线视频| 久久青草综合色| 狠狠婷婷综合久久久久久88av| 麻豆成人av在线观看| 欧美日韩一级在线毛片| 国产xxxxx性猛交| 9色porny在线观看| 久久中文字幕一级| 一级黄色大片毛片| 日本wwww免费看| 日韩 欧美 亚洲 中文字幕| 人人澡人人妻人| 叶爱在线成人免费视频播放| 麻豆乱淫一区二区| 色婷婷av一区二区三区视频| 欧美午夜高清在线| 国产成人精品久久二区二区免费| 母亲3免费完整高清在线观看| 999久久久精品免费观看国产| 波多野结衣一区麻豆| 久久久国产欧美日韩av| xxx96com| 91av网站免费观看| 在线观看www视频免费| 色婷婷av一区二区三区视频| 国产不卡av网站在线观看| 国产成人精品久久二区二区免费| 日韩人妻精品一区2区三区| 欧美 日韩 精品 国产| 99久久综合精品五月天人人| 亚洲熟妇中文字幕五十中出 | 亚洲伊人色综图| 免费av中文字幕在线| 日韩熟女老妇一区二区性免费视频| 亚洲一区中文字幕在线| 在线观看日韩欧美| 久久精品熟女亚洲av麻豆精品| 99精品在免费线老司机午夜| 99国产精品一区二区蜜桃av | 精品欧美一区二区三区在线| 亚洲久久久国产精品| 超碰成人久久| 亚洲熟女精品中文字幕| 国产真人三级小视频在线观看| 国产成人av激情在线播放| 免费女性裸体啪啪无遮挡网站| 丰满饥渴人妻一区二区三| 欧美黄色淫秽网站| 欧美黄色片欧美黄色片| 老司机深夜福利视频在线观看| 久久九九热精品免费| 国产一区有黄有色的免费视频| 精品一区二区三区av网在线观看| 日本a在线网址| 黄色 视频免费看| 老司机福利观看| 精品人妻熟女毛片av久久网站| 欧美国产精品va在线观看不卡| 成年人黄色毛片网站| 欧美av亚洲av综合av国产av| 亚洲精品久久午夜乱码| 99国产精品一区二区蜜桃av | 欧美日韩亚洲国产一区二区在线观看 | 黄片大片在线免费观看| 女人高潮潮喷娇喘18禁视频| 黑人欧美特级aaaaaa片| 男女免费视频国产| 999久久久国产精品视频| 色尼玛亚洲综合影院| 日韩视频一区二区在线观看| 国产亚洲av高清不卡| 欧美丝袜亚洲另类 | 91麻豆精品激情在线观看国产 | 纯流量卡能插随身wifi吗| 免费在线观看亚洲国产| 国产成人系列免费观看| 99国产精品免费福利视频| 国产精品一区二区在线不卡| 在线永久观看黄色视频| 国产精品一区二区免费欧美| av有码第一页| 脱女人内裤的视频| 国产亚洲欧美在线一区二区| 亚洲欧美一区二区三区久久| 少妇裸体淫交视频免费看高清 | 精品国产乱码久久久久久男人| 757午夜福利合集在线观看| 视频在线观看一区二区三区| 老汉色∧v一级毛片| 三级毛片av免费| 搡老岳熟女国产| 国产片内射在线| 国产欧美日韩精品亚洲av| 在线国产一区二区在线| 男女之事视频高清在线观看| 一级a爱视频在线免费观看| tocl精华| 丝瓜视频免费看黄片| 亚洲第一欧美日韩一区二区三区| 人妻 亚洲 视频| 久久人人爽av亚洲精品天堂| 久久久国产精品麻豆| 国产精品乱码一区二三区的特点 | 精品视频人人做人人爽| 在线国产一区二区在线| 国产精品成人在线| 国产淫语在线视频| 9191精品国产免费久久| 欧美日韩视频精品一区| 久久久精品免费免费高清| 午夜福利影视在线免费观看| 亚洲熟女精品中文字幕| 久久精品aⅴ一区二区三区四区| 欧美乱码精品一区二区三区| 欧美色视频一区免费| 少妇被粗大的猛进出69影院| 国产亚洲精品一区二区www | 亚洲欧美色中文字幕在线| 国产1区2区3区精品| 在线观看免费日韩欧美大片| 精品视频人人做人人爽| 咕卡用的链子| 在线观看日韩欧美| 国产91精品成人一区二区三区| 757午夜福利合集在线观看| 精品久久久久久电影网| 国产欧美日韩精品亚洲av| 侵犯人妻中文字幕一二三四区|