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

    石津灌區(qū)冬小麥關(guān)鍵生育期逐時降水的變化特征

    2016-11-10 05:18:38任慶福翁白莎裴宏偉
    關(guān)鍵詞:乳熟期中雨歷時

    任慶福,翁白莎,裴宏偉,袁 喆

    (1 北京林業(yè)大學(xué) 水土保持學(xué)院,北京 100083;2 中國水利水電科學(xué)研究院 流域水循環(huán)模擬與調(diào)控國家重點(diǎn)實(shí)驗(yàn)室,北京 100038;3 中國科學(xué)院遺傳與發(fā)育生物學(xué)研究所 農(nóng)業(yè)資源研究中心,河北 石家莊 050022)

    ?

    石津灌區(qū)冬小麥關(guān)鍵生育期逐時降水的變化特征

    任慶福1,2,翁白莎2,裴宏偉3,袁喆2

    (1 北京林業(yè)大學(xué) 水土保持學(xué)院,北京 100083;2 中國水利水電科學(xué)研究院 流域水循環(huán)模擬與調(diào)控國家重點(diǎn)實(shí)驗(yàn)室,北京 100038;3 中國科學(xué)院遺傳與發(fā)育生物學(xué)研究所 農(nóng)業(yè)資源研究中心,河北 石家莊 050022)

    【目的】 研究氣候變化背景下石津灌區(qū)冬小麥關(guān)鍵生育期(拔節(jié)期、抽穗期、乳熟期)內(nèi)有效降水過程的變化,為灌區(qū)雨水的高效利用、制定科學(xué)合理的灌溉制度提供科學(xué)依據(jù)。【方法】 采用灌區(qū)內(nèi)部及周邊15個氣象站1983-2012年4-6月的逐時降水資料,基于降水歷時與降水強(qiáng)度的分類標(biāo)準(zhǔn),運(yùn)用Mann-Kendall趨勢檢驗(yàn),借助ArcGIS平臺的反距離插值(IDW)工具,分析了4類有效降水(短歷時小雨、長歷時小雨、短歷時中雨、長歷時中雨)在冬小麥3個生育期內(nèi)發(fā)生時刻的變化規(guī)律、逐站逐年的累計(jì)降水次數(shù)以及總降雨量的時空變化規(guī)律。【結(jié)果】 (1)午后至夜間(14:00-06:00)為4類有效降水的高發(fā)時段,其中長歷時中雨在拔節(jié)期的發(fā)生時刻集中在10:00-15:00時段。(2)4類有效降水的累計(jì)降水次數(shù)以及總降雨量30年的變化趨勢一致:①短歷時小雨在冬小麥拔節(jié)期、抽穗期呈增加趨勢,在乳熟期表現(xiàn)出減少趨勢;②長歷時小雨在拔節(jié)期呈增加趨勢,抽穗期、乳熟期為減少趨勢;③短歷時中雨、長歷時中雨在1998年以前呈減少趨勢,1998年以后短歷時中雨在拔節(jié)期為減少趨勢,其他時期表現(xiàn)為增加趨勢,長歷時中雨在抽穗期為減少趨勢,其他時期則為增加趨勢。(3)4類有效降水的累計(jì)降水次數(shù)以及總降雨量的空間分布一致:拔節(jié)期4種類型降水均表現(xiàn)為北多南少;抽穗期不同歷時小雨表現(xiàn)為中間少四圍多,不同歷時中雨則表現(xiàn)為西南多東北少;乳熟期不同歷時小雨的空間分布與抽穗期一致,短歷時中雨呈北多南少的分布,長歷時中雨則呈東多西少的分布格局。【結(jié)論】 氣候變化背景下冬小麥拔節(jié)期有效降水的發(fā)生次數(shù)和降雨量逐漸增加;乳熟期的有效降水次數(shù)和降雨量呈減少趨勢。

    石津灌區(qū);冬小麥生育期;逐時降水;有效降水;降水歷時

    充足的水分是保證作物高產(chǎn)穩(wěn)產(chǎn)的關(guān)鍵因子,在氣候變化背景下,華北平原降雨量減少已是不爭的事實(shí)[1-3]。盡管近50年來華北平原糧食產(chǎn)量持續(xù)增加,但這是以超采地下水為代價,華北平原灌溉用水占到水資源消耗總量的70%~75.2%[4-5],作物的水分供需矛盾日趨緊張。從長遠(yuǎn)來看,在灌溉用水保持零增長或負(fù)增長的條件下,實(shí)現(xiàn)糧食產(chǎn)量的提高,必須提高農(nóng)田水分利用效率。農(nóng)田水分利用效率主要包括灌溉水利用率、降雨利用率和作物水分利用率[6],但以往研究偏重于灌溉水利用效率[7]以及作物水分利用效率[8],對于降雨利用率的研究也只集中在降雨量上[9-10]。就農(nóng)作物而言,關(guān)鍵生育期內(nèi)的降水過程對作物產(chǎn)量的影響至關(guān)重要,降水的發(fā)生時刻、發(fā)生頻率、歷時、強(qiáng)度等指標(biāo)在時空上的變化很大程度上決定著降水對作物的有效性[11],然而對降水過程的這些指標(biāo)特征研究多見于氣象[12-14]、極端水文事件[15-16]以及城市暴雨[17-18]等領(lǐng)域,在作物的降雨利用方面卻鮮有報(bào)道,因而深入了解這些降雨特性對設(shè)計(jì)及管理農(nóng)田灌溉系統(tǒng)、提高農(nóng)田水分利用效率具有重要的理論指導(dǎo)意義。

    石津灌區(qū)是全國大型灌區(qū),是華北平原典型的井渠結(jié)合灌區(qū)。存在著灌溉水源不足、水分利用效率較低等問題,農(nóng)業(yè)用水矛盾突出。深入研究該灌區(qū)的降水過程,可以充分挖掘作物對降雨的利用潛力、緩解當(dāng)前灌區(qū)的用水矛盾。本研究利用石津灌區(qū)及周邊15個氣象站1983-2012年冬小麥關(guān)鍵生育期的逐時降水資料,以降水歷時、降水強(qiáng)度為分類指標(biāo),深入分析4類降水在灌區(qū)冬小麥3個關(guān)鍵生育期內(nèi)的發(fā)生時刻、發(fā)生頻次及雨量的時空變化規(guī)律,為灌區(qū)降雨的高效利用、制定科學(xué)合理的灌溉制度、提高農(nóng)田水分利用效率、實(shí)現(xiàn)作物的穩(wěn)產(chǎn)高產(chǎn)提供必要的理論依據(jù)。

    1 研究區(qū)概況

    石津灌區(qū)位于東經(jīng)114°19′~116°30′、北緯37°30′~38°18′,涉及冀縣、寧晉、晉縣、武強(qiáng)、武邑、深縣、深澤、藁城、衡水、趙縣、辛集等11個縣市,平均海拔30 m(圖1)。灌區(qū)屬于暖溫帶大陸性季風(fēng)氣候,春季多風(fēng)干旱、夏季高溫多雨,年均降水量為488 mm,且降水年內(nèi)分布不均,70%的全年降水主要集中在6-8月;多年平均水面蒸發(fā)量為1 000~1 200 mm。多年平均氣溫在12~13 ℃;年均日照時數(shù)為2 626 h,無霜期190~200 d;灌區(qū)土壤以褐土、潮土為主。灌區(qū)內(nèi)主要作物為冬小麥、夏玉米。

    圖 1 石津灌區(qū)地理位置及周邊氣象站點(diǎn)

    2 資料與方法

    2.1數(shù)據(jù)來源

    逐時降水?dāng)?shù)據(jù)采用國家氣象信息中心提供的《中國國家級地面氣象站逐小時降水?dāng)?shù)據(jù)集》,作物生育期數(shù)據(jù)采用中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)提供的《中國農(nóng)作物生長發(fā)育和農(nóng)田土壤濕度旬值數(shù)據(jù)集》與《中國生態(tài)系統(tǒng)定位觀測與研究數(shù)據(jù)集——農(nóng)田生態(tài)系統(tǒng)卷》(河北欒城站1998-2008)[19]。選取石津灌區(qū)內(nèi)部及其周邊共15個氣象站點(diǎn),為保證灌區(qū)降水?dāng)?shù)據(jù)的連續(xù)性及完整性,將1983-2012年4-6月作為研究時段(圖1)。

    2.2研究方法

    (1)小麥生育期的確定。作物物候期的觀測通常以植株高度、葉片數(shù)量和葉片長度以及開花時間等為主要指標(biāo),考慮到在灌區(qū)尺度上的觀測誤差,忽略氣溫的變化對冬小麥生育期的影響,本研究假定研究區(qū)在研究時段(1983-2012年)灌區(qū)內(nèi)冬小麥生育期不發(fā)生變化。根據(jù)藁城、欒城、深縣、南宮4個站點(diǎn)的作物生長發(fā)育日期及欒城試驗(yàn)站的觀測結(jié)果,確定石津灌區(qū)冬小麥關(guān)鍵生育期為拔節(jié)期(4月10日-5月1日)、抽穗期(5月1日-5月21日)以及乳熟期(5月21-6月10日)。

    (2)降水過程的確定及分類。確定降水間歇長度為2 h作為2次降水過程的分界,即當(dāng)一次降水過程連續(xù)2 h 沒有降水時,劃分為2次降水過程;將一次降水過程的全部時間定義為降水歷時(h);將該降水過程的雨量定義為次降雨量(mm);平均雨強(qiáng)(mm/h)定義為次降雨量與降水歷時之比,并將峰值雨強(qiáng)(mm/h)定義為一次降水過程中最大的單位降水雨強(qiáng);同時將文中的總降雨量定義為冬小麥某生育期的全部雨量(mm);降水次數(shù)定義為冬小麥某生育期內(nèi)發(fā)生降水的總次數(shù)(次)。

    根據(jù)降水平均雨強(qiáng)將降水過程分為小雨、中雨、大雨、暴雨4個等級;根據(jù)降水歷時將降水過程分為長歷時降水(≥6 h)和短歷時降水(<6 h)[13,15],因此根據(jù)降水歷時、降水平均雨強(qiáng)將降水過程分為8類(表1)。

    (3)冬小麥生育期有效降水類型的確定。殷水清等[15]認(rèn)為,降水強(qiáng)度≥10 mm/h的降水過程為極端降水事件,根據(jù)表1分類標(biāo)準(zhǔn),認(rèn)為短歷時大雨、短歷時暴雨、長歷時大雨、長歷時暴雨均為極端降水。極端降水可能會產(chǎn)生大量的徑流、引起深層滲漏,使得冬小麥對降水的利用率下降。而降水強(qiáng)度低、持續(xù)時間短的降水通常有效性最高,但這類降水存在被冬小麥冠層截留率高、入滲到土壤中易蒸發(fā)、持續(xù)性低等問題。由此可見,短歷時小雨對作物的有效性較高,但受到冬小麥冠層、蒸發(fā)、土壤入滲等因子的影響;而長歷時小雨、短歷時中雨、長歷時中雨盡管易發(fā)生深層滲漏,但入滲到土壤中的水分充足,消退時間較長。因此本研究重點(diǎn)分析這4類降水。

    冬小麥冠層截留量受到小麥品種、葉面積指數(shù)、降雨強(qiáng)度、降水歷時等因素的影響,本研究依據(jù)王迪等[20]的試驗(yàn)結(jié)果,取冬小麥在拔節(jié)期、抽穗期、乳熟期的最大冠層截留量分別為0.7,1.5,1.2 mm,認(rèn)為在冬小麥拔節(jié)期、抽穗期、乳熟期分別大于0.68 mm、1.47 mm、1.20 mm的降水為有效降水。

    (4)灌區(qū)降水特征指標(biāo)的選取及處理。本研究分別對冬小麥拔節(jié)、抽穗、乳熟等3個生育期的短歷時小雨、短歷時中雨、長歷時小雨、長歷時中雨的發(fā)生時刻、發(fā)生頻率、降水量等3個指標(biāo)進(jìn)行分析。各過程降雨量通過反距離權(quán)重(Inverse Distance Weighted,IDW)空間插值法求得;降水發(fā)生時刻、發(fā)生頻率按站點(diǎn)統(tǒng)計(jì)。運(yùn)用Mann-Kendall非參數(shù)檢驗(yàn)法對降水過程各指標(biāo)序列進(jìn)行趨勢分析;采用最小二乘法計(jì)算各降水指標(biāo)的氣候傾向率。

    應(yīng)用Origin9.0、ArcGIS10.0等軟件以及Fortran語言對試驗(yàn)數(shù)據(jù)進(jìn)行整理分析。

    表 1 降水過程的分類結(jié)果Table 1 Classification of precipitation

    3 結(jié)果與分析

    3.1各降水過程的統(tǒng)計(jì)特征

    對灌區(qū)內(nèi)部及周邊15個典型站點(diǎn)(晉縣、趙縣、欒城、寧晉、饒陽、深縣、安平、深澤、新河、武強(qiáng)、辛集、衡水、武邑、冀縣、棗強(qiáng))在冬小麥拔節(jié)期、抽穗期、乳熟期的降水進(jìn)行分類,并統(tǒng)計(jì)各降雨過程的總降雨量、總降水次數(shù)、次降雨量、降雨歷時、平均雨強(qiáng)、峰值雨強(qiáng),結(jié)果見表2和表3。由表2可以看出,石津灌區(qū)1983-2012年冬小麥關(guān)鍵生育期15個站點(diǎn)逐年逐站統(tǒng)計(jì),共發(fā)生 6 151 次降水,總降雨量為29 151.9 mm,平均每站點(diǎn)每年降雨量為64.8 mm。從總降水次數(shù)來看,1983-2012年灌區(qū)在冬小麥拔節(jié)期、抽穗期、乳熟期15個站點(diǎn)分別發(fā)生了1 695,2 222,2 234次降水,其中短歷時降水發(fā)生了1 238,1 601,1 759次,分別占到總降水次數(shù)的73.0%,72.1%,78.7%。3個生育期小雨發(fā)生次數(shù)分別為1 608,2 055,1 990次,分別占到總降水次數(shù)的94.9%,92.5%,89.1%。短歷時小雨在3個生育期發(fā)生次數(shù)分別為1 179,1 482,1 562次,逐漸增多,相應(yīng)的總降雨量也逐漸增加,每次降雨量都在1.25~1.28 mm。3個生育期內(nèi)均未發(fā)生長歷時暴雨,而短歷時暴雨在拔節(jié)、抽穗2個時期內(nèi)各發(fā)生1次,在乳熟期發(fā)生了8次;長歷時大雨在抽穗期發(fā)生了1次,僅在乳熟期發(fā)生了2次,但平均每次降雨量可達(dá)68.75 mm,由此可見灌區(qū)冬小麥在乳熟期有極端降水的發(fā)生。

    表 2 各降水類型降雨量與降雨次數(shù)逐站統(tǒng)計(jì)結(jié)果(1983-2012年)Table 2 Amount and frequency of different precipitation at each station (1983-2012)

    注:“-”表示無此數(shù)據(jù)。表3同。

    Note:“-” means no data.The same for Table 3.

    表 3 各降水類型降水歷時和降水強(qiáng)度逐站統(tǒng)計(jì)結(jié)果(1983-2012年)Table 3 Duration and intensity of different precipitation at each station (1983-2012)

    由表3可以看出,短歷時降雨在冬小麥3個生育期內(nèi)的平均歷時為2.2 h,長歷時降雨平均歷時為9.7 h;短歷時降雨在3個生育期之間的降雨歷時、平均雨強(qiáng)、降雨峰值變化不大,其中短歷時小雨平均歷時為2.22 h、平均雨強(qiáng)為0.51 mm/h、降雨峰值為0.86 mm/h;短歷時中雨也具有同樣的規(guī)律,平均歷時為2.66 h、平均雨強(qiáng)為4.31 mm/h、降雨峰值為8.14 mm/h;長歷時小雨在乳熟期降水歷時有所減少,但降雨峰值可達(dá)3.16 mm/h;長歷時中雨在拔節(jié)期平均雨強(qiáng)為3.45 mm/h。極端降水在3個生育期內(nèi)的降雨歷時、平均雨強(qiáng)、降雨峰值變化較大,短歷時暴雨在拔節(jié)、抽穗期平均歷時均為1.00 h,但在乳熟期平均歷時達(dá)到2.00 h;短歷時大雨的平均雨強(qiáng)在抽穗期為9.42 mm/h,在乳熟期達(dá)到10.53 mm/h。

    將次降雨量超過冬小麥最大冠層截留量的4類降水(短歷時小雨、短歷時中雨、長歷時小雨、長歷時中雨)進(jìn)行重新統(tǒng)計(jì),結(jié)果見表4。由表4可知,灌區(qū)長、短歷時小雨平均發(fā)生次數(shù)為3.0次/(a·站),而長、短歷時中雨平均發(fā)生次數(shù)平均為0.2~0.7次/(a·站),因此本研究對長、短歷時小雨2類降水的時間尺度為1年;對長、短歷時中雨2類降水的時間尺度為5年。

    表 4 冬小麥不同生育期4類有效降水次數(shù)及有效降水量逐站統(tǒng)計(jì)結(jié)果(1983-2012年)Table 4 Amount and frequency of the 4 effective precipitation types at different wheat growth stages at each station (1983-2012)

    3.24類有效降水過程的發(fā)生時刻

    3.2.1不同歷時小雨發(fā)生時刻逐年變化由圖2-a可以看出,短歷時小雨在冬小麥3個生育期發(fā)生總次數(shù)整體上呈現(xiàn)“午后”單峰型變化,主要集中在15:00-20:00。其中拔節(jié)期短歷時小雨主要發(fā)生在16:00-18:00,在17:00發(fā)生了40次,在20:00達(dá)到一個峰值;乳熟期的發(fā)生時刻略有提前,在16:00發(fā)生了39次;在抽穗期短歷時小雨主要發(fā)生在早晨08:00以前、晚上21:00以后;3個生育期在08:00-13:00時段發(fā)生短歷時小雨頻率低,拔節(jié)期在13:00只發(fā)生過7次。由圖2-b可以看出,長歷時小雨整體上呈現(xiàn)“夜間-午后”雙峰型變化,但3個生育期的“峰值”時刻各不相同。其中夜間時段拔節(jié)期集中在02:00-06:00,抽穗期集中在00:00-08:00,乳熟期則集中在01:00-11:00;午后時段3個生育期都在15:00以后。

    圖 2 不同歷時小雨在各時刻的發(fā)生次數(shù)(1983-2012年)

    不同歷時小雨在各時刻發(fā)生次數(shù)的逐年變化(1983-2012年)見圖3。

    圖 3不同歷時小雨在各時刻發(fā)生次數(shù)的逐年變化(1983-2012年)

    Ⅰ.拔節(jié)期;Ⅱ.抽穗期;Ⅲ.乳熟期;a.短歷時小雨;b.長歷時小雨

    Fig.3Annual variation of frequency of short and long light rain types at different hours (1983-2012)

    Ⅰ.Jointing stage;Ⅱ.Heading stage;Ⅲ.Milk ripe stage;a.Short light rain;b.Long light rain

    圖3-Ⅰ-a、Ⅱ-a、Ⅲ-a顯示了短歷時小雨在冬小麥3個生育期各時刻發(fā)生次數(shù)的逐年變化。從變化趨勢上看,在1992年之前,拔節(jié)期短歷時小雨的發(fā)生主要集中在上午時段,1992年之后,趨向于午后15:00以后及夜間時段;抽穗期這種趨勢更加明顯,且2002年以后,短歷時小雨主要發(fā)生在午后13:00-20:00;在乳熟期,短歷時小雨的發(fā)生時刻呈現(xiàn)從“午后”型向“白天”型變化,表現(xiàn)為在1992年之前發(fā)生時刻主要集中在15:00-21:00時段,1992年以后逐漸向中午時刻集中,2006年13:00達(dá)7次。

    圖3-Ⅰ-b、Ⅱ-b、Ⅲ-b顯示了長歷時小雨在冬小麥3個生育期各時刻發(fā)生次數(shù)的逐年變化。從變化趨勢來看,1994年之前,長歷時小雨在冬小麥拔節(jié)、乳熟期發(fā)生時刻分布于05:00-20:00時段,在抽穗期發(fā)生時刻集中在00:00-08:00;1994年以后,降雨發(fā)生時刻趨向于夜間,拔節(jié)期發(fā)生時刻呈“夜間-傍晚”型,抽穗期集中在15:00以后,乳熟期向深夜00:00-06:00時段集中。

    3.2.2不同歷時中雨發(fā)生時刻的變化圖4-a中短歷時中雨在冬小麥3個生育期呈單峰型,發(fā)生時刻主要集中在14:00-01:00,在拔節(jié)期30年15個站點(diǎn)短歷時中雨未在01:00、03:00-08:00、10:00、12:00-13:00時段發(fā)生過,但在22:00發(fā)生過8次;17:00是乳熟期短歷時中雨發(fā)生的高頻時段,發(fā)生了24次;08:00-12:00是短歷時中雨發(fā)生低頻時段,在冬小麥3個生育時期共發(fā)生了18次。圖4-b中長歷時中雨在冬小麥拔節(jié)期的發(fā)生時刻主要集中在午間10:00-15:00時段,呈單峰型;長歷時中雨在抽穗期主要集中在18:00-06:00,屬于夜間型降雨;在乳熟期降水發(fā)生時刻呈“凌晨-下午”雙峰型,主要發(fā)生在04:00-06:00、17:00-18:00時段。

    圖 4 不同歷時中雨在各時刻的發(fā)生次數(shù)(1983-2012年)

    3.34類有效降水過程的降水次數(shù)及雨量

    3.3.1不同歷時小雨降水次數(shù)及降雨量的時空變化灌區(qū)冬小麥3個生育期不同歷時小雨的降水次數(shù)、降雨量1983-2012年的變化趨勢檢驗(yàn)結(jié)果見表5。表5表明,降水次數(shù)與降雨量變化趨勢一致。在拔節(jié)期短歷時小雨、長歷時小雨的總降水次數(shù)與總降雨量均呈增加的趨勢,且短歷時小雨的降水次數(shù)與總降雨量增加趨勢顯著(α=0.05),15個站點(diǎn)短歷時小雨的總降水次數(shù)以6.08次/10a的速率增加,總降雨量的氣候傾向率為16.16 mm/10a;長歷時小雨的總降水次數(shù)則以1.18次/10a的速率增加,總降雨量的增加速率為12.29 mm/10a。在抽穗期,短歷時小雨的總降水次數(shù)與總降雨量仍然呈增加趨勢,但增加趨勢不顯著,而長歷時小雨的總降水次數(shù)與總降雨量則呈微弱下降趨勢,降水次數(shù)下降速率為2.00次/10a,降雨量氣候傾向率為-23.26 mm/10a。在乳熟期,短歷時小雨、長歷時小雨的總降水次數(shù)與總降雨量均呈減少趨勢,且長歷時小雨的降雨次數(shù)與降雨量顯著減少(α=0.05),長歷時小雨的降雨量在乳熟期以49.23 mm/10a的速率減少,降水次數(shù)的減少速率為4.63次/10a,且在α=0.10水平上降雨量也以1.47 mm/次的速率顯著減少。

    表 5 不同歷時小雨降水次數(shù)及降雨量的Mann-Kendall趨勢檢驗(yàn)結(jié)果(1983-2012年)Table 5 Frequency and amount of short and long light rain types by Mann-Kendall test (1983-2012)

    注:**表示顯著水平為0.05;*表示顯著水平為0.10。

    Note:** indicates the significant level is 0.05,and * indicates the significant level is 0.10.

    不同歷時小雨降水次數(shù)與降雨量的空間分布見圖5。

    圖 5不同歷時小雨降水次數(shù)與降雨量的空間分布

    Ⅰ.拔節(jié)期;Ⅱ.抽穗期;Ⅲ.乳熟期;a.短歷時小雨;b.長歷時小雨

    Fig.5Spatial distribution of frequency and amount of short and long light rain types

    Ⅰ.Jointing stage;Ⅱ.Heading stage;Ⅲ.Milk ripe stage;a.Short light rain;b.Long light rain

    圖5-Ⅰ-a、Ⅰ-b顯示了1983-2012年灌區(qū)在冬小麥拔節(jié)期短歷時小雨、長歷時小雨降雨量的空間分布:短歷時小雨30年的總降雨量為90.5 mm,平均3.0 mm/a,長歷時小雨30年的總降雨量為 315.6 mm,平均10.5 mm/a。降水發(fā)生次數(shù)與降雨量空間分布一致,短歷時小雨表現(xiàn)為沿深澤、辛集、衡水一帶分別向東北、西南方向遞減,長歷時小雨表現(xiàn)為灌區(qū)西北部的藁城、晉縣向東南方向遞減;短歷時小雨發(fā)生次數(shù)最多的為衡水站,發(fā)生了47次,平均1.5次/a,降雨量也達(dá)到109.2 mm,平均2.3 mm/次;但在衡水站長歷時小雨30年只發(fā)生了22次,降雨量在290 mm以下。

    圖5-Ⅱ-a、Ⅱ-b顯示了1983-2012年灌區(qū)在冬小麥抽穗期短歷時小雨、長歷時小雨降雨量的空間分布,其中短歷時小雨30年的總降雨量為83.9 mm,平均2.8 mm/a,長歷時小雨30年的總降雨量為342.3 mm,平均11.4 mm/a。與拔節(jié)期不同,此期短歷時小雨、長歷時小雨的降雨次數(shù)與降雨量的空間分布均表現(xiàn)為由灌區(qū)四圍向中北部遞減,沿晉縣站、深縣站一帶為灌區(qū)降雨量最低值,短歷時小雨在深縣站30年降雨次數(shù)為19次,降雨量在77 mm以下;長歷時小雨在冀縣站30年的降雨次數(shù)達(dá)到45次,降雨量394.6 mm,平均13.2 mm/次。

    圖5-Ⅲ-a、Ⅲ-b顯示了1983-2012年灌區(qū)在冬小麥乳熟期短歷時小雨、長歷時小雨降雨量的空間分布,其中短歷時小雨30年的總降雨量為100.2 mm,平均3.3 mm/a,長歷時小雨30年的總降雨量為270.8 mm,平均9.0 mm/a。乳熟期短歷時小雨的降雨次數(shù)及降雨量空間分布表現(xiàn)為由深縣、辛集、寧晉一帶分別向東、西方向遞增,在藁城、晉縣一帶達(dá)到降雨量最大值,在107.2 mm以上,晉縣站降雨次數(shù)為34次,深縣站降雨次數(shù)為24次,相差10次;長歷時小雨則是沿深縣、辛集、趙縣一帶分別向南、北兩方向遞增,分別在灌區(qū)南部的新河、冀州、衡水一帶,北部地區(qū)的深澤達(dá)到降雨量最大,在286 mm以上,但深澤站降雨次數(shù)為25次,降雨量為323.4 mm,平均13.0 mm/次。

    3.3.2不同歷時中雨降雨量的時空變化灌區(qū)冬小麥3個生育期不同歷時中雨降雨量在5年尺度上的變化見圖6。

    圖 6 不同歷時中雨5年尺度降雨量的變化(1983-2012年)

    由圖6可以看出,長歷時中雨降雨量整體上要大于短歷時中雨,從1983-1997年除乳熟期的短歷時中雨降雨量有所增加(1988-1992年降雨量達(dá)到516 mm)外,短(長)歷時中雨在冬小麥3個生育期均呈下降趨勢,且在1993-1997年,冬小麥拔節(jié)期、抽穗期并未發(fā)生長歷時中雨。1998年以后,拔節(jié)期短歷時中雨降雨量表現(xiàn)出減少的趨勢,長歷時中雨則呈增加趨勢;在抽穗期、乳熟期,這種趨勢相反,短歷時中雨降雨量表現(xiàn)出增加趨勢,而長歷時中雨降雨量則表現(xiàn)出降低的趨勢。

    由圖7可以看出,整體上看短歷時中雨、長歷時中雨在冬小麥3個生育期的降水次數(shù)與降水量的空間分布一致。在冬小麥拔節(jié)期,短歷時中雨降雨量為37.8 mm,多年平均1.3 mm,其空間分布整體上自灌區(qū)北向南方向遞減,灌區(qū)北部的深澤降雨量在65.5 mm以上,深澤站降雨次數(shù)為7次,是15個站點(diǎn)中短歷時中雨發(fā)生次數(shù)最多的,灌區(qū)降雨量最小值位于東北部的武強(qiáng)、武邑一帶,降雨量在32.9 mm以下,只發(fā)生了2~3次降水(圖7-Ⅰ-a)。長歷時中雨降雨量為82.0 mm,多年平均2.7 mm,空間分布整體上自灌區(qū)東北向西南方向遞減。灌區(qū)東北部的深縣、武強(qiáng)一帶降雨量在112.3~148.5 mm,與短歷時中雨相比,盡管降雨量大,但降水次數(shù)整體上較少,其中欒城、趙縣、冀縣站只發(fā)生過1次,晉縣、深縣站無長歷時中雨發(fā)生(圖7-Ⅰ-b)。

    圖 7 不同歷時中雨降水次數(shù)與降雨量的空間分布

    在冬小麥抽穗期,短歷時中雨降雨量為71.6 mm,多年平均2.4 mm,其空間分布整體上自灌區(qū)西南向東北方向遞減,在灌區(qū)西北部的深澤、晉縣一帶降雨量在56.3 mm以下,2個站點(diǎn)的降雨次數(shù)分別為3次、4次,灌區(qū)東南部的寧晉、辛集降雨量在87.4 mm以上,降水次數(shù)也達(dá)到9次(圖7-Ⅱ-a)。長歷時中雨降雨量為145.4 mm,多年平均4.8 mm,其空間分布整體上與短歷時中雨分布一致,但降雨量高值區(qū)位于灌區(qū)中西部的辛集,降雨量在222.3 mm以上,降水次數(shù)為6次,寧晉站長歷時中雨只發(fā)生過2次,因而在灌區(qū)西南端的寧晉地區(qū)降雨量也偏低(圖7-Ⅱ-b)。

    在冬小麥乳熟期,短歷時中雨降雨量為123.5 mm,多年平均4.1 mm,其空間分布整體上與拔節(jié)期相似,均為自北向南遞減,灌區(qū)北部的深縣降雨量在147.0 mm以上,降水次數(shù)達(dá)到14次,降雨量少的地區(qū)位于灌區(qū)東南部的武邑、衡水一帶,降雨量在91.5 mm以下,降水次數(shù)為6次(圖7-Ⅲ-a)。長歷時中雨降雨量為90.0 mm,多年平均3.0 mm,空間分布表現(xiàn)為自東向西遞減,灌區(qū)最東端的武強(qiáng)地區(qū)降雨量在163.9 mm以上,武強(qiáng)站降水次數(shù)為8次,到灌區(qū)西部的晉縣、趙縣、寧晉一帶,降雨量在 41.9~81.5 mm,寧晉站30年來只發(fā)生過1次短歷時中雨(圖7-Ⅲ-b)。

    4 討 論

    (1)不同歷時、不同雨型降水的發(fā)生時刻,對冬小麥生育期內(nèi)降水的吸收更為重要。短歷時小雨發(fā)生于18:00-05:00,此時處于冬小麥蒸散作用微弱時段,根系吸收少,加之降雨歷時短,降水很快入滲到土壤中,且降水在土壤中的消退時間短,因此該時段冬小麥對降水的利用率低;而短歷時小雨發(fā)生于白天,冬小麥根系吸水作用增強(qiáng),該時段的降水對冬小麥最為有效;同樣,短歷時中雨、長歷時中雨可使降水充分入滲到土壤中,土壤水分達(dá)到飽和,土壤水分消退時間長,冬小麥對該類降雨利用率高。

    本研究結(jié)果表明,短歷時小雨在冬小麥拔節(jié)期、乳熟期主要發(fā)生在16:00-18:00,在抽穗期則發(fā)生在早晨08:00以前、晚上21:00以后;短歷時中雨在冬小麥3個生育期的發(fā)生時刻主要集中在14:00-01:00;長歷時小雨在冬小麥3個生育期的發(fā)生均呈現(xiàn)“夜間-午后”型,長歷時中雨在冬小麥拔節(jié)期的發(fā)生時刻主要集中在午間10:00-15:00時段,在抽穗期主要集中在18:00-06:00,屬于夜間型降雨,在乳熟期降水發(fā)生時刻主要集中在04:00-06:00、17:00-18:00時段。Yu等[12]研究指出,中國中東部地區(qū),長歷時降水的峰值雨強(qiáng)出現(xiàn)在清晨,短歷時降水易出現(xiàn)在下午和傍晚;李建等[13]研究顯示,北京地區(qū)夏季(6-8月)長歷時降水的峰值多出現(xiàn)在后半夜至清晨,而短歷時降水峰值多集中在下午至前半夜;殷水清等[15]研究認(rèn)為,海河流域夏季長歷時降水峰值發(fā)生在早上,短歷時降水峰值發(fā)生在下午。由此對于短歷時小雨、短歷時中雨的發(fā)生,與前人研究結(jié)果一致,均發(fā)生于午后至翌日凌晨01:00。而對于長歷時小雨與長歷時中雨,本研究結(jié)果顯示這2種降雨的發(fā)生時刻更趨于白天、或“后半夜-午后”,與前人研究結(jié)果并不一致,主要原因:①前人研究時段主要集中在夏季(6-8月),而本研究時段主要集中在春季(4-5月);②前人研究基于降水歷時描述了降水量的日變化特征,曲線峰值表示了降雨量的大小[12]或某時刻有無降水發(fā)生[15],但并未考慮不同的降水強(qiáng)度等級,歷時短的降水發(fā)生時刻與峰值強(qiáng)度出現(xiàn)時刻相差不大,但歷時長的降水發(fā)生時刻與峰值強(qiáng)度出現(xiàn)時刻間隔時間較長,因此不同的降水強(qiáng)度下,長歷時降雨的發(fā)生時刻與長歷時降雨的日變化規(guī)律并不一致。

    (2)研究也指出,海河流域1960-2010年春季降水呈增加趨勢[21-22]。本研究時段主要集中在春季(4-5月),短歷時小雨降雨量在冬小麥拔節(jié)期、抽穗期均呈增加趨勢,進(jìn)入乳熟期之后呈減小趨勢,長歷時小雨降雨量在乳熟期顯著減少。整體上與前人研究基本一致,由于本研究將春季(4-5月)劃分為3個階段(拔節(jié)期、抽穗期、乳熟期),且在冬小麥乳熟期包括了初夏,因此在冬小麥乳熟期,短歷時小雨、長歷時小雨降雨量的變化呈減少的趨勢。1998年為短歷時中雨、長歷時中雨變化趨勢的轉(zhuǎn)折點(diǎn)。

    5 結(jié) 論

    (1)灌區(qū)1983-2012年4類降水發(fā)生時刻:短歷時小雨在冬小麥拔節(jié)期、乳熟期主要發(fā)生在16:00-18:00,在抽穗期則發(fā)生在早晨08:00以前、晚上21:00以后,短歷時中雨在冬小麥3個生育期的發(fā)生時刻主要集中在14:00-01:00;長歷時小雨在冬小麥3個生育期的發(fā)生均呈現(xiàn)“夜間-午后”型,長歷時中雨在冬小麥拔節(jié)期的發(fā)生時刻主要集中在午間10:00-15:00時段,在抽穗期主要集中在18:00-06:00,在乳熟期降水發(fā)生時刻主要集中在04:00-06:00、17:00-18:00時段。

    (2)灌區(qū)1983-2012年4類降水在冬小麥3個生育期內(nèi)累計(jì)次數(shù)與總降雨量的時間變化一致:①短歷時小雨在冬小麥拔節(jié)期、抽穗期呈增加趨勢,在乳熟期表現(xiàn)出減少趨勢;②長歷時小雨在拔節(jié)期呈增加趨勢,在抽穗期、乳熟期為減少趨勢;③短歷時中雨、長歷時中雨在1998年以前呈減少趨勢,1998年以后短歷時中雨在拔節(jié)期為減少趨勢,其他時期表現(xiàn)為增加趨勢,長歷時中雨在抽穗期為減少趨勢,其他時期則為增加趨勢。

    (3)灌區(qū)1983-2012年冬小麥3個生育期內(nèi)4類降水累計(jì)降水次數(shù)與總降雨量的空間分布一致:①拔節(jié)期短歷時小雨表現(xiàn)為沿深澤、辛集、衡水一帶分別向東北、西南方向遞減;長歷時小雨表現(xiàn)為由灌區(qū)西北部的藁城、晉縣向東南方向遞減;短歷時中雨整體上自北向南遞減;長歷時中雨則表現(xiàn)為自東北向西南方向遞減。②抽穗期短歷時小雨、長歷時小雨均表現(xiàn)為由灌區(qū)四圍向中北部遞減;短歷時中雨、長歷時中雨空間分布一致,整體上均表現(xiàn)為自西南向東北方向遞減。③乳熟期短歷時小雨表現(xiàn)為由深縣、辛集、寧晉一帶分別向東、西方向遞增;長歷時小雨則是沿深縣、辛集、趙縣一帶分別向南、北兩方向遞增;短歷時中雨空間分布整體上與拔節(jié)期相似,均為自北向南遞減;長歷時中雨則表現(xiàn)為自東向西遞減。

    (4)本研究假定了灌區(qū)在1983-2012年冬小麥物候期并未發(fā)生變化,忽略了氣溫的影響。而研究指出,氣溫的升高使得河北省冬小麥在1998年以后,起身期和抽穗期較1991-1997年間各提早了5 d,拔節(jié)期提早了4 d,返青期和成熟期提早了3 d[23];90年代以后冬小麥播種期比80年代推遲2~5 d,成熟期提早1~2 d,生育期縮短了2~7 d[24]。因此需要對溫度影響下的冬小麥關(guān)鍵生育期的降水變化做進(jìn)一步探討。

    [1]Zhai P M,Zhang X B,Wan H,et al.Trends in total precipitation and frequency of daily precipitation extremes over China [J].Journal of Climate,2005,18(7):1096-1108.

    [2]Qian W,Lin X.Regional trends in recent precipitation indices in China [J].Meteorology and Atmospheric Physics,2005,90(3/4):193-207.

    [3]姚惠明,吳永祥,關(guān)鐵生.中國降水演變趨勢診斷及其新事實(shí) [J].水科學(xué)進(jìn)展,2013,24(1):1-10.

    Yao H M,Wu Y X,Guan T S.Diagnose of precipitation evolution trend in China and new facts [J].Advances in Water Science,2013,24(1):1-10.

    [4]姚治君,林耀明,高迎春,等.華北平原分區(qū)適宜性農(nóng)業(yè)節(jié)水技術(shù)與潛力 [J].自然資源學(xué)報(bào),2000,15(3):259-264.

    Yao Z J,Lin Y M,Gao Y C,et al.The water-saving technologies and potentials for regional suitable agriculture in North China Plain [J].Journal of Natural Resources,2000,15(3):259-264.

    [5]張光輝,劉中培,費(fèi)宇紅,等.華北平原區(qū)域水資源特征與作物布局結(jié)構(gòu)適應(yīng)性研究 [J].地球?qū)W報(bào),2010,31(1):17-22.

    Zhang G H,Liu Z P,Fei Y H,et al.The relationship between the distribution of irrigated crops and the supply capability of regional water resources in North China Plain [J].Acta Geoscientic Sinica,2010,31(1):17-22.

    [6]王會肖,劉昌明.作物水分利用效率內(nèi)涵及研究進(jìn)展 [J].水科學(xué)進(jìn)展,2000,11(1):99-104.

    Wang H X,Liu C M.Advances in crop water use efficiency research [J].Advances in Water Science,2000,11(1):99-104.

    [7]陳皓銳,黃介生,伍靖偉,等.井渠結(jié)合灌區(qū)用水效率指標(biāo)尺度效應(yīng)研究框架 [J].農(nóng)業(yè)工程學(xué)報(bào),2009,25(8):1-7.

    Chen H R,Huang J S,Wu J W,et al.Scale effect research framework of irrigation water use efficiency indices in well-canal combined irrigation area [J].Transactions of the Chinese Society of Agricultural Engineering,2009,25(8):1-7.

    [8]張玉翠,沈彥俊,裴宏偉,等.冬小麥生長季蒸散與水分利用效率變化特征 [J].南水北調(diào)與水利科技,2010,8(5):39-41.

    Zhang Y C,Shen Y J,Pei H W,et al.Characteristics of evapotranspiration and water use efficiency of typical irrigated winter wheat in NCP [J].South to North Water Transfers and Water Science & Technology,2010,8(5):39-41.

    [9]Sun H Y,Shen Y J,Yu Q,et al.Effect of precipitation change on water balance and WUE of the winter wheat-summer maize rotation in the North China Plain [J].Agricultural Water Management,2010,97(8):1139-1145.

    [10]邵立威,張喜英,陳素英,等.降水、灌溉和品種對玉米產(chǎn)量和水分利用效率的影響 [J].灌溉排水學(xué)報(bào),2009(1):48-51.

    Shao L W,Zhang X Y,Chen S Y,et al.Yield and water use efficiency affected by rain fall,irrigation and maize varieties [J].Journal of Irrigation and Drainage,2009(1):48-51.

    [11]徐鳳琴.有效降水量淺析 [J].氣象水文海洋儀器,2009,26(1):96-100.

    Xu F Q.Concise analysis of effective precipitation [J].Meteorological Hydrological and Marine Instruments,2009,26(1):96-100.

    [12]Yu R,Zhou T,Xiong A.Diurnal variations of summer precipitation over contiguous China [J].Geophysical Research Letters,2007,34(1):L1704.

    [13]李建,宇如聰,王建捷.北京市夏季降水的日變化特征 [J].科學(xué)通報(bào),2008,53(7):829-832.

    Li J,Yu R C,Wang J J.The diurnal variation of summer precipitation in Beijing [J].Chinese Science Bulletion,2008,53(7):829-832.

    [14]宇如聰,李建,陳昊明,等.中國大陸降水日變化研究進(jìn)展 [J].氣象學(xué)報(bào),2014(5):948-968.

    Yu R C,Li J,Chen H M,et al.Progress in studies of the precipitation diurnal variation over contiguous China [J].Acta Meteorologica Sinica,2014(5):948-968.

    [15]殷水清,高歌,李維京,等.1961~2004年海河流域夏季逐時降水變化趨勢 [J].中國科學(xué):地球科學(xué),2012(2):256-266.

    Yin S Q,Gao G,Li W J,et al.Long-term precipitation changeby hourly data in Haihe River Basin during 1961-2004 [J].Sci China:Earth Sci,2012(2):256-266.

    [16]李書嚴(yán),馬京津,軒春怡,等.1951~2008年北京極端天氣事件分析 [J].氣候與環(huán)境研究,2012,17(2):244-250.

    Li S Y,Ma J J,Xuan C Y,et al.Analysis of extreme weather events in Beijing during 1951-2008 [J].Climatic and Environmental Research,2012,17(2):244-250.

    [17]孫嫻,魏娜,王式功,等.西安夏季降水的日變化特征研究 [J].冰川凍土,2013,35(5):1176-1182.

    Sun X,Wei N,Wang S G,et al.Study of the daily variation of summer precipitation in Xi’an [J].Journal of Glaciology and Geocryology,2013,35(5):1176-1182.

    [18]劉偉東,尤煥苓,任國玉,等.北京地區(qū)精細(xì)化的降水變化特征 [J].氣候與環(huán)境研究,2014,19(1):61-68.

    Liu W D,You H L,Ren G Y,et al.Subtle precipitation characteristics in Beijing area [J].Climatic and Environmental Research,2014,19(1):61-68.

    [19]胡春勝,程一松.中國生態(tài)系統(tǒng)定位觀測與研究數(shù)據(jù)集:農(nóng)田生態(tài)系統(tǒng)卷:河北欒城站(1998-2008) [M].北京:中國農(nóng)業(yè)出版社,2011.

    Hu C S,Cheng Y S.The datasets of the Chinese ecosystem research network:volume of agricultural ecosystem:Luancheng station in Hebei Province 1998-2008 [M].Beijing:Chinese Agriculture Press,2011.

    [20]王迪,李久生,饒敏杰.噴灌冬小麥冠層截留試驗(yàn)研究 [J].中國農(nóng)業(yè)科學(xué),2006,39(9):1859-1864.

    Wang D, Li J S, Rao M J.Winter wheat canopy interception under sprinkler irrigation [J].Scientia Agriculture Sinica,2006,39(9):1859-1864.

    [21]王利娜,朱厚華,魯帆,等.海河流域近50年降水量時空變化特征分析 [J].干旱地區(qū)農(nóng)業(yè)研究,2012,30(2):242-246.

    Wang L N,Zhu H H,Lu F,et al.Characteristics of temporal and spatial variation of precipitation in Haihe River Basin during recent 50 years [J].Agricultural Research in the Arid Areas,2012,30(2):242-246.

    [22]Chu J L,Xia J,Xu C Y,et al.Spatial and temporal variability of daily precipitation in Haihe River Basin,1958-2007 [J].Journal of Geographical Sciences,2010(2):248-260.

    [23]李元華,車少靜.河北省溫度和降水變化對農(nóng)業(yè)的影響 [J].中國農(nóng)業(yè)氣象,2005,26(4):224-228.

    Li Y H,Che S J.Impacts of changes of temperature and precipitation on agriculture in Hebei Regions [J].Chinese Journal of Agrometeogrogy,2005,26(4):224-228.

    [24]谷永利,林艷,李元華.氣溫變化對河北省冬小麥主要發(fā)育期的影響分析 [J].干旱區(qū)資源與環(huán)境,2007,21(12):141-145.

    Gu Y L,Lin Y,Li Y H.The analysis on the effect of temperature fluctuation on primary growing periods of winter wheat in Hebei [J].Jounal of Arid Land Resources and Environment,2007,21(12):141-145.

    Hourly precipitation during key growth periods of winter wheat in Shijin irrigation district

    REN Qingfu1,2,WENG Baisha2,PEI Hongwei3,YUAN Zhe2

    (1CollegeofSoilandWaterConservation,BeijingForestryUniversity,Beijing100083,China;2StateKeyLaboratoryofSimulationandRegulationofWaterCycleinRiverBasin,ChinaInstituteofWaterResourcesandHydropowerResearch,Beijing100038,China;3CenterforAgriculturalResourcesResearch,InstituteofGeneticsandDevelopmentalBiology,ChineseAcademyofSciences,Shijiazhuang,Hebei050022,China)

    【Objective】 This study investigated effective precipitation process on key growth periods of winter wheat (jointing,heading,and milk ripe) under climate change in Shijin irrigation district to provide scientific basis for rainfall utilization and reasonable irrigation schedule.【Method】 Based on the hourly precipitation data from 15 meteorological stations in study area on April,May,and June from 1983 to 2012,the precipitation was classified according to duration and intensity.The occurrence time,frequency and amount of the four effective precipitation types (short light rain,long light rain,short moderate rain,and long moderate rain) in three growth stages of winter wheat were analyzed with the Mann-Kendall trend test and the inverse distance interpolation (IDW) method.【Result】 (1) The occurrence time of the effective precipitation types was mainly from afternoon to morning (14:00-06:00),and long moderate rain mainly occurred at jointing stage during 10:00-15:00.(2) The trend and amount of all precipitation types were same.The short light rain increased at jointing and heading stages while decreased at milk ripe stage.Long light rain increased at jointing stage while decreased at heading and milk ripe stages.Short and long moderate rain types decreased during the three growth stages before 1998.After then,short moderate rain increased at heading and milk ripe stages and long moderate rain increased at jointing and milk ripe stages.(3) The spatial distribution of frequency and amount of the precipitation types was also identical.All precipitation types at jointing stage were high in north and low in south.At heading stage,the amounts of the short and long light rain types decreased from border to center,while the short and long moderate rain types increased from northeast to southwest.Distribution of short and long light rain types at milk ripe stage was consistent at heading stage.The amount of short moderate rain at milk ripe stage was high in north and low in south while that of long moderate rain was high in east and low in west. 【Conclusion】 The frequency and amount of effective precipitation would increase at jointing stage and decrease at milk ripe stage of winter wheat under climate change.

    Shijin irrigation district;winter wheat growth periods;hourly precipitation;effective precipitation;precipitation duration

    時間:2016-08-0909:40DOI:10.13207/j.cnki.jnwafu.2016.09.013

    2015-02-09

    國家“十二五”科技支撐計(jì)劃項(xiàng)目“氣候變化對水資源影響與風(fēng)險(xiǎn)評估技術(shù)”(2012BAC19B03);國家自然科學(xué)基金面上項(xiàng)目“基于水資源系統(tǒng)的廣義干旱風(fēng)險(xiǎn)評價與風(fēng)險(xiǎn)區(qū)劃研究”(51279207)

    任慶福(1983-),男,山西汾陽人,博士,主要從事農(nóng)業(yè)水文水資源和氣候變化研究。E-mail:wxws.2008@163.com

    翁白莎(1986-),女,福建晉江人,工程師,博士,主要從事變化環(huán)境下干旱風(fēng)險(xiǎn)評價與應(yīng)對研究。

    E-mail:wengbs@iwhr.com

    S271

    A

    1671-9387(2016)09-0092-13

    網(wǎng)絡(luò)出版地址:http://www.cnki.net/kcms/detail/61.1390.S.20160809.0940.026.html

    猜你喜歡
    乳熟期中雨歷時
    起落架用30CrMnSiNi2A鋼表面Cd-Ti鍍層在海水和鹽霧環(huán)境中的腐蝕行為研究
    下 雨
    秸稈還田與灌溉方式對春玉米根冠抗氧化能力的影響
    量詞“只”的形成及其歷時演變
    常用詞“怠”“惰”“懶”的歷時演變
    淺談跨座式單軌交通建設(shè)中雨污水管線遷改工作
    飼用小黑麥在濟(jì)源地區(qū)的刈割期及刈割次數(shù)研究
    無題(6)
    對《紅樓夢》中“不好死了”與“……好的”的歷時考察
    攀西地區(qū)冬閑田種植晚熟型燕麥的最佳刈割期研究
    中國奶牛(2019年1期)2019-02-15 07:19:42
    久久久精品免费免费高清| 久久久久久久久久久免费av| 在线a可以看的网站| 夜夜看夜夜爽夜夜摸| 尾随美女入室| 黄色配什么色好看| 三级男女做爰猛烈吃奶摸视频| 国产精品不卡视频一区二区| 国产白丝娇喘喷水9色精品| 色综合色国产| 亚洲精品亚洲一区二区| 亚洲成人精品中文字幕电影| 亚洲欧洲国产日韩| 午夜久久久久精精品| 国产高清三级在线| 国产探花在线观看一区二区| 国产精品不卡视频一区二区| .国产精品久久| 国产亚洲精品av在线| 日韩,欧美,国产一区二区三区| 看十八女毛片水多多多| 亚洲av中文av极速乱| 成人毛片60女人毛片免费| freevideosex欧美| 亚洲成人中文字幕在线播放| 亚洲欧美成人综合另类久久久| 97超视频在线观看视频| 免费大片18禁| 久久午夜福利片| 综合色av麻豆| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 一夜夜www| 久久久久久久久大av| 五月伊人婷婷丁香| 美女国产视频在线观看| 91久久精品国产一区二区三区| 嫩草影院新地址| 丰满少妇做爰视频| 国产免费视频播放在线视频 | 大香蕉97超碰在线| 久久精品夜色国产| 久久久久久久亚洲中文字幕| 久久久久国产网址| 日韩在线高清观看一区二区三区| 搞女人的毛片| 黄色一级大片看看| 国产欧美日韩精品一区二区| 国产男女超爽视频在线观看| 亚洲欧美一区二区三区黑人 | 人妻系列 视频| 干丝袜人妻中文字幕| 亚洲人成网站在线播| 中文在线观看免费www的网站| 丝瓜视频免费看黄片| 极品少妇高潮喷水抽搐| 天堂影院成人在线观看| 人妻一区二区av| www.av在线官网国产| 在线a可以看的网站| 永久网站在线| av免费观看日本| 亚洲av成人av| 久久久久久久午夜电影| 午夜精品一区二区三区免费看| 尾随美女入室| 久久精品夜夜夜夜夜久久蜜豆| 国产高清有码在线观看视频| 日本三级黄在线观看| 国产一区二区三区综合在线观看 | 高清午夜精品一区二区三区| 久久人人爽人人爽人人片va| 亚洲精品第二区| 少妇的逼水好多| 日韩一本色道免费dvd| 国产av不卡久久| 日韩亚洲欧美综合| 看黄色毛片网站| 免费高清在线观看视频在线观看| 久久精品久久精品一区二区三区| 国产精品蜜桃在线观看| 欧美日韩综合久久久久久| www.色视频.com| 国产精品久久视频播放| 国产 一区 欧美 日韩| 国产熟女欧美一区二区| 国产国拍精品亚洲av在线观看| 国产人妻一区二区三区在| 亚洲丝袜综合中文字幕| 国产精品福利在线免费观看| 久久97久久精品| 晚上一个人看的免费电影| 日韩国内少妇激情av| 日本午夜av视频| 女人十人毛片免费观看3o分钟| 高清毛片免费看| 精品不卡国产一区二区三区| 2018国产大陆天天弄谢| 午夜爱爱视频在线播放| 内射极品少妇av片p| 日本色播在线视频| 一级毛片我不卡| 十八禁国产超污无遮挡网站| 亚洲真实伦在线观看| 一级二级三级毛片免费看| 国产69精品久久久久777片| 久久久久久久久久成人| 高清日韩中文字幕在线| 国产精品美女特级片免费视频播放器| 国产精品99久久久久久久久| 国产成人91sexporn| 久久国内精品自在自线图片| 欧美高清成人免费视频www| 免费高清在线观看视频在线观看| 国产av国产精品国产| 亚洲久久久久久中文字幕| 九九在线视频观看精品| 成人特级av手机在线观看| 国产乱来视频区| 国产探花极品一区二区| 美女黄网站色视频| 麻豆乱淫一区二区| 网址你懂的国产日韩在线| 最新中文字幕久久久久| 高清欧美精品videossex| 中文字幕亚洲精品专区| 国产精品久久久久久久电影| 最近视频中文字幕2019在线8| 亚洲成人一二三区av| 免费看光身美女| 亚洲伊人久久精品综合| 老司机影院成人| 国产又色又爽无遮挡免| www.色视频.com| 好男人在线观看高清免费视频| 日韩欧美一区视频在线观看 | 午夜福利成人在线免费观看| 午夜福利在线观看免费完整高清在| 国产亚洲91精品色在线| 国产单亲对白刺激| 天堂中文最新版在线下载 | 亚洲内射少妇av| 亚洲精品成人av观看孕妇| 国产探花极品一区二区| 国产亚洲91精品色在线| 久久草成人影院| 亚洲最大成人中文| 国产一区亚洲一区在线观看| 麻豆av噜噜一区二区三区| 看免费成人av毛片| 亚洲欧洲日产国产| 中文字幕免费在线视频6| 国产av码专区亚洲av| 亚洲熟妇中文字幕五十中出| 亚洲成人一二三区av| 黄色日韩在线| 亚洲精品国产成人久久av| 国产精品伦人一区二区| 日本黄大片高清| a级毛片免费高清观看在线播放| 男人狂女人下面高潮的视频| 国产一级毛片七仙女欲春2| or卡值多少钱| 国产大屁股一区二区在线视频| 午夜福利网站1000一区二区三区| 一夜夜www| 老女人水多毛片| 国产永久视频网站| 日韩av在线大香蕉| 免费人成在线观看视频色| 日韩欧美精品v在线| 性插视频无遮挡在线免费观看| 十八禁国产超污无遮挡网站| 日本色播在线视频| 能在线免费观看的黄片| 在线a可以看的网站| 久久久久久国产a免费观看| 嘟嘟电影网在线观看| 久久久久久久久久人人人人人人| 搡老乐熟女国产| 久久久久久久大尺度免费视频| 99视频精品全部免费 在线| a级毛片免费高清观看在线播放| 午夜福利成人在线免费观看| 婷婷色麻豆天堂久久| 久久97久久精品| 99热这里只有是精品50| 精品人妻视频免费看| videossex国产| 美女被艹到高潮喷水动态| 91aial.com中文字幕在线观看| 丝袜美腿在线中文| 亚洲人成网站在线观看播放| 色综合亚洲欧美另类图片| 水蜜桃什么品种好| 人妻制服诱惑在线中文字幕| 综合色av麻豆| 老师上课跳d突然被开到最大视频| ponron亚洲| 国产黄片视频在线免费观看| 日韩视频在线欧美| 五月伊人婷婷丁香| 国产人妻一区二区三区在| 少妇的逼水好多| 久久精品国产自在天天线| 成年人午夜在线观看视频 | 日本猛色少妇xxxxx猛交久久| 男人狂女人下面高潮的视频| 中文字幕av成人在线电影| 午夜福利视频精品| 久久亚洲国产成人精品v| 日韩制服骚丝袜av| 成人鲁丝片一二三区免费| 精品国产一区二区三区久久久樱花 | 乱码一卡2卡4卡精品| 午夜福利在线在线| 国产精品熟女久久久久浪| 不卡视频在线观看欧美| 97在线视频观看| 成人亚洲精品av一区二区| 亚洲av日韩在线播放| 午夜亚洲福利在线播放| 日韩欧美一区视频在线观看 | 国产淫语在线视频| 精品久久久久久电影网| 干丝袜人妻中文字幕| 国产一级毛片七仙女欲春2| 国产高清有码在线观看视频| kizo精华| 床上黄色一级片| 女的被弄到高潮叫床怎么办| 婷婷色麻豆天堂久久| 性色avwww在线观看| 国产午夜精品论理片| 日韩大片免费观看网站| 欧美极品一区二区三区四区| 中文欧美无线码| 午夜精品在线福利| 国产精品1区2区在线观看.| 久久久久久久午夜电影| 天堂网av新在线| 91av网一区二区| 国产老妇伦熟女老妇高清| 黄色一级大片看看| 午夜亚洲福利在线播放| 亚洲精品一二三| 七月丁香在线播放| 亚洲不卡免费看| 亚洲av国产av综合av卡| 成人av在线播放网站| 波野结衣二区三区在线| 国产黄色小视频在线观看| 国产免费又黄又爽又色| 色吧在线观看| 久久久久久久久久久丰满| 自拍偷自拍亚洲精品老妇| 最近手机中文字幕大全| 女人久久www免费人成看片| 久久久精品94久久精品| 日韩电影二区| 国产精品久久久久久精品电影| 天堂俺去俺来也www色官网 | 午夜免费激情av| 欧美激情国产日韩精品一区| 成人特级av手机在线观看| 老师上课跳d突然被开到最大视频| 少妇人妻一区二区三区视频| 99热6这里只有精品| 九九在线视频观看精品| 亚洲高清免费不卡视频| 亚洲欧美日韩东京热| 色综合站精品国产| 夜夜爽夜夜爽视频| 99热网站在线观看| 麻豆乱淫一区二区| 国产乱人视频| 亚洲精品日本国产第一区| 国产精品1区2区在线观看.| 美女xxoo啪啪120秒动态图| 久久久a久久爽久久v久久| 国产精品一区二区性色av| 精品久久国产蜜桃| 国产精品伦人一区二区| .国产精品久久| 丰满人妻一区二区三区视频av| 一本一本综合久久| 十八禁网站网址无遮挡 | av卡一久久| 人妻系列 视频| .国产精品久久| xxx大片免费视频| 三级毛片av免费| 亚洲无线观看免费| 亚洲av成人精品一二三区| 99久国产av精品国产电影| 亚洲色图av天堂| 在线观看美女被高潮喷水网站| 日韩欧美一区视频在线观看 | 91精品一卡2卡3卡4卡| 国产成人免费观看mmmm| 亚洲国产精品国产精品| 少妇丰满av| 三级国产精品片| av网站免费在线观看视频 | 久久久久久久久久久免费av| 69av精品久久久久久| 97超碰精品成人国产| av又黄又爽大尺度在线免费看| 成年免费大片在线观看| 午夜福利网站1000一区二区三区| 好男人在线观看高清免费视频| 人体艺术视频欧美日本| 午夜老司机福利剧场| 国产成人freesex在线| 久久久久久国产a免费观看| 高清毛片免费看| 五月伊人婷婷丁香| 久久久久久久大尺度免费视频| 亚洲真实伦在线观看| 国产亚洲av嫩草精品影院| 久久久久久久久久久免费av| 男女视频在线观看网站免费| h日本视频在线播放| 观看美女的网站| 日韩成人伦理影院| 欧美一区二区亚洲| 熟妇人妻久久中文字幕3abv| 亚洲精品久久午夜乱码| 日韩,欧美,国产一区二区三区| 成人漫画全彩无遮挡| 久久久精品免费免费高清| 一级片'在线观看视频| 国产精品麻豆人妻色哟哟久久 | 国产欧美日韩精品一区二区| 夫妻午夜视频| 国产成人午夜福利电影在线观看| 97在线视频观看| 一级毛片aaaaaa免费看小| 国产男人的电影天堂91| 成年女人在线观看亚洲视频 | 久久久色成人| 永久网站在线| 免费观看的影片在线观看| 人体艺术视频欧美日本| 午夜福利视频1000在线观看| 日本黄大片高清| 亚洲最大成人手机在线| 午夜福利在线观看免费完整高清在| 一级毛片我不卡| 国产男女超爽视频在线观看| 18禁裸乳无遮挡免费网站照片| 亚洲色图av天堂| 亚洲国产精品成人综合色| 日韩国内少妇激情av| 视频中文字幕在线观看| 成人高潮视频无遮挡免费网站| 欧美性感艳星| 天天一区二区日本电影三级| 两个人视频免费观看高清| 高清视频免费观看一区二区 | 欧美最新免费一区二区三区| 18禁在线无遮挡免费观看视频| 国产黄片美女视频| 亚洲高清免费不卡视频| 全区人妻精品视频| 亚洲最大成人手机在线| 男的添女的下面高潮视频| 亚洲欧美中文字幕日韩二区| 亚洲精品乱码久久久久久按摩| 久久久久精品性色| 久久久亚洲精品成人影院| 男人和女人高潮做爰伦理| 欧美三级亚洲精品| 国产精品熟女久久久久浪| 寂寞人妻少妇视频99o| 国产精品嫩草影院av在线观看| av.在线天堂| 亚洲四区av| 欧美三级亚洲精品| 麻豆成人午夜福利视频| 日韩精品青青久久久久久| 99热6这里只有精品| 久久久午夜欧美精品| a级毛色黄片| 久久99蜜桃精品久久| 干丝袜人妻中文字幕| 精品人妻视频免费看| 国产v大片淫在线免费观看| 少妇人妻一区二区三区视频| 亚洲久久久久久中文字幕| 日本wwww免费看| 尤物成人国产欧美一区二区三区| 我的女老师完整版在线观看| 国产成人a区在线观看| 日韩三级伦理在线观看| 国产色婷婷99| 国国产精品蜜臀av免费| 久久久精品免费免费高清| 色综合站精品国产| 麻豆成人av视频| 午夜亚洲福利在线播放| 成人亚洲精品一区在线观看 | 高清日韩中文字幕在线| 国产高清三级在线| 亚洲国产欧美人成| 亚洲国产日韩欧美精品在线观看| 国产高清有码在线观看视频| 成人亚洲欧美一区二区av| 欧美xxxx黑人xx丫x性爽| 熟女电影av网| 乱系列少妇在线播放| 久久久亚洲精品成人影院| 麻豆乱淫一区二区| 看非洲黑人一级黄片| 亚洲精品国产av蜜桃| www.色视频.com| 国产精品麻豆人妻色哟哟久久 | 不卡视频在线观看欧美| 美女内射精品一级片tv| 身体一侧抽搐| 麻豆乱淫一区二区| 3wmmmm亚洲av在线观看| 国产大屁股一区二区在线视频| 国产精品一区二区三区四区久久| 国产伦精品一区二区三区四那| 搡老乐熟女国产| 久久精品国产亚洲av涩爱| 日日干狠狠操夜夜爽| 伦理电影大哥的女人| 一级毛片我不卡| 人妻夜夜爽99麻豆av| 久久精品夜色国产| 精品久久久久久久久亚洲| 亚洲精品一区蜜桃| 日韩av免费高清视频| 国产欧美日韩精品一区二区| 欧美日韩视频高清一区二区三区二| 搞女人的毛片| 亚洲国产最新在线播放| 久久久亚洲精品成人影院| 三级国产精品片| 精品久久久久久久久亚洲| 伊人久久国产一区二区| 国产亚洲91精品色在线| 国产伦理片在线播放av一区| 男女视频在线观看网站免费| 身体一侧抽搐| 天堂俺去俺来也www色官网 | 亚洲精品第二区| 亚洲精品日韩av片在线观看| 直男gayav资源| 亚洲熟妇中文字幕五十中出| 尾随美女入室| 午夜福利在线观看吧| 三级国产精品欧美在线观看| 边亲边吃奶的免费视频| 一二三四中文在线观看免费高清| 全区人妻精品视频| av女优亚洲男人天堂| 美女脱内裤让男人舔精品视频| 丝袜美腿在线中文| 三级国产精品欧美在线观看| 国产69精品久久久久777片| 亚洲欧美成人综合另类久久久| 国产大屁股一区二区在线视频| 久久人人爽人人爽人人片va| 伦精品一区二区三区| 少妇猛男粗大的猛烈进出视频 | 亚洲国产精品专区欧美| 99久久精品一区二区三区| 我要看日韩黄色一级片| 少妇高潮的动态图| 能在线免费看毛片的网站| 中文天堂在线官网| 国产精品伦人一区二区| 高清日韩中文字幕在线| 赤兔流量卡办理| 天美传媒精品一区二区| 国产精品.久久久| 亚洲精品国产av成人精品| 国产永久视频网站| 亚洲国产成人一精品久久久| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美日本视频| 久久韩国三级中文字幕| 身体一侧抽搐| 久久草成人影院| 国产 亚洲一区二区三区 | 免费人成在线观看视频色| 欧美一区二区亚洲| 精品不卡国产一区二区三区| 午夜激情福利司机影院| 国产精品av视频在线免费观看| 青春草视频在线免费观看| 一级毛片久久久久久久久女| 国产有黄有色有爽视频| 国产乱来视频区| 草草在线视频免费看| 午夜精品一区二区三区免费看| 成人漫画全彩无遮挡| 老司机影院毛片| 国产成人aa在线观看| 欧美bdsm另类| 可以在线观看毛片的网站| 一级毛片久久久久久久久女| 高清午夜精品一区二区三区| 人妻系列 视频| 菩萨蛮人人尽说江南好唐韦庄| 免费观看精品视频网站| 啦啦啦中文免费视频观看日本| 免费观看av网站的网址| 人妻系列 视频| 精品一区二区免费观看| 久久久欧美国产精品| 欧美zozozo另类| 午夜福利在线观看免费完整高清在| 99热这里只有是精品在线观看| 国产精品久久久久久精品电影| 欧美变态另类bdsm刘玥| 91久久精品国产一区二区成人| 亚洲第一区二区三区不卡| 亚洲一区高清亚洲精品| 免费不卡的大黄色大毛片视频在线观看 | 亚洲国产精品专区欧美| 看黄色毛片网站| 好男人视频免费观看在线| 91久久精品电影网| 晚上一个人看的免费电影| 中文精品一卡2卡3卡4更新| 亚洲性久久影院| 午夜老司机福利剧场| 国产精品国产三级专区第一集| 国产久久久一区二区三区| 国产精品日韩av在线免费观看| 白带黄色成豆腐渣| 人妻少妇偷人精品九色| 搡女人真爽免费视频火全软件| av专区在线播放| 91精品伊人久久大香线蕉| 中文字幕av成人在线电影| 九九爱精品视频在线观看| 久久久精品94久久精品| 三级男女做爰猛烈吃奶摸视频| 一级黄片播放器| 欧美bdsm另类| 欧美变态另类bdsm刘玥| 亚洲国产最新在线播放| 亚洲乱码一区二区免费版| 国产精品麻豆人妻色哟哟久久 | 免费观看av网站的网址| 国产一级毛片在线| 亚洲av福利一区| www.av在线官网国产| 亚洲欧美中文字幕日韩二区| 亚洲欧美成人综合另类久久久| 小蜜桃在线观看免费完整版高清| 一区二区三区免费毛片| 人人妻人人澡欧美一区二区| 国产大屁股一区二区在线视频| 免费看美女性在线毛片视频| 大话2 男鬼变身卡| 亚洲av电影不卡..在线观看| 欧美另类一区| 免费av不卡在线播放| a级毛色黄片| 国产一区亚洲一区在线观看| 免费大片黄手机在线观看| 精品一区二区三区视频在线| 精品国产三级普通话版| 日韩在线高清观看一区二区三区| 午夜福利视频1000在线观看| 97人妻精品一区二区三区麻豆| 欧美潮喷喷水| 男女视频在线观看网站免费| 日产精品乱码卡一卡2卡三| videos熟女内射| 深爱激情五月婷婷| 麻豆久久精品国产亚洲av| 精品久久久久久成人av| 激情 狠狠 欧美| 国产午夜福利久久久久久| 亚洲av电影不卡..在线观看| 综合色av麻豆| 天堂中文最新版在线下载 | 啦啦啦韩国在线观看视频| 亚洲aⅴ乱码一区二区在线播放| 男女边摸边吃奶| 亚洲国产精品国产精品| 国产综合精华液| 一区二区三区高清视频在线| 精品亚洲乱码少妇综合久久| 噜噜噜噜噜久久久久久91| 插逼视频在线观看| 成人二区视频| 久久久久久国产a免费观看| 欧美区成人在线视频| 亚洲一区高清亚洲精品| 麻豆国产97在线/欧美| 人人妻人人澡人人爽人人夜夜 | 边亲边吃奶的免费视频| 精品亚洲乱码少妇综合久久| 少妇裸体淫交视频免费看高清| 一区二区三区高清视频在线| 精品人妻一区二区三区麻豆| 一级毛片久久久久久久久女| 纵有疾风起免费观看全集完整版 | 少妇人妻一区二区三区视频| 久久这里有精品视频免费| 亚洲va在线va天堂va国产| 婷婷色综合大香蕉| kizo精华| 国产淫片久久久久久久久| 美女国产视频在线观看| 亚洲欧美日韩东京热| 秋霞在线观看毛片| 亚洲av成人av|