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

    基于密切值法和水氮管理模型的華北平原農田水氮優(yōu)化管理

    2017-11-24 06:08:18李周晶胡克林李保國
    農業(yè)工程學報 2017年14期
    關鍵詞:水氮夏玉米施肥量

    徐 強,李周晶,胡克林,李保國

    (中國農業(yè)大學資源與環(huán)境學院,農業(yè)部華北耕地保育重點實驗室,北京100193)

    基于密切值法和水氮管理模型的華北平原農田水氮優(yōu)化管理

    徐 強,李周晶,胡克林※,李保國

    (中國農業(yè)大學資源與環(huán)境學院,農業(yè)部華北耕地保育重點實驗室,北京100193)

    優(yōu)化農田水氮管理措施可為實現糧食高產、資源高效及環(huán)境友好的目標提供科學依據。該研究以華北平原泰安地區(qū)為例,利用農田生態(tài)系統(tǒng)水熱碳氮過程耦合模型(soil water heat carbon and nitrogen simulator,WHCNS)分別對冬小麥季設置的165個水肥組合和夏玉米季設置的55個水肥組合進行了情景模擬分析,在綜合考慮農學、環(huán)境和經濟效益的基礎上,采用密切值法優(yōu)化了農田水肥管理方案。結果表明:受到華北地區(qū)年內降雨分配不均的影響,冬小麥產量隨著灌水量的增加呈先增加后穩(wěn)定的趨勢;而 夏玉米產量與灌水量沒有明顯的關系。冬小麥和夏玉米產量均隨著施肥量的增加而增加,后保持穩(wěn)定。水分滲漏和氮素淋洗量均隨著灌水量或施肥量的增加而顯著增加。在研究區(qū)作物秸稈全部還田及高累積氮的條件下,冬小麥季灌水240 mm和施肥60 kg/hm2(以N計,下同),夏玉米季不灌溉和施肥90 kg/hm2分別為研究區(qū)當年冬小麥季和夏玉米季最佳的水肥管理方案。在所有水肥組合情景中,優(yōu)化的水肥管理方案不僅能保證冬小麥-夏玉米最大周年產量的97%、具有較高的水氮利用效率和最佳的產投比,而且氮素淋洗和氣體損失分別比最大值降低了77%和71%。因此,該方法可以用來優(yōu)化華北平原農田的水肥管理措施。

    灌溉;氮;優(yōu)化;冬小麥-夏玉米輪作;WHCNS模型;情景分析;密切值法

    0 引 言

    華北平原是中國冬小麥和夏玉米主要種植區(qū),同時也是中國水肥高投入區(qū)。多年來長期過度的開采地下水進行灌溉,已導致該地區(qū)成為世界上最大的“漏斗區(qū)”,造成了嚴重的生態(tài)問題[1]。另一方面,氮肥的過量投入使得該地區(qū)生態(tài)環(huán)境問題進一步加劇,如土壤酸化、地表水富營養(yǎng)化和地下水硝酸鹽污染等[2-3]。不僅如此,水氮的大量投入還造成了資源的巨大浪費。在華北平原冬小麥-夏玉米輪作區(qū),冬小麥季廣泛采用大水漫灌,全生育期灌水多達4~6次,灌水總量達500 mm,水分利用效率卻較低[4]。而農民為了追求高產普遍過量施用氮肥,華北地區(qū)氮肥平均用量為550 kg/(hm2·a),而作物產量增加有限[5-6]。由于中國耕地日益減少、人口不斷增加,未來中國糧食的需求還會繼續(xù)增長,只有持續(xù)增產才能滿足國家糧食安全的需求[7]。因此,減少該地區(qū)的水肥投入,制定科學合理的水肥管理措施對于該地區(qū)的農業(yè)可持續(xù)發(fā)展十分必要[8]。

    利用不同水肥組合的田間試驗來提高作物產量及水氮利用效率的研究已有大量報道,然而由于田間試驗費時費力,缺乏系統(tǒng)性并受制于氣候的多變性,如何優(yōu)化水肥用量始終是田間試驗面臨的難題。一些研究指出,土壤-作物系統(tǒng)模型能夠有效地模擬水氮運移及作物生長過程,從而在保證環(huán)境友好的條件下給出符合實際情況的最佳農田管理方案[9-10]。Hu等[11]用RZWQM模型評估了一年兩作體系上的不同氮肥管理效果,發(fā)現灌溉和施肥量均減半的情況下可獲得很好的經濟和環(huán)境效益。陳研等[12]在土壤-作物系統(tǒng)過程模型的基礎上,用動態(tài)規(guī)劃的方法對河北曲周田間水、氮資源管理措施進行了優(yōu)化,發(fā)現優(yōu)化后的水、氮利用效率分別比對照高出 13%和16%。Zhang等[9]利用DNDC模型優(yōu)化了華北平原夏玉米季的施氮量,結果表明當施氮180 kg/hm2(以N計,下同)時,產量可達到最大值,同時硝態(tài)氮的淋洗低于18.4 kg/hm2。其他一些學者的研究表明,優(yōu)化施肥可以減少土壤剖面殘留硝態(tài)氮的積累,極大地降低氮素淋失的環(huán)境風險[13-14]。然而以上研究僅僅考慮了農學、環(huán)境或經濟效益的某一方面,具有一定的片面性。目前,對于農田水肥管理優(yōu)化方案缺乏一個統(tǒng)一的優(yōu)化評價方法。

    密切值方法是一種多目標決策的選優(yōu)方法,被廣泛地應用于評價農業(yè)經濟項目、地下水污染和水質的研究中[15-16],目前還未在農田水肥優(yōu)化管理中進行應用。最近,本研究小組針對中國集約化的農業(yè)生產模式開發(fā)了農田水氮管理模型(soil water heat carbon and nitrogen simulator,WHCNS),Liang等[10]應用該模型對中國西北綠洲地區(qū)不同水肥管理模式的灌溉制度進行了優(yōu)化;Li等[17]應用WHCNS模型對華北平原小麥-玉米典型輪作區(qū)的不同水肥與栽培管理組合模式進行了校驗和評估。然而上述研究在優(yōu)化目標選取上比較單一,僅僅考慮了硝酸鹽淋失和作物產量,缺乏一個能綜合考慮農學、環(huán)境和經濟效益的擇優(yōu)方法。因此,本研究以華北平原泰安地區(qū)為例,將 WHCNS模型與密切值方法相結合對農田水肥組合模式進行綜合評價,在綜合考慮農學、環(huán)境和經濟效益的基礎上,得到最佳水肥管理方案,旨在為華北平原同時實現糧食高產、資源高效及環(huán)境友好的目標提供技術支持。

    1 材料與方法

    1.1 研究區(qū)概況

    本研究在山東省泰安市大汶口(35°58′N, 117°03′E)進行為期2 a(2009年10月—2011年9月)的田間試驗。該試驗點屬溫帶大陸性季風氣候,年均日照時數為2 627 h。年均氣溫12 ℃,年均降水量為697 mm,多集中在7—9月。試驗田主要為冬小麥-夏玉米輪作,土壤類型為沖積始成土。表土有機質為16.7 g/kg,全氮為1.01 g/kg,速效磷為52.65 mg/kg,速效鉀為96.15 mg/kg。該試驗田土壤基本理化性質見表1。

    1.2 試驗設計

    試驗設置4種不同水肥及栽培管理模式:1)農民傳統(tǒng)模式(FP);2)高產高效模式(產量和氮肥利用效率比FP提高15%~20%,OPT-1);3)再高產模式(HY),即不計資源投入成本使產量最大化;4)再高產高效模式(產量和氮肥利用效率比FP提高30%~50%,OPT-2)。小區(qū)設計采用4個重復的完全隨機區(qū)組設計,每個小區(qū)6 m × 40 m。試驗品種冬小麥為泰農18,夏玉米為鄭單958。FP模式代表了該地區(qū)的傳統(tǒng)水氮管理和栽培模式。OPT-1和OPT-2模式均建立在FP模式基礎上,其灌溉日期和數量以及施肥日期和數量等均根據作物需求來調整。HY模式指不計資源投入成本,通過優(yōu)化田間管理措施使產量最大化(表 2)。詳細的田間水肥管理措施、測定過程和方法見文獻[17]。

    表1 研究區(qū)土壤剖面基本理化性質及水力學特性參數Table 1 Basic physicochemical and hydraulic properties for soil profile in study area

    表2 2009—2011年不同處理的具體田間管理措施Table 2 Field management recorded for different treatments from 2009 to 2011

    1.3 WHCNS模型及情景設置

    1.3.1 WHCNS模型簡介

    農田生態(tài)系統(tǒng)水熱碳氮過程耦合模型(Soil Water Heat Carbon and Nitrogen Simulator,WHCNS),主要包括氣象、土壤水運動、土壤熱傳導、氮素運移及轉化、有機質周轉、作物生長和田間管理等模塊。該模型以天為步長,由氣象數據和作物生物學參數驅動。在模型中,采用Penman-Monteith公式估算參考作物蒸散量[18]。土壤水分入滲和再分布過程分別采用 Green-Ampt模型和Richard’s方程進行模擬。根系吸水采用 van Genuchten模型,并引入了補償性吸水機制。土壤熱運動采用對流-傳導方程來描述。土壤無機氮的運移采用對流-擴散方程,源匯項中考慮了碳氮循環(huán)各過程(有機質礦化、生物固持、尿素水解、氨揮發(fā)、硝化和反硝化等)和作物吸收,土壤有機質周轉動態(tài)直接來源于Daisy模型[19]。作物生長發(fā)育進程、干物質生產和分配及作物產量的模擬采用荷蘭的PS123作物模型[20],通過水氮脅迫校準因子來實現水氮限制下作物產量的模擬。該模型詳細原理見文獻[21]。

    本研究中WHCNS模型參數的確定采用“試錯法”。利用HY模式田間試驗實測數據對模型參數進行校準,然后用其他處理(FP、OPT-1和OPT-2)的實測數據驗證模型。模型模擬的土壤含水率、土壤硝態(tài)氮、作物產量和葉面積指數(leaf area index,LAI)與實測值之間的線性相關系數分別為0.90、0.84、0.88和0.91,均高于0.80[17],說明校準后的 WHCNS模型可用于該地區(qū)不同水氮管理模式下土壤含水率、土壤硝態(tài)氮、作物產量和LAI的模擬。具體的模型輸入參數和模型校驗過程見文獻[17]。

    1.3.2 水氮情景設置

    為了優(yōu)化水氮管理措施,以2009—2010年OPT-2模式為例進行情景分析。根據當地農民習慣水肥投入量,在冬小麥生長季設置了15個灌水量和11個施氮量處理,共得到 165個水肥組合情景。在夏玉米生長季設置了 5個灌水量和11個施肥量處理,共得到55個水肥組合情景。水氮的投入情景如下:

    冬小麥季:1)灌水量范圍為0~420 mm,依次遞增30 mm。灌水次數和時間同OPT-2模式,每次灌水量相同;2)施氮量范圍為0~300 kg/hm2,依次遞增30 kg/hm2。

    施肥時間不變,基肥與追肥比為4∶6。

    夏玉米季:1)灌水量范圍為0~120 mm,依次遞增30 mm,灌水時間為播種后;2)施氮量范圍為 0~300 kg/hm2,依次遞增30 kg/hm2。施肥時間不變,基肥與追肥比為1∶5。

    1.4 密切值法

    密切值法的基本原理就是找出一個方案或者一個指標的最優(yōu)點和最劣點,而最接近最優(yōu)點,離最劣點最遠的方案為最佳的方案[15-16]。對于某一水肥管理方案,假定有m個水肥管理措施(Q1,Q2,…Qm)和n個綜合評價指標(A1,A2,…An)??梢韵冉?1個有量綱的矩陣,然后將每個元素進行歸一化,可以得到新的無量綱矩陣C:

    式中Cij為第i個水氮管理措施(Qi)下的第j個綜合評價指標(Aj)歸一化后的值。令:

    在水肥管理評價指標中由于有正向指標和逆向指標,因此在式(2)中“+”表示正向指標;“-”表示逆向指標。從而建立無量綱樣本矩陣:

    因此,一個樣本集的“最優(yōu)點”應為評價指標最小值的集合QG,一個樣本集的“最劣點”為評價指標最大值的集合QB,分別為

    其次,最優(yōu)點QG和最劣點QB為所有樣本中評價指標的虛擬點的集合,求出各個樣本點與這些虛擬點(最優(yōu)或者最劣)的距離就可以為水肥管理的評價提供一個定量的依據。在應用時可以增加各項指標的權重系數,使結果更合理。第i個評價方案樣本與“最優(yōu)點”和“最劣點”集合的歐氏距離di-G和di-B分別為

    其中ωj為第j個評價指標的權重,

    最后,計算密切值Ej。該值反映某個樣本點與某一端點的接近程度,它的計算過程實質上是將樣本點Qi與最優(yōu)點集合和最劣點集合的距離進一步無量綱化處理。根據端點的不同,密切值分為最優(yōu)密切值Ei-G和最劣密切值Ei-B。

    當評價方案的最優(yōu)密切值越小,最劣密切值越大,表明該水肥管理方案與最優(yōu)點越近,離最劣點越遠,該水肥管理方案越好。把水肥管理方案按照最優(yōu)密切值或者最劣密切值進行排序,就可以得到不同水肥管理方案的優(yōu)劣排序結果。

    2 結果與分析

    2.1 不同水氮管理組合情景下的作物產量及水氮利用效率

    利用校驗后的 WHCNS模型分別對冬小麥設置的165個水肥組合情景和夏玉米設置的55個水肥組合情景進行模擬分析,得到了不同水肥組合情景下冬小麥和夏玉米的作物產量、水分滲漏、氮素淋洗和水氮利用效率(圖1)。由圖1a可知,冬小麥的產量與灌水量密切相關,冬小麥的產量隨著灌水量的增加而增加,當灌水量達到251 mm時,冬小麥的產量達到最大值。再繼續(xù)增加灌水量,作物產量基本保持不變。而夏玉米產量和灌水量之間并沒有明顯的關系,這主要是由于華北平原約 70%的降雨主要集中在7-9月的夏玉米生育期,夏玉米僅靠降雨就可滿足其水分需求,因此灌水量與其產量之間沒有明顯的規(guī)律。由圖 1b可知,當冬小麥季灌水量小于150 mm時,水分滲漏量幾乎為 0;但當灌水量大于150 mm 時,水分滲漏量迅速增加。而夏玉米季的水分滲漏量隨著灌水量增加而增加,由于該年份夏玉米的降雨量較多(447 mm),即使不灌水,水分滲漏也可達到31 mm。另外,土體N的淋洗量與水分滲漏量成正比(圖1c)。在冬小麥季當灌水量小于150 mm時,N素淋洗量幾乎為0;當超過150 mm時,N素淋洗量迅速增加。在夏玉米季N素淋洗與水分滲漏的動態(tài)變化趨勢基本一致,隨著灌水量的增加而增加。圖1d顯示了不同水氮投入情景下模型模擬計算得到的冬小麥和夏玉米的水分利用效率(water use efficiency,WUE)。在冬小麥季,WUE隨著灌水量增加呈增加的趨勢,當灌水量達到270 mm時,WUE達到最大值,繼續(xù)增加灌水量,WUE幾乎保持不變;夏玉米季,WUE和灌水量之間的關系并不明顯。圖 1e是氮素利用效率(nitrogen use efficiency,NUE)與灌水量和施氮量的關系。在冬小麥季,NUE隨著灌水量的增加而增加,當灌水量增加到90 mm時,NUE達到最大值。隨著灌水量的繼續(xù)增加,NUE基本保持不變;在夏玉米季,NUE和灌水量之間的關系不明顯。

    圖1 不同施氮水平下作物產量、水氮損失、水分利用效率(WUE)和氮素利用效率(NUE)隨灌水量變化Fig. 1 Change in crop yield, drainage, water-N loss, water use efficiency (WUE) and nitrogen use efficiency (NUE) with irrigation amount under different fertilizer levels

    作物產量與氮肥投入量的關系如圖 2所示,冬小麥產量隨著氮肥投入的增加先增加后保持不變。夏玉米產量隨著氮肥投入的增加而增加,當氮肥投入量達到70 kg/hm2時,夏玉米產量達到最大值,為9 005 kg/hm2;繼續(xù)增加施肥量,作物產量基本保持不變。

    圖2 不同灌水水平下作物產量與施氮量的關系Fig. 2 Relationship between fertilizer-N application rate and crop yield under different irrigation levels

    2.2 基于密切值法的農田水氮管理措施優(yōu)化

    為了獲得最佳水氮管理方案,需要將不同水氮組合的輸出結果進行分類并附以權重系數。在計算密切值時,本研究考慮了農學效應(作物產量、WUE和NUE)、環(huán)境效應(N淋洗和氣體損失)和經濟效益(產投比),其權重系數分別設為0.7、-0.2和0.5[17]。由于不施肥或不灌溉時的經濟效益和環(huán)境效應最好,因此在選擇最佳水肥管理方案的過程中添加了 1個限制條件:要求冬小麥和夏玉米產量分別高于7 500和8 500 kg/hm2。按照密切值排序,最終評價結果見表3。

    對于冬小麥來說,灌水量420 mm結合施肥300 kg/hm2的管理方案,具有最低的產投比8.2,其最優(yōu)密切值最大且最劣密切值最小,分別為10.96和0.04。該方案下N淋洗和N的氣體損失量分別高達65.7和57.9 kg/hm2,為所有情景中最高,而NUE最低,為22.9 kg/kg,說明該水肥組合管理方案是最差的。當灌水量降低到240 mm,同時施肥量減少到60 kg/hm2時,最優(yōu)密切值最小且最劣密切值最大,分別為9.51和0.90。此時冬小麥產量僅比最大值降低了約5%,而N淋洗和N的氣體損失量分別減少到14.9和17.4 kg/hm2,NUE顯著提高至32.9 kg/kg,并且此方案具有最高的產投比 21.6,說明該情景能夠在保持較高產量的前提下有效降低環(huán)境影響,同時還節(jié)約了水氮資源,可見該水肥管理方案最佳。

    對于夏玉米來說,灌水量 120 mm結合施肥量300 kg/hm2的管理方案,具有最大的最優(yōu)密切值和最小的最劣密切值,此時產投比最低,為 10.9。該方案下氮素淋洗和氣體損失量分別達到77.6和56.1 kg/hm2,為所有情景中最高,而NUE為最低(24.1 kg/kg),說明該方案是最差的。在不灌水和施肥量為90 kg/hm2的管理方案下,夏玉米季的產量達到 8 909 kg/hm2,僅比最大產量下降了1.3%。氮素的淋洗和氣體損失量分別為17.5和15.2 kg/hm2,同時WUE和NUE為所有情景中最高,分別為2.66 kg/m3和39.1 kg/kg,而且此情景具有最大的產投比44.0,說明該方案能夠在保證較高產量的前提下維持較高的經濟效益,同時還能降低環(huán)境影響,可見該方案最佳。綜上,冬小麥季灌水量240 mm配合施氮肥60 kg/hm2,夏玉米季不灌水和施氮量為90 kg/hm2分別是本研究區(qū)該年份冬小麥和夏玉米季最佳的水肥管理方案,能保證冬小麥-夏玉米最大周年產量的 97%和最佳的產投比,而氮素淋洗和氣體損失分別比最高值降低了77%和71%。

    表3 2009-2010年冬小麥和夏玉米產量分別高于7 500和8 500 kg·hm-2的田間水肥管理Table 3 Management of irrigation (Ir) and fertilizer-N application rate (FN) in winter wheat with yield of 7 500 kg·hm-2 and summer maize with yield of 8 500 kg·hm-2 during period of 2009 to 2010

    3 討 論

    大量研究表明,合理地減少灌溉量或施肥量,不僅能提高糧食產量和水氮利用效率,而且能顯著減少水氮的淋失[4,9-10,22-23]。Sexton等[22]研究了美國文代爾地區(qū)水氮投入對玉米產量和硝酸鹽淋洗的影響,認為最佳施氮量可保證最大產量的95%,而硝酸鹽淋洗量可降低30%~40%。Fang等[23]在華北平原的研究表明:當灌水量由300 kg/hm2減少至200 kg/hm2時,不僅能保證作物不減產,還能使硝酸鹽淋洗量減少60%。Liang等[10]在內蒙古荒漠綠洲區(qū)的研究表明,優(yōu)化后的水氮管理措施可使水分滲漏和氮素淋洗分別降低 65%和59%。本研究表明,冬小麥產量隨灌水量增加而升高,當產量達到最大值時,再繼續(xù)增加灌水量,作物產量保持不變。由于夏季充沛的降雨,夏玉米產量與灌水量的關系并不明顯。冬小麥和夏玉米產量均隨著施肥量的增加而增加,達到平臺后保持不變,而且兩季作物的水分滲漏和氮素淋洗量均隨著灌水量或施肥量的增加而顯著增加,這與前人的研究結果是一致的[24-27]。

    由于影響水氮優(yōu)化結果的田間管理因素眾多[17,28-30],如灌溉制度、施肥模式、種植制度、秸稈還田及耕作方式等,不同學者在不同的地區(qū)得到的結果不盡相同。Wang等[31]建議在北京地區(qū)冬小麥和夏玉米季的最佳施肥量分別為110和120 kg/hm2。Zhang等[9]應用DNDC模型研究得到夏玉米季最佳施肥量為180 kg/hm2時,產量可達到最大值,而硝態(tài)氮的淋洗<18.4 kg/hm2。Zhao等[4]利用APSIM 模型對華北平原吳橋試驗點冬小麥-夏玉米水氮管理優(yōu)化的結果表明:冬小麥灌水總量為225 mm分3次灌溉,配合施氮330 kg/hm2(小麥150 kg/hm2,玉米180 kg/hm2)可使周年產量達到最高,且對環(huán)境的影響最小。本研究發(fā)現冬小麥季灌水240 mm和施肥60 kg/hm2,夏玉米季不灌溉和施肥90 kg/hm2為研究區(qū)最佳的水肥管理方案。本研究結果與上述研究結果不同,主要是以上研究僅僅考慮了作物產量和硝酸鹽淋洗的某個方面,而本研究綜合考慮了作物產量、水氮利用效率、氮素損失及產投比等多方面。另外,Zhao等[32]在北京地區(qū)推薦的冬小麥季施肥量為 55~104 kg/hm2,夏玉米季的施肥量為 30~122 kg/hm2。本研究結果與該結果比較一致,但是施肥量要低于其他地區(qū)的研究結果[4,9,31],這主要與本試驗地作物秸稈全部還田及冬小麥和夏玉米播種前土體無機氮累積量較高(分別為150和187 kg/hm2)有關[17]。

    然而,本研究也存在一些不足,首先在使用密切值法進行選優(yōu)時,各指標權重系數的設定具有一定的經驗性,在今后研究中可嘗試利用層次分析等方法更為客觀地選取權重系數。其次,由于在選優(yōu)過程中并沒能考慮降雨年型對優(yōu)化結果的影響,導致得到的結果存在一定的局限性,在這方面也需要完善。

    4 結 論

    1)不同水肥組合的情景分析結果表明,冬小麥產量隨著灌水量增加呈先增加后保持不變的趨勢,而夏玉米產量與灌水量沒有明顯的關系。冬小麥和夏玉米產量均隨著施肥量的增加先增加后保持不變。水分滲漏和氮素淋洗量均隨著灌水量的增加而顯著增加。冬小麥的WUE和NUE隨灌水量的增加而增加,當灌水量達到一定程度后基本保持不變,而夏玉米的WUE和NUE與灌水量之間的關系不明顯。

    2)在試驗地作物秸稈全部還田及土體多年高累積氮的條件下,灌水240 mm和施肥量60 kg/hm2為研究區(qū)冬小麥季的最佳水肥管理方案,產量可達到 7 813 kg/hm2;對于夏玉米,不灌溉和施肥量為90 kg/hm2是研究區(qū)的最佳水肥管理方案,產量可達到8 909 kg/hm2。

    3)在所有情景中,優(yōu)化的水肥管理方案不僅能保證冬小麥-夏玉米最大周年產量的97%和最佳的產投比,而氮素淋洗和氣體損失分別比最高值降低了77%和71%。該方法不僅能保證較高的農學和經濟效益,而且大大減少了氮素的環(huán)境影響。因此,密切值方法結合農田水氮管理模型WHCNS可以用來優(yōu)化華北平原農田的水肥管理措施。

    [1] Hu C, Delgado, J A, Zhang X, et al. Assessment of groundwater use by wheat (TriticumaestivumL.) in the Luancheng Xian Region and potential implications for water conservation in the Northwestern North China Plain[J].Journal of Soil and Water Conservation, 2005, 60(2): 80-88.

    [2] Guo J H, Liu X J, Zhang Y, et al. Significant acidification in major Chinese croplands[J]. Science, 2010, 327: 1008-1010.

    [3] Cui Zhenling, Zhang Fusuo, MiaoYuxin, et al. Soil nitrate-N levels required for high yield maize production in the North China Plain[J]. Nutrient Cycling in Agro-ecosystems, 2008,82(2): 187-196.

    [4] Zhao Zhigang, Qin Xin, Wang Enli, et al. Modelling to increase the eco-efficiency of a wheat-maize double cropping system[J]. Agriculture, Ecosystems and Environment, 2015,210: 36-46.

    [5] 陳新平,張福鎖,崔振嶺,等. 小麥-玉米輪作體系養(yǎng)分資源綜合管理理論與實踐[M]. 北京:中國農業(yè)大學出版社,2006:27-41.

    [6] Zhang Weifeng, Cao Guoxin, Li Xiaolin, et al. Closing yield gaps in China by empowering smallholder farmers[J]. Nature,2016, 537(7622): 671

    [7] Zhu Zhaoliang, Chen Deli. Nitrogen fertilizer use in China Contributions to food production, impacts on the environment and best management strategies[J]. Nutrient Cycling in Agroecosystems, 2002, 63: 117-127.

    [8] 王敬國,林杉,李保國. 氮循環(huán)與中國農業(yè)氮管理[J]. 中國農業(yè)科學,2016,49(3):503-517.Wang Jingguo, Lin Shan, Li Baoguo. Nitrogen cycling and management strategies in Chinese agriculture[J]. Scientia Agricultura Sinica, 2016, 49(3): 503-517. (in Chinese with English abstract)

    [9] Zhang Yitao, Wang Hongyuan, Liu Shen, et al. Identifying critical nitrogen application rate for maize yield and nitrate leaching in a Haplic Luvisol soil using the DNDC model[J].Science of the Total Environment, 2015, 514: 388-398.

    [10] Liang Hao, Qi Zhiming, Hu Kelin, et al. Can nitrate contaminated groundwater be remediated by optimizing flood irrigation rate with high nitrate water in a desert oasis using the WHCNS model?[J]. Journal of Environmental Management, 2016, 181: 16-25.

    [11] Hu C S, Saseendran S, Green T, et al. Evaluating nitrogen and water management in a double-cropping system using RZWQM[J]. Vadose Zone Journal, 2006, 5: 493-505.

    [12] 陳研,胡克林,馮凌,等. 基于土壤-作物系統(tǒng)模擬模型的冬小麥田間水氮優(yōu)化管理[J]. 農業(yè)工程學報,2007,23(6):55-60.Chen Yan, Hu Kelin, Feng Ling, et al. Optimal management of water and nitrogen for winter wheat based on simulation model in soil-plant system in agricultural field[J].Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2007, 23(6): 55-60. (in Chinese with English abstract)

    [13] 王歡元,胡克林,李保國,等. 不同管理模式下農田水氮利用效率及其環(huán)境效應[J]. 中國農業(yè)科學,2011,13:2701-2710.Wang Huanyuan, Hu Kelin, Li Baoguo, et al. Analysis of water and nitrogen use efficiencies and their environmental impact under different water and nitrogen management practices[J]. Scientia Agricultura Sinica, 2011, 13: 2701-2710. (in Chinese with English abstract)

    [14] 王立剛,李虎,楊黎,等. 冬小麥/夏玉米輪作系統(tǒng)不同施氮量的長期環(huán)境效應及區(qū)域氮調控模擬[J]. 中國農業(yè)科學,2013,14:2932-2941.Wang Ligang, Li Hu, Yang Li, et al. Simulation of long-term and regional environmental effects of different N applications in the winter wheat/summer maize system[J]. Scientia Agricultura Sinica, 2013, 14: 2932-2941. (in Chinese with English abstract)

    [15] 王登瀛. 多目標決策方案優(yōu)選的密切值法[J]. 系統(tǒng)工程,1989,7(1):33-35.Wang Dengying. Close value method for optimization of multi-objective decision[J]. Systems Engineering, 1989, 7(1):33-35. (in Chinese with English abstract)

    [16] 樓文高. 用改進的密切值法綜合評價農業(yè)技術經濟方案[J]. 農業(yè)系統(tǒng)科學與綜合研究,2002,18(2):92-95.Lou Wengao. Comprehensive assessment of agricultural economic projects using osculating method[J]. SSCSA, 2002,18(2): 92-95. (in Chinese with English abstract)

    [17] Li Zhoujing, Hu Kelin, Li Baoguo, et al. Evaluation of water and nitrogen use efficiencies in a double cropping system under different integrated management practices based on a model approach[J]. Agricultural Water Management, 2015,159: 19-34.

    [18] Allen R G, Pereira L S, Raes D, et al. Crop evapotranspiration guidelines for computing crop water requirements. FAO Irrigation and drainage paper No. 56[R].Rome: FAO, 1998, 65-78.

    [19] Hansen S, Jcnsen H E, Niclsen N E, et al. NPo-research, A10:DAISY: Soil plant atmosphere system model[M]. Copenhagen:The National Agency for Environmental Protection, 1990.

    [20] Driessen P M, Konjin N T. Land use system analysis[M].Netherlands: Wageningen Agricultural University, 1992.

    [21] 梁浩,胡克林,李保國,等. 土壤-作物-大氣系統(tǒng)水熱碳氮過程耦合模型構建[J]. 農業(yè)工程學報,2014,30(24):54-66.Liang Hao, Hu Kelin, Li Baoguo, et al. Coupled simulation of soil water-heat-carbon-nitrogen process and crop growth at soil-plant-atmosphere continuum system and its validation[J].Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(24): 54-66. (in Chinese with English abstract)

    [22] Sexton B T, Moncrief J F, Rosen C J, et al. Optimizing Nand irrigation inputs for corn based on N leaching and yield on a coarsetextured soil[J]. Journal of environmental quality, 1996,25: 982-992.

    [23] Fang Quanxiao, Ma Liwang, Yu Qiang, et al. Modeling nitrogen and water management effects in a wheat-maize double-cropping system[J]. Journal of Environmental Quality,2008, 37:2232-2242.

    [24] Michalczyk A, Kersebaum K C, Roelcke M, et al. Modelbased optimisation of nitrogen and water management for wheat-maize systems in the North China Plain[J]. Nutrient Cycling in Agro-ecosystems. 2014, 98(2): 203-222.

    [25] 王小彬,代快,趙全勝,等. 農田水氮關系及其協(xié)同管理[J]. 生態(tài)學報,2010,30(24):7001-7015.Wang Xiaobin, Dai Kuai, Zhao Quansheng, et al. Opinions on water-nitrogen relations and their synergism management[J]. Acta Ecologica Sinica, 2010, 30(24): 7001-7015. (in Chinese with English abstract)

    [26] 呂麗華,董志強,張經廷,等. 水氮對冬小麥-夏玉米產量及氮利用效應研究[J]. 中國農業(yè)科學,2014,47(19):3839-3849.Lü Lihua, Dong Zhiqiang, Zhang Jingting, et al. Effect of water and nitrogen on yield and nitrogen utilization of winter wheat and summer maize[J]. Scientia Agricultura Sinica,2014, 47(19): 3839-3849. (in Chinese with English abstract)

    [27] Fu Qiuping, Wang Quanjiu, Shen Xinlei, et al. Optimizing water and nitrogen inputs for winter wheat cropping system on the Loess Plateau, China[J]. Journal of Arid Land, 2014,6(2): 230-242.

    [28] Chen Deli, Suter H, Islam A, et al. Prospects of improving efficiency of fertiliser nitrogen in Australian agriculture: A review of enhanced efficiency fertilisers[J]. Australian Journal of Soil Research, 2008, 46, 289-301.

    [29] Wang Dong, Xu Zhenzhu, Zhao Junye, et al. Excessive nitrogen application decreases grain yield and increases nitrogen loss in a wheat-soil system[J]. Acta Agriculturae Scandinavica Section B-Soil & Plant Science, 2011, 61:681-692.

    [30] Jin Libin, Cui Haiyan, Li Bo, et al. Effects of integrated agronomic management practices on yield and nitrogen efficiency of summer maize in North China[J]. Field Crops Research, 2012, 134: 30-35.

    [31] Wang Xiangping, Huang Guanhua. Evaluation on the irrigation and fertilization management practices under the application of treated sewage water in Beijing, China[J].Agricultural Water Management, 2008, 95: 1011-1027.

    [32] Zhao Rongfang, Chen Xinping, Zhang Hailin, et al.Fertilization and nitrogen balance in a wheat/maize rotation system in North China[J]. Agronomy Journal, 2006, 98:938-945.

    Optimal management of water and nitrogen for farmland in North China Plain based on osculating value method and WHCNS model

    Xu Qiang, Li Zhoujing, Hu Kelin※, Li Baoguo
    (College of Resources and Environmental Sciences, China Agricultural University, Key Laboratory of Arable Land Conservation in North China, Ministry of Agriculture, Beijing100193,China)

    Optimizing water and nitrogen (N) management is important for high grain yield, resource-efficient and environment-friendly targets in North China Plain (NCP). In this study, the Soil Water Heat Carbon and Nitrogen Simulation(WHCNS) model was combined with the osculating value method to achieve this objective. A 2-year experiment involving 4 integrated management practices was conducted in Tai’an City, Shandong Province in the NCP. These management practices were designed as follows: 1) traditional farming practice (FP); 2) optimized combination of cropping and fertilization (OPT-1);3) practice for high yield (HY), which did not consider the cost of resource inputs to maximize grain yield; and 4) further optimized combination of cropping and fertilization (OPT-2), which was based on the HY practice. The dynamics of soil water content and soil nitrate concentration in different soil depths were monitored, crop dry matter and leaf area index at the key crop growth stages and yield were measured. These dataset was used to test the WHCNS model. After calibration and validation, the scenario analysis was conducted by the validated WHCNS model, which included 165 scenarios of winter wheat and 55 scenarios of summer maize with different water and N management combinations. The osculating value method was then employed to obtain the best management practices (BMPs). This method took into account agronomic, environmental and economic benefits. The results showed that the winter wheat yield increased firstly with the amount of irrigation and then kept steady when reached the maximum. There was no significant relationship between maize yields with irrigation due to the summer maize growth was synchronized with rains in NCP. The yield of winter wheat and summer maize increased with the amount of fertilizer firstly, and then remained stable when reached the maximum. The amount of water drainage and nitrogen leaching increased significantly with the irrigation amount and fertilizer application rates. The optimal irrigation amount and fertilization rates were 240 mm and 60 kg/hm2(based on N) for winter wheat, respectively. The fertilizer-N application at a rate of 90 kg/hm2with no irrigation was the BMPs for summer maize. Of all the scenarios, the BMPs could guarantee 97% of the maximum annual yield, higher water and N use efficiency and the highest value-cost ratio (VCR), whilst N leaching and gas loss reduced by 77% and 71%, respectively. Thus, this method can be used to optimize management of water and nitrogen in farmland effectively.

    irrigation; nitrogen; optimum; winter wheat-summer maize double-cropping system; WHCNS model;scenario analysis; osculating value method

    10.11975/j.issn.1002-6819.2017.14.021

    S152

    A

    1002-6819(2017)-14-0152-07

    徐 強,李周晶,胡克林,李保國. 基于密切值法和水氮管理模型的華北平原農田水氮優(yōu)化管理[J]. 農業(yè)工程學報,2017,33(14):152-158.

    10.11975/j.issn.1002-6819.2017.14.021 http://www.tcsae.org

    Xu Qiang, Li Zhoujing, Hu Kelin, Li Baoguo. Optimal management of water and nitrogen for farmland in North China Plain based on osculating value method and WHCNS model[J]. Transactions of the Chinese Society of Agricultural Engineering(Transactions of the CSAE), 2017, 33(14): 152-158. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2017.14.021 http://www.tcsae.org

    2017-01-13

    2017-06-10

    國家重點研發(fā)計劃項目資助(2016YFD0800102);長江學者和創(chuàng)新團隊發(fā)展計劃項目(IRT0412)

    徐 強,男,黑龍江哈爾濱人,博士生,主要從事資源環(huán)境系統(tǒng)模型應用。北京 中國農業(yè)大學資源與環(huán)境學院,100193。

    Email: qiangxu2014@163.com

    ※通信作者:胡克林,男,湖北鐘祥人,博士,教授,主要從事土壤空間變異及溶質運移的研究。北京 中國農業(yè)大學資源與環(huán)境學院,100193。

    Email: hukel@cau.edu.cn

    猜你喜歡
    水氮夏玉米施肥量
    蔬菜施肥量控制技術指南
    小麥收割之后 如何種植夏玉米才能高產
    夏玉米高產的關鍵栽培技術措施
    現代園藝(2017年19期)2018-01-19 02:50:04
    水氮耦合對煙株生長、產量及主要化學成分的影響
    水氮交互作用對棉花產量的影響研究
    水氮耦合對膜下滴灌棉花干物質積累的影響研究
    水果紅薯無土高產栽培的最佳施肥量研究
    省工棉2號最佳種植密度及最適施肥量的研究
    不同施肥量對煙草生長發(fā)育及經濟產量的影響
    膜下滴灌水氮耦合對棉花干物質積累和氮素吸收及水氮利用效率的影響
    欧美少妇被猛烈插入视频| 亚洲自偷自拍三级| 男女边摸边吃奶| 免费av不卡在线播放| 亚洲色图av天堂| 麻豆成人av视频| 赤兔流量卡办理| tube8黄色片| 国产欧美日韩一区二区三区在线 | 一边亲一边摸免费视频| 少妇丰满av| 成人毛片a级毛片在线播放| 亚洲成色77777| 在线精品无人区一区二区三 | av女优亚洲男人天堂| 一二三四中文在线观看免费高清| av在线天堂中文字幕| 国产亚洲午夜精品一区二区久久 | 国产成人91sexporn| 大陆偷拍与自拍| 国产探花极品一区二区| 内地一区二区视频在线| 热re99久久精品国产66热6| 亚洲精品日本国产第一区| 国产亚洲av嫩草精品影院| 一区二区三区乱码不卡18| 免费大片18禁| 青春草亚洲视频在线观看| 国产成年人精品一区二区| 国产伦精品一区二区三区视频9| 国产精品偷伦视频观看了| 在线观看av片永久免费下载| 亚洲欧美精品自产自拍| 欧美最新免费一区二区三区| 网址你懂的国产日韩在线| 免费大片黄手机在线观看| 一级毛片aaaaaa免费看小| 精品人妻一区二区三区麻豆| 久久久久久久大尺度免费视频| 好男人在线观看高清免费视频| 一级毛片我不卡| 啦啦啦中文免费视频观看日本| 少妇裸体淫交视频免费看高清| a级毛片免费高清观看在线播放| 亚洲内射少妇av| 国产精品国产三级国产av玫瑰| 亚洲欧美成人综合另类久久久| 少妇人妻一区二区三区视频| 日日啪夜夜撸| 我的老师免费观看完整版| 国产精品一区二区三区四区免费观看| 精品人妻熟女av久视频| 欧美性感艳星| av又黄又爽大尺度在线免费看| 日韩不卡一区二区三区视频在线| 日本色播在线视频| 99精国产麻豆久久婷婷| 尾随美女入室| 成人欧美大片| 国产精品国产av在线观看| 熟女电影av网| 亚洲精品日韩av片在线观看| 2018国产大陆天天弄谢| 插阴视频在线观看视频| 久久精品国产亚洲av天美| 人人妻人人看人人澡| 人体艺术视频欧美日本| 久久久久久伊人网av| 久久99热这里只频精品6学生| 丝袜脚勾引网站| 久热久热在线精品观看| 小蜜桃在线观看免费完整版高清| 91久久精品电影网| 爱豆传媒免费全集在线观看| 能在线免费看毛片的网站| 一级av片app| 久久这里有精品视频免费| 欧美区成人在线视频| 91aial.com中文字幕在线观看| 日韩av免费高清视频| 一二三四中文在线观看免费高清| 性色av一级| 精品视频人人做人人爽| 久久女婷五月综合色啪小说 | 一级片'在线观看视频| 好男人在线观看高清免费视频| 亚洲久久久久久中文字幕| 亚洲,欧美,日韩| 久久久亚洲精品成人影院| 丰满人妻一区二区三区视频av| 亚洲,一卡二卡三卡| 建设人人有责人人尽责人人享有的 | 精品久久久久久久久av| 最近2019中文字幕mv第一页| 久热这里只有精品99| 波野结衣二区三区在线| 日韩精品有码人妻一区| av在线播放精品| 少妇 在线观看| 天堂网av新在线| 久久久久久久久久人人人人人人| 你懂的网址亚洲精品在线观看| 国产高清不卡午夜福利| 欧美激情久久久久久爽电影| 白带黄色成豆腐渣| 久久鲁丝午夜福利片| 白带黄色成豆腐渣| 黄色日韩在线| 九九在线视频观看精品| av在线蜜桃| 久久热精品热| 国产一区二区三区av在线| av播播在线观看一区| 国产永久视频网站| 久久99精品国语久久久| 五月伊人婷婷丁香| 亚洲国产精品专区欧美| 黑人高潮一二区| 夫妻午夜视频| 国产综合精华液| 18禁在线无遮挡免费观看视频| 欧美成人a在线观看| 春色校园在线视频观看| a级毛色黄片| 亚洲成人一二三区av| 国产中年淑女户外野战色| 香蕉精品网在线| 亚洲精品一区蜜桃| 又黄又爽又刺激的免费视频.| 精品一区二区三卡| 91狼人影院| 亚洲色图av天堂| 在线天堂最新版资源| 久久久午夜欧美精品| 亚洲精品,欧美精品| 精品国产乱码久久久久久小说| av免费观看日本| 久久久久久久久大av| 国国产精品蜜臀av免费| 校园人妻丝袜中文字幕| 一个人观看的视频www高清免费观看| 丰满少妇做爰视频| av女优亚洲男人天堂| 干丝袜人妻中文字幕| 国产久久久一区二区三区| 别揉我奶头 嗯啊视频| 可以在线观看毛片的网站| 高清视频免费观看一区二区| .国产精品久久| 寂寞人妻少妇视频99o| 免费看a级黄色片| 亚洲国产av新网站| 一边亲一边摸免费视频| 啦啦啦中文免费视频观看日本| 久久国内精品自在自线图片| 最近的中文字幕免费完整| 久久精品夜色国产| www.色视频.com| 大码成人一级视频| 我要看日韩黄色一级片| 国产高清国产精品国产三级 | 午夜免费男女啪啪视频观看| 岛国毛片在线播放| 精品国产一区二区三区久久久樱花 | 18禁裸乳无遮挡动漫免费视频 | 日日啪夜夜爽| 国产日韩欧美在线精品| 国产探花极品一区二区| 美女脱内裤让男人舔精品视频| 汤姆久久久久久久影院中文字幕| 国产伦精品一区二区三区视频9| 国产精品三级大全| 91久久精品国产一区二区成人| 一本一本综合久久| 免费看av在线观看网站| 如何舔出高潮| 亚洲人成网站高清观看| 日韩成人伦理影院| 欧美日韩综合久久久久久| 日韩av在线免费看完整版不卡| 亚洲国产av新网站| av卡一久久| 香蕉精品网在线| 内射极品少妇av片p| 2021天堂中文幕一二区在线观| 欧美少妇被猛烈插入视频| 国产淫片久久久久久久久| 国产 一区 欧美 日韩| 欧美日韩综合久久久久久| 亚洲欧美成人综合另类久久久| 亚洲色图综合在线观看| 亚洲成人一二三区av| 亚洲欧美成人精品一区二区| 国产av码专区亚洲av| 午夜福利在线在线| 亚洲国产精品专区欧美| 久久ye,这里只有精品| 大片免费播放器 马上看| 日韩 亚洲 欧美在线| 国产爽快片一区二区三区| 亚洲丝袜综合中文字幕| 看十八女毛片水多多多| 欧美区成人在线视频| freevideosex欧美| 最近最新中文字幕免费大全7| 搡老乐熟女国产| 亚洲精品久久午夜乱码| 国产精品熟女久久久久浪| 干丝袜人妻中文字幕| 免费少妇av软件| 国产成人精品福利久久| 边亲边吃奶的免费视频| 亚洲欧洲国产日韩| 激情五月婷婷亚洲| 欧美成人一区二区免费高清观看| 国产在视频线精品| 亚洲成人一二三区av| 男女无遮挡免费网站观看| 色综合色国产| 只有这里有精品99| 自拍欧美九色日韩亚洲蝌蚪91 | 内地一区二区视频在线| 久热这里只有精品99| 国产有黄有色有爽视频| av播播在线观看一区| 欧美精品国产亚洲| 91久久精品国产一区二区成人| 下体分泌物呈黄色| 内地一区二区视频在线| 亚洲国产精品999| 欧美精品人与动牲交sv欧美| 国产极品天堂在线| 日韩中字成人| 久久久久久国产a免费观看| 国产免费一区二区三区四区乱码| 国产69精品久久久久777片| 99热这里只有是精品50| 成人免费观看视频高清| 两个人的视频大全免费| 日韩伦理黄色片| 综合色av麻豆| 国产老妇伦熟女老妇高清| 99热这里只有是精品在线观看| 真实男女啪啪啪动态图| 亚洲国产精品999| 国内精品美女久久久久久| av网站免费在线观看视频| 在线观看一区二区三区| 欧美国产精品一级二级三级 | 亚洲精品第二区| 国产精品嫩草影院av在线观看| 麻豆乱淫一区二区| 亚洲人成网站在线播| 一级爰片在线观看| 亚洲色图av天堂| 成人国产av品久久久| 久久99蜜桃精品久久| 国产国拍精品亚洲av在线观看| av播播在线观看一区| 欧美xxxx黑人xx丫x性爽| 成人美女网站在线观看视频| 国产高潮美女av| 青青草视频在线视频观看| 人妻系列 视频| 亚洲精华国产精华液的使用体验| 日韩精品有码人妻一区| 少妇 在线观看| 免费大片黄手机在线观看| 丰满人妻一区二区三区视频av| 狂野欧美激情性bbbbbb| 亚洲色图av天堂| 一级片'在线观看视频| 麻豆乱淫一区二区| 波野结衣二区三区在线| 日韩av免费高清视频| 99热这里只有是精品在线观看| 免费av观看视频| av网站免费在线观看视频| 国产成人免费观看mmmm| 青春草亚洲视频在线观看| 亚洲国产最新在线播放| 美女高潮的动态| 女的被弄到高潮叫床怎么办| 午夜免费观看性视频| 中文欧美无线码| 一个人观看的视频www高清免费观看| 国产日韩欧美在线精品| 国产黄频视频在线观看| 18禁在线无遮挡免费观看视频| 国产探花在线观看一区二区| 亚洲熟女精品中文字幕| 成年人午夜在线观看视频| 亚洲成人av在线免费| 亚洲色图综合在线观看| 亚洲欧美日韩东京热| 日韩强制内射视频| 国产成人午夜福利电影在线观看| 18禁裸乳无遮挡动漫免费视频 | 一区二区av电影网| 97精品久久久久久久久久精品| 日韩成人伦理影院| 亚洲av免费在线观看| 欧美成人精品欧美一级黄| 三级经典国产精品| 青青草视频在线视频观看| 小蜜桃在线观看免费完整版高清| 亚洲国产欧美人成| 国产探花在线观看一区二区| 久久久国产一区二区| 丝袜喷水一区| 黑人高潮一二区| 97在线人人人人妻| 国产伦精品一区二区三区视频9| 国产精品伦人一区二区| av在线观看视频网站免费| 人妻少妇偷人精品九色| 久久热精品热| 最近的中文字幕免费完整| av在线蜜桃| 日本一二三区视频观看| 全区人妻精品视频| 久久久久九九精品影院| av.在线天堂| 国产成人a∨麻豆精品| 亚洲国产日韩一区二区| 国产av国产精品国产| 欧美激情国产日韩精品一区| 精品少妇久久久久久888优播| 国产精品国产av在线观看| 久久久国产一区二区| av在线老鸭窝| 精品午夜福利在线看| 国产综合懂色| 我的女老师完整版在线观看| av在线app专区| 国产精品av视频在线免费观看| 免费看日本二区| 街头女战士在线观看网站| 日韩一区二区三区影片| 国产黄频视频在线观看| 亚洲精品中文字幕在线视频 | 久久久久久久亚洲中文字幕| 久久99蜜桃精品久久| 好男人视频免费观看在线| 午夜精品一区二区三区免费看| 欧美亚洲 丝袜 人妻 在线| 免费观看无遮挡的男女| 夜夜爽夜夜爽视频| 久久99蜜桃精品久久| 午夜亚洲福利在线播放| 伊人久久国产一区二区| 神马国产精品三级电影在线观看| 在线观看免费高清a一片| 精品久久久噜噜| 中文欧美无线码| 丰满少妇做爰视频| 欧美日本视频| 在线观看一区二区三区| 丝袜脚勾引网站| 亚洲国产日韩一区二区| 国产亚洲91精品色在线| 亚洲成人一二三区av| 日本免费在线观看一区| 国产精品精品国产色婷婷| 国产乱人偷精品视频| 日韩人妻高清精品专区| 亚洲aⅴ乱码一区二区在线播放| 成人美女网站在线观看视频| 国产午夜精品久久久久久一区二区三区| 一二三四中文在线观看免费高清| 亚洲精品乱码久久久v下载方式| av网站免费在线观看视频| 久久97久久精品| 夜夜爽夜夜爽视频| 黄片wwwwww| 亚洲成人久久爱视频| 人妻制服诱惑在线中文字幕| 爱豆传媒免费全集在线观看| 免费看av在线观看网站| 免费观看av网站的网址| 青春草亚洲视频在线观看| 国产亚洲av片在线观看秒播厂| 一区二区三区精品91| 午夜福利视频1000在线观看| 免费av观看视频| 成人国产av品久久久| 日韩一区二区三区影片| 精品国产露脸久久av麻豆| 偷拍熟女少妇极品色| 又黄又爽又刺激的免费视频.| 日本与韩国留学比较| 69人妻影院| 成年女人在线观看亚洲视频 | 亚洲av不卡在线观看| 在线天堂最新版资源| 校园人妻丝袜中文字幕| 大又大粗又爽又黄少妇毛片口| h日本视频在线播放| 国产高清国产精品国产三级 | 日本一本二区三区精品| 3wmmmm亚洲av在线观看| 亚洲欧美中文字幕日韩二区| 三级经典国产精品| 亚洲最大成人av| 日本一本二区三区精品| 亚洲成人中文字幕在线播放| 老司机影院成人| 99热6这里只有精品| 大香蕉久久网| 精品久久久久久久久av| 尤物成人国产欧美一区二区三区| 免费高清在线观看视频在线观看| 一级毛片我不卡| 内射极品少妇av片p| 精品久久国产蜜桃| 中文字幕免费在线视频6| 永久网站在线| 一个人观看的视频www高清免费观看| 波多野结衣巨乳人妻| 免费av毛片视频| 国产精品伦人一区二区| 国产乱人视频| 18禁裸乳无遮挡免费网站照片| 久久久久久久久久成人| 国产黄a三级三级三级人| 亚洲图色成人| 国产精品福利在线免费观看| 在现免费观看毛片| 免费大片黄手机在线观看| 久久精品综合一区二区三区| 中文在线观看免费www的网站| 中文字幕人妻熟女乱码| 在线免费观看不下载黄p国产| 免费观看av网站的网址| 在线观看免费高清a一片| 熟女av电影| 国产男女内射视频| 天天躁狠狠躁夜夜躁狠狠躁| 一区在线观看完整版| 九色亚洲精品在线播放| 男女免费视频国产| 亚洲欧美一区二区三区久久| 人人妻人人爽人人添夜夜欢视频| 国产激情久久老熟女| 在线观看人妻少妇| 国产国语露脸激情在线看| 超碰97精品在线观看| 久久久亚洲精品成人影院| 老司机深夜福利视频在线观看 | 一级黄片播放器| 亚洲欧美精品自产自拍| 咕卡用的链子| 国产精品.久久久| 久久久久久久久免费视频了| 国产精品国产av在线观看| av又黄又爽大尺度在线免费看| 亚洲美女搞黄在线观看| 观看av在线不卡| 老司机亚洲免费影院| 高清视频免费观看一区二区| 男女边吃奶边做爰视频| 国产精品麻豆人妻色哟哟久久| 精品一区二区三区四区五区乱码 | 国产老妇伦熟女老妇高清| 午夜福利在线免费观看网站| 精品少妇久久久久久888优播| 免费看不卡的av| 免费久久久久久久精品成人欧美视频| 乱人伦中国视频| 国产精品一国产av| 国产日韩欧美视频二区| 亚洲国产精品一区二区三区在线| 毛片一级片免费看久久久久| 在线精品无人区一区二区三| 免费av中文字幕在线| 国产xxxxx性猛交| 一级毛片 在线播放| 999久久久国产精品视频| 成人手机av| 亚洲欧美精品综合一区二区三区| 黄频高清免费视频| 啦啦啦视频在线资源免费观看| a级片在线免费高清观看视频| 久久久精品94久久精品| avwww免费| 亚洲av福利一区| 无遮挡黄片免费观看| 日韩一卡2卡3卡4卡2021年| 狠狠精品人妻久久久久久综合| 中文欧美无线码| 丰满饥渴人妻一区二区三| 精品福利永久在线观看| 99久国产av精品国产电影| 国产av精品麻豆| 卡戴珊不雅视频在线播放| 国产成人精品福利久久| 成人漫画全彩无遮挡| 亚洲第一av免费看| 18禁裸乳无遮挡动漫免费视频| 亚洲精品国产av蜜桃| 国产精品偷伦视频观看了| 亚洲熟女毛片儿| 日韩av不卡免费在线播放| 日韩免费高清中文字幕av| 精品酒店卫生间| 啦啦啦在线免费观看视频4| 亚洲精华国产精华液的使用体验| 久久久精品94久久精品| 在线观看免费视频网站a站| 成年女人毛片免费观看观看9 | 亚洲国产中文字幕在线视频| 国产成人系列免费观看| 色婷婷av一区二区三区视频| 日韩中文字幕欧美一区二区 | 国产精品成人在线| 一区二区三区激情视频| 视频在线观看一区二区三区| 97人妻天天添夜夜摸| 国产精品久久久久久人妻精品电影 | 欧美人与性动交α欧美软件| 欧美激情高清一区二区三区 | 精品国产乱码久久久久久男人| 老司机深夜福利视频在线观看 | 晚上一个人看的免费电影| 中文字幕人妻熟女乱码| 欧美成人精品欧美一级黄| 婷婷色综合大香蕉| 久久人人97超碰香蕉20202| 操出白浆在线播放| 999精品在线视频| 国产精品av久久久久免费| 亚洲欧美一区二区三区久久| 欧美激情高清一区二区三区 | 黄频高清免费视频| 国产精品嫩草影院av在线观看| av.在线天堂| 久久精品亚洲熟妇少妇任你| 51午夜福利影视在线观看| 国产精品嫩草影院av在线观看| 人成视频在线观看免费观看| 欧美日韩精品网址| 老鸭窝网址在线观看| 九色亚洲精品在线播放| 男女边摸边吃奶| 一边摸一边做爽爽视频免费| 91国产中文字幕| 蜜桃国产av成人99| av免费观看日本| 成年av动漫网址| 亚洲伊人久久精品综合| 国产精品麻豆人妻色哟哟久久| 国产男人的电影天堂91| 人妻 亚洲 视频| 国产日韩一区二区三区精品不卡| 亚洲,欧美,日韩| 国产1区2区3区精品| 亚洲人成77777在线视频| 国产免费一区二区三区四区乱码| 丝袜美腿诱惑在线| 色婷婷久久久亚洲欧美| 久久免费观看电影| 麻豆精品久久久久久蜜桃| 91国产中文字幕| 十分钟在线观看高清视频www| 丝袜喷水一区| 国产成人啪精品午夜网站| 亚洲视频免费观看视频| 天堂8中文在线网| 国产在线免费精品| 好男人视频免费观看在线| 欧美最新免费一区二区三区| 高清欧美精品videossex| 亚洲一码二码三码区别大吗| 少妇人妻精品综合一区二区| 18在线观看网站| 超碰97精品在线观看| 日韩精品有码人妻一区| 国产亚洲av片在线观看秒播厂| 国产一区亚洲一区在线观看| 日韩一区二区视频免费看| 久久ye,这里只有精品| 欧美日本中文国产一区发布| 午夜福利视频精品| 你懂的网址亚洲精品在线观看| 久久久久网色| 亚洲av中文av极速乱| 晚上一个人看的免费电影| 菩萨蛮人人尽说江南好唐韦庄| 老熟女久久久| 国产成人欧美在线观看 | 久热爱精品视频在线9| 中文欧美无线码| 九九爱精品视频在线观看| 18禁裸乳无遮挡动漫免费视频| 日韩精品免费视频一区二区三区| 久久ye,这里只有精品| 国产一区二区三区综合在线观看| 91精品伊人久久大香线蕉| 中文字幕另类日韩欧美亚洲嫩草| 五月天丁香电影| 在线观看免费日韩欧美大片| 纯流量卡能插随身wifi吗| 亚洲美女视频黄频| 搡老岳熟女国产| 欧美精品av麻豆av| 19禁男女啪啪无遮挡网站| 国产成人a∨麻豆精品| 免费女性裸体啪啪无遮挡网站| 9热在线视频观看99| 18在线观看网站| 老汉色∧v一级毛片| 国产欧美日韩综合在线一区二区| 极品少妇高潮喷水抽搐| 欧美中文综合在线视频| 国产精品二区激情视频|