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

    華北兩熟區(qū)農(nóng)田磷流失風(fēng)險預(yù)警研究*

    2019-08-31 02:22:30冀宏杰張懷志張認(rèn)連徐愛國田昌玉
    關(guān)鍵詞:投入量磷素輪作

    冀宏杰, 張懷志, 張認(rèn)連, 徐愛國, 田昌玉

    華北兩熟區(qū)農(nóng)田磷流失風(fēng)險預(yù)警研究*

    冀宏杰, 張懷志, 張認(rèn)連, 徐愛國, 田昌玉

    (中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所/農(nóng)業(yè)農(nóng)村部植物營養(yǎng)與養(yǎng)分循環(huán)重點開放實驗室 北京 100081)

    輪作是長期以來在各地形成的相對固定的作物組合模式, 同時也是農(nóng)田養(yǎng)分平衡狀況計算的基本單元, 多年累積形成的養(yǎng)分盈虧則直接影響環(huán)境。在植物必需的3種大量營養(yǎng)元素中, 對環(huán)境影響最大的營養(yǎng)元素是磷和氮, 與氮素相比, 磷素循環(huán)途徑則相對穩(wěn)定, 沒有氣態(tài)損失, 磷的表觀養(yǎng)分平衡與環(huán)境污染的相關(guān)性高, 且計算方便, 因而本文選擇磷素作為環(huán)境風(fēng)險預(yù)警指示元素。為探究簡易、有效的輪作環(huán)境效應(yīng)預(yù)警方法, 從篩選適宜預(yù)警的關(guān)鍵元素和環(huán)境風(fēng)險產(chǎn)生的關(guān)鍵環(huán)節(jié)入手, 提出適合我國小農(nóng)戶現(xiàn)實國情的、基于農(nóng)田地塊級別的預(yù)警指標(biāo)體系, 并以華北兩熟區(qū)5種輪作類型、38個地塊的調(diào)查數(shù)據(jù)為例進行了驗證。該預(yù)警體系包括3級評價指標(biāo)(周年磷平衡、耕層土壤質(zhì)地、周年磷肥運籌)共6種預(yù)警級別并用3種色系進行標(biāo)識(深綠、淺綠、淺黃、深黃、淺紅、深紅)。該預(yù)警方法的特點為: 研究對象為整個輪作周期而非單茬作物; 磷平衡既考慮周年的總量均衡(即整個輪作周期磷素投入、產(chǎn)出平衡狀況), 又考慮不同茬次間的科學(xué)分配(即磷素統(tǒng)籌), 且二者均采用相對平衡值而非絕對平衡值, 從而使得不同輪作類型之間的磷平衡可以進行統(tǒng)一量化和比較; 土壤質(zhì)地簡化為沙土、壤土、黏土3種。驗證結(jié)果表明, 磷肥用量、運籌合理、環(huán)境友好(深綠)的地塊僅占10.53%; 而高環(huán)境風(fēng)險(深紅、淺紅)的輪作地塊達(dá)57.89%, 即一半以上的地塊目前的輪作方式存在高或極高的環(huán)境風(fēng)險。該預(yù)警方法具有科學(xué)、實用、簡便, 所需數(shù)據(jù)容易獲得的特點, 可用于農(nóng)戶對自己地塊磷流失風(fēng)險進行判斷, 也可用于國家、區(qū)域或地方農(nóng)業(yè)主管部門進行種植業(yè)結(jié)構(gòu)調(diào)整時參考。

    兩熟區(qū); 輪作; 農(nóng)田養(yǎng)分; 磷; 流失風(fēng)險; 預(yù)警系統(tǒng); 華北

    兩熟區(qū)是我國糧食作物和經(jīng)濟作物生產(chǎn)供應(yīng)的重要區(qū)域, 如華北平原、黃土高原、江淮丘陵等區(qū)域[1]。在漫長的歷史發(fā)展過程中, 兩熟區(qū)形成了眾多的種植類型, 特別是輪作, 栽培面積廣泛。其中既有冬小麥(L.)-夏玉米(L.)這樣典型的糧食作物組合, 也有冬小麥-夏大豆[(L.) Merr.]、小麥-花生(L.)、小麥-棉花(spp.)等糧食作物與經(jīng)濟作物組合, 甚至還有大蒜(L.)-棉花等具有區(qū)域特色的作物組合。多樣的輪作類型充分利用了當(dāng)?shù)毓鉁責(zé)豳Y源, 對于國家大宗糧油作物、特色作物供給, 保障國家糧食安全、促進農(nóng)民增收、農(nóng)業(yè)增效發(fā)揮了不可替代的作用。

    從養(yǎng)分與環(huán)境的關(guān)系來看, 當(dāng)某種養(yǎng)分的投入大于支出時, 就會造成該養(yǎng)分盈余, 反之則造成養(yǎng)分虧缺, 短期、適度盈余或虧缺不會影響環(huán)境, 也不會影響作物生長, 而長期大量盈余或虧缺則不僅會影響產(chǎn)量、品質(zhì), 降低肥料利用效率, 還會引起諸多環(huán)境問題[2]。尤其是磷素, 磷盈余引起的土壤磷含量提高會增大通過地表徑流和滲漏量進入環(huán)境的量[3-4], 從而會對土壤、水環(huán)境造成負(fù)面影響, 如地表水富營養(yǎng)化等問題[5-7]。因此, 建立基于地塊的預(yù)警系統(tǒng), 進行農(nóng)田磷素流失預(yù)警對優(yōu)化施肥措施, 提高肥料利用效率及保護環(huán)境意義重大。

    一個區(qū)域輪作制度是多年逐漸形成的相對穩(wěn)定的作物組合, 由于累積效應(yīng)往往造成不同輪作系統(tǒng)之間的養(yǎng)分平衡狀況差異巨大。輪作是作物種植的基本單元, 也是決定周年養(yǎng)分投入的基本單元。由于輪作對于周年養(yǎng)分乃至多年養(yǎng)分累積(盈余或虧缺)的決定性作用, 故地塊的多年養(yǎng)分平衡狀況可通過周年輪作的養(yǎng)分盈余評估計算[2]。因此以周年輪作系統(tǒng)為研究的基本單元。一方面, 由于組成某種輪作系統(tǒng)的作物是固定的, 而作物對養(yǎng)分吸收的偏好及農(nóng)民生產(chǎn)習(xí)慣等因素, 長期固定輪作會導(dǎo)致土壤中某種養(yǎng)分的過?;騾T乏, 以及病蟲害加劇、土壤鹽漬化、酸化等問題, 不利于農(nóng)田生態(tài)系統(tǒng)穩(wěn)定, 最終造成產(chǎn)量降低。另一方面, 即使同一種輪作系統(tǒng), 由于不同農(nóng)戶管理水平的差異, 不同地塊間養(yǎng)分平衡狀況有極大差異。因而預(yù)警系統(tǒng)既需要考慮輪作的共性特征, 又需要考慮農(nóng)戶個體極端風(fēng)險, 基于地塊進行預(yù)警系統(tǒng)研究對于環(huán)境風(fēng)險管控更具科學(xué)價值。而我國尚缺乏針對地塊的從環(huán)境風(fēng)險角度考慮的預(yù)警系統(tǒng)。

    目前的環(huán)境風(fēng)險預(yù)警往往針對的重金屬指標(biāo)[8-9]、畜禽養(yǎng)殖[10], 或從社會經(jīng)濟等宏觀角度出發(fā)[11], 尚缺乏農(nóng)業(yè)環(huán)境風(fēng)險預(yù)警系統(tǒng)。據(jù)農(nóng)業(yè)調(diào)查, 農(nóng)業(yè)面源污染在總污染中占比嚴(yán)重, 據(jù)2010年《第一次全國污染源普查公報》結(jié)果, 我國農(nóng)業(yè)污染源(不包括典型地區(qū)農(nóng)村生活源)的總氮(TN)和總磷(TP)排放分別為270.46萬t和28.47萬t, 占地表水體污染總負(fù)荷的57.2%和67.4%。可見, 農(nóng)業(yè)面源污染已經(jīng)成為當(dāng)前地表水體污染的主要來源[12-13]。

    在植物必需的3種大量營養(yǎng)元素中, 對環(huán)境影響最大的營養(yǎng)元素是磷和氮, 這可以從歷年《中國環(huán)境狀況公報》公布的地表水主要污染物類型看出[14-16]。然而由于在農(nóng)田生態(tài)系統(tǒng)中氮素的循環(huán)途徑非常多, 氮素?fù)]發(fā)與徑流、淋溶等損失的影響因素多, 造成氮表觀平衡值往往變化比較大; 與氮素相比, 磷素循環(huán)途徑則相對穩(wěn)定, 沒有氣態(tài)損失, 磷的表觀養(yǎng)分平衡與環(huán)境污染的相關(guān)性高, 且計算方便, 因而本文選擇磷素作為環(huán)境風(fēng)險預(yù)警指示元素。

    施肥量既關(guān)系到農(nóng)民生產(chǎn)成本, 又關(guān)系到農(nóng)業(yè)環(huán)境健康。近些年來, 我國農(nóng)業(yè)面源污染形勢日趨嚴(yán)峻, 對農(nóng)業(yè)環(huán)境保護空前重視, “綠水青山就是金山銀山”理念日益深入人心, 亟需開發(fā)一種針對地塊種植磷素流失風(fēng)險進行預(yù)警的簡便易行的方法, 既便于國家、區(qū)域或地方農(nóng)業(yè)主管部門進行種植業(yè)結(jié)構(gòu)調(diào)整, 也利于農(nóng)戶自身判斷與決策。本研究從分析農(nóng)田養(yǎng)分平衡關(guān)鍵元素入手, 通過篩選農(nóng)田地塊適宜預(yù)警的關(guān)鍵元素(磷)、環(huán)境風(fēng)險產(chǎn)生的關(guān)鍵環(huán)節(jié)入手, 提出適合我國小農(nóng)戶現(xiàn)實國情的、基于農(nóng)田地塊級別的預(yù)警指標(biāo)體系, 并以華北兩熟區(qū)5種輪作類型、38個地塊的調(diào)查數(shù)據(jù)為例進行了驗證。

    1 研究方法

    1.1 輪作系統(tǒng)農(nóng)田磷流失風(fēng)險預(yù)警評價指標(biāo)體系的建立

    農(nóng)田磷素流失風(fēng)險首先取決于輪作周期內(nèi)磷素投入與輸出量的差異(即: 數(shù)量因素), 也取決于周年內(nèi)不同茬次間的比例(即: 運籌因素), 還要考慮流失過程的控制因素(即: 過程因素)。因而本文著重從以上3個方面選擇關(guān)鍵因素形成預(yù)警評價指標(biāo)體系。

    數(shù)量因素: 農(nóng)田磷素平衡是廣為應(yīng)用的指標(biāo), 主要用于土壤肥力評估、肥料利用效率、肥料量推薦、農(nóng)學(xué)閾值等方面[17-20], 而磷素平衡與磷素流失關(guān)系、環(huán)境閾值的研究相對較少。

    過程因素: 磷在土壤中的移動與土壤質(zhì)地有關(guān), 土壤質(zhì)地是影響磷流失的重要因素之一。一般來說, 質(zhì)地細(xì)的土壤的磷吸附能力強。在美國兩種不同黏粒含量的土壤研究發(fā)現(xiàn), 黏粒含量32%的土壤造成的磷素流失遠(yuǎn)遠(yuǎn)小于黏粒含量為5%的土壤[21]; 中國紫色土研究發(fā)現(xiàn)0~15 cm層中0.001 mm黏粒團聚度與徑流泥沙中全磷、速效磷含量顯著正相關(guān)[22]。土壤質(zhì)地主要取決于母質(zhì), 也受到地形等因素影響。

    運籌因素: 磷肥總量和磷肥在不同茬次間分配是影響磷素流失的重要因子。根據(jù)土壤磷有效性變化規(guī)律和輪作作物的磷養(yǎng)分吸收規(guī)律, 對輪作周期的每種作物的磷進行科學(xué)分配。研究表明, 一年兩熟地區(qū)的旱地, 磷肥重點施用在冬季作物上(冬重夏輕)[23], 而這恰恰在當(dāng)前肥料運籌中被忽視(詳見本文農(nóng)戶調(diào)查結(jié)果)。磷素在不同季節(jié)(或茬次)間分配的科學(xué)性與否, 在一定程度上決定了作物對土壤磷素吸收利用的程度, 并進一步?jīng)Q定其環(huán)境污染概率。

    基于以上3類因素, 確定了預(yù)警系統(tǒng)采用3級指標(biāo)評價體系, 分別為: 周年磷平衡、土壤質(zhì)地、周年磷肥運籌(秋播與夏播作物磷肥投入之比, 即“P秋/P夏”)。一級指標(biāo)為周年磷平衡狀況, 該指標(biāo)決定了不同的色系[綠色(環(huán)境相對友好)、紅色(環(huán)境風(fēng)險預(yù)警)、黃色(環(huán)境風(fēng)險中等)], 在此基礎(chǔ)上再依據(jù)土壤質(zhì)地(二級指標(biāo))與磷肥在不同茬次間分配情況(三級指標(biāo))細(xì)分, 共分為6種預(yù)警級別(深綠、淺綠、淺黃、深黃、淺紅、深紅, 深和淺表示各級別的程度)。根據(jù)不同土壤質(zhì)地對磷盈余的保持能力(黏土>壤土>沙土)分別將沙土、壤土、黏土的最適宜磷素平衡定為0.9~1.1、0.8~1.2和0.7~1.3。

    1.2 風(fēng)險預(yù)警實施步驟

    首先分別測定、計算地塊的周年磷平衡、耕層土壤質(zhì)地與周年磷肥運籌, 然后將所得結(jié)果填入系統(tǒng)預(yù)警表相應(yīng)位置即可自行得到該地塊、當(dāng)前輪作類型的預(yù)警級別。具體步驟如下:

    1)計算地塊的周年磷平衡, 即農(nóng)田磷素周年投入與產(chǎn)出之比:

    周年磷平衡=兩茬作物的磷投入量/兩茬作物的磷帶走量 (1)

    投入包括一年兩茬作物投入的基肥和追肥中有機肥、化肥的磷含量, 以肥料施用量乘以該肥料磷含量(P2O5)計算; 帶走量包括一年兩茬作物收獲攜出量, 以單位產(chǎn)量帶走磷養(yǎng)分量(P2O5)乘以產(chǎn)量計算。

    2)測定目標(biāo)地塊的耕層土壤質(zhì)地

    采用土壤質(zhì)地速測方法, 即手測法判斷耕層土壤的質(zhì)地, 分為砂質(zhì)、壤質(zhì)、黏質(zhì)3種[24]。

    3)計算秋播作物與夏播作物磷肥投入比

    計算秋播作物與夏播作物磷肥投入比(P秋/P夏), 該步驟的科學(xué)依據(jù)為根據(jù)土壤磷循環(huán)的生物地球化學(xué)特征, 周年中磷肥應(yīng)重點施用在秋播作物上以提高磷肥利用效率[23], 計算公式如下:

    P秋/P夏=秋播作物磷投入量/夏播作物磷投入量 (2)

    秋播作物磷投入量和夏播作物投入量分別包括一年兩茬作物所投入的基肥和追肥中有機肥、化肥的磷含量, 以肥料施用量乘以該肥料磷含量(P2O5)計算。

    4)根據(jù)以上步驟計算結(jié)果, 依據(jù)表1判斷地塊農(nóng)田的磷預(yù)警級別。

    表1 一年兩作輪作農(nóng)田磷流失預(yù)警指標(biāo)體系表1)

    1)用綠、黃、紅3種色系代表預(yù)警級別的高低, 紅色代表風(fēng)險等級最高, 綠色代表風(fēng)險最低, 黃色風(fēng)險居中; 深淺表示級別的程度。1) The alert level is shown by the colors. The red means high phosphorus loss risk; the green means low phosphorus loss risk; the yellow means medium phosphorus loss risk. Deep and light mean degree of the level.

    用綠、黃、紅3種色系代表預(yù)警級別, 紅色代表磷流失風(fēng)險等級最高, 綠色有環(huán)保之含義代表磷流失風(fēng)險最低, 而黃色表示磷流失風(fēng)險居中。色系的選擇由周年磷平衡值決定, 每種土壤質(zhì)地均包括此3種色系: 綠色系列代表周年磷投入量比較合理; 黃色系列代表周年磷投入量偏少, 需增加磷肥投入量; 紅色系列代表周年磷投入量偏多, 需減少磷肥投入量。

    每種色系又根據(jù)周年磷肥運籌細(xì)分為深淺兩級, 共計6種顏色, 分別代表6種預(yù)警級別, 各顏色所代表的含義如下: 深綠代表磷肥施用合理、茬次間分配合理、環(huán)境友好。淺綠代表磷肥周年投入總量合理, 但茬次間的分配不合理, 需調(diào)整磷肥的茬次間分配; 環(huán)境風(fēng)險較低。淺黃代表磷肥周年投入偏少, 但茬次間的分配合理, 需增加磷肥的周年總投入量; 環(huán)境風(fēng)險較低。深黃代表磷肥周年投入偏少, 需調(diào)整季節(jié)分配, 需增加磷肥的周年總投入量; 環(huán)境風(fēng)險中等。淺紅代表磷肥投入過量, 季節(jié)分配合理, 需減少磷肥的周年投入量; 環(huán)境風(fēng)險高。深紅代表磷投入過量, 且茬次間分配不合理, 需減少磷投入量, 且調(diào)整茬次間的分配; 環(huán)境風(fēng)險極高。

    1.3 應(yīng)用驗證

    采用本課題組調(diào)查數(shù)據(jù)對該方法體系進行說明, 并運用得到的結(jié)果對該區(qū)域的輪作風(fēng)險進行分析。數(shù)據(jù)來源為課題組2015年在山東、河北兩個省6個地區(qū)的11個縣市共38個農(nóng)戶地塊的調(diào)查結(jié)果, 主要包括華北地區(qū)5種輪作類型: 冬小麥-夏玉米、冬小麥-夏大豆、小麥-花生、小麥-棉花、大蒜-棉花; 調(diào)查內(nèi)容包括地塊位置、輪作情況、作物產(chǎn)量、施肥、耕層土壤質(zhì)地等信息。

    調(diào)查取樣的原則采用“隨機等比例”調(diào)查, 即首先按照各縣統(tǒng)計數(shù)據(jù)中的各種作物播種面積確定各縣的取樣點數(shù)量, 然后在各縣該種作物的集中種植區(qū)重點取樣, 而具體取樣位置采用隨機原則, 以最大限度地保證取樣的代表性。編號采用“縣名+順序號”命名法。

    2 結(jié)果與分析

    在華北兩熟區(qū)38個農(nóng)戶地塊上的周年磷平衡、秋播與夏播磷肥投入比、耕層土壤質(zhì)地及根據(jù)本文建立的預(yù)警方法得到預(yù)警結(jié)果如表2所示。

    2.1 不同輪作類型的周年磷平衡

    結(jié)果表明, 調(diào)查的38個地塊中, 兩茬作物的周年磷平衡值大于1的有30個, 占總樣本量的79%, ≤1的有8個, 占總樣本量的21%, 說明絕大部分地塊磷投入量超過帶走量, 為盈余狀態(tài); 按照輪作類型統(tǒng)計, 調(diào)查的5種輪作類型中, 周年磷平衡平均值從大到小依次為: 冬小麥-夏玉米(1.66)、冬小麥-花生(1.60)、大蒜-棉花(1.00)、冬小麥-大豆(0.85)、冬小麥-棉花(0.75)??梢钥闯? 種植面積最大、調(diào)查樣本量也最多的冬小麥-夏玉米輪作類型的周年磷平衡平均值恰恰也是最大, 表明當(dāng)前華北地區(qū)冬小麥-玉米輪作總體上磷肥投入過量, 普遍超過了其需求量造成磷盈余。這也與土壤速效磷含量持續(xù)上升的諸多證據(jù)相互印證[25-26]。廣大區(qū)域內(nèi)種植面積最大的大宗糧食作物磷肥盈余如此嚴(yán)重, 不僅造成嚴(yán)重了肥料資源浪費, 也會有潛在的日益增加的環(huán)境風(fēng)險, 亟需加強磷肥管理。

    表2 華北不同輪作系統(tǒng)的周年磷平衡、秋播作物與夏播作物磷肥投入比、土壤質(zhì)地及預(yù)警結(jié)果

    1)周年磷肥運籌指秋播與夏播磷肥投入比; 2)“#DIV/0!”代表夏播作物磷投入量為零, 此種情況按秋播作物與夏播作物磷肥投入比大于1計。1) The annual phosphorus allocation is the phosphorus input ratio of autumn-sowing crop to summer-sowing crop. 2) “#DIV/0!” means that the amount of phosphorus input of summer-sowing crop is zero, the phosphorus input ratio of autumn-sowing crop to summer-sowing crop is greater than 1.

    統(tǒng)計分析還表明, 同一輪作類型內(nèi)部, 不同農(nóng)戶之間的周年磷平衡值差異巨大, 以樣本量最大的兩種輪作類型為例: 冬小麥-夏玉米輪作最大值為2.80, 最小值僅為0.70, 兩者相差3倍, 標(biāo)準(zhǔn)差為0.58; 冬小麥-花生最大值為2.30, 最小值僅為1.00, 兩者相差1.3倍, 標(biāo)準(zhǔn)差為0.44(表3)。說明華北地區(qū)同一輪作類型內(nèi), 不同農(nóng)戶之間磷投入差異巨大, 農(nóng)戶對磷肥施用缺乏科學(xué)足夠認(rèn)識, 從另一側(cè)面反映了加強磷肥管理的必要性。

    表3 華北幾種典型輪作類型的周年磷平衡

    2.2 不同輪作類型的周年磷肥運籌(秋播與夏播磷肥投入比)情況

    調(diào)查的38個地塊中, 秋播與夏播磷肥投入比(周年磷肥運籌)>1和≤1的地塊數(shù)量正好相等(#DIV/0!代表夏播作物磷投入量為零, 此種情況按“P秋/P夏>1”計), 說明從周年磷肥統(tǒng)籌的角度來衡量(“P秋/P夏>1”為合理, “P秋/P夏≤1”為不合理), 所調(diào)查的50%地塊磷肥運籌不合理; 若按照輪作類型統(tǒng)計, 調(diào)查的5種輪作類型中, 1種(冬小麥-夏玉米)平均值>1, 為1.28(由于#DIV/0!無法計算, 此處統(tǒng)計平均值時去掉), 整體比較合理; 1種(冬小麥-花生)平均值為0.87, 整體表現(xiàn)為不合理; 另外3種(冬小麥-大豆、大蒜-棉花、冬小麥-棉花)均為夏播作物不施磷肥, 雖然按照此評價方法結(jié)果為合理, 但是由于夏播作物完全不施磷是否能夠保障該茬作物全生育期供應(yīng), 尚需更多田間試驗驗證。

    統(tǒng)計分析還表明, 同一輪作類型內(nèi)部, 不同農(nóng)戶之間的周年磷肥運籌差異巨大, 以樣本量最大的麥玉輪作類型為例(由于#DIV/0!無法計算, 此處統(tǒng)計平均值時去掉): 秋播與夏播磷肥投入比最大值為2.70, 最小值僅為0.90, 兩者相差2倍, 標(biāo)準(zhǔn)差為0.61。

    2.3 預(yù)警結(jié)果情況

    38個地塊中, 磷肥施用合理、茬次間分配合理、環(huán)境友好(深綠)的地塊僅4個, 占總樣本量10.53%; 磷肥周年投入總量合理, 但茬次間的分配不合理, 需調(diào)整磷肥的茬次間分配, 環(huán)境風(fēng)險較低(淺綠)的地塊5個, 占總樣本量13.16%; 磷肥周年投入偏少, 但茬次間的分配合理, 需增加磷肥的周年總投入量, 環(huán)境風(fēng)險較低(淺黃)的地塊5個, 占總樣本量13.16%; 磷肥周年投入偏少, 需調(diào)整季節(jié)分配, 增加磷肥的周年總投入量, 環(huán)境風(fēng)險中等(深黃)的地塊2個, 占總樣本量5.26%; 磷肥投入過量, 季節(jié)分配合理, 需減少磷肥的周年投入量, 環(huán)境風(fēng)險高(淺紅)的地塊9個, 占總樣本量23.68%; 磷肥投入過量, 且茬次間分配不合理, 需減少磷投入量, 且調(diào)整茬次間的分配, 環(huán)境風(fēng)險極高(深紅)的地塊13個, 占總樣本量34.21%。由以上結(jié)果可知, 磷流失高環(huán)境風(fēng)險(包括高和極高)的輪作地塊達(dá)57.89%, 即一半以上地塊的輪作方式存在高或極高的環(huán)境風(fēng)險, 亟需引起足夠的重視, 后續(xù)可通過調(diào)整輪作類型的方式加以改善, 綜合考慮經(jīng)濟效益和環(huán)境風(fēng)險, 尋找到最優(yōu)輪作方式。

    3 討論與結(jié)論

    3.1 討論

    關(guān)于磷平衡的計算對象, 本文采用了一個輪作周期中所有作物的磷素平衡, 而非單茬次作物的磷平衡, 這是由于對于有長期種植歷史的地塊來說, 一個輪作周期是一個基本的種植單元。而該地塊多年磷素的累積(或盈虧狀況)直接取決于輪作(而非單茬)磷平衡值的正負(fù)與絕對值大小。

    關(guān)于周年磷平衡的計算方法, 本文采用周年磷投入量與帶走量的比值(而非二者的差值)。即: 周年磷平衡=兩茬作物磷投入總量/兩茬作物磷帶走總量。通常磷平衡值的計算是投入與產(chǎn)出的差值, 即磷平衡的絕對值[2,27]。但此種方法得到的磷平衡值, 由于不同作物(輪作)類型之間對磷素吸收量的差異較大, 從而難以用統(tǒng)一的標(biāo)準(zhǔn)判斷結(jié)果的合理性, 而只能分別對不同輪作類型采用不同的數(shù)據(jù)單獨判斷, 限制了應(yīng)用的靈活性。本研究中磷平衡的計算采用投入與產(chǎn)出的比值(商)而非差值, 使不同輪作類型之間的磷平衡可以進行統(tǒng)一量化和比較; 當(dāng)然具體值還可根據(jù)將來試驗進一步優(yōu)化。

    為了操作方便, 本方法周年磷肥運籌簡化為以P秋/P夏=1(即兩茬作物的磷投入量相同)為界限, 僅分成>1和≤1兩種情況。將來如果某地有不同輪作類型茬次最佳比例數(shù)據(jù), 可依據(jù)數(shù)據(jù)對該地區(qū)該指標(biāo)進行細(xì)化更新。

    土壤質(zhì)地分類是指按照土壤顆粒組成的比例對土壤進行的分類。最常見的土壤質(zhì)地分類有國際制、美國農(nóng)業(yè)部制、蘇聯(lián)卡慶斯基制3種。這3種質(zhì)地分類制都是與其粒組分級標(biāo)準(zhǔn)和顆粒大小分析時的土粒分散方法相配套的[28]。但往往需要進行實驗室測定黏粒粒級數(shù)據(jù), 且分類級別都較多, 如1987年《中國土壤》(第2版)中公布了中國的質(zhì)地分類制, 分為3組12種質(zhì)地名稱[29]。本系統(tǒng)從實用性、簡便性角度將土壤質(zhì)地粗分為砂土、壤土和黏土3種, 未進一步細(xì)化。

    本文建立的預(yù)警方法是為單個地塊的輪作風(fēng)險給出相對合理的評價, 即應(yīng)用對象為地塊, 但是當(dāng)用戶數(shù)量達(dá)到一定規(guī)模后, 可進一步對區(qū)域性的所有輪作類型進行分類匯總, 找出該研究區(qū)域各輪作類型的共性風(fēng)險及差異性評價, 進而通過種植結(jié)構(gòu)調(diào)整等措施加以政策引導(dǎo), 形成區(qū)域農(nóng)業(yè)良性循環(huán)。因此, 該方法既可作為區(qū)域土壤養(yǎng)分管理預(yù)警依據(jù), 又可作為主動調(diào)控的方法。

    3.2 結(jié)論

    本文從影響輪作環(huán)境風(fēng)險的關(guān)鍵元素、關(guān)鍵環(huán)節(jié)入手, 提出適合我國小農(nóng)戶現(xiàn)實國情的、基于農(nóng)田地塊的輪作風(fēng)險預(yù)警指標(biāo)體系, 并以華北平原的農(nóng)戶調(diào)查數(shù)據(jù)為例具體說明了預(yù)警系統(tǒng)的使用方法。結(jié)果表明: 該預(yù)警方法科學(xué)、實用、簡便, 所需數(shù)據(jù)容易獲得; 3種色系的預(yù)警結(jié)果易于普通農(nóng)戶識別推廣。華北38個地塊預(yù)警結(jié)果表明, 磷肥施用合理、茬次間分配合理、環(huán)境友好(深綠)的地塊僅4個, 占總樣本量10.53%; 高環(huán)境風(fēng)險(包括高和極高)的輪作地塊達(dá)57.89%, 即一半以上的地塊目前的輪作方式存在高或極高的環(huán)境風(fēng)險。該預(yù)警方法既可用于農(nóng)戶對自己地塊磷流失風(fēng)險進行判斷, 進而選擇更優(yōu)化的輪作模式, 也可用于國家、區(qū)域或地方農(nóng)業(yè)主管部門進行種植業(yè)結(jié)構(gòu)調(diào)整時參考。

    [1] 李鵬, 關(guān)強, 鄔夢成, 等. 我國旱地兩熟區(qū)土壤養(yǎng)分時空變異的調(diào)查研究[J]. 土壤通報, 2017, 48(2): 380–386 LI P, GUAN Q, WU M C, et al. Investigation of temporal and spatial variability of soil nutrients in double cropping area of the upland, China[J]. Chinese Journal of Soil Science, 2017, 48(2): 380–386

    [2] 冀宏杰, 張懷志, 張維理, 等. 我國農(nóng)田磷養(yǎng)分平衡研究進展[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報, 2015, 23(1): 1–8 JI H J, ZHANG H Z, ZHANG W L, et al. Research progress on cropland phosphorus balance in China[J]. Chinese Journal of Eco-Agriculture, 2015, 23(1): 1–8

    [3] 宋科, 張維理, 徐愛國, 等. 太湖水網(wǎng)地區(qū)不同種植類型農(nóng)田磷素滲漏流失研究[J]. 植物營養(yǎng)與肥料學(xué)報, 2009, 15(6): 1288–1294 SONG K, ZHANG W L, XU A G, et al. Phosphorus leaching losses in different planting farmlands in the riverine plain area of Taihu Lake[J]. Plant Nutrition and Fertilizer Science, 2009, 15(6): 1288–1294

    [4] 張燕, 李永梅, 張懷志, 等. 滇池流域農(nóng)田徑流磷素流失的土壤影響因子[J]. 水土保持學(xué)報, 2011, 25(4): 41–45 ZHANG Y, LI Y M, ZHANG H Z, et al. The soil factors influencing phosphorus loss by field surface runoff in Dianchi watershed[J]. Journal of Soil and Water Conservation, 2011, 25(4): 41–45

    [5] 閆興成, 張重乾, 季銘, 等. 富營養(yǎng)化湖泊夏季表層水體溫室氣體濃度及其影響因素[J]. 湖泊科學(xué), 2018, 30(5): 1420–1428 YAN X C , ZHANG C Q, JI M, et al. Concentration of dissolved greenhouse gas and its influence factors in the summer surface water of eutrophic lake[J]. Journal of Lake Sciences, 2018, 30(5): 1420–1428

    [6] 鄒銀洪, 張潤宇, 陳敬安, 等. 黏土礦物在富營養(yǎng)化水體和底泥磷污染控制中的應(yīng)用研究進展[J]. 地球科學(xué)進展, 2018, 33(6): 578–589 ZOU Y H, ZHANG R Y, CHEN J A, et al. Research advance in the application of clay minerals to phosphorus pollution control in eutrophic water bodies and sediments[J]. Advances in Earth Science, 2018, 33(6): 578–589

    [7] 劉玉萍, 陳西, 王延華, 等. 滇池流域土壤養(yǎng)分分布及其對水體富營養(yǎng)化的影響[J]. 南京師大學(xué)報: 自然科學(xué)版, 2017, 40(4): 129–136 LIU Y P, CHEN X, WANG Y H, et al. Soil nutrients distribution in the small catchment and their effects on the eutrophication of Dianchi Lake[J]. Journal of Nanjing Normal University: Natural Science Edition, 2017, 40(4): 129–136

    [8] 陳思萱. 土壤重金屬污染時空模擬與環(huán)境風(fēng)險預(yù)警研究[D]. 長沙: 中南大學(xué), 2014 CHEN S X. Spatial-temporal modeling of heavy metal pollution in soil and its environmental risk warning[D]. Changsha: Central South University, 2014

    [9] 侯彥林. 中國農(nóng)田重金屬污染預(yù)警系統(tǒng)研究[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2012, 31(4): 697–705 HOU Y L. The warning system of heavy metal pollution in the farmland of China[J]. Journal of Agro-Environment Science, 2012, 31(4): 697–705

    [10] YAN B J, PAN Y C, WU W Y, et al. WebGIS-based early warning system of farmland environmental pollution caused by livestock waste[J]. Agricultural Science & Technology, 2011, 12(7): 1062–1066

    [11] 李維新, 張永春, 張海平, 等. 太湖流域水環(huán)境風(fēng)險預(yù)警系統(tǒng)構(gòu)建[J]. 生態(tài)與農(nóng)村環(huán)境學(xué)報, 2010, 26(z1): 4–8 LI W X, ZHANG Y C, ZHANG H P, et al. Construction of water environmental risk early warning system for Taihu Lake basin[J]. Journal of Ecology and Rural Environment, 2010, 26(z1): 4–8

    [12] 國家統(tǒng)計局. 第一次全國污染源普查公報[R/OL]. (2010-02-11). http://www.stats.gov.cn/tjsj/tjgb/qttjgb/qgqttjgb/201002/t20100211_30641.html National Bureau of Statistics. The First National Pollution Source Census Bulletin[R/OL]. (2010-02-11). http://www.stats. gov.cn/tjsj/tjgb/qttjgb/qgqttjgb/201002/t20100211_30641.html

    [13] 楊林章, 吳永紅. 農(nóng)業(yè)面源污染防控與水環(huán)境保護[J]. 中國科學(xué)院院刊, 2018, 33(2): 168–176 YANG L Z, WU Y H. Prevention and control of agricultural non-point source pollution and aquatic environmental protection[J]. Bulletin of Chinese Academy of Sciences, 2018, 33(2): 168–176

    [14] 中國人民共和國生態(tài)環(huán)境部. 2017中國生態(tài)環(huán)境狀況公報[R/OL]. (2018-05-31). http://www.mee.gov.cn/hjzl/zghjzkgb/ lnzghjzkgb/Ministry of Ecology and Environment of the People’s Republic of China. 2017 China Ecological Environment Bulletin[R/OL]. (2018-05-31). http://www.mee.gov.cn/hjzl/zghjzkgb/ lnzghjzkgb/

    [15] 中國人民共和國生態(tài)環(huán)境部. 2016中國生態(tài)環(huán)境狀況公報[R/OL]. (2017-06-05). http://www.mee.gov.cn/hjzl/zghjzkgb/ lnzghjzkgb/Ministry of Ecology and Environment of the People’s Republic of China. 2016 China Ecological Environment Bulletin[R/OL]. (2017-06-05). http://www.mee.gov.cn/hjzl/zghjzkgb/lnzghjzkgb/

    [16] 中國人民共和國生態(tài)環(huán)境部. 2015中國生態(tài)環(huán)境狀況公報[R/OL]. (2016-06-01). http://www.mee.gov.cn/hjzl/zghjzkgb/ lnzghjzkgb/Ministry of Ecology and Environment of the People’s Republic of China. 2015 China Ecological Environment Bulletin[R/OL]. (2016-06-01). http://www.mee.gov.cn/hjzl/zghjzkgb/lnzghjzkgb/

    [17] 吳啟華, 劉曉斌, 張淑香, 等. 施用常規(guī)磷水平的80%可實現(xiàn)玉米高產(chǎn)、磷素高效利用和土壤磷平衡[J]. 植物營養(yǎng)與肥料學(xué)報, 2016, 22(6): 1468–1476 WU Q H, LIU X B, ZHANG S X, et al. Application of 80% of routine phosphorus rate to keep high yield and P efficiency of maize and P balance in soil[J]. Journal of Plant Nutrition and Fertilizer, 2016, 22(6): 1468–1476

    [18] 劉彥伶, 李渝, 張雅蓉, 等. 長期施肥對黃壤性水稻土磷平衡及農(nóng)學(xué)閾值的影響[J]. 中國農(nóng)業(yè)科學(xué), 2016, 49(10): 1903–1912 LIU Y L, LI Y, ZHANG Y R, et al. Effect of long-term fertilization on the P balance and critical value of soil Olsen-P in paddy soil from yellow earth[J]. Scientia Agricultura Sinica, 2016, 49(10): 1903–1912

    [19] 展曉瑩, 任意, 張淑香, 等. 中國主要土壤有效磷演變及其與磷平衡的響應(yīng)關(guān)系[J]. 中國農(nóng)業(yè)科學(xué), 2015, 48(23): 4728–4737 ZHAN X Y, REN Y, ZHANG S X, et al. Changes in Olsen phosphorus concentration and its response to phosphorus balance in the main types of soil in China[J]. Scientia Agricultura Sinica, 2015, 48(23): 4728–4737

    [20] 李娟, 章明清, 孔慶波. 不同施肥模式對菜稻輪作產(chǎn)量和菜田氮磷平衡的影響[J]. 中國農(nóng)學(xué)通報, 2016, 32(3): 146–150 LI J, ZHANG M Q, KONG Q B. Effect of different fertilization patterns on yield of vegetable-rice rotation and nitrogen and phosphorus balance in vegetable field[J]. Chinese Agricultural Science Bulletin, 2016, 32(3): 146–150

    [21] COX F R, HENDRICKS S E. Soil test phosphorus and clay content effects on runoff water quality[J]. Journal of Environment Quality, 2000, 29(5): 1582

    [22] 劉玉民, 龍偉, 劉亞敏, 等. 不同種植模式下紫色土養(yǎng)分流失影響因子研究[J]. 水土保持學(xué)報, 2005, 19(5): 81–84 LIU Y M, LONG W, LIU Y M, et al. Study on influencing factors of nutrient loss under different land use patterns in purple soils[J]. Journal of Soil and Water Conservation, 2005, 19(5): 81–84

    [23] 魯如坤, 劉鴻翔, 聞大中, 等. 我國典型地區(qū)農(nóng)業(yè)生態(tài)系統(tǒng)養(yǎng)分循環(huán)和平衡研究Ⅳ. 農(nóng)田養(yǎng)分平衡的評價方法和原則[J]. 土壤通報, 1996, 27(5): 197–199 LU R K, LIU H X, WEN D Z, et al. Research on nutrient cycling and balance of agricultural ecosystem in typical areas of China. Ⅳ. Evaluation methods and principles of farmland nutrient balance[J]. Chinese Journal of Soil Science, 1996, 27(5): 197–199

    [24] 張甘霖, 李德成. 野外土壤描述與采樣手冊[M]. 北京: 科學(xué)出版社, 2017 ZHANG G L, LI D C. Manual of Field Soil Description and Sampling[M]. Beijing: Science Press, 2017

    [25] 全國農(nóng)業(yè)技術(shù)推廣服務(wù)中心, 中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與區(qū)劃研究所. 耕地質(zhì)量演變趨勢研究——國家級耕地土壤監(jiān)測數(shù)據(jù)整編[M]. 北京: 中國農(nóng)業(yè)科學(xué)技術(shù)出版社, 2008 National Agro-Tech Extension and Service Center, Institute of Agricultural Resources and Regional Planning of Chinese Academy of Agricultural Sciences. Research of the Evolution of Arable Land Quality — Date Set of National Soil Monitoring of Farmland[M]. Beijing: China Agricultural Science and Technology Press, 2008

    [26] 黃欣欣, 廖文華, 劉建玲, 等. 河北省典型土壤磷庫狀況及對長期施磷和秸稈還田的響應(yīng)[J]. 中國土壤與肥料, 2017, (6): 30–36HUANG X X, LIAO W H, LIU J L, et al. Response of long-term phosphorus fertilization and straw return to various fractions of phosphorus in typical soil of Hebei Province[J]. Soils and Fertilizers Sciences in China, 2017, (6): 30–36

    [27] 李書田, 金繼運. 中國不同區(qū)域農(nóng)田養(yǎng)分輸入、輸出與平衡[J]. 中國農(nóng)業(yè)科學(xué), 2011, 44(20): 4207–4229 LI S T, JIN J Y. Characteristics of nutrient input/output and nutrient balance in different regions of China[J]. Scientia Agricultura Sinica, 2011, 44(20): 4207–4229

    [28] 周健民, 沈仁芳. 土壤學(xué)大辭典[M]. 北京: 科學(xué)出版社, 2013 ZHOU J M, SHEN R F. Dictionary of Soil Science[M]. Beijing: Science Press, 2013

    [29] 熊毅, 李慶奎. 中國土壤[M]. 北京: 科學(xué)出版社, 1987 XIONG Y, LI Q K. Soils of China[M]. Beijing: Science Press, 1987

    Early warning system of field phosphorus loss risk for double cropping area in North China*

    JI Hongjie, ZHANG Huaizhi, ZHANG Renlian, XU Aiguo, TIAN Changyu

    (Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences / Key Laboratory of Plant Nutrient and Nutrient Cycle, Ministry of Agriculture and Rural Affairs, Beijing 100081, China)

    Rotation is a relatively fixed crop combination pattern that forms in various places over a long period. It is also the basic unit for the calculation of farmland nutrient balance. The accumulation of nutrient surplus or deficiency over the years directly affects the environment. Among the three essential nutrients for plants, the most important for the environment are phosphorus and nitrogen. Compared with nitrogen, the phosphorus cycle is relatively stable, has no gas loss, and reveals a high correlation between apparent nutrient balance and environmental pollution of phosphorus, and the calculation for which is convenient. Therefore, phosphorus was selected as the indicator for environment risk in this study. In order to establish a simple and effective method to detect early warnings of environmental effects, starting with the screening of key elements suitable for early warning and key links to phosphorus loss in farmland, this article proposed an early warning indicator system based on farmland plots suitable for small-holder farmers in China. The study also performed trial calculations with survey data from 38 plots, including five typical rotation systems in North China. The early warning system included three evaluation indicators — annual phosphorus balance, plough layer soil texture, and annual phosphate fertilizer operation. A total of six early warning levels of phosphorus loss was set in the early warning system — dark green, light green, light yellow, deep yellow, light red, and deep red, in which the red meant high phosphorus loss risk, the green meant low phosphorus loss risk, the yellow meant medium phosphorus loss risk, deep and light mean degree of the level. The characteristics of the early warning system were 1) measuring the entire rotation cycle rather than a single crop; and 2) the phosphorus balance considered the total balance of the anniversary (i.e., the phosphorus input and output balance of the entire rotation cycle) and considered the difference between different crops. Scientific allocation (i.e., phosphorus co-ordination) used relative equilibrium values rather than absolute equilibrium values so that the phosphorus balance between different rotation types could be uniformly quantified and compared. Soil texture was simplified to sand, loam, and clay species. The trial results showed that a low early warning level (dark green) which was rational in application of phosphate fertilizer with reasonable application rate and co-ordination between rotation crops only occupied 10.53% of all plots analyzed, while the high environmental risk (dark red, light red) rotation plots accounted for 57.89%, i.e., more than half of the plots using the current rotation methods had high or extremely high environmental risks. The early warning system was scientific, practical, and simple, and the required data were easy to be obtained. It can be used for farmers to judge the risk of phosphorus loss in their own plots and for national, regional, or local agricultural authorities to reference when adjusting planting structures.

    Double cropping area; Crop rotation system; Farmland nutrient; Phosphorus; Loss risk; Early warning system; North China

    , JI Hongjie, E-mail: jihongjie@caas.cn

    Apr. 19, 2019;

    Jun. 27, 2019

    S19

    2096-6237(2019)09-1394-08

    10.13930/j.cnki.cjea.190229

    2019-04-19

    2019-06-27

    冀宏杰, 主要從事土壤養(yǎng)分庫調(diào)控增效技術(shù)研究。E-mail: jihongjie@caas.cn

    * This study was supported by the Special Fund for Agro-scientific Research in the Public Interest of China (201503121-13) and the National Key R&D Program of China (2016YFC0501302).

    * 農(nóng)業(yè)農(nóng)村部公益性行業(yè)專項(201503121-13)和國家重點研發(fā)計劃項目(2016YFC0501302)資助

    冀宏杰, 張懷志, 張認(rèn)連, 徐愛國, 田昌玉. 華北兩熟區(qū)農(nóng)田磷流失風(fēng)險預(yù)警研究[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報(中英文), 2019, 27(9): 1394-1401

    JI H J, ZHANG H Z, ZHANG R L, XU A G, TIAN C Y.Early warning system of field phosphorus loss risk for double cropping area in North China[J]. Chinese Journal of Eco-Agriculture, 2019, 27(9): 1394-1401

    猜你喜歡
    投入量磷素輪作
    磷素添加對土壤水分一維垂直入滲特性的影響
    農(nóng)民瓜菜“接茬輪作”增效益
    詞匯加工深度理論對二語課堂詞匯附帶習(xí)得的影響
    山東青年(2020年10期)2020-11-23 02:01:12
    影響英語閱讀中詞匯附帶習(xí)得的學(xué)習(xí)者因素
    輪作休耕,讓土地歇一會兒
    耕地輪作效益調(diào)查報告
    投入量三要素對詞匯習(xí)得的影響探究
    文教資料(2014年23期)2014-02-23 07:06:16
    作物高效利用土壤磷素的研究進展
    土壤與作物(2013年4期)2013-03-11 18:47:35
    磷素營養(yǎng)對大豆磷素吸收及產(chǎn)量的影響
    磷素子粒生產(chǎn)效率不同的小麥品種磷素吸收利用差異
    欧美日本亚洲视频在线播放| 亚洲精品久久国产高清桃花| 久久久久性生活片| 级片在线观看| 欧美不卡视频在线免费观看| www日本黄色视频网| 乱人视频在线观看| 日韩 亚洲 欧美在线| 国产 一区 欧美 日韩| 婷婷亚洲欧美| 国产真实伦视频高清在线观看| eeuss影院久久| 日韩欧美一区二区三区在线观看| www.色视频.com| 我的女老师完整版在线观看| 色吧在线观看| 97在线视频观看| 日韩亚洲欧美综合| 可以在线观看毛片的网站| 亚洲精华国产精华液的使用体验 | 亚洲三级黄色毛片| 在线a可以看的网站| 无遮挡黄片免费观看| 有码 亚洲区| 男人和女人高潮做爰伦理| 久久人妻av系列| av福利片在线观看| 免费av毛片视频| 搡老妇女老女人老熟妇| 成熟少妇高潮喷水视频| 国产黄片美女视频| 日日摸夜夜添夜夜添小说| 国产欧美日韩精品一区二区| 悠悠久久av| 深夜精品福利| 给我免费播放毛片高清在线观看| 国产色婷婷99| 一区二区三区高清视频在线| 天堂动漫精品| 亚洲国产精品成人久久小说 | 国产精品一区二区性色av| 久久中文看片网| 亚洲高清免费不卡视频| 亚洲天堂国产精品一区在线| 国产极品精品免费视频能看的| 99久久中文字幕三级久久日本| 在线观看午夜福利视频| 日韩欧美国产在线观看| 晚上一个人看的免费电影| 免费黄网站久久成人精品| 特级一级黄色大片| 久久精品国产自在天天线| а√天堂www在线а√下载| 性欧美人与动物交配| 成人永久免费在线观看视频| 一a级毛片在线观看| 亚洲av熟女| 晚上一个人看的免费电影| 亚洲欧美日韩东京热| 久久精品人妻少妇| 五月伊人婷婷丁香| 黑人高潮一二区| 激情 狠狠 欧美| a级毛片a级免费在线| 欧美在线一区亚洲| 男人的好看免费观看在线视频| 高清毛片免费观看视频网站| 精品欧美国产一区二区三| 国产欧美日韩一区二区精品| 国产人妻一区二区三区在| 韩国av在线不卡| 欧美成人精品欧美一级黄| 亚洲国产精品国产精品| 99久久成人亚洲精品观看| 国产探花在线观看一区二区| 看片在线看免费视频| 老熟妇乱子伦视频在线观看| 国产高清不卡午夜福利| 亚洲国产欧洲综合997久久,| 午夜爱爱视频在线播放| 国产伦一二天堂av在线观看| 少妇丰满av| 99久久无色码亚洲精品果冻| 久久久久久久亚洲中文字幕| 亚洲激情五月婷婷啪啪| 一级毛片久久久久久久久女| 亚洲图色成人| 国产高潮美女av| 国产黄a三级三级三级人| av女优亚洲男人天堂| 日日摸夜夜添夜夜添小说| 日韩人妻高清精品专区| 国产黄色小视频在线观看| 老司机福利观看| 大香蕉久久网| 深夜精品福利| 麻豆av噜噜一区二区三区| 国产免费男女视频| avwww免费| 在线观看午夜福利视频| 啦啦啦韩国在线观看视频| 国产乱人视频| 偷拍熟女少妇极品色| 一级av片app| 免费无遮挡裸体视频| 精品国产三级普通话版| 高清毛片免费观看视频网站| 国产精品人妻久久久久久| 亚洲真实伦在线观看| 亚洲av一区综合| 久久久国产成人精品二区| 精品少妇黑人巨大在线播放 | 老熟妇乱子伦视频在线观看| 国产一区亚洲一区在线观看| 日本撒尿小便嘘嘘汇集6| 国产午夜精品论理片| 午夜影院日韩av| 最近中文字幕高清免费大全6| 亚洲久久久久久中文字幕| 亚洲欧美日韩东京热| 在线观看av片永久免费下载| 亚洲乱码一区二区免费版| 国产男人的电影天堂91| 在线看三级毛片| 久久鲁丝午夜福利片| 成人国产麻豆网| 一级毛片我不卡| 黄色欧美视频在线观看| 亚洲成人中文字幕在线播放| 高清毛片免费观看视频网站| 午夜精品一区二区三区免费看| 三级经典国产精品| 麻豆久久精品国产亚洲av| 国产 一区精品| 国产黄a三级三级三级人| 国产淫片久久久久久久久| 看片在线看免费视频| 人人妻人人澡欧美一区二区| 免费搜索国产男女视频| 夜夜夜夜夜久久久久| 女的被弄到高潮叫床怎么办| 成年版毛片免费区| 国产亚洲91精品色在线| 女的被弄到高潮叫床怎么办| 老熟妇仑乱视频hdxx| 18禁在线播放成人免费| 日本一本二区三区精品| 97人妻精品一区二区三区麻豆| 六月丁香七月| 亚洲精华国产精华液的使用体验 | 91久久精品国产一区二区成人| 久久久久精品国产欧美久久久| 国产精品一区二区性色av| 久久久精品欧美日韩精品| 中出人妻视频一区二区| av天堂中文字幕网| 国产一区二区三区在线臀色熟女| 99久久成人亚洲精品观看| 一区二区三区高清视频在线| av卡一久久| 亚洲三级黄色毛片| 亚洲国产精品sss在线观看| 黑人高潮一二区| 国产片特级美女逼逼视频| 日韩成人伦理影院| 国产精品一区二区三区四区久久| 午夜福利成人在线免费观看| 久久久久久九九精品二区国产| 床上黄色一级片| 色视频www国产| 国语自产精品视频在线第100页| 亚洲av免费在线观看| 麻豆乱淫一区二区| 日本成人三级电影网站| 成人漫画全彩无遮挡| 亚洲久久久久久中文字幕| 身体一侧抽搐| 岛国在线免费视频观看| 淫妇啪啪啪对白视频| 国产av在哪里看| 性插视频无遮挡在线免费观看| 免费人成在线观看视频色| 亚洲欧美日韩无卡精品| 国产精品一区www在线观看| 久久久久久久久久成人| 免费高清视频大片| 精品久久久久久成人av| 国产av在哪里看| 欧美xxxx黑人xx丫x性爽| 亚洲人成网站高清观看| 伦理电影大哥的女人| 欧美bdsm另类| 亚洲欧美清纯卡通| 国产探花极品一区二区| 国国产精品蜜臀av免费| 亚洲天堂国产精品一区在线| 91av网一区二区| 亚洲精品国产av成人精品 | 观看美女的网站| 亚洲国产精品久久男人天堂| 免费看a级黄色片| 日韩欧美精品免费久久| 一个人看视频在线观看www免费| 一进一出好大好爽视频| 97超级碰碰碰精品色视频在线观看| 国产黄片美女视频| 欧美日韩乱码在线| 日日摸夜夜添夜夜添小说| 在线观看一区二区三区| 国产日本99.免费观看| 黄色日韩在线| 一个人免费在线观看电影| 国产在线精品亚洲第一网站| 国产一区二区三区av在线 | 国产色婷婷99| 最好的美女福利视频网| 1024手机看黄色片| 欧美激情久久久久久爽电影| 舔av片在线| 一本精品99久久精品77| 一a级毛片在线观看| 精品人妻熟女av久视频| 亚洲第一区二区三区不卡| 欧美日韩乱码在线| 久久久欧美国产精品| 国产高清激情床上av| 精品人妻熟女av久视频| 日本黄色视频三级网站网址| 成人午夜高清在线视频| 色5月婷婷丁香| www.色视频.com| 99热这里只有是精品在线观看| 全区人妻精品视频| 69av精品久久久久久| 国产精品电影一区二区三区| 午夜a级毛片| 三级毛片av免费| 国产精品亚洲一级av第二区| 岛国在线免费视频观看| 麻豆一二三区av精品| 久久这里只有精品中国| 欧美成人精品欧美一级黄| 久久久久久九九精品二区国产| 夜夜爽天天搞| 亚洲国产日韩欧美精品在线观看| 有码 亚洲区| 直男gayav资源| 人人妻,人人澡人人爽秒播| 俺也久久电影网| 两性午夜刺激爽爽歪歪视频在线观看| 噜噜噜噜噜久久久久久91| 老师上课跳d突然被开到最大视频| 亚洲自拍偷在线| 国内精品宾馆在线| 天天一区二区日本电影三级| 国产一区二区三区在线臀色熟女| 国产欧美日韩精品一区二区| 亚洲国产欧洲综合997久久,| 老女人水多毛片| 欧美日韩精品成人综合77777| 国产精品一区www在线观看| 一级黄色大片毛片| 国产成人a∨麻豆精品| 看黄色毛片网站| 亚州av有码| 校园春色视频在线观看| 97人妻精品一区二区三区麻豆| 能在线免费观看的黄片| 1000部很黄的大片| 亚洲四区av| 国产精品一区二区三区四区免费观看 | 熟女人妻精品中文字幕| 精品一区二区三区视频在线观看免费| 日韩欧美精品免费久久| 婷婷精品国产亚洲av| 18禁在线播放成人免费| 99久久精品一区二区三区| 国产精品无大码| 久久久久国内视频| 超碰av人人做人人爽久久| 久久久久国内视频| 免费av观看视频| 亚洲综合色惰| 国产精品美女特级片免费视频播放器| 黄色配什么色好看| 精品乱码久久久久久99久播| 国产麻豆成人av免费视频| 免费看光身美女| 久久精品国产自在天天线| 国产精华一区二区三区| 真人做人爱边吃奶动态| 无遮挡黄片免费观看| 久久午夜亚洲精品久久| 搞女人的毛片| 亚洲av美国av| 人妻久久中文字幕网| 国产一级毛片七仙女欲春2| 欧美一区二区亚洲| 极品教师在线视频| 国产成人freesex在线 | 99久久精品国产国产毛片| 国产精品99久久久久久久久| 99热只有精品国产| 国产精品一二三区在线看| 一个人看的www免费观看视频| 国产av不卡久久| 99热这里只有是精品在线观看| 婷婷亚洲欧美| 九九热线精品视视频播放| 成人综合一区亚洲| 欧美精品国产亚洲| 国产探花在线观看一区二区| 国产国拍精品亚洲av在线观看| 国产大屁股一区二区在线视频| 久久亚洲精品不卡| 国产精品伦人一区二区| 村上凉子中文字幕在线| 国产精品一二三区在线看| 91午夜精品亚洲一区二区三区| 男人和女人高潮做爰伦理| 亚洲人成网站在线播放欧美日韩| 欧美在线一区亚洲| 免费av观看视频| 床上黄色一级片| 日本一本二区三区精品| 国产精品一区二区三区四区久久| 成人综合一区亚洲| 国产黄a三级三级三级人| 天堂网av新在线| 久久精品国产清高在天天线| 乱码一卡2卡4卡精品| 亚洲美女视频黄频| 尤物成人国产欧美一区二区三区| www.色视频.com| 亚洲美女搞黄在线观看 | 午夜久久久久精精品| 婷婷六月久久综合丁香| 久久精品国产亚洲av天美| 色噜噜av男人的天堂激情| 成人鲁丝片一二三区免费| 免费观看在线日韩| 最近视频中文字幕2019在线8| 欧美国产日韩亚洲一区| 亚洲精品在线观看二区| 久久人人爽人人片av| 免费在线观看影片大全网站| 国产亚洲av嫩草精品影院| 中文亚洲av片在线观看爽| 成人无遮挡网站| 黄色日韩在线| 高清毛片免费观看视频网站| 欧美日韩综合久久久久久| 99久久无色码亚洲精品果冻| 网址你懂的国产日韩在线| 欧美激情在线99| 又黄又爽又免费观看的视频| 性插视频无遮挡在线免费观看| 99久久九九国产精品国产免费| 女人十人毛片免费观看3o分钟| 国产成人freesex在线 | 身体一侧抽搐| 国产精品一区www在线观看| 1000部很黄的大片| 热99re8久久精品国产| ponron亚洲| 亚洲成人久久爱视频| 精品人妻偷拍中文字幕| 欧美色欧美亚洲另类二区| 97碰自拍视频| 最新在线观看一区二区三区| 午夜福利18| 久久久久性生活片| 91午夜精品亚洲一区二区三区| 亚洲国产色片| 国产乱人视频| 非洲黑人性xxxx精品又粗又长| 久久热精品热| 日韩 亚洲 欧美在线| 麻豆久久精品国产亚洲av| 给我免费播放毛片高清在线观看| 亚洲国产高清在线一区二区三| 久久久久国内视频| 成人特级黄色片久久久久久久| 中国美女看黄片| 国产黄片美女视频| 在线观看午夜福利视频| 亚洲电影在线观看av| 久久久久九九精品影院| 亚洲国产高清在线一区二区三| 午夜老司机福利剧场| 黄色日韩在线| 午夜免费激情av| 亚洲第一区二区三区不卡| 男人的好看免费观看在线视频| 99视频精品全部免费 在线| 18禁在线无遮挡免费观看视频 | 午夜免费激情av| 高清午夜精品一区二区三区 | 啦啦啦啦在线视频资源| 日本色播在线视频| 特大巨黑吊av在线直播| 少妇被粗大猛烈的视频| 日本三级黄在线观看| 少妇丰满av| 亚洲av不卡在线观看| 看非洲黑人一级黄片| 色av中文字幕| 搡女人真爽免费视频火全软件 | 日韩欧美 国产精品| 国产淫片久久久久久久久| 色哟哟哟哟哟哟| 麻豆国产av国片精品| 免费在线观看影片大全网站| 可以在线观看的亚洲视频| 国产高潮美女av| 别揉我奶头~嗯~啊~动态视频| 我的老师免费观看完整版| 99久久精品国产国产毛片| 你懂的网址亚洲精品在线观看 | 精品久久国产蜜桃| 午夜日韩欧美国产| 天美传媒精品一区二区| 国产白丝娇喘喷水9色精品| 99久久无色码亚洲精品果冻| 国产高清不卡午夜福利| 两个人视频免费观看高清| 国产在线男女| 亚洲一级一片aⅴ在线观看| 日本在线视频免费播放| 无遮挡黄片免费观看| 成人性生交大片免费视频hd| 插逼视频在线观看| av在线播放精品| 小蜜桃在线观看免费完整版高清| 久久人人精品亚洲av| 能在线免费观看的黄片| 在线天堂最新版资源| 女人十人毛片免费观看3o分钟| 91在线观看av| 精品久久久久久久久av| 一本精品99久久精品77| 亚洲成人精品中文字幕电影| 午夜爱爱视频在线播放| 插逼视频在线观看| 久久6这里有精品| 日韩成人伦理影院| 亚洲内射少妇av| 精品久久久久久久久久久久久| 精品午夜福利视频在线观看一区| 亚洲人成网站在线播放欧美日韩| 欧美高清成人免费视频www| 日韩强制内射视频| 日韩精品中文字幕看吧| 午夜福利视频1000在线观看| 日韩精品青青久久久久久| 亚洲精品456在线播放app| 九九久久精品国产亚洲av麻豆| 成人美女网站在线观看视频| 舔av片在线| 色综合色国产| 99在线视频只有这里精品首页| 日本五十路高清| 国产91av在线免费观看| 波多野结衣高清作品| 91在线观看av| 六月丁香七月| 一区二区三区高清视频在线| 免费看日本二区| 香蕉av资源在线| 国产精品电影一区二区三区| 天美传媒精品一区二区| 全区人妻精品视频| 久久久久性生活片| 久久婷婷人人爽人人干人人爱| 国内精品久久久久精免费| 亚洲性夜色夜夜综合| 看十八女毛片水多多多| 亚洲欧美中文字幕日韩二区| av国产免费在线观看| 一本精品99久久精品77| 日韩强制内射视频| h日本视频在线播放| 国产精品嫩草影院av在线观看| 蜜臀久久99精品久久宅男| 99热这里只有精品一区| 亚洲内射少妇av| 久久精品国产亚洲av涩爱 | 亚洲国产精品国产精品| 久久久久久久午夜电影| 老熟妇乱子伦视频在线观看| 亚洲国产色片| eeuss影院久久| 成人美女网站在线观看视频| 亚洲av一区综合| 插逼视频在线观看| 少妇熟女欧美另类| 亚洲精品日韩av片在线观看| 国产高清三级在线| 欧美一区二区亚洲| 全区人妻精品视频| 韩国av在线不卡| 久久婷婷人人爽人人干人人爱| 成人一区二区视频在线观看| 综合色丁香网| 亚洲欧美精品自产自拍| 国产精品伦人一区二区| 日本精品一区二区三区蜜桃| 特级一级黄色大片| 亚洲av免费高清在线观看| 国产精品久久久久久亚洲av鲁大| 国产久久久一区二区三区| 色哟哟·www| 国产精品国产高清国产av| 国产精品嫩草影院av在线观看| 色在线成人网| 99热这里只有是精品50| 校园人妻丝袜中文字幕| 日本黄大片高清| 干丝袜人妻中文字幕| 精品乱码久久久久久99久播| 久久久久免费精品人妻一区二区| 亚洲精品国产av成人精品 | 国产在视频线在精品| 97在线视频观看| 网址你懂的国产日韩在线| 人妻少妇偷人精品九色| 成人国产麻豆网| 国产伦精品一区二区三区四那| 欧美在线一区亚洲| 中文字幕av成人在线电影| 亚洲美女搞黄在线观看 | 国产精品一区二区三区四区免费观看 | 国产午夜精品论理片| 插阴视频在线观看视频| 久久人人精品亚洲av| 国产高清三级在线| 搡老妇女老女人老熟妇| 九九在线视频观看精品| 亚洲最大成人手机在线| 亚洲欧美精品自产自拍| 国产人妻一区二区三区在| 波多野结衣高清无吗| 国产69精品久久久久777片| 成人三级黄色视频| 国产久久久一区二区三区| 国产精品免费一区二区三区在线| 99热这里只有是精品50| 麻豆国产97在线/欧美| 看免费成人av毛片| 国产片特级美女逼逼视频| 最新中文字幕久久久久| 中文字幕av在线有码专区| av卡一久久| 日日啪夜夜撸| 亚洲天堂国产精品一区在线| 在线观看美女被高潮喷水网站| 日本-黄色视频高清免费观看| 女人被狂操c到高潮| 免费高清视频大片| 三级经典国产精品| 久久久久性生活片| 亚洲精品色激情综合| 少妇熟女aⅴ在线视频| 好男人在线观看高清免费视频| 国产一区二区在线观看日韩| 国产在视频线在精品| 天堂动漫精品| 看片在线看免费视频| 日韩成人av中文字幕在线观看 | 两个人视频免费观看高清| av福利片在线观看| 午夜福利在线观看吧| 国产精品一区二区性色av| aaaaa片日本免费| 国产不卡一卡二| 少妇的逼水好多| 一区二区三区免费毛片| 2021天堂中文幕一二区在线观| 十八禁国产超污无遮挡网站| av视频在线观看入口| 国产爱豆传媒在线观看| 久久久久精品国产欧美久久久| 一级毛片久久久久久久久女| 亚洲成人久久性| 欧美国产日韩亚洲一区| 超碰av人人做人人爽久久| 12—13女人毛片做爰片一| 别揉我奶头 嗯啊视频| 俄罗斯特黄特色一大片| 国产免费一级a男人的天堂| 校园春色视频在线观看| 激情 狠狠 欧美| 国产一级毛片七仙女欲春2| 99热这里只有精品一区| 99九九线精品视频在线观看视频| 免费观看在线日韩| 精品99又大又爽又粗少妇毛片| 丰满乱子伦码专区| 精品熟女少妇av免费看| 免费黄网站久久成人精品| 成人漫画全彩无遮挡| 亚洲成人久久性| 精品人妻一区二区三区麻豆 | 日韩欧美在线乱码| 午夜免费激情av| 少妇被粗大猛烈的视频| 国产亚洲精品综合一区在线观看| 久久国内精品自在自线图片| 成人国产麻豆网| 欧美日韩国产亚洲二区| 午夜爱爱视频在线播放| 国产精品免费一区二区三区在线| 又黄又爽又免费观看的视频| 在线看三级毛片| 国产国拍精品亚洲av在线观看|