• <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免费观看日本| 自拍欧美九色日韩亚洲蝌蚪91| 18+在线观看网站| 日韩 亚洲 欧美在线| 综合色丁香网| 99热国产这里只有精品6| 久久av网站| 亚洲欧美成人综合另类久久久| 99热全是精品| 亚洲精品乱久久久久久| 久久精品久久久久久噜噜老黄| 欧美日韩成人在线一区二区| 亚洲精品一区蜜桃| 99热这里只有是精品在线观看| 又黄又爽又刺激的免费视频.| 我要看黄色一级片免费的| 老司机影院毛片| 飞空精品影院首页| 免费黄网站久久成人精品| 成年人午夜在线观看视频| 国产综合精华液| 国产精品一国产av| 亚洲一区二区三区欧美精品| 亚洲精品美女久久av网站| 亚洲欧美成人精品一区二区| 亚洲三级黄色毛片| 免费观看性生交大片5| 一边亲一边摸免费视频| 亚洲精品乱码久久久久久按摩| 黑丝袜美女国产一区| 国产在线免费精品| 宅男免费午夜| 亚洲图色成人| 晚上一个人看的免费电影| 久久午夜综合久久蜜桃| 久久久久人妻精品一区果冻| 日日摸夜夜添夜夜爱| 国产熟女欧美一区二区| 国产日韩欧美亚洲二区| 侵犯人妻中文字幕一二三四区| 在线 av 中文字幕| 99九九在线精品视频| 国产综合精华液| 欧美日韩国产mv在线观看视频| 欧美日韩视频高清一区二区三区二| 精品久久国产蜜桃| 精品国产一区二区久久| 久久99精品国语久久久| 永久免费av网站大全| 久久久亚洲精品成人影院| 欧美人与善性xxx| 国产一区有黄有色的免费视频| 美女国产高潮福利片在线看| 精品卡一卡二卡四卡免费| 多毛熟女@视频| 欧美亚洲 丝袜 人妻 在线| 精品少妇黑人巨大在线播放| 韩国精品一区二区三区 | √禁漫天堂资源中文www| 日本wwww免费看| 亚洲综合色网址| 午夜免费鲁丝| 美女国产高潮福利片在线看| 亚洲av在线观看美女高潮| 在线观看一区二区三区激情| 亚洲精品美女久久av网站| 国产成人精品无人区| 一级a做视频免费观看| 亚洲高清免费不卡视频| 一级毛片我不卡| 中国国产av一级| 在线亚洲精品国产二区图片欧美| 男男h啪啪无遮挡| 黑人猛操日本美女一级片| 日本午夜av视频| 精品国产一区二区三区久久久樱花| 涩涩av久久男人的天堂| 韩国av在线不卡| 国语对白做爰xxxⅹ性视频网站| 视频在线观看一区二区三区| 欧美精品高潮呻吟av久久| 麻豆精品久久久久久蜜桃| av在线播放精品| 国产一区亚洲一区在线观看| 国产福利在线免费观看视频| 精品久久久精品久久久| 又大又黄又爽视频免费| 国产av精品麻豆| 精品一品国产午夜福利视频| 狠狠精品人妻久久久久久综合| 欧美成人午夜免费资源| 女的被弄到高潮叫床怎么办| 精品少妇内射三级| 一区在线观看完整版| 国产精品.久久久| 日韩在线高清观看一区二区三区| 国产在视频线精品| 国产精品一区二区在线不卡| 少妇被粗大猛烈的视频| 国产白丝娇喘喷水9色精品| 日韩电影二区| 国产成人精品一,二区| 亚洲图色成人| 美女国产高潮福利片在线看| 91精品伊人久久大香线蕉| 99热网站在线观看| 香蕉丝袜av| 80岁老熟妇乱子伦牲交| 午夜免费观看性视频| 日韩精品免费视频一区二区三区 | 亚洲美女视频黄频| 国产免费一级a男人的天堂| 在线观看一区二区三区激情| 精品一品国产午夜福利视频| 亚洲成人av在线免费| 汤姆久久久久久久影院中文字幕| 国产成人精品无人区| 久久女婷五月综合色啪小说| 尾随美女入室| 国产高清三级在线| videossex国产| 免费大片18禁| 国产精品三级大全| 男女啪啪激烈高潮av片| 大陆偷拍与自拍| 女人久久www免费人成看片| 欧美3d第一页| 99久久精品国产国产毛片| 观看美女的网站| 一边亲一边摸免费视频| 国产精品欧美亚洲77777| 久久ye,这里只有精品| 制服人妻中文乱码| 精品久久蜜臀av无| 日韩av在线免费看完整版不卡| 成人免费观看视频高清| 亚洲欧美色中文字幕在线| 国产日韩欧美在线精品| 老司机影院成人| 丰满乱子伦码专区| 777米奇影视久久| 99热网站在线观看| 日本wwww免费看| 成人国产麻豆网| 成人国产av品久久久| 又黄又爽又刺激的免费视频.| 欧美成人午夜精品| 久久久久久久大尺度免费视频| av在线app专区| 极品少妇高潮喷水抽搐| 国产成人欧美| 丝袜美足系列| 十分钟在线观看高清视频www| 90打野战视频偷拍视频| 国产又色又爽无遮挡免| 精品一区二区三卡| 少妇被粗大的猛进出69影院 | 亚洲精品久久久久久婷婷小说| 亚洲综合色惰| 亚洲,一卡二卡三卡| 亚洲国产色片| 激情视频va一区二区三区| 韩国av在线不卡| 99九九在线精品视频| 免费久久久久久久精品成人欧美视频 | 十分钟在线观看高清视频www| 国产一区二区在线观看av| 日韩成人av中文字幕在线观看| 侵犯人妻中文字幕一二三四区| 日本欧美国产在线视频| 男男h啪啪无遮挡| 91久久精品国产一区二区三区| 一本—道久久a久久精品蜜桃钙片| 欧美人与性动交α欧美精品济南到 | 三上悠亚av全集在线观看| 成年av动漫网址| 国产精品人妻久久久久久| 亚洲精品乱码久久久久久按摩| 成人国产av品久久久| 一区二区三区乱码不卡18| 一边亲一边摸免费视频| 亚洲一级一片aⅴ在线观看| 久久久久视频综合| 一级毛片电影观看| 免费久久久久久久精品成人欧美视频 | 久久久久精品人妻al黑| 欧美精品一区二区免费开放| 国产高清国产精品国产三级| 国产日韩欧美视频二区| 内地一区二区视频在线| 熟妇人妻不卡中文字幕| 国产精品久久久久久久久免| 毛片一级片免费看久久久久| 日产精品乱码卡一卡2卡三| 99精国产麻豆久久婷婷| 成年女人在线观看亚洲视频| 免费人成在线观看视频色| 免费在线观看完整版高清| 亚洲精品美女久久av网站| 十八禁网站网址无遮挡| 女性生殖器流出的白浆| 精品少妇黑人巨大在线播放| 久久av网站| 亚洲久久久国产精品| 成人免费观看视频高清| 成人亚洲欧美一区二区av| 成人免费观看视频高清| 成人亚洲欧美一区二区av| 美女主播在线视频| 啦啦啦视频在线资源免费观看| 国产精品蜜桃在线观看| 看十八女毛片水多多多| 色婷婷av一区二区三区视频| 亚洲国产精品一区三区| 一二三四中文在线观看免费高清| 少妇猛男粗大的猛烈进出视频| 久久久久人妻精品一区果冻| 国产精品欧美亚洲77777| 一级a做视频免费观看| √禁漫天堂资源中文www| 久久热在线av| xxxhd国产人妻xxx| 咕卡用的链子| 91成人精品电影| 久久久精品免费免费高清| 九色成人免费人妻av| 少妇的逼水好多| 国产色婷婷99| 国产成人精品一,二区| 精品亚洲成国产av| 亚洲欧美成人精品一区二区| 美女脱内裤让男人舔精品视频| 69精品国产乱码久久久| 久久影院123| 欧美日韩亚洲高清精品| 性高湖久久久久久久久免费观看| www.色视频.com| 色婷婷久久久亚洲欧美| 国产色婷婷99| 日韩av在线免费看完整版不卡| 国产男女内射视频| 欧美丝袜亚洲另类| 国产毛片在线视频| 制服诱惑二区| 男男h啪啪无遮挡| 久久人人97超碰香蕉20202| 色94色欧美一区二区| 哪个播放器可以免费观看大片| 国产日韩一区二区三区精品不卡| a级毛色黄片| 精品亚洲成a人片在线观看| 久久精品人人爽人人爽视色| 亚洲国产精品专区欧美| 美国免费a级毛片| 久久国产精品男人的天堂亚洲 | 婷婷色麻豆天堂久久| 免费大片18禁| 黄色 视频免费看| 成人免费观看视频高清| 亚洲精品乱久久久久久| 国产精品久久久久久精品电影小说| 国产精品不卡视频一区二区| 日本色播在线视频| 国产1区2区3区精品| 亚洲av在线观看美女高潮| 日本免费在线观看一区| 一级片'在线观看视频| 黄片播放在线免费| 香蕉精品网在线| a级片在线免费高清观看视频| 黄片无遮挡物在线观看| av卡一久久| 丁香六月天网| 91在线精品国自产拍蜜月| 妹子高潮喷水视频| 赤兔流量卡办理| 一区在线观看完整版| 男男h啪啪无遮挡| 一级爰片在线观看| 亚洲一码二码三码区别大吗| www.熟女人妻精品国产 | 中国美白少妇内射xxxbb| 哪个播放器可以免费观看大片| 久久午夜福利片| 日本欧美视频一区| 蜜臀久久99精品久久宅男| 人人澡人人妻人| 九九爱精品视频在线观看| 熟女电影av网| 少妇的逼好多水| 国产成人欧美| 久久女婷五月综合色啪小说| 国产精品免费大片| 国产无遮挡羞羞视频在线观看| 高清视频免费观看一区二区| 免费在线观看黄色视频的| av国产久精品久网站免费入址| 免费人成在线观看视频色| 亚洲一码二码三码区别大吗| 欧美性感艳星| 日韩制服丝袜自拍偷拍| 男女免费视频国产| 午夜久久久在线观看| 好男人视频免费观看在线| 亚洲人与动物交配视频| 国产一级毛片在线| 国产高清三级在线| 国产探花极品一区二区| 国产成人精品无人区| 免费播放大片免费观看视频在线观看| 搡老乐熟女国产| 国产精品嫩草影院av在线观看| 韩国高清视频一区二区三区| 赤兔流量卡办理| 最近最新中文字幕大全免费视频 | 亚洲国产色片| 少妇高潮的动态图| 大码成人一级视频| 秋霞在线观看毛片| 大香蕉久久网| 最近最新中文字幕大全免费视频 | 精品熟女少妇av免费看| 亚洲精品乱久久久久久| 黄色视频在线播放观看不卡| 黄色配什么色好看| 国产精品三级大全| 高清毛片免费看| 精品人妻偷拍中文字幕| 精品少妇内射三级| 国产精品久久久av美女十八| 婷婷色综合大香蕉| 男人爽女人下面视频在线观看| 午夜免费男女啪啪视频观看| 免费播放大片免费观看视频在线观看| 51国产日韩欧美| 久久久久久久久久久免费av| 午夜久久久在线观看| 妹子高潮喷水视频| 免费人妻精品一区二区三区视频| 亚洲婷婷狠狠爱综合网| 久久久欧美国产精品| 草草在线视频免费看| 国产精品99久久99久久久不卡 | 最黄视频免费看| av不卡在线播放| 色5月婷婷丁香| 99热网站在线观看| 国产乱来视频区| 一级毛片黄色毛片免费观看视频| 久久热在线av| 一个人免费看片子| 成人黄色视频免费在线看| 一本大道久久a久久精品| 国产成人一区二区在线| 女人精品久久久久毛片| 国产精品一二三区在线看| 波野结衣二区三区在线| 男女高潮啪啪啪动态图| 极品少妇高潮喷水抽搐| 成人黄色视频免费在线看| 妹子高潮喷水视频| 男人操女人黄网站| 香蕉精品网在线| av免费观看日本| 香蕉精品网在线| 国产熟女午夜一区二区三区| 亚洲精品456在线播放app| 日本与韩国留学比较| 激情视频va一区二区三区| 91国产中文字幕| 欧美3d第一页| 美女主播在线视频| 18禁观看日本| 精品一区二区免费观看| 乱人伦中国视频| 日本wwww免费看| 狂野欧美激情性bbbbbb| 色网站视频免费| 国产视频首页在线观看| 日韩电影二区| 日本色播在线视频| 美女国产高潮福利片在线看| 亚洲成色77777| 国产激情久久老熟女| 免费黄网站久久成人精品| 亚洲中文av在线| 一级毛片电影观看| 乱码一卡2卡4卡精品| 国产黄频视频在线观看| 中文乱码字字幕精品一区二区三区| 99re6热这里在线精品视频| av免费观看日本| 久久精品夜色国产| 在线天堂最新版资源| 日韩成人av中文字幕在线观看| 青春草视频在线免费观看| 亚洲熟女精品中文字幕| 欧美日韩av久久| 美女脱内裤让男人舔精品视频| 免费观看性生交大片5| 日韩欧美一区视频在线观看| 99热网站在线观看| 婷婷色av中文字幕| 天天操日日干夜夜撸| 男男h啪啪无遮挡| 亚洲成av片中文字幕在线观看 | 夫妻性生交免费视频一级片| av免费在线看不卡| 超碰97精品在线观看| 亚洲欧美精品自产自拍| 久久久久久久精品精品| 97在线人人人人妻| 99香蕉大伊视频| 美女国产高潮福利片在线看| 色网站视频免费| 两个人看的免费小视频| 亚洲精品久久成人aⅴ小说| av免费在线看不卡| 蜜桃国产av成人99| 人成视频在线观看免费观看| 黑人猛操日本美女一级片| 国产乱人偷精品视频| 欧美精品一区二区免费开放| 九色亚洲精品在线播放| 免费女性裸体啪啪无遮挡网站| 国产精品久久久久久av不卡| 宅男免费午夜| 免费av中文字幕在线| 日本欧美视频一区| 亚洲av.av天堂| 婷婷成人精品国产| 熟女电影av网| 国产av国产精品国产| 午夜福利,免费看| 九色成人免费人妻av| 精品亚洲成a人片在线观看| 观看av在线不卡| 色婷婷久久久亚洲欧美| 亚洲欧美成人精品一区二区| 激情五月婷婷亚洲| 精品一区二区三卡| 国产成人免费观看mmmm| 亚洲国产成人一精品久久久| 亚洲精品乱码久久久久久按摩| 久久热在线av| 欧美激情 高清一区二区三区| 中文精品一卡2卡3卡4更新| 精品亚洲乱码少妇综合久久| 观看av在线不卡| 丰满少妇做爰视频| 亚洲精品国产av蜜桃| 亚洲av中文av极速乱| 五月天丁香电影| 久久人妻熟女aⅴ| 咕卡用的链子| 22中文网久久字幕| 大香蕉97超碰在线| a级毛片在线看网站| 9热在线视频观看99| 男女免费视频国产| 老司机影院成人| 国产成人精品无人区| 美女国产高潮福利片在线看| 亚洲欧美精品自产自拍| 国产精品熟女久久久久浪| 亚洲人成77777在线视频| 免费观看无遮挡的男女| 80岁老熟妇乱子伦牲交| 一区二区三区四区激情视频| 久久鲁丝午夜福利片| 日本av手机在线免费观看| 中国美白少妇内射xxxbb| 少妇人妻 视频| 国产精品久久久久久精品电影小说| 人成视频在线观看免费观看| 黄色 视频免费看| 搡女人真爽免费视频火全软件| 亚洲图色成人| 日本av免费视频播放| 久久久a久久爽久久v久久| 一级毛片黄色毛片免费观看视频| 午夜激情久久久久久久| av在线老鸭窝| 一本—道久久a久久精品蜜桃钙片| 国产精品嫩草影院av在线观看| 全区人妻精品视频| 午夜福利影视在线免费观看| 久久人妻熟女aⅴ| videossex国产| 久久久久久久久久久免费av| 亚洲 欧美一区二区三区| 亚洲精品美女久久av网站| 黄片播放在线免费| 国产成人午夜福利电影在线观看| 高清av免费在线| 亚洲国产欧美日韩在线播放| 久久久久国产网址| 亚洲国产精品专区欧美| 26uuu在线亚洲综合色| 国产精品无大码| 丝袜美足系列| 国产不卡av网站在线观看| 中文天堂在线官网| 久久国产亚洲av麻豆专区| 亚洲精品久久成人aⅴ小说| 一区二区三区精品91| 大片免费播放器 马上看| 午夜免费男女啪啪视频观看| 日韩精品免费视频一区二区三区 | 97在线视频观看| 成人综合一区亚洲| √禁漫天堂资源中文www| a级片在线免费高清观看视频| 免费大片18禁| 久久久久久久久久久免费av| 一本久久精品| 国产福利在线免费观看视频| 中文乱码字字幕精品一区二区三区| 成人毛片a级毛片在线播放| 一二三四在线观看免费中文在 | 亚洲精品国产色婷婷电影| 欧美激情 高清一区二区三区| 久久久久久久精品精品| 色94色欧美一区二区| 高清欧美精品videossex| 伊人亚洲综合成人网| 欧美精品亚洲一区二区| 肉色欧美久久久久久久蜜桃| 大香蕉97超碰在线| 亚洲av日韩在线播放| a 毛片基地| 亚洲精品视频女| 美女视频免费永久观看网站| 亚洲欧美成人综合另类久久久| 乱码一卡2卡4卡精品| 精品午夜福利在线看| 精品人妻在线不人妻| 亚洲av电影在线观看一区二区三区| 18+在线观看网站| 国产精品一区二区在线不卡| 亚洲欧美成人精品一区二区| 美女xxoo啪啪120秒动态图| 亚洲国产精品999| 中国国产av一级| 伦精品一区二区三区| 激情五月婷婷亚洲| 另类亚洲欧美激情| 男女国产视频网站| 成年动漫av网址| xxx大片免费视频| 国产精品久久久av美女十八| 少妇人妻久久综合中文| 91成人精品电影| 亚洲精品,欧美精品| 久久久久国产精品人妻一区二区| 国产黄色免费在线视频| 99热这里只有是精品在线观看| freevideosex欧美| 国产乱来视频区| 2022亚洲国产成人精品| 免费黄色在线免费观看| 边亲边吃奶的免费视频| 97精品久久久久久久久久精品| 蜜臀久久99精品久久宅男| 美女国产高潮福利片在线看| 91午夜精品亚洲一区二区三区| 亚洲欧美一区二区三区国产| 国产欧美日韩综合在线一区二区| 国产精品国产av在线观看| 丰满少妇做爰视频| 国产成人aa在线观看| 少妇被粗大猛烈的视频| 少妇的逼水好多| 国语对白做爰xxxⅹ性视频网站| 国产精品无大码| 精品视频人人做人人爽| 王馨瑶露胸无遮挡在线观看| a级毛片在线看网站| 精品亚洲成a人片在线观看| 色视频在线一区二区三区| 国产福利在线免费观看视频| 欧美成人午夜免费资源| 丰满迷人的少妇在线观看| 国产免费现黄频在线看| 51国产日韩欧美| 成人国语在线视频| 国产欧美亚洲国产| 在线免费观看不下载黄p国产| 日韩大片免费观看网站| 亚洲在久久综合| 国产精品99久久99久久久不卡 | 五月伊人婷婷丁香| 国产av一区二区精品久久| videossex国产| 欧美日韩视频精品一区| 久久精品aⅴ一区二区三区四区 | 一边摸一边做爽爽视频免费| 爱豆传媒免费全集在线观看| 国产在线免费精品| 日韩一区二区三区影片| 国产成人精品婷婷| 欧美最新免费一区二区三区| 亚洲av综合色区一区| √禁漫天堂资源中文www| 国产av码专区亚洲av| 女的被弄到高潮叫床怎么办| 国产亚洲最大av| 你懂的网址亚洲精品在线观看| 中文字幕亚洲精品专区| 久久久精品94久久精品| 熟女电影av网| 纯流量卡能插随身wifi吗| 青春草亚洲视频在线观看| 在线看a的网站| 97在线人人人人妻|