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

    中國(guó)化肥面源污染環(huán)境風(fēng)險(xiǎn)時(shí)空變化

    2017-08-11 08:46:42劉欽普
    關(guān)鍵詞:低度面源化肥

    劉欽普

    (南京曉莊學(xué)院環(huán)境科學(xué)學(xué)院,南京 211171)

    中國(guó)化肥面源污染環(huán)境風(fēng)險(xiǎn)時(shí)空變化

    劉欽普

    (南京曉莊學(xué)院環(huán)境科學(xué)學(xué)院,南京 211171)

    為了合理施肥,防治化肥面源污染,對(duì)中國(guó)化肥施用的環(huán)境風(fēng)險(xiǎn)進(jìn)行研究。通過(guò)建立化肥施用環(huán)境安全閾值模型和環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)方法,計(jì)算中國(guó)各省區(qū)時(shí)空維度化肥施用環(huán)境安全閾值,評(píng)價(jià)化肥施用環(huán)境風(fēng)險(xiǎn)程度及變化。研究結(jié)果表明.中國(guó)化肥施用環(huán)境安全閾值近25年來(lái)在200~300 kg·hm-2之間變化,平均值為243 kg·hm-2,近15年來(lái)的平均值為251 kg·hm-2,與國(guó)家環(huán)保部生態(tài)鄉(xiāng)鎮(zhèn)建設(shè)規(guī)定的化肥施用強(qiáng)度標(biāo)準(zhǔn)基本吻合;總化肥施用在1994年由環(huán)境安全轉(zhuǎn)為低度風(fēng)險(xiǎn),氮、磷、鉀單質(zhì)化肥施用分別在1988、1999、2008年由環(huán)境安全轉(zhuǎn)為低度風(fēng)險(xiǎn)狀態(tài);2014年中國(guó)總化肥施用環(huán)境安全閾值平均為285 kg·hm-2,化肥施用環(huán)境風(fēng)險(xiǎn)指數(shù)平均為0.54,氮、磷、鉀化肥風(fēng)險(xiǎn)指數(shù)平均分別為0.54、0.53、0.49,化肥施用總體處于低度環(huán)境風(fēng)險(xiǎn);2014年西藏等11個(gè)省區(qū)化肥施用處于環(huán)境安全狀態(tài),遼寧等16個(gè)省區(qū)處于低度環(huán)境風(fēng)險(xiǎn)狀態(tài),河南、海南和天津三省區(qū)處于中度環(huán)境風(fēng)險(xiǎn),陜西處于嚴(yán)重環(huán)境風(fēng)險(xiǎn)??傊?中國(guó)化肥施用環(huán)境風(fēng)險(xiǎn)自1994年進(jìn)入低度風(fēng)險(xiǎn)狀態(tài),并呈現(xiàn)上升的趨勢(shì),自2006年以后有所降低。2014年全國(guó)各省區(qū)的化肥施用環(huán)境風(fēng)險(xiǎn)程度總體為低度環(huán)境風(fēng)險(xiǎn),但省區(qū)之間差異較大,各類(lèi)風(fēng)險(xiǎn)呈現(xiàn)出聚集分布的特點(diǎn)。

    化肥施用強(qiáng)度;環(huán)境安全閾值;環(huán)境風(fēng)險(xiǎn)指數(shù);時(shí)空變化;中國(guó)

    隨著點(diǎn)源污染環(huán)境風(fēng)險(xiǎn)的不斷控制,面源污染環(huán)境風(fēng)險(xiǎn)已成為世界范圍內(nèi)地表水與地下水污染的主要風(fēng)險(xiǎn)來(lái)源[1]。為了避免農(nóng)業(yè)面源污染現(xiàn)場(chǎng)監(jiān)測(cè)耗時(shí)費(fèi)力的困難,人們?cè)噲D用模型方法對(duì)面源污染進(jìn)行模擬和環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)。當(dāng)前研究的熱點(diǎn)主要集中在1980年代后國(guó)外出現(xiàn)的SWAT、AnnAGNPS、APEX等大型的面源污染機(jī)理模型的改進(jìn)和應(yīng)用[2-3]。汪朝輝等[4]基于GIS和SWAT模型對(duì)清江流域面源污染進(jìn)行模擬,揭示了面源污染環(huán)境風(fēng)險(xiǎn)的主要影響因素是降水、土地利用方式和耕作制度等。田耀武等[5]基于AnnAGNPS模型對(duì)三峽庫(kù)區(qū)秭歸縣面源污染輸出進(jìn)行評(píng)價(jià),認(rèn)為該模型對(duì)秭歸縣域泥沙和營(yíng)養(yǎng)物質(zhì)等模擬方式有待進(jìn)一步改進(jìn),才能更好地成為面源污染環(huán)境風(fēng)險(xiǎn)防治管理和調(diào)控工具。尹黎明等[6]認(rèn)為APEX模型在淮河中上游坡面尺度進(jìn)行不同土地管理措施下的水土流失評(píng)價(jià)較為適宜。國(guó)內(nèi)學(xué)者也構(gòu)建了相關(guān)面源污染模型。劉曼蓉等、李懷恩等、李定強(qiáng)等[7]分別建立了暴雨徑流污染的概化模型與統(tǒng)計(jì)相關(guān)模型、降雨量-徑流量-污染負(fù)荷之間的數(shù)學(xué)統(tǒng)計(jì)模型和面源污染物遷移機(jī)理模型等。陳會(huì)等[8]基于排水過(guò)程分析,利用模擬方法,對(duì)水稻灌區(qū)中地表排水和稻田滲漏排水的面源污染過(guò)程進(jìn)行研究。以上的面源污染模擬和環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)模型體現(xiàn)了對(duì)面源污染的綜合研究,模擬過(guò)程也較為復(fù)雜,但對(duì)化肥施用的單項(xiàng)污染建立的環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)模型還不多。劉欽普[9]提出了化肥污染環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)經(jīng)驗(yàn)?zāi)P?對(duì)中國(guó)的化肥施用環(huán)境風(fēng)險(xiǎn)進(jìn)行了初步研究。但是,該模型的環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)標(biāo)準(zhǔn)主要是采用國(guó)家環(huán)保部門(mén)制定的生態(tài)縣和生態(tài)鄉(xiāng)鎮(zhèn)化肥施用強(qiáng)度不低于250 kg·hm-2的常量標(biāo)準(zhǔn)[10-11],沒(méi)有考慮評(píng)價(jià)標(biāo)準(zhǔn)的時(shí)空變化。本文利用筆者建立的化肥施用環(huán)境安全閾值計(jì)算模型,對(duì)中國(guó)化肥施用的環(huán)境風(fēng)險(xiǎn)時(shí)空變化做深入研究,期望完善化肥施用環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)方法,并對(duì)我國(guó)生態(tài)農(nóng)業(yè)建設(shè)的化肥面源污染控制提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 化肥污染環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)方法

    化肥施用面源污染環(huán)境風(fēng)險(xiǎn)是農(nóng)業(yè)生產(chǎn)活動(dòng)中施用化肥對(duì)生態(tài)環(huán)境造成污染的可能性,是一種非突發(fā)性環(huán)境風(fēng)險(xiǎn)。這種風(fēng)險(xiǎn)是客觀存在的,具有很大的不確定性。筆者設(shè)計(jì)的環(huán)境風(fēng)險(xiǎn)指數(shù)模型得出的數(shù)值是介于0和1之間的小數(shù),有風(fēng)險(xiǎn)程度的上限和下限,便于對(duì)風(fēng)險(xiǎn)程度的認(rèn)識(shí)和比較[9]。其計(jì)算公式為.

    式中:Rt為化肥污染環(huán)境風(fēng)險(xiǎn)總指數(shù);Ri為單質(zhì)肥料(氮、磷或鉀等)污染環(huán)境風(fēng)險(xiǎn)指數(shù);Wi為單質(zhì)肥料氮、磷、鉀環(huán)境污染風(fēng)險(xiǎn)權(quán)重,這里依據(jù)前期研究,分別取0.648、0.230、0. 122;Ti為單質(zhì)肥料環(huán)境安全閾值,是指在為獲得某一季作物目標(biāo)產(chǎn)量而不危害環(huán)境的某種化肥的最大施用量(即單位播種面積化肥施用量上限),計(jì)算方法見(jiàn)下列公式(3)、(4)、(5);Fi為某化肥施用強(qiáng)度,是指本年內(nèi)實(shí)際用于農(nóng)業(yè)生產(chǎn)的單位播種面積化肥施用量,kg·hm-2。

    由公式(1)、(2)可見(jiàn),Ri(或Rt)介于0和1之間,當(dāng)Ri等于0.5時(shí),則Fi和Ti兩者相等,是施肥環(huán)境安全的臨界點(diǎn)。通過(guò)比較化肥施用強(qiáng)度與環(huán)境安全閾值的大小,即可表征化肥施用的環(huán)境風(fēng)險(xiǎn)大小。

    1.2 化肥施用環(huán)境安全閾值測(cè)度方法

    化肥施用環(huán)境風(fēng)險(xiǎn)指數(shù)模型公式(2)中的環(huán)境安全閾值是一個(gè)關(guān)鍵參數(shù)。假設(shè)化肥施用環(huán)境安全閾值與作物目標(biāo)產(chǎn)量有關(guān),作物目標(biāo)產(chǎn)量綜合反映自然生態(tài)條件、社會(huì)經(jīng)濟(jì)條件和農(nóng)田管理水平等因素的綜合作用。作物目標(biāo)產(chǎn)量可以根據(jù)某地區(qū)近幾年的產(chǎn)量的平均值向上浮動(dòng)10%得到。由于糧食作物是各地主要的農(nóng)作物,其產(chǎn)量高低可以代表該地區(qū)自然經(jīng)濟(jì)狀況。一般來(lái)講,化肥施用環(huán)境安全閾值應(yīng)小于當(dāng)?shù)刈魑镒罡弋a(chǎn)量時(shí)的化肥施用量。例如,崔玉亭等的研究表明,在某一產(chǎn)量要求下,考慮生態(tài)為主的化肥施用生態(tài)適宜量是考慮經(jīng)濟(jì)為主的化肥施用經(jīng)濟(jì)適宜量的85%左右[12]。參考巨曉棠提出的理論施氮量的計(jì)算方法[13],依據(jù)糧食作物目標(biāo)產(chǎn)量確定施氮量,在此施氮量的基礎(chǔ)上,向下浮動(dòng)10%~20%,得出氮肥施用環(huán)境安全閾值。根據(jù)大田作物氮磷鉀施肥比例一般為1∶0.5∶0.5的施肥要求,筆者提出了氮、磷、鉀總肥及單質(zhì)化肥施用的環(huán)境安全閾值測(cè)度模型.

    式中:Tt為總化肥施用環(huán)境安全閾值;ρ為化肥施用環(huán)境安全閾值調(diào)節(jié)系數(shù),一般取0.8或0. 9;A為單位產(chǎn)量的作物需氮量;Yi為某地區(qū)近期n年中某一年的作物產(chǎn)量;TN為氮肥環(huán)境安全閾值;TP為磷肥環(huán)境安全閾值;TK為鉀肥環(huán)境安全閾值。

    1.3 化肥施用環(huán)境風(fēng)險(xiǎn)分級(jí)方法

    根據(jù)化肥施用強(qiáng)度(Fi)超過(guò)環(huán)境安全閾值(Ti)的倍數(shù),把化肥施用環(huán)境風(fēng)險(xiǎn)程度從環(huán)境安全到極嚴(yán)重風(fēng)險(xiǎn)分為5個(gè)不同的風(fēng)險(xiǎn)等級(jí)類(lèi)型,如表1所示。將表1分類(lèi)依據(jù)中的有關(guān)數(shù)據(jù),代入方程(2)進(jìn)行計(jì)算,得出相應(yīng)的風(fēng)險(xiǎn)指數(shù)分級(jí)。

    表1 化肥施用環(huán)境風(fēng)險(xiǎn)指數(shù)(Rior Rt)分級(jí)及環(huán)境風(fēng)險(xiǎn)分類(lèi)Table 1 Classification of fertilization environmental risk index(Rior Rt)

    本研究使用的基礎(chǔ)資料主要來(lái)自國(guó)家統(tǒng)計(jì)局網(wǎng)站統(tǒng)計(jì)資料和有關(guān)年份中國(guó)統(tǒng)計(jì)年鑒[14],對(duì)無(wú)法直接通過(guò)年鑒或調(diào)查獲得的部分?jǐn)?shù)據(jù)用相關(guān)數(shù)據(jù)替代或估算。由于各統(tǒng)計(jì)年鑒中都沒(méi)有給出復(fù)合肥中氮磷鉀的含量,根據(jù)市場(chǎng)調(diào)查情況,參考有關(guān)文獻(xiàn),復(fù)合肥中的氮磷鉀含量統(tǒng)一按1∶1∶1處理[15]。

    2 結(jié)果與分析

    2.1 中國(guó)化肥施用環(huán)境風(fēng)險(xiǎn)時(shí)間變化特征

    根據(jù)作物播種面積(含果園種植面積)和化肥施用總量計(jì)算1980—2014年中國(guó)各年份單位播種面積化肥施用總強(qiáng)度。按照三年滑動(dòng)平均的方法,計(jì)算各年份糧食作物的目標(biāo)產(chǎn)量,用公式(1)~(5)計(jì)算出各年份的環(huán)境風(fēng)險(xiǎn)指數(shù)。35年來(lái),我國(guó)化肥施用總強(qiáng)度及其環(huán)境風(fēng)險(xiǎn)指數(shù)總體上呈現(xiàn)增長(zhǎng)的趨勢(shì)(表2)。

    計(jì)算結(jié)果表明,中國(guó)化肥施用總強(qiáng)度及氮、磷、鉀單質(zhì)化肥施用強(qiáng)度分別從1980年的86.4、64.2、19.2、3.0 kg·hm-2增加到2014年的333.7、174.2、87.2、75.8 kg·hm-2,35年來(lái)分別增長(zhǎng)到3.9、2.7、4.5、25.5倍。從2005年我國(guó)普遍實(shí)行配方施肥以來(lái),氮、磷化肥施用強(qiáng)度增長(zhǎng)緩慢,鉀肥施用強(qiáng)度增加較快(圖1)。N∶P2O5∶K2O(簡(jiǎn)寫(xiě)為N∶P∶K)的比例1980年為1∶0.3∶0.05,2005年為1∶0.44∶0.35,2014年為1∶0.50∶0.43,逐漸趨于1∶0.5∶0.5的發(fā)達(dá)國(guó)家水平[16]。同時(shí),我國(guó)化肥施用環(huán)境安全閾值近25年來(lái)在200~300 kg·hm-2之間變化,平均值是243 kg·hm-2,近15年來(lái)的平均值是251 kg· hm-2,特別是在2000—2010年間,化肥施用環(huán)境安全閾值平均值為250 kg·hm-2,正好與我國(guó)環(huán)境保護(hù)部門(mén)在2010年制定的生態(tài)鄉(xiāng)鎮(zhèn)建設(shè)化肥施用強(qiáng)度小于250 kg·hm-2的評(píng)價(jià)指標(biāo)相吻合??偦适┯铆h(huán)境風(fēng)險(xiǎn)指數(shù)及氮磷鉀單質(zhì)肥料的風(fēng)險(xiǎn)指數(shù)分別從1980年的0.39、0.47、0.35、0.08增加到2014年的0.54、0.55、0.55、0.51,總體處于低度環(huán)境風(fēng)險(xiǎn)(圖2)。由圖2不難發(fā)現(xiàn).①氮、磷、鉀化肥施用風(fēng)險(xiǎn)分別在1988、1999、2008年由環(huán)境安全轉(zhuǎn)為低度風(fēng)險(xiǎn),由于三個(gè)單質(zhì)化肥的平衡作用,總化肥風(fēng)險(xiǎn)自1994年由環(huán)境安全轉(zhuǎn)為低度風(fēng)險(xiǎn)。②自2006年以來(lái),氮肥的風(fēng)險(xiǎn)有所下降,磷肥的風(fēng)險(xiǎn)保持穩(wěn)定,鉀肥風(fēng)險(xiǎn)卻由安全增加為低度風(fēng)險(xiǎn)。

    總的來(lái)說(shuō),我國(guó)化肥施用環(huán)境風(fēng)險(xiǎn)自2006年以后呈降低趨勢(shì)。其主要原因是,2005年以后,我國(guó)大部分地區(qū)實(shí)行測(cè)土配方施肥,注重了氮、磷、鉀養(yǎng)分的合理搭配,雖然總的化肥施用量有所增加,但糧食產(chǎn)量增加,提高了化肥施用環(huán)境安全閾值,加上氮肥的施用風(fēng)險(xiǎn)有所下降、磷肥的施用風(fēng)險(xiǎn)保持穩(wěn)定、鉀肥風(fēng)險(xiǎn)稍微上升,使得總的化肥施用環(huán)境風(fēng)險(xiǎn)有所降低。

    表2 1980—2014年中國(guó)化肥施用總強(qiáng)度(Ft/kg·hm-2)及其環(huán)境風(fēng)險(xiǎn)指數(shù)(Rt)的變化Table 2 Total fertilization intensity(Ft/kg·hm-2)and fertilization risk(Rt)from 1980 to 2014 in China

    圖1 1980—2014年中國(guó)化肥施用強(qiáng)度的變化Figure 1 Changes of fertilization intensity in China from 1980 to 2014

    圖2 1980—2014年中國(guó)化肥施用環(huán)境風(fēng)險(xiǎn)變化Figure 2 Environmental risk changes of fertilization in China from 1980 to 2014

    2.2 中國(guó)化肥施用環(huán)境風(fēng)險(xiǎn)空間變化特征

    利用國(guó)家統(tǒng)計(jì)局編制的《2015中國(guó)統(tǒng)計(jì)年鑒》[17]數(shù)據(jù),根據(jù)公式(3)、(4)、(5),分別計(jì)算出2014年中國(guó)各省區(qū)的氮磷鉀單質(zhì)肥料和總化肥施用環(huán)境安全閾值,然后根據(jù)公式(1)、(2)分別計(jì)算出各省的化肥施用環(huán)境風(fēng)險(xiǎn)指數(shù)。結(jié)果表明,2014年中國(guó)總化肥施用環(huán)境安全閾值平均為285 kg·hm-2,環(huán)境風(fēng)險(xiǎn)指數(shù)為0.54,為低度環(huán)境風(fēng)險(xiǎn);氮磷鉀化肥施用環(huán)境風(fēng)險(xiǎn)指數(shù)分別為0.54、0.53、0.49。各省區(qū)的化肥施用環(huán)境風(fēng)險(xiǎn)差異較大,呈現(xiàn)出聚集分布的特點(diǎn),如表3和圖3所示。

    圖3a表明,2014年西藏等11個(gè)省區(qū)處于化肥施用環(huán)境安全狀態(tài),除黑龍江、吉林和上海外,多分布在西部和西南部地區(qū);遼寧等16個(gè)省區(qū)處于低度環(huán)境風(fēng)險(xiǎn)狀態(tài),除新疆外,多分布在華北地區(qū)和華南沿海地區(qū);河南、海南和天津三省市處于中度環(huán)境風(fēng)險(xiǎn),陜西處于嚴(yán)重環(huán)境風(fēng)險(xiǎn)。對(duì)于單質(zhì)化肥氮肥來(lái)說(shuō),其環(huán)境風(fēng)險(xiǎn)分布除個(gè)別省外基本與總肥分布一致,說(shuō)明當(dāng)前中國(guó)的化肥施用環(huán)境風(fēng)險(xiǎn)主要還是氮肥施用風(fēng)險(xiǎn)(圖3b)。磷肥施用環(huán)境風(fēng)險(xiǎn)相對(duì)較平衡,除東北三省和滬、浙、閩、湘、貴、渝、川、藏這11個(gè)省區(qū)為環(huán)境安全區(qū)域和晉陜?cè)ソ颦?個(gè)省區(qū)為中度風(fēng)險(xiǎn)外,其余15個(gè)省區(qū)屬于低度環(huán)境風(fēng)險(xiǎn)(圖3c)。鉀肥施用環(huán)境風(fēng)險(xiǎn)總體最低,但差別較大,從環(huán)境安全到極嚴(yán)重風(fēng)險(xiǎn)都有,海南省為極嚴(yán)重風(fēng)險(xiǎn)區(qū),廣西為嚴(yán)重風(fēng)險(xiǎn)區(qū),廣東、福建和陜西為中風(fēng)險(xiǎn)區(qū),云南、安徽、河南、山西、山東、吉林、北京、天津?yàn)榈投蕊L(fēng)險(xiǎn)區(qū),其余18個(gè)省區(qū)為環(huán)境安全區(qū)域(圖3d)。

    表3 2014年中國(guó)化肥施用環(huán)境風(fēng)險(xiǎn)指數(shù)及風(fēng)險(xiǎn)類(lèi)型Table 3 Environmental risk indexes and types of fertilization environmental risk in China in 2014

    圖3 2014年中國(guó)化肥施用環(huán)境風(fēng)險(xiǎn)空間分布Figure 3 Distribution of environmental risk of fertilization in China in 2014

    總之,氮肥施用是目前我國(guó)化肥面源污染的主要風(fēng)險(xiǎn),其次為磷肥,相比來(lái)說(shuō),鉀肥施用環(huán)境風(fēng)險(xiǎn)還不突出。盡管人們還沒(méi)有認(rèn)識(shí)到施鉀對(duì)環(huán)境有什么危害,但是區(qū)域差別大,有上升的潛勢(shì),應(yīng)引起警惕??偟膩?lái)說(shuō),各類(lèi)風(fēng)險(xiǎn)呈現(xiàn)集中連片分布,地域聯(lián)系比較明顯,對(duì)其內(nèi)在規(guī)律性尚需進(jìn)一步探討。

    3 討論

    3.1 化肥施用環(huán)境安全閾值動(dòng)態(tài)模型比單一標(biāo)準(zhǔn)更適合不同區(qū)域情況

    化肥施用環(huán)境安全閾值是農(nóng)田化肥面源污染環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)和管理的關(guān)鍵參數(shù),關(guān)系到化肥面源污染環(huán)境風(fēng)險(xiǎn)控制的成敗。目前普遍認(rèn)為,發(fā)達(dá)國(guó)家為防止化肥污染環(huán)境風(fēng)險(xiǎn),設(shè)置的化肥施用環(huán)境安全上限值為225 kg·hm-2[18],然而亦有文章認(rèn)為這是氮肥施用的環(huán)境安全上限值[19-20]。我國(guó)環(huán)境保護(hù)部制定的生態(tài)鄉(xiāng)鎮(zhèn)建設(shè)指標(biāo)中給出了單位播種面積化肥施用強(qiáng)度不超過(guò)250 kg·hm-2的指標(biāo)。但這些上限值或指標(biāo)都是一個(gè)單一的常數(shù),沒(méi)有考慮指標(biāo)使用的具體環(huán)境條件。其實(shí),化肥施用的環(huán)境安全閾值與具體的生態(tài)環(huán)境和社會(huì)經(jīng)濟(jì)條件有著密切的關(guān)系,采取一刀切的方式,對(duì)不同環(huán)境條件用同一數(shù)值作為化肥施用安全標(biāo)準(zhǔn)顯然是粗糙的,甚至是不合理的。本文提出的化肥施用環(huán)境安全閾值模型是基于當(dāng)?shù)丶Z食作物目標(biāo)產(chǎn)量和相應(yīng)需肥量設(shè)計(jì)的。依據(jù)糧食作物目標(biāo)產(chǎn)量確定施肥安全閾值,既不會(huì)因?yàn)殚撝颠^(guò)低,限制肥料的施用,影響糧食生產(chǎn),也不會(huì)因?yàn)殚撝颠^(guò)高,誤導(dǎo)過(guò)量施肥而污染環(huán)境。建立符合區(qū)域特點(diǎn)的不同自然經(jīng)濟(jì)條件下化肥施用環(huán)境安全閾值測(cè)度模型和環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)方法,不僅是目前化肥施用面源污染風(fēng)險(xiǎn)控制迫切需要解決的一個(gè)理論問(wèn)題,也將對(duì)化肥施用環(huán)境風(fēng)險(xiǎn)田間管理起到重要的指導(dǎo)作用。

    3.2 2014年中國(guó)各省區(qū)環(huán)境風(fēng)險(xiǎn)趨于平衡

    本研究與筆者前期對(duì)2012年中國(guó)化肥施用環(huán)境風(fēng)險(xiǎn)的研究結(jié)果相比[9],全國(guó)化肥施用風(fēng)險(xiǎn)總體平均都是在低度風(fēng)險(xiǎn)范圍,但是各省區(qū)的環(huán)境風(fēng)險(xiǎn)程度趨于平均,環(huán)境安全的省區(qū)由5個(gè)增加到11個(gè),低度風(fēng)險(xiǎn)省區(qū)由9個(gè)增加到16個(gè),中度風(fēng)險(xiǎn)由13個(gè)減少到3個(gè),僅有一個(gè)嚴(yán)重風(fēng)險(xiǎn)省區(qū)。其中最主要的原因是,前期研究采用的環(huán)境安全閾值是國(guó)家關(guān)于生態(tài)縣建設(shè)化肥施用強(qiáng)度不小于250 kg·hm-2常數(shù)值,本研究應(yīng)用環(huán)境安全閾值模型計(jì)算各省區(qū)的環(huán)境安全閾值,把環(huán)境安全閾值與當(dāng)?shù)氐沫h(huán)境條件相結(jié)合,更接近各省的實(shí)際情況,且2014年中國(guó)總體平均環(huán)境安全閾值為285 kg·hm-2,提高了閾值標(biāo)準(zhǔn),風(fēng)險(xiǎn)程度有所降低。

    3.3 化肥施用強(qiáng)度計(jì)算方法急需統(tǒng)一

    化肥施用強(qiáng)度是我國(guó)生態(tài)縣和生態(tài)鄉(xiāng)鎮(zhèn)建設(shè)及評(píng)價(jià)化肥施用環(huán)境風(fēng)險(xiǎn)的一個(gè)重要指標(biāo),但其計(jì)算方法經(jīng)常出現(xiàn)混亂。國(guó)家環(huán)??偩?007年發(fā)布的“關(guān)于印發(fā)《生態(tài)縣、生態(tài)市、生態(tài)省建設(shè)指標(biāo)(修訂稿)》的通知”(環(huán)發(fā)[2007]195號(hào))中關(guān)于化肥施用強(qiáng)度的指標(biāo)計(jì)算是按照耕地面積[10],而后來(lái)的國(guó)家環(huán)保部“關(guān)于印發(fā)《國(guó)家級(jí)生態(tài)鄉(xiāng)鎮(zhèn)申報(bào)及管理規(guī)定(試行)》的通知(環(huán)發(fā)[2010]75號(hào))關(guān)于化肥施用強(qiáng)度指標(biāo)的解釋是“指鄉(xiāng)鎮(zhèn)轄區(qū)內(nèi)實(shí)際用于農(nóng)業(yè)生產(chǎn)的化肥施用量(包括氮肥、磷肥、鉀肥和復(fù)合肥)與播種面積之比”[11]。按照播種面積計(jì)算與按耕地面積計(jì)算的結(jié)果有很大不同。因?yàn)槲覈?guó)各地的復(fù)種指數(shù)差別很大,在復(fù)種指數(shù)高的地區(qū)播種面積與耕地面積會(huì)相差2到3倍[21],所以,在一些地區(qū)用耕地面積計(jì)算出的化肥施用強(qiáng)度與用播種面積計(jì)算的結(jié)果相比,差別會(huì)很大。這表明,化肥施用強(qiáng)度指標(biāo)的計(jì)算問(wèn)題應(yīng)引起研究者和引用者的注意,避免出現(xiàn)錯(cuò)誤的結(jié)論。

    4 結(jié)論與建議

    4.1 結(jié)論

    (1)我國(guó)化肥施用環(huán)境安全閾值近25年來(lái)在200~300 kg·hm-2之間變化,平均值為243 kg·hm-2,近15年來(lái)的平均值為251 kg·hm-2,與國(guó)家環(huán)保部生態(tài)鄉(xiāng)鎮(zhèn)建設(shè)規(guī)定的化肥施用強(qiáng)度標(biāo)準(zhǔn)基本吻合。在1994年我國(guó)化肥施用由環(huán)境安全轉(zhuǎn)為低度風(fēng)險(xiǎn),氮磷鉀單質(zhì)化肥施用分別在1988、1999、2008年由環(huán)境安全轉(zhuǎn)為低度風(fēng)險(xiǎn)??傮w來(lái)說(shuō),自2006年以來(lái),我國(guó)化肥施用環(huán)境風(fēng)險(xiǎn)有所降低。

    (2)2014年我國(guó)化肥施用環(huán)境安全閾值平均為285 kg·hm-2,化肥施用環(huán)境風(fēng)險(xiǎn)指數(shù)為0.54,為低度環(huán)境風(fēng)險(xiǎn);氮磷鉀化肥風(fēng)險(xiǎn)指數(shù)分別為0.54、0.53、0.49,氮磷施用處于低度風(fēng)險(xiǎn),鉀肥處于環(huán)境安全臨界狀態(tài)。

    (3)2014年各省區(qū)的化肥施用環(huán)境風(fēng)險(xiǎn)差異較大,呈現(xiàn)出聚集分布的特點(diǎn)。西藏等11個(gè)省區(qū)處于化肥施用環(huán)境安全狀態(tài),遼寧等16個(gè)省區(qū)處于低度環(huán)境風(fēng)險(xiǎn)狀態(tài),河南、海南和天津三省區(qū)處于中度環(huán)境風(fēng)險(xiǎn),陜西處于嚴(yán)重環(huán)境風(fēng)險(xiǎn);氮肥施用環(huán)境風(fēng)險(xiǎn)分布基本與總肥分布一致,磷肥相對(duì)較平衡,鉀肥差別較大,從環(huán)境安全到極嚴(yán)重風(fēng)險(xiǎn)都有出現(xiàn)。

    4.2 建議

    (1)加快化肥面源污染環(huán)境風(fēng)險(xiǎn)監(jiān)測(cè)監(jiān)控體系建設(shè),建立農(nóng)業(yè)生態(tài)環(huán)境預(yù)警機(jī)制?;拭嬖次廴驹谝恍┑貐^(qū)相當(dāng)嚴(yán)重的原因之一就是對(duì)面源污染的監(jiān)測(cè)預(yù)警能力不足,不能及時(shí)掌握化肥生產(chǎn)、銷(xiāo)售、施用、流失和污染等各方面的信息,及時(shí)地采取防治措施。因此,政府要加大科技、資金、設(shè)備、人員等的投入,建立省、市、縣等各級(jí)面源污染風(fēng)險(xiǎn)監(jiān)測(cè)信息專(zhuān)網(wǎng),特別是化肥面源污染風(fēng)險(xiǎn)大的省區(qū),實(shí)現(xiàn)覆蓋全區(qū)域的農(nóng)業(yè)面源污染數(shù)據(jù)的快速采集、審核、傳輸、反饋,以及對(duì)農(nóng)業(yè)面源污染風(fēng)險(xiǎn)的預(yù)測(cè)預(yù)報(bào),有效掌握農(nóng)業(yè)面源污染風(fēng)險(xiǎn)現(xiàn)狀、來(lái)源和發(fā)展趨勢(shì),為全面防治農(nóng)業(yè)面源污染提供科學(xué)依據(jù)。

    (2)增加科技和教育的投入,提高農(nóng)民的科學(xué)素養(yǎng)和環(huán)境保護(hù)意識(shí)。農(nóng)民是農(nóng)業(yè)生產(chǎn)的主體,也是控制農(nóng)業(yè)化肥過(guò)量施用的重要力量。特別是化肥施用環(huán)境風(fēng)險(xiǎn)較大的地區(qū),應(yīng)高度重視研究農(nóng)戶(hù)行為與化肥面源污染的關(guān)系。對(duì)農(nóng)民進(jìn)行農(nóng)業(yè)生產(chǎn)知識(shí)和環(huán)境保護(hù)知識(shí)的教育和培訓(xùn),使農(nóng)民認(rèn)識(shí)到合理施肥和控制化肥污染的好處。通過(guò)一定的法規(guī)、政策、教育、管理和技術(shù)等措施指導(dǎo)和優(yōu)化農(nóng)戶(hù)施肥行為。農(nóng)業(yè)和土壤科學(xué)的研究要與生產(chǎn)實(shí)踐緊密結(jié)合,研究不同土壤在不同耕作制度下的合理施肥技術(shù)并及時(shí)推廣。只有使面源污染風(fēng)險(xiǎn)控制成為農(nóng)戶(hù)的自覺(jué)行為,才能真正從源頭控制農(nóng)田化肥面源污染環(huán)境風(fēng)險(xiǎn)。

    (3)加快制定土壤污染防治法、土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)和化肥施用技術(shù)規(guī)范及相關(guān)政策等,加強(qiáng)農(nóng)田管理。目前我國(guó)保護(hù)土壤環(huán)境的專(zhuān)項(xiàng)法律基本上還是空白,土壤環(huán)境標(biāo)準(zhǔn)體系也很不健全。國(guó)家應(yīng)該盡快出臺(tái)有關(guān)法律制度和政策,對(duì)土壤環(huán)境保護(hù)做出總體規(guī)劃與部署,對(duì)化肥農(nóng)藥等化學(xué)物質(zhì)的使用也應(yīng)有相應(yīng)的規(guī)定和要求。規(guī)定關(guān)于防治農(nóng)業(yè)化肥污染的基本原則、管理體制、主要權(quán)利和義務(wù)、法律責(zé)任等。探索有機(jī)種植方式補(bǔ)貼機(jī)制及濫用化肥懲處機(jī)制,以制度保障面源污染控制和農(nóng)業(yè)可持續(xù)發(fā)展。

    [1]周亮,徐建剛,孫東琪,等.淮河流域農(nóng)業(yè)面源污染空間特征解析及分類(lèi)控制[J].環(huán)境科學(xué),2013,34(2).547-554.

    ZHOU Liang,XU Jian-gang,SUN Dong-qi,et al.Spatial heterogeneity and classified control of agricultural non-point source pollution in Huaihe River Basin[J].Environmental Science,2013,34(2).547-554.

    [2]Chaplot V.Impact of DEM mesh size and soil map scale on SWAT runoff,sediment,and NO3-N loads predictions[J].J Hydrol,2005,315. 207-222.

    [3]起曉星.基于糧食安全和環(huán)境風(fēng)險(xiǎn)控制目標(biāo)的可持續(xù)農(nóng)業(yè)土地利用模式研究[D].北京.中國(guó)農(nóng)業(yè)大學(xué),2015.6-11.

    QI Xiao-xing.Study on sustainable agricultural land use patterns under the premise of ensuring food security with lower environmental cost[D]. Beijing.China Agricultural University,2015.6-11.

    [4]Wang Z H,Zhao D Z,Cao B,et al.Research on simulation of non-point source pollution in Qingjiang River Basin based on SWAT model and GIS[J].Journal of Yangtze River Scientific Research Institute,2010,27 (1).57-61.

    [5]田耀武,黃志霖,肖文發(fā).基于AnnAGNPS模型的三峽庫(kù)區(qū)秭歸縣非點(diǎn)源污染輸出評(píng)價(jià)[J].生態(tài)學(xué)報(bào),2011,31(16).4568-4578.

    TIAN Yao-wu,HUANG Zhi-lin,XIAO Wen-fa.Assessment of nonpoint source pollution export from Zigui county in the Three Gorges Reservoir area using the AnnAGNPS model[J].Acta Ecologica Sinica, 2011,31(16).4568-4578.

    [6]尹黎明,潘劍軍,王秀英.APEX模型在淮河中上游坡面尺度的適宜性研究[J].中國(guó)水土保持科學(xué),2012,10(2).1-7.

    YIN Li-ming,PAN Jian-jun,WANG Xiu-ying.Applicability of APEX model for slope scale in the middle and upland Huaihe River[J].Science of Soil and Water Conservation,2012,10(2).1-7.

    [7]苑韶峰,呂軍.流域農(nóng)業(yè)非點(diǎn)源污染研究概況[J].土壤通報(bào),2004, 35(4).507-511.

    YUAN Shao-feng,Lü Jun.A general introduction of agricultural nonpoint source pollution in watersheds[J].Chinese Journal of Soil Science, 2004,35(4).507-511.

    [8]陳會(huì),王康周,祖昊.基于排水過(guò)程分析的水稻灌區(qū)農(nóng)田面源污染模擬[J].農(nóng)業(yè)工程學(xué)報(bào),2012,28(6).112-117.

    CHEN Hui,WANG Kang-zhou,ZU Hao.Simulation of agricultural non-point source pollution from paddy rice irrigation district based on analyses of drainage processes[J].Transactions of the Chinese Society of Agricultural Engineering,2012,28(6).112-117.

    [9]劉欽普.中國(guó)化肥投入?yún)^(qū)域差異及環(huán)境風(fēng)險(xiǎn)分析[J].中國(guó)農(nóng)業(yè)科學(xué), 2014,47(18).3596-3605.

    LIU Qin-pu.Distribution of fertilizer application and its environmental risk in different provinces of China[J].ScientiaAgriculturaSinica,Agricultural Science,2014,47(18).3596-3605.

    [10]國(guó)家環(huán)境保護(hù)總局.關(guān)于印發(fā)《生態(tài)縣、生態(tài)市、生態(tài)省建設(shè)指標(biāo)(修訂稿)》的通知[EB/OL].[2016-03-01]http://www.zhb.gov.cn/ gkml/zj/wj/200910/t20091022_172492.htm.

    Ministry of Environmental Protection of P.R.C.Notice on document distribution of“indicators about ecological county,ecological city,ecological province construction(revised)”[EB/OL].[2016-03-01].http:// www.zhb.gov.cn/gkml/zj/wj/200910/t20091022_172492.htm.

    [11]環(huán)境保護(hù)部.關(guān)于印發(fā)《國(guó)家級(jí)生態(tài)鄉(xiāng)鎮(zhèn)申報(bào)及管理規(guī)定(試行)》的通知[EB/OL].[2016-03-01].http://www.wyepb.gov.cn/html/news5/ show_news_w2_1_606.html,

    Ministry of Environmental Protection of P.R.C.Notice on document distribution of“application and management regulation for ecological villages and towns at state level”[EB/OL].[2016-03-01].http://www. wyepb.gov.cn/html/news5/show_news_w2_1_606.html.

    [12]崔玉亭,程序,韓純?nèi)?等.蘇南太湖流域水稻經(jīng)濟(jì)生態(tài)適宜施氮量研究[J].生態(tài)學(xué)報(bào),2000,20(4).659-662.

    CUI Yu-ting,CHENG Xu,HAN Chun-ru,et al.The economic and ecological satisfactory amount of nitrogen fertilizer using on rice in Tai Lake Watershed[J].ActaEcologicaSinica,2000,20(4).659-662.

    [13]巨曉棠.理論施氮量的改進(jìn)及驗(yàn)證.兼論確定作物氮肥推薦量的方法[J].土壤學(xué)報(bào),2015,52(2).1-13.

    JU Xiao-tang.Improvement and validation of theoretical N rate(TNR). Discussing the methods for N fertilizer recommendation[J].Acta PedologicaSinica,2015,52(2).1-13.

    [14]國(guó)家統(tǒng)計(jì)局.國(guó)家數(shù)據(jù)[EB/OL].[2016-03-01].http://data.stats.gov. cn/

    National Bureau of Statistics of China.National data.[EB/OL].[2016-03-01].http://data.stats.gov.cn.

    [15]趙建勛,程燚.安徽省肥料施用現(xiàn)狀與對(duì)策[J].安徽農(nóng)學(xué)通報(bào),2011, 17(1).102-103.

    ZHAO Jian-xun,CHENG Yan.Current situation and countermeasures of fertilization in Anhui Province[J].Anhui Agri Sci Bull,2011,17(1). 102-103.

    [16]宋秀杰,程大軍,張?chǎng)?等.北京種植業(yè)結(jié)構(gòu)調(diào)整及化肥面源污染控制[C].2010中國(guó)環(huán)境科學(xué)學(xué)會(huì)學(xué)術(shù)年會(huì)論文集,2010.第一卷.

    SONG Xiu-jie,CHENG Da-jun,ZHANG Xin,et al.Adjustment of planting industry and control of fertilizer no-point pollution in Beijing[C].Conference proceeding of Chinese Environmental Science Association,2010.first volume.

    [17]國(guó)家統(tǒng)計(jì)局.中國(guó)統(tǒng)計(jì)年鑒2015[EB/OL].[2016-03-01]http://www. stats.gov.cn/tjsj/ndsj/2015/indexch.htm.

    National Bureau of Statistics of China.China statistical yearbook 2015[EB/OL].[2016-03-01]http://www.stats.gov.cn/tjsj/ndsj/2015/indexch.htm.

    [18]孔箐鋅,陳阜.我國(guó)農(nóng)田污染現(xiàn)狀與防治對(duì)策[J].農(nóng)業(yè)工程技術(shù).新能源產(chǎn)業(yè),2010(4).23-26.

    KONG Qing-xin,CHEN Fu.Current pollution situation and countermeasures of cropland in China[J].Agricultural Engineering Technology(Renewable Energy Industry),2010(4).23-26.

    [19]農(nóng)業(yè)部.打好農(nóng)業(yè)面源污染防治攻堅(jiān)戰(zhàn)促進(jìn)農(nóng)業(yè)可持續(xù)發(fā)展.中華人民共和國(guó)國(guó)土資源部(網(wǎng)站)[EB/OL].2015-08-16[2016-03-01].http://search.mlr.gov.cn/was/search

    Ministry of Agriculture of the People’s of Republic of China.Fight against non-point pollution for sustainable development of agriculture of China,Ministry of Land and Resources of the People’s Republic of China(website)[EB/OL].2015-08-16[2016-03-01].http://search. mlr.gov.cn/was/search.

    [20]王奇,詹賢達(dá),王會(huì).我國(guó)糧食安全與水安全之間的關(guān)系初探——基于糧食產(chǎn)量與化肥施用的定量關(guān)系[J].中國(guó)農(nóng)業(yè)資源區(qū)劃, 2013,34(1).81-86.

    WANG Qi,ZHAN Xian-da,WANG Hui.The relationship between food and water environment security in China[J].Chinese Journal of A-gricultural Resources and Regional Planning,2013,34(1).81-86.

    [21]丁明軍,陳倩,辛良杰,等.1999—2013年中國(guó)耕地復(fù)種指數(shù)的時(shí)空演變格局[J].地理學(xué)報(bào),2015,70(7).1080-1090.

    DING Ming-jun,CHEN Qian,XIN Liang-jie,et al.Spatial and temporal variations of multiple cropping index in China based on SPOT-NDVI during 1999—2013[J].Acta Geographica Sinica,2015,70(7). 1080-1090.

    Spatio-temoral changes of fertilization environmental risk of China

    LIU Qin-pu
    (School of Environmental Science,Nanjing Xiaozhuang University,Nanjing 211171,China)

    It is very important for non-point source pollution control to evaluate the risk of non-point source pollution from fertilization. Studies on the spatio-temporal changes of environmental risk of fertilization in China were carried out in order to fully recognize the differences in environmental risk of fertilization at national and provincial levels to accelerate eco-agriculture construction of China.Models of fertilization environmental threshold and risk assessment were established to calculate the fertilization thresholds and assess the fertilization environmental risk of China.Results showed that,the thresholds of fertilization safety were between 200~300 kg·hm-2over past 25 years with an average of 243 kg·hm-2.The fertilization risk changed from safety to low risk degree in 1994;the risks of nitrogen,phosphorus and potash changed to low risk from safety in 1988,1999,2008,respectively;the environmental risk of fertilization has been decreasing since 2006.The environmental risk was at low degree in 2014 with averages of fertilization environmental safety thresholds and risk index of 285 kg·hm-2and 0.54,respectively.The risk indexes of nitrogen,phosphorus and potash were 0.54,0.53,0.49,respectively.The environmental risks at provincial level showed great differences,and the types of risks were aggregately distributed in 2014.There were 11 provincial areas at environmental safety degree,16 areas at low risk degree,three areas of Henan,Hainan and Tianjin at medium risk degree,and the area of Shaanxi at serious risk degree.

    intensity of fertilization;environmental safety threshold;environmental risk index;spatio-temporal change;China

    X820.4

    A

    1672-2043(2017)07-1247-07

    10.11654/jaes.2016-1378

    劉欽普.中國(guó)化肥面源污染環(huán)境風(fēng)險(xiǎn)時(shí)空變化[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2017,36(7).1247-1253.

    LIU Qin-pu.Spatio-temoral changes of fertilization environmental risk of China[J].Journal of Agro-Environment Science,2017,36(7).1247-1253.

    2016-10-29

    劉欽普(1957—),男,河南許昌人,教授,博士,從事土地資源評(píng)價(jià)和生態(tài)經(jīng)濟(jì)方面的研究。E-mail:liuqinpu@163.com

    國(guó)家自然科學(xué)基金項(xiàng)目(31470519);江蘇省自然科學(xué)基金項(xiàng)目(BK20131399)

    Project supported:The National Natural Science Foundation of China(31470519).The Natural Science Foundation of the Jiangsu of China(BK20131399)

    猜你喜歡
    低度面源化肥
    青少年中低度近視控制中低濃度阿托品和角膜塑形鏡的聯(lián)合應(yīng)用效果
    7月我國(guó)化肥進(jìn)口量增加
    中藥日鐘陰陽(yáng)方控制低度青少年近視的效果分析
    農(nóng)業(yè)面源污染的危害與治理
    現(xiàn)代農(nóng)業(yè)離不開(kāi)化肥
    澄江市農(nóng)業(yè)面源污染成因及對(duì)策
    化肥農(nóng)藥減量,怎么減?
    化肥需求增長(zhǎng)不均 有人歡喜有人憂(yōu)
    基于SWAT模型的漳河流域面源污染模擬研究
    農(nóng)業(yè)面源污染對(duì)水質(zhì)的影響及防治對(duì)策
    黄色日韩在线| 欧美成人一区二区免费高清观看| av国产免费在线观看| 一a级毛片在线观看| 日日撸夜夜添| 色噜噜av男人的天堂激情| 男女做爰动态图高潮gif福利片| 99热网站在线观看| 少妇的逼好多水| 午夜亚洲福利在线播放| 日韩高清综合在线| 久99久视频精品免费| 亚洲av电影不卡..在线观看| 亚洲男人的天堂狠狠| 99热精品在线国产| 亚洲专区国产一区二区| 我的老师免费观看完整版| 国产女主播在线喷水免费视频网站 | 最近中文字幕高清免费大全6 | 黄色丝袜av网址大全| 窝窝影院91人妻| 国产成人一区二区在线| 国产午夜精品论理片| 久久久久国产精品人妻aⅴ院| 午夜福利高清视频| 国产午夜精品久久久久久一区二区三区 | 久久国产精品人妻蜜桃| 中国美白少妇内射xxxbb| 欧美区成人在线视频| 日韩欧美 国产精品| 午夜影院日韩av| av国产免费在线观看| 国产精品一区二区三区四区久久| 久久6这里有精品| 久9热在线精品视频| 中文字幕免费在线视频6| 成人美女网站在线观看视频| 国产午夜精品论理片| 成人国产综合亚洲| 人妻夜夜爽99麻豆av| 美女免费视频网站| 变态另类丝袜制服| 日本免费a在线| 亚洲最大成人手机在线| 午夜爱爱视频在线播放| 麻豆精品久久久久久蜜桃| 亚洲人成网站在线播| 深夜精品福利| 色在线成人网| 亚州av有码| 国产精品日韩av在线免费观看| 日韩欧美免费精品| 亚洲最大成人手机在线| 免费看av在线观看网站| 淫妇啪啪啪对白视频| 成人国产综合亚洲| 欧美又色又爽又黄视频| 免费不卡的大黄色大毛片视频在线观看 | 给我免费播放毛片高清在线观看| 一级毛片久久久久久久久女| 久久午夜福利片| bbb黄色大片| 国产男人的电影天堂91| 亚洲精品国产成人久久av| 99国产精品一区二区蜜桃av| 在现免费观看毛片| 99久国产av精品| 国产午夜精品论理片| 男女啪啪激烈高潮av片| 国产精品免费一区二区三区在线| 99久久精品国产国产毛片| 日本黄色片子视频| 免费观看精品视频网站| 免费av不卡在线播放| 亚洲av中文av极速乱 | 天天躁日日操中文字幕| 亚洲五月天丁香| 国产真实乱freesex| 老司机午夜福利在线观看视频| 别揉我奶头 嗯啊视频| 亚洲三级黄色毛片| 特级一级黄色大片| 白带黄色成豆腐渣| 亚洲电影在线观看av| av国产免费在线观看| 99热只有精品国产| 露出奶头的视频| 最近最新免费中文字幕在线| 国产69精品久久久久777片| 欧美潮喷喷水| 免费无遮挡裸体视频| 国产午夜精品久久久久久一区二区三区 | 精品一区二区免费观看| 在现免费观看毛片| 91狼人影院| 精品久久久久久,| 国产黄色小视频在线观看| 麻豆成人午夜福利视频| 久久久久精品国产欧美久久久| 搡老熟女国产l中国老女人| 赤兔流量卡办理| 亚洲欧美日韩卡通动漫| 国产黄色小视频在线观看| 男人的好看免费观看在线视频| 国产精品永久免费网站| 日韩欧美在线乱码| 亚洲欧美日韩东京热| 中文亚洲av片在线观看爽| 人人妻,人人澡人人爽秒播| 他把我摸到了高潮在线观看| 淫秽高清视频在线观看| 男女视频在线观看网站免费| 国产人妻一区二区三区在| 婷婷丁香在线五月| 黄色欧美视频在线观看| 久久精品国产亚洲av涩爱 | 欧美xxxx黑人xx丫x性爽| 性欧美人与动物交配| 日韩欧美在线乱码| 亚洲综合色惰| 国产中年淑女户外野战色| 中国美白少妇内射xxxbb| av国产免费在线观看| 午夜福利在线观看吧| 九色国产91popny在线| 亚洲一区二区三区色噜噜| 国产一区二区在线av高清观看| 非洲黑人性xxxx精品又粗又长| 一个人看视频在线观看www免费| 久久国产乱子免费精品| 国产精华一区二区三区| 成人欧美大片| 一a级毛片在线观看| 嫁个100分男人电影在线观看| 天堂网av新在线| 大型黄色视频在线免费观看| 亚洲aⅴ乱码一区二区在线播放| 久久精品国产亚洲av天美| 国产一级毛片七仙女欲春2| 亚洲精品粉嫩美女一区| 成熟少妇高潮喷水视频| av在线蜜桃| 听说在线观看完整版免费高清| 中出人妻视频一区二区| 不卡一级毛片| 男女之事视频高清在线观看| 色哟哟哟哟哟哟| 婷婷六月久久综合丁香| 午夜精品在线福利| 老司机深夜福利视频在线观看| 久久久久久伊人网av| 国产成人影院久久av| 99久久精品一区二区三区| 国产精品久久久久久亚洲av鲁大| 国产伦人伦偷精品视频| 成人一区二区视频在线观看| 黄片wwwwww| 99riav亚洲国产免费| 天堂av国产一区二区熟女人妻| 国产 一区精品| 欧美xxxx黑人xx丫x性爽| 亚洲欧美日韩高清在线视频| 午夜福利欧美成人| 成人午夜高清在线视频| 特大巨黑吊av在线直播| 成人精品一区二区免费| 熟女电影av网| 亚洲美女视频黄频| 国产伦人伦偷精品视频| 久99久视频精品免费| 内地一区二区视频在线| 成人国产麻豆网| 波多野结衣高清作品| 在线观看午夜福利视频| 又爽又黄无遮挡网站| 在线国产一区二区在线| 夜夜看夜夜爽夜夜摸| 少妇被粗大猛烈的视频| av在线天堂中文字幕| 国产淫片久久久久久久久| 老师上课跳d突然被开到最大视频| 午夜福利在线观看免费完整高清在 | 国产精品1区2区在线观看.| 床上黄色一级片| 久久99热6这里只有精品| 网址你懂的国产日韩在线| 国产国拍精品亚洲av在线观看| 国产亚洲欧美98| 国产真实乱freesex| 亚洲,欧美,日韩| 国内揄拍国产精品人妻在线| 久久九九热精品免费| 日韩在线高清观看一区二区三区 | 好男人在线观看高清免费视频| 国产精品人妻久久久久久| 欧美bdsm另类| 无人区码免费观看不卡| 日日摸夜夜添夜夜添av毛片 | 国产又黄又爽又无遮挡在线| 久久精品人妻少妇| 99久久精品国产国产毛片| 尾随美女入室| 亚洲自偷自拍三级| 欧美3d第一页| 精品人妻偷拍中文字幕| 国模一区二区三区四区视频| av在线亚洲专区| 99久久中文字幕三级久久日本| 国产麻豆成人av免费视频| 熟妇人妻久久中文字幕3abv| 国产成人福利小说| 乱码一卡2卡4卡精品| 色哟哟哟哟哟哟| 精品免费久久久久久久清纯| 国产在线男女| 在线免费观看的www视频| 亚洲欧美激情综合另类| 成人特级黄色片久久久久久久| 成人特级av手机在线观看| 成人二区视频| 神马国产精品三级电影在线观看| 1024手机看黄色片| 亚洲经典国产精华液单| 最近视频中文字幕2019在线8| 黄色丝袜av网址大全| av天堂在线播放| 国模一区二区三区四区视频| 亚洲国产精品合色在线| av在线天堂中文字幕| 国产在线精品亚洲第一网站| 日韩精品有码人妻一区| 成人二区视频| av天堂在线播放| 男插女下体视频免费在线播放| 欧美日本视频| 变态另类丝袜制服| 男女那种视频在线观看| 久久精品夜夜夜夜夜久久蜜豆| 午夜激情福利司机影院| 两个人的视频大全免费| 两个人视频免费观看高清| av专区在线播放| 精品国产三级普通话版| 成人鲁丝片一二三区免费| 亚洲电影在线观看av| 日本a在线网址| 联通29元200g的流量卡| 99热只有精品国产| 国产av一区在线观看免费| 亚洲美女黄片视频| 一进一出抽搐动态| 我的老师免费观看完整版| 久久亚洲真实| 黄色配什么色好看| 国产精品精品国产色婷婷| 九色国产91popny在线| 黄片wwwwww| 国产中年淑女户外野战色| 亚洲中文字幕日韩| 精品久久久久久,| 性色avwww在线观看| 韩国av一区二区三区四区| 午夜a级毛片| 欧美色欧美亚洲另类二区| 九九在线视频观看精品| 少妇猛男粗大的猛烈进出视频 | 尾随美女入室| 色噜噜av男人的天堂激情| 午夜免费激情av| 美女cb高潮喷水在线观看| 天堂影院成人在线观看| 成人av一区二区三区在线看| 国产精品爽爽va在线观看网站| 欧美色欧美亚洲另类二区| 搡老岳熟女国产| 精品乱码久久久久久99久播| 桃红色精品国产亚洲av| 国产aⅴ精品一区二区三区波| 国产熟女欧美一区二区| 欧美黑人欧美精品刺激| 一本一本综合久久| eeuss影院久久| 99久久无色码亚洲精品果冻| 九九爱精品视频在线观看| 国产老妇女一区| 天堂动漫精品| 久久久久性生活片| 成人高潮视频无遮挡免费网站| 亚洲精品成人久久久久久| 舔av片在线| 午夜福利成人在线免费观看| 国产高清不卡午夜福利| 淫秽高清视频在线观看| 午夜福利成人在线免费观看| ponron亚洲| 极品教师在线免费播放| 老司机午夜福利在线观看视频| 亚洲av五月六月丁香网| 国产精品免费一区二区三区在线| 99在线视频只有这里精品首页| 少妇的逼水好多| 色噜噜av男人的天堂激情| 亚洲四区av| 国产精品不卡视频一区二区| 精品国产三级普通话版| 国产真实乱freesex| 国产伦在线观看视频一区| 亚洲美女视频黄频| 亚洲成人久久爱视频| 国产毛片a区久久久久| 丰满人妻一区二区三区视频av| 日韩亚洲欧美综合| 真人做人爱边吃奶动态| 免费看a级黄色片| 99久久无色码亚洲精品果冻| 男女那种视频在线观看| 午夜日韩欧美国产| 最好的美女福利视频网| 99国产精品一区二区蜜桃av| 日本一本二区三区精品| 久久热精品热| av视频在线观看入口| 中文在线观看免费www的网站| 女同久久另类99精品国产91| 久久久久久大精品| 99久久成人亚洲精品观看| 97人妻精品一区二区三区麻豆| 国产日本99.免费观看| 又粗又爽又猛毛片免费看| 日本欧美国产在线视频| 欧美三级亚洲精品| 一本久久中文字幕| 精品一区二区三区视频在线观看免费| 高清毛片免费观看视频网站| 国产精品亚洲一级av第二区| 久久久久国内视频| 夜夜看夜夜爽夜夜摸| 最后的刺客免费高清国语| 国语自产精品视频在线第100页| 两个人视频免费观看高清| 欧美+日韩+精品| 黄色视频,在线免费观看| 麻豆一二三区av精品| 五月伊人婷婷丁香| 1024手机看黄色片| 熟妇人妻久久中文字幕3abv| 女生性感内裤真人,穿戴方法视频| 精品久久久久久久久久免费视频| 波多野结衣高清作品| 久久久久久久亚洲中文字幕| 国产亚洲91精品色在线| 日韩大尺度精品在线看网址| 精品久久国产蜜桃| 小蜜桃在线观看免费完整版高清| 中文字幕熟女人妻在线| 亚洲成av人片在线播放无| 亚洲精品国产成人久久av| 麻豆一二三区av精品| 成年女人永久免费观看视频| 国产欧美日韩精品一区二区| 国产精品日韩av在线免费观看| 美女被艹到高潮喷水动态| 亚洲熟妇熟女久久| 国产av麻豆久久久久久久| 88av欧美| 天堂网av新在线| 国产精品一区二区三区四区久久| 欧美潮喷喷水| 狂野欧美激情性xxxx在线观看| 中文字幕高清在线视频| 国国产精品蜜臀av免费| 亚洲国产日韩欧美精品在线观看| 亚洲性夜色夜夜综合| 中文字幕高清在线视频| 国产精品女同一区二区软件 | 国产成人福利小说| 精品久久久久久久久久久久久| 亚洲美女黄片视频| 男人和女人高潮做爰伦理| 国产午夜精品论理片| 亚洲国产欧洲综合997久久,| 免费看av在线观看网站| 久久精品国产亚洲av香蕉五月| 成人性生交大片免费视频hd| 国产白丝娇喘喷水9色精品| 国产成人福利小说| 欧洲精品卡2卡3卡4卡5卡区| 午夜免费成人在线视频| 久久九九热精品免费| 久久久久国产精品人妻aⅴ院| 变态另类丝袜制服| 亚洲不卡免费看| 97热精品久久久久久| 蜜桃亚洲精品一区二区三区| 成年女人看的毛片在线观看| 亚洲七黄色美女视频| 国产精品野战在线观看| 一区福利在线观看| 久9热在线精品视频| 国产乱人伦免费视频| 亚洲最大成人av| 又爽又黄a免费视频| 亚洲自拍偷在线| 国产精品一区www在线观看 | 精品免费久久久久久久清纯| 日韩一区二区视频免费看| 久久精品国产亚洲av香蕉五月| 99热这里只有是精品在线观看| 日韩中文字幕欧美一区二区| 欧美激情在线99| 中国美白少妇内射xxxbb| 国产亚洲91精品色在线| 欧洲精品卡2卡3卡4卡5卡区| 亚洲真实伦在线观看| 国产淫片久久久久久久久| 波野结衣二区三区在线| 99九九线精品视频在线观看视频| 亚洲国产日韩欧美精品在线观看| 欧美成人性av电影在线观看| 一个人看的www免费观看视频| 久久久久国内视频| 婷婷精品国产亚洲av| 舔av片在线| 窝窝影院91人妻| 亚洲18禁久久av| 美女被艹到高潮喷水动态| 99久国产av精品| 欧美日韩黄片免| 亚洲欧美日韩高清专用| 亚洲avbb在线观看| 久久这里只有精品中国| 久久人妻av系列| 欧美精品国产亚洲| 看黄色毛片网站| 欧美一区二区亚洲| 国产精品1区2区在线观看.| 久久久久九九精品影院| 国产高潮美女av| 亚洲人与动物交配视频| 97人妻精品一区二区三区麻豆| 日韩欧美在线二视频| 成人三级黄色视频| 精品欧美国产一区二区三| 亚洲成人久久性| 午夜激情欧美在线| 麻豆精品久久久久久蜜桃| 99国产极品粉嫩在线观看| 女人被狂操c到高潮| 色噜噜av男人的天堂激情| 99久久成人亚洲精品观看| 亚洲av免费在线观看| 国内精品久久久久精免费| 九九热线精品视视频播放| 国产白丝娇喘喷水9色精品| 悠悠久久av| xxxwww97欧美| 久久精品久久久久久噜噜老黄 | 男女做爰动态图高潮gif福利片| 久久国产精品人妻蜜桃| 91久久精品国产一区二区三区| 亚洲性久久影院| 岛国在线免费视频观看| 欧美三级亚洲精品| 韩国av在线不卡| 亚洲精品一区av在线观看| 国产美女午夜福利| av天堂中文字幕网| 精品一区二区三区视频在线观看免费| 亚洲七黄色美女视频| 不卡视频在线观看欧美| 免费看美女性在线毛片视频| 午夜免费男女啪啪视频观看 | 综合色av麻豆| 老女人水多毛片| 亚洲欧美日韩无卡精品| 18禁黄网站禁片午夜丰满| 久久精品91蜜桃| 亚洲精华国产精华精| 亚洲中文日韩欧美视频| 国产亚洲91精品色在线| 精品久久久噜噜| 在线免费观看的www视频| 老师上课跳d突然被开到最大视频| 久久99热这里只有精品18| 99热这里只有精品一区| 一级a爱片免费观看的视频| 国产亚洲欧美98| 精品人妻1区二区| 国产精品一区www在线观看 | 欧美成人a在线观看| 国产精品亚洲美女久久久| 国产精品精品国产色婷婷| 91在线观看av| 嫩草影院精品99| 色综合婷婷激情| 床上黄色一级片| 欧美黑人巨大hd| 国产在视频线在精品| 香蕉av资源在线| 校园春色视频在线观看| 国产不卡一卡二| 别揉我奶头~嗯~啊~动态视频| 国产午夜福利久久久久久| 性色avwww在线观看| 国语自产精品视频在线第100页| 少妇猛男粗大的猛烈进出视频 | 色综合亚洲欧美另类图片| 欧美中文日本在线观看视频| 日本与韩国留学比较| 51国产日韩欧美| 欧美一区二区精品小视频在线| 69av精品久久久久久| 日韩精品中文字幕看吧| 国产精品精品国产色婷婷| 欧美+日韩+精品| 嫩草影院新地址| 国产精品不卡视频一区二区| 长腿黑丝高跟| 麻豆一二三区av精品| 亚洲内射少妇av| 国产精品野战在线观看| 国产一区二区三区在线臀色熟女| 在线播放无遮挡| 一夜夜www| 91狼人影院| 69人妻影院| 午夜日韩欧美国产| 99九九线精品视频在线观看视频| 成人精品一区二区免费| a级毛片免费高清观看在线播放| 成人永久免费在线观看视频| 欧美性猛交╳xxx乱大交人| av在线亚洲专区| 欧美性感艳星| av在线亚洲专区| 在线观看美女被高潮喷水网站| 亚洲无线观看免费| 日本免费a在线| 一个人看的www免费观看视频| 日韩高清综合在线| 动漫黄色视频在线观看| 国产亚洲精品久久久com| 一个人看的www免费观看视频| 日韩欧美精品免费久久| 亚洲一区二区三区色噜噜| 亚洲人与动物交配视频| 欧美+亚洲+日韩+国产| aaaaa片日本免费| 最好的美女福利视频网| 国产69精品久久久久777片| 少妇被粗大猛烈的视频| 精品99又大又爽又粗少妇毛片 | 欧洲精品卡2卡3卡4卡5卡区| 亚洲美女视频黄频| 九色国产91popny在线| 久久久久久久午夜电影| 国产精品野战在线观看| 俄罗斯特黄特色一大片| 尾随美女入室| 在线观看美女被高潮喷水网站| 国产精品国产高清国产av| 国产精品久久视频播放| 精品一区二区三区视频在线观看免费| 天天躁日日操中文字幕| 看十八女毛片水多多多| 久久天躁狠狠躁夜夜2o2o| 国产熟女欧美一区二区| 日本一二三区视频观看| 国内精品美女久久久久久| xxxwww97欧美| 午夜福利欧美成人| 男女那种视频在线观看| 男女边吃奶边做爰视频| 赤兔流量卡办理| 日本a在线网址| 国语自产精品视频在线第100页| 久久国内精品自在自线图片| 一区二区三区激情视频| 亚洲人与动物交配视频| 午夜免费成人在线视频| 一卡2卡三卡四卡精品乱码亚洲| 国产免费av片在线观看野外av| 十八禁国产超污无遮挡网站| 亚洲熟妇中文字幕五十中出| 91在线观看av| 欧美激情在线99| 精品不卡国产一区二区三区| 免费观看的影片在线观看| 在线观看av片永久免费下载| 午夜激情福利司机影院| 男女下面进入的视频免费午夜| 午夜久久久久精精品| 国产成人福利小说| 日韩国内少妇激情av| 色精品久久人妻99蜜桃| 在线观看av片永久免费下载| 免费看av在线观看网站| 丰满人妻一区二区三区视频av| 久久久久久久久中文| 国产高清有码在线观看视频| 永久网站在线| 国产视频内射| 一区福利在线观看| 国产在线精品亚洲第一网站| 成人特级av手机在线观看| 亚洲最大成人中文| 亚洲四区av| av在线观看视频网站免费| 国内久久婷婷六月综合欲色啪| 亚洲国产欧洲综合997久久,| 亚洲中文日韩欧美视频| 亚洲国产精品成人综合色| 亚洲国产高清在线一区二区三| 久9热在线精品视频| 人人妻人人澡欧美一区二区|