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

    稻田徑流易發(fā)期不同類型肥料的氮素流失風(fēng)險(xiǎn)

    2017-08-11 08:46:28侯朋福薛利祥俞映倞江瑜汪吉東薛利紅楊林章
    關(guān)鍵詞:田面穗肥易發(fā)

    侯朋福,薛利祥,俞映倞,江瑜,汪吉東,薛利紅*,楊林章

    (1.江蘇省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與環(huán)境研究所,南京 210014;2.中國(guó)農(nóng)業(yè)科學(xué)院作物科學(xué)研究所,北京 100081)

    稻田徑流易發(fā)期不同類型肥料的氮素流失風(fēng)險(xiǎn)

    侯朋福1,薛利祥1,俞映倞1,江瑜2,汪吉東1,薛利紅1*,楊林章1

    (1.江蘇省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與環(huán)境研究所,南京 210014;2.中國(guó)農(nóng)業(yè)科學(xué)院作物科學(xué)研究所,北京 100081)

    依據(jù)60年典型自然降雨資料分析了太湖地區(qū)稻田徑流易發(fā)期,結(jié)合不同肥料減量運(yùn)籌長(zhǎng)期定位試驗(yàn)(2009年開(kāi)始)近3年的徑流監(jiān)測(cè)數(shù)據(jù),比較了徑流易發(fā)期內(nèi)不同處理的實(shí)際徑流氮素控制效果。結(jié)果表明,太湖地區(qū)水稻生長(zhǎng)前期降雨概率較大,6月21日—7月6日單日降雨概率均超過(guò)50%。從階段統(tǒng)計(jì)結(jié)果來(lái)看,基肥期和蘗肥期(含基肥-蘗肥階段)降雨概率和降雨量均明顯高于穗肥期、蘗肥-穗肥階段和穗肥后至成熟階段,降雨概率分別達(dá)48.15%和49.81%,降雨量分別為12.81 mm和12.84 mm,均超過(guò)實(shí)際監(jiān)測(cè)到徑流的同期最低降雨量11.5 mm,產(chǎn)生徑流可能性較大,是稻季徑流易發(fā)期。從定位試驗(yàn)的實(shí)際產(chǎn)生徑流和田面水氮素濃度結(jié)果來(lái)看,易發(fā)期徑流和田面水氮素以銨態(tài)氮為主,硝態(tài)氮差異不顯著。與常規(guī)用量分次施肥處理(CN)相比,化肥減量?jī)?yōu)化處理(RF)和有機(jī)無(wú)機(jī)減量配施處理(OCN)徑流易發(fā)期徑流氮素濃度分別較CN處理平均降低8.83%和19.18%。緩控釋肥減量替代處理(SCU)的基肥期和基肥-蘗肥階段徑流氮素濃度和田面水氮素濃度明顯高于其他處理,徑流氮素濃度較CN處理分別增加了20.3%和11.72%,但蘗肥期徑流氮素濃度減少30.72%。全有機(jī)肥減量替代處理(OF)肥期徑流和田面水氮素濃度降低,其中基肥期徑流氮素濃度較CN處理降低9.04%、蘗肥期降低28.53%,但基肥-蘗肥階段氮素濃度較CN處理增加了19.7%,增加了氮素徑流損失風(fēng)險(xiǎn)。不同氮肥減量措施能夠降低易發(fā)期內(nèi)不同階段徑流氮素濃度,但在徑流易發(fā)期的徑流氮素?fù)p失控制效果不能一概而論。

    稻田;徑流;易發(fā)期;肥料類型;氮素流失

    隨著工業(yè)廢水和城市生活污水等點(diǎn)源污染得到控制,農(nóng)村面源污染控制逐步成為現(xiàn)代農(nóng)業(yè)和社會(huì)可持續(xù)發(fā)展的重大課題,農(nóng)田氮肥面源污染作為其中的重要組成部分備受關(guān)注[1]。太湖地區(qū)水資源豐富,是我國(guó)農(nóng)業(yè)集約化種植區(qū)域之一,水稻種植面積占耕地總面積的75%[2-3]。王桂良等[4]基于文獻(xiàn)數(shù)據(jù)的Meta分析發(fā)現(xiàn),氮素徑流損失是長(zhǎng)江流域單季稻生長(zhǎng)期的主要損失途徑。課題組前期研究表明,太湖地區(qū)不同氮肥管理模式下稻田氮素環(huán)境排放總量為40.4~85.8 kg·hm-2,其中徑流氮損失為6.8~9.2 kg·hm-2,占氮素環(huán)境排放的10.7%~23%,是氮素?fù)p失的主要途徑[5]。為了控制稻田氮素?fù)p失,研究人員通過(guò)系列研究逐步形成了調(diào)整肥料用量及運(yùn)籌方式、新型緩控釋肥料替代和有機(jī)肥料替代等系列源頭減量控制技術(shù)[6]。稻田田面在田埂保護(hù)下形成封閉徑流體系,只有在特殊情況如降雨發(fā)生時(shí),田面水才會(huì)溢出形成機(jī)會(huì)徑流,因此稻田徑流的發(fā)生具有突發(fā)性偶然性的特征,年際間隨降雨波動(dòng)較大[6-7]。因此,不同減量技術(shù)由于不同時(shí)期肥料投入差異較大,季節(jié)氮素徑流損失可能呈現(xiàn)不同的階段損失特征[8-10]。目前對(duì)徑流氮素?fù)p失控制效果的評(píng)價(jià)主要依靠徑流池或徑流收集管等對(duì)徑流收集測(cè)定后比較隨機(jī)發(fā)生徑流的氮素?fù)p失強(qiáng)度,對(duì)不同時(shí)期徑流發(fā)生風(fēng)險(xiǎn)關(guān)注較少[10-13]。不同類型肥料的養(yǎng)分釋放特性差異較大[14]。尿素等無(wú)機(jī)化肥施入稻田后迅速溶于水,短時(shí)間內(nèi)田面水氮素濃度升至峰值,有機(jī)肥由于本身成分的復(fù)雜性,其養(yǎng)分釋放通常呈現(xiàn)一個(gè)長(zhǎng)期的波動(dòng)過(guò)程,而緩控釋肥料由于肥料的緩慢釋放,田面水氮素濃度變化也呈現(xiàn)出長(zhǎng)期的變化過(guò)程[9,15-17]。因此不同減量技術(shù)氮素?fù)p失主要與徑流-施肥時(shí)間間隔及作物的不同生育期有關(guān),短期監(jiān)測(cè)數(shù)據(jù)的分析可能由于年際間降雨的差異而無(wú)法真實(shí)反映或掩蓋單項(xiàng)技術(shù)的實(shí)際流失特征及控制效果。太湖地區(qū)地處長(zhǎng)江下游,水稻種植季降雨頻發(fā),本研究擬利用歷史氣象數(shù)據(jù),分析區(qū)域多年份降雨數(shù)據(jù),尋求太湖地區(qū)徑流發(fā)生的易發(fā)期,并結(jié)合長(zhǎng)期定位監(jiān)測(cè)試驗(yàn)結(jié)果,對(duì)不同源頭減量技術(shù)的徑流氮素?fù)p失控制效果進(jìn)行綜合評(píng)價(jià),提出適宜該區(qū)域的徑流氮素減排技術(shù)。

    1 材料與方法

    1.1 歷史氣象資料分析

    徑流是多種因素綜合發(fā)生的,受降雨量、降雨強(qiáng)度、田間植被覆蓋度/LAI、土壤水分狀況或水深等影響[18]。而稻田田面在田埂保護(hù)下形成封閉徑流體系,只有在特殊情況如降雨發(fā)生時(shí),田面水才會(huì)溢出形成機(jī)會(huì)徑流[6],因此作者在本文中根據(jù)徑流發(fā)生受驅(qū)動(dòng)因子的先后將降雨劃分為前驅(qū)動(dòng)因子,將田間植被和土壤水分劃分為后驅(qū)動(dòng)因子。而有研究表明只要有徑流發(fā)生,徑流初始期就是氮素流失的高風(fēng)險(xiǎn)期。由于田間和土壤條件時(shí)空變化差異較大,且屬后驅(qū)動(dòng)因子,因此本文中僅以前驅(qū)動(dòng)因子降雨為依據(jù)來(lái)定性分析稻田徑流易發(fā)期(即徑流潛在發(fā)生風(fēng)險(xiǎn))。

    1.1.1 田間管理及統(tǒng)計(jì)區(qū)間劃分

    根據(jù)蘇州歷年水稻栽插實(shí)際,一般年份水稻栽插時(shí)間在6月5號(hào)左右開(kāi)始,6月25號(hào)栽插結(jié)束,按基肥期與栽插同步及施肥后7 d作為肥期,將基肥期區(qū)間劃分為6月6日至7月2日,共27日歷天。此外,一般栽插后10 d左右于水稻緩苗后追施一次分蘗肥,并將施肥后7 d時(shí)間劃分為蘗肥期,因此蘗肥期區(qū)間劃分為6月16日至7月12日,共27日歷天。太湖地區(qū)水稻孕穗期一般集中在7月底8月初,按生產(chǎn)實(shí)踐并考慮7 d肥期,將7月26日至8月12日劃分為穗肥期,共18日歷天。其余時(shí)間均為非肥期,并按生育進(jìn)程分為蘗肥-穗肥階段(7月13日至7月25日)和穗肥后至成熟階段(8月13日至11月10日)。需說(shuō)明的是,由于地區(qū)水稻栽插時(shí)間歷時(shí)較長(zhǎng),因此在統(tǒng)計(jì)區(qū)間上基肥和蘗肥期有時(shí)間重合,基肥-蘗肥階段無(wú)法統(tǒng)計(jì)體現(xiàn)。

    此外,該地區(qū)田間水分管理以干濕交替灌溉為主,前期栽插后保持5~8 cm淺水層促進(jìn)水稻活棵并抑制田間雜草發(fā)生,施肥前和孕穗揚(yáng)花期田間灌水保持5~8 cm水層,夠苗期曬田,其他時(shí)間灌水后至自然落干再補(bǔ)充水分。大面積田塊田埂高度一般在20 cm左右,但水稻生長(zhǎng)前期為了防止淹苗,田埂排水口高度一般控制在10 cm左右,水層較淺,拔節(jié)后水稻株高增加,排水口高度和田埂高度一致。

    1.1.2 統(tǒng)計(jì)分析方法

    將1956—2015年60年歷史氣象資料(中國(guó)氣象局氣象數(shù)據(jù)中心)按上述統(tǒng)計(jì)區(qū)間分別統(tǒng)計(jì)區(qū)間內(nèi)降雨頻次(天數(shù))、累積降雨量,并折算降雨概率(降雨頻次÷區(qū)間天數(shù))和單次平均降雨量(累積降雨量÷降雨頻次)。

    1.2 田間試驗(yàn)

    1.2.1 試驗(yàn)設(shè)計(jì)

    試驗(yàn)于2009年稻季開(kāi)始在太湖東岸蘇州市相城區(qū)望亭鎮(zhèn)迎湖村進(jìn)行,供試土壤類型為黃泥土。稻季采用大田育苗移栽,水稻收獲后作冬閑處理。試驗(yàn)設(shè)置5個(gè)處理,3次重復(fù),隨機(jī)區(qū)組排列,小區(qū)面積40 m2,小區(qū)間作埂隔離防滲。每個(gè)小區(qū)單設(shè)進(jìn)、排水口。水稻栽插規(guī)格為30 cm×13.3 cm,每穴定植2株。近3年(2014—2016年)水稻栽插及施基肥、施分蘗肥和施穗肥時(shí)間分別為2014-06-24、2014-07-10和2014-08-10,2015-06-21、2015-07-01和2015-08-13,2016-06-21、2016-07-06和2016-08-10。田間水分管理同當(dāng)?shù)匾恢?。按高產(chǎn)栽培要求進(jìn)行病蟲(chóng)草害防控。試驗(yàn)處理如下.

    當(dāng)?shù)爻R?guī)施肥處理(CN).總氮量270 kg·hm-2,分次施用,基肥∶分蘗肥∶穗肥分配比例為3∶3∶4。磷肥用量為80 kg·hm-2,一次性基施,鉀肥用量90 kg·hm-2,分兩次使用,基施50%、穗肥施用50%。

    化肥減量?jī)?yōu)化處理(RF).總氮量200 kg·hm-2,分次施用,基肥∶分蘗肥∶穗肥分配比例為3∶3∶4。磷肥用量為0 kg·hm-2,鉀肥用量90 kg·hm-2,分兩次使用,基施50%、穗肥施用50%。

    緩控釋肥減量替代處理(SCU).總氮量200 kg· hm-2,70%的氮采用水稻專用硫包衣尿素(N37%,漢楓集團(tuán)),基施,剩余30%氮采用尿素作為穗肥追施。磷肥用量為65 kg·hm-2,一次性基施,鉀肥用量90 kg· hm-2,分兩次使用,基施50%、穗肥施用50%。

    有機(jī)無(wú)機(jī)減量配施處理(OCN).總氮量200 kg· hm-2,其中有機(jī)肥占20%,采用商品有機(jī)肥(基施),其余80%為尿素?;省梅痔Y肥∶穗肥分配比例為3∶3∶4。磷肥用量為65 kg·hm-2,一次性基施,鉀肥用量90 kg·hm-2,分兩次使用,基施50%、穗肥施用50%。

    全有機(jī)肥減量替代處理(OF).總氮量200 kg· hm-2,采用商品有機(jī)肥,一次性基施。

    1.2.2 采樣及測(cè)定方法

    田面水.分別于基、追肥后一周時(shí)間內(nèi),每天連續(xù)采集各田塊田面水樣,采集部位按S型每小區(qū)采集3個(gè)點(diǎn)位樣品混合;采集后的混合水樣過(guò)濾后用AA3流動(dòng)分析儀分析水樣的總氮、銨態(tài)氮和硝態(tài)氮含量。

    徑流樣品.每次降雨產(chǎn)流后記錄同期降雨量并采集各小區(qū)預(yù)埋徑流收集管內(nèi)的徑流水樣,隨即帶回實(shí)驗(yàn)室過(guò)濾后利用AA3流動(dòng)分析儀分析徑流水樣的總氮、銨態(tài)氮和硝態(tài)氮含量。

    1.2.3 統(tǒng)計(jì)分析方法

    按1.1區(qū)間劃分方法,對(duì)徑流易發(fā)期不同區(qū)間3年發(fā)生徑流進(jìn)行均值計(jì)算,來(lái)評(píng)價(jià)不同肥料類型該期徑流發(fā)生風(fēng)險(xiǎn)。常規(guī)數(shù)據(jù)處理和作圖在Excel 2007中進(jìn)行,方差分析在SPSS 16.0中采用Duncan法對(duì)數(shù)據(jù)進(jìn)行差異顯著性分析(P<0.05)。

    2 結(jié)果與分析

    2.1 徑流易發(fā)期頻次統(tǒng)計(jì)

    2.1.1 1956—2015年歷史年降雨量和稻季降雨量

    圖1為1956—2015年60年歷史降雨數(shù)據(jù)統(tǒng)計(jì)分析。結(jié)果表明,近60年太湖地區(qū)年降雨量在680.1~1 699.7 mm之間,平均年降雨量為1 132.6 mm。水稻季(6月6日至11月10日)降雨量約為215.3~1 285.1 mm,平均約為629.1 mm。其中稻季降雨最低和最高年份分別出現(xiàn)在1967年和1999年。稻季降雨占全年降雨量的比例平均為54.68%(27.94%~75.61%),占比較高。從近10年降雨數(shù)據(jù)來(lái)看,年均降雨量約為1 182.3 mm(914~1364 mm),波動(dòng)較小。稻季降雨量均值和占比分別為650.8 mm和55.17%。近10年降雨統(tǒng)計(jì)數(shù)據(jù)與60年統(tǒng)計(jì)數(shù)據(jù)均值接近,說(shuō)明利用近60年降雨數(shù)據(jù)預(yù)測(cè)和指征該區(qū)域降雨具有一定代表性。

    2.1.2 稻季不同肥期區(qū)間降雨概率

    近60年水稻生長(zhǎng)季單日降雨概率列于圖2。結(jié)果可以看出,水稻生長(zhǎng)前期降雨概率較大,除個(gè)別天次,6月21日—7月6日單日降雨概率均超過(guò)50%。從肥期劃分統(tǒng)計(jì)結(jié)果來(lái)看,基肥期(BF)和蘗肥期(TF)降雨概率明顯高于其他時(shí)期,是稻田徑流的易發(fā)期,降雨概率分別達(dá)到48.15%和49.81%。蘗肥-穗肥階段(T-H)、穗肥期(HF)和穗肥后至成熟階段(HM)降雨概率較小,分別為35.13%、30.28%和32.74%。五個(gè)時(shí)期降雨概率表現(xiàn)為.TF>BF>T-H>H-M>HF。

    2.1.3 稻季不同肥期區(qū)間單次平均降雨量

    圖1 近60年蘇州地區(qū)年降雨量和稻季降雨量(1956—2015)Figure 1 Annual and rice season rainfall in Suzhou in recent 60 years(1956—2015)

    近60年日均降雨量數(shù)據(jù)表明(圖3),與降雨概率相一致,基肥期和蘗肥期日均降雨量較高,分別為12.81 mm和12.84 mm。此外,穗肥期日均降雨量也較高,為11.73 mm。蘗肥-穗肥階段和穗肥后-成熟階段降雨量較低,分別為10.11 mm和9.43 mm。五個(gè)時(shí)期日均降雨量整體表現(xiàn)為.TF>BF>HF>T-H>H-M。

    2.2 徑流易發(fā)期實(shí)際發(fā)生徑流

    2.2.1 徑流易發(fā)期實(shí)際徑流發(fā)生頻次

    圖2 水稻生長(zhǎng)季單日降雨概率統(tǒng)計(jì)(1956—2015)Figure 2 Day rainfall probability in rice growing season(1956—2015)

    圖3 水稻生長(zhǎng)季有降雨天次日平均降雨量(1956—2015)Figure 3 The average rainfall for rainfall day in rice growing season(1956—2015)

    前述分析表明,基肥期和蘗肥期(含基肥-蘗肥階段)是徑流發(fā)生的高風(fēng)險(xiǎn)期和徑流減排的關(guān)鍵時(shí)期。表1為長(zhǎng)期定位監(jiān)測(cè)試驗(yàn)2014、2015年和2016年在基肥期、蘗肥期和基肥至蘗肥階段三個(gè)徑流易發(fā)期實(shí)際發(fā)生徑流頻次的統(tǒng)計(jì)結(jié)果。從結(jié)果可以看出,近3年徑流易發(fā)期共有9次徑流事件發(fā)生,降雨量集中在11.5~282.3 mm。其中,基肥期共有4次徑流發(fā)生,其中2014年1次(降雨量80 mm),2015年1次(降雨量24.3 mm),2016年2次(降雨量分別為70.2 mm和157.5 mm);蘗肥期共有2次徑流發(fā)生,分別發(fā)生在2014年(降雨量73.9mm)和2016年(降雨量11.5mm);基肥至蘗肥階段共有3次徑流發(fā)生,3年各產(chǎn)生1次徑流事件(降雨量分別為19、31.3和282.3 mm)。

    表1 徑流易發(fā)期實(shí)際發(fā)生徑流頻次及同期降雨量Table 1 The actual frequency of runoff and precipitation in the same period in runoff prone period

    2.2.2 徑流氮素濃度

    表2為三年監(jiān)測(cè)期實(shí)際徑流發(fā)生不同形態(tài)氮素質(zhì)量濃度。為便于統(tǒng)計(jì)分析和分期比較,對(duì)同一時(shí)期徑流事件實(shí)測(cè)氮素質(zhì)量濃度進(jìn)行了均值計(jì)算(表2)。結(jié)果表明,稻田徑流易發(fā)期徑流氮素?fù)p失以銨態(tài)氮為主,基肥期和蘗肥期的徑流總氮和銨態(tài)氮質(zhì)量濃度明顯高于基肥-蘗肥階段。從總氮質(zhì)量濃度結(jié)果來(lái)看,與氮素常規(guī)用量CN處理相比,年際間變化較一致,化肥減量?jī)?yōu)化處理RF和有機(jī)肥部分替代處理OCN在三個(gè)時(shí)期的徑流氮素濃度均較低,RF分別減少了10.64%、2.49%和13.36%,OCN分別減少了18.61%、14.71%和24.22%,均減少了徑流損失的風(fēng)險(xiǎn)。硫包衣尿素減量替代處理SCU基施后的基肥期和基肥-蘗肥階段的徑流氮素濃度較高,2014年和2016年基肥期總氮濃度均顯著高于當(dāng)?shù)爻R?guī)施肥處理CN,增加了氮素徑流流失風(fēng)險(xiǎn)(基肥期增加20.3%、基肥-蘗肥階段增加11.72%),但蘗肥期氮素濃度較低,其中2014年差異達(dá)到顯著水平,降低了蘗肥期氮素流失風(fēng)險(xiǎn)(減少30.72%)。而全有機(jī)肥減量替代處理OF肥期的徑流氮素濃度較低(基肥期降低9.04%、蘗肥期降低28.53%),但非肥期(基肥-蘗肥階段)徑流氮素濃度較高,增加了19.7%。與無(wú)機(jī)化肥減量?jī)?yōu)化處理RF相比,相同氮素用量下,SCU和OF處理均表現(xiàn)為基肥期和基肥-蘗肥階段徑流氮素濃度較高,OCN處理徑流氮素濃度較低。徑流銨態(tài)氮濃度變化趨勢(shì)與總氮濃度一致。不同階段硝態(tài)氮濃度年際間變化差異較大,其中基肥期2016年6月24發(fā)生的徑流硝態(tài)氮濃度明顯高于其他3次,基肥-蘗肥階段的2014年7月2日徑流硝態(tài)氮濃度明顯高于其他2次,這可能與干濕交替水分管理有關(guān)。處理間硝態(tài)氮濃度比較來(lái)看,當(dāng)?shù)爻R?guī)化肥處理的硝態(tài)氮濃度三個(gè)時(shí)期均較高, 2014年6月26(基肥期)和2014年7月2日(基肥-蘗肥階段)顯著高于OCN處理,2016年6月28(基肥期)顯著高于SCU處理,2014年7月13(蘗肥期)顯著高于OF處理,其他處理間差異均未達(dá)到顯著水平。

    表2 徑流易發(fā)期實(shí)際發(fā)生徑流氮素質(zhì)量濃度/mg·L-1Table 2 Nitrogen concentration in runoff water in runoff prone period

    2.3 基肥期和蘗肥期田面水氮素濃度

    太湖地區(qū)稻季田面水氮素徑流損失是太湖水體富營(yíng)養(yǎng)化的重要原因。因此,田面水氮素質(zhì)量濃度在一定程度上可以反映稻田氮素徑流損失風(fēng)險(xiǎn)。為進(jìn)一步實(shí)證不同減量替代處理在徑流易發(fā)期徑流流失風(fēng)險(xiǎn),對(duì)定位試驗(yàn)稻田施肥后基肥期和蘗肥期田面水氮素質(zhì)量濃度進(jìn)行了分析(表3)。結(jié)果表明,稻田田面水氮素以銨態(tài)氮為主,硝態(tài)氮濃度較低。田面水銨態(tài)氮濃度變化與總氮濃度變化一致。從總氮濃度結(jié)果來(lái)看,不同處理下田面水氮素濃度均在施肥后1 d達(dá)到峰值,此后逐漸降低。從田面水氮素濃度均值來(lái)看,除基肥期SCU處理的氮素濃度高于CN處理,其他處理基肥期和蘗肥期氮素濃度均低于CN處理。此外,從基肥期SCU處理的單日田面水氮素濃度變化來(lái)看, SCU處理的田面水峰值濃度低于CN處理,但后期CN處理濃度下降較快。

    表3 不同施肥處理基肥期和蘗肥期田面水氮素質(zhì)量濃度(2016)(mg·L-1)Table 3 Nitrogen concentration in surface water in basal fertilization and tillering fertilization stages under different fertilization treatments(mg·L-1)

    與減量?jī)?yōu)化施肥處理RF相比,等氮量的SCU處理的基肥期田面水氮素濃度較高,但蘗肥期較低。全有機(jī)肥替代OF處理的田面水氮素濃度在基肥期和蘗肥期均低于RF處理。而有機(jī)肥部分替代處理OCN則表現(xiàn)為基肥期較低,但蘗肥期較高,這可能與蘗肥期無(wú)機(jī)氮素投入相同,但基肥期施入的有機(jī)肥緩慢釋放有關(guān)。

    田面水硝態(tài)氮濃度結(jié)果表明,與CN處理相比,基肥期施用有機(jī)肥處理(OCN和OF)的硝態(tài)氮濃度較高,減量?jī)?yōu)化施肥處理RF和硫包衣尿素替代處理SCU的硝態(tài)氮濃度較低,蘗肥期施用無(wú)機(jī)尿素處理(CN、RF和OCN)的硝態(tài)氮濃度較高,CN處理的濃度最高,而在該期不施肥的SCU和OF處理硝態(tài)氮濃度較低,但處理間差異未達(dá)到顯著水平。

    3 討論

    對(duì)于稻田徑流發(fā)生流失風(fēng)險(xiǎn)的研究,多集中于徑流濃度和損失量的比較,對(duì)于不同生育期發(fā)生風(fēng)險(xiǎn)的關(guān)注較少[10-13]。稻田田面在田埂保護(hù)下形成封閉徑流體系,只有在特殊情況如降雨發(fā)生時(shí),田面水才會(huì)溢出形成機(jī)會(huì)徑流[6]。有研究指出,增加田埂排水口高度是減少稻田氮素流失的重要途徑,但水稻生長(zhǎng)前期為了防止淹苗,田埂一般較低,水層較淺,這更增加了水稻生長(zhǎng)前期的流失風(fēng)險(xiǎn)[14,19]。此外,焦瑞峰等[18]研究發(fā)現(xiàn),水稻生長(zhǎng)后期水稻冠層對(duì)于降雨能量有一定的緩沖作用,這在一定程度上減緩了徑流發(fā)生的風(fēng)險(xiǎn)。本研究結(jié)果表明,水稻生長(zhǎng)前期降雨概率較大,6月21日—7月6日單日降雨概率均超過(guò)50%(圖2)。從階段統(tǒng)計(jì)結(jié)果來(lái)看,基肥期和蘗肥期(含基肥-蘗肥階段)降雨概率和降雨量均明顯高于其他階段,降雨概率分別達(dá)48.15%和49.81%,降雨量分別為12.81 mm和12.84 mm,均超過(guò)實(shí)際監(jiān)測(cè)到徑流的同期最低降雨量11.5 mm(表1),產(chǎn)生徑流可能性較大。這說(shuō)明,由降雨驅(qū)動(dòng)的水稻生長(zhǎng)前期降雨概率和降雨量本身高于其他生育期。而當(dāng)?shù)氐那捌谒畬庸芾泶胧?栽插后保持5~8 cm淺水層促進(jìn)水稻活棵并抑制田間雜草發(fā)生,此外前期為了防止淹苗,田埂排水口高度一般控制在10 cm左右)及較小的水稻冠層葉面積加重了徑流發(fā)生的風(fēng)險(xiǎn)。此外,水稻高產(chǎn)氮肥運(yùn)籌中基肥和蘗肥的氮肥投入約占施氮總量的50%~60%[20-22],這更增加了前期氮素流失風(fēng)險(xiǎn)。

    新型緩控釋肥和有機(jī)肥等替代肥料的養(yǎng)分釋放特性與無(wú)機(jī)化肥差異較大[14]。尿素等化肥施入稻田后迅速溶于水,短時(shí)間內(nèi)田面水氮素濃度升至峰值,有機(jī)肥由于本身成分的復(fù)雜性及受微生物作用,其氮素釋放通常呈現(xiàn)一個(gè)長(zhǎng)期的波動(dòng)過(guò)程,而緩控釋肥料由于肥料的緩慢釋放,田面水氮素濃度變化也呈現(xiàn)出長(zhǎng)期的變化過(guò)程[9,15-17]。為了明確不同肥料類型及施用方式在徑流易發(fā)期徑流流失風(fēng)險(xiǎn),本研究利用長(zhǎng)期定位試驗(yàn)從實(shí)際觀測(cè)徑流和施肥后田面水氮素質(zhì)量濃度兩方面進(jìn)行了分析。結(jié)果表明,徑流和田面水氮素均以銨態(tài)氮為主,硝態(tài)氮濃度較低,不同年際間徑流氮素濃度變化處理間較一致?;蕼p量處理較常規(guī)用量處理的徑流和田面水氮素濃度均降低,徑流易發(fā)期3個(gè)階段分別較CN處理降低了10.64%、2.49%和13.36%,有利于降低徑流氮素流失。而硫包衣尿素由于一次施用較大,基肥期和基肥-蘗肥階段實(shí)測(cè)徑流氮素濃度和基肥期田面水氮素濃度明顯高于其他處理(基肥期較CN增加了20.3%、基肥-蘗肥階段較CN增加了11.72),增加了氮素徑流流失風(fēng)險(xiǎn),但蘗肥期徑流和田面水氮素濃度均較低,徑流總氮濃度較CN降低了30.72%,流失風(fēng)險(xiǎn)較小。這可能與硫包衣尿素的釋放特性有關(guān)。谷佳林等[23]對(duì)硫包衣尿素肥料釋放特性研究發(fā)現(xiàn),硫包衣尿素養(yǎng)分釋放速率受溫度和水分影響較大,肥料所處環(huán)境溫度越高,水分含量越大,養(yǎng)分釋放速率越快。稻田溫度較高,移栽時(shí)土壤水分含量處于飽和狀態(tài),這加快了硫包衣尿素的釋放[24-25]。有機(jī)肥氮素分解是一個(gè)復(fù)雜的過(guò)程,受其化學(xué)成分和土壤溫濕度、通氣狀況、C/N等多因素影響[26-27]。本研究結(jié)果表明,有機(jī)肥部分替代(20%替代)和全量替代處理徑流易發(fā)期氮素流失風(fēng)險(xiǎn)表現(xiàn)不一致,基肥期和基肥-蘗肥階段OF處理的氮素濃度較高,蘗肥期OCN處理的徑流和田面水氮素濃度較高,這可能是由于添加量和施肥方式差異導(dǎo)致的。此外,與無(wú)機(jī)化肥處理相比,有機(jī)肥全量替代處理的肥期(基肥期和蘗肥期)徑流氮素濃度和田面水氮素濃度均低于高量和等量無(wú)機(jī)化肥處理(CN、RF),較CN處理基肥期降低了9.04%、蘗肥期降低了28.53%,但非肥期(基肥-蘗肥階段)徑流氮素濃度較高,較CN增加了19.7%,氮素徑流損失增加。與有機(jī)肥全量添加不同的是,20%有機(jī)肥添加處理OCN的徑流易發(fā)期的徑流氮素濃度低于高量和等量無(wú)機(jī)化肥處理(CN、RF),三個(gè)階段分別較CN處理減少了18.61%、14.71%和24.22%,但等氮量下蘗肥期田面水峰值濃度較無(wú)機(jī)化肥RF處理高,這可能與有機(jī)肥的釋放密切相關(guān)。與單施有機(jī)肥相比,配施化肥可以降低初始碳氮比,并為異養(yǎng)微生物提供充足的無(wú)機(jī)氮源,促進(jìn)有機(jī)肥氮素礦化,增加氮素礦化量[28-30]。此外,從徑流硝態(tài)氮濃度結(jié)果來(lái)看,除CN處理的徑流硝態(tài)氮濃度較高外,其他處理均不顯著,但值得注意的是,不同階段硝態(tài)氮濃度年際間變化差異較大,其中基肥期2016年6月24日發(fā)生的徑流硝態(tài)氮濃度明顯高于其他3次,基肥-蘗肥階段的2014年7月2日徑流硝態(tài)氮濃度明顯高于其他2次。有研究表明,在一定土壤水分含量下,硝化速率隨水分含量的增加而增加,在氧供應(yīng)受限時(shí),硝化速率開(kāi)始下降[31]。因此,不同年際間硝態(tài)氮濃度的變化可能與干濕交替水分管理導(dǎo)致的土壤含水量及氧化還原電位變化有關(guān)。由于田間徑流監(jiān)測(cè)僅采用徑流收集管對(duì)徑流氮素濃度進(jìn)行了評(píng)估,因此對(duì)不同時(shí)期徑流損失量未進(jìn)行量化比較,這需要在今后試驗(yàn)中進(jìn)一步補(bǔ)充量化,以此對(duì)不同階段氮素實(shí)際損失量進(jìn)行比較分析。

    4 結(jié)論

    (1)基肥期和蘗肥期(含基肥-蘗肥階段)降雨概率和降雨量較高,降雨概率分別達(dá)48.15%和49.81%,降雨量分別為12.81 mm和12.84 mm,均超過(guò)實(shí)際監(jiān)測(cè)到徑流的同期最低降雨量11.5 mm,產(chǎn)生徑流可能性較大,是稻田徑流的易發(fā)期。

    (2)不同類型肥料能夠降低易發(fā)期內(nèi)不同階段徑流氮素濃度,但在徑流易發(fā)期的徑流氮素?fù)p失控制效果不能一概而論。與常規(guī)化肥用量處理相比,化肥減量和有機(jī)肥減量替代處理能夠降低徑流易發(fā)期徑流氮素?fù)p失;但硫包衣尿素減量替代處理的基肥期和基肥-蘗肥階段徑流氮素?fù)p失風(fēng)險(xiǎn)增加,但蘗肥期徑流氮素濃度減少;而全有機(jī)肥減量替代處理的肥期徑流和田面水氮素濃度降低,其中基肥期徑流氮素濃度較CN處理降低9.04%、蘗肥期降低28.53%,但基肥-蘗肥階段徑流氮素濃度較CN處理增加了19.7%,氮素徑流損失風(fēng)險(xiǎn)較高。

    [1]楊林章,施衛(wèi)明,薛利紅,等.農(nóng)村面源污染治理的“4R”理論與工程實(shí)踐——總體思路與“4R”治理技術(shù)[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2013, 32(1).1-8.

    YANG Lin-zhang,SHI Wei-ming,XUE Li-hong,et al.Reduce-retainreuse-restore technology for the controlling the agricultural non-point source pollution in countryside in China.General countermeasures and technologies[J].Journalof Agro-Environment Science,2013,32(1).1-8.

    [2]Yang S,Peng S,Xu J,et al.Nitrogen loss from paddy field with different water and nitrogen managements in Taihu Lake region of China[J].Communications in Soil Science and Plant Analysis,2013,44(16).2393-2407.

    [3]Zhao X,Zhou Y,Min J,et al.Nitrogen runoff dominates water nitrogen pollution from rice-wheat rotation in the Taihu Lake region of China[J].Agriculture Ecosystems&Environment,2012,156(4).1-11.

    [4]王桂良,崔振嶺,陳新平,等.南方稻田活性氮損失途徑及其影響因素[J].應(yīng)用生態(tài)學(xué)報(bào),2015,26(8).2337-2345.

    WANG Gui-liang,CUI Zhen-ling,CHEN Xin-ping,et al.Reactive nitrogen loss pathways and their effective factors in paddy field in Southern China[J].Chinese Journal of Applied Ecology,2015,26(8).2337-2345.

    [5]薛利紅,俞映倞,楊林章.太湖流域稻田不同氮肥管理模式下的氮素平衡特征及環(huán)境效應(yīng)評(píng)價(jià)[J].環(huán)境科學(xué),2011,32(4).1133-1138.

    XUE Li-hong,YU Ying-liang,YANG Lin-zhang.Nitrogen balance and environmental impact of paddy field under different N management methods in Taihu Lake region[J].Environmental Science,2011,32(4). 1133-1138.

    [6]薛利紅,楊林章,施衛(wèi)明,等.農(nóng)村面源污染治理的“4R”理論與工程實(shí)踐——源頭減量技術(shù)[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2013,32(5).881-888. XUE Li-hong,YANG Lin-zhang,SHI Wei-ming,et al.Reduce-retainreuse-restore technology for controlling the agricultural non-point pollution in countryside in China.Source reduction technology[J].Journal of Agro-Environment Science,2013,32(5).881-888.

    [7]Ye Y S,Liang X Q,Chen Y X,et al.Alternate wetting and drying irrigation and controlled-release nitrogen fertilizer in late-season rice.Effects on dry matter accumulation,yield,water and nitrogen use[J].Field Crops Research,2013,144(6).212-224.

    [8]侯朋福,薛利祥,范立慧,等.氮肥類型及施用方式對(duì)節(jié)水抗旱稻田面水氮?jiǎng)討B(tài)及產(chǎn)量的影響[J].中國(guó)稻米,2016,22(3).43-48.

    HOU Peng-fu,XUE Li-xiang,FAN Li-hui,et al.Effects of nitrogen fertilizer type and application method on surface water nitrogen dynamics and yield of drought-resistance rice[J].China Rice,2016,22(3).43-48.

    [9]Shaviv A,Mikkelsen R L.Controlled-release fertilizers to increase efficiency of nutrient use and minimize environmental degradation.A review[J].Nutrient Cyclingin Agroecosystems,1993,35(1).1-12.

    [10]Xue L H,Yu Y L,Yang L Z.Maintaining yields and reducing nitrogen loss in rice-wheat rotation system in Taihu Lake region with proper fertilizer management[J].EnvironmentalResearch Letters,2014,9(11). 115010.

    [11]Yu Y L,Xue L H,Yang L Z.Winter legumes in rice crop rotations reduces nitrogen loss,and improves rice yield and soil nitrogen supply[J].Agronomy for Sustainable Development,2014,34(3).633-640.

    [12]Qiao J,Yang L Z,Yan T M,et al.Nitrogen fertilizer reduction in riceproduction for two consecutive years in the Taihu Lake area[J].Agriculture Ecosystems&Environment,2012,146(1).103-112.

    [13]紀(jì)雄輝,鄭圣先,魯艷紅,等.施用尿素和控釋氮肥的雙季稻田表層水氮素動(dòng)態(tài)及其徑流損失規(guī)律[J].中國(guó)農(nóng)業(yè)科學(xué),2006,39(12). 2521-2530.

    JI Xiong-hui,ZHENG Sheng-xian,LU Yan-hong,et al.Dynamics of floodwater nitrogen and its runoff loss,urea and controlled release nitrogen fertilizer application regulation in rice[J].Scientia Agricultura Sinica,2006,39(12).2521-2530.

    [14]傅朝棟.基于ArcGIS的區(qū)域尺度稻田降雨徑流氮磷流失負(fù)荷估算[D].杭州.浙江大學(xué),2015.

    FU Chao-dong.ArcGIS based regional estimation of nitrogen and phosphorus loss loads in surface runoff from rice paddy fields[D]. Hangzhou.Zhejiang University,2015.

    [15]郭智,周煒,陳留根,等.施用豬糞有機(jī)肥對(duì)稻麥兩熟農(nóng)田稻季養(yǎng)分徑流流失的影響[J].水土保持學(xué)報(bào),2013,27(6).21-25.

    GUO Zhi,ZHOU Wei,CHEN Liu-gen,et al.Effect of pig manure application on surface runoff losses of soil nitrogen and phosphorus during the paddy season in intersive rice-wheat rotation field[J].Journal of Soil and Water Conservation,2013,27(6).21-25.

    [16]Guo H Y,Zhu J G,Wang X R,et al.Case study on nitrogen and phosphorus emissions from paddy field in Taihu region[J].Environmental Geochemistry and Health,2004,26(2).209-219.

    [17]金潔,楊京平,施洪鑫,等.水稻田面水中氮磷素的動(dòng)態(tài)特征研究[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2005,24(2).357-361.

    JIN Jie,YANG Jing-ping,SHI Hong-xin,et al.Variations of nitrogen and phosphorus in surface water body of a paddy field[J].Journal of A-gro-Environment Science,2005,24(2).357-361.

    [18]焦瑞鋒,張鴻睿,卞新民.稻田降雨徑流擾動(dòng)系數(shù)的解析分析[J].浙江農(nóng)業(yè)學(xué)報(bào),2012,24(5).898-903.

    JIAO Rui-feng,ZHANG Hong-rui,BIAN Xin-min.Analytical analysis of paddy runoff disturbance coefficient[J].Acta Agriculturae Zhejiangensis,2012,24(5).898-903.

    [19]Mishra A,Ghorai A K,Singh S R.Rainwater,soil and nutrient conservation in rainfed rice lands in Eastern India[J].Agricultural Water Management,1998,38(1).45-57.

    [20]薛利紅,李剛?cè)A,侯朋福,等.太湖地區(qū)稻田持續(xù)高產(chǎn)的減量施氮技術(shù)體系研究[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2016,35(4).729-736.

    XUE Li-hong,LI Gang-hua,HOU Peng-fu,et al.Nitrogen reduction technique system for sustaining high yield of paddy fields[J].Journal of Agro-Environment Science,2016,35(4).729-736.

    [21]林晶晶.水稻對(duì)基蘗穗肥氮素的吸收利用差異及適宜比例的研究[D].南京.南京農(nóng)業(yè)大學(xué),2014.

    LIN Jing-jing.Studies on absorption and utilization and optimum proportion of nitrogen of basic,tillering and panicle fertilizer in rice[D]. Nanjing.Nanjing Agricultural University,2014.

    [22]范立慧,徐珊珊,侯朋福,等.不同地力下基蘗肥運(yùn)籌比例對(duì)水稻產(chǎn)量及氮肥吸收利用的影響[J].中國(guó)農(nóng)業(yè)科學(xué),2016,49(10). 1872-1884.

    FAN Li-hui,XU Shan-shan,HOU Peng-fu,et al.Effect of different ratios of basal to tiller nitrogen on rice yield and nitrogen utilization under different soil fertility[J].Scientia Agricultura Sinica,2016,49 (10).1872-1884.

    [23]谷佳林,徐凱,付鐵梅,等.不同密閉材料硫包衣尿素氮素釋放特性及對(duì)夏玉米生長(zhǎng)的影響[J].植物營(yíng)養(yǎng)與肥料學(xué)報(bào),2011,17(3). 630-637.

    GU Jia-lin,XU Kai,FU Tie-mei,et al.Nitrogen release characteristics of different hermetic material sulfur coated urea and their effects on summer maize[J].Plant Nutrition and Fertilizer Science,2011,17(3). 630-637.

    [24]Miao X K,Xing X M,Ding Y F,et al.Yield and nitrogen uptake of bowl-seedling machine-transplanted rice with slow-release nitrogen fertilizer[J].Agronomy Journal,2015,108(1).313-320.

    [25]李敏,郭熙盛,葉舒婭,等.硫膜和樹(shù)脂膜控釋尿素對(duì)水稻產(chǎn)量、光合特性及氮肥利用率的影響[J].植物營(yíng)養(yǎng)與肥料學(xué)報(bào),2013,19 (4).808-815.

    LI Min,GUO Xi-sheng,YE Shu-ya,et al.Effects of sulfur-and polymer-coated controlled release urea on yield,photosynthetic characteristics and nitrogen fertilizer efficiency of rice[J].Journal of Plant Nutrition and Fertilizer,2013,19(4).808-815.

    [26]劉更另,金維續(xù).中國(guó)有機(jī)肥料[M].北京.農(nóng)業(yè)出版社,1991.

    LIU Geng-ling,JIN Wei-xu.Chinese organic fertilizer[M].Beijing.A-gricultural Press,1991.

    [27]朱大威,黃耀,盧燕宇.有機(jī)肥氮素釋放動(dòng)態(tài)模型的初步研究[J].南京農(nóng)業(yè)大學(xué)學(xué)報(bào),2006,29(3).146-150.

    ZHU Da-wei,HUANG Yao,LU Yan-yu.A primary study on modeling nitrogen release from organic manure[J].Journal of Nanjing Agricultural University,2006,29(3).146-150.

    [28]張璐,石林,蔡澤江,等.有機(jī)肥與化肥不同比例配施下水稻土銨態(tài)氮釋放特征[J].中國(guó)土壤與肥料,2015(4).15-22.

    ZHANG Lu,SHI Lin,CAI Ze-jiang,et al.The ammonium nitrogen release characteristic of paddy soil with mixed application of different ratios of organic and inorganic fertilizer[J].Soil and Fertilizer Sciences in China,2015(4).15-22.

    [29]Chantigny M H,Angers D A,Prévost D,et al.Dynamics of soluble organic C and C mineralization in cultivated soils with varying N fertilization[J].Soil Biology&Biochemistry,1999,31(4).543-550.

    [30]Sikora L J,Enkiri N K.Efficiency of compost-fertilizer blends compared with fertilizer alone[J].Soil Science,2000,165(5).444-451.

    [31]張威,張旭東,何紅波,等.干濕交替條件下土壤氮素轉(zhuǎn)化及其影響研究進(jìn)展[J].生態(tài)學(xué)雜志,2010,29(4).783-789.

    ZHANG Wei,ZHANG Xu-dong,HE Hong-bo,et al.Research advances in soil nitrogen transformation as related to drying/wetting cycles[J].Chinese Journal of Ecology,2010,29(4).783-789.

    Nitrogen loss risk of paddy field under different fertilizations in runoff prone period

    HOU Peng-fu1,XUE Li-xiang1,YU Ying-liang1,JIANG Yu2,WANG Ji-dong1,XUE Li-hong1*,YANG Lin-zhang1
    (1.Institute of Agricultural Resources and Environment,Jiangsu Academy of Agricultural Sciences,Nanjing 210014,China;2.Institute of Crop Science,Chinese Academy of Agricultural Sciences,Beijing 100081,China)

    The runoff prone period of paddy field in Taihu Lake region was analyzed according to the typical natural rainfall data of 60 years in this study.Then,the actual runoff nitrogen control effect from paddy field in runoff prone period of different fertilization reduction treatment were compared using the runoff monitoring data in recent 3 years of a long-term fertilization reduction experiment that initiated in 2009.The results showed that,the rainfall probability was larger in early rice growth stage in Taihu Lake region and daily rainfall probability between June 21st and July 6th were more than 50%.The rainfall probability of basal fertilization stage and tillering fertilization stage were 48.15%and 49.81%,respectively,obviously higher than other stages.The mean precipitation of these two stages were 12.81 mm and 12.84 mm,both exceed the minimum rainfall for runoff(11.5 mm),implying high probability for paddy runoff in these two periods.The nitrogen in runoff water and surface water in the runoff prone period were mainly ammonium nitrogen as the long-term experiment study indicated.Compared with conventional nitrogen rate nitrogen rate(CN)treatment,the nitrogen concentrations of reduced nitrogen rate treatment(RF)and organic combined with chemical nitrogen treatment(OCN)were decreased 8.83%and 19.18%,respectively.The nitrogenconcentrations in runoff water and surface water of sulfur coating urea treatment(SCU)in basal fertilization and basal fertilization-tillering fertilization stage were higher than those of other treatments,which increased 20.3%and 11.72%respectively incomparasion with CN treatment.The nitrogen concentration in runoff water of organic fertilizer treatment(OF)decreased 9.04%and 28.53%in basal fertilization and tillering fertilization stage,respectively,but increased 19.7%in tillering fertilization-panicle fertilization stage compared with CN treatment,which increased the risk of runoff nitrogen loss.The nitrogen reduction measures referred in this study can reduce the runoff water nitrogen concentrations in the runoff prone period of paddy field,but the effect on nitrogen loss control in the runoff prone period are different.

    paddy field;runoff;prone period;fertilizer type;nitrogen loss

    X592

    A

    1672-2043(2017)07-1353-09

    10.11654/jaes.2017-0695

    侯朋福,薛利祥,俞映倞,等.稻田徑流易發(fā)期不同類型肥料的氮素流失風(fēng)險(xiǎn)[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2017,36(7).1353-1361.

    HOU Peng-fu,XUE Li-xiang,YU Ying-liang,,et al.Nitrogen loss risk of paddy field under different fertilizations in runoff prone period[J].Journal of Agro-Environment Science,2017,36(7).1353-1361.

    2017-05-12

    侯朋福(1986—),男,山東高密人,博士,助理研究員,主要從事農(nóng)田養(yǎng)分管理與面源污染控制研究。E-mail:pengfuhou100smby@163.com

    *通信作者:薛利紅E-mail:njxuelihong@gmail.com

    公益性行業(yè)(農(nóng)業(yè))科研專項(xiàng)(201503106);國(guó)家自然科學(xué)基金青年科學(xué)基金項(xiàng)目(41601319);國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2017YFD0300100)

    Project supported:The Special Fund for Agro-scientific Research on Public Interests(201503106);The Young Scientists Fund of the National Natural Science Foundation of China(41601319);The National Key Research and Development Program China(2017YFD0300100)

    猜你喜歡
    田面穗肥易發(fā)
    長(zhǎng)期秸稈還田對(duì)水稻產(chǎn)量與田面水環(huán)境的影響
    機(jī)用鎳鈦銼在乳磨牙根管治療中的應(yīng)用
    貴州省地質(zhì)災(zāi)害易發(fā)分區(qū)圖
    夏季羊易發(fā)疾病及防治方法
    春耕稻田滯水減排控制面源污染效果研究
    冬季雞腸炎易發(fā) 科學(xué)防治有方法
    不同穗肥用量和運(yùn)籌對(duì)小麥生長(zhǎng)的影響
    摻混控釋肥側(cè)深施對(duì)稻田田面水氮素濃度的影響
    水稻全程機(jī)械化灌溉技術(shù)模式應(yīng)用
    穗肥施用時(shí)期對(duì)喀斯特區(qū)域雜交稻產(chǎn)量及農(nóng)藝性狀的影響
    国产不卡一卡二| av福利片在线观看| 中文字幕av在线有码专区| 午夜日韩欧美国产| 欧美激情久久久久久爽电影| 亚洲熟女毛片儿| 神马国产精品三级电影在线观看| 日韩欧美一区二区三区在线观看| 俄罗斯特黄特色一大片| 禁无遮挡网站| www.精华液| 91麻豆av在线| 在线国产一区二区在线| 亚洲电影在线观看av| 最好的美女福利视频网| 变态另类丝袜制服| 中亚洲国语对白在线视频| 熟女人妻精品中文字幕| 日韩欧美国产在线观看| 国产精品99久久久久久久久| 1000部很黄的大片| 丰满人妻一区二区三区视频av | 身体一侧抽搐| 丁香欧美五月| 国产免费av片在线观看野外av| 免费看a级黄色片| av在线蜜桃| 一二三四在线观看免费中文在| 成人欧美大片| 丁香六月欧美| 午夜久久久久精精品| 欧美色欧美亚洲另类二区| 成年人黄色毛片网站| 欧美激情在线99| 欧美性猛交黑人性爽| 国产91精品成人一区二区三区| 成人国产一区最新在线观看| 此物有八面人人有两片| 久久欧美精品欧美久久欧美| 亚洲精品在线观看二区| 999久久久国产精品视频| 精品不卡国产一区二区三区| 亚洲国产精品合色在线| 日韩 欧美 亚洲 中文字幕| 神马国产精品三级电影在线观看| 噜噜噜噜噜久久久久久91| 亚洲国产色片| 一级毛片精品| 伊人久久大香线蕉亚洲五| 波多野结衣高清作品| 国内精品美女久久久久久| 无遮挡黄片免费观看| 亚洲av电影在线进入| 亚洲色图 男人天堂 中文字幕| 丰满人妻熟妇乱又伦精品不卡| 99久久99久久久精品蜜桃| 每晚都被弄得嗷嗷叫到高潮| 国产精品,欧美在线| 亚洲黑人精品在线| 99热6这里只有精品| 啦啦啦观看免费观看视频高清| 亚洲熟女毛片儿| 日韩欧美一区二区三区在线观看| 亚洲国产欧美一区二区综合| 久久热在线av| ponron亚洲| 亚洲熟妇中文字幕五十中出| 成人av在线播放网站| 亚洲人成网站高清观看| 国产亚洲av高清不卡| 草草在线视频免费看| 两个人视频免费观看高清| 日本黄色片子视频| 欧美激情久久久久久爽电影| av欧美777| 亚洲av第一区精品v没综合| 亚洲欧美激情综合另类| 亚洲国产欧美网| 国内少妇人妻偷人精品xxx网站 | 精品久久久久久久人妻蜜臀av| 国产一级毛片七仙女欲春2| 免费av毛片视频| 亚洲国产欧美一区二区综合| 桃色一区二区三区在线观看| 中文在线观看免费www的网站| 成人精品一区二区免费| 国产亚洲精品久久久久久毛片| 国产高清激情床上av| 91在线精品国自产拍蜜月 | 午夜日韩欧美国产| 97人妻精品一区二区三区麻豆| 国产精品自产拍在线观看55亚洲| 国产成人一区二区三区免费视频网站| 久久久久免费精品人妻一区二区| 美女午夜性视频免费| 日日干狠狠操夜夜爽| www国产在线视频色| 亚洲专区字幕在线| 久久亚洲精品不卡| 男插女下体视频免费在线播放| 九色国产91popny在线| 亚洲熟女毛片儿| 国产精品,欧美在线| 桃红色精品国产亚洲av| 午夜精品一区二区三区免费看| 全区人妻精品视频| 亚洲av免费在线观看| 免费观看的影片在线观看| 黄色丝袜av网址大全| 一二三四在线观看免费中文在| 九九久久精品国产亚洲av麻豆 | 无限看片的www在线观看| 久久久国产欧美日韩av| 老鸭窝网址在线观看| 亚洲男人的天堂狠狠| 激情在线观看视频在线高清| 久久精品国产清高在天天线| 欧美日韩福利视频一区二区| 真人做人爱边吃奶动态| 午夜福利欧美成人| 亚洲av熟女| 欧美日韩亚洲国产一区二区在线观看| 无限看片的www在线观看| 在线a可以看的网站| 黄色 视频免费看| 国产av在哪里看| 久久精品aⅴ一区二区三区四区| 美女高潮喷水抽搐中文字幕| 国产精品自产拍在线观看55亚洲| 十八禁人妻一区二区| 精品一区二区三区四区五区乱码| 亚洲国产中文字幕在线视频| 99久久99久久久精品蜜桃| 亚洲av五月六月丁香网| 国产高潮美女av| 亚洲欧美日韩无卡精品| 亚洲av美国av| 欧美成人性av电影在线观看| 麻豆国产av国片精品| 国产单亲对白刺激| 中文字幕av在线有码专区| 黄色日韩在线| 97人妻精品一区二区三区麻豆| 欧美黑人欧美精品刺激| 亚洲国产色片| 又黄又爽又免费观看的视频| 日本精品一区二区三区蜜桃| 成人鲁丝片一二三区免费| 老熟妇乱子伦视频在线观看| 88av欧美| 亚洲一区高清亚洲精品| 国产亚洲精品久久久久久毛片| 久久精品影院6| 很黄的视频免费| 不卡av一区二区三区| 久久亚洲精品不卡| 女警被强在线播放| 激情在线观看视频在线高清| 麻豆国产97在线/欧美| 日本撒尿小便嘘嘘汇集6| 国产免费男女视频| 在线观看一区二区三区| 别揉我奶头~嗯~啊~动态视频| 欧美日本视频| 97超视频在线观看视频| 欧美日韩一级在线毛片| 色老头精品视频在线观看| 国产欧美日韩一区二区三| 亚洲成av人片免费观看| 老司机深夜福利视频在线观看| 精品国产三级普通话版| 久9热在线精品视频| 亚洲av中文字字幕乱码综合| 男人舔女人下体高潮全视频| 可以在线观看的亚洲视频| 国产伦一二天堂av在线观看| 97超视频在线观看视频| 免费看日本二区| 亚洲色图av天堂| 亚洲精品色激情综合| 黄色片一级片一级黄色片| 搞女人的毛片| 青草久久国产| 欧美高清成人免费视频www| 亚洲国产精品999在线| a级毛片在线看网站| 亚洲av成人精品一区久久| 日本黄大片高清| 久久精品综合一区二区三区| 国产精品永久免费网站| 久久久久亚洲av毛片大全| 搡老妇女老女人老熟妇| 日日夜夜操网爽| 亚洲自拍偷在线| 成人无遮挡网站| 黑人欧美特级aaaaaa片| 最好的美女福利视频网| 香蕉国产在线看| 亚洲国产日韩欧美精品在线观看 | 观看美女的网站| 亚洲第一电影网av| 日本五十路高清| 天天躁日日操中文字幕| 看黄色毛片网站| 少妇人妻一区二区三区视频| 国产欧美日韩一区二区精品| 国产精品一区二区精品视频观看| 国产精品自产拍在线观看55亚洲| 国产精品一区二区免费欧美| 亚洲自拍偷在线| 黄片大片在线免费观看| 日韩国内少妇激情av| 国产精品一及| 精品熟女少妇八av免费久了| 国产亚洲av高清不卡| 亚洲av电影在线进入| 国产精品永久免费网站| 欧美av亚洲av综合av国产av| 给我免费播放毛片高清在线观看| 国产精品乱码一区二三区的特点| 亚洲中文av在线| 91麻豆av在线| 午夜福利在线在线| 后天国语完整版免费观看| 搡老熟女国产l中国老女人| 久久精品aⅴ一区二区三区四区| 日韩有码中文字幕| 黑人巨大精品欧美一区二区mp4| 午夜福利免费观看在线| 久久人妻av系列| 亚洲真实伦在线观看| 中国美女看黄片| 中文字幕高清在线视频| 88av欧美| 欧美日韩一级在线毛片| 99久久精品一区二区三区| 最新在线观看一区二区三区| 欧美成人一区二区免费高清观看 | 欧美激情久久久久久爽电影| 精品国产三级普通话版| 日本免费a在线| 午夜激情欧美在线| 午夜久久久久精精品| 女同久久另类99精品国产91| 久久婷婷人人爽人人干人人爱| 99国产精品一区二区蜜桃av| 亚洲第一欧美日韩一区二区三区| 一个人看视频在线观看www免费 | 巨乳人妻的诱惑在线观看| 免费看美女性在线毛片视频| 两人在一起打扑克的视频| 99久久99久久久精品蜜桃| 2021天堂中文幕一二区在线观| 天天添夜夜摸| 在线看三级毛片| 韩国av一区二区三区四区| 男女午夜视频在线观看| 一夜夜www| 日韩欧美在线二视频| 国产精品一区二区精品视频观看| 少妇裸体淫交视频免费看高清| 国产成人av激情在线播放| 美女午夜性视频免费| 国产精品美女特级片免费视频播放器 | 亚洲国产中文字幕在线视频| 在线视频色国产色| 好男人在线观看高清免费视频| 真实男女啪啪啪动态图| 国产成人系列免费观看| 午夜激情欧美在线| 国产午夜福利久久久久久| 69av精品久久久久久| 亚洲国产中文字幕在线视频| 国产一区二区三区视频了| 少妇裸体淫交视频免费看高清| 国产单亲对白刺激| 国产极品精品免费视频能看的| 51午夜福利影视在线观看| 国产激情欧美一区二区| 国产高清视频在线观看网站| www.熟女人妻精品国产| 巨乳人妻的诱惑在线观看| 狂野欧美白嫩少妇大欣赏| 精品不卡国产一区二区三区| 99久久无色码亚洲精品果冻| 亚洲国产欧美网| 五月玫瑰六月丁香| 久久久久精品国产欧美久久久| 久久国产乱子伦精品免费另类| 99在线视频只有这里精品首页| 老司机在亚洲福利影院| 观看美女的网站| 国产伦精品一区二区三区四那| 草草在线视频免费看| 国产亚洲欧美98| 久久久久久人人人人人| 国产亚洲精品久久久com| 亚洲最大成人中文| 国产精品av久久久久免费| 亚洲在线自拍视频| 精品乱码久久久久久99久播| 日本三级黄在线观看| 午夜福利18| av女优亚洲男人天堂 | 亚洲欧美日韩高清专用| 国产精品亚洲一级av第二区| 嫩草影院入口| 欧美日韩瑟瑟在线播放| 舔av片在线| 国内精品久久久久精免费| 国产主播在线观看一区二区| 欧美另类亚洲清纯唯美| 亚洲国产精品999在线| 黄色片一级片一级黄色片| 精品国产亚洲在线| 午夜福利在线观看吧| 亚洲va日本ⅴa欧美va伊人久久| 久久久久久国产a免费观看| 国产精品一区二区免费欧美| 久久中文看片网| 国产激情久久老熟女| av天堂在线播放| 99精品久久久久人妻精品| 黑人巨大精品欧美一区二区mp4| 日日干狠狠操夜夜爽| 中文资源天堂在线| 露出奶头的视频| 精品一区二区三区av网在线观看| 色吧在线观看| 欧美日韩亚洲国产一区二区在线观看| www日本黄色视频网| 久久这里只有精品19| 午夜福利在线观看免费完整高清在 | 国产成人福利小说| 两个人看的免费小视频| 97超视频在线观看视频| 两个人看的免费小视频| 精品国产美女av久久久久小说| 国内精品久久久久久久电影| 久久久久亚洲av毛片大全| 久99久视频精品免费| 国产乱人伦免费视频| 日本a在线网址| 亚洲成av人片在线播放无| 国产精品一区二区三区四区久久| 精品久久久久久久久久免费视频| 色吧在线观看| 黑人欧美特级aaaaaa片| 色综合亚洲欧美另类图片| 中文字幕精品亚洲无线码一区| 午夜福利在线观看免费完整高清在 | 久久人人精品亚洲av| 操出白浆在线播放| 三级国产精品欧美在线观看 | 又粗又爽又猛毛片免费看| 久久久色成人| 99久久久亚洲精品蜜臀av| 午夜免费激情av| 国产精品98久久久久久宅男小说| 最近最新免费中文字幕在线| 婷婷精品国产亚洲av| 51午夜福利影视在线观看| 成人永久免费在线观看视频| 听说在线观看完整版免费高清| 亚洲国产中文字幕在线视频| 九九久久精品国产亚洲av麻豆 | 2021天堂中文幕一二区在线观| 在线免费观看不下载黄p国产 | 久99久视频精品免费| 欧美不卡视频在线免费观看| 欧美在线黄色| 成人特级黄色片久久久久久久| 丰满人妻一区二区三区视频av | 叶爱在线成人免费视频播放| 制服人妻中文乱码| 国产一区在线观看成人免费| 女警被强在线播放| 久久香蕉精品热| 欧美日韩乱码在线| 啦啦啦韩国在线观看视频| 日韩欧美一区二区三区在线观看| 色视频www国产| 国内精品美女久久久久久| 99国产综合亚洲精品| 99精品在免费线老司机午夜| 久久精品夜夜夜夜夜久久蜜豆| 精品一区二区三区四区五区乱码| 熟女电影av网| 舔av片在线| 久久热在线av| 国产精品香港三级国产av潘金莲| 国产精品美女特级片免费视频播放器 | 精品电影一区二区在线| а√天堂www在线а√下载| 啦啦啦观看免费观看视频高清| 男人舔女人下体高潮全视频| 两人在一起打扑克的视频| 两个人的视频大全免费| 久久亚洲真实| 国产熟女xx| 美女黄网站色视频| 久久久久国产一级毛片高清牌| 日韩人妻高清精品专区| 又大又爽又粗| 欧美日本视频| 欧美黑人欧美精品刺激| 欧美另类亚洲清纯唯美| 大型黄色视频在线免费观看| 久久久久国内视频| 亚洲成av人片在线播放无| 国产午夜精品久久久久久| 高清在线国产一区| 欧美日本亚洲视频在线播放| 国产淫片久久久久久久久 | 亚洲av电影在线进入| 亚洲欧美日韩高清在线视频| 日本 av在线| 精品日产1卡2卡| 2021天堂中文幕一二区在线观| 久久这里只有精品19| 18禁观看日本| 麻豆成人午夜福利视频| 欧美日韩国产亚洲二区| 精品国产美女av久久久久小说| 欧美中文综合在线视频| 好男人电影高清在线观看| 欧美激情在线99| 一本综合久久免费| 成在线人永久免费视频| 国产探花在线观看一区二区| 久久中文看片网| 国产 一区 欧美 日韩| 男女做爰动态图高潮gif福利片| 老司机午夜福利在线观看视频| 国产野战对白在线观看| 美女 人体艺术 gogo| 一区二区三区高清视频在线| 免费观看的影片在线观看| 变态另类成人亚洲欧美熟女| 日韩av在线大香蕉| 看黄色毛片网站| 国产一区二区在线av高清观看| 国产aⅴ精品一区二区三区波| 男女做爰动态图高潮gif福利片| 亚洲欧美日韩卡通动漫| e午夜精品久久久久久久| 国内精品美女久久久久久| 在线观看午夜福利视频| 国产毛片a区久久久久| 久久久国产精品麻豆| 一级毛片女人18水好多| 1024香蕉在线观看| 日韩欧美在线二视频| 2021天堂中文幕一二区在线观| 久久性视频一级片| 色老头精品视频在线观看| or卡值多少钱| 日韩欧美在线乱码| 欧美成人性av电影在线观看| 国产精品亚洲一级av第二区| 精品熟女少妇八av免费久了| 久久国产乱子伦精品免费另类| 精品99又大又爽又粗少妇毛片 | 色尼玛亚洲综合影院| 久久中文看片网| 成人欧美大片| 国产黄色小视频在线观看| 亚洲午夜精品一区,二区,三区| 国产精品久久视频播放| 精品一区二区三区视频在线观看免费| 无人区码免费观看不卡| 成人三级黄色视频| 久久伊人香网站| 99精品欧美一区二区三区四区| 日日夜夜操网爽| 国产精品av久久久久免费| 亚洲午夜精品一区,二区,三区| 男人的好看免费观看在线视频| 每晚都被弄得嗷嗷叫到高潮| 性欧美人与动物交配| 国产不卡一卡二| 久久热在线av| 亚洲狠狠婷婷综合久久图片| 三级男女做爰猛烈吃奶摸视频| 免费观看的影片在线观看| 大型黄色视频在线免费观看| 欧美一级毛片孕妇| 久久精品国产清高在天天线| 可以在线观看毛片的网站| 国产精品自产拍在线观看55亚洲| h日本视频在线播放| 亚洲色图 男人天堂 中文字幕| 999精品在线视频| 99久久无色码亚洲精品果冻| 悠悠久久av| 国产69精品久久久久777片 | 国内精品久久久久精免费| 久久久精品欧美日韩精品| 婷婷亚洲欧美| 最近最新中文字幕大全免费视频| 特级一级黄色大片| 国产真人三级小视频在线观看| 国产精品一及| 久久久国产成人精品二区| 国产高潮美女av| 国模一区二区三区四区视频 | 久久亚洲精品不卡| 亚洲av电影不卡..在线观看| 亚洲国产精品sss在线观看| 亚洲欧美激情综合另类| 国产乱人视频| 给我免费播放毛片高清在线观看| 黑人巨大精品欧美一区二区mp4| 色播亚洲综合网| 色老头精品视频在线观看| 18禁美女被吸乳视频| 一个人看视频在线观看www免费 | 两个人视频免费观看高清| 无限看片的www在线观看| 日韩欧美免费精品| 亚洲国产欧美网| 听说在线观看完整版免费高清| 国产亚洲精品综合一区在线观看| 九九热线精品视视频播放| 啦啦啦韩国在线观看视频| 老汉色av国产亚洲站长工具| 啦啦啦观看免费观看视频高清| 十八禁网站免费在线| 天堂网av新在线| 哪里可以看免费的av片| 久久久精品大字幕| 香蕉av资源在线| 亚洲精品乱码久久久v下载方式 | 亚洲av电影在线进入| 一个人看的www免费观看视频| 后天国语完整版免费观看| 亚洲欧美日韩无卡精品| 久久人妻av系列| 极品教师在线免费播放| 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产色片| 一a级毛片在线观看| 亚洲在线观看片| 黄频高清免费视频| 日韩中文字幕欧美一区二区| 免费在线观看成人毛片| 99久久精品热视频| 亚洲国产欧美网| 一边摸一边抽搐一进一小说| 人人妻,人人澡人人爽秒播| 午夜成年电影在线免费观看| 国产午夜精品久久久久久| 成人国产综合亚洲| 麻豆成人午夜福利视频| 男人舔女人下体高潮全视频| 国产伦精品一区二区三区四那| 欧美zozozo另类| 久久久久久九九精品二区国产| 精品99又大又爽又粗少妇毛片 | 老司机深夜福利视频在线观看| 午夜福利视频1000在线观看| 亚洲欧美日韩无卡精品| 精品国产乱子伦一区二区三区| 精品国产乱码久久久久久男人| 全区人妻精品视频| 久久精品国产99精品国产亚洲性色| 深夜精品福利| 亚洲自偷自拍图片 自拍| 后天国语完整版免费观看| 日韩成人在线观看一区二区三区| 亚洲天堂国产精品一区在线| 免费看美女性在线毛片视频| 不卡一级毛片| 欧美日韩亚洲国产一区二区在线观看| 成年女人永久免费观看视频| av福利片在线观看| 国产视频内射| 丰满人妻一区二区三区视频av | 日本成人三级电影网站| xxxwww97欧美| 国产高清视频在线观看网站| 色老头精品视频在线观看| 亚洲人与动物交配视频| 久久精品aⅴ一区二区三区四区| 午夜日韩欧美国产| 黄片大片在线免费观看| 丁香欧美五月| 国产激情偷乱视频一区二区| 18美女黄网站色大片免费观看| 欧美日本亚洲视频在线播放| 久久久精品大字幕| 亚洲国产精品成人综合色| 此物有八面人人有两片| 国产aⅴ精品一区二区三区波| 黄频高清免费视频| 老熟妇仑乱视频hdxx| 国产精品永久免费网站| 亚洲九九香蕉| 国产极品精品免费视频能看的| 成人三级做爰电影| 国产午夜精品久久久久久| 亚洲专区国产一区二区| 国产97色在线日韩免费| 五月玫瑰六月丁香| 精品福利观看| 夜夜夜夜夜久久久久| 丰满人妻一区二区三区视频av | 夜夜看夜夜爽夜夜摸| 性欧美人与动物交配| 亚洲精品国产精品久久久不卡| 小说图片视频综合网站| 成年女人毛片免费观看观看9| 成年女人永久免费观看视频| 久久香蕉精品热| 老司机福利观看|