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

    基于降水量的黑龍江省玉米種植適宜性研究

    2017-03-16 08:14:42孫彥坤高松陽曹義娜許佳琦王振華
    關(guān)鍵詞:距平百分率降雨量

    孫彥坤,高松陽,曹義娜,許佳琦,王振華

    (東北農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院,哈爾濱 150030)

    基于降水量的黑龍江省玉米種植適宜性研究

    孫彥坤,高松陽,曹義娜,許佳琦,王振華

    (東北農(nóng)業(yè)大學(xué)資源與環(huán)境學(xué)院,哈爾濱 150030)

    利用黑龍江省63個氣象站點(diǎn)1965~2014年降水資料,運(yùn)用降水距平百分率、Mann-Kendall突變檢驗(yàn)和標(biāo)準(zhǔn)化降水指數(shù)等研究方法,分析黑龍江省近50年來不同區(qū)域間降水變化趨勢和特點(diǎn),結(jié)合玉米生長需水特征判斷玉米易缺水時(shí)段。結(jié)果表明,近50年全省玉米生長季降水量平均值約457.7 mm,可滿足玉米生長發(fā)育對水分基本需求;黑龍江省各區(qū)域中,中南部地區(qū)、三江平原區(qū)和北部地區(qū)生長季降水量高于全省平均值,西部地區(qū)低于全省平均值,特別是泰來和肇源,生長季降水量平均值低于375 mm,必須依靠灌溉維持玉米正常生長;黑龍江省中南部地區(qū)和北部地區(qū)生長季降水量年際間變化較小,三江平原地區(qū)和西部地區(qū)則波動較大,易發(fā)生水旱災(zāi)害;黑龍江省近50年生長季降水量主要發(fā)生2次突變,分別在20世紀(jì)80年代初期和90年代中后期;黑龍江省西部地區(qū)玉米生長期基本均處于缺水狀態(tài),需灌溉補(bǔ)充水分,其他地區(qū)缺水期主要在6月和7月。研究可為黑龍江省降水研究和玉米種植區(qū)水資源配置等提供參考。

    黑龍江?。唤邓?;突變;玉米

    水分是限制玉米氣候生產(chǎn)潛力主要因素,玉米生育期耗水量大,且對水分脅迫反應(yīng)敏感,干旱對玉米生長發(fā)育、形態(tài)特征、生理生化及產(chǎn)量品質(zhì)均有重要影響[1]。研究表明,春玉米全生育期需水量約為500~600 mm;如低于375 mm,滿足底墑條件下播種,籽粒和秸稈產(chǎn)量顯著減少[2-3]。黑龍江省是重要玉米生產(chǎn)基地,但降水時(shí)空分布不均,氣象災(zāi)害多發(fā),隨生長季平均氣溫升高,玉米光溫生產(chǎn)潛力總體呈上升趨勢;但由于降水量減少,玉米氣候生產(chǎn)潛力總體呈下降趨勢[4]。黑龍江省玉米生產(chǎn)以雨養(yǎng)為主,降水量和時(shí)空分布決定干旱發(fā)生程度。研究表明,黑龍江省東部地區(qū)和中部地區(qū)降水量較充足,西部地區(qū)干旱少雨,常出現(xiàn)旱澇災(zāi)害[5]。黑龍江省玉米種植面積不斷擴(kuò)大,但農(nóng)業(yè)排灌設(shè)施建設(shè)相對滯后,亟待準(zhǔn)確、定量評估黑龍江省降水特點(diǎn),為該地區(qū)農(nóng)業(yè)生產(chǎn)及水旱災(zāi)害防御提供依據(jù)。本文采用黑龍江省63個氣象站點(diǎn)50年逐日降水量數(shù)據(jù)資料,研究分析黑龍江省近50年降水分布和變化情況,為黑龍江省玉米生產(chǎn)及旱澇災(zāi)害防御提供參考。

    1 材料與方法

    1.1 數(shù)據(jù)處理

    為便于統(tǒng)計(jì)與比較,根據(jù)黑龍江省玉米分布區(qū)域范圍、地理位置和種植制度特點(diǎn),結(jié)合各產(chǎn)區(qū)農(nóng)業(yè)自然資源狀況,將各氣象站點(diǎn)分4個區(qū)域分析:①中南部玉米產(chǎn)區(qū)(以下簡稱“中南部玉米區(qū)”),包括哈爾濱市所轄市縣(除依蘭縣)、牡丹江市所轄市縣、綏化市東部肇東市、慶安等共20個站點(diǎn);②松嫩平原西部干旱半干旱玉米產(chǎn)區(qū)(以下簡稱“西部玉米區(qū)”),包括齊齊哈爾市所轄市縣、大慶市所轄市縣、綏化市西部安達(dá)縣、明水縣、望奎縣、青岡縣等共17個站點(diǎn);③三江平原溫和半濕潤玉米產(chǎn)區(qū)(以下簡稱“三江平原玉米區(qū)”),包括佳木斯市所轄市縣、鶴崗市所轄市縣、雙鴨山市所轄市縣、雞西市所轄市縣、依蘭縣等共15個站點(diǎn)[6];④黑龍江省北部地區(qū),包括黑河、伊春及大興安嶺地區(qū)等共11個站點(diǎn)。分區(qū)情況如圖1所示。黑龍江省北部地區(qū)非玉米主產(chǎn)區(qū),但隨全球變暖,玉米種植界限北移,故本文引入黑龍江省北部地區(qū)部分市縣資料一并分析。

    圖1 黑龍江省玉米產(chǎn)區(qū)分區(qū)Fig.1Partition of the maize growing areas in Heilongjiang Province

    本文以黑龍江省63個氣象站點(diǎn)1965~2014年逐日降水量數(shù)據(jù)為研究對象,資料來源于黑龍江省氣象臺,利用Excel軟件對數(shù)據(jù)作標(biāo)準(zhǔn)化處理及統(tǒng)計(jì)分析,計(jì)算各市縣月降水量、年降水量及生長季降水量,計(jì)算各區(qū)域月降水量、年降水量及生長季降水量平均值。黑龍江省春玉米播種期一般為4月下旬至5月中旬,由于在雨養(yǎng)環(huán)境中播種前段降水影響播種期和苗期土壤水分,受全球氣候變暖等因素影響,黑龍江省農(nóng)業(yè)界限溫度起始日期提前[7],所以文中將4~9月作為黑龍江省玉米生長季分析。

    1.2 方法

    1.2.1 降水量距平百分率

    降水量距平百分率是表征某時(shí)段降水量較常年值偏差程度指標(biāo)之一,直觀反映降水異常引起氣象干旱[8]。降水量距平百分率因其計(jì)算簡單、物理意義明確,被廣泛應(yīng)用于氣象干旱監(jiān)測和評估[9],其表達(dá)式為:

    1.2.2 Mann-Kendall突變檢驗(yàn)

    氣候突變是指兩個穩(wěn)定氣候階段之間統(tǒng)計(jì)特征量發(fā)生顯著變化現(xiàn)象。Mann-Kendall法是一種非參數(shù)統(tǒng)計(jì)檢驗(yàn)方法[10],不需樣本遵從一定分布,也不受少數(shù)異常值影響[11]。本文應(yīng)用Mann-Kendall法突變檢驗(yàn),通過構(gòu)造秩序列:

    秩序列Sk是第i時(shí)刻數(shù)值大于j時(shí)刻數(shù)值個數(shù)累計(jì)數(shù)。

    在時(shí)間序列隨機(jī)獨(dú)立假定下,定義統(tǒng)計(jì)量

    式中,UF1=0,E(Sk)、var(Sk)是累計(jì)數(shù)Sk均值和方差,在x1,x2,…,xn相互獨(dú)立,且有相同連續(xù)分布時(shí),由下式算出:

    UFi為標(biāo)準(zhǔn)正態(tài)分布,是按時(shí)間序列x順序x1,x2,…,xn計(jì)算出統(tǒng)計(jì)量序列,給定顯著性水平α,若|UFi|>UFα,則表明序列存在明顯趨勢變化。

    按照時(shí)間序列x逆序xn,xn-1,…,x1,重復(fù)上述過程得到逆序統(tǒng)計(jì)量值

    UB1=-UFi(k=n,n-1,…,1)

    其中,UB1=0。

    繪制UFk和UBk曲線圖,若UFk>0,則表明序列呈上升趨勢,UFk<0則呈下降趨勢。當(dāng)其超過臨界線時(shí),表明上升或下降趨勢顯著。如果UFk和UBk兩條曲線出現(xiàn)交點(diǎn),交點(diǎn)在臨界線之間,交點(diǎn)對應(yīng)時(shí)刻即突變發(fā)生時(shí)間。

    1.2.3 標(biāo)準(zhǔn)化降水指數(shù)

    標(biāo)準(zhǔn)化降水指數(shù)(SPI)是通過計(jì)算某時(shí)段內(nèi)降水量Г分布概率,正態(tài)標(biāo)準(zhǔn)化處理,最終用標(biāo)準(zhǔn)化降水累計(jì)頻率分布劃分干旱等級研究方法;標(biāo)準(zhǔn)化降水指數(shù)在不同地區(qū)和時(shí)段均具有良好計(jì)算穩(wěn)定性,可有效反映旱澇狀況[12],適于月以上尺度氣候狀況的干旱監(jiān)測與評估[13]。

    SPI計(jì)算公式為:

    G(x)為與Г函數(shù)相關(guān)降水分布概率;x為降水量樣本值;S為概率密度正負(fù)系數(shù)。

    當(dāng)G(x)>0.5時(shí),S=1;當(dāng)G(x)≤0.5時(shí),S=-1。G (x)可由以下Г分布函數(shù)概率密度積分公式得到

    式中,β和γ分別為Г分布函數(shù)尺度和形狀參數(shù);c0、c1、c2和d1、d2、d3分別為Г分布函數(shù)轉(zhuǎn)換為累計(jì)頻率簡化近似求解公式計(jì)算參數(shù),取值如下:c0=2.515517,c1=0.802853,c2=0.010328,d1=1.432788,d2=0.189269,d3=0.001308。

    SPI計(jì)算涉及較復(fù)雜Г分布函數(shù)[14-15]。

    2 結(jié)果與分析

    2.1 黑龍江省全省及各區(qū)域降水基本情況

    黑龍江省近50年(1965~2014年)全省年平均降水量約519.9 mm,最大值為2013年690.9 mm,最小值為2001年385.5 mm;全省生長季(4~9月)降雨量平均值為457.7 mm,低于玉米生長季正常需水量,其中最大值為603.0 mm(1994年),最小值為333.7 mm(2001年)。近50年全省玉米生長季平均降水量低于500 mm年份多達(dá)36年,其中7年少于375 mm;大于600 mm僅1年(1994年)。

    2.1.1 中南部玉米區(qū)降水基本情況

    神華大準(zhǔn)鐵路田家石畔特大橋位于大準(zhǔn)線點(diǎn)岱溝至龍王渠區(qū)間,中心里程K3+332,橋梁全長1254m,采用2101-Ⅱ-26《人行道支架構(gòu)造圖》制作的寬度為1.55m角鋼支架人行道,在多年的運(yùn)營過程中,由于列車震動、銹蝕等原因,造成橋上部分拖架螺栓折斷。從人行道托架螺栓的折斷形式、托架支撐結(jié)構(gòu)的失效形態(tài)以及托架螺栓銹蝕、外觀質(zhì)量較差的情況可作出如下推斷:托架承擔(dān)了較大的人行道靜荷載,而螺栓銹蝕導(dǎo)致其自身強(qiáng)度和承載能力不斷降低,造成托架螺栓折斷的主要原因是在長期交變動應(yīng)力作用下產(chǎn)生的疲勞破壞。

    中南部玉米區(qū)玉米種植面積大、產(chǎn)量高,是全省玉米高產(chǎn)區(qū)。該區(qū)域地勢平坦,土壤多為黑鈣土,土質(zhì)肥沃,是全省溫度、土壤條件和生產(chǎn)水平最優(yōu)越地區(qū)。中南部玉米種植區(qū)年降雨量平均值約547.3 mm,最高值793.6 mm(1994年),最低值409.4 mm(1999年)。其中年降水量最大值1994年出現(xiàn)在尚志市,為1 024.6 mm,同年木蘭縣年降水量也超過1 000 mm,達(dá)1 017 mm;全區(qū)最小年降水量1989年出現(xiàn)在雙城區(qū),當(dāng)年年降水量為246.9 mm。全區(qū)生長季降水量平均值約為475.8 mm,最大值為704.5 mm(1994年),最小值為333.8 mm(1999年)。其中生長季最大降水量出現(xiàn)在尚志市,為934.7 mm(1994年),同年木蘭縣和賓縣生長季降雨量也超過900 mm,分別為919.4 mm和909.5 mm;最小值在雙城1989年,為215.1 mm。全區(qū)平均生長季降雨量低于500 mm年份為31年,多于500 mm年份為19年,其中僅1年超過600 mm(1994年)。本區(qū)域內(nèi)近50年生長季降水量平均值不少于500 mm市縣有5個,分別為:五常519.4 mm、延壽504.7 mm、木蘭527.0 mm、巴彥518.4 mm、尚志547.2 mm;區(qū)域內(nèi)無生長季降水量平均值低于375 mm市縣。

    2.1.2 松嫩平原西部玉米區(qū)降水基本情況

    松嫩平原西部玉米種植區(qū)是黑龍江省玉米主產(chǎn)區(qū)之一,該區(qū)域南部降水量少,土壤多為風(fēng)沙土,熱量資源充足;中北部地區(qū)土壤多為黑土,土質(zhì)肥沃。西部玉米區(qū)年降雨量平均值約459.9 mm,最高值為692.4 mm(1998年),最低值為2001年407.4 mm。其中年降水量最大值為1998年甘南縣981.3 mm,最小值在2001年泰來縣,僅189.2 mm。全區(qū)生長季降水量平均值約為421.6 mm,最大值為631.8 mm(1998年),最小值為275.4 mm(2004年)。其中生長季最大降雨量出現(xiàn)在甘南縣,為926.7 mm(1998年),最小值在2004年龍江縣,僅為158.7 mm,且西區(qū)各市縣生長季降水量最低值均低于300 mm。全區(qū)平均生長季降雨量低于500 mm年份多達(dá)39年,其中低于375 mm年份多達(dá)16年;超過600 mm年份為1年(1998年)。本區(qū)域內(nèi)無近50年生長季降雨量平均值達(dá)500 mm市縣,其中泰來356.9 mm和肇源374.1 mm生長季降雨量平均值低于375 mm。

    2.1.3 三江平原玉米區(qū)降水基本情況

    三江平原玉米種植區(qū)土壤多為草甸土,土質(zhì)肥沃;由于該區(qū)域地勢較低,易受澇災(zāi)。三江平原玉米區(qū)年降雨量平均值約539.7 mm,最大值為796.8 mm(1994年),最小值為351.7 mm(1975年)。其中年降水量最大值為1994年鶴崗市965 mm,且全區(qū)各市縣年降水量最大值均超過700 mm;最小值1975年出現(xiàn)在樺南縣,為276.2 mm。全區(qū)生長季降水量平均值約464.0 mm,最大值為702.8 mm(1994年),最小值為313.2 mm(1977年)。其中生長季最大降雨量出現(xiàn)在鶴崗市,為859.4 mm(1994年),最小值2003年出現(xiàn)在寶清,為226.7 mm,且全區(qū)各市縣最低生長季降水量除虎林市為301.1 mm外均少于300 mm。本區(qū)域內(nèi)生長季降雨量低于500 mm年份有30年,降雨量500~600 mm有18年,有2年高于600 mm。本區(qū)域內(nèi)鶴崗(553.2 mm)和蘿北(505.8 mm)近50年生長季降雨量平均值不少于500 mm,且其余13個市縣生長季降雨量平均值均不少于400 mm。

    2.1.4 黑龍江省北部地區(qū)降水基本情況

    黑龍江省北部地區(qū)年降雨量平均值約535.7 mm,最大值為715.3 mm(2013年),最小值399.2 mm(2007年)。本區(qū)域內(nèi)降水量最大值為1985年伊春市995.5 mm,最小值同樣出現(xiàn)在伊春市,為1979年274.3 mm。本區(qū)域生長季降水量平均值約為471.8 mm,最大值為625.4 mm(2013年),最小值為327.6 mm(2007年)。其中生長季最大降雨量1985年出現(xiàn)在伊春市,為875.9 mm;最小降雨量出現(xiàn)在漠河,為203.7 mm(2002年)。本區(qū)域內(nèi)生長季降雨量少于500 mm年份為35年,少于375 mm年份為2年,多于600 mm年份為2年。區(qū)域內(nèi)伊春546.9 mm、鐵力546.4 mm、五營532.7 mm近50年生長季降雨量平均值超過500 mm,無低于375 mm市縣。

    2.2 全省及各區(qū)域近50年降水變化情況

    2.2.1 全省及各區(qū)域年降水量變化情況

    近50年黑龍江省全省及各區(qū)域年降水量變化形式見圖2,由圖2可見,全省年降雨量均值和各分區(qū)年降雨量值近50年中大多數(shù)年份在400~600 mm,全省年降雨量1965~1975年在500 mm線附近浮動,1976~1980年降雨量減少,2009~2014年又呈波動上升趨勢,2013年達(dá)最大值。省內(nèi)各區(qū)域年降水量變化趨勢與全省年降水量變化趨勢一致,中南部玉米區(qū)年降水量曲線變化情況與全省年降水量曲線接近,在1985~1989年有略高于全省線時(shí)段;西部玉米區(qū)在近50年中大部分年降水量均低于全省平均值,僅1998年出現(xiàn)高于全省平均值高值點(diǎn)(692.4 mm),有15年年降水量低于400mm;三江平原玉米區(qū)年降水量波動范圍較大,兩個年降水量高值點(diǎn)均出現(xiàn)在三江平原玉米區(qū)(1981和1994年),年降水量接近800 mm,而該區(qū)域年降水量低值則低至351.7 mm(1975年);北部地區(qū)年降水量總體上略高于全省年降水量,僅在2003年達(dá)到高值(660.6 mm),高于全省及其他各區(qū)域年降水量值。

    2.2.2 全省及各區(qū)域生長季降水量變化情況

    近50年黑龍江省全省及各分區(qū)玉米生長季降水量變化形式如圖3所示,各分區(qū)玉米生長季降水量變化趨勢與全省范圍玉米生長季降水量總體上保持一致,全省玉米生長季降水量在1965~1980年處在降水量較低時(shí)期,低于500 mm;1981~1998年全省玉米生長季降水量值在500 mm浮動;在1999~2011年又經(jīng)歷降水量低于500 mm相對少雨期,從2012年開始又呈上升趨勢。最高值在三江平原區(qū)(702.8mm),最低值在西部玉米區(qū)(275.4 mm)。近50年生長季降水量波動變化幅度?。蝗皆衩讌^(qū)在1975~1979年生長季降水量明顯下降。

    圖2 1965~2014年黑龍江省全省及各分區(qū)年降水量變化特征Fig.2 Annual precipitation of Heilongjiang Province and each region in 1965-2014

    圖3 1965~2014年黑龍江省全省及各分區(qū)4~9月降水量變化特征Fig.3Precipitation in April to September of Heilongjiang Province and each region in 1965-2014

    2.3各區(qū)域生長季降水距平百分率變化情況

    如圖4所示,中南部玉米區(qū)生長季降水量平穩(wěn),降水距平百分率波動不大,僅1年距平百分率超40%(1994年),有10年距平百分率絕對值略高于20%,其余均低于20%??梢姡心喜坑衩讌^(qū)在1974~1980年處于生長季降水量下降時(shí)期,1983~1988年處于生長季降水量上升時(shí)期,1994年降水量達(dá)高值后距平百分率不斷波動,且負(fù)值多于正值,說明降水量下降年份多于上升年份。西部玉米區(qū)降水距平百分率波動幅度較大,該區(qū)域前15年(1965~1979年)間僅4年降水距平百分率為正值,說明該時(shí)期西部玉米區(qū)處于生長季降水量下降時(shí)期;1983~1988年,距平百分率均為正,降水量處于上升階段;在1998年達(dá)到高值后,距平百分率連續(xù)4年為負(fù)值,降水量再次下降,此后距平百分率不斷波動,2012~2014年再次出現(xiàn)降水量上升期。三江平原玉米區(qū)總體降水距平百分率波動不大,但該區(qū)域高值和低值絕對值很大,高值約為51.5%、低值約為-32.5%,說明該區(qū)域在氣候異常時(shí)易出現(xiàn)降水,導(dǎo)致旱澇災(zāi)害發(fā)生,該區(qū)域降水不斷波動,僅1975~1980年降水量持續(xù)下降。北部地區(qū)近50年生長季降水量變化較為平穩(wěn),距平百分率無絕對值超過40%極值,2013和2007年分別超過30%高值和低值(其余大多年份降水距平百分率絕對值小于20%;該區(qū)域1973~1980年降水距平百分率為負(fù)值,生長季降水量處于下降期,此后降水距平百分率不斷波動變化,多為負(fù)值。

    圖4 各區(qū)域生長季降水距平百分率變化趨勢Fig.4Variation trends of each region's growth season precipitation anomaly percentage

    2.4 各區(qū)域多年生長季降水量突變分析

    利用Mann-Kendall突變檢驗(yàn)法研究黑龍江省各玉米種植區(qū)多年降水量突變性,結(jié)果表明,除西部玉米區(qū)個別年份有顯著上升趨勢外,四個區(qū)域近50年生長季降水量未出現(xiàn)顯著變化。各區(qū)域Mann-Kendall突變檢驗(yàn)結(jié)果如圖5所示。

    由于M-K方法不適用于多點(diǎn)突變檢驗(yàn),故采用滑動T檢驗(yàn)法確定突變點(diǎn)。結(jié)果結(jié)合Mann-Kendall檢驗(yàn)結(jié)果,可判斷各區(qū)域生長季降水量突變點(diǎn)分別為:中南部玉米區(qū)突變點(diǎn)在1984和1995年附近,西部玉米區(qū)突變點(diǎn)在1981和1998年附近,三江平原玉米區(qū)突變點(diǎn)在1981和1994年附近,北部地區(qū)突變點(diǎn)在1981年和1997年附近。黑龍江省近50年玉米生長季降水量突變主要發(fā)生在20世紀(jì)80年代初期和90年代中后期,與靳春香[16]等研究結(jié)果一致。

    圖5 各區(qū)域Mann-Kendall突變檢驗(yàn)Fig.5Results of Mann-Kendall tests in each region

    表1 各區(qū)域突變點(diǎn)Table 1 Abrupt changes of precipitation in each region

    2.5 各區(qū)域年尺度標(biāo)準(zhǔn)化降水指數(shù)變化情況

    黑龍江省各區(qū)域以年為尺度標(biāo)準(zhǔn)化降水指數(shù)(SPI12)變化情況如圖6所示??傮w上,各區(qū)域SPI值與其年降水量變化趨勢一致。-0.5<SPI<0.5為正常;-1.0<SPI≤-0.5、-1.5<SPI≤-1.0、 -2.0<SPI≤-1.5和SPI≤-2.0分別代表輕度干旱、中度干旱、重度干旱和特別干旱;0.5≤SPI<1.0、1.0≤SPI<1.5、1.5≤SPI<2.0和SPI≥2.0分別代表輕度濕潤、中度濕潤、重度濕潤和特別濕潤。由圖6可知,4個區(qū)域各出現(xiàn)1次特別干旱和特別濕潤情況,三江平原玉米區(qū)特別干旱年出現(xiàn)在1975年,中南部玉米區(qū)、西部玉米區(qū)和北部地區(qū)特別干旱年則比較接近,分別出現(xiàn)在1999、2001和2007年;另一方面,中南部玉米區(qū)和三江平原玉米區(qū)特別濕潤年份均出現(xiàn)在1994年,西部玉米區(qū)為1998年,北部地區(qū)則出現(xiàn)在2013年。上述情況與各區(qū)域年降水量變化情況吻合,說明標(biāo)準(zhǔn)化降水指數(shù)體現(xiàn)不同地區(qū)年際間干濕程度變化。特別干旱年和特別濕潤年除三江平原區(qū)特別干旱年出現(xiàn)在20世紀(jì)70年代外,其余均出現(xiàn)在20世紀(jì)90年代后,表明20世紀(jì)90年代中期以來黑龍江省極端氣候頻繁。

    由于SPI計(jì)算特性,即假定所有地點(diǎn)發(fā)生極端旱澇時(shí)間概率相同,標(biāo)準(zhǔn)化降水指數(shù)無法標(biāo)識旱澇頻發(fā)地區(qū);但由于其本身注重標(biāo)準(zhǔn)化正態(tài)分布,可體現(xiàn)旱澇年際變化。由圖5結(jié)合SPI旱澇等級標(biāo)準(zhǔn)可見,黑龍江省在20世紀(jì)70年代處于比較干旱時(shí)期,80年代呈濕潤趨勢,90年代末至21世紀(jì)初處較干旱階段,2012年始再次進(jìn)入相對濕潤時(shí)期。

    圖6 各區(qū)域年度標(biāo)準(zhǔn)化降水指數(shù)變化情況Fig.6SPI12of each region

    2.6 各區(qū)域生長季降水量月均變化特征

    黑龍江省各區(qū)域多年生長季降水量逐月平均值曲線見圖7。由圖7可知,4個區(qū)域月均降水量變化趨勢基本一致,均為4~6月逐漸增加,在7月時(shí)達(dá)最大值,8和9月逐漸下降。西部玉米區(qū)除7月外其他各月降水量均少于另外3個區(qū)域,特別是5、8和9月;中南部玉米區(qū)和北部地區(qū)在生長季中各月降水量較接近;三江平原玉米區(qū)7月降水量明顯少于另外3個區(qū)域,其余5個月降水量與中南部玉米區(qū)和北部地區(qū)接近。

    玉米在不同生長期需水量不同,生育期干旱導(dǎo)致玉米減產(chǎn),其中抽雄吐絲期水分脅迫造成減產(chǎn)最嚴(yán)重[17]。黑龍江省西部半干旱區(qū)玉米4~9月需水量依次為:25~40.8、45~51.5、115~136.5、185~215.8、130~158.7和65~72 mm,其他地區(qū)玉米生長季需水量接近或略低于該參考值[18-19]。

    圖7 黑龍江省各區(qū)域生長季降水量逐月平均值Fig.7Average value of monthly precipitation in growth season of each region in Heilongjiang Province

    由圖5可知,黑龍江省各區(qū)域6和7月降水量明顯低于玉米當(dāng)月需水量,說明各區(qū)域在正常年份兩月降水無法滿足玉米正常生長發(fā)育需要;中南部玉米區(qū)、三江平原玉米區(qū)和北部地區(qū)在4、5、8和9月降水量均可達(dá)到或接近需水量參考值,說明這3個地區(qū)4、5、8、9月降水量可滿足玉米正常生長需要;西部玉米區(qū)在整個生長季降水量均低于需水量參考值,說明西部玉米區(qū)在正常年份玉米生長整體上處于缺水狀態(tài)。各區(qū)域生長季每月玉米灌溉水量參考值見表2。由于黑龍江省各地區(qū)玉米播種時(shí)間不同,統(tǒng)一自5月始計(jì)算各區(qū)域玉米生長季灌溉需水量。

    表2 黑龍江省各區(qū)域玉米生長季灌溉需水量Table 2Irrigation water requirement of each region in Heilongjiang Province in maize growth season

    3 討論

    作物產(chǎn)量與氣候條件相關(guān),玉米生長發(fā)育受水分、溫度、光照、營養(yǎng)物質(zhì)等自然條件變化和干旱和洪澇等農(nóng)業(yè)氣象災(zāi)害影響。本文重點(diǎn)分析降水因素對黑龍江省玉米生產(chǎn)影響,但未考慮玉米品種差異、研究區(qū)地理特點(diǎn)和人為因素等影響,有待后續(xù)研究深入探討。

    前人研究通常未將黑龍江省北部地區(qū)作為玉米種植區(qū)影響,本文重點(diǎn)討論降水對玉米種植適宜性影響,考慮全球氣候變暖、黑龍江省玉米種植界限北擴(kuò)等原因,將北部地區(qū)納入研究范圍,全面分析黑龍江省降水資源分布情況。

    4 結(jié)論

    a.近50年黑龍江省范圍年降水量未出現(xiàn)大幅波動,總體上在400~600 mm變動;玉米生長季降水量在400~500 mm波動,平均值約457.7 mm。在降水量正常年份,可滿足玉米生長基本需要,但無法滿足玉米高產(chǎn)穩(wěn)產(chǎn)所需水分,需灌溉補(bǔ)充。

    b.從各分區(qū)看,中南部玉米區(qū)、三江平原玉米區(qū)和黑龍江省北部地區(qū)生長季降水量高于全省平均值,分別為475.8、464.0和471.8 mm。其中,中南部玉米區(qū)和北部地區(qū)生長季降水量變化趨勢較平穩(wěn);而三江平原地區(qū)年際間波動幅度較大,最大值達(dá)702.8 mm,最小值為313.2 mm,說明該地區(qū)降水量年際分布不均,出現(xiàn)降水異常時(shí)易發(fā)生水旱災(zāi)害。松嫩平原西部玉米區(qū)是4個區(qū)域中生長季降水量最少地區(qū),生長季降水量平均值約421.6 mm,低于全省平均值;特別是泰來和肇源生長季降水量平均值低于375 mm,無法滿足玉米生長對水分基本需求,須灌溉維持玉米正常生長。

    c.根據(jù)降水距平百分率和Mann-Kendall突變分析結(jié)果,西部玉米區(qū)生長季降水量波動強(qiáng)烈,其次為三江平原地區(qū),中南部玉米區(qū)和黑龍江省北部地區(qū)生長季降水量年際間分布均勻。Mann-Kendall突變分析表明,4個區(qū)域均存在2次突變,發(fā)生時(shí)間接近,分別在20世紀(jì)80年代初期和90年代中后期,第一次突變由降轉(zhuǎn)升,第二次突變由升轉(zhuǎn)降。標(biāo)準(zhǔn)化降水指數(shù)變化情況則顯示20世紀(jì)90年代以來黑龍江省極端氣候頻繁出現(xiàn)。

    d.各區(qū)域生長季逐月降水量分析顯示,黑龍江省內(nèi)4個區(qū)域玉米生長季月均降水量變化趨勢一致,均從4月始逐漸增加,至7月達(dá)最大值,8和9月下降。與玉米生長季需水量參考值比較說明,中南部玉米區(qū)、三江平原玉米區(qū)和北部地區(qū)玉米缺水發(fā)生時(shí)段主要集中在6和7月,5、8、9月自然降水量可滿足玉米正常生長需要;西部玉米區(qū)整個生長季處于缺水狀態(tài),須灌溉維持玉米正常生長。

    [1]楊重一,龐士力,孫彥坤.黑龍江省作物氣候生產(chǎn)潛力估算[J].東北農(nóng)業(yè)大學(xué)學(xué)報(bào),2010,41(3):75-78.

    [2]董秋婷,李茂松,劉江,等.近50年東北地區(qū)春玉米干旱時(shí)空演變特征[J].自然災(zāi)害學(xué)報(bào),2011,20(4):52-59.

    [3]劉鈺,汪林,倪廣恒,等.中國主要作物灌溉需水量空間分布特征[J].農(nóng)業(yè)工程學(xué)報(bào),2009,25(12):6-12.

    [4]孫彥坤,田寶星,高見,等.氣候變化對黑龍江省黑土區(qū)玉米氣候生產(chǎn)潛力影響[J].東北農(nóng)業(yè)大學(xué)學(xué)報(bào),2013,44(11):44-49.

    [5]楊若子,周廣勝.東北三省玉米主要農(nóng)業(yè)氣象災(zāi)害綜合危險(xiǎn)性評估[J].氣象學(xué)報(bào),2015,73(6):1141-1153.

    [6]張研,錢春榮.黑龍江省玉米單產(chǎn)水平現(xiàn)狀與發(fā)展策略[J].東北農(nóng)業(yè)大學(xué)學(xué)報(bào),2010,41(7):155-160.

    [7]康蕾,張紅旗.我國五大糧食主產(chǎn)區(qū)農(nóng)業(yè)干旱態(tài)勢綜合研究[J].中國生態(tài)農(nóng)業(yè)學(xué)報(bào),2014,22(8):928-937.

    [8]馬建勇,許吟隆,潘婕.東北地區(qū)農(nóng)業(yè)氣象災(zāi)害趨勢變化及其對糧食產(chǎn)量影響[J].中國農(nóng)業(yè)氣象,2012,33(2):283-288.

    [9]云雅如,方修琦,王媛,等.黑龍江省過去20年糧食作物種植格局變化及其氣候背景[J].自然資源學(xué)報(bào),2005,20(5):697-705.

    [10]魏湜,金益,張樹權(quán),等.黑龍江玉米生理生態(tài)與栽培[M].北京:中國農(nóng)業(yè)出版社,2014.

    [11]孫彥坤,高璐迪,雷德福,等.黑龍江省農(nóng)業(yè)界限溫度時(shí)空變化特征分析[J].東北農(nóng)業(yè)大學(xué)學(xué)報(bào),2016,47(5):44-53.

    [12]國家氣候中心,中國氣象科學(xué)研究院,國家氣象中心.氣象干旱等級(GB/T20481-2006)[S].北京,中國標(biāo)準(zhǔn)出版社,2006.

    [13]高曉容.東北地區(qū)玉米主要?dú)庀鬄?zāi)害風(fēng)險(xiǎn)評估研究[D].南京:南京信息工程大學(xué),2012.

    [14]袁文平,周廣勝.標(biāo)準(zhǔn)化降水指標(biāo)與Z指數(shù)在我國應(yīng)用對比分析[J].植物生態(tài)學(xué)報(bào),2004,28(4):523-529.

    [15]黃晚華,楊曉光,李茂松,等.基于標(biāo)準(zhǔn)化降水指數(shù)的中國南方季節(jié)性干旱近58a演變特征[J].農(nóng)業(yè)工程學(xué)報(bào),2010,26(7): 50-59.

    [16]靳春香,王秀茹,王希,等.黑龍江省近50年降水變化趨勢及空間分布特征[J].中國水土保持科學(xué),2015,13(1):76-83.

    [17]白向歷,孫世賢,楊國航,等.不同剩余時(shí)期水分脅迫對玉米產(chǎn)量及生長發(fā)育影響[J].玉米科學(xué),2009,17(2):60-63.

    [18]司振江.黑龍江省西部半干旱區(qū)干旱特征及預(yù)測模型研究[D].長春:吉林大學(xué),2006.

    [19]張建平,王春乙,楊曉光,等.近26a來東北三省玉米生育期內(nèi)需水量時(shí)空分布特征[C].第26屆中國氣象學(xué)會年會農(nóng)業(yè)氣象防災(zāi)減災(zāi)與糧食安全分會場論文集,2009.

    Study on planting suitability of maize in Heilongjiang Province based

    on precipitation

    /SUN Yankun,GAO Songyang,CAO Yina,XU Jiaqi,WANG Zhenhua(School of

    Resources and Environmental Sciences,NortheastAgricultural University,Harbin 150030,China)

    This paper was based on the precipitation data of 63 meteorological stations in Heilongjiang Province,from 1965 to 2014,the variation trend and characteristics of precipitation in recent 50 years among different regions in Heilongjiang Province by using the methods of precipitation anomaly percentage,Mann-Kendall method and SPI,and estimate the short periods of water in combination with the water demand characteristics of maize.Results showed that,In recent 50 years,the average precipitation of Heilongjiang Province during maize growth season was about 457.7 mm,that could satisfy the basic needs of maize growth and development on water.In Heilongjiang Province,the growing season precipitations of the central and southern region,the Sanjiang Plain region and the northern region were higher than the average precipitation of the whole province,but the precipitation of the western region was lower,especially in Tailai and Zhaoyuan,the average growing season precipitations were lower than 375 mm,irrigation must be used to maintain normal growth of maize.The inter-annual change of growing season precipitations in thenorthern region and the central and southern region were small;but in the Sanjiang Plain region and the western region were unstable,floods and droughts were likely to occur.There were two abrupt changes of growing season precipitations in Heilongjiang Province,in early 1980s and mid and late 1990s.In the western region,almost the whole maize growth period was in lack of water,irrigation was needed to add moisture;in the other regions,the main water shortage period was in June and July.The results provided reference for the precipitation research and water resources allocation in maize planting areas in Heilongjiang Province.

    Heilongjiang Province;precipitation;abrupt change;maize

    S513

    A

    1005-9369(2017)01-0033-09

    2016-11-07

    聯(lián)合國環(huán)境規(guī)劃署綠色發(fā)展項(xiàng)目(CPR/06/209-02)

    孫彥坤(1966-),男,教授,博士,博士生導(dǎo)師,研究方向?yàn)檗r(nóng)業(yè)氣象學(xué)。E-mail:511064380@qq.com

    時(shí)間2017-1-9 15:46:12[URL]http://www.cnki.net/kcms/detail/23.1391.S.20170109.1546.020.html

    孫彥坤,高松陽,曹義娜,等.基于降水量的黑龍江省玉米種植適宜性研究[J].東北農(nóng)業(yè)大學(xué)學(xué)報(bào),2017,48(1):33-41.

    Sun Yankun,Gao Songyang,Cao Yina,et al.Study on planting suitability of maize in Heilongjiang Province based on precipitation[J].Journal of Northeast Agricultural University,2017,48(1):33-41.(in Chinese with English abstract)

    猜你喜歡
    距平百分率降雨量
    簡圖法理解標(biāo)記有絲分裂百分率法測定細(xì)胞周期時(shí)間
    颶風(fēng)Edouard(2014)暖心結(jié)構(gòu)的多資料對比分析
    降雨量與面積的關(guān)系
    近40年阿里地區(qū)云量和氣溫的年際變化
    西藏科技(2018年9期)2018-10-17 05:51:30
    甘肅省降水和冰雹天氣氣候分析
    洞庭湖區(qū)降雨特性分析
    基于距平百分率的萬源市氣象干旱灰色預(yù)測
    羅甸縣各鄉(xiāng)鎮(zhèn)實(shí)測降雨量分析及應(yīng)用研究
    美國“富人稅”最重的十個州
    海外星云 (2012年13期)2012-04-29 00:44:03
    全美“最胖”城市
    海外星云 (2012年11期)2012-04-29 00:44:03
    久久久久久人人人人人| 午夜福利乱码中文字幕| 午夜免费观看性视频| 人人妻人人爽人人添夜夜欢视频| 天天添夜夜摸| 99热国产这里只有精品6| 亚洲情色 制服丝袜| 天堂8中文在线网| av国产久精品久网站免费入址| 亚洲欧美清纯卡通| 巨乳人妻的诱惑在线观看| 天天添夜夜摸| 少妇的丰满在线观看| 亚洲国产欧美网| 欧美人与性动交α欧美软件| 免费黄频网站在线观看国产| 看免费成人av毛片| 在线观看免费午夜福利视频| 精品国产一区二区久久| 精品卡一卡二卡四卡免费| 中文字幕制服av| 国产伦理片在线播放av一区| 国产精品成人在线| av天堂久久9| 国产精品久久久久久精品电影小说| 亚洲伊人色综图| 五月开心婷婷网| 国产精品无大码| 亚洲欧洲精品一区二区精品久久久 | 亚洲精品国产一区二区精华液| 夜夜骑夜夜射夜夜干| 日日啪夜夜爽| 免费观看av网站的网址| 啦啦啦在线观看免费高清www| 国语对白做爰xxxⅹ性视频网站| 国产在线一区二区三区精| 女性生殖器流出的白浆| 51午夜福利影视在线观看| 国产97色在线日韩免费| 国产精品嫩草影院av在线观看| 91成人精品电影| 精品国产一区二区三区久久久樱花| 亚洲国产av新网站| 精品久久久精品久久久| 日日啪夜夜爽| 你懂的网址亚洲精品在线观看| 女人精品久久久久毛片| 久久精品亚洲熟妇少妇任你| 国产亚洲av片在线观看秒播厂| 人人妻人人澡人人爽人人夜夜| 少妇被粗大猛烈的视频| 91精品国产国语对白视频| 999精品在线视频| 亚洲久久久国产精品| 中文欧美无线码| 国产日韩一区二区三区精品不卡| 老汉色∧v一级毛片| 亚洲,欧美,日韩| 国产精品成人在线| 卡戴珊不雅视频在线播放| 女的被弄到高潮叫床怎么办| 五月开心婷婷网| 亚洲欧洲国产日韩| 欧美 亚洲 国产 日韩一| 日本91视频免费播放| 99精国产麻豆久久婷婷| 一边亲一边摸免费视频| 国产乱来视频区| 亚洲视频免费观看视频| 美女主播在线视频| 日韩欧美精品免费久久| 一区福利在线观看| 国产精品偷伦视频观看了| 人人澡人人妻人| 亚洲精品国产区一区二| 九草在线视频观看| 精品久久久久久电影网| 久久精品亚洲熟妇少妇任你| 久久天堂一区二区三区四区| 看非洲黑人一级黄片| 久久 成人 亚洲| 熟妇人妻不卡中文字幕| 久久精品国产亚洲av高清一级| 大片电影免费在线观看免费| 侵犯人妻中文字幕一二三四区| 国产午夜精品一二区理论片| 成人黄色视频免费在线看| 久久青草综合色| 三上悠亚av全集在线观看| 伦理电影大哥的女人| 国产麻豆69| 亚洲婷婷狠狠爱综合网| 波多野结衣av一区二区av| 中文字幕人妻丝袜制服| 免费av中文字幕在线| 亚洲欧美激情在线| 国产精品.久久久| 久久这里只有精品19| 国语对白做爰xxxⅹ性视频网站| 国产精品麻豆人妻色哟哟久久| 亚洲欧美精品自产自拍| 丰满乱子伦码专区| 国产精品一二三区在线看| 在线精品无人区一区二区三| av片东京热男人的天堂| 国产黄色视频一区二区在线观看| 久久亚洲国产成人精品v| 国产一区二区三区av在线| 精品少妇久久久久久888优播| 中文精品一卡2卡3卡4更新| 亚洲四区av| 男人添女人高潮全过程视频| 亚洲美女黄色视频免费看| 久久青草综合色| 久久精品国产a三级三级三级| 国产熟女欧美一区二区| 激情五月婷婷亚洲| 妹子高潮喷水视频| 丝袜美腿诱惑在线| 国产亚洲午夜精品一区二区久久| 一区二区三区激情视频| 亚洲成人av在线免费| 各种免费的搞黄视频| 国语对白做爰xxxⅹ性视频网站| 国产色婷婷99| 亚洲色图 男人天堂 中文字幕| 久久 成人 亚洲| 亚洲成人免费av在线播放| 中文字幕人妻丝袜制服| 亚洲情色 制服丝袜| 999久久久国产精品视频| 黑丝袜美女国产一区| 中文精品一卡2卡3卡4更新| 少妇人妻精品综合一区二区| 女性生殖器流出的白浆| 2018国产大陆天天弄谢| 国产一级毛片在线| 免费高清在线观看日韩| 欧美少妇被猛烈插入视频| 丰满饥渴人妻一区二区三| 九草在线视频观看| 国产亚洲欧美精品永久| 丝袜在线中文字幕| 国产高清不卡午夜福利| 国产男女内射视频| 美女视频免费永久观看网站| 大话2 男鬼变身卡| 日本爱情动作片www.在线观看| 亚洲精品一区蜜桃| 亚洲av国产av综合av卡| 一级毛片我不卡| h视频一区二区三区| 免费人妻精品一区二区三区视频| 一级毛片我不卡| 日韩一区二区视频免费看| 国产精品成人在线| 精品少妇内射三级| 国产成人午夜福利电影在线观看| 操美女的视频在线观看| 日韩中文字幕欧美一区二区 | 中文字幕人妻熟女乱码| 国产 精品1| 精品少妇黑人巨大在线播放| 日韩大片免费观看网站| 欧美亚洲日本最大视频资源| 高清黄色对白视频在线免费看| 亚洲av电影在线观看一区二区三区| 女人被躁到高潮嗷嗷叫费观| 国产精品久久久久久人妻精品电影 | 久久精品国产a三级三级三级| 国产成人精品久久二区二区91 | 一区福利在线观看| 国产麻豆69| 老司机在亚洲福利影院| av有码第一页| 亚洲国产中文字幕在线视频| 久久久久久人妻| 久久99一区二区三区| 亚洲情色 制服丝袜| 日日啪夜夜爽| 青草久久国产| av电影中文网址| 国产午夜精品一二区理论片| 最近最新中文字幕免费大全7| 国产精品久久久久久精品电影小说| 一本一本久久a久久精品综合妖精| 亚洲国产av新网站| 免费黄色在线免费观看| 国产有黄有色有爽视频| 我的亚洲天堂| h视频一区二区三区| 观看av在线不卡| 老司机亚洲免费影院| 亚洲成人国产一区在线观看 | 热re99久久国产66热| 精品一区二区三卡| 国产精品秋霞免费鲁丝片| 女性被躁到高潮视频| 国产精品二区激情视频| 美女扒开内裤让男人捅视频| 九九爱精品视频在线观看| 交换朋友夫妻互换小说| 久久精品久久精品一区二区三区| 精品人妻在线不人妻| 美女主播在线视频| 国产日韩欧美视频二区| 天堂中文最新版在线下载| 两个人看的免费小视频| 十八禁高潮呻吟视频| 日韩av在线免费看完整版不卡| 国产1区2区3区精品| 欧美97在线视频| 免费观看a级毛片全部| 超碰成人久久| 午夜福利免费观看在线| 制服人妻中文乱码| 精品视频人人做人人爽| 日韩一卡2卡3卡4卡2021年| 一级黄片播放器| 亚洲成色77777| 亚洲国产中文字幕在线视频| 亚洲av成人不卡在线观看播放网 | 欧美日韩亚洲综合一区二区三区_| 久久午夜综合久久蜜桃| 久久久国产欧美日韩av| av福利片在线| 韩国高清视频一区二区三区| av网站免费在线观看视频| 中文字幕av电影在线播放| avwww免费| 精品国产乱码久久久久久小说| 久久久久网色| 国产av精品麻豆| 久久人人爽av亚洲精品天堂| 丝袜喷水一区| 精品国产一区二区久久| 国产精品国产三级国产专区5o| 精品国产露脸久久av麻豆| 中国国产av一级| 欧美av亚洲av综合av国产av | 亚洲成人免费av在线播放| 中文欧美无线码| 欧美国产精品va在线观看不卡| 91精品三级在线观看| 成人三级做爰电影| 久久久国产欧美日韩av| 精品国产一区二区久久| 日韩中文字幕欧美一区二区 | 精品午夜福利在线看| 如何舔出高潮| 在线观看免费日韩欧美大片| 国产精品99久久99久久久不卡 | 亚洲欧美一区二区三区国产| 成年美女黄网站色视频大全免费| 王馨瑶露胸无遮挡在线观看| 看免费成人av毛片| 91老司机精品| 天堂中文最新版在线下载| 国产片内射在线| 丝袜喷水一区| av片东京热男人的天堂| www.自偷自拍.com| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲欧美激情在线| 亚洲激情五月婷婷啪啪| 老司机深夜福利视频在线观看 | 欧美老熟妇乱子伦牲交| 色综合欧美亚洲国产小说| 日韩制服丝袜自拍偷拍| 婷婷色综合大香蕉| 最近的中文字幕免费完整| 欧美日韩一区二区视频在线观看视频在线| 好男人视频免费观看在线| 人妻一区二区av| 日本猛色少妇xxxxx猛交久久| 亚洲精品一区蜜桃| 少妇人妻 视频| 亚洲精品成人av观看孕妇| 精品免费久久久久久久清纯 | 国产精品国产三级国产专区5o| 成年美女黄网站色视频大全免费| 少妇猛男粗大的猛烈进出视频| 亚洲第一青青草原| 美女扒开内裤让男人捅视频| 日本欧美视频一区| 亚洲av男天堂| 伊人久久国产一区二区| 国产极品粉嫩免费观看在线| 如何舔出高潮| 日本色播在线视频| 亚洲精品中文字幕在线视频| 亚洲精品,欧美精品| 曰老女人黄片| 亚洲国产精品一区二区三区在线| 我要看黄色一级片免费的| www.av在线官网国产| 国产日韩一区二区三区精品不卡| 亚洲精华国产精华液的使用体验| 精品午夜福利在线看| 国产在线一区二区三区精| 亚洲天堂av无毛| 久久人人爽人人片av| 欧美在线一区亚洲| 在现免费观看毛片| 51午夜福利影视在线观看| 亚洲国产欧美网| 成年人免费黄色播放视频| 欧美中文综合在线视频| 男女无遮挡免费网站观看| 最近最新中文字幕免费大全7| 成年av动漫网址| 91aial.com中文字幕在线观看| 欧美日韩av久久| av.在线天堂| 好男人视频免费观看在线| 国产一区二区三区综合在线观看| 日本爱情动作片www.在线观看| 亚洲精品一二三| 国产男女内射视频| 日本欧美视频一区| 色综合欧美亚洲国产小说| 国产亚洲av片在线观看秒播厂| 人人妻人人澡人人看| 亚洲成人一二三区av| 亚洲av国产av综合av卡| 一边摸一边做爽爽视频免费| 色精品久久人妻99蜜桃| 国产精品久久久久久精品电影小说| av.在线天堂| 亚洲国产最新在线播放| 国产有黄有色有爽视频| 午夜福利影视在线免费观看| 99国产综合亚洲精品| 99久久综合免费| 岛国毛片在线播放| 久久 成人 亚洲| 日本黄色日本黄色录像| 青春草亚洲视频在线观看| 久久婷婷青草| 国产精品三级大全| 国产av一区二区精品久久| 国产男女内射视频| 黑人巨大精品欧美一区二区蜜桃| 男女边摸边吃奶| 亚洲一区中文字幕在线| 黑丝袜美女国产一区| 免费av中文字幕在线| 国产成人系列免费观看| 又粗又硬又长又爽又黄的视频| 涩涩av久久男人的天堂| 色婷婷av一区二区三区视频| 午夜福利视频在线观看免费| 日韩大片免费观看网站| 免费少妇av软件| 久久久久视频综合| 久久国产精品男人的天堂亚洲| 国产成人欧美在线观看 | 国产av精品麻豆| 大片免费播放器 马上看| 女人久久www免费人成看片| 久久99热这里只频精品6学生| 欧美人与善性xxx| 国产一区二区 视频在线| 亚洲成人国产一区在线观看 | 黄网站色视频无遮挡免费观看| 又大又爽又粗| 性色av一级| 亚洲精品,欧美精品| 国产精品一二三区在线看| 国产精品国产三级专区第一集| 国产亚洲av高清不卡| 又大又爽又粗| 国产精品久久久久久精品古装| 日本av免费视频播放| 亚洲精品中文字幕在线视频| 日韩伦理黄色片| 中文字幕制服av| 99热国产这里只有精品6| 亚洲精品国产av蜜桃| 别揉我奶头~嗯~啊~动态视频 | 国产亚洲av高清不卡| 亚洲精品久久久久久婷婷小说| 久热爱精品视频在线9| 色精品久久人妻99蜜桃| 性色av一级| 亚洲色图 男人天堂 中文字幕| 老汉色∧v一级毛片| 69精品国产乱码久久久| 亚洲美女搞黄在线观看| 成人国语在线视频| 精品国产一区二区三区四区第35| 美女扒开内裤让男人捅视频| 蜜桃在线观看..| 成人影院久久| 91精品三级在线观看| av国产精品久久久久影院| 在线天堂中文资源库| 91精品国产国语对白视频| 午夜福利免费观看在线| 我要看黄色一级片免费的| 亚洲国产精品一区二区三区在线| 国产 精品1| 亚洲av国产av综合av卡| 国产成人欧美| 亚洲欧洲国产日韩| 国产一区有黄有色的免费视频| 久久久久视频综合| 两个人免费观看高清视频| 国产高清国产精品国产三级| 日本黄色日本黄色录像| 亚洲成人手机| 你懂的网址亚洲精品在线观看| 久久精品亚洲熟妇少妇任你| 亚洲激情五月婷婷啪啪| 999精品在线视频| 欧美日韩福利视频一区二区| 精品亚洲成a人片在线观看| 成人国语在线视频| 一区二区日韩欧美中文字幕| 久久久久久久久免费视频了| 成人免费观看视频高清| 夜夜骑夜夜射夜夜干| 永久免费av网站大全| 国产黄色视频一区二区在线观看| 香蕉国产在线看| 999久久久国产精品视频| 精品一区二区免费观看| 国产日韩欧美在线精品| 国产亚洲最大av| 国产成人一区二区在线| 麻豆av在线久日| 久久人人爽人人片av| 制服诱惑二区| 国产欧美日韩一区二区三区在线| 色婷婷av一区二区三区视频| 国产又爽黄色视频| 欧美精品av麻豆av| 在线天堂中文资源库| 亚洲国产欧美网| 欧美xxⅹ黑人| 国产日韩欧美视频二区| 日本午夜av视频| 99久久人妻综合| 18禁裸乳无遮挡动漫免费视频| 国产精品香港三级国产av潘金莲 | 国产日韩欧美视频二区| 免费在线观看黄色视频的| 97在线人人人人妻| 日韩av免费高清视频| 另类亚洲欧美激情| 久久99热这里只频精品6学生| 波野结衣二区三区在线| 曰老女人黄片| 男人操女人黄网站| 精品亚洲成国产av| 亚洲欧洲日产国产| 中文字幕人妻丝袜一区二区 | 亚洲国产欧美日韩在线播放| 人人妻人人澡人人看| 免费黄频网站在线观看国产| 黑人猛操日本美女一级片| 亚洲久久久国产精品| 色播在线永久视频| 交换朋友夫妻互换小说| 蜜桃国产av成人99| 日韩大片免费观看网站| 午夜免费男女啪啪视频观看| 啦啦啦在线免费观看视频4| 色网站视频免费| 丝袜脚勾引网站| 国产黄频视频在线观看| 日日啪夜夜爽| 免费av中文字幕在线| 亚洲久久久国产精品| 99九九在线精品视频| 成人毛片60女人毛片免费| 一区二区三区激情视频| 日本vs欧美在线观看视频| 国产精品香港三级国产av潘金莲 | 999精品在线视频| 色婷婷av一区二区三区视频| 成人午夜精彩视频在线观看| 中文字幕av电影在线播放| av.在线天堂| 一级黄片播放器| bbb黄色大片| 国产又爽黄色视频| 91老司机精品| 欧美激情高清一区二区三区 | 人人妻人人澡人人看| 免费少妇av软件| 久久热在线av| 夫妻性生交免费视频一级片| 亚洲欧美精品综合一区二区三区| 免费观看a级毛片全部| 日本91视频免费播放| 亚洲精品国产色婷婷电影| 五月开心婷婷网| 成年人午夜在线观看视频| 91aial.com中文字幕在线观看| 高清在线视频一区二区三区| 欧美中文综合在线视频| 国产一区二区激情短视频 | av卡一久久| 在线看a的网站| 在线观看www视频免费| 精品一区二区免费观看| 国产精品国产av在线观看| 午夜福利视频精品| 美女大奶头黄色视频| 国产成人欧美| 777米奇影视久久| 亚洲男人天堂网一区| 欧美乱码精品一区二区三区| 婷婷色麻豆天堂久久| 999久久久国产精品视频| 大片电影免费在线观看免费| 久久鲁丝午夜福利片| 高清欧美精品videossex| 久久99精品国语久久久| 国产日韩欧美亚洲二区| 中国国产av一级| 最近中文字幕2019免费版| 日韩 亚洲 欧美在线| 日本一区二区免费在线视频| avwww免费| 日日摸夜夜添夜夜爱| 交换朋友夫妻互换小说| 亚洲专区中文字幕在线 | 91精品三级在线观看| 我要看黄色一级片免费的| 国产黄频视频在线观看| tube8黄色片| 青青草视频在线视频观看| 丝袜美腿诱惑在线| 久久精品久久久久久久性| 最新在线观看一区二区三区 | 久久鲁丝午夜福利片| 中文字幕精品免费在线观看视频| 欧美精品一区二区免费开放| 搡老岳熟女国产| 精品一区二区三区av网在线观看 | 少妇猛男粗大的猛烈进出视频| 亚洲av欧美aⅴ国产| 亚洲专区中文字幕在线 | 大香蕉久久网| 国产 精品1| 亚洲av成人精品一二三区| 久热这里只有精品99| 丁香六月欧美| 亚洲激情五月婷婷啪啪| 极品少妇高潮喷水抽搐| 黄片播放在线免费| 亚洲欧美一区二区三区黑人| 老司机靠b影院| 美女扒开内裤让男人捅视频| 女性生殖器流出的白浆| 少妇人妻精品综合一区二区| 操美女的视频在线观看| 人人妻人人爽人人添夜夜欢视频| 99久国产av精品国产电影| 国产精品无大码| 美女午夜性视频免费| 啦啦啦在线免费观看视频4| 欧美精品人与动牲交sv欧美| 久久久久久人妻| 国产视频首页在线观看| 国产av码专区亚洲av| 97精品久久久久久久久久精品| 色婷婷av一区二区三区视频| 国产成人一区二区在线| 五月天丁香电影| 又大又黄又爽视频免费| 免费看不卡的av| 最近手机中文字幕大全| 亚洲熟女精品中文字幕| 精品一品国产午夜福利视频| 亚洲成人手机| 亚洲精品,欧美精品| 国产免费又黄又爽又色| 久久 成人 亚洲| 国产黄色免费在线视频| 免费黄频网站在线观看国产| 最近中文字幕高清免费大全6| 伊人亚洲综合成人网| 自线自在国产av| 青春草亚洲视频在线观看| 国产免费福利视频在线观看| 18禁观看日本| 久久午夜综合久久蜜桃| 亚洲欧美一区二区三区黑人| 成年av动漫网址| 制服人妻中文乱码| 国产成人精品久久久久久| 天天躁夜夜躁狠狠躁躁| 制服人妻中文乱码| 亚洲国产精品一区二区三区在线| 熟妇人妻不卡中文字幕| 丰满迷人的少妇在线观看| 最近最新中文字幕大全免费视频 | 黑丝袜美女国产一区| 亚洲精品国产区一区二| 男女边吃奶边做爰视频| 欧美日韩亚洲综合一区二区三区_| 国产成人欧美| 男女高潮啪啪啪动态图| 久久人妻熟女aⅴ| 久久人人97超碰香蕉20202| 制服人妻中文乱码| 久久99一区二区三区| 老熟女久久久| 人人妻人人添人人爽欧美一区卜| 一级片免费观看大全| 男女午夜视频在线观看| 久久青草综合色|