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

    基于SPI指數(shù)的農(nóng)作物生長期干旱時間變化研究

    2014-09-10 10:37:13申海鳳商彥蕊劉公英
    湖北農(nóng)業(yè)科學(xué) 2014年11期
    關(guān)鍵詞:邢臺縣干旱夏玉米

    申海鳳+商彥蕊+劉公英

    摘要:邢臺縣作為河北省重要糧食生產(chǎn)縣之一,受旱災(zāi)影響日趨嚴重。根據(jù)1954-2010年逐月降水?dāng)?shù)據(jù),采用標準化降水指數(shù)(SPI)進行趨勢分析和M-K趨勢突變檢驗,研究邢臺縣冬小麥和夏玉米不同生育階段干旱發(fā)生的規(guī)律。結(jié)果表明,邢臺縣干旱發(fā)生強度和頻率都在增大,以中旱和重旱為主,20世紀50年代后期干旱趨勢顯著,60年代和90年代干旱趨勢相對明顯。冬小麥SPI1指數(shù)下降速率在灌漿-成熟階段最大、拔節(jié)階段次之、再之是苗生長階段,此3個階段干旱加劇較其它生育階段明顯,除播種前和抽穗-開花階段,其他各生育階段SPI指數(shù)自1980年初均有不同程度下降,干旱趨勢相對明顯。夏玉米生育期的早中期階段SPI1指數(shù)呈下降的趨勢,播種-拔節(jié)階段、拔節(jié)-抽雄階段是干旱的頻發(fā)期,除播種前和灌漿-成熟階段,其他各生育階段SPI1指數(shù)自1970年后期開始下降,呈現(xiàn)干旱趨勢。

    關(guān)鍵詞:邢臺縣;冬小麥;夏玉米;生育階段;干旱

    中圖分類號:S423文獻標識碼:A文章編號:0439-8114(2014)11-2536-06

    The Time Changing of Drought in Crop Growth Period Based on the SPI Index

    ——Taking Xingtai County in Hebei Province as an Example

    SHEN Hai-feng1,SHANG Yan-rui1,2, LIU Gong-ying1

    (1.School of Resources and Environmental Science,Hebei Normal University,Shijiazhuang 050024,China;

    2.Hebei Key Laboratory of Environmental Change and Ecological Construction, Shijiazhuang 050016, China)

    Abstract: Xingtai county, as one of the major grain production counties in Hebei province, is seriously affected by drought. Based on the monthly precipitation data of 1954-2010 years, standardized precipitation index (SPI) for trend fitting and Mann-Kendall trend mutation testing was used to study the inter-annual variation of meteorological drought of the winter wheat and summer maize growing in Xingtai county. More attention was payed to analyze the regularity of drought in its different growth stages. The results showed that intensity and frequency of droughts was increasing in Xingtai county. Moderate drought and severe drought was the main trend. The trend of drought was significant in the late 1950s, relatively more obvious in late 1960s and 1990s. Winter wheat SPI1 index was decreased fast in filling to maturing stage of winter wheat, then came to the jointing stage. The last was the further growth of the seedling stage. The three stages was meteorological drought intensified significantly over other growth stages before sowing and heading-flowering stage. The SPI1 Index of other growth stages was declined in varying degrees since the early 1980s, but the trend was relatively obvious. Summer corn SPI1 index showed a downward trend in summer maize growth period. Drought was frequent in planting-jointing stage and jointing-heading stage, in addition to before planting and filling-mature stage. Other growth stages of the SPI1 index started to decline since the late 1970s, with the same trend as that of the drought.

    Key words: Xingtai country; winterwheat; summermaize; growth period; drought

    基金項目:國家自然科學(xué)基金項目(NO41171402)

    干旱是一種嚴重威脅人類生存的自然現(xiàn)象。干旱出現(xiàn)的頻率高、持續(xù)時間長、波及范圍廣,對國民經(jīng)濟特別是對農(nóng)業(yè)生產(chǎn)有重大的影響。為了監(jiān)測和評價干旱的發(fā)生強度及趨勢,國內(nèi)外的學(xué)者已提出了多種干旱指數(shù),包括降水距平指數(shù)、降水距平百分率、降水異常指數(shù)(RAI)、表層水供應(yīng)指數(shù)(SWSI)、Palmer干旱指數(shù)(PDSI)、標準化降水指數(shù)(SPI)、Z指數(shù)、土壤濕度指數(shù)等[1],其中,基于降水的標準化降水指數(shù)(SPI)具有數(shù)據(jù)獲取、計算簡便的優(yōu)點,能夠較好地反映出不同時間尺度下干旱發(fā)生特征,而且在不同地區(qū)和不同時段均具有良好的穩(wěn)定性和連續(xù)性,因而被廣泛應(yīng)用[2],如車少靜等[3]利用SPI指數(shù)分析了近41年河北省旱澇時空變化分布特征;郁凌華等[4]利用SPI指數(shù)分析了黃淮海地區(qū)夏玉米生長季內(nèi)的旱澇災(zāi)害特征;陳莉等[5]利用SPI指數(shù)對東北地區(qū)近50年農(nóng)作物生長季干旱趨勢進行了研究。應(yīng)用邢臺縣1954-2010年逐月降水資料,對邢臺縣SPI指數(shù)進行了測算和分析[6],研究對冬小麥-夏玉米不同生育階段干旱發(fā)生規(guī)律,為政府和農(nóng)戶進行科學(xué)的灌溉提供依據(jù)。

    邢臺縣位于河北省西南部,太行山南端東麓,屬暖溫帶半干旱季風(fēng)氣候區(qū),大陸性季風(fēng)特征顯著[7]。該縣年平均氣溫為12~13 ℃,年平均降水量為534.86 mm,年平均蒸發(fā)量為1 093.4 mm,蒸發(fā)量遠大于降水量,且降水季節(jié)分布不勻,夏季平均降水量341.5 mm,冬季平均降水量僅為14.7 mm。邢臺縣農(nóng)業(yè)以種植糧食作物為主,經(jīng)濟作物為輔,其中,糧食作物以種植冬小麥和夏玉米為主。據(jù)河北統(tǒng)計年鑒數(shù)據(jù)統(tǒng)計知,2010年邢臺縣農(nóng)作物播種面積達40 099 hm2,糧食播種面積達到30 939 hm2,占農(nóng)作物總播種面積的77.2%,其中冬小麥播種面積為12 133 hm2,夏玉米播種面積為15 322 hm2,分別占糧食總播種面積的39.2%和49.5%。邢臺縣農(nóng)業(yè)災(zāi)害頻發(fā),2010年受災(zāi)面積達16 581 hm2,成災(zāi)面積達14 764 hm2,其中,旱災(zāi)受災(zāi)面積11 845 hm2,旱災(zāi)成災(zāi)面積11 271 hm2。在成災(zāi)面積中糧食作物減產(chǎn)三至五成面積達8 596 hm2,減產(chǎn)五至八成面積達826 hm2,減產(chǎn)8成以上達460 hm2,旱情十分嚴峻。

    1材料與方法

    1.1SPI指數(shù)

    SPI是由McKee等[8]在1993年研究不同時間尺度干旱的發(fā)生頻率及持續(xù)時間之間的關(guān)系時首次提出。一般情況下,降水呈現(xiàn)偏態(tài)分布,在統(tǒng)計計算中常用Gamma密度函數(shù)來擬合不同時間尺度的降水分布規(guī)律[9]。SPI是基于降水累積概率分布的一個有效分析干濕周期的指標,它將降水的偏態(tài)分布轉(zhuǎn)換為標準的正態(tài)分布,進行統(tǒng)計計算。

    假設(shè)某一時段的降水量為x,則其伽馬分布的概率密度函數(shù)為:

    g(x)=xa-1e-x/?茁 (1)

    Γ(α)=yα-1e-ydy (2)

    式中:Γ(α)為伽馬函數(shù);α為形狀參數(shù),α> 0;β為尺度參數(shù),β>0;x為降水量,x>0。

    α,β可以采用最大似然估計法估算:

    = (3)

    A=ln()-ln(x) (4)

    式中n為計算系列的長度。于是給定時間尺度的累積概率為:

    G(X)=g(x)dx=xα-1e-x/?茁dx (5)

    令t=x/β,上式可變?yōu)椴煌耆馁ゑR函數(shù),

    G(x)=xe-x-?茁dx(6)

    由于伽馬方程不包含x=0的情況,而實際的降水量x可以為0,所以累計概率可以表示為:

    H(x)=q+(1-q)G(x) (7)

    式中,q是降水為0的概率。如果m表示時間系列中降水量為零值的個數(shù),則q=m/n。累積概率H(x)可以通過下式轉(zhuǎn)換為標準正態(tài)分布函數(shù):

    Z=SPI=-t- (8)

    其中,t= 0≤H(x)≤0.5 0.5≤H(x)<1(9)

    c0=2.515 517,c1=0.802 853,c2=0.010 328,

    d1=1.432 788,d2=0.189 269,d3=0.001 308。

    在干旱發(fā)生的時間特征上,計算1個月和12個月尺度的SPI值。1個月時間尺度的SPI1值可以表征短時間內(nèi)降水條件,它直接影響土壤的短期水分條件和農(nóng)作物不同生長階段水分的供給條件[7]。12個月尺度SPI12值則反應(yīng)了一個區(qū)域長期的降水條件,從側(cè)面揭示一個區(qū)域的地下水補給情況、水庫庫存水量的狀況,進而分析一個區(qū)域的農(nóng)業(yè)發(fā)展所需的水量條件。

    SPI值隨著時間序列的變化在0值附近頻繁的上下波動。SPI值大于0,表示降水量超出了同一時期平均降水量,則可能發(fā)生洪澇;SPI值小于0,表示降水量低于同時期平均降水量,則可能發(fā)生干旱。SPI的絕對值越大表示發(fā)生旱澇的可能越大。參照Mckee等[10]干旱等級劃分干旱等級:1.0>SPI>-1.0,表示降水正常;-1.00≥SPI>-1.5,表示發(fā)生中等干旱;-1.5≥SPI>-2.0,表示發(fā)生重度干旱;SPI≤-2.0,則表示發(fā)生極端干旱。

    1.2氣象統(tǒng)計方法

    使用線性回歸方程y=at+b(y為SPI值,a為斜率,t為降水時間序列)來擬合SPI指數(shù)隨時間序列的變化,可以反映出SPI指數(shù)隨時間的線性變化規(guī)律[11-13]。

    運用M-K非參數(shù)檢驗法[11-13]對冬小麥和夏玉米不同生育階段SPI指數(shù)做趨勢及突變檢驗:UFk值>0,表示原序列SPI指數(shù)呈現(xiàn)上升趨勢,反之則呈現(xiàn)下降趨勢;若UFk值超過置信度線(u0.05=±1.96)時,表示原序列SPI指數(shù)呈上升或下降趨勢明顯;若UFk和UBk兩條曲線的交叉點位于置信區(qū)間內(nèi),則該點對應(yīng)的時間便是序列突變開始的時間。

    2干旱時間變化規(guī)律

    2.1邢臺縣冬小麥生長期干旱年際變化

    圖1 SPI12顯示了1954-2011年邢臺縣冬小麥-夏玉米干旱發(fā)生的年際變化。數(shù)據(jù)顯示,58年間邢臺縣共發(fā)生中旱程度以上旱災(zāi)92次,其中中旱56次、重旱26次、極端干旱10次。20世紀80年代以來,干旱事件發(fā)生頻次明顯增加,其中1986、1987、1999和2000年出現(xiàn)了極端干旱。整體擬合趨勢表明,隨時間向后推移SPI指數(shù)整體呈現(xiàn)遞減趨勢,說明邢臺縣干旱整體呈現(xiàn)加劇趨勢。

    2.2冬小麥不同生育階段干旱變化

    冬小麥生育期為10月初至次年6月上旬,不同發(fā)育階段起止時間見表1,冬小麥SPI變化規(guī)律見圖2。干旱對于冬小麥不同生育階段影響不同。播種前,發(fā)生干旱,土壤底墑不足,影響冬小麥出苗;苗生長和越冬階段,冬小麥處于營養(yǎng)生長階段,小麥耗水較少,受干旱影響小,但苗生長初期發(fā)生干旱將影響幼苗生長,極端干旱可能導(dǎo)致幼苗枯死,影響后期作物產(chǎn)量;返青階段,冬小麥生長迅速,小麥對水分需求速度和需求量明顯增大,干旱發(fā)生影響冬小麥返青拔節(jié);拔節(jié)階段,干旱影響其生殖生長,造成冬小麥減產(chǎn);抽穗-開花階段,冬小麥葉面積達到最大,水分虧缺會造成冬小麥減產(chǎn)愈加明顯;灌漿成熟階段,干旱影響小麥粒數(shù)和粒徑,質(zhì)量下降[16]。故分析冬小麥不同生育階段干旱發(fā)生規(guī)律對用戶選擇適應(yīng)模式意義重大。

    由圖2可以看出,邢臺縣冬小麥整個生育期SPI值呈現(xiàn)下降趨勢,說明降水量不斷下降,氣象干旱加劇。各生育階段SPI趨勢斜率大小比較:灌漿-成熟階段(-0.011)<拔節(jié)階段(-0.008)<苗生長階段(-0.007)<返青階段(-0.003)<越冬階段(-0.002)<播種前、抽穗-開花階段(0.010),表明降水減少速率在灌漿-成熟期最大,拔節(jié)期次之,再之是苗期,返青期和越冬期最小。冬小麥在苗生長階段、拔節(jié)階段、成熟階段氣象干旱加劇較其它生育階段明顯。統(tǒng)計不同生育階段干旱發(fā)生次數(shù)可知,58年間發(fā)生中旱及以上程度干旱,越冬期8次,苗初期9次,苗中期8次,苗后期0次,越冬期平均2次,返青期6次,拔節(jié)期7次,抽穗-開花期5次,灌漿-成熟期6次。除越冬期之外,其他各個生育期均受到干旱威脅。

    運用M-K非參數(shù)檢驗方法對邢臺縣冬小麥各生育階段58年SPI指數(shù)作具體趨勢突變檢驗(圖3),檢驗結(jié)果可以總結(jié)出:

    圖3a冬小麥播種前,主要突變年份是1963年、1978年,SPI指數(shù)在1956-1960年和1995-2003年呈下降趨勢,處于相對干旱階段;圖3b苗生長階段,突變年份是1963年、1969年、1982年和1993年。SPI指數(shù)在1954-1960年、1963-1967年、1969-1973年和1981-2011年呈下降態(tài)勢,表明干旱趨勢加強;圖3c越冬階段,除1979年外,UFk值均小于0,表明冬小麥越冬階段在整個時間序列中持續(xù)呈現(xiàn)干旱趨勢,說明邢臺縣受冬旱影響加重;圖3d返青階段,突變年份是1959年、1979年和2000年,SPI指數(shù)在1960-1962年、1968-1988年和2001-2011年呈下降趨勢,表明干旱趨勢加強;圖3e拔節(jié)階段,主要突變年份為1959年、1961、1978年、2005年,SPI指數(shù)在1974-2011年呈下降趨勢,表明冬小麥拔節(jié)階段在整個時間序列中呈現(xiàn)干旱趨勢;圖3f抽穗-開花階段,主要突變年份為1956年、1958年、1961年、1962年和1964年,在1955-1963年、1966-1969年和1972-1984年,SPI指數(shù)下降,干旱趨勢增強;圖3g灌漿-成熟階段,主要突變年份為1956年、1959年、1963年、1973年、1983年和1985年,1960-1970年和1974-2011年SPI指數(shù)下降,干旱趨勢增強。

    2.3夏玉米不同生育階段干旱變化

    夏玉米生育期為6月下旬至10月初,不同生育期起止時間見表2,夏玉米SPI變化規(guī)律見圖4。夏玉米主要生長季為邢臺縣的雨季。夏玉米全生育期需水量在350~400 mm,而平均降水量約為440 mm,但是各生育階段降水差異大,影響其生長發(fā)育。播種前,干旱發(fā)生,土壤底墑不足,影響夏玉米出苗;播種-拔節(jié)期,此階段為幼苗生長發(fā)育階段,需水量占全生育期需水量的20%,相對較少,但此時干旱發(fā)生影響幼苗生長,影響后期生殖生長;拔節(jié)-抽雄期,干旱對夏玉米的株高和葉面積生長影響最大;抽雄-灌漿期,夏玉米進入生殖生長階段,干旱發(fā)生將對夏玉米的穗數(shù)產(chǎn)生一定的影響;灌漿-成熟期,干旱發(fā)生,對夏玉米顆數(shù)和粒徑大小都將產(chǎn)生影響,加劇夏玉米減產(chǎn)[14]。因此要著重關(guān)注邢臺縣歷年夏玉米不同生育期的干旱規(guī)律,以便于保證夏糧生產(chǎn)的順利進行。

    運用M-K非參數(shù)檢驗方法對邢臺縣夏玉米各生育階段58年SPI指數(shù)做具體趨勢突變檢驗,見圖5。由圖5可以看出,夏玉米不同生育期SPI斜率比較,播種前、灌漿-成熟階段(0.010)>苗生長階段、抽雄-灌漿階段(-0.011)>拔節(jié)-抽雄階段(-0.018)。SPI值在夏玉米生長的早中期呈下降的趨勢,后期的灌漿成熟階段階段反而呈現(xiàn)出上升趨勢,說明在此階段干旱對夏玉米生長的危險有所下降。統(tǒng)計58年間夏玉米不同生育期的中旱以上干旱共發(fā)生35次,其中,播種前一月中旱5次,播種-拔節(jié)階段中旱5次、重旱1次,拔節(jié)-抽雄階段中旱3次、重旱3次,抽雄-灌漿階段中旱8次、重旱2次,灌漿-成熟階段中旱5次、重旱3次。由統(tǒng)計數(shù)據(jù)可知,播種-拔節(jié)階段、拔節(jié)-抽雄階段是干旱的頻發(fā)期,降水量少,影響夏玉米穩(wěn)定需水量的供給。故在拔節(jié)-抽雄開花生長階段要重點監(jiān)測和防范。

    圖5a和圖5b分別是播種前和播種-拔節(jié)階段,該兩個階段同時也是冬小麥抽雄-成熟階段,已在冬小麥干旱變化突變檢驗中進行分析,故在此不做具體分析;圖5c拔節(jié)-抽雄階段,主要突變年份有1972年、1975年、1990年、2000年、2002年和2006年,SPI指數(shù)在1970-2011年呈下降趨勢,夏玉米在該階段干旱化持續(xù)加強;圖5d抽雄-灌漿階段,主要的突變年份是1957年,1958-1962年、1964-1974年和1977-2011年,干旱趨勢加強,尤其在2001-2008年,干旱趨勢通過了0.05臨界線的顯著性水平,表明降水減少,氣候變干的趨勢十分顯著;圖5e灌漿-成熟階段,主要突變年份是1961年、1964年、1974年和1978年,1956-1960年、1995-2003年,SPI指數(shù)呈下降趨勢,有干旱化傾向。

    3小結(jié)與討論

    1)根據(jù)SPI12變化規(guī)律,總結(jié)出邢臺縣農(nóng)作物氣象干旱逐漸加重,以中旱和重旱為主,1980年代后期以來干旱發(fā)生的頻次明顯增加。不同年代干旱變化趨勢不同,20世紀50年代后期下降最快,干旱趨勢明顯;60年代和90年代呈微下降趨勢;20世紀80年代上升最快,濕澇較明顯;20世紀70年代和21世紀初呈微上升趨勢。

    2)邢臺縣冬小麥不同生育階段干旱發(fā)生規(guī)律:冬小麥整個生育期SPI指數(shù)呈下降趨勢,表明氣象干旱加劇,其中,灌漿-成熟期下降速率最大,拔節(jié)期次之,再之是苗期。冬小麥在苗期、拔節(jié)期,成熟期氣象干旱加劇較其它生育期明顯。突變分析表明,除播種前和抽穗-開花階段,其他各生育階段SPI指數(shù)自1980年代初均有不同程度下降,干旱趨勢相對明顯,與邢臺縣整體干旱變化相一致。

    3)邢臺縣夏玉米不同生育階段干旱發(fā)生規(guī)律:SPI顯示夏玉米生育期的早中期階段呈下降的趨勢,播種-拔節(jié)階段、拔節(jié)-抽雄階段是干旱的頻發(fā)期,降水少,影響夏玉米穩(wěn)定需水量的供給。突變分析表明,除播種前和灌漿-成熟階段,其他各生育階段SPI指數(shù)自20世紀70年代后期呈下降趨勢,干旱發(fā)生相對增加與與邢臺縣整體干旱變化大致一致。

    4)本研究根據(jù)1954-2011年58年降水?dāng)?shù)據(jù)運用SPI指數(shù)進行了冬小麥-夏玉米生育期干旱發(fā)生規(guī)律研究,為農(nóng)戶適時進行灌溉提供了科學(xué)依據(jù)。河北農(nóng)村統(tǒng)計資料表明2010年邢臺縣發(fā)生大范圍的干旱,而年尺度的SPI顯示2010分析2010年以年為尺度的各月累積降水可知,各月降水量較多,而SPI12值是以年為單位的累積降水進行標準化得到的,所以2010年SPI12顯示未有明顯干旱。同時本研究采用的是縣域尺度研究區(qū),氣象站點較少,沒有系統(tǒng)進行邢臺縣冬小麥-夏玉米生育期干旱的空間變化研究。同時研究一個地區(qū)農(nóng)作物生育期干旱變化,不僅要分析降水因子,還要進一步分析降水之外的其它氣象因子、下墊面因子等。

    參考文獻:

    [1] 張強,鞠笑生,李淑華.Palmer干旱指數(shù)、地表濕潤指數(shù)與降水距平指數(shù)的比較[J].地理學(xué)報,2003,58:117-124.

    [2] GUTTMAN N B. Comparing the Palmer drought index and the standardized precipitation index[J]. Journal of the American Water Resources Association,1998,34(1):113-121.

    [3] 車少靜,李春強,申雙和. 基于SPI的近41年河北省旱澇時空特征分析[J].中國農(nóng)業(yè)氣象,2010,31(1):137-143.

    [4] 郁凌華,趙艷霞.黃淮海地區(qū)夏玉米生長季內(nèi)的旱澇災(zāi)害分析[J].災(zāi)害學(xué),2013,28(2):71-75.

    [5] 陳莉,方麗娟,李帥.東北地區(qū)近50年農(nóng)作物生長季干旱趨勢研究[J].災(zāi)害學(xué),2010,25(4):5-10.

    [6] 楊麗娜,吳智杰,許新路,等.基于標準化降水指數(shù)的邢臺地區(qū)的干旱特征[J].干旱區(qū)資源與環(huán)境,2011,25(11):120-124.

    [7] 商彥蕊,史培軍.自然災(zāi)害系統(tǒng)脆弱性研究[M].西安:西安地圖出版社,2004.191-198.

    [8] MCKEE T B, DOESKEN N J, KLEIST J. The relationship of drought fre-quency and duration to time scales [A]. In: Proc.8th conference onapplied climatology[C]. Anaheim,California,1993: 179-184.

    [9] EDWARDS D C, MCKEE T B. Characteristics of 20th century drought in the United States atmultiple scale[R]. Atmospheric Science PaperNo.634,1997:1-30.

    [10] MCKEE T B,DOESKEN N J,KLEIST J.The relationship of drought frequency and duration to time scales[R].California: Eighth Conference on Applied Climatology,1993.

    [11] 陳瑩,陳興偉.福建省近50年干旱時空特征演變_基于標準化降水指數(shù)分析J].2011,20(3):57-63.

    [12] 王鶯,趙福年,姚玉璧,等.基于Z指數(shù)的石羊河流域干旱特征分析[J].災(zāi)害學(xué),2013,28(2):100-106.

    [13] 魏鳳英.現(xiàn)代氣象統(tǒng)計診斷與預(yù)測技術(shù)[M].北京:氣象出版社,2011:69-72.

    [14] 肖俊夫,劉占東,劉祖貴,等.不同時期干旱和干旱程度對夏玉米生長發(fā)育及耗水特性的影響[J].玉米科學(xué),2011,19(4):54-58,64.

    [15] 吳冬麗,王春已,薛紅喜,等.華北地區(qū)冬小麥干旱時空分布特征[J].自然災(zāi)害學(xué)報,2012,21(1):18-25.

    [16] 吳東麗,王春乙,張雪芬,等.華北冬小麥作物氣候干旱指數(shù)研究[J].科技導(dǎo)報,2009,27(7):32-36.

    運用M-K非參數(shù)檢驗方法對邢臺縣夏玉米各生育階段58年SPI指數(shù)做具體趨勢突變檢驗,見圖5。由圖5可以看出,夏玉米不同生育期SPI斜率比較,播種前、灌漿-成熟階段(0.010)>苗生長階段、抽雄-灌漿階段(-0.011)>拔節(jié)-抽雄階段(-0.018)。SPI值在夏玉米生長的早中期呈下降的趨勢,后期的灌漿成熟階段階段反而呈現(xiàn)出上升趨勢,說明在此階段干旱對夏玉米生長的危險有所下降。統(tǒng)計58年間夏玉米不同生育期的中旱以上干旱共發(fā)生35次,其中,播種前一月中旱5次,播種-拔節(jié)階段中旱5次、重旱1次,拔節(jié)-抽雄階段中旱3次、重旱3次,抽雄-灌漿階段中旱8次、重旱2次,灌漿-成熟階段中旱5次、重旱3次。由統(tǒng)計數(shù)據(jù)可知,播種-拔節(jié)階段、拔節(jié)-抽雄階段是干旱的頻發(fā)期,降水量少,影響夏玉米穩(wěn)定需水量的供給。故在拔節(jié)-抽雄開花生長階段要重點監(jiān)測和防范。

    圖5a和圖5b分別是播種前和播種-拔節(jié)階段,該兩個階段同時也是冬小麥抽雄-成熟階段,已在冬小麥干旱變化突變檢驗中進行分析,故在此不做具體分析;圖5c拔節(jié)-抽雄階段,主要突變年份有1972年、1975年、1990年、2000年、2002年和2006年,SPI指數(shù)在1970-2011年呈下降趨勢,夏玉米在該階段干旱化持續(xù)加強;圖5d抽雄-灌漿階段,主要的突變年份是1957年,1958-1962年、1964-1974年和1977-2011年,干旱趨勢加強,尤其在2001-2008年,干旱趨勢通過了0.05臨界線的顯著性水平,表明降水減少,氣候變干的趨勢十分顯著;圖5e灌漿-成熟階段,主要突變年份是1961年、1964年、1974年和1978年,1956-1960年、1995-2003年,SPI指數(shù)呈下降趨勢,有干旱化傾向。

    3小結(jié)與討論

    1)根據(jù)SPI12變化規(guī)律,總結(jié)出邢臺縣農(nóng)作物氣象干旱逐漸加重,以中旱和重旱為主,1980年代后期以來干旱發(fā)生的頻次明顯增加。不同年代干旱變化趨勢不同,20世紀50年代后期下降最快,干旱趨勢明顯;60年代和90年代呈微下降趨勢;20世紀80年代上升最快,濕澇較明顯;20世紀70年代和21世紀初呈微上升趨勢。

    2)邢臺縣冬小麥不同生育階段干旱發(fā)生規(guī)律:冬小麥整個生育期SPI指數(shù)呈下降趨勢,表明氣象干旱加劇,其中,灌漿-成熟期下降速率最大,拔節(jié)期次之,再之是苗期。冬小麥在苗期、拔節(jié)期,成熟期氣象干旱加劇較其它生育期明顯。突變分析表明,除播種前和抽穗-開花階段,其他各生育階段SPI指數(shù)自1980年代初均有不同程度下降,干旱趨勢相對明顯,與邢臺縣整體干旱變化相一致。

    3)邢臺縣夏玉米不同生育階段干旱發(fā)生規(guī)律:SPI顯示夏玉米生育期的早中期階段呈下降的趨勢,播種-拔節(jié)階段、拔節(jié)-抽雄階段是干旱的頻發(fā)期,降水少,影響夏玉米穩(wěn)定需水量的供給。突變分析表明,除播種前和灌漿-成熟階段,其他各生育階段SPI指數(shù)自20世紀70年代后期呈下降趨勢,干旱發(fā)生相對增加與與邢臺縣整體干旱變化大致一致。

    4)本研究根據(jù)1954-2011年58年降水?dāng)?shù)據(jù)運用SPI指數(shù)進行了冬小麥-夏玉米生育期干旱發(fā)生規(guī)律研究,為農(nóng)戶適時進行灌溉提供了科學(xué)依據(jù)。河北農(nóng)村統(tǒng)計資料表明2010年邢臺縣發(fā)生大范圍的干旱,而年尺度的SPI顯示2010分析2010年以年為尺度的各月累積降水可知,各月降水量較多,而SPI12值是以年為單位的累積降水進行標準化得到的,所以2010年SPI12顯示未有明顯干旱。同時本研究采用的是縣域尺度研究區(qū),氣象站點較少,沒有系統(tǒng)進行邢臺縣冬小麥-夏玉米生育期干旱的空間變化研究。同時研究一個地區(qū)農(nóng)作物生育期干旱變化,不僅要分析降水因子,還要進一步分析降水之外的其它氣象因子、下墊面因子等。

    參考文獻:

    [1] 張強,鞠笑生,李淑華.Palmer干旱指數(shù)、地表濕潤指數(shù)與降水距平指數(shù)的比較[J].地理學(xué)報,2003,58:117-124.

    [2] GUTTMAN N B. Comparing the Palmer drought index and the standardized precipitation index[J]. Journal of the American Water Resources Association,1998,34(1):113-121.

    [3] 車少靜,李春強,申雙和. 基于SPI的近41年河北省旱澇時空特征分析[J].中國農(nóng)業(yè)氣象,2010,31(1):137-143.

    [4] 郁凌華,趙艷霞.黃淮海地區(qū)夏玉米生長季內(nèi)的旱澇災(zāi)害分析[J].災(zāi)害學(xué),2013,28(2):71-75.

    [5] 陳莉,方麗娟,李帥.東北地區(qū)近50年農(nóng)作物生長季干旱趨勢研究[J].災(zāi)害學(xué),2010,25(4):5-10.

    [6] 楊麗娜,吳智杰,許新路,等.基于標準化降水指數(shù)的邢臺地區(qū)的干旱特征[J].干旱區(qū)資源與環(huán)境,2011,25(11):120-124.

    [7] 商彥蕊,史培軍.自然災(zāi)害系統(tǒng)脆弱性研究[M].西安:西安地圖出版社,2004.191-198.

    [8] MCKEE T B, DOESKEN N J, KLEIST J. The relationship of drought fre-quency and duration to time scales [A]. In: Proc.8th conference onapplied climatology[C]. Anaheim,California,1993: 179-184.

    [9] EDWARDS D C, MCKEE T B. Characteristics of 20th century drought in the United States atmultiple scale[R]. Atmospheric Science PaperNo.634,1997:1-30.

    [10] MCKEE T B,DOESKEN N J,KLEIST J.The relationship of drought frequency and duration to time scales[R].California: Eighth Conference on Applied Climatology,1993.

    [11] 陳瑩,陳興偉.福建省近50年干旱時空特征演變_基于標準化降水指數(shù)分析J].2011,20(3):57-63.

    [12] 王鶯,趙福年,姚玉璧,等.基于Z指數(shù)的石羊河流域干旱特征分析[J].災(zāi)害學(xué),2013,28(2):100-106.

    [13] 魏鳳英.現(xiàn)代氣象統(tǒng)計診斷與預(yù)測技術(shù)[M].北京:氣象出版社,2011:69-72.

    [14] 肖俊夫,劉占東,劉祖貴,等.不同時期干旱和干旱程度對夏玉米生長發(fā)育及耗水特性的影響[J].玉米科學(xué),2011,19(4):54-58,64.

    [15] 吳冬麗,王春已,薛紅喜,等.華北地區(qū)冬小麥干旱時空分布特征[J].自然災(zāi)害學(xué)報,2012,21(1):18-25.

    [16] 吳東麗,王春乙,張雪芬,等.華北冬小麥作物氣候干旱指數(shù)研究[J].科技導(dǎo)報,2009,27(7):32-36.

    運用M-K非參數(shù)檢驗方法對邢臺縣夏玉米各生育階段58年SPI指數(shù)做具體趨勢突變檢驗,見圖5。由圖5可以看出,夏玉米不同生育期SPI斜率比較,播種前、灌漿-成熟階段(0.010)>苗生長階段、抽雄-灌漿階段(-0.011)>拔節(jié)-抽雄階段(-0.018)。SPI值在夏玉米生長的早中期呈下降的趨勢,后期的灌漿成熟階段階段反而呈現(xiàn)出上升趨勢,說明在此階段干旱對夏玉米生長的危險有所下降。統(tǒng)計58年間夏玉米不同生育期的中旱以上干旱共發(fā)生35次,其中,播種前一月中旱5次,播種-拔節(jié)階段中旱5次、重旱1次,拔節(jié)-抽雄階段中旱3次、重旱3次,抽雄-灌漿階段中旱8次、重旱2次,灌漿-成熟階段中旱5次、重旱3次。由統(tǒng)計數(shù)據(jù)可知,播種-拔節(jié)階段、拔節(jié)-抽雄階段是干旱的頻發(fā)期,降水量少,影響夏玉米穩(wěn)定需水量的供給。故在拔節(jié)-抽雄開花生長階段要重點監(jiān)測和防范。

    圖5a和圖5b分別是播種前和播種-拔節(jié)階段,該兩個階段同時也是冬小麥抽雄-成熟階段,已在冬小麥干旱變化突變檢驗中進行分析,故在此不做具體分析;圖5c拔節(jié)-抽雄階段,主要突變年份有1972年、1975年、1990年、2000年、2002年和2006年,SPI指數(shù)在1970-2011年呈下降趨勢,夏玉米在該階段干旱化持續(xù)加強;圖5d抽雄-灌漿階段,主要的突變年份是1957年,1958-1962年、1964-1974年和1977-2011年,干旱趨勢加強,尤其在2001-2008年,干旱趨勢通過了0.05臨界線的顯著性水平,表明降水減少,氣候變干的趨勢十分顯著;圖5e灌漿-成熟階段,主要突變年份是1961年、1964年、1974年和1978年,1956-1960年、1995-2003年,SPI指數(shù)呈下降趨勢,有干旱化傾向。

    3小結(jié)與討論

    1)根據(jù)SPI12變化規(guī)律,總結(jié)出邢臺縣農(nóng)作物氣象干旱逐漸加重,以中旱和重旱為主,1980年代后期以來干旱發(fā)生的頻次明顯增加。不同年代干旱變化趨勢不同,20世紀50年代后期下降最快,干旱趨勢明顯;60年代和90年代呈微下降趨勢;20世紀80年代上升最快,濕澇較明顯;20世紀70年代和21世紀初呈微上升趨勢。

    2)邢臺縣冬小麥不同生育階段干旱發(fā)生規(guī)律:冬小麥整個生育期SPI指數(shù)呈下降趨勢,表明氣象干旱加劇,其中,灌漿-成熟期下降速率最大,拔節(jié)期次之,再之是苗期。冬小麥在苗期、拔節(jié)期,成熟期氣象干旱加劇較其它生育期明顯。突變分析表明,除播種前和抽穗-開花階段,其他各生育階段SPI指數(shù)自1980年代初均有不同程度下降,干旱趨勢相對明顯,與邢臺縣整體干旱變化相一致。

    3)邢臺縣夏玉米不同生育階段干旱發(fā)生規(guī)律:SPI顯示夏玉米生育期的早中期階段呈下降的趨勢,播種-拔節(jié)階段、拔節(jié)-抽雄階段是干旱的頻發(fā)期,降水少,影響夏玉米穩(wěn)定需水量的供給。突變分析表明,除播種前和灌漿-成熟階段,其他各生育階段SPI指數(shù)自20世紀70年代后期呈下降趨勢,干旱發(fā)生相對增加與與邢臺縣整體干旱變化大致一致。

    4)本研究根據(jù)1954-2011年58年降水?dāng)?shù)據(jù)運用SPI指數(shù)進行了冬小麥-夏玉米生育期干旱發(fā)生規(guī)律研究,為農(nóng)戶適時進行灌溉提供了科學(xué)依據(jù)。河北農(nóng)村統(tǒng)計資料表明2010年邢臺縣發(fā)生大范圍的干旱,而年尺度的SPI顯示2010分析2010年以年為尺度的各月累積降水可知,各月降水量較多,而SPI12值是以年為單位的累積降水進行標準化得到的,所以2010年SPI12顯示未有明顯干旱。同時本研究采用的是縣域尺度研究區(qū),氣象站點較少,沒有系統(tǒng)進行邢臺縣冬小麥-夏玉米生育期干旱的空間變化研究。同時研究一個地區(qū)農(nóng)作物生育期干旱變化,不僅要分析降水因子,還要進一步分析降水之外的其它氣象因子、下墊面因子等。

    參考文獻:

    [1] 張強,鞠笑生,李淑華.Palmer干旱指數(shù)、地表濕潤指數(shù)與降水距平指數(shù)的比較[J].地理學(xué)報,2003,58:117-124.

    [2] GUTTMAN N B. Comparing the Palmer drought index and the standardized precipitation index[J]. Journal of the American Water Resources Association,1998,34(1):113-121.

    [3] 車少靜,李春強,申雙和. 基于SPI的近41年河北省旱澇時空特征分析[J].中國農(nóng)業(yè)氣象,2010,31(1):137-143.

    [4] 郁凌華,趙艷霞.黃淮海地區(qū)夏玉米生長季內(nèi)的旱澇災(zāi)害分析[J].災(zāi)害學(xué),2013,28(2):71-75.

    [5] 陳莉,方麗娟,李帥.東北地區(qū)近50年農(nóng)作物生長季干旱趨勢研究[J].災(zāi)害學(xué),2010,25(4):5-10.

    [6] 楊麗娜,吳智杰,許新路,等.基于標準化降水指數(shù)的邢臺地區(qū)的干旱特征[J].干旱區(qū)資源與環(huán)境,2011,25(11):120-124.

    [7] 商彥蕊,史培軍.自然災(zāi)害系統(tǒng)脆弱性研究[M].西安:西安地圖出版社,2004.191-198.

    [8] MCKEE T B, DOESKEN N J, KLEIST J. The relationship of drought fre-quency and duration to time scales [A]. In: Proc.8th conference onapplied climatology[C]. Anaheim,California,1993: 179-184.

    [9] EDWARDS D C, MCKEE T B. Characteristics of 20th century drought in the United States atmultiple scale[R]. Atmospheric Science PaperNo.634,1997:1-30.

    [10] MCKEE T B,DOESKEN N J,KLEIST J.The relationship of drought frequency and duration to time scales[R].California: Eighth Conference on Applied Climatology,1993.

    [11] 陳瑩,陳興偉.福建省近50年干旱時空特征演變_基于標準化降水指數(shù)分析J].2011,20(3):57-63.

    [12] 王鶯,趙福年,姚玉璧,等.基于Z指數(shù)的石羊河流域干旱特征分析[J].災(zāi)害學(xué),2013,28(2):100-106.

    [13] 魏鳳英.現(xiàn)代氣象統(tǒng)計診斷與預(yù)測技術(shù)[M].北京:氣象出版社,2011:69-72.

    [14] 肖俊夫,劉占東,劉祖貴,等.不同時期干旱和干旱程度對夏玉米生長發(fā)育及耗水特性的影響[J].玉米科學(xué),2011,19(4):54-58,64.

    [15] 吳冬麗,王春已,薛紅喜,等.華北地區(qū)冬小麥干旱時空分布特征[J].自然災(zāi)害學(xué)報,2012,21(1):18-25.

    [16] 吳東麗,王春乙,張雪芬,等.華北冬小麥作物氣候干旱指數(shù)研究[J].科技導(dǎo)報,2009,27(7):32-36.

    猜你喜歡
    邢臺縣干旱夏玉米
    小麥收割之后 如何種植夏玉米才能高產(chǎn)
    美術(shù)作品
    夏玉米高產(chǎn)的關(guān)鍵栽培技術(shù)措施
    基于距平的白城地區(qū)干旱時間分布特征分析
    臨夏地區(qū)干旱特征及干濕氣候區(qū)劃
    夏季高溫干旱時節(jié)高山蔬菜種植管理策略
    小說里的世界 虛擬里的真實
    種植食用玫瑰 畝產(chǎn)收入六千
    歷史文化名村保護的可持續(xù)發(fā)展策略——以河北省邢臺縣英談村為例
    天津市2011年良種補貼推介品種目錄
    曰老女人黄片| www.精华液| 亚洲18禁久久av| 精品福利观看| 99久久成人亚洲精品观看| 18禁裸乳无遮挡免费网站照片| 国产av不卡久久| 村上凉子中文字幕在线| 91在线观看av| av片东京热男人的天堂| 亚洲av电影在线进入| 久久久国产欧美日韩av| 白带黄色成豆腐渣| 日日夜夜操网爽| 日本黄色片子视频| 午夜精品在线福利| 露出奶头的视频| 午夜a级毛片| 18禁观看日本| 岛国在线观看网站| 怎么达到女性高潮| 首页视频小说图片口味搜索| 91在线精品国自产拍蜜月 | 一个人看视频在线观看www免费 | 亚洲成a人片在线一区二区| 亚洲七黄色美女视频| 亚洲国产欧洲综合997久久,| 麻豆久久精品国产亚洲av| 国产1区2区3区精品| 国产v大片淫在线免费观看| 免费人成视频x8x8入口观看| 亚洲,欧美精品.| 国产三级在线视频| 国产精品影院久久| 免费看光身美女| 亚洲第一电影网av| 久久久国产成人免费| 色综合亚洲欧美另类图片| av福利片在线观看| 韩国av一区二区三区四区| 真实男女啪啪啪动态图| 97超级碰碰碰精品色视频在线观看| 搡老妇女老女人老熟妇| 97超视频在线观看视频| 搡老岳熟女国产| 国产午夜福利久久久久久| 国产亚洲精品一区二区www| 国产成人av激情在线播放| 久久久久性生活片| 一本久久中文字幕| 国产成人精品久久二区二区91| 精品国产乱码久久久久久男人| 国产欧美日韩精品亚洲av| 日韩 欧美 亚洲 中文字幕| 久久精品人妻少妇| 99久久久亚洲精品蜜臀av| 国产精品久久久久久久电影 | 欧美中文综合在线视频| 亚洲欧美激情综合另类| 欧美xxxx黑人xx丫x性爽| 亚洲国产欧洲综合997久久,| 九色成人免费人妻av| 日本撒尿小便嘘嘘汇集6| 国产亚洲精品久久久久久毛片| a级毛片在线看网站| 男人舔女人下体高潮全视频| 黑人欧美特级aaaaaa片| 禁无遮挡网站| 亚洲成人久久性| 一级毛片精品| 久久久国产成人精品二区| 成人午夜高清在线视频| 97超视频在线观看视频| 久久天堂一区二区三区四区| 国产真人三级小视频在线观看| 亚洲真实伦在线观看| 一区二区三区高清视频在线| 丁香六月欧美| 黄色丝袜av网址大全| 久久伊人香网站| 国产久久久一区二区三区| 国内精品久久久久久久电影| 99riav亚洲国产免费| 男人的好看免费观看在线视频| 中文字幕熟女人妻在线| 亚洲,欧美精品.| 久久精品国产清高在天天线| 熟女电影av网| 香蕉久久夜色| 99热这里只有是精品50| 精品一区二区三区视频在线 | 久久精品国产99精品国产亚洲性色| 露出奶头的视频| 人人妻人人澡欧美一区二区| 国产成+人综合+亚洲专区| 看片在线看免费视频| 久久久国产成人精品二区| 两性夫妻黄色片| 俄罗斯特黄特色一大片| 人人妻人人澡欧美一区二区| 日本黄大片高清| 真人做人爱边吃奶动态| 成人亚洲精品av一区二区| 久久性视频一级片| 国产高清有码在线观看视频| 国产精品久久久av美女十八| 12—13女人毛片做爰片一| 婷婷精品国产亚洲av| 9191精品国产免费久久| 国产久久久一区二区三区| 亚洲欧美激情综合另类| 国产精品久久视频播放| 91av网站免费观看| 久久久水蜜桃国产精品网| 精品久久久久久,| 99视频精品全部免费 在线 | 欧美又色又爽又黄视频| 成人高潮视频无遮挡免费网站| 亚洲欧美日韩高清在线视频| 免费观看的影片在线观看| 一个人观看的视频www高清免费观看 | 在线十欧美十亚洲十日本专区| 欧美中文综合在线视频| 夜夜夜夜夜久久久久| 男插女下体视频免费在线播放| 国产精品精品国产色婷婷| 香蕉av资源在线| 免费看a级黄色片| 白带黄色成豆腐渣| 亚洲中文日韩欧美视频| 亚洲欧洲精品一区二区精品久久久| av片东京热男人的天堂| 亚洲成av人片免费观看| 精品国产超薄肉色丝袜足j| 色av中文字幕| 欧美色欧美亚洲另类二区| 欧美中文日本在线观看视频| 美女高潮的动态| 一级作爱视频免费观看| 国产精品 国内视频| 叶爱在线成人免费视频播放| 欧美日韩中文字幕国产精品一区二区三区| 中文字幕熟女人妻在线| 在线永久观看黄色视频| 黄色 视频免费看| 亚洲国产日韩欧美精品在线观看 | 亚洲七黄色美女视频| 国产精品av久久久久免费| 亚洲熟妇熟女久久| 1024香蕉在线观看| 国产又色又爽无遮挡免费看| 在线视频色国产色| 亚洲成人久久性| 久久精品91蜜桃| 日韩精品中文字幕看吧| 啦啦啦韩国在线观看视频| 亚洲男人的天堂狠狠| 国产视频内射| 精品无人区乱码1区二区| 国产精品 国内视频| 久久久久国内视频| 免费大片18禁| 99国产综合亚洲精品| 一个人免费在线观看的高清视频| 琪琪午夜伦伦电影理论片6080| 国产精品久久久人人做人人爽| 中文字幕人成人乱码亚洲影| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲专区中文字幕在线| 成人国产一区最新在线观看| 91字幕亚洲| 欧美一区二区精品小视频在线| 99久久久亚洲精品蜜臀av| 欧美日韩精品网址| 欧美精品啪啪一区二区三区| 国产伦一二天堂av在线观看| 欧美高清成人免费视频www| 精品国内亚洲2022精品成人| 国产精品一区二区免费欧美| 亚洲天堂国产精品一区在线| 欧美中文综合在线视频| svipshipincom国产片| 免费在线观看成人毛片| av片东京热男人的天堂| 欧美中文日本在线观看视频| svipshipincom国产片| 蜜桃久久精品国产亚洲av| 男女那种视频在线观看| 亚洲av成人精品一区久久| 精品国产乱子伦一区二区三区| 国产又黄又爽又无遮挡在线| 久久人妻av系列| 国产v大片淫在线免费观看| 精品熟女少妇八av免费久了| 三级毛片av免费| 国产单亲对白刺激| 一级作爱视频免费观看| 久久这里只有精品19| 日本撒尿小便嘘嘘汇集6| 在线免费观看不下载黄p国产 | 淫秽高清视频在线观看| 国产精品免费一区二区三区在线| 女生性感内裤真人,穿戴方法视频| 亚洲熟妇熟女久久| 国产亚洲欧美在线一区二区| 丰满的人妻完整版| 男女下面进入的视频免费午夜| 国产一区二区在线观看日韩 | 亚洲成av人片在线播放无| 欧美一级a爱片免费观看看| 日韩大尺度精品在线看网址| 窝窝影院91人妻| 韩国av一区二区三区四区| 国产成人精品久久二区二区免费| 啦啦啦韩国在线观看视频| 中文字幕高清在线视频| 国产精品99久久99久久久不卡| 国产毛片a区久久久久| 欧美激情久久久久久爽电影| 欧美日韩瑟瑟在线播放| 国产精品亚洲一级av第二区| 亚洲中文字幕一区二区三区有码在线看 | 国产午夜福利久久久久久| 国产综合懂色| 久久久久国产一级毛片高清牌| 美女扒开内裤让男人捅视频| 精品国产乱子伦一区二区三区| 欧美精品啪啪一区二区三区| 国产极品精品免费视频能看的| 亚洲电影在线观看av| 日韩欧美国产一区二区入口| 91av网一区二区| 久久久久久久久久黄片| 一夜夜www| 老鸭窝网址在线观看| 日韩高清综合在线| 精品熟女少妇八av免费久了| 免费在线观看亚洲国产| 亚洲avbb在线观看| 国产成人av激情在线播放| 韩国av一区二区三区四区| 欧美不卡视频在线免费观看| 精品一区二区三区视频在线观看免费| ponron亚洲| av中文乱码字幕在线| 在线观看一区二区三区| 久久天堂一区二区三区四区| 淫秽高清视频在线观看| 久久久久久国产a免费观看| avwww免费| 婷婷六月久久综合丁香| 免费在线观看影片大全网站| 波多野结衣高清无吗| 无遮挡黄片免费观看| 97超视频在线观看视频| 国产精品自产拍在线观看55亚洲| 久久天躁狠狠躁夜夜2o2o| 国产亚洲精品久久久com| 午夜精品久久久久久毛片777| 两人在一起打扑克的视频| 欧美极品一区二区三区四区| www.精华液| 午夜精品一区二区三区免费看| 日韩欧美在线乱码| 国产精品99久久99久久久不卡| 精品电影一区二区在线| 少妇人妻一区二区三区视频| АⅤ资源中文在线天堂| 亚洲精品粉嫩美女一区| 久久精品综合一区二区三区| 99国产综合亚洲精品| 午夜福利在线观看吧| 国产欧美日韩一区二区精品| 又爽又黄无遮挡网站| 日韩欧美国产在线观看| 久久久久久久久中文| 日韩中文字幕欧美一区二区| 亚洲美女黄片视频| 在线视频色国产色| 亚洲精品美女久久av网站| 国产精品av久久久久免费| 免费看十八禁软件| 欧美中文日本在线观看视频| 一本精品99久久精品77| 日日夜夜操网爽| 天天一区二区日本电影三级| 中文亚洲av片在线观看爽| av国产免费在线观看| 淫秽高清视频在线观看| 免费看十八禁软件| 日本撒尿小便嘘嘘汇集6| 久久伊人香网站| 亚洲成av人片在线播放无| 国产精品98久久久久久宅男小说| 欧美+亚洲+日韩+国产| 搡老岳熟女国产| 国产精品免费一区二区三区在线| 午夜福利免费观看在线| 色播亚洲综合网| 国产69精品久久久久777片 | 国产精品亚洲一级av第二区| 欧美又色又爽又黄视频| 国产精品久久久人人做人人爽| 麻豆国产97在线/欧美| 国模一区二区三区四区视频 | 麻豆成人av在线观看| 12—13女人毛片做爰片一| 国产不卡一卡二| 国产精品日韩av在线免费观看| 亚洲精华国产精华精| 99久久精品国产亚洲精品| 精品99又大又爽又粗少妇毛片 | 欧美黄色淫秽网站| 国产又黄又爽又无遮挡在线| 国产亚洲精品综合一区在线观看| 亚洲欧美精品综合一区二区三区| 久久香蕉精品热| 男人舔女人下体高潮全视频| 亚洲av日韩精品久久久久久密| 免费人成视频x8x8入口观看| 在线观看日韩欧美| 两个人看的免费小视频| 国产精品乱码一区二三区的特点| 12—13女人毛片做爰片一| www日本黄色视频网| or卡值多少钱| 操出白浆在线播放| 精品久久蜜臀av无| 精品久久久久久久人妻蜜臀av| 色哟哟哟哟哟哟| 亚洲国产精品sss在线观看| or卡值多少钱| 久久久色成人| 日本黄色片子视频| 中文字幕最新亚洲高清| 日韩欧美在线乱码| 亚洲色图av天堂| 人人妻人人澡欧美一区二区| 亚洲精品色激情综合| 国产伦人伦偷精品视频| 黄色片一级片一级黄色片| 亚洲在线观看片| 成人特级黄色片久久久久久久| 久久精品夜夜夜夜夜久久蜜豆| 国内毛片毛片毛片毛片毛片| 操出白浆在线播放| 人人妻人人看人人澡| 999精品在线视频| 国产蜜桃级精品一区二区三区| 狠狠狠狠99中文字幕| 又爽又黄无遮挡网站| 国产av麻豆久久久久久久| 国产一区二区三区视频了| 免费av毛片视频| 欧美黑人巨大hd| 婷婷丁香在线五月| 亚洲最大成人中文| 亚洲国产精品sss在线观看| 国内毛片毛片毛片毛片毛片| 免费在线观看亚洲国产| 亚洲国产欧美网| 国内少妇人妻偷人精品xxx网站 | 嫁个100分男人电影在线观看| 亚洲国产中文字幕在线视频| 搞女人的毛片| 日本免费一区二区三区高清不卡| 十八禁网站免费在线| 日本 欧美在线| 成人鲁丝片一二三区免费| 日日夜夜操网爽| 精品电影一区二区在线| 国产精品一及| 嫁个100分男人电影在线观看| av黄色大香蕉| 日韩中文字幕欧美一区二区| 国产精品久久久人人做人人爽| 国产精品 国内视频| 黄色丝袜av网址大全| 黄色日韩在线| avwww免费| 色播亚洲综合网| 精品日产1卡2卡| 搡老熟女国产l中国老女人| 国产视频一区二区在线看| 国产伦精品一区二区三区视频9 | 一区二区三区激情视频| 精品国产三级普通话版| 狠狠狠狠99中文字幕| 成年女人毛片免费观看观看9| 国产伦在线观看视频一区| 久久国产乱子伦精品免费另类| 日本免费a在线| 国产精品久久电影中文字幕| 不卡av一区二区三区| av在线蜜桃| 国产精品av视频在线免费观看| 欧美激情在线99| 日韩欧美国产在线观看| 黄色日韩在线| 88av欧美| 成人性生交大片免费视频hd| 亚洲 欧美一区二区三区| 欧美+亚洲+日韩+国产| 嫩草影院精品99| 欧美中文日本在线观看视频| 色老头精品视频在线观看| 国产不卡一卡二| 最近最新免费中文字幕在线| 日韩欧美在线二视频| 熟女电影av网| 老汉色av国产亚洲站长工具| 国内精品久久久久久久电影| 亚洲国产日韩欧美精品在线观看 | 网址你懂的国产日韩在线| 久久这里只有精品中国| 午夜福利在线观看吧| 久久午夜综合久久蜜桃| 一区二区三区高清视频在线| 国产一区二区三区在线臀色熟女| 国产美女午夜福利| 特级一级黄色大片| 国内揄拍国产精品人妻在线| 亚洲精华国产精华精| xxx96com| 男人舔奶头视频| 国产高清有码在线观看视频| 亚洲欧洲精品一区二区精品久久久| 亚洲中文字幕日韩| 深夜精品福利| 在线看三级毛片| 又黄又粗又硬又大视频| 久久精品aⅴ一区二区三区四区| 俺也久久电影网| 人妻丰满熟妇av一区二区三区| 伊人久久大香线蕉亚洲五| 日日干狠狠操夜夜爽| 欧洲精品卡2卡3卡4卡5卡区| av在线蜜桃| 一夜夜www| 亚洲第一电影网av| 日本黄大片高清| 嫩草影院精品99| 伦理电影免费视频| 久久久久久人人人人人| 叶爱在线成人免费视频播放| 一卡2卡三卡四卡精品乱码亚洲| 成人特级黄色片久久久久久久| 亚洲精品美女久久久久99蜜臀| 国产真人三级小视频在线观看| 成人午夜高清在线视频| 国产 一区 欧美 日韩| 久久精品人妻少妇| 99久久无色码亚洲精品果冻| 久久久精品大字幕| 日韩欧美 国产精品| 亚洲无线在线观看| 成人亚洲精品av一区二区| 亚洲国产高清在线一区二区三| 成年女人永久免费观看视频| 狂野欧美激情性xxxx| 熟女电影av网| 久久久精品欧美日韩精品| 欧美成人免费av一区二区三区| 女同久久另类99精品国产91| 国产午夜精品久久久久久| 成人亚洲精品av一区二区| 欧美三级亚洲精品| 欧美日韩乱码在线| 国产精品一区二区免费欧美| 亚洲美女黄片视频| 日韩高清综合在线| 国产熟女xx| 宅男免费午夜| 亚洲一区高清亚洲精品| 亚洲电影在线观看av| 免费大片18禁| 欧美国产日韩亚洲一区| 两个人的视频大全免费| 99国产精品一区二区蜜桃av| 国产精品电影一区二区三区| 中文字幕人成人乱码亚洲影| 精品一区二区三区视频在线 | 国产97色在线日韩免费| 麻豆成人av在线观看| 亚洲狠狠婷婷综合久久图片| 国产av一区在线观看免费| 成人国产一区最新在线观看| 国产精品久久视频播放| 国产伦精品一区二区三区四那| 欧美av亚洲av综合av国产av| 亚洲成av人片免费观看| 搞女人的毛片| 老汉色av国产亚洲站长工具| 欧美日韩亚洲国产一区二区在线观看| 黄色 视频免费看| 成人午夜高清在线视频| 欧美最黄视频在线播放免费| 亚洲国产看品久久| 国产成人av教育| 国产精品亚洲一级av第二区| 国内毛片毛片毛片毛片毛片| 男女那种视频在线观看| 午夜a级毛片| 久久草成人影院| 久久久久亚洲av毛片大全| 久久这里只有精品19| 亚洲av成人不卡在线观看播放网| 欧美成狂野欧美在线观看| 又紧又爽又黄一区二区| 国产黄色小视频在线观看| 天堂网av新在线| 国产免费av片在线观看野外av| 久久草成人影院| 男女午夜视频在线观看| 欧美日韩福利视频一区二区| 一区二区三区激情视频| 天天躁日日操中文字幕| 老司机在亚洲福利影院| 男女之事视频高清在线观看| 成人精品一区二区免费| АⅤ资源中文在线天堂| 亚洲第一电影网av| 国产精品香港三级国产av潘金莲| 欧美丝袜亚洲另类 | 美女黄网站色视频| 亚洲第一欧美日韩一区二区三区| 精品一区二区三区视频在线观看免费| 久久久久久久精品吃奶| 中文资源天堂在线| 亚洲精品美女久久av网站| 日本在线视频免费播放| 久久久国产欧美日韩av| 在线观看日韩欧美| 成人特级黄色片久久久久久久| 成人午夜高清在线视频| 亚洲国产欧美网| 十八禁网站免费在线| 日韩欧美国产在线观看| 色噜噜av男人的天堂激情| 99久久精品国产亚洲精品| av在线蜜桃| 亚洲五月天丁香| 搞女人的毛片| 噜噜噜噜噜久久久久久91| 亚洲在线自拍视频| 国产蜜桃级精品一区二区三区| 狂野欧美激情性xxxx| 波多野结衣高清无吗| 麻豆成人午夜福利视频| 亚洲专区国产一区二区| 日本免费a在线| 小蜜桃在线观看免费完整版高清| 1024香蕉在线观看| 一卡2卡三卡四卡精品乱码亚洲| 午夜免费成人在线视频| 一本一本综合久久| 在线a可以看的网站| 亚洲aⅴ乱码一区二区在线播放| 99国产精品99久久久久| 国产高清视频在线观看网站| 国产精品永久免费网站| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲人成伊人成综合网2020| 99久久久亚洲精品蜜臀av| 男插女下体视频免费在线播放| 老司机在亚洲福利影院| 免费观看人在逋| 在线观看舔阴道视频| 亚洲中文av在线| 一级毛片女人18水好多| 九九久久精品国产亚洲av麻豆 | 一本一本综合久久| 搡老熟女国产l中国老女人| 久久精品国产亚洲av香蕉五月| 精品久久久久久成人av| 欧美av亚洲av综合av国产av| 久久久久性生活片| 男人舔女人的私密视频| 久久久国产成人免费| 亚洲熟女毛片儿| 天天躁日日操中文字幕| 激情在线观看视频在线高清| 18禁美女被吸乳视频| 欧美乱码精品一区二区三区| 九九热线精品视视频播放| 丁香六月欧美| 亚洲精品国产精品久久久不卡| 亚洲av成人精品一区久久| 国产精品,欧美在线| 美女扒开内裤让男人捅视频| 国产1区2区3区精品| 麻豆国产97在线/欧美| 窝窝影院91人妻| 波多野结衣高清作品| 淫妇啪啪啪对白视频| 一二三四社区在线视频社区8| 小蜜桃在线观看免费完整版高清| 美女 人体艺术 gogo| 国产精品一区二区免费欧美| 午夜a级毛片| 日韩欧美在线二视频| 久久精品aⅴ一区二区三区四区| 色综合欧美亚洲国产小说| 国产一区在线观看成人免费| 婷婷丁香在线五月| 男人舔女人下体高潮全视频| 亚洲avbb在线观看| 亚洲美女视频黄频| 亚洲国产精品久久男人天堂| 久久精品国产亚洲av香蕉五月| 午夜成年电影在线免费观看| 亚洲国产看品久久| 欧美另类亚洲清纯唯美| 国产高清激情床上av| 老司机福利观看|